Human gene expression profiles of susceptibility and ... - Nature

177 downloads 69 Views 528KB Size Report
Sep 23, 2010 - Sciences, Stellenbosch University, Cape Town, South Africa. Tuberculosis (TB) still ... airways, rendering the disease contagious.12 Sensing of.
Genes and Immunity (2011) 12, 15–22 & 2011 Macmillan Publishers Limited All rights reserved 1466-4879/11 www.nature.com/gene

ORIGINAL ARTICLE

Human gene expression profiles of susceptibility and resistance in tuberculosis J Maertzdorf1, D Repsilber2, SK Parida1, K Stanley3, T Roberts3, G Black3, G Walzl3 and SHE Kaufmann1 Department of Immunology, Max Planck Institute for Infection Biology, Berlin, Germany; 2Leibniz Institute for Farm Animal Biology, Genetics and Biometry, Dummerstorf, Germany and 3Division of Molecular Biology and Human Genetics, MRC Centre for Molecular and Cellular Biology, DST and NRF Centre of Excellence for Biomedical TB Research, Faculty of Health Sciences, Stellenbosch University, Cape Town, South Africa 1

Tuberculosis (TB) still poses a profound burden on global health, owing to significant morbidity and mortality worldwide. Although a fully functional immune system is essential for the control of Mycobacterium tuberculosis infection, the underlying mechanisms and reasons for failure in part of the infected population remain enigmatic. Here, whole-blood microarray gene expression analyses were performed in TB patients and in latently as well as uninfected healthy controls to define biomarkers predictive of susceptibility and resistance. Fc gamma receptor 1B (FCGRIB)was identified as the most differentially expressed gene, and, in combination with four other markers, produced a high degree of accuracy in discriminating TB patients and latently infected donors. We determined differentially expressed genes unique for active disease and identified profiles that correlated with susceptibility and resistance to TB. Elevated expression of innate immune-related genes in active TB and higher expression of particular gene clusters involved in apoptosis and natural killer cell activity in latently infected donors are likely to be the major distinctive factors determining failure or success in controlling M. tuberculosis infection. The gene expression profiles defined in this study provide valuable clues for better understanding of progression from latent infection to active disease and pave the way for defining predictive correlates of protection in TB. Genes and Immunity (2011) 12, 15–22; doi:10.1038/gene.2010.51; published online 23 September 2010 Keywords: tuberculosis; microarray; biomarkers

Introduction Tuberculosis (TB) is an ancient disease caused by the bacterial pathogen Mycobacterium tuberculosis, which has co-evolved with the human population for millennia.1 Infection with M. tuberculosis is controlled in most infected individuals, while in a minority of cases containment of infection fails over time, and active TB disease develops.2 Although only a minority of exposed individuals develop active disease, morbidity and mortality in humans worldwide remain high, with nearly 2 million deaths annually.3,4 In the face of this major impact on global health, remarkably little progress has been made in the past towards the development of new tools and elucidation of key features involved in protection and pathogenesis.5–7 It is generally accepted that a competent immune response is key to control of M. tuberculosis infection.8 Although the adaptive immune response induces development of solid granulomas, which contain the bacteria,9 Correspondence: Dr J Maertzdorf and Professor SHE Kaufmann, Department of Immunology, Max Planck Institute for Infection Biology, Chariteplatz 1, Berlin 10117, Germany. E-mail: [email protected] or [email protected] Received 18 May 2010; revised and accepted 12 August 2010; published online 23 September 2010

control of latent M. tuberculosis infection also requires robust innate immunity. The precise mechanisms by which these innate processes contribute to the control of latent infection remain elusive.10,11 In part of the infected population, latent infection progresses to active disease. Once immune control of infection wanes, granulomas liquefy and containment breaks down. As a result, actively replicating M. tuberculosis are released into the airways, rendering the disease contagious.12 Sensing of M. tuberculosis by macrophages and dendritic cells via pattern recognition receptors such as the Toll-like receptor family trigger signalling cascades, which initiate cytokine secretion, suggesting a role in protection against M. tuberculosis.13 However, the exact role of Toll-like receptor and other pattern recognition receptors in immunity against M. tuberculosis remains controversial.11,14 Bacterial killing is a key function of macrophages, which are activated by cytokines, notably tumour necrosis factor alpha and interferon gamma, which are mainly produced by T cells and natural killer (NK) cells. NK cells are wellknown key players in innate immunity, which actively induce death of M. tuberculosis-infected monocytes.15 It has been proposed to reduce TB incidence to o1 per 106 inhabitants by 2050.4 To reach this ambitious goal, better intervention measures are needed.7 It has been calculated that a combination of tools for rapid diagnosis, efficacious vaccines and novel drugs will be needed to reduce current TB incidence by 490%.16

