Background

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.

Methods

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.

Results

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.

Conclusions

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 [1]. 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 [4], and others [5]. 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[6], M. avium[7], 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 [13].

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 [14] and lipopeptides [15], 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 [18], 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 [20]). 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) [23]. 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 [24], 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.

Table 1

Characteristics of study population

Status

Control

Case

Total

P value*

Gender

Male

712

40.6%

1043

59.4%

1755

<0.001

Female

784

63.3%

455

36.7%

1239

Ethnicity

Mandinka

454

47.8%

495

52.2%

949

0.250

Jola

287

56.9%

217

43.1%

504

0.004

Fula

296

58.7%

208

41.3%

504

<0.001

Wollof

263

71.3%

106

28.7%

369

<0.001

Other

190

28.8%

470

71.2%

660

<0.001

Unknown

6

75.0%

2

25.0%

8

Total

1496

1498

2994

*Standard Chi-Square tests with 1 degree of freedom were performed for comparing groups.

Genotyping

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

Bioinformatic analysis

Potential alternative splice sites of MARCO were investigated using the Alternative Splice Site Predictor (ASSP) tool (http://wangcomputing.com/assp/index.html) [27]. Alignments of intron 1 and regions of high homology between multiple species were determined using ESPERR (evolutionary and sequence pattern extraction through reduced representations) [28]. 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) [29].

Statistical analysis

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 [30]. Haplotype blocks were defined as regions demonstrating strong evidence of historical recombination between less than 5% of single-nucleotide polymorphism (SNP)–pair comparisons [31]. All control genotype distributions were in Hardy-Weinberg equilibrium (0.05 level).

Figure 1

Haploview analysis of SNPs inMARCO[GeneID: 8685], located on chromosome 2q12. Linkage Disequilibrium (LD) in the African population is visualized across the MARCO locus, including upstream and downstream regions (chr2: 119414167–119468337). Polymorphisms are identified by their dbSNP rs numbers, and their position relative to the gene structure is marked. SNPs marked in bold were typed in this study. SNPs marked in green were included in the Wellcome Trust Case Control consortium [32] analysis and were used in validation studies. Empty squares indicate a high degree of LD (LD coefficient D' = 1) between pairs of markers. Numbers indicate the D' value expressed as a percentile. Black squares indicate pairs in strong LD with logarithm of odds (LOD) scores for LD ≥2; grey squares, D' < 1 with LOD ≥2; white squares, D' < 1.0 and LOD < 2.

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.

Table 2

Association analysis ofMARCOSNPs in a Gambian tuberculosis case control study

CHR

SNP

Position

Genotypic

Dominant

2

rs17180481

119414167

5.2

7.60E-02

4.2

4.10E-02

2

rs6748401

119414287

3.9

1.40E-01

0

8.30E-01

2

rs6748401

119414527

2.6

2.80E-01

0.5

4.90E-01

2

rs17009242

119415253

1.6

4.50E-01

1.6

2.10E-01

2

rs17009268

119415479

NA

NA

NA

NA

2

rs1073529

119421449

1.8

4.20E-01

1.7

1.90E-01

2

rs1985119

119426870

1.3

5.20E-01

1.2

2.70E-01

2

rs4491733

119434498

12.3

2.10E-03

6.5

1.10E-02

2

rs7559955

119441199

7.8

2.00E-02

6.9

8.50E-03

2

rs12998782

119441610

7.3

2.60E-02

6.5

1.10E-02

2

rs2119112

119442389

1

6.10E-01

0.4

5.10E-01

2

rs2278589

119445346

0

9.90E-01

0

9.60E-01

2

rs11693199

119448272

5.9

5.10E-02

3.1

7.90E-02

2

rs1371562

119449418

0.1

9.30E-01

0.1

7.10E-01

2

rs10204875

119453881

2.9

2.40E-01

1.4

2.30E-01

2

rs12997087

119454047

0.6

7.60E-01

0.4

5.20E-01

2

rs1371565

119454355

0.5

7.80E-01

0.5

4.90E-01

2

rs1371566

119454417

0.5

7.60E-01

0.4

5.50E-01

2

rs6761637

119455533

1.1

5.80E-01

1

3.20E-01

2

rs7599352

119456760

1.5

4.70E-01

1.1

2.90E-01

2

rs13389814

119458420

8

1.90E-02

6.9

8.50E-03

2

rs3731611

119459081

1.7

4.30E-01

1.6

2.10E-01

2

rs3731612

119459335

1.7

4.30E-01

1.6

2.10E-01

2

rs12987402

119465026

1.6

4.40E-01

1.3

2.50E-01

2

rs11678719

119468337

5.1

7.70E-02

0.7

4.00E-01

Table 3

The T allele of rs7559955 shows heterozygote protection

Control

Tuberculosis

Total

C/C

542

581

1123

C/T

651

553

1204

T/T

156

150

306

total

1349

1284

2633

40.2%

45.2%

48.3%

43.1%

11.6%

11.7%

HW P

0.060

0.295

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) [29] (Figure 2B).

Figure 2

TESS Analysis of potential transcription factor binding. The T variant allele creates potential for EFII and C/EBPalpha binding. A) The rs7559955 allele in intron 1 is found in a region of high conservation. Alignments of the homologous region from human, chimpanzee, macaque, mouse, rat, dog, and cow were performed, and the region of high conservation, and thus high regulatory potential as determined by ESPERR, is marked in blue. B) Transcription factor binding site analysis using TESS indicates that the presence of the T allele may confer the addition of EFII or C/EBPalpha transcription factor binding sites.

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 [35]). 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 [37] 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 [9]. 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 [23]. 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 [42]. Unfortunately, the intronic SNPs identified in our study and in the Chinese Han population [23] 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) [46]. 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 [23]. 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.

Acknowledgements

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.

Pre-publication history

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.