Evaluation of the relationship between MARCO and CD36 single-nucleotide polymorphisms and susceptibility to pulmonary tuberculosis in a Chinese Han population

Background 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. Methods 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. Results 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. Conclusions 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. Electronic supplementary material The online version of this article (doi:10.1186/s12879-017-2595-2) contains supplementary material, which is available to authorized users.


Background
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% [1]. 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 [2]. Polymorphisms in genes for the Toll-like receptor (TLR), interleukin (IL) [3,4], vitamin D receptor (VDR) [5] and interferongamma (IFN-γ) [6] 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 [7]. 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 [8]. Mice deficient in MARCO have shown an impaired ability to eliminate bacteria in the lungs [9]. In human alveolar macrophages, a specific protein that bound inhaled particles and bacteria was established to be MARCO [10]. Later, people used Mycobacterium marinum, a close relative of Mtb, to infect zebrafish, and found MARCO knockdown zebrafish suffered an increased bacterial burden [11]. Another study described how mice without MARCO showed a delay in the organization of marginal zone macrophages within the spleen [12]. Furthermore, MARCO recognized trehalose 6,6′-dimycolate (TDM; cord factor), which is one of the most nosogenetic components of Mtb [13].
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 [14]. Notably, those specific lipids and lipoproteins that make up the main components of mycobacterial organism cell walls are recognized by CD36 [15], [16]. However, CD36 −/− mice were found to suffer a reduced mycobacterial burden, indicating CD36 is an important determinant of host susceptibility to mycobacterial infection [17].
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.

DNA isolation
Genomic DNA was extracted from 200 μL of each whole blood sample using QIAamp a DNA Micro Kit (QIA-GEN, 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. [18]. The SNP, rs12998782, is associated with susceptibility to pulmonary TB in the Gambian population [19], 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 [20]. Rs10499859 was included in a study of left ventricular hypertrophy in a Korean population but did not show any significance [21]. As for rs3211956, a study of metabolic syndrome showed it had no association with this SNP [22]. 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 highresolution 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).

Statistical analysis
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.

Results
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 r 2 for the two MARCO SNPs and three CD36 SNPs. The two CD36 variants (rs1194182 and rs10499859) were in strong LD (D' = 0.646; r 2 = 0.302) (Fig. 1). Also the two MARCO variants (rs12998782 and rs17009726) were in strong LD (D' = 0.877; r 2 = 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.

Discussion
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 [19]. However, genotypes or gene alleles of rs17009726 showed no significant difference between cases and controls, which did     [18]. 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 [13]. 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 [25]. 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 [13], 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 [26]. 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 [17]. 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.
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. [18] 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) [25]. Whether MARCO and CD36 are associated with extra-pulmonary tuberculosis, such as TBM, or intestinal, bone and urinary tuberculosis, needs further investigation.

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

Additional files
Additional file 1: Table S1. PTB patients information. Some laboratory examination results and clinical findings were listed. (XLSX 32 kb) Additional file 2: Table S2.