Host gene expression profiles in tuberculosis J Maertzdorf et al

16

Biomarkers will be of critical importance for this approach as they can facilitate the identification of targets for new drugs and vaccines, allow better monitoring of clinical trials and provide guidelines for rapid diagnosis. In the present study, we evaluated transcriptional signatures in TB patients, healthy donors with latent M. tuberculosis infection (LTBI) and healthy non-infected donors (NIDs) in an attempt to define biomarkers predictive of susceptibility for, or resistance against, active TB.

Results Gene expression profiling in TB patients and household contacts (LTBI and NID) Genome-wide transcription profiles in whole blood from 33 TB patients, 34 LTBI and from 9 NIDs from a TBendemic area were generated using microarray chips comprising B45 000 unique features. Our strategy is aimed at analysing three separate comparisons to identify differentially expressed genes between study groups, that is, TB versus LTBI, TB versus NID and LTBI versus NIDs (Figure 1). Our ultimate goal is the identification of unique profiles that are associated with active TB. For our first goal, that is, identifying diseaserelated gene profiles, we focussed on comparisons of TB with LTBI as well as TB with NID. We proposed that these comparisons define patterns that distinguish between active TB disease and health. For our second goal, that is, identifying profiles that correlate with resistance and susceptibility, we focussed on genes that are differentially expressed in TB versus LTBI, but not in TB versus NID. This approach filtered out general disease-related patterns. A total of 2048 differentially expressed transcripts (representing 1935 genes) were identified between TB patients and LTBI. These included 988 transcripts (927 genes) that were also differentially expressed between TB patients and NIDs. In this analysis, a fold-change cutoff of 40.2 or o0.2 (log 2 scale) with qo0.01 was used (q-value equals the raw P-value corrected for multiple testing). A clustering analysis was performed to interrogate whether donors within this study could be divided into distinct groups based on their gene expression profiles. The resulting two-way clustered heat map, showing intensity levels in all donors, is depicted in Supplementary Figure 1. This analysis revealed clustering of most of the TB donors, whereas no obvious clusters of LTBI and NIDs were identified. Thus, there was a discriminative pattern between TB versus healthy individuals with or without M. tuberculosis infection, but no apparent difference between LTBI and NID.

Validation of differential gene expression by reverse transcriptase (RT)-PCR From the genes with highest degree of differential expression between TB and LTBI, 10 were selected to validate differential expression levels by quantitative RT-PCR. Beta-2-microglobulin was chosen as reference gene as this gene showed the least variation in expression within and between each study group (data not shown). A single sample from the TB group was not included because of insufficient amounts of RNA. Although most of the high-ranking genes showed distinct levels of expression between study groups, not all microarray results could be confirmed by RT-PCR (Supplementary Figure S2). The most strongly differentially expressed gene that was identified in TB versus control groups in this study was Fc gamma receptor 1B (FCGR1B). Expression validation by RT-PCR confirmed significant differences in expression of this gene between study groups (Figure 2a). An earlier study by our group had identified a minimal set of three genes (CD64, RAB33A and LTF (lactoferrin)), whose combined expression patterns showed a discriminating power between TB and LTBI in Caucasians.17 These candidate biomarker genes were also included in our RT-PCR validation to investigate whether this subset of genes could discriminate between TB and LTBI individuals from the South African cohort studied herein as well. Similar to those previous findings, CD64 and LTF showed a significant difference in gene expression between the study groups, whereas RAB33A expression levels were similar between TB and LTBI (Figure 2b). These findings underline differences and communalities in candidate TB biomarkers between study populations of different ethnicities. Universal biomarkers shared by a broad variety of ethnicities would be preferable. Random forest analysis was applied to calculate the discriminating power of these three genes (CD64, LTF and RAB33A) to distinguish between TB donors and LTBI. With CD64 being the most powerful discriminating gene, the combination of the three markers achieved a sensitivity and specificity of 88 and 91%, respectively, in identifying TB patients (Table 1a). Running a similar analysis on all 13 genes that had been subjected to RT-PCR expression evaluation, a combination of five most prominently differentiating genes turned out to give the highest accuracy in discriminating between TB and LTBI, that is, FCGR1B, CD64, LTF, guanylate binding protein 5 and Granzyme A. This intriguing subset of genes formed a biosignature capable of increasing sensitivity and specificity to 94% (30/32) and 97% (33/ 34), respectively (see Table 1b).

Figure 1 (a) Characteristics of study groups and group wise comparisons of gene expression profiles. LTBI, healthy donors latently infected with M. tuberculosis; NID, healthy non-infected donors; TB, tuberculosis patients. (b) Pairwise comparisons of gene expression between study groups. Genes and Immunity

