Int J Biol Sci 2012; 8(3):375-382. doi:10.7150/ijbs.4047
Polymorphic Allele of Human MRC1 Confer Protection against Tuberculosis in a Chinese Population
1. Institute of Cell Biology, Zhejiang University, Hangzhou 310058, P.R.China;
This is an open access article distributed under the terms of the Creative Commons Attribution (CC BY-NC) License. See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Zhang X, Jiang F, Wei L, Li F, Liu J, Wang C, Zhao M, Jiang T, Xu D, Fan D, Sun X, Li JC. Polymorphic Allele of Human MRC1 Confer Protection against Tuberculosis in a Chinese Population. Int J Biol Sci 2012; 8(3):375-382. doi:10.7150/ijbs.4047. Available from http://www.ijbs.com/v08p0375.htm
Mannose receptor is a member of the C-type lectin receptor family involved in pathogen molecular-pattern recognition, and plays a critical role in shaping host immune response. Single nucleotide polymorphisms (SNPs) in the MRC1 gene may affect expression levels and differences in the structure and function of proteins in different individuals, thereby affecting individual susceptibility to pulmonary tuberculosis. However, to date, MRC1 polymorphisms associated with susceptibility to pulmonary tuberculosis have not yet been reported. The present study aimed to investigate potential associations of SNPs in the MRC1 gene with pulmonary tuberculosis in a Chinese population. Six SNPs (G1186A, G1195A, T1212C, C1221G, C1303T and C1323T) in exon 7 of the MRC1 gene were genotyped using the PCR and DNA sequencing methods in the pulmonary tuberculosis patients and the healthy controls. Linkage disequilibrium analysis was performed between polymorphic sites. The study found that the allele frequency of G1186A (rs34039386) of the MRC1 gene in a Chinese population was higher in the pulmonary tuberculosis group than the healthy control group. There was a significant difference in frequency distribution between the two groups (P = 0.037; OR = 0.76; 95% CI, 0.58-0.98). Genotypic analysis also indicated that the AG genotypes in a Chinese population were significantly correlated with pulmonary tuberculosis (P < 0.01; OR = 0.57; 95% CI, 0.37-0.87). After adjustment for age and gender, G1186A sites were found to be dominant (P < 0.01; OR = 0.59; 95% CI, 0.40-0.87), over-dominant (P = 0.045; OR = 0.69; 95% CI, 0.47-0.99) and additive models (P = 0.041; OR = 0.76; 95% CI, 0.59-0.99) in association with pulmonary tuberculosis. But, no association was found between the other 5 SNPs (G1195A, T1212C, C1221G, C1303T and C1323T) and tuberculosis (P > 0.05). This study is the first to report that genetic variants in the MRC1 gene can be associated with pulmonary tuberculosis in a Chinese population, and may reduce the risk of infecting pulmonary tuberculosis. This also provides a new experimental basis to clarify the pathogenesis of pulmonary tuberculosis.
Keywords: Mannose receptor, MRC1 gene, Tuberculosis, Single-nucleotide polymorphism, Chinese.
Tuberculosis (TB) is one of the world's most serious public health threats. Over the past several years, China is most profoundly affected by TB. Each year, the new cases of TB account for 18% of the world's population . Although about one third of the world's population is thought to be infected with Mycobacterium tuberculosis (MTB), only 5-15% of people develop clinically active TB during their lifetime . Some evidence suggests that certain genetic factors may be involved in innate immunity and play important roles in susceptibility to TB at the individual level. Genetic studies showed that both genes and environmental factors are associated with the pathogenesis of TB [3, 4]. Therefore, it is important to identify genes that mediate susceptibility to TB.
TB is an infectious disease caused by MTB, which mainly lives on the monocyte/macrophage system. Cellular immunity is involved in resistance to infectious disease caused by MTB. Activated macrophages perceive the invasion of MTB and lead to active and passive immune response, such as, antigen presenting, T cell activation, B cell activation, the production of interleukin (IL), interferon (IFN) and transformation growth factor (TGF) . So, monocyte/macrophage system plays a key role in the early identification of MTB and the incidence of pulmonary TB . Pattern recognition receptors (PRRs) are located on the surface of the macrophages and dendritic cells, which belong to the body's natural immune system and are the core of the receptor molecules identifying the pathogen. Many classes of PRRs have been described, including Toll-like receptors (TLR), NOD-like receptors (NLR) and C-type lectin receptors (CLR) [6-8]. Recently, polymorphisms in the TLR [3, 9-18] and NLR [19, 20] genes have been shown to be associated with susceptibility to pulmonary TB. The mannose receptor (MR) belongs to CLR, and the predisposition of MRC1 gene variants to pulmonary TB have not been reported yet.
MR is a member of the CLR family, which plays an important role in innate immunity . MR is predominantly present on alveolar macrophages and dendritic cells and recognizes glycan structures containing mannose, fucose and N-acetylglucoasmine, which are commonly found on the cell walls of pathogenic micro-organism such as mycobacteria, fungus, parasites, and yeast [6, 22, 23]. MR binds to mannose-capped lipoarabinomannans (ManLAM), a cell wall component of MTB [24, 25], loaded to the antigen-presenting cells (APC), and then presented to T cells which play a role in immune response . MR can also help macrophages to phagocytize MTB [27, 28] and plays an important role in innate immunity [25, 26]. The MRC1 gene, encoding the human MR, is located on chromosome 10p12 and consist of 30 exons. Several reports have shown that the MRC1 gene is associated with the susceptibility to some diseases, such as asthma, sarcoidosis and leprosy [29-31]. We speculated that the MRC1 gene may be an excellent candidate for susceptibility to TB, and tested for associations between MRC1 polymorphisms and the pulmonary TB in our Chinese case-control analysis.
Materials and Methods
Patients and Control Subjects
A total of 222 Chinese subjects with pulmonary TB, aged 18-66 years (mean age 37.5 ± 8.6 years) were recruited from the Sixth Hospital of Shaoxing and Hangzhou Red Cross Hospital. The diagnosis of pulmonary TB was made according to established criteria , 98 cases were confirmed by the M. tuberculosis cultures and 124 cases were diagnosed by chest X-ray and computed tomography (CT) revealing evidence of typical active tuberculosis. A total of 232 healthy controls, aged 20-64 years (mean age 39.6 ± 9.4 years) were recruited from the Zhejiang Hospital between October 2008 and July 2011 (Table 1).
The early morning, fasting blood samples from the patients and controls were collected in 3.0 mL EDTA tubes and stored at -70 °C. For each gene polymorphisms, it was ensured that the final selected samples of TB patients were similar to controls in the ratio of gender, age, presence of TB history of relatives, BCG vaccination and Tuberculin skin test (>10 mm). This study was approved by the Ethics Committee of the Faculty of Medicine (Zhejiang University, China), and informed consents were obtained from all subjects before blood sampling.
We acquired the MRC1 gene sequence by using the NCBI Gene database (http://www.ncbi.nlm.nih.gov/gene) and the Ensemble Genome Browser Database (http://asia.ensembl.org). We also obtained the information of single nucleotide polymorphisms (SNPs) in this gene by using dbSNP database (http://www.ncbi.nlm.nih.gov/snp/), and screened the frequency of SNPs by using HapMap database (http://www.hapmap.org).
Finally, we selected six SNPs in exon 7 of the MRC1 gene, G1186A (rs34039386, Gly396Ser), G1195A (rs71497223, Gly396Ser), T1212C (rs71497224, Ile404Ile), C1221G (rs34284571, Leu407Phe), C1303T (Leu435 Phe) and C1323T (Asn441Asn).
Characteristics of healthy controls and TB patients.
Data are reported as number with percent in parentheses, unless otherwise stated. TB: tuberculosis; ND: not determined. aDifference between total patients and controls (t test). bDifference between total patients and controls (χ2 test).
SNPs in the MRC1 gene were analyzed by polymerase chain reaction (PCR) and direct sequencing. Genomic DNAs were extracted from the peripheral blood leukocytes with DNA extraction kit (QIAamp® DNA Blood Min Kit, Germany) according to the manufacturer's instruction.
The PCR primer included forward: TTGAGGCTGCAATGAGACAT and reverse: AGTGTAAGGTAGACTGCTCT. The PCR was performed by denaturing at 95 °C for 5 min, followed by 35 cycles at 95 °C for 50 s, 63 °C for 50 s, and 72 °C for 30 s, and a final extension at 72 °C for 10 min. The amplified products were purified with PCR cleaning kit (AxyPrep PCR cleaning kit, USA) and then identified by scanning with the ABI 3100 sequencer (Applied Biosystems, Carlsbad, CA, USA).
The chi-square test was used to compare allele and genotype distribution in the TB patients and control subjects by the two-tailed Fisher exact method with the GraphPad Prism version 5.01 software. The odd ratio (OR) and 95% confidence interval (CI) were calculated by Miettinen's method. The value of OR is used to estimate the relative risk. If the value of OR is equal to 1, it indicates that the factor has no effect on the incidence of disease. If the value of OR is greater than 1, it indicates that the factor is a risk to the incidence of disease. If the value of OR is less than 1, it indicates that the factor is protective.
Pairwise linkage disequilibrium was calculated with SNPstas (http://bioinfo.iconcologia.net/SNPstats), which was estimated by calculating pairwise based on D′(D′value is the relative size of the LD linkage disequilibrium, which is on behalf of the differences between the expected haplotype frequency and the observed haplotype frequency and scaled in [−1,1] range) and R (R is the correlation coefficient between alleles) that are on behalf of the degree of linkage disequilibrium between each pair of SNPs . Values of D′> 0.75 were considered to be strong pairwise linkage disequilibrium and we construct haplotypes between each pair of SNPs. Haplotype frequencies and associations were calculated with expectation-maximization algorithm. Using logistic regression method, we calculated the allele combinations of OR values and 95% CI referring to the highest frequency of haplotype allele and the additive genetic model. P value was calculated to estimate haplotype association with pulmonary TB.
The MRC1 gene was amplified and sequenced by PCR in all samples. Six polymorphisms were found in this region. All MRC1 SNPs were in Hardy-Weinberg equilibrium in the control group and the pulmonary TB group (P > 0.05) (Table 2). The frequencies of the G/G, A/G and A/A genotypes for the G1186A site in the control group were 0.267, 0.517 and 0.216, respectively, while, the frequencies of the three genotypes in the pulmonary TB group were 0.383, 0.423 and 0.194, respectively. The frequency of A/G genotype was lower in the pulmonary TB group than the control group and there was a significant difference between the two groups (P < 0.01; OR = 0.57; 95% CI, 0.37 - 0.87), but the frequency of A/A genotype was not significantly different between the two groups (P > 0.05). No significant differences between patients and controls were observed in the genotype distribution of G1195A, T1212C, C1221G, C1303T and C1323T (P > 0.05) (Table 3).
The allele frequency of G for the G1186A was higher in the pulmonary TB group than the control group and there was a significant difference between the two groups (P = 0.037; OR = 0.76; 95% CI, 0.58-0.98) (Table 2). The allele frequencies of G1195A, T1212C, C1221G, C1303T and C1323T were not significantly different between patients and controls (P > 0.05) (Table 2).
Distribution of the MRC1 SNP allele frequencies in the TB group (n=222) and the control group (n=232) in a Chinese population.
SNP: single nucleotide polymorphism; HWE: Hardy-Weinberg equilibrium; OR: odd ratio; CI: confidence interval; N: number of alleles; Freq: frequency. P value and OR were obtained by Chi-square test. Two-tailed Fisher's exact test was used because of small sample numbers. All of the MRC1 SNPs were in Hardy-Weinberg equilibrium.
Distribution of the MRC1 SNP genotype frequencies in the TB group (n=222) and the control group (n=232) in a Chinese population.
SNP: single nucleotide polymorphism; OR: odds ratios; CI: confidence intervals; N: number of alleles; Freq: frequency. P-value and OR were obtained by Chi-square test. Two-tailed Fisher's exact test was used because of small sample numbers.
We performed logistic repression analysis with all SNPs using a recessive model adjusted for age and sex and observed that G1186A SNP was associated with the pulmonary TB in co-dominant, dominant (P < 0.01; OR = 0.59; 95% CI, 0.40-0 .87), over-dominant (P = 0.045; OR = 0.69; 95% CI, 0.47-0 .99) and additive model (P = 0.041; OR = 0.76; 95% CI ,0.59-0 .99) (Table 4, Table 5). According to the minimum Akaike information criterion (AIC), this genetic model was the dominant mode. However, the genetic model of other SNPs did not correlate with the pulmonary TB (P > 0.05).
Pairwise LD between the 6 SNPs of the MRC1 gene was calculated for the cases and controls in a Chinese population (Figure 1). We found strong LD (D′> 0.75, P < 0.0001) between some pairs of markers in the MRC1 gene including G1186A and C1221G, G1186A and C1323T, G1195A and T1212C, G1195A and G1221C, G1195A and C1303T, T1212C and C1221G, T1212C and C1303T, C1221G and C1303T. Therefore, based on logistic-regression analyses, we observed that the haplotype of each SNP pair was associated with pulmonary TB (Table 6). We found that the frequency of the haplotype AGTCCT in the pulmonary TB group was higher than the control group, whereas, the frequency of the haplotype GGTCCT, GGCCCT, GGTCCC in the pulmonary TB group was lower than the control group. All haplotypes with the MRC1 gene polymorphisms showed no significant association with pulmonary TB (Table 6).
Linkage disequilibrium analysis of single-nucleotide polymorphisms (SNPs) in the MRC1 gene in a Chinese population. The number at the intersection of each pair of SNPs represents the pairwise D, D′, r2 and P values between two SNPs. Values of D′ > 0.75 were considered to be strong pairwise linkage disequilibrium.(Click on the image to enlarge.)
Association of MRC1 polymorphisms with TB in a Chinese population with logistic regression.
aAdjusted for age and sex. Freq: frequency; OR: odd ratio; CI: confidence interval.
The association of G1186A with TB in a Chinese population using logistic regression
aAdjusted for age and sex. OR: odd ratio; 95% CI: 95% confidence interval; AIC: Akaike's Information Criterion; BIC: Bayesian Information Criterion.
Haplotype frequencies of polymorphisms variants of the MRC1 gene in patients with TB (n=222) and healthy controls (n=232) in a Chinese population
aAdjusted for age and sex. Freq: frequency of haplotype; OR: odds ratio; CI: confidence interval.
Recent studies have reported that MR plays an important role in innate and adaptive immunity . These receptors activate immune response by extracellular carbohydrate recognition that binds the glycan structures of microbial pathogens, and may play an important role in TB infection . However, polymorphisms in the MRC1 gene, encoding MR, have not been reported to be associated with pulmonary TB. We postulated that MRC1 polymorphisms may be associated with the development of pulmonary TB. Six SNPs in exon 7 of the MRC1 gene were analyzed. The results showed that both the frequency of genotype (P = 0.037; OR = 0.76; 95% CI, 0.58-0.98) and allele (P < 0.01; OR = 0.57; 95% CI, 0.37-0.87) for the G1186A (rs34039386) site were significantly different between the pulmonary TB group and the control group. The OR values of genotype and allele for G1186A were less than 1, indicating that MRC1 polymorphisms have effects on the development of pulmonary TB, and may reduce the risk of being infected with pulmonary TB.
The G1186A site of MRC1 is non-synonymous mutation that changes Gly to Ser. The mutation of this site may change the function of MR, which may also affect the binding of MR protein to polysaccharide structure on the MTB surface, thus, contributing to the secretion of anti-inflammatory factors. Meanwhile, we also analyzed association between the established haplotype of MRC1 gene and the incidence of pulmonary TB. The results showed that haplotype AGTCCT, GGTCCT, GGCCCT and GGTCCC in the MRC1 gene had no significant association with pulmonary TB (P > 0.05). In addition, the other 5 polymorphism sites of MRC1 gene [G1195A (rs71497223, Gly396Ser), T1212C (rs71497224, Ile404Ile), C1221G (rs34284571, Leu407Phe), C1303T (Leu435Phe) and C1323T (Asn441Asn)] also showed no significant association with pulmonary TB (P > 0.05), suggesting that these polymorphism sites were not associated with susceptibility to pulmonary TB.
Andrea et al.  found that the allele frequency of G1186A site in the exon 7 of MRC1 gene had significant disparity between leprosy per se cases and healthy controls in the Vietnamese population (P = 0.036, OR = 0.76; 95% CI, 0.60-0.96). The OR value was less than 1, and the results showed that this site has significant protective role for leprosy and may reduce the risk of leprosy infection. Similarly, the allele frequency of G1186A site was also found to be significantly different between leprosy per se cases and healthy controls (P = 0.016, OR = 1.34; 95% CI, 1.06-1.70) in the Brazilian population . The OR value was more than 1, suggesting that this polymorphism is associated with susceptibility to leprosy. Hattori et al.  reported that the allele frequency of G1186A site in the exon 7 of MRC1 gene had no significant association with asthma in the Japanese population. We proposed that the individual differences caused by genetic factors, may be the key factor in the difference of disease susceptibility and the gene polymorphisms . In the present study, we analyzed alleles, genotypes and haplotypes, and concluded that the G1186A polymorphism of MRC1 gene was associated with TB, and may have a protective function.
This work was supported by grants from Zhejiang Province Special Sci-Tech Projects (No.2009C03011-3), National Natural Science Foundation of China (No.81072724), and National Special Sci-Tech Projects (No.2012ZX10005001-006; No.2008ZX10005-010).
Conflict of Interests
The authors have declared that no conflict of interest exists.
1. Vashishtha VM. WHO Global Tuberculosis Control Report 2009: Tuberculosis elimination is a distant dream. Indian Pediatr. 2009;46:401-2
2. Rosman MD O-EA. Clinical presentation and treatment of tuberculosis. In: (ed.) Fishman AP. Fishman's Pulmonary Diseases and Disorders. New York, USA: McGraw-Hill. 1998:2483-502
3. Xue Y, Jin L, Li AZ. et al. Microsatellite polymorphisms in intron 2 of the toll-like receptor 2 gene and their association with susceptibility to pulmonary tuberculosis in Han Chinese. Clin Chem Lab Med. 2010;48:785-9
4. Delgado JC, Baena A, Thim S. et al. Ethnic-specific genetic associations with pulmonary tuberculosis. J Infect Dis. 2002;186:1463-8
5. Tsolaki AG. Innate immune recognition in tuberculosis infection. Adv Exp Med Biol. 2009;653:185-97
6. Brown GD. Dectin-1: a signalling non-TLR pattern-recognition receptor. Nat Rev Immunol. 2006;6:33-43
7. Jo EK. Mycobacterial interaction with innate receptors: TLRs, C-type lectins, and NLRs. Curr Opin Infect Dis. 2008;21:279-86
8. Hill AV. Aspects of genetic susceptibility to human infectious diseases. Annu Rev Genet. 2006;40:469-86
9. Ma X, Liu Y, Gowen BB. et al. Full-exon resequencing reveals toll-like receptor variants contribute to human susceptibility to tuberculosis disease. PLoS One. 2007;2:e1318
10. Jin L, Zhu XF, Zhong JP. et al. The Arg753Gln Polymorphism of the Human Toil-like Receptor 2 Gene and Its Association with Tuberculosis Disease in Zhejiang Han Population Chinese. Chinese Journal of Cell Biology. 2007;29:229-31
11. Newport MJ, Allen A, Awomoyi AA. et al. The toll-like receptor 4 Asp299Gly variant: no influence on LPS responsiveness or susceptibility to pulmonary tuberculosis in The Gambia. Tuberculosis (Edinb). 2004;84:347-52
12. Davila S, Hibberd ML, Hari Dass R. et al. Genetic association and expression studies indicate a role of toll-like receptor 8 in pulmonary tuberculosis. PLoS Genet. 2008;4:e1000218
13. Schroder NW, Schumann RR. Single nucleotide polymorphisms of Toll-like receptors and susceptibility to infectious disease. Lancet Infect Dis. 2005;5:156-64
14. Thuong NT, Hawn TR, Thwaites GE. et al. A polymorphism in human TLR2 is associated with increased susceptibility to tuberculous meningitis. Genes Immun. 2007;8:422-8
15. Zhang YX, Xue Y, Liu JY. et al. Association of TIRAP (MAL) gene polymorhisms with susceptibility to tuberculosis in a Chinese population. Genet Mol Res. 2011;10:7-15
16. Xue Y, Zhao ZQ, Hong D. et al. Lack of association between MD-2 promoter gene variants and tuberculosis. Genet Mol Res. 2010;9:1584-90
17. Xue Y, Zhao ZQ, Wang HJ. et al. Toll-like receptors 2 and 4 gene polymorphisms in a southeastern Chinese population with tuberculosis. Int J Immunogenet. 2010;37:135-8
18. Ogus AC, Yoldas B, Ozdemir T. et al. The Arg753GLn polymorphism of the human toll-like receptor 2 gene in tuberculosis disease. Eur Respir J. 2004;23:219-23
19. Proell M, Riedl SJ, Fritz JH. et al. The Nod-like receptor (NLR) family: a tale of similarities and differences. PLoS One. 2008;3:e2119
20. Austin CM, Ma X, Graviss EA. Common nonsynonymous polymorphisms in the NOD2 gene are associated with resistance or susceptibility to tuberculosis disease in African Americans. J Infect Dis. 2008;197:1713-6
21. Torrelles JB, Azad AK, Schlesinger LS. Fine discrimination in the recognition of individual species of phosphatidyl-myo-inositol mannosides from Mycobacterium tuberculosis by C-type lectin pattern recognition receptors. J Immunol. 2006;177:1805-16
22. East L, Isacke CM. The mannose receptor family. Biochim Biophys Acta. 2002;1572:364-86
23. McGreal EP, Miller JL, Gordon S. Ligand recognition by antigen-presenting cell C-type lectin receptors. Curr Opin Immunol. 2005;17:18-24
24. Beharka AA, Crowther JE, McCormack FX. et al. Pulmonary surfactant protein A activates a phosphatidylinositol 3-kinase/calcium signal transduction pathway in human macrophages: participation in the up-regulation of mannose receptor activity. J Immunol. 2005;175:2227-36
25. Schlesinger LS, Kaufman TM, Iyer S. et al. Differences in mannose receptor-mediated uptake of lipoarabinomannan from virulent and attenuated strains of Mycobacterium tuberculosis by human macrophages. J Immunol. 1996;157:4568-75
26. Kang PB, Azad AK, Torrelles JB. et al. The human macrophage mannose receptor directs Mycobacterium tuberculosis lipoarabinomannan-mediated phagosome biogenesis. J Exp Med. 2005;202:987-99
27. Schlesinger LS. Macrophage phagocytosis of virulent but not attenuated strains of Mycobacterium tuberculosis is mediated by mannose receptors in addition to complement receptors. J Immunol. 1993;150:2920-30
28. Nigou J, Zelle-Rieser C, Gilleron M. et al. Mannosylated lipoarabinomannans inhibit IL-12 production by human dendritic cells: evidence for a negative signal delivered through the mannose receptor. J Immunol. 2001;166:7477-85
29. Alter A, de Leseleuc L, Van Thuc N. et al. Genetic and functional analysis of common MRC1 exon 7 polymorphisms in leprosy susceptibility. Hum Genet. 2010;127:337-48
30. Hattori T, Konno S, Hizawa N. et al. Genetic variants in the mannose receptor gene (MRC1) are associated with asthma in two independent populations. Immunogenetics. 2009;61:731-8
31. Hattori T, Konno S, Takahashi A. et al. Genetic variants in mannose receptor gene (MRC1) confer susceptibility to increased risk of sarcoidosis. BMC Med Genet. 2010;11:151
32. Diagnostic criteria for pulmonary tuberculosis. Ministry of Health of China. http://www.moh.gov.cn/publicfiles/business/cmsresources/zwgkzt/cmsrsdocument/doc3242.pdf
33. Liu CP, Li XG, Lou JT. et al. Association analysis of the PHOX2B gene with Hirschsprung disease in the Han Chinese population of Southeastern China. J Pediatr Surg. 2009;44:1805-11
Corresponding author: Ji-Cheng Li, Zhejiang University School of Medicine, No. 866, Yuhangtang Road, Hangzhou 310058, China. Tel: +86 571-88208088; Fax: +86 571-88208088. Email: lijichenedu.cn.