We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Evaluation of the relationship between MARCO and CD36 single-nucleotide polymorphisms and susceptibility to pulmonary tuberculosis in a Chinese Han population
BMC Infectious Diseasesvolume 17, Article number: 488 (2017)
Gene polymorphisms impact greatly on a person’s susceptibility to pulmonary tuberculosis (PTB). Macrophage receptor with collagenous structure (MARCO) and CD36 are two scavenger receptors (SRs) that can recognize Mycobacterium tuberculosis (Mtb) and play a key role in tuberculosis infection. Gene polymorphisms of MARCO and CD36 may contribute to tuberculosis risk.
To investigate whether genetic polymorphisms of MARCO and CD36 are associated with susceptibility to PTB, genomic DNA samples from patients (n = 202) and healthy controls (n = 216) were collected and analyzed by polymerase chain reaction with high-resolution melting analysis.
We studied two single nucleotide polymorphisms (SNPs) in MARCO (rs12998782 and rs17009726) and three SNPs in CD36 (rs1194182, rs3211956 and rs10499859). Rs12998782 (P = 0.018) might be associated with susceptibility to PTB. Rs1194182 (P < 0.01) and rs10499859 (P < 0.001) might be associated with resistance to PTB. Rs17009726 and rs3211956 were not associated with susceptibility/resistance to PTB.
These data showed that MARCO rs12998782 may increase PTB risk while two SNPs of CD36, rs1194182 and rs10499859 may reduce the risk, indicating MARCO and CD36 as important receptors in response to PTB.
Tuberculosis (TB) has been a major health problem within human civilizations for thousands of years. After many years struggling against TB, an efficacious regimen has been established that has saved thousands of lives. However, with the appearance of drug resistance as well as an increase in the number of HIV-infected individuals who are susceptible to TB infection, this disease remains a threat. Globally in 2015, an estimated 10.4 million new cases of TB occurred, of which China, India, Nigeria, Indonesia, Pakistan and South Africa accounted for 60% . Though environmental factors and infection status have a major influence on the incidence of TB, evidence indicates that susceptibility or resistance of TB is also determined by genetic factors . Polymorphisms in genes for the Toll-like receptor (TLR), interleukin (IL) [3, 4], vitamin D receptor (VDR)  and interferon-gamma (IFN-γ)  are associated with susceptibility to tuberculosis.
TB is caused by Mycobacterium tuberculosis (Mtb), a pathogenic bacterium that can invade many organs, causing different types of TB, with PTB the most common. Mtb is an intracellular bacterium that mostly infects alveolar macrophages of the host . The cell membrane of macrophages contains SRs that were once thought to be related to the internalization of modified low-density lipoprotein (LDL). However, several studies showed that SRs also recognize several pathogens and play a crucial role in the body’s defense system. MARCO is a class A scavenger receptor that was thought to be upregulated in response to bacterial invasion . Mice deficient in MARCO have shown an impaired ability to eliminate bacteria in the lungs . In human alveolar macrophages, a specific protein that bound inhaled particles and bacteria was established to be MARCO . Later, people used Mycobacterium marinum, a close relative of Mtb, to infect zebrafish, and found MARCO knockdown zebrafish suffered an increased bacterial burden . Another study described how mice without MARCO showed a delay in the organization of marginal zone macrophages within the spleen . Furthermore, MARCO recognized trehalose 6,6′-dimycolate (TDM; cord factor), which is one of the most nosogenetic components of Mtb .
CD36 is another scavenger receptor and belongs to the class B family of scavenger receptors. CD36 binds many ligands such as oxidized low density lipoprotein, thrombospondin, and certain bacteria, as well as other ligands . Notably, those specific lipids and lipoproteins that make up the main components of mycobacterial organism cell walls are recognized by CD36 , . However, CD36−/− mice were found to suffer a reduced mycobacterial burden, indicating CD36 is an important determinant of host susceptibility to mycobacterial infection .
MARCO and CD36 are both pattern recognition receptors expressed on macrophages that can recognize Mtb; they connect innate and adaptive immunities. A few studies have investigated the relationship between genetic polymorphisms in MARCO and PTB susceptibility, but similar studies on CD36 are lacking. The aim of our study was to investigate whether gene polymorphisms of MARCO and CD36 are associated with susceptibility to PTB in a Chinese Han population.
We chose 202 PTB cases from the Tuberculosis Research Institute of Shenyang Tenth People’s Hospital from September 2015 to May 2016. Cases were diagnosed after radiographic examination and were positive for TB cultures or smears. Strain identification was performed after smear positive using gene chip. Patients with cancer, autoimmune disease, HIV and other pulmonary diseases were excluded. The healthy control group consisted of 216 uninfected participants from the Medical Examination Center of The First Affiliated Hospital of China Medical University, which matched cases in the same cohort according to age, sex, race and geographical location. All controls had no history of tuberculosis and were tested normal after blood routine examination, urine routine test, liver function test and radiographic examination. All patients were of Chinese Han descent. The demographic characteristics of study population were listed in Table 1. The gender ratios and mean ages were not significantly different between cases and controls. Details of clinical findings and laboratory examination results about patients and controls were listed in Additional files 1 and 2.
Genomic DNA was extracted from 200 μL of each whole blood sample using QIAamp a DNA Micro Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Extracted DNA was diluted to 30 ng/μL and was stored at −80 °C (ThermoFisher Scientific, Bartlesville, OK, USA).
Single nucleotide polymorphism selection and genotyping
We chose two SNPs of MARCO (rs12998782 and rs17009726) and three SNPs of CD36 (rs1194182, rs3211956 and rs10499859). Rs17009726 was studied in a Chinese Han population by Ma et al. . The SNP, rs12998782, is associated with susceptibility to pulmonary TB in the Gambian population , but a similar study has not been conducted in a Chinese Han population. The SNP, rs1194182, was reported to be related to cerebral malaria syndrome . Rs10499859 was included in a study of left ventricular hypertrophy in a Korean population but did not show any significance . As for rs3211956, a study of metabolic syndrome showed it had no association with this SNP . All SNPs were selected from an NCBI dbSNP database (http://www.ncbi.nlm.nih.gov/projects/SNP/) with a minor allele frequency (MAF)> 0.05 in a Chinese Han population.
We used polymerase chain reaction with high-resolution melting analysis (HRM-PCR; Roche Applied Science, Mannheim, Germany) for SNP genotyping. Primer Premier 6 was used to design primers with annealing temperatures ranging from 60 to 65 °C and amplicons were limited to 100–200 bp (Table 2). Real-time PCR cycling and HRM analyses were conducted using a Light Cycler® 480 System (Roche Applied Science). Samples were spiked with wild-type DNA as a standard. Wild-type and mutant homozygotes were distinguished by spiking samples with a known genotype sequence before PCR. Each sample of unknown genomic DNA (30 ng) was used as a template, along with an additional 3.0 ng (rs12998782, rs1194182, rs3211956, rs10499859) of known wild-type DNA. Because only two known genotypes of rs17009726 (wild-type and mutant heterozygous) exist, a supplement of wild-type DNA standard was not needed. PCRs were conducted in 96-well plates using touchdown PCR cycling. HRM curve data were obtained by melting over the range 65–95 °C at a rate of 25 data acquisitions per 1 °C. Results were analyzed using Light Cycler® 480 Gene Scanning software (Roche Applied Science). Forty randomly selected samples were verified by direct sequencing on an ABI7000 sequence detection system (Applied Biosystems, Carlsbad, CA, USA).
We performed statistics using SPSS 20.0. A Chi-square test was used to analyze genotype and gene allele frequencies in patient and control groups. Odd ratios (ORs) and 95% confidence intervals (CIs) were added to analyze the association between polymorphisms and the risk of tuberculosis. A two-tailed P value less than 0.05 was considered statistically significant. All control genotype distributions were in Hardy–Weinberg equilibrium.
As observed in Table 3, for the rs12998782 polymorphism of the MARCO gene, the presence of a TT homozygous mutation was higher in the patient group than in controls. While for rs17009726, neither allele nor genotype frequencies showed differences between the two groups.
As for rs1194182, the frequency of the G allele was slightly higher in the control group than in patients, while the GG homozygous mutation was significantly higher in the control group (CC + CG vs. GG). Furthermore, both the G allele and the GG homozygous mutation of rs10499859 showed a significant increase in the control group (P< 0.001) than in patients. However, a difference was not observed in rs3211956 between the two groups.
Taken together, our results suggest that polymorphisms of MARCO at rs12998782 might be associated with susceptibility to PTB, while polymorphisms of CD36 at rs1194182 and rs10499859, respectively, might be associated with protection against PTB.
We then combined 12,998,782 (from MARCO) with rs1194182 (from CD36), and rs12998782 (from MARCO) with 10,499,859 (from CD36; Table 4), to obtain nine combination genotypes, respectively. Interestingly, we found when two homozygous mutations appeared together (TT of rs12998782 and GG of rs1194182, or TT of rs12998782 and GG of rs10499859), a significant difference was not observed between patients and controls, thus indicating the opposite function of these polymorphisms in MARCO and CD36.
Then we estimated the LD block using a Haploview version 4.2. Pairwise LD was calculated by both D’ and r2 for the two MARCO SNPs and three CD36 SNPs. The two CD36 variants (rs1194182 and rs10499859) were in strong LD (D’ = 0.646; r2 = 0.302) (Fig. 1). Also the two MARCO variants (rs12998782 and rs17009726) were in strong LD (D’ = 0.877; r2 = 0.335) (Fig. 2). Further, we examined whether haplotypes of CD36 and MARCO were associated with PTB. For CD36, frequencies of haplotype CA were significantly higher in PTB group as compared to the control group and it may confer a higher risk of susceptibility effect for PTB. While the haplotype GG and CG played protective effect for PTB, the frequencies of haplotype GG and CG were significantly higher in control group (Table 5). However, none of the haplotypes of block1 in MARCO showed any significant difference between PTB and control group (datas not shown).
As nutrition status can affect the susceptibility to PTB, we performed stratified analysis according to albumin (ALB) levels (Table 6). We found in normal-ALB group (ALB> 40g/L), still TT of rs12998782 was significantly higher in PTB group than control group. Also GG of rs1194182 and rs10499859 showed an significant increase in control group. While in the low-ALB group, there was no difference in distribution of any genotypes or alleles between PTB group and control group.
Our results indicate that polymorphisms in the human class A scavenger receptor, MARCO, and the class B scavenger receptor, CD36, are associated with susceptibility/resistance to PTB in a Chinese Han population. According to our study, the TT genotype of rs12998782 in the MARCO gene can increase the risk of tuberculosis. This is consistent with findings from a recent study of a Gambian population that revealed rs12998782 was associated with susceptibility to PTB . However, genotypes or gene alleles of rs17009726 showed no significant difference between cases and controls, which did not coincide with the results of a previous study on a Chinese Han population performed by Ma et al. . Then we discussed separately according to age and sex, there was still no significant difference between PTB group and control group. Also we discussed according to treatment times (Initial treatment of PTB and retreatment of PTB), no difference was observed.
The MARCO gene is located on 2q14.2, and rs12998782 is located within intron 1. Though intron variants are not responsible for changes in amino acids, they may influence the transcription and splicing process. Studies on certain genes have revealed that intron 1 is essential for the full expression of rat genes [23, 24]. Thus, a variant in rs12998782 may impact the expression of MARCO on human alveolar macrophages, which may fail to recognize Mtb efficiently and may be defective in cooperating with TLR2/TLR4 . A recent study found that it is the reduced phagocytotic ability of pathogens, rather than reduced cytokine production, that explains the association between MARCO polymorphisms and susceptibility to PTB . A decreased ability to phagocytose Mtb will lead to a reduction of its clearance, which would increase susceptibility to PTB. However, due to a disparity between the experimental model (mice) and our investigated patients (human), mice with a MARCO knockout may represent a more notable defect in responding to Mtb , while for humans, MARCO polymorphisms may only have a slight impact.
On the other hand, we found the GG genotypes of rs1194182 and rs10499859 in the CD36 gene are associated with resistance to PTB. The CD36 gene is located on 7q21.11; rs1194182 is located in the 5′ untranslated region (5′ UTR), while rs10499859 is in intron 1. The 5′ UTR is a major site of translational regulation, a key step in eukaryotic gene expression . Single nucleotide polymorphisms in the 5′ UTR of the CD36 gene may affect its expression and the function of the protein. Thus, the GG homozygous mutation of rs1194182 may lead to the dysfunction of CD36. Interestingly, this result is consistent with a study performed by Hawkes et al., which demonstrated that CD36−/− mice suffered fewer mycobacterial infections than CD36+/+ mice . CD36 participates in the process of apoptosis and Mtb infected macrophages undergoing apoptosis can be phagocytosed by new macrophages induced by intracellular Mtb, thus leading to the expansion of infection [27, 28]. In other words, CD36 may cooperate with Mtb in infection. Hence, this may explain the homozygous mutations of rs1194182 and rs10499859 acting as protective factors against PTB.
We used Mutation Taster (http://www.mutationtaster.org/) to predict possible function changes of rs12998782 and rs10499859. We found these two SNPs were low conservative (PhyloP = −0.559, PhastCons = 0.008; PhyloP = −0.656, PhastCons = 0) and would lead to splice site changes and protein features might be affected.
As for the combination of MARCO and CD36, no significant difference was observed between patients and controls when two homozygous mutations from MARCO and CD36, respectively, coexisted. This indicates that polymorphisms of these two genes may have opposite effects on the susceptibility to PTB.
To date, this is the first time that the relationship between CD36 SNPs (rs1194182 and rs10499859) and PTB has been analyzed. For MARCO, Ma et al.  identified several SNPs that correlated with TB in a Chinese Han population; however, rs12998782 was not included in these studies.
One limitation of our study was the insufficient information available concerning patients and controls: More stratified analyses are required to make our study more enriched. In addition, extra-pulmonary tuberculosis also deserves attention. A recent study revealed two SNPs of MARCO, rs2278589 and rs6751745, were associated with PTB but not with TB meningitis (TBM) . Whether MARCO and CD36 are associated with extra-pulmonary tuberculosis, such as TBM, or intestinal, bone and urinary tuberculosis, needs further investigation.
In summary, our findings demonstrated that genetic variations in the scavenger receptors, MARCO and CD36, might be associated with susceptibility/resistance to pulmonary tuberculosis in a Chinese Han population. It remains to be determined whether genetic variations in MARCO and CD36 alter the expression of these receptors and whether genetic polymorphisms impact the severity of infection and the persistence of Mtb over a longer time period.
Cluster of differentiation 36
Polymerase chain reaction with high-resolution melting analysis.
Minor allele frequency
Macrophage receptor with collagenous structure
Single nucleotide polymorphism
Global Tuberculosis Report 2016. World Health Organization (WHO) 2016. Available from: http://www.who.int/tb/publications/global_report/en/
Bellamy R. Susceptibility to mycobacterial infections: the importance of host genetics. Genes Immun. 2003;4:4–11.
Wang M, Xu G, Lu L, Xu K, Chen Y, Pan H, et al. Genetic polymorphisms of IL-17A, IL-17F, TLR4 and miR-146a in association with the risk of pulmonary tuberculosis. Sci Rep. 2016;6:28586.
Salie M, Daya M, Lucas LA, Warren RM, van der Spuy GD, van Helden PD, et al. Association of toll-like receptors with susceptibility to tuberculosis suggests sex-specific effects of TLR8 polymorphisms. Infect Genet Evol. 2015;34:221–9.
Hu Q, Chen Z, Liang G, Mo F, Zhang H, Xu S, et al. Vitamin D receptor gene associations with pulmonary tuberculosis in a Tibetan Chinese population. BMC Infect Dis. 2016;16:469.
Shamsi M, Zolfaghari MR, Farnia P. Association of IFN-gamma and P2X7 receptor gene polymorphisms in susceptibility to tuberculosis among Iranian patients. Acta Microbiol Immunol Hung. 2016;63:93–101.
Killick KE, Ni Cheallaigh C, O'Farrelly C, Hokamp K, MacHugh DE, Harris J. Receptor-mediated recognition of mycobacterial pathogens. Cell Microbiol. 2013;15:1484–95.
Luc JW, van Laan d, Kangas M, Do pp EA, Broug-Holub E, Elomaa O, Tryggvason K, et al. Macrophage scavenger receptor MARCO: in vitro and in vivo regulation and involvement in the anti-bacterial host defense. Immunol Lett. 1997;57:203–8.
Arredouani M, Yang Z, Ning Y, Qin G, Soininen R, Tryggvason K, et al. The scavenger receptor MARCO is required for lung defense against pneumococcal pneumonia and inhaled particles. J Exp Med. 2004;200:267–72.
Arredouani MS, Palecanda A, Koziel H, Huang YC, Imrich A, Sulahian TH, et al. MARCO is the major binding receptor for Unopsonized particles and bacteria on human alveolar macrophages. J Immunol. 2005;175:6058–64.
Benard EL, Roobol SJ, Spaink HP, Meijer AH. Phagocytosis of mycobacteria by zebrafish macrophages is dependent on the scavenger receptor Marco, a key control factor of pro-inflammatory signalling. Dev Comp Immunol. 2014;47:223–33.
Chen Y, Pikkarainen T, Elomaa O, Soininen R, Kodama T, Kraal G, et al. Defective microarchitecture of the spleen marginal zone and impaired response to a thymus-independent type 2 antigen in mice lacking scavenger receptors MARCO and SR-A. J Immunol. 2005;175:8173–80.
Bowdish DM, Sakamoto K, Kim MJ, Kroos M, Mukhopadhyay S, Leifer CA, et al. MARCO, TLR2, and CD14 are required for macrophage cytokine responses to mycobacterial trehalose dimycolate and mycobacterium tuberculosis. PLoS Pathog. 2009;5:e1000474.
Silverstein RL, Febbraio M. CD36, a scavenger receptor involved in immunity, metabolism, angiogenesis, and behavior. Sci Signal. 2009;2:re3.
Hoebe K, Georgel P, Rutschmann S, Du X, Mudd S, Crozat K, et al. CD36 is a sensor of diacylglycerides. Nature. 2005;433
Philips JA, Rubin EJ, Perrimon N. Drosophila RNAi screen reveals CD36 family member required for Mycobaterial infection. Sci. 2005;309
Hawkes M, Li X, Crockett M, Diassiti A, Finney C, Min-Oo G, et al. CD36 deficiency attenuates experimental mycobacterial infection. BMC Infect Dis. 2010;10:299.
Ma MJ, Wang HB, Li H, Yang JH, Yan Y, Xie LP, et al. Genetic variants in MARCO are associated with the susceptibility to pulmonary tuberculosis in Chinese Han population. PLoS One. 2011;6:e24069.
Bowdish DM, Sakamoto K, Lack NA, Hill PC, Sirugo G, Newport MJ, et al. Genetic variants of MARCO are associated with susceptibility to pulmonary tuberculosis in a Gambian population. BMC Med Genet. 2013;14:47.
Sambo MR, Trovoada MJ, Benchimol C, Quinhentos V, Goncalves L, Velosa R, et al. Transforming growth factor beta 2 and heme oxygenase 1 genes are risk factors for the cerebral malaria syndrome in Angolan children. PLoS One. 2010;5:e11141.
Park JK, Kim MK, Choi BY, Jung Y, Song K, Kim YM, et al. Validation study of candidate single nucleotide polymorphisms associated with left ventricular hypertrophy in the Korean population. BMC Med Genet. 2015;16:13.
Farook VS, Puppala S, Schneider J, Fowler SP, Chittoor G, Dyer TD, et al. Metabolic syndrome is linked to chromosome 7q21 and associated with genetic variants in CD36 and GNAT3 in Mexican Americans. Obesity. 2012;20:2083–92.
R. CP, Klett, Bonner TI. Identification and characterization of the rat M1 Muscarinic receptor promoter. J Neurochem 1999;72.
Kim C-H, Kim H-S, Cubells JF, Kim K-S. A previously Undescribed Intron and extensive 5′ upstream sequence, but not Phox2a-mediated Transactivation are necessary for high level cell type-specific expression of the human Norepinephrine transporter gene. J Biol Chem. 1999;274:6507–18.
Thuong NT, Tram TT, Dinh TD, Thai PV, Heemskerk D, Bang ND, et al. MARCO variants are associated with phagocytosis, pulmonary tuberculosis susceptibility and Beijing lineage. Genes Immun. 2016;17(7):419-25.
Pickering BM, Willis AE. The implications of structured 5′ untranslated regions on translation and disease. Semin Cell Dev Bio. 2005;16:39–47.
Voll RE, Herrmann M, Roth EA, Stach C, Kalden JR. Immunosuppressive effects of apoptotic cells. Nature. 1997;39
Davis JM, Ramakrishnan L. The role of the granuloma in expansion and dissemination of early tuberculous infection. Cell. 2009;136:37–49.
We’d like to thank all of the patients and individuals for their participation.
This research was supported by Shenyang science and technology plan project (F13–220–9-04).
Availability of data and materials
All data generated or analyzed during this study are included in this published article. Consent for publication of raw data not obtained but dataset is fully anonymous in a manner that can easily be verified by any user of the dataset. Publication of the dataset clearly and obviously presents minimal risk to confidentiality of study participants. The Ethics Committee of the First Affiliated Hospital of China Medical University have approved the publication of raw data.
WL designed the study, carried out the experiment, participated in the statistical analysis and drafted the manuscript. HK participated in the design and directed the experiment and analysis. GJ directed the experiment operation, participated in statistical analysis and assisted in drafting the manuscript. JS and BS participated in the sample collections and DNA isolation. LC and YC assisted in SNPs selection, primer design and experiment operation. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Our study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University. Written informed consent was obtained from patients as delineated by the protocol, which was also approved by the Ethics Committee of the First Affiliated Hospital of China Medical University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.