Host gene expression profiles in tuberculosis J Maertzdorf et al

Table 1 Random forest analysis of gene subsets discriminating

17

between TB and LTBI

a

Predicted Gene subset LTF RAB33A

Actual

CD64

TB

LTBI

TB

28 (88%)

4 (12%)

LTBI

3 (9%)

31 (91%)

b

Predicted

FCGR1B CD64 LTF GBP5

Actual

Gene subset

TB

LTBI

TB

30 (94%)

2 (6%)

LTBI

1 (3%)

33 (97%)

compared with healthy individuals with (LTBI) or without (NID) M. tuberculosis infection.

Figure 2 Differences in gene expression levels between TB, LTBI and NID in the South African cohort studied in this report. (a) Scatter dot plot for FCGR1B expression for each individual sample, as PCR cycle differences in the internal reference gene beta-2microglobulin , indicating the mean and 95% confidence interval. (b) Box plots indicating expression variation of indicated genes for each study group. P value o0.01. Shown are the median with the boundaries representing 25th and 75th percentiles, and whiskers (error bars) above and below indicating the 90th and 10th percentiles.

Functional categorization of differentially expressed genes To categorize differentially regulated genes, we used the web-based tool Database for Annotation, Visualization and Integrated Discovery.18,19 Functional annotation clustering of genes common in comparisons of TB versus LTBI and TB versus NID revealed a high proportion of genes involved in protein binding, cell communication and signal transduction (Figure 3). Also, genes involved in distinct immune responses, apoptosis and stress responses were enriched in this gene ontology analysis (enrichment of gene ontology terms depicted in Figure 3, Po0.05). Genes differentially expressed in both comparisons behaved similarly, 792/988 transcripts (744/927 genes) showed increased or decreased expression in both TB versus LTBI and TB versus NID analyses. These genes, which showed highly significant enrichment in gene ontology terms, were mainly involved in immune response regulation and in apoptosis. Our results verify that systemic changes in immune responses are a useful indicator for active TB disease. The similarities in gene expression in both comparisons are likely to reflect a systemic inflammatory response in TB patients

Unique expression profiles in TB disease and LTBI In an attempt to identify patterns of gene expression indicative of progression to active TB disease from LTBI, we next focussed our analyses on genes differentially expressed in TB versus LTBI only. For this, we filtered out all genes from the TB versus LTBI comparison that were differentially expressed in TB versus NID, as well. This filtering revealed 1060 transcripts (representing 1008 genes) between TB and LTBI, but not between TB and NID. A comparative analysis led to the finding that a gene cluster is involved in apoptosis regulation with reduced expression in TB compared with LTBI (Table 2). Because genes present in both comparisons had been filtered out, this particular gene cluster showed similar levels of expression in TB and NID. Similarly, a group of genes involved in host defence responses and mainly active in granulocyte and macrophage effector function and differentiation was expressed at lower levels in LTBI compared with TB and NID (Table 2). These subsets of genes indicating decreased apoptotic activity and concomitantly increased innate host defence are apparently indicative for progression of TB disease from LTBI. Functional and cell-type-associated expression profiles Functional categorization of differentially expressed genes unique for TB versus LTBI revealed profound upregulation of Toll-like receptor-associated genes including MyD88 and IRF7. Simultaneously, the data reveal an increased expression of several interferoninducible genes, including guanylate binding protein 1 and 2, interferon alpha-inducible protein 6 and 27, ISG15 ubiquitin-like modifier and 20 ,50 -oligoadenylate synthetase 1 (see Table 3), which might illustrate the induction of interferon responses upon activation of pattern recognition receptors pathways.20,21 Intriguingly, a distinct gene expression profile between TB and LTBI was observed when analysing expression of well-defined subsets of macrophage- and NK cell-associated genes. Macrophage-associated genes were upregulated in TB Genes and Immunity

Host gene expression profiles in tuberculosis J Maertzdorf et al

18

Figure 3 Enriched biological properties and molecular functions of differentially expressed genes based on gene ontology (GO) classifications for genes shared by comparisons of TB versus LTBI and TB versus NID. Individual genes can be found under multiple GO annotations. Percentages indicate the number of genes sharing a certain GO term relative to the complete list of differentially expressed genes in both comparisons.

patients, whereas NK cell-associated genes were downregulated in TB patients as compared with LTBI patients (Table 3, list of cell-type-associated genes in part derived from Gaucher et al.22). Similarly, expression of complement- and inflammasome-associated genes were elevated in TB patients, whereas B cell-associated genes were, for the most part, decreased. It is to be noted that for all these functional and cell type-associated genes, gene expression levels did not differ between TB and NID. However, we did not perform a comparative phenotypic analysis of cellular composition of macrophages and NK cells in our samples. Thus, we cannot formally exclude that the observed differences in gene expression reflect divergent total cell counts of these cell types in LTBI patients compared with both TB patients and NIDs.

