Keywords
Vitamin D Receptor, Polymorphism, Fok I genotypes, Tuberculosis, Uganda Introduction
Vitamin D Receptor, Polymorphism, Fok I genotypes, Tuberculosis, Uganda Introduction
According to the World Health Organization (WHO) report on the ‘Use of high burden country lists for TB by WHO in the post-2015 era (WHO/HTM/TB/2015.29), Uganda was removed from the 22 high TB burden countries. However Uganda is still among the 41 high TB/human immunodeficiency virus (HIV) burden countries1 and TB remains a major public health problem. Although it is still unclear how individuals develop active TB2 there is emerging evidence of multi drug resistant TB (MDRTB). Genetic predisposition to TB has been suggested in several studies3,4 but the immunological associations with genetic polymorphisms are still unclear5. Vitamin D insufficiency is common worldwide and recent studies have shown that vitamin D insufficiency is associated with a higher risk of active TB. The vitamin D receptor gene (VDR) has been identified as a candidate gene for TB susceptibility, however; studies reporting from different ethnic groups have been inconsistent6. Various VDR polymorphisms have been identified and associated with TB susceptibility or resistance2,7,8. Similarly the widely studied and functional single nucleotide polymorphism (SNP) rs2228570 of VDR fok I gene has shown variations. This polymorphism is due to the alteration in one of the start codons in which the amino acid thymine is replaced with cytosine (T/C). The dominant homozygous FF variant has high transcriptional activity with three amino acids less since translation starts at the second codon, while the homozygous recessive (ff) has 427 amino acids and is the longer form. The absence or presence of the restriction site is designated as F or f respectively7,11. This study investigated the frequency distribution of VDR fok I gene among TB patients and compared it with controls in a Ugandan population.
After obtaining permission and informed consent (see ethical considerations) a pilot study in newly diagnosed smear positive TB patients and healthy controls was conducted between the months of April and June 2013 at the Mulago National Referral and Teaching Hospital located in North of Kampala, Uganda. By consecutive sampling adult patients who presented with persistent cough for more than 3 weeks had a positive Ziehl-Nielsen smear test for the first time and signed a written consent were considered eligible. Health workers and medical trainees who worked at the TB out patient’s clinic and other wards exposed to TB were matched for sex and age with the control group. All subjects were screened for human immunodeficiency virus (HIV) and data and consent forms were filled. Blood was then collected into tubes with ethylenediaminetetraacetic acid (EDTA) anticoagulant (Becton, Dickinson and company New Jersey USA) and whole blood samples were stored at 2–8°C before DNA column extraction.
Genotyping of VDR fok I gene was performed by PCR–direct sequencing method on ABI 310 sequencer. Human genomic DNA was extracted using the Genotype DNA isolation kit version 2.0 (HAIN Life Science, Germany) according to the manufacturer’s instructions. To check the purity genomic DNA was run on a 1.0% agarose gel in a 1 × Tris Acetate EDTA (TAE) buffer containing ethidium bromide for 1 hour at 120 volts and visualized under an UV transilluminator. The primer sequences used in our study were forward A 5’-GC TGG CCC TGG CAC TGA CTC TG TCT -3’ and reverse 5’-ATG GAA ACA CCT TGC TTC TTC TCC CTC as -3’ described by Harris et al. (1997). PCR was performed in 15.1 µl reaction volumes which contained 7 µl of PCR water, 1 µl of 25 mM Mgcl2, 1 µl of 10X master mix (Fisher biotec company Australia), 1 µl of 0.22 mg of each of the forward and reverse primers (Integrated DNA technologies company USA), 0.1 µl of 5U of Pre heated Taq polymerase and finally 5 µl of pure DNA were added and the reaction was thoroughly mixed. The PCR GTQ Thermocycler programme involved an initial denaturation at 95°C for 5 minutes, followed by 30 cycles 70°C annealing for 45 seconds, elongation at 72°C for 1 minute and a final elongation at 72°C, for 10 minutes. Amplified DNA was purified using an extraction kit (JET quick company, USA) according to manufacturer’s instructions. Cycle sequencing PCR was performed using the Big Dye terminator KIT version 3, in the Thermocycler Gene Amp PCR system 9700. The cycling programme involved denaturation at 96°C for 1 minute for 1 cycle ,annealing at 96°C for 10 minutes 30 cycles, elongation at 70°C for 30 seconds and a final elongation at 60°C for 5 minutes. The Dye EX 2.0 spin kit (250) (QIAGEN, Germany) was used to remove dye and other PCR products following the manufacturer’s instructions. Sequencing of the PCR products was done using the ABI Big Dye Termination kit (Applied Biosystems, USA) and the ABI prism 310 Genetic analyzer (Applied Biosystems). Sequences obtained were compared to those in the Gene Bank database by applying BLAST; NCBI tool. Finally DNA baser sequence assembly software version 4.7.0 bioinformatics tool was used for the Identification of fok I gene polymorphism. Polymorphism of FF genotype was identified on sequences with a start codon of ACG instead of an ATG followed by an ATG downstream. The homozygous ff had a start codon of ATG and consequently another ATG in the sequence. The heterozygous Ff genotype had an ACG, ACG and an ATG.
The sample size was calculated based on data of a previous study9 with mean of 78.3 for TB and 83.5 for healthy controls and a total of 82 samples at a power of 80% would detect an effect. We therefore divided the samples into cases and controls for a pilot study. Data were summarized into odds ratios, 95% confidence interval (95% CI), and alpha of p<0.05 was considered significant using STATA software (Stata Corp. STATA 12.0, College Station, Texas, USA). The frequency of vitamin D fok I gene was determined using percentages. A Chi square test was used to determine the allele and genotype frequency distribution and potential deviation from Hardy Weinberg equilibrium. Fok I genotypes were categorized as FF for dominant, Ff heterozygous and ff for homozygous variant.
A total of 82 participants were enrolled in the study. Of these, 41 had pulmonary TB while 41 were healthy subjects with no history of the disease. Majority of participants were males (66.0%). Among the participants with TB 73.2% were males while the females were 26.8%. The percentage of HIV/TB patients in this study of 24.4% was lower than the non HIV/TB patients. Thirteen (15.9%) of the 82 participants were HIV seropositive, of these three were from the control group. Most of the TB patients were drivers, farmers, teachers and others. Other lifestyle factors did not show significant variation. Details of the above results are shown in Table 1. Only four (10%) participants were aware of exposure to TB (data not shown).
Following DNA extraction and amplification of the fok I gene in all 82 samples the electrophoresis results showed a band between 200 and 300 bps (Figure 1). The PCR products were then sequenced and a BLAST on the query sequences was done against those in the NCBI Gene bank to confirm the sequences of fok I. The gene product was named vitamin D3 receptor isoform VDRB1 vitamin D3 receptor isoform VDRA. It showed sequences between 238 bp and 260 bp with a percentage identity between 98–99% (see Dataset 1).
The frequencies of Fok I genotypes, FF, Ff, and ff in TB patients and in the controls were assessed. Table 2 shows the genotype and allele frequency distribution of this study population with their odds ratio, confidence intervals and p values. No significant difference was observed in the distribution of the homozygous genotypes in the study population and the two alleles. FF genotype was used as the reference and other genotypes were compared while the f allele was used as reference to the F allele (Table 2). The heterozygous Ff genotype was associated with TB, (OR 3.2) since it occurred more in the TB patients than the healthy subjects as shown in Table 2. There was a predominance in the allele distribution of F dominant verse the f recessive allele in both the TB and control group. Both HIV and non HIV individuals predominantly had the FF genotype with only one HIV/TB patient having the recessive ff genotype. Hardy Weinberg equilibrium was tested to establish if the genotype and allele distribution among the TB patients and the control population was in equilibrium. The Chi square test estimated a deviation coefficient of 0.041 in the TB patients and 0.045 in the healthy subjects.
An analysis of fok I gene was performed to determine the frequency distribution of the fok I gene polymorphism among male and female subjects (Table 3). The frequency distribution of fok I FF homozygous dominant genotype was male was 57.3% for male and 32.9 % for females. There was one female with heterozygous Ff genotype while homozygous recessive ff genotype was only found male subjects. These results showed a significant difference in the distribution of the fok I genotypes among male and female subjects given a p value of 0.02. When a multivariate analysis was further performed against fok I, gender and other variables a p 0.01 was obtained.
Fok I Genotypes | Male, n, (%) | Female, n, (%) | P value |
---|---|---|---|
FF
Ff ff | 47 (57.3%) 3 (3.7%) 4 (4.9%) | 27 (32.9%) 1 (1.2 %) 0 (0.0%) |
0.02 |
VDR polymorphisms have been studied in different populations and have shown various results. We document for the first time the frequency distribution of VDR polymorphism in the Ugandan population. This study reports a frequency of 87.8% of FF genotype in the Pulmonary Tuberculosis (PTB) patients; however previous studies have reported a 65% frequency of FF genotype among Indian PTB patients, 72% among the Venda of South Africa and 62% among west-Africans. The low frequency of f allele in this population is consistent with another study10 showing that the f allele of Fok I occurs less frequently in Africans compared to Caucasians and Asians, (Africans 24%, Caucasians 34%, and Asians 51%). The ff genotype frequency of 4.8% was similar to the one described for PTB patients in the African case control study of Gambia Guinea and Guinea Bissau. The ff genotype was reported to be 5% in an Indian healthy population11, while it was 6% in the African American population and 3% in the Venda of South Africa12. The ff genotype was only found among the males in this study. This observation differs from an Indian study showing the FF genotype only in males and a significant difference in the genotypic distribution between males and females. The VDR polymorphism was not in equilibrium with Hardy Weinberg and this could be attributed to the low frequency of the ff genotype. Among other studies that have also reported disequilibrium in the VDR are the studies in the three West African countries6 and a study on the Paraguay population8. Disequilibrium observed in these studies may be due to genotyping errors but was overcome by use of high quality DNA which was achieved by secondary PCR. Other factors are genetic drift and mutations. There was no significant difference in the distribution of fok I genotype in different ethnic groups of our study population, (p value 0.30). However a difference was observed when the analysis was based on gender. The low numbers of females in this study may probably be explained by the female sex being a protection against TB according to a study13. Among the 5 % to 10 % of individuals who develop active disease about 70% are males.
Although the frequency distribution of the homozygous fok I genotypes was not significantly different in the Ugandan TB patients and controls the Ff heterozygous genotype appears to be associated with TB in the Ugandan population. Therefore further studies are needed to elucidate this relationship.
F1000Research: Dataset 1. Human VDR fok I gene nucleotide sequences, 10.5256/f1000research.9109.d13042014
Written informed consent to participate in the study and publish these data was obtained from all participants.
Permission to carry out the study was obtained from the Research and Ethics Committee Mulago Hospital, (ref MREC: 329), the Institutional Review Board of the School of Biomedical Sciences Higher Degrees Research and Ethics committee, (ref SBS 108), and from the Uganda National Council of Science and Technology (ref No.1431). Before investigations, informed consent was obtained from the study subjects. Patients’ personal information was kept confidential by using serial codes instead of names on the questionnaire.
AE: conceived the study, AE and EJ: designed the experiments, AE, WW and EJ: developed proposal and was involved in research, MP contributed to study design and other statistics expertise AE WW, MP, KA and EJ prepared the first draft of manuscript. All authors were involved in the revision of the draft manuscript and have agreed to the final content.
This research was sponsored by the Makerere University and the Swedish International Development Cooperation Agency Bilateral under the Gender Main streaming Directorate. Support for manuscript writing was provided by Fogarty International Center, National Institutes of Health (grant # D43TW009771 “HIV co-infections in Uganda: TB, Cryptococcus, and Viral Hepatitis.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We would like to acknowledge the study participants, Mulago hospital and MBN laboratories for analysis.
We thank the Swedish International Development Cooperation Agency through Gender Main streaming, Makerere University for funding this research.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Competing Interests: No competing interests were disclosed.
Competing Interests: No competing interests were disclosed.
References
1. Wigginton JE, Cutler DJ, Abecasis GR: A note on exact tests of Hardy-Weinberg equilibrium.Am J Hum Genet. 2005; 76 (5): 887-93 PubMed Abstract | Publisher Full TextCompeting Interests: No competing interests were disclosed.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | |||
---|---|---|---|
1 | 2 | 3 | |
Version 1 29 Jul 16 |
read | read | read |
Click here to access the data.
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)