- Research article
- Open Access
Genetic variants of MARCO are associated with susceptibility to pulmonary tuberculosis in a Gambian population
BMC Medical Genetics volume 14, Article number: 47 (2013)
The two major class A scavenger receptors are scavenger receptor A (SRA), which is constitutively expressed on most macrophage populations, and macrophage receptor with collagenous structure (MARCO), which is constitutively expressed on a more restricted subset of macrophages, (e.g. alveolar macrophages) but whose expression increases on most macrophages during the course of infection. Although the primary role of SRA appears to be clearance of modified host proteins and lipids, mice defective in expression of either MARCO or SRA are immunocompromised in multiple models of infection and in vitro assays, the scavenger receptors have been demonstrated to bind bacteria and to enhance pro-inflammatory signalling to many bacterial lung pathogens; however their importance in Mycobacterium tuberculosis infection, is less clear.
To determine whether polymorphisms in either SRA or MARCO were associated with tuberculosis, a case–control study of was performed. DNA samples from newly-detected, smear-positive, pulmonary tuberculosis cases were collected from The Gambia. Controls for this study consisted of DNA from cord bloods obtained from routine births at local Gambian health clinics. Informed written consent was obtained from patients or their parents or guardians. Ethical approval was provided by the joint The Gambian Government/MRC Joint Ethics Committee.
We studied the frequencies of 25 polymorphisms of MSR1 (SRA) and 22 in MARCO in individuals with tuberculosis (n=1284) and matched controls (n=1349). No SNPs within the gene encoding or within 1 kb of the promoter sequence of MSR1 were associated with either susceptibility or resistance to tuberculosis. Three SNPs in MARCO (rs4491733, Mantel-Haenszel 2x2 χ2 = 6.5, p = 0.001, rs12998782, Mantel-Haenszel 2x2 χ2 = 6.59, p = 0.001, rs13389814 Mantel-Haenszel 2x2 χ2 = 6.9, p = 0.0009) were associated with susceptibility to tuberculosis and one (rs7559955, Mantel-Haenszel 2x2 χ2 = 6.9, p = 0.0009) was associated with resistance to tuberculosis.
These findings identify MARCO as a potentially important receptor in the host response to tuberculosis.
Mycobacterium tuberculosis, the causative agent of tuberculosis, is one of few pathogens that infects and persists in host macrophages. Consequently it has successfully colonized up to a third of the world’s population and of those, many millions develop active disease . Phagocytosis of M. tuberculosis is mediated by a number of receptors including the mannose receptor and DC-SIGN, which recognise mannose-capped lipoarabinomannan (Man-lam) [2, 3], complement receptor via recognition of opsonin-coated bacteria , and others . Although the class A scavenger receptors, SRA (class A scavenger receptor) and MARCO (macrophage receptor with collagenous structure), are broadly classified as “phagocytic receptors”, and have been demonstrated to internalize mycobacterial species such as M. leprae, M. avium, M. bovis Bacille Calmette-Guérin [8, 9], and M. tuberculosis [10, 11], it has been proposed that there is functional redundancy in mycobacterial uptake by macrophages so that any individual receptor is dispensable [10, 12]. Intriguingly, whereas uptake of M. tuberculosis into alveolar macrophages during the course of acute infection may be mediated by many equivalent receptors, long-term extra-pulmonary persistence (e.g. in the adipose tissue), may be mediated through scavenger receptor uptake .
Although the importance of phagocytic receptors in uptake of M. tuberculosis is not entirely clear, their role in induction of pro-inflammatory responses appears to be more straightforward. Class A and B scavenger receptors are required for maximal cytokine responses to mycobacterial lipoarabinomannans  and lipopeptides , and both the class A scavenger receptors and the C-type lectin, Mincle, are required for optimal toll-like receptor (TLR) and Syk/Card9 signaling responses to mycobacterial trehalose dimycolate (cord factor), respectively [16, 17]. The biological importance of this enhancement of cytokine responses remains to be fully elucidated; however, in lung infection models of M. tuberculosis, the absence of SRA is protective , but in models of disseminated disease, its absence is fatal [8, 11, 19].
To some degree, the lack of clarity surrounding whether these receptors are of key importance in host defence towards tuberculosis is probably due to deficiencies in the mouse model. The lung pathology of tuberculosis is sufficiently different in mice, in that many of the hallmark features of disease (e.g. granulomas) do not occur, and genes found to be associated in human studies are not necessarily associated with murine susceptibility and vice versa (reviewed in ). In order to determine whether the class A scavenger receptors are crucial to protection against human disease, we have performed a case–control study of single nucleotide polymorphisms (SNPs) using samples from a well-described Gambian population [21, 22]. We demonstrate that polymorphisms in MARCO, but not MSR1, the gene encoding SRA, are associated with tuberculosis infection in the Gambian population. Encouragingly, these results are consistent with a recently published case–control study in a Chinese Han population (n= 923 cases and 1033 controls) . In both the Chinese Han and Gambian populations, SNPs in intron 1, which we identify as a putative alternative promoter site, are associated with susceptibility to tuberculosis, implying that changes in MARCO function or expression contribute to host defence against tuberculosis.
Patient DNA Samples
DNA samples from newly-detected, smear-positive, pulmonary tuberculosis cases and healthy controls were collected from The Gambia as described [21, 22]. Samples were comprised of cases and controls (please see Table 1 for details). Gambian pulmonary tuberculosis cases (n=1,498) presented with a compatible clinical picture of tuberculosis and were diagnosed with culture or smear positivity; all cases that were smear positive but culture negative had radiographic confirmation. Exclusion criteria included presentation of autoimmune, cancer, or other diseases, such as HIV-1, which are known to impact host immunity. The majority (>95%) of the tuberculosis cohort was screened for HIV-1, with positive cases excluded from the study because HIV infection increases the risk of tuberculosis. Gambian controls (n=1,496) were recruited from routine births at local Gambian health clinics. Some samples were removed from the study due to quality control issues, including low genotype frequency and relatedness, leaving 1284 cases and 1349 controls for the final analysis. These cases and controls have been previously analyzed studying tuberculosis (46) and malaria , and we have a detailed understanding of the underlying ethnic stratification through genome-wide SNP data. Logistic regression based upon 6-axis of variation (to reduce the impact of ethnicity on the final association analysis) did not affect our final conclusions. Informed written consent was obtained from patients or their parents or guardians. Ethical approval was provided by the joint The Gambian Government/MRC Joint Ethics Committee.
Polymorphisms in MSR1 and MARCO genes were identified from the National Center for Biotechnology Information dbSNP database (http://www.ncbi.nlm.nih.gov/projects/SNP/). SNPs were genotyped using the Sequenom (San Diego) MassARRAY system, and the primer extension products were analyzed using MALDI-TOF mass spectrometry [25, 26].
Potential alternative splice sites of MARCO were investigated using the Alternative Splice Site Predictor (ASSP) tool (http://wangcomputing.com/assp/index.html) . Alignments of intron 1 and regions of high homology between multiple species were determined using ESPERR (evolutionary and sequence pattern extraction through reduced representations) . To investigate whether there might be conserved transcriptional elements in intron 1 of the gene encoding MARCO, the Transcription Element Search System (TESS) program was used (http://www.cbil.upenn.edu/tess) .
Statistical analysis of genotype associations was performed using SPSS version 12.0 (SPSS, Inc., Chicago, IL). Analysis of linkage disequilibrium (LD) and haplotypes (Figure 1) was performed using the Haploview version 3.2 program . Haplotype blocks were defined as regions demonstrating strong evidence of historical recombination between less than 5% of single-nucleotide polymorphism (SNP)–pair comparisons . All control genotype distributions were in Hardy-Weinberg equilibrium (0.05 level).
SNPs in MARCO are associated with resistance/susceptibility to tuberculosis
We chose 25 tag SNPs in MARCO and 22 in MSR1; preliminary analysis indicated that these SNPs were polymorphic and found at reasonable frequencies in the Gambian population (e.g. > 5%). None of the SNPs genotyped in MSR1 were associated with tuberculosis susceptibility (Additional file 1 Table S1). Three SNPs in MARCO showed some evidence of association with susceptibility to tuberculosis (p < 0.02, Table 2), and for one of these, rs7559955, heterozygotes showed reduced risk of tuberculosis (p = 0.008, Table 3). These SNPs were located within two regions of the MARCO gene. The first cluster of three SNPs is located in intron 1 of the gene (rs4491733, rs7559955, rs12998782), and 1 SNP (rs13389814) was found in intron 16.
The T variant in rs7559955 may confer a regulatory element
The presence of a polymorphism in rs7559955 is associated with resistance to tuberculosis and may confer heterozygote protection (Table 3). Bioinformatic analysis using ESPERR indicates that the 200 bp region surrounding rs7559955 contains a relatively short sequence that is highly conserved amongst seven species (human, chimpanzee, macaque, mouse, rat, dog, and cow) and is thus a potential regulatory sequence (Figure 2A). Score values above the 0.1 default upper limit indicate very marked resemblance to alignment patterns typical of regulatory elements [33, 34]. We thus hypothesized that this region might be a transcription factor binding site. Analysis using Transcription Element Search System (TESS) indicated that there was a weak match in both the ancestral and variant allele for the transcription factor NF-1 (score 8.0). The T allele appears to result in the addition of a transcription factor binding site that has homology to that of EFII (score 9.155) or C/EBPalpha(score 10)  (Figure 2B).
Herein, we demonstrate that polymorphisms within the human class A scavenger receptor MARCO correlate with susceptibility/resistance to tuberculosis in a Gambian population. Interestingly, none of the SNPs identified in this study encoded non-synonymous mutations, but rather were within introns (i.e. rs13389814 within intron 16) and putative regulatory regions (i.e. rs7559955 within a potential promoter site in intron 1). This may imply that changes in MARCO expression rather than overt changes in structure are key to conferring resistance or susceptibility to tuberculosis. In fact, although SRA (encoded by MSR1) shares high amino acid similarity with MARCO (76% amino acid similarity) and has many overlapping ligands, the two receptors differ considerably at the level of regulation of expression. Whereas SRA expression is regulated primarily by factors associated with macrophage differentiation (e.g. GM-CSF) and lipid accumulation, MARCO is primarily regulated by inflammation and bacterial infection (reviewed in ). Unlike SRA, which is expressed on virtually all macrophages, MARCO is constitutively expressed on some subsets of macrophages (i.e. alveolar, peritoneal, marginal zone of the spleen, and medullary cords of the lymph nodes), and is rapidly up-regulated on others during the course of infection. Indeed, in vivo MARCO expression is quickly (<1 h) increased upon challenge with bacteria or bacterial products, and this increase in expression occurs even at sites distal to bacterial challenge (e.g. on Kupffer cells post lung infection) [36, 37]. Changes in MARCO expression during the course of experimental tuberculosis infections have not been performed, but expression of MARCO increases upon systemic challenge with BCG  and is higher on macrophages within and proximal to BCG-containing granulomas. MARCO-expressing macrophages in the marginal zone appear to phagocytose more BCG than neighbouring macrophages that do not express MARCO . Whether the presence of a SNP within intron 1 (rs7559955) results in enhanced or decreased expression of MARCO, constitutively or after exposure to pathogens or their products, and whether changes in expression alter binding, uptake, or cytokine production, requires further investigation.
Recently Ma et al. reported that SNPs in MARCO were associated with susceptibility to tuberculosis in the Chinese Han population . Interestingly, they identified that rs17009716 (chromosomal position 119441930), was significantly associated with susceptibility to tuberculosis. This SNP is also within intron 1, less than 400 base pairs from and in high LD with rs12998782 and rs7559955 in Chinese (CHB), European (CEU), and African (YRI) populations. This provides support for potential consistency of a MARCO association across diverse populations and suggests a role for MARCO, in host defence against tuberculosis, possibly at the level of gene regulation, although further functional studies are warranted. Although animal models using MARCO knockout mice strongly associate MARCO with many pulmonary infections, such as pneumonia and influenza [38–41], whether polymorphisms in MARCO are associated with susceptibility to infection in humans remains to be determined. A recent Danish study found that non-synonymous SNPs within MARCO were found at very low frequency (0.005-5%) within a Danish population, and in a study investigating whether MARCO was associated with chronic obstructive pulmonary disease (COPD) and lung infection in COPD patients, only one of these was associated with an increased risk of sepsis and none with pneumonia or COPD itself . Unfortunately, the intronic SNPs identified in our study and in the Chinese Han population  were not investigated. It will be essential to validate these studies in more diverse populations, and perhaps to investigate more SNPs, in order to confirm the relevance of these results.
The observed association of MARCO variants with sepsis is consistent with in vivo[36, 37, 43, 44] and in vitro[44, 45] observations that MARCO regulates pro-inflammatory cytokine production to whole bacteria and TLR ligands. The class A scavenger receptors in general, and MARCO in particular, have been demonstrated to enhance TLR signalling to mycobacterial cell wall components, such as the major immunogenic lipid on the mycobacterial cell wall, trehalose dimycolate (TDM) . Whether this increased regulation of cytokines is of benefit or is detrimental to the host during the course of either pulmonary or systemic tuberculosis infection remains to be determined.
In summary, we have demonstrated that genetic variation in a class A scavenger receptor, MARCO, is associated with susceptibility/resistance to tuberculosis in a Gambian population, consistent with previously reported data from a Chinese Han population . Further studies are warranted to determine whether genetic variation in MARCO may alter expression of the receptor, and whether regulation of MARCO expression contributes to protection or susceptibility to tuberculosis at the level of macrophage infection, initiation of a pro-inflammatory response, or long-term persistence.
Single nucleotide polymorphism
Macrophage receptor with collagenous structure
Wellcome Trust Case Control Consortium
Scavenger receptor cysteine-rich domain
Peripheral blood mononuclear cell
Organization, W.H. Global tuberculosis control: WHO report 2011. 2011, cited 2011; Available from: http://apps.who.int/iris/bitstream/10665/75938/1/9789241564502_eng.pdf
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 (7): 2920-2930.
Maeda N: The cell surface receptor DC-SIGN discriminates between mycobacterium species through selective recognition of the mannose caps on lipoarabinomannan. J Biol Chem. 2003, 278 (8): 5513-5516. 10.1074/jbc.C200586200.
Schlesinger LS: Phagocytosis of mycobacterium tuberculosis is mediated by human Monocyte complement receptors and complement component C3. J Immunol. 1990, 144 (7): 2771-2780.
Rothfuchs AG: Dectin-1 interaction with mycobacterium tuberculosis leads to enhanced IL-12p40 production by splenic dendritic cells. J Immunol. 2007, 179 (6): 3463-3471.
Montoya D: Divergence of macrophage phagocytic and antimicrobial programs in leprosy. Cell Host Microbe. 2009, 6 (4): 343-353. 10.1016/j.chom.2009.09.002.
Bermudez LE, Parker A, Goodman JR: Growth within macrophages increases the efficiency of mycobacterium avium in invading other macrophages by a complement receptor-independent pathway. Infect Immun. 1997, 65 (5): 1916-1925.
Haworth R: The macrophage scavenger receptor type a is expressed by activated macrophages and protects the host against lethal endotoxic shock. J Exp Med. 1997, 186 (9): 1431-1439. 10.1084/jem.186.9.1431.
Ito S: Roles of a macrophage receptor with collagenous structure (MARCO) in host defense and heterogeneity of splenic marginal zone macrophages. Arch Histol Cytol. 1999, 62 (1): 83-95. 10.1679/aohc.62.83.
Zimmerli S, Edwards S, Ernst JD: Selective receptor blockade during phagocytosis does not alter the survival and growth of mycobacterium tuberculosis in human macrophages. Am J Respir Cell Mol Biol. 1996, 15 (6): 760-770. 10.1165/ajrcmb.15.6.8969271.
Pedroza-Gonzalez A: In situ analysis of lung antigen-presenting cells during murine pulmonary infection with virulent mycobacterium tuberculosis. Int J Exp Pathol. 2004, 85 (3): 135-145. 10.1111/j.0959-9673.2004.00381.x.
Court N: Partial redundancy of the pattern recognition receptors, scavenger receptors, and C-type lectins for the long-term control of mycobacterium tuberculosis infection. J Immunol. 2010, 184 (12): p 7057-p 7070.
Neyrolles O: Is adipose tissue a place for mycobacterium tuberculosis persistence?. PLoS One. 2006, 1: e43-10.1371/journal.pone.0000043.
Jozefowski S: Mycobacterium tuberculosis lipoarabinomannan enhances LPS-induced TNF-alpha production and inhibits NO secretion by engaging scavenger receptors. Microb Pathog. 2011, 50 (6): p.350-p.359.
Drage MG: TLR2 And its co-receptors determine responses of macrophages and dendritic cells to lipoproteins of mycobacterium tuberculosis. Cell Immunol. 2009, 258 (1): 29-37. 10.1016/j.cellimm.2009.03.008.
Schoenen H: Cutting edge: mincle is essential for recognition and adjuvanticity of the mycobacterial cord factor and its synthetic analog trehalose-dibehenate. J Immunol. 2010, 184 (6): p 2756-p 2760.
Bowdish DM: MARCO, TLR2, and CD14 are required for macrophage cytokine responses to mycobacterial trehalose dimycolate and mycobacterium tuberculosis. PLoS Pathog. 2009, 5 (6): e1000474-10.1371/journal.ppat.1000474.
Sever-Chroneos Z: Prolonged survival of scavenger receptor class A-deficient mice from pulmonary mycobacterium tuberculosis infection. Tuberculosis (Edinb). 2011, 91 Suppl 1: p. S69-p. S74.
Suzuki H: A role for macrophage scavenger receptors in atherosclerosis and susceptibility to infection. Nature. 1997, 386 (6622): 292-10.1038/386292a0.
Apt A, Kramnik I: Man and mouse TB: contradictions and solutions. Tuberculosis (Edinb). 2009, 89 (3): 195-198. 10.1016/j.tube.2009.02.002.
Bennett S: Investigation of environmental and host-related risk factors for tuberculosis in Africa II. Investigation of host genetic factors. Am J Epidemiol. 2002, 155 (11): 1074-1079. 10.1093/aje/155.11.1074.
Lienhardt C: Investigation of environmental and host-related risk factors for tuberculosis in africa I. Methodological aspects of a combined design. Am J Epidemiol. 2002, 155 (11): 1066-1073. 10.1093/aje/155.11.1066.
Ma MJ: Genetic variants in MARCO are associated with the susceptibility to pulmonary tuberculosis in Chinese Han population. PLoS One. 2011, 6 ((8): p. e24069-
Jallow M: Genome-wide and fine-resolution association analysis of malaria in west Africa. Nat Genet. 2009, 41 (6): 657-665. 10.1038/ng.388.
Jurinke C: The use of MassARRAY technology for high throughput genotyping. Adv Biochem Eng Biotechnol. 2002, 77: 57-74.
Jurinke C: Automated genotyping using the DNA MassArray technology. Methods Mol Biol. 2002, 187: 179-192.
Wang M, Marin A: Characterization and prediction of alternative splice sites. Gene. 2006, 366 (2): 219-227. 10.1016/j.gene.2005.07.015.
Taylor J: ESPERR: learning strong and weak signals in genomic sequence alignments to identify functional elements. Genome Res. 2006, 16 (12): 1596-1604. 10.1101/gr.4537706.
Schug J: Using TESS to predict transcription factor binding sites in DNA sequence. Curr Protoc Bioinformatics. 2008, Chapter 2: p. Unit 2-p. Unit 6.
Barrett JC: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Gabriel SB: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.
Thye T: Genome-wide association analyses identifies a susceptibility locus for tuberculosis on chromosome 18q11.2. Nat Genet. 2010, 42 (9): p 739-p 741.
King DC: Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences. Genome Res. 2005, 15 (8): 1051-1060. 10.1101/gr.3642605.
Kolbe D: Regulatory potential scores from genome-wide three-way alignments of human, mouse, and rat. Genome Res. 2004, 14 (4): 700-707. 10.1101/gr.1976004.
Bowdish DM, Gordon S: Conserved domains of the class a scavenger receptors: evolution and function. Immunol Rev. 2009, 227 (1): 19-31. 10.1111/j.1600-065X.2008.00728.x.
van der Laan LJ: Macrophage scavenger receptor MARCO: in vitro and in vivo regulation and involvement in the anti-bacterial host defense. Immunol Lett. 1997, 57 (1–3): 203-208.
van der Laan LJ: Regulation and functional involvement of macrophage scavenger receptor MARCO in clearance of bacteria in vivo. J Immunol. 1999, 162 (2): 939-947.
Arredouani M: The scavenger receptor MARCO is required for lung defense against pneumococcal pneumonia and inhaled particles. J Exp Med. 2004, 200 (2): 267-272. 10.1084/jem.20040731.
Arredouani MS: The macrophage scavenger receptor SR-AI/II and lung defense against pneumococci and particles. Am J Respir Cell Mol Biol. 2006, 35 (4): 474-478. 10.1165/rcmb.2006-0128OC.
Palecanda A, Kobzik L: Receptors for unopsonized particles: the role of alveolar macrophage scavenger receptors. Curr Mol Med. 2001, 1 (5): 589-595. 10.2174/1566524013363384.
Ghosh S: MARCO regulates early inflammatory responses against influenza: a useful macrophage function with adverse outcome. Am J Respir Cell Mol Biol. 2011, 45 (5): p 1036-p 1044.
Thomsen M: Genetic variation in the scavenger receptor MARCO and its association with chronic obstructive pulmonary disease and lung infection in 10,604 individuals. Respiration.
Elomaa O: Structure of the human macrophage MARCO receptor and characterization of its bacteria-binding region. J Biol Chem. 1998, 273 (8): 4530-8. 10.1074/jbc.273.8.4530.
Dorrington M, Roche A, Chauvin SE, Tu Z, Mossman K, Weiser JN, Bowdish DME: MARCO is required for TLR2- and NOD2-mediated responses to streptococcus pneumoniae and clearance of pneumococcal colonization in the murine nasopharynx. J Immunol. 2013, 190 (1): p 250-8.
Jozefowski S: Disparate regulation and function of the class a scavenger receptors SR-AI/II and MARCO. J Immunol. 2005, 175 (12): 8032-41.
Bowdish DME, Sakamoto K, Mukhopadhyay S, Leifer CA, Tryvaggson K, Gordon S, Russell DG: MARCO is a binding receptor for mycobacterial trehalose 6,6’-dimycolate and is required to initiate TLR2 induced signaling. In preparation. 2008
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/14/47/prepub
This work was funded by the Wellcome Trust. A.V.S.H. was supported by a Wellcome Trust Principal Fellowship. This work was funded in part by NIH Grant Number 1 R15 AI094436-01A1 to KS. DMEB would like to thank Dr. Chris Verschoor for critical reading of the manuscript. Work in the Bowdish laboratory is supported in part by the McMaster Immunology Research Centre and the M. G. DeGroote Institute for Infectious Disease Research.
The authors have no competing interests.
DMEB, KS, NL and FOV performed experiments. PCH, GS, MN, AVSH provided infrastructure and support for sample collection and processing. DMEB, SG, AVSH, FOV designed the experiments. FOV performed the statistical analysis. KS, SG, and AVSH provided funding. All authors read, edited, commented and approved on the manuscript.
Electronic supplementary material
About this article
Cite this article
Bowdish, D.M., Sakamoto, K., Lack, N.A. et al. Genetic variants of MARCO are associated with susceptibility to pulmonary tuberculosis in a Gambian population. BMC Med Genet 14, 47 (2013). https://doi.org/10.1186/1471-2350-14-47
- Scavenger receptors
- Mycobacterium tuberculosis
- Single nucleotide polymorphisms
- Case control study