Discussion Although several key components involved in M. tuberculosis infection and active TB disease have been described, the mechanisms that determine successful resistance or susceptibility to TB disease remain illdefined. Definition of correlates of protection in TB will most likely require a multitude of biomarkers, termed biosignature.23,24 The complexity of the infection process that comprises LTBI and active TB disease necessitates a broader strategy aimed at identifying unique gene expression profiles. In the present study, we applied Genes and Immunity

whole-genome microarray expression analysis on wholeblood samples from TB patients and their controls (LTBI and NID) from a high TB-endemic area to identify such biomic profiles. Previous studies using microarray analysis in TB patients have identified gene patterns, which identify subjects at risk for recurrent disease among patients with cured TB25 and a distinct subset of genes to classify TB patients and LTBI in a Caucasian population.17 TB patients and healthy household control subjects (both latently infected and uninfected) included in this study were recruited from a cohort with high incidence of infection and disease. Accordingly, the recruitment of NID proved difficult, as almost all individuals in this cohort were found to give a positive result in a tuberculin skin test. FCGR1B (Fc gamma receptor 1B) was identified as the gene with the highest degree of differential expression between the three groups of donors in this study. An earlier study by our group had identified CD64 (alias FCGR1A; Fc gamma receptor 1A) as one of the highest ranking differentially expressed genes in TB-diseased versus latently infected Caucasian subjects.17 Despite of their different annotations, these genes have remarkably similar sequences. The oligonucleotide probes on the microarray chips used in both studies, were 100% complementary to the sequences of the probes used in either study. The microarray chips in both studies thus targeted transcripts from both genes. In contrast, in our RT-PCR evaluation, we designed primers to specifically target the CD64 encoding mRNA. Primers for FCGR1B

Host gene expression profiles in tuberculosis J Maertzdorf et al

Table 2

Table 3

Unique genes related to apoptosis and defence response differentially expressed between TB and LTBI

Functional and cell-type-associated gene expression in TB versus LTBI

Transcript ID

Gene symbol

Gene symbol

GOTERM: apoptosis NM_000024 NM_004208 NM_000633 NM_014739 NM_001165 NM_004052 NM_001731 NM_012115 NM_006139 NM_019604 NM_001316 NM_001008895 NM_001008540 NM_005494 NM_004944 NM_004836 NM_000123 NM_001992 NM_001033030 NM_000639 NM_006144 NM_004134 NM_014061 NM_006785 NM_002520 NM_153719 NM_024065 AK097958 NM_024900 NM_004703 NM_006265 NM_016542 NM_003031 NM_012238 NM_032195 NM_004619 NM_000553

ADRB2 AIFM1 BCL2 BCLAF1 BIRC3 BNIP3 BTG1 CASP8AP2 CD28 CRTAM CSE1L CUL4A CXCR4 DNAJB6 DNASE1L3 EIF2AK3 ERCC5 F2R FAIM FASLG GZMA HSPA9 MAGEH1 MALT1 NPM1 NUP62 PDCL3 PERP PHF17 RABEP1 RAD21 RP6-213H19.1 SIAH1 SIRT1 SON TRAF5 WRN

GOTERM: defence response NM_001637 NM_145641 NM_021160 NM_004331 NM_001725 NM_006077 NM_006678 NM_001250 NM_005194 NM_001805 NM_203339 NM_172313 NM_021175 NM_000634 NM_006863 NM_006864 NM_006840 NM_181657 NM_022468 NM_000265 NM_032206 NM_003005 NM_001002236 NM_015136 NM_003152 NM_003190 NM_006019 NM_005423

AOAH APOL3 BAT5 BNIP3 L BPI CBARA1 CD300C CD40 CEBPB CEBPE CLU CSF3R HAMP IL8RA LILRA1 LILRB3 LILRB5 LTB4R MMP25 NCF1 NLRC5 SELP SERPINA1 STAB1 STAT5A TAPBP TCIRG1 TFF2

M (log 2)

0.51 0.26 0.31 0.30 0.63 0.22 0.50 0.31 0.58 0.57 0.27 0.38 0.38 0.25 0.22 0.27 0.36 0.27 0.32 0.32 0.77 0.32 0.41 0.20 0.43 0.31 0.23 0.28 0.29 0.29 0.63 0.81 0.21 0.32 0.23 0.60 0.24 0.33 0.22 0.32 0.40 0.30 0.31 0.39 0.27 0.21 0.44 0.78 0.51 0.38 0.35 0.38 0.44 0.22 0.26 0.36 0.52 0.26 0.38 0.51 0.25 0.27 0.20 0.46 0.21

