The frequency distribution of vitamin D Receptor fok I gene polymorphism among Ugandan pulmonary TB patients

Background: Mycobacterium tuberculosis (TB) is still a major problem globally and especially in Africa. Vitamin D deficiency has been linked to TB in the past and studies have found vitamin D deficiency to be common among Ugandan TB patients. The functional activity of vitamin D is dependent on the genotype of the vitamin D receptor (VDR) polymorphic genes. Recent findings have indicated that VDR polymorphisms may cause increased resistance or susceptibility to TB. The vitamin D ligand and its receptor play a pivotal role in innate immunity by eliciting antimicrobial activity, which is important in prevention of TB. The fok I vitamin D receptor gene has extensively been examined in TB patients but findings so far have been inconclusive. Objectives: This study sought to investigate the frequency distribution of the VDR fok I gene polymorphisms in pulmonary TB patients and controls. Methods: A pilot case control study of 41 newly diagnosed TB patients and 41 healthy workers was set up. Vitamin D receptor fok I gene was genotyped. Results: The frequency distribution of fok I genotype in Ugandan TB patients was 87.8% homozygous-dominant (FF), 7.3% (Ff) heterozygous and 4.8% (ff) homozygous recessive. For normal healthy subjects the frequencies were (FF) 92.6%, (Ff) 2.4% and (ff) 4.8%. No significant difference was observed in the FF and ff genotypes among TB patients and controls. The Ff heterozygous genotype distribution appeared more in TB patients than in controls. A significant difference was observed in the fok I genotype among gender p value 0.02. No significant difference was observed in ethnicity, p value 0.30. Conclusions: The heterozygous Ff fok I genotype may be associated with TB in the Ugandan population.


Abstract
(TB) is still a major problem globally Background: Mycobacterium tuberculosis and especially in Africa. Vitamin D deficiency has been linked to TB in the past and studies have found vitamin D deficiency to be common among Ugandan TB patients. The functional activity of vitamin D is dependent on the genotype of the vitamin D receptor (VDR) polymorphic genes. Recent findings have indicated that VDR polymorphisms may cause increased resistance or susceptibility to TB. The vitamin D ligand and its receptor play a pivotal role in innate immunity by eliciting antimicrobial activity, which is important in prevention of TB. The vitamin D receptor gene has extensively been fok I examined in TB patients but findings so far have been inconclusive.
: This study sought to investigate the frequency distribution of the Objectives VDR gene polymorphisms in pulmonary TB patients and controls.

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 countries 1 and TB remains a major public health problem. Although it is still unclear how individuals develop active TB 2 there is emerging evidence of multi drug resistant TB (MDRTB). Genetic predisposition to TB has been suggested in several studies 3,4 but the immunological associations with genetic polymorphisms are still unclear 5 . 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 inconsistent 6 . Various VDR polymorphisms have been identified and associated with TB susceptibility or resistance 2,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 respectively 7,11 . This study investigated the frequency distribution of VDR fok I gene among TB patients and compared it with controls in a Ugandan population.

Methods
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 fok I gene 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 Mgcl 2 , 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.

Statistical analysis
The sample size was calculated based on data of a previous study 9 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.

Results
Description of social demographic characteristics of participants 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).
Detection of fok I gene PCR amplicons and BLAST identification of fok I nucleotides 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 D 3 receptor isoform VDRA. It showed sequences between 238 bp and 260 bp with a percentage identity between 98-99% (see Dataset 1).

Frequency distribution of fok I genotypes and alleles among TB patients and controls
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.
Genotype fok I distribution by gender An analysis of fok I gene was performed to determine the frequency distribution of the fok I gene polymorphism among male and female subjects (  10 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 population 11 , while it was 6% in the African American population and 3% in the Venda of South Africa 12 . 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 countries 6 and a study on the Paraguay population 8 . 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 study 13 . Among the 5 % to 10 % of individuals who develop active disease about 70% are males.

Conclusion
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.

Consent
Written informed consent to participate in the study and publish these data was obtained from all participants.

Ethical considerations
Permission to carry out the study was obtained from the Research The vitamin D receptor plays an important role in the immune system and the roles have widely been studied. Several studies have been conducted to assess the frequency distribution as well as roles of Vitamin D's in various diseases conditions; including; multiple sclerosis, cancer, diabetes and TB. However, the findings have been varying and inconclusive, imposing the need of carrying further investigations in different populations and settings to determine the roles and frequency distrition. It was on this varying state of Vitamin D that the authors, Ester , conducted a study in Uganda aiming to et al investigate the frequency distribution of the VDR gene ( I) polymorphisms in pulmonary TB patients fok compared to controls in Ugandan population.
The article though short and with relatively small sample size of 82 (TB patients and health controls) is well written, appropriate description is given covering all sections; introduction, objectives methods, results, discussions and conclusion. The study found no significance difference in the frequency distribution of homozygous -dominant (FF) and homozygous receive (ff) genotype among TB patients and controls. However, heterozygous (Ff) genotype was more frequent among TB patients compared to controls (7.3% vs 2.4%). The author concluded that Ff gene may be associated with TB because it occurred more in TB patients than in health controls.

Specific comment:
In the discussion section the author gave references to other studies done elsewhere on the frequency of FF and ff genotypes, however, references on the frequency of heterozygous (Ff) is missing and should be added in the discussion section. In this study the difference was observed on this gene for the Ugandan studied population hence reference to other studies findings is important.

I have read this submission. I believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.
No competing interests were disclosed.