Abbreviations: LTBI, healthy donors latently infected with M. tuberculosis; TB, tuberculosis. Cutoff values used: M40.2 and o–0.2; qo0.01 (M ¼ fold change (log 2 scale); q equals the P-value corrected for multiple testing). Positive values (in bold) indicate higher expression; negative values (in italics) indicate lower expression in TB compared with LTBI.

M (log 2)

q

IFN-induced and TLR-associated genes AIM2 1.41 GBP1 1.56 GBP2 0.79 IFI27 0.63 IFI6 0.24 IFIT3 0.92 IRF7 0.51 ISG15 0.56 MYD88 0.35 OAS1 0.51 OASL 0.39 SPP1 0.23 TLR2 0.43 TLR4 0.30 TLR5 1.02 TLR6 0.32 TLR8 0.39 TRIM21 0.47 TRIM22 0.68 WARS 0.83

1.8E08 1.1E07 2.2E08 3.5E02 2.5E02 6.9E04 9.1E06 1.1E02 2.5E04 3.3E03 8.3E03 1.0E02 8.5E04 1.8E02 2.3E09 7.2E05 2.6E03 3.0E09 5.6E04 1.3E06

Complement-associated genes C1QA C1QB C3AR1 SERPING1

1.02990 1.35988 0.37221 1.13444

1.0E06 7.2E08 1.2E03 1.5E09

Inflammasome-associated genes CASP5 PYCARD IL1B

0.55520 0.27709 0.56420

9.5E06 8.5E03 5.0E04

0.34 0.23 0.26 0.85 0.79 0.25

1.7E02 1.9E02 4.9E02 6.3E05 8.1E04 3.7E05

0.26 0.45 0.43 0.89 0.25 0.72 to 0.25 0.57 0.45

4.4E06 4.0E02 1.5E02 9.5E08 4.3E02 o2.0E02 2.4E05 7.8E04

0.57388 0.82985 0.45887 0.57109

7.7E03 1.8E07 3.2E04 4.9E06

Macrophage-associated genes CD36 CD68 IFI16 MARCO NOD2 STAB1 NK cell-associated genes CD160 GNLY GZMB GZMK KIR2DL2 KLRsa NCR3 PRF1 B cell-associated genes POU2AF1 BANK1 CD19 PIK3AP1

19

Abbreviations: IFN, interferon; NK, natural killer; TLR, Toll-like receptor. M ¼ fold change (log2 scale). Positive values (in bold) indicate higher expression, negative values (in italics) indicate lower expression in TB compared with LTBI. q ¼ significance level (P-value corrected for multiple testing). a KLRs includes genes KLRB1, KLRC1, KLRC2, KLRC3, KLRD1, KLRF1, KLRK1.

Genes and Immunity

Host gene expression profiles in tuberculosis J Maertzdorf et al

20

were cross-reactive with CD64 owing to the high sequence similarity, thus allowing expression analysis of both genes. Although flow cytometry analyses revealed a higher percentage of monocytes—the main cell type expressing CD64—among peripheral blood mononuclear cells from TB donors compared with LTBI, the B15-fold difference in mRNA expression (estimate from RT-PCR analysis) cannot solely be explained by this slightly higher proportion of monocytes. Increased expression of the CD64 gene in monocytes from TB patients was also confirmed at the protein level in our previous study.17 A major function of the family of Fc receptors for IgG (FcgRs) is binding of antibodies by their constant domain. The role of antibodies has long been considered to be of minor importance in TB,26 although more recent studies indicate that antibodies may have a role in control against M. tuberculosis infections.27 In addition to their antibody binding capacity, the members within the FcgR family comprise a group of molecules that can simultaneously trigger activating and inhibitory signalling pathways to set thresholds for cell activation and thus generate a well-balanced immune response.28 The strong upregulation of the activating receptor FCGR1 in TB patients as identified here could reflect a crucial role of this molecule in the sustenance of chronic immune activation in these individuals, which might have a detrimental effect on disease outcome. Marked upregulation of FCGR1 in TB patients of different ethnicity warrants further elucidation of its precise role in TB, notably whether FCGR1 upregulation indicates an immune modulatory function, or an ‘attempt in vain’ of a protective response. The profiles defined in this study comprise a significant number of differentially expressed genes between individuals in all three study groups. These differences were most pronounced between diseased donors (TB) versus both groups of healthy controls (LTBI and NID), independent of M. tuberculosis infection. Functional annotation clustering of these genes indicated a high prevalence of factors involved in immune modulation, signal transduction, gene expression activity, as well as protein and metal ion binding factors. Additionally, distinct clusters of genes involved in apoptosis regulation and stress responses were differentially expressed between TB patients and both control groups, independent of infection with M. tuberculosis. Although these clusters can be considered to be disease susceptibility-related, these profiles may not be specific for TB. In general, we assume that many alterations in gene expression patterns reflect an overall activation and differentiation of a sustained inflammatory response. To define TB-specific susceptibility- and resistance-related genes within these clusters, we filtered out all genes that were differentially expressed in both TB versus LTBI and TB versus NID comparisons to exclude general pathology-related differences in genes expression. Differentially expressed gene profiles unique for TB versus LTBI comprised an intriguing subset of genes involved in apoptosis regulation, which were found to be less active in TB compared with LTBI. In addition, a group of gene clusters involved in defence responses of granulocytes and in macrophage effector functions and differentiation, was expressed at lower levels in LTBI. These expression levels were likely to be modulated by M. tuberculosis infection. Although bacterial

Genes and Immunity

infections often induce an increase in granulocyte numbers in the blood, this subset of genes showed no differences in expression levels between TB and NID. Thus, these patterns most likely indicate a downregulation in expression of these genes in LTBI. The most straightforward conclusion from our findings would rate increased expression of apoptosis-related genes and simultaneous downregulated expression of genes involved in granulocyte activation in LTBI beneficial for preventing outbreak of active disease. Alternatively, higher gene expression levels and activation of granulocytes and monocytes in TB patients could reflect mobilization of professional phagocytes during active TB in an attempt to strengthen defence in affected lung tissue despite the risk of heightened tissue damage. Control of M. tuberculosis infection depends on a delicate balance of immune control mainly by macrophages and T lymphocytes. Disturbance of this balance results in progression of latent infection to active TB disease. This imbalance in immune control apparently correlates with increased innate immunity, as indicated by elevated defence response gene expression, Toll-like receptor signalling and macrophage effector functions, as observed in TB patients in this study. Simultaneous reduced expression of particular subsets of genes involved in apoptosis and NK cell activity could affect eradication of replicating M. tuberculosis and thus have a key role in failure of the immune defence to successfully control latent infection. Reciprocally, fine-tuning of apoptosis and the innate immune defence responses appears critical for sustained control of M. tuberculosis in LTBI. The gene expression profiles identified in this study will need further validation by wide-scale gene expression analyses in different geographical and ethnical cohorts. Such studies are currently ongoing in cohorts from other African countries to validate current biomarker gene profiles and to identify TB-specific-gene expression signatures, as well as expression patterns common to TB and other chronic infectious diseases. More detailed analyses of expression profiles that may be crucial for protection and susceptibility in M. tuberculosisinfected individuals will thus be required to ultimately define a unique biosignature for TB.

Materials and methods Subject enrollment and sample collection The study presented here was approved by ethical committees in both Stellenbosch (South Africa) and Berlin (Germany) and written informed consent has been obtained from all study participants. From a cohort of TB patients and household contacts recruited at Stellenbosch University (South Africa), 33 TB patients, 34 healthy donors with LTBI and 9 healthy NIDs were included in this study (Figure 1a depicts the gender and age distribution of donors within groups). Household contacts are defined herein as immediate family members of a TB case (acid-fast bacilli, sputum-smear positive) sharing at least 3 h contact per day within the previous 3 months. TB patients and controls were characterized based on chest radiography, M. tuberculosis sputum culture results and tuberculin skin testing. All subjects were HIV– and samples from TB patients were taken before chemotherapy. From every donor, 2.5 ml of peripheral whole blood was collected in PAXgene tubes

Host gene expression profiles in tuberculosis J Maertzdorf et al

Table 4

21

qRT-PCR target genes and primer sequences

Gene name

Accession no.

Primer

Sequence (50 –30 )

B2M

NM_004048

CASP1

NM_033292

CD64

NM_000566

Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse Forward Reverse

GCTCGCGCTACTCTCTCTTT CTCTGCTGGATGACGTGAGT TACAGAGCTGGAGGCATTTG TTCCCGAATACCATGAGACA AGGCCTGGTTTGCAGCTTT CTGCCTCGCAGGGTCTTG TGTGCCAGGAAACAAGAGAG GACAGTGGACAAAGGTGTGG TGCAGGAAATGCAAAGAAAG AACTGGACCCTGTCGTTCTC TGGCAAATCCTACCTGATGA CCATATCCAAATTCCCTTGG CCAGAAGTGGTTGTTTCCCT TAGATGCAGGTTCCATAGCG GGAAGGGGTGCACCGGAAGG CACGGGGAGCAAGTGGGCAG GAAGAGACTCGTGCAATGGA AAGGCCAAAGGAAGTGACC GGTGACTTCCATCCAGGACT CCACATCTTATCGCATCCAC ACTGGGAGCAACTGGACTCT TAGAACCTCGCAATGACAGC AGAGAGACTCCCCCATCCAGT CCATCAAGGGTCACAGCATC AGATCCAGGTGCCCTCCAA GAGCATGTTGTGGGCATCAG GGGTTAACATTAGGCTGGGA CCGAACTGAGTATTGGTGGA

CEACAM1

NM_001024912

GBP1

NM_002053

GBP5

NM_052942

DEFA3

NM_005217

FCGR1B

NM_001017986

GZMA

NM_006144

HP

NM_005143

IFI27

NM_005532

LTF

NM_002343

Rab33A

NM_004794

S100A12

NM_005621

(PreAnalytix, Europe BD, Erembodegem, Belgium) and stored at 801C before processing. RNA extraction and microarray procedure RNA was extracted from 2.5 ml of peripheral whole blood collected in PAXgene tubes, using the PAXgene Blood RNA Kit following the manufacturer’s protocol (both tubes and kits from PreAnalytix). RNA concentration and integrity were determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Foster City, CA, USA). RNA labelling was done with the Fluorescent Linear Amplification Kit (Agilent Technologies) according to the manufacturer’s instructions. Quantity and labelling efficiency were verified before hybridization of the samples to whole-genome oligonucleotide microarray (Agilent Technologies). In this study, we used two-colour microarray slides on which Cy3- and Cy5- labelled RNA from each study group (TB, LTBI or NID) was co-hybridized with RNA from a gender- and age-matched individual from the other groups. Dye colours were balanced within group comparisons to avoid any dye-specific bias. Because of the lower number of NID available in this study, each RNA sample from this group was co-hybridized with a pool of RNA from three matched individuals from the other two groups. Data analysis Microarray chips were processed and scanned at 5 mm using an Agilent scanner. Image analysis was performed with Feature Extraction software (version 6.1.1, Agilent Technologies). Raw microarray data were processed using the function read.maimages of the Bioconductor29 R package limma30 and intensities of the spots not flagged out were background corrected using the method norm-exp.31 Background-corrected values were

Amplicon size 82 88 59 139 109 97 134 98 80 148 87 70 56 85

lowess normalized to obtain as unbiased red/green ratios as possible. For global comparability, the data of all arrays were quantile normalized.32 Further, features with mean log intensities o7 were also flagged out, as for these low-intensity spots random 1 noise in the data would lead to spurious signals. To statistically verify differential expression in the three pairwise comparisons of LTBI, NID and TB cases, we applied a combination of t-test with subsequent false discovery rate correction based on the approach of Efron and Tibshirani.33 Data were log transformed and differentially expressed genes were identified based on log-fold changes (M-values) in average gene expression with a q value o0.01 (q equals the P-value corrected for multiple testing) by two-sample t-test. Functional annotation analysis and clustering was performed with the help of the Database for Annotation, Visualization and Integrated Discovery Bioinformatics Resources 2008 (http:// david.abcc.ncifcrf.gov).18,19 Random forest analysis was done according to the method of Liaw and Wiener.34 Quantitative RT-PCR Differential expression of several genes was validated by quantitative RT-PCR (qRT-PCR). cDNA was generated by reverse transcription using oligo-dT primers and Superscript II (Invitrogen GmbH, Darmstadt, Germany). Transcripts were quantified by qRT-PCR based on SYBR Green incorporation (Applied Biosystems GmbH, Darmstadt, Germany) on an ABI PRISM 7900 thermocycler. Target genes and primer sequences are listed in Table 4.

Conflict of interest The authors declare no conflict of interest. Genes and Immunity

Host gene expression profiles in tuberculosis J Maertzdorf et al

22

Acknowledgements We thank ML Grossman for thoroughly revising the paper, J Weiner for biostatistical support and the Microarray Core Facility unit at the Max Planck Institute for Infection Biology in Berlin for sample processing. This study was funded by Grant number 37772 from the Bill & Melinda Gates Foundation through the Grand Challenges in Global Health Initiative, to all authors except DR

References 1 Smith NH, Hewinson RG, Kremer K, Brosch R, Gordon SV. Myths and misconceptions: the origin and evolution of Mycobacterium tuberculosis. Nat Rev Microbiol 2009; 7: 537–544. 2 Chan J, Flynn J. The immunological aspects of latency in tuberculosis. Clin Immunol 2004; 110: 2–12. 3 Kaufmann SH, Cole ST, Mizrahi V, Rubin E, Nathan C. Mycobacterium tuberculosis and the host response. J Exp Med 2005; 201: 1693–1697. 4 World Health Organization Annual Report 2009. Global Tuberculosis Control—Epidemiology, Strategy, Financing. WHO: Geneva, Switzerland, 2009. 5 Comas I, Gagneux S. The past and future of tuberculosis research. PLoS Pathog 2009; 5: e1000600. 6 Chaisson RE, Harrington M. How research can help control tuberculosis. Int J Tuberc Lung Dis 2009; 13: 558–568. 7 Marais BJ, Raviglione M, Donald PR, Harries AD, Kritski AL, Graham SM et al. Scale up services and research priorities for TB diagnosis, management and control. Call to Action. Lancet 2010; 375: 2179–2191. 8 Flynn FL, Chan J. Immunology of tuberculosis. Annu Rev Immunol 2001; 19: 93–129. 9 Ulrichs T, Kaufmann SH. New insights into the function of granulomas in human tuberculosis. J Pathol 2006; 208: 261–269. 10 Dorhoi A, Kaufmann SH. Fine-tuning of T cell responses during infection. Curr Opin Immunol 2009; 21: 367–377. 11 Korbel DS, Schneider BE, Schaible UE. Innate immunity in tuberculosis: myths and truth. Microb Infect 2008; 10: 995–1004. 12 Russell DG, Cardona P-J, Kim M-J, Allain S, Altare F. Foamy macrophages and the progression of the human tuberculosis granuloma. Nat Immunol 2009; 10: 943–948. 13 Quesniaux V, Fremond C, Jacobs M, Parida S, Nicolle D, Yeremeev V et al. Toll-like receptor pathways in the immune responses to mycobacteria. Microbes Infect 2004; 6: 946–959. 14 Kaufmann SH. Tuberculosis: deadly combination. Nature 2008; 15: 295–296. 15 Brill KJ, Li Q, Larkin R, Canaday DH, Kaplan DR, Boom WH et al. Human natural killer cells mediate killing of intracellular Mycobacterium tuberculosis H37Rv via granule-independent mechanisms. Infect Immun 2001; 69: 1755–1765.

16 World Health Organization. The Stop TB Strategy, Building on and Enhancing DOTS to Meet the TB-Related Millennium Development Goals. WHO: Geneva, 2006 (WHO/HTM/TB/ 2006.368). 17 Jacobsen M, Repsilber D, Gutschmidt A, Neher A, Feldmann K, Mollenkopf HJ et al. Candidate biomarkers for discrimination between infection and disease caused by Mycobacterium tuberculosis. J Mol Med 2007; 85: 613–621. 18 Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC et al. DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol 2003; 4: P3. 19 Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc 2009; 4: 44–57. 20 Brikos C, O’Neill LA. Signalling of toll-like receptors. Handb Exp Pharmacol 2008; 183: 21–50. 21 Takeuchi O, Akira S. Pattern recognition receptors and inflammation. Cell 2010; 140: 805–820. 22 Gaucher D, Therrien R, Kettaf N, Angermann BR, Boucher G, Filali-Mouhim A et al. Yellow fever vaccine induces integrated multilineage and polyfunctional immune responses. J Exp Med 2008; 205: 3119–3131. 23 Jacobsen M, Mattow J, Repsilber D, Kaufmann SH. Novel strategies to identify biomarkers in tuberculosis. Biol Chem 2008; 389: 487–495. 24 Parida SK, Kaufmann SH. The quest for biomarkers in tuberculosis. Drug Discov Today 2010; 15: 148–157. 25 Mistry R, Cliff JM, Clayton CL, Beyers N, Mohamed YS, Wilson PA et al. Gene-expression patterns in whole blood identify subjects at risk for recurrent tuberculosis. J Infect Dis 2007; 195: 357–365. 26 Glatman-Freedman A, Casadevall A. Serum therapy for tuberculosis revisited: reappraisal of the role of antibodymediated immunity against Mycobacterium tuberculosis. Clin Microbiol Rev 1998; 11: 514–532. 27 Abebe F, Bjune G. The protective role of antibody responses during Mycobacterium tuberculosis infection. Clin Exp Immunol 2009; 157: 235–243. 28 Nimmerjahn F, Ravetch JV. Fcg receptors as regulators of immune responses. Nat Rev Immunol 2008; 8: 34–47. 29 Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004; 5: R80. 30 Smyth GK, Michaud J, Scott HS. Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics 2005; 21: 2067–2075. 31 Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A et al. A comparison of background correction methods for two-colour microarrays. Bioinformatics 2007; 23: 2700–2707. 32 Smyth GK, Speed T. Normalization of cDNA microarray data. Methods 2003; 31: 265–273. 33 Efron B, Tibshirani R. Empirical Bayes methods and false discovery rates for microarrays. Genet Epidemiol 2002; 23: 70–86. 34 Liaw A, Wiener M. Classification and regression by randomforest. R News 2002; 2: 18–22.

Supplementary Information accompanies the paper on Genes and Immunity website (http://www.nature.com/gene)

Genes and Immunity