A predictive signature gene set for discriminating active from latent tuberculosis in Warao Amerindian children
© Verhagen et al; licensee BioMed Central Ltd. 2013
Received: 21 October 2012
Accepted: 29 January 2013
Published: 1 February 2013
Tuberculosis (TB) continues to cause a high toll of disease and death among children worldwide. The diagnosis of childhood TB is challenged by the paucibacillary nature of the disease and the difficulties in obtaining specimens. Whereas scientific and clinical research efforts to develop novel diagnostic tools have focused on TB in adults, childhood TB has been relatively neglected. Blood transcriptional profiling has improved our understanding of disease pathogenesis of adult TB and may offer future leads for diagnosis and treatment. No studies applying gene expression profiling of children with TB have been published so far.
We identified a 116-gene signature set that showed an average prediction error of 11% for TB vs. latent TB infection (LTBI) and for TB vs. LTBI vs. healthy controls (HC) in our dataset. A minimal gene set of only 9 genes showed the same prediction error of 11% for TB vs. LTBI in our dataset. Furthermore, this minimal set showed a significant discriminatory value for TB vs. LTBI for all previously published adult studies using whole blood gene expression, with average prediction errors between 17% and 23%. In order to identify a robust representative gene set that would perform well in populations of different genetic backgrounds, we selected ten genes that were highly discriminative between TB, LTBI and HC in all literature datasets as well as in our dataset. Functional annotation of these genes highlights a possible role for genes involved in calcium signaling and calcium metabolism as biomarkers for active TB. These ten genes were validated by quantitative real-time polymerase chain reaction in an additional cohort of 54 Warao Amerindian children with LTBI, HC and non-TB pneumonia. Decision tree analysis indicated that five of the ten genes were sufficient to classify 78% of the TB cases correctly with no LTBI subjects wrongly classified as TB (100% specificity).
Our data justify the further exploration of our signature set as biomarkers for potential childhood TB diagnosis. We show that, as the identification of different biomarkers in ethnically distinct cohorts is apparent, it is important to cross-validate newly identified markers in all available cohorts.
KeywordsBiomarker Children Mycobacterium tuberculosis Transcriptomics
It is estimated that one third of the world’s population is infected with Mycobacterium tuberculosis and that each year about nine million people develop tuberculosis (TB), one million (11%) of whom are children under 15 years of age . A unique aspect of TB in children is the rapid progression to disease, typically within the first year following infection, unlike in adults, where TB infection can persist for decades without progression into an active infection . Bacteriological confirmation in the diagnosis of childhood TB is the exception rather than the rule with culture remaining negative in around 70% of cases with probable TB . Using blood transcriptional profiling, several signature gene sets have been identified in adult cohorts from South Africa, The Gambia and The United Kingdom [4–6]. However, a significant overlap was shown with a biomarker set for sarcoidosis, suggesting the need for more specific biomarker sets . To statistically verify differential expression between active TB, latent TB infection (LTBI) and healthy controls (HC) different methods have been used, varying from statistical tests [4, 6] to prediction models using the k-nearest neighbours algorithm . Correlation analysis, a method selecting genes that are correlated with a single differentially expressed gene, was used to identify a biomarker set in a Gambian cohort . No studies applying gene expression profiling of children with TB have been published, and it is unknown whether the existing signature gene sets are applicable to childhood cohorts.
In Venezuela, a high TB incidence rate (3190 per 100,000) has been reported in Warao Amerindian children living in the Orinoco Delta in northeastern Venezuela . In this study, we identified new gene signatures in childhood TB by comparing gene expression profiles of Warao Amerindian children with TB, LTBI and HC. We validated the identified gene signatures from this study in an independent cohort of children with LTBI, HC or non-TB pneumonia. Furthermore, we estimated the predictive value of our gene signatures in previously performed adult studies and we compared the discriminatory power of the literature signature gene sets with our gene set.
Identification of signature genes
Characteristics of children with TB, LTBI and HC in which microarray analyses were performed
Number of donors
1.1 – 14.5
2.2 – 14.6
1.3 – 11.5
Class errors of the 116 signature gene set
Class error (%)
Validation of signature gene sets in independent sample sets
Performance of signature gene sets: cross prediction matrix showing prediction errors
Signature gene set
Berry et al.[]
Maertzdorf et al.[]
116 gene set
Minimal TB-LTBI set
Minimal TB-LTBI-HC set
Robust 10 gene q-PCR set
86 gene set
393 gene set
11 gene set
5 gene set
TB vs. LTBI
Berry et al.
Maertzdorf et al.
Maertzdorf et al.
average prediction error
TB vs. LTBI vs. HC
Berry et al.
Maertzdorf et al.
Maertzdorf et al.
TB vs. LTBI vs. HC vs. other disease
Berry et al.
Identification of the minimal discriminatory signature gene sets
A minimal discriminatory gene set to discriminate between TB and LTBI was identified by variable selection random forest (VarSelRF) . This procedure produced four genes for the set described in this manuscript, three genes for the dataset described by Berry et al., three genes for the dataset described by Maertzdorf et al. and two genes for the other dataset described by Maertzdorf et al. that are required to distinguish TB from LTBI samples (Additional file 3: Table S3). After removal of the overlaps, i.e. the genes that were selected in more than one dataset, nine genes remained (Additional file 3: Table S3). Similarly, a minimal gene set to distinguish TB from LTBI and from HC was identified. In this comparison, six genes were identified for the set described in this manuscript and respectively 35, three, and 25 genes were identified for the three other datasets [4–6] (Additional file 3: Table S3). After removal of the overlaps, 42 genes remained (Additional file 3: Table S3). Performance of these minimal sets was comparable to performance of the 116 gene set (Table 3). However, as only a few genes included in the minimal sets were present in more than one of the datasets, the selected minimal sets seem to be a summary of four small sets that perform well on each of the included datasets rather than a robust representative set that would perform well in ethnically different populations.
Identification and validation of the optimal signature gene set
Set of 10 signature genes with their role in TB, lung disease or inflammatory processes
Possible role in TB, lung disease or inflammatory processes
Cholinergic muscarine 2 receptor
cAMP regulation on airway smooth muscle.
· Loss of muscarine receptor function is associated with airway hyperreactivity .
Phagocytosis, clathrin-mediated endocytosis in alveolar macrophages .
· Inhibition of clathrin-mediated endocytosis led to inhibition of lipopolysaccharide (LPS) internalization and cytokine/chemokine release from macrophages stimulated by LPS .
Sorting nexin 17
Intracellular binding protein for the adhesion molecule P-selectin .
· P-selectin is important in the early phase of cell migration in TB infection and increased P-selectin serum levels are found in TB patients .
Phosphatidylinositol glycan anchor biosynthesis class C
Biosynthesis of glycosylphosphatidylinositol .
· Incorporation of the mycobacterial cell wall component lipoarabinomannan (LAM) into the macrophage cell membrane, a process that is dependent on successful insertion of a glycosylphosphatidylinositol anchor, is one of the key virulence factors for M. tuberculosis.
S100 calcium binding protein P
Calcium-binding protein involved in intracellular and extracellular calcium sensing and signal transduction .
· M. tuberculosis-mediated inhibition of a cytosolic rise in calcium is one of the essential steps in phagosome maturation .
Taste receptor type 2 member 46
Regulation of ciliary beat frequency through modulation of intracellular calcium concentration .
· M. tuberculosis-mediated inhibition of a cytosolic rise in calcium is one of the essential steps in phagosome maturation . Decreased expression of TAS2R receptors has been shown to lead to a decrease in intracellular calcium concentration .
· TAS2Rs are expressed on human airway smooth muscle where they cause bronchodilation through a localized calcium response .
Inhibition of formation of stress granules.
· Stress granules are host RNA cytoplasmic granules formed in response to infections by a pathway involving phosphorylation of the translation initiation factor eIF2α .
Encodes for the delta globin chain of HbA2.
· Involved in oxygen transport from the lung to the peripheral tissues.
Glycine dehydrogenase (decarboxylating)
Metabolic enzyme promoting cellular transformation.
· Altered GLDC expression has been correlated with survival time in lung cancer patients .
Acyl-CoA thioesterase 7
Expressed in macrophages, plays a role in inflammation through production of arachidonic acid.
· The molecular and cellular functions of ACOT7 have identified the enzyme as a candidate drug target in inflammatory diseases .
Receiver operating characteristic analysis of selected genes
Down- or upregulation in TB
TB (n = 9) vs. LTBI (n = 29)
TB (n = 9) vs. HC (n = 25)
TB (n = 9) vs. pneumonia (n = 18)
Although peripheral blood transcriptional signatures discriminating between TB, LTBI and HC subjects have been identified in adult studies [4–6], concerns about the specificity of these signature sets have been raised . Furthermore, the performance of these signature sets in children, that show high rates of progressive tuberculosis due to immaturity of the immune response, has not been investigated so far. In this study, we identified a 116 signature gene set that discriminated TB from LTBI and HC with class errors of 11%, 22%, and 0% for the respective classes of TB, LTBI and HC (Table 2). While this 116 gene signature set also showed a good discriminative value between TB and LTBI in adults from South Africa, The Gambia and The United Kingdom, signature sets that were identified in those adult cohorts were unable to discriminate TB from LTBI in our childhood cohort (Table 3).
Gene clusters that were enriched in our signature set included genes in the categories of (programmed) cell death and calcium binding (Additional file 2: Table S2). Both the Gambian as well as the South-African study of Maertzdorf et al.[5, 6] also described enrichment of genes involved in cell death. Other similarities between the functional annotations in the South African study  and our study are the enrichment of genes involved in regulation of cell proliferation, regulation of caspase activity and protein kinase activity. Specifically, CD64 was identified as the most powerful discriminating gene seperating TB from LTBI cases in the South African study . As CD64 has also been identified as a marker for general innate immune response activity and sepsis, this marker may not be specific to TB . Berry et al. observed that genes downstream of type I interferon-αβ receptor signaling were over-respresented in patients with active TB . However, type I inferferon signaling is also induced in response to respiratory viruses  and Streptococcus pneumoniae, questioning the specificity of genes involved in type I interferon receptor signaling as biomarkers for active TB.
The enrichment of genes involved in calcium signaling in our TB biomarker set has not been described before in adult studies using whole-blood gene expression [4–6], nor in studies based on transcriptional profiling of peripheral blood mononuclear cells (PBMCs) [31, 32]. A close relation between abnormal calcium metabolism and radiological extent of disease has been described in pulmonary TB patients [33, 34]. Alterations in serum calcium, particularly cases of hypercalcemia, have been observed in adult TB patients [33–35]. Hypercalcemia in pediatric TB patients is an infrequently recognized and poorly understood phenomenon . In lung tissue, several processes related to calcium homeostasis are thought to contribute to M. tuberculosis persistence and the aggregation of macrophages in granulomas. Over-production of 1,2-dihydroxyvitamin D3, which plays a traditional role in calcium metabolism, in alveolar macrophages in granulomas has a protective effect against oxidative injuries from the nitric oxide burst from granulomatous macrophages [37–39]. Furthermore, M. tuberculosis inhibits a calcium-dependent phagolysosome formation pathway which leads to the prevention of maturation of M. tuberculosis-containing phagosomes into phagolysosomes. This process, referred to as the M. tuberculosis phagosome maturation arrest, is critical for M. tuberculosis persistence in the human host . S100P, which significantly discriminated TB from LTBI in our study children, and TAS2R46, which significantly distinguished TB cases from HC, are genes involved in calcium signaling [21, 23, 24]. Possibly, altered expression of these genes in TB patients reflects M. tuberculosis-mediated changes in calcium metabolism in lung tissue that can be measured in peripheral whole blood.
Although the groups of TB, LTBI and HC were reasonably well age-matched (Table 1), we cannot exclude the possibility that age-dependent differences in immune responses have influenced gene expression profiles. Age-related differences in both innate as well as antigen-specific responses to M. tuberculosis are well recognised [40, 41]. Alveolar macrophage antimicrobial activity and recruitment of monocytes as well as the production of cytokines and certain aspects of antigen presentation appear to be less efficient in young children. This is particularly true in children younger than one year of age . Therefore, the exclusion of children less than one year of age in our study is likely to have prevented a significant influence of age-related immune differences on gene expression results. Furthermore, the signature set that we identified showed a good discriminative value between TB and LTBI in adults from several regions [4–6]. This is an indication that the genes that were selected in our analysis make up a signature set that performs well in individuals of all ages.
We identified a minimal gene set of 42 genes that was able to separate TB cases from LTBI and HC in all previously described (adult) cohorts [4–6] as well as in our childhood cohort. However, as this minimal set was possibly over-optimized to fit exactly those sets that were used for its composition, this set might not perform well in a newly identified cohort from a different geographic region. As the datasets used for the composition of the minimal set were based on European, African and South American populations the minimal set may not be applicable to individuals from Asia, while this region carries almost two-third of the global TB burden . Furthermore, this signature set could be only indicative of damage to the lung epithelium, similar to what has been described for the overlap of the gene set determined by Berry et al. and the biosignature characteristic for sarcoidosis . Therefore, we used bootstrapping procedures to select a robust set of ten genes that had a high discriminative value in our population, in the two populations described by Maertzdorf et al.[5, 6] and in the comparison between TB, LTBI, HC and other inflammatory and infectious diseases in the dataset of Berry et al.. Although this approach probably leads to less overfitting of the selected set towards the source databases used and less overlap with other infectious diseases in comparison with the minimal gene set we identified, the discriminatory power of this ten gene set is less than that of the minimal set (Table 3). Future cohorts can be of help in the reduction of the 116 gene set to a dataset with similar performance in discriminating TB from LTBI, HC and other inflammatory diseases as the minimal gene set without overfitting the dataset to the source datasets.
From the ten gene set, a combination of five (S100P, HBD, PIGC, CHRM2 and ACOT7) could be used in decision tree analysis to differentiate TB from LTBI, HC and non-TB pneumonia with 78% sensitivity and 96% specificity in our dataset (Figure 2). Additionally, the expression profile of children that were treated for TB shifted from an active TB classification (oval in Figure 2) towards a classification as not suffering from active TB (hexagon in Figure 2) at five months post treatment initiation. This indicates that these biomarkers reflect a dynamic response that changes as mycobactericidal activity diminishes.
The discriminatory value of the 116 gene signature set for the classification of cases in the cohort described by Berry and colleagues  was significantly better in people from London compared to people from South Africa. An explanation for the greater similarity between our study population with people from London than with people from South Africa comes from population-genetic studies in which a decrease in the level of genetic variation between populations is observed with increasing geographic distance from Africa, consistent with the out-of-Africa spread of human populations . The finding that previously published signature sets based on individuals from South Africa [4, 6] do not provide a good discriminatory value between TB, LTBI and HC in The Gambia  points towards a high heterozygosity in TB immune response between different African countries. A high-resolution survey of genotype variation based on single-nucleotide polymorphisms, copy-number variants and haplotype analysis of a worldwide sample of 29 populations revealed that the genetic distance between individuals from Asia and Native American or Colombian individuals is significantly less than the genetic distance between Asian and South African populations . Bayesian cluster analysis clustered individuals from East Asia together with Native American or Colombian individuals, indicating their close phylogenetic relationship . Clustering of Native American individuals with Asian individuals based on their genetic similarities was also observed in a recently published quantitative assessment of human genetic variation worldwide . Therefore, we speculate that the applicability of our signature set in Asian populations might be better than the applicability of sets identified in African or European populations.
This study provides a signature gene set that was demonstrated to be instrumental for the diagnosis of childhood TB. As the identification of different biomarkers in ethnically distinct cohorts is apparent, it is important to cross-validate newly identified markers in all available cohorts. Especially, more childhood cohorts should be investigated as TB diagnosis based on traditional methods is less sensitive and specific in children compared to adults.
The Warao Amerindians are an indigenous population living in wooden houses raised on stilts along the Orinoco river banks. With a population of approximately 30,000, the Warao people are the second most important Native American group in Venezuela. In this study, 27 HIV-negative children 1 to 15 years of age with TB (n = 9), LTBI (n = 9) and HC (n = 9) were recruited between May 2010 and December 2010. General characteristics of the study subjects are given in Table 1. Detailed information of the study subjects that was recorded on inclusion is given in Additional file 1: Table S1. Tuberculin skin test (TST) and QuantiFERON-TB Gold In-Tube assay (QFT-GIT) were performed in all children. Sputum samples were collected from all children with expectoration and gastric aspirates were taken from all children under 6 years of age. Children with active TB were diagnosed based on culture of M. tuberculosis (n = 2) or on the basis of clinical, epidemiological and radiological features (n = 7). The latter group were children with a TST ≥ 10 mm or a positive QFT-GIT result who presented all of: persistent fever > 38°C recorded daily for at least two weeks, persistent cough for more than three weeks, weight loss (> 5% reduction in weight compared with the highest weight recorded in last three months) or failure to thrive (documented crossing of percentile lines in the preceding three months), persistent lethargy or decrease in playfulness/activity reported by the parent and absence of clinical response on broad-spectrum antibiotics. Standard antero-posterior and lateral chest radiographs (CXRs) were taken from all children. Two independent experts, blinded to all clinical information, evaluated the CXRs and documented their findings on a standard report form. Where the two objective experts disagreed, a third expert was consulted and final consensus was achieved. A diagnosis of TB was only made when the CXR was consistent with TB  and the child showed a positive clinical response to anti-TB treatment. Children were followed up clinically, radiologically and, in case of a negative TST at inclusion, by means of TST at six and 12 months after inclusion. LTBI was defined as a TST ≥ 10 mm and a positive QFT-GIT with a negative culture result on inclusion in the absence of radiological and clinical evidence of TB disease on inclusion as well as on six and 12 months after inclusion. HC were children with a TST = 0 mm on inclusion and on six and 12 months after inclusion. The HC had a negative QFT-GIT and a negative culture result on inclusion in the absence of radiological and clinical evidence of TB disease on inclusion as well as on six and 12 months after inclusion. TB patients were sampled before initiation of anti-TB treatment. Of three of the nine TB patients, a follow-up sample was taken when the patient had taken anti-TB treatment for five months.
RNA isolation and microarray procedures
From every child, 2.5 ml of peripheral whole blood was collected in PAXgene RNA tubes (PreAnalytix, Hombrechtikon, Switzerland) and stored at −80°C prior to processing. RNA was isolated using the PAXgene Blood RNA kit (PreAnalytix) following the manufacturer’s protocol. RNA quality was checked using the RNA 6000 Nano Kit on a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA). Gene expression data were obtained using GeneChip Human Exon 1.0 ST Arrays (Affymetrix, Santa Clara, CA). The probe preparation and hybridization was done following Affymetrix protocol. Arrays were washed and stained according to the protocol on a GeneChip Fluidics Station 450 (Affymetrix) and scanned on a Genechip Scanner 3000 + autoloader (Affymetrix). Scan results were converted to CEL with the Affymetric scanning software. Spot intensity data were processed with Affymetrix Power Tools (version 1.14-4.1.1) using robust multi-array average (RMA) and plier-gcbg normalization on the core gene set (version HuEx-1 0-st-v2.r2.dt1.hg18). Gene expression values were log2-transformed and differentially expressed genes were identified based on log2 fold changes (M-values). P-values were calculated with a Bayes-regularized one-way ANOVA  followed by multiple testing correction of the p-values (q-values) according to the method of Storey and Tibshirani . Microarray data have been deposited in the GEO database under accession GSE41055 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=ltmhxwsmskeyyte&acc=GSE41055).
Random forest based identification of TB biomarker genes
A random forest classification was performed to identify signature genes for discrimination of TB, LTBI and HC. This classification model, consisting of 5000 decision trees trained on random subsets of samples and variables was trained based on log2 gene expression data as a function of individuals belonging to either TB, LTBI or HC classes. Irrelevant genes were removed using the random forest-based local importance measure as described in PhenoLink  where genes with a negative or neutral contribution of correctly classifying the samples were removed. This process was repeated until fewer than three genes could be removed per iteration. Next, genes were selected provided that they had a contribution to the correct classification of at least three samples of the same class. Again, this process was repeated until fewer than three genes were removed. The resulting set was used for classification and determination of the classification error, both as OOB error and as average error, on all gene expression sets.
Literature datasets describing gene expression studies of TB, LTBI, HC and other infectious diseases were obtained from GEO (http://nar.oxfordjournals.org/content/39/suppl_1/D1005.full), accession numbers GSE19491 , GSE19492 , GSE28623  and GSE25534 . For GSE19491, GSE19492 and GSE28623 available normalized data were used, while for GSE25534 the 2 dye array data was normalized using RMA in Arraystar (DNAStar, Madison, WI) to allow inter-slide comparison of gene expression data. First, the TB biomarker set identified in this manuscript was applied to the literature cohorts [4–6] and the biomarker sets identified in the literature cohorts were applied to each other for estimation of the predictive value of each identified set in the other populations. Next, to determine the minimal discriminatory gene sets based on the dataset described in this study as well as on the three previously published datasets [4–6], expression data from probes corresponding with the TB biomarker set determined in this manuscript were subjected to VarSelRF, as described by Diaz-Uriarte et al., which progressively eliminates genes with the lowest random forest-based local importance measure until no further improvements in the OOB error rate are reported. This procedure is designed to identify small, non-redundant sets of genes that have good predictive performance. However, over fitting of the selected minimum discriminatory gene sets on the source databases used could lead to a poor performance of these minimal sets when applying them to classify other patient populations. Therefore, .632+ bootstrapping procedures were performed to select genes with a high discriminative power in all datasets. Bootstrapping repeatedly analyzes subsamples rather than subsets of the data. As each subsample is a random sample with replacement from the full sample, this procedure performs well in the selection of classifiers even when there is over fitting .
To test possible associations between the class probability of being wrongly classified as LTBI (dependent variable) and age, gender and geographical region (independent variables) of TB cases included in the GSE19491/GSE19492 dataset, multivariable linear regression analysis was performed.
Enrichment analysis of signature genes was performed using the web-based DAVID bioinformatics tool (http://david.abcc.ncifcrf.gov) . Q-values for enriched annnotations were determined by a modified Fisher’s exact test (EASE Score ) and corrected for multiple testing by the Benjamini-Hochberg approach.
Quantitative RT-PCR on microarray and additional samples
Differential expression of several genes was validated by qRT-PCR. cDNA was generated by reverse transcription using Superscript III Reverse Transcriptase and Random Primers (Invitrogen) following manufacturer’s protocol after DNAse treatment of the RNA using TURBO DNA-free (Ambion). GAPDH was chosen as reference gene. qRT-PCR was performed on the 27 microarray samples as well as on a validation cohort (n = 54). These were RNA samples collected from an additional 54 children during the same period and in the same manner as the microarray samples. Of the 54 children of which RNA was collected for qRT-PCR, 20 were diagnosed with LTBI, 16 were HC and 18 children were diagnosed with a radiologically confirmed pneumonia , of which six were LTBI and 12 were HC.
ROC methodology was applied to evaluate the discriminatory ability of signature genes. All expression values were normally distributed (Kolmogorov-Smirnov’s test, p > 0.05). One-way analysis of variance and unpaired Student’s t tests were performed to compare expression values of these genes in children with TB, LTBI, HC and non-TB pneumonia. When the variances across groups were not equal (Levene’s test p < 0.05), Welch correction for nonhomogeneity of variance was applied.
The nature and objectives of the study were explained to the parents of children in Spanish or were simultaneously translated to their native language by Spanish-Warao bilingual native interpreters. The study was approved by the ethical committee of the Instituto de Biomedicina, the Regional Health Services, and the Delta Amacuro Indigenous Health Office (Servicio de Atención y Orientación al Indígena). Children were enrolled if their parents or primary caregivers provided written informed consent. Illiterate parents or caregivers signed by means of a thumb print.
The authors thank Mercedes España, coordinator of the Venezuelan National TB Control Program, Maria Eugenia Melendez, head of the pediatric pulmonology department of the Caracas Children’s Hospital Â´J.M. de los RíosÂ´, and Heriberto Perez, pulmonologist at the Caracas Hospital Â´José Maria VargasÂ´, for evaluation of CXRs. Furthermore, we thank the participating families and the field workers involved in the recruitment of participants, in particular the medical students of the Escuela de Medicina ‘José Maria Vargas’ of the Universidad Central de Venezuela and the personnel of the Centros de Diagnostico Medico Cubano (CDIs) in Curiapo and Pedernales. We also thank the staff of the Laboratorio de Tuberculosis, Instituto de Biomedicina (Caracas, Venezuela) and the Laboratory of Pediatric Infectious Diseases, Radboud University Medical Centre (Nijmegen, The Netherlands) for their clinical assistance and technical support. The study was supported by FUNDAIM, Fundación para la Investigación en Micobacterias, Caracas, Venezuela. The ‘Stichting VSB Fonds’ (Utrecht, The Netherlands) is greatly acknowledged for providing a stipend to LV. AZ was supported by the European Commission FP7 Marie Curie IEF Action 274586.
- World Health Organization (WHO): Guidance for national tuberculosis programmes on the management of tuberculosis in children. 2006, Geneva, Switzerland: WHOGoogle Scholar
- Newton SM, Brent AJ, Anderson S, Whittaker E, Kampmann B: Paediatric tuberculosis. Lancet Infect Dis. 2008, 8: 498-510.View ArticlePubMedGoogle Scholar
- Zar HJ, Hanslo D, Apolles P, Swingler G, Hussey G: Induced sputum versus gastric lavage for microbiological confirmation of pulmonary tuberculosis in infants and young children: a prospective study. Lancet. 2005, 365: 130-134.View ArticlePubMedGoogle Scholar
- Berry MP, Graham CM, McNab FW, Xu Z, Bloch SA, Oni T, Wilkinson KA, Banchereau R, Skinner J, Wilkinson RJ, Quinn C, Blankenship D, Dhawan R, Cush JJ, Mejias A, Ramilo O, Kon OM, Pascual V, Banchereau J, Chaussabel D, O’Garra A: An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature. 2010, 466: 973-977.PubMed CentralView ArticlePubMedGoogle Scholar
- Maertzdorf J, Ota M, Repsilber D, Mollenkopf HJ, Weiner J, Hill PC, Kaufmann SH: Functional correlations of pathogenesis-driven gene expression signatures in tuberculosis. PLoS One. 2011, 6: e26938-PubMed CentralView ArticlePubMedGoogle Scholar
- Maertzdorf J, Repsilber D, Parida SK, Stanley K, Roberts T, Black G, Walzl G, Kaufmann SH: Human gene expression profiles of susceptibility and resistance in tuberculosis. Genes Immun. 2011, 12: 15-22.View ArticlePubMedGoogle Scholar
- Maertzdorf J, Weiner J, Mollenkopf HJ, Network T, Bauer T, Prasse A, Muller-Quernheim J, Kaufmann SH: Common patterns and disease-related signatures in tuberculosis and sarcoidosis. Proc Natl Acad Sci USA. 2012, 109: 7853-7858.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernández De Larrea C, Fañdino C, López D, Del Nogal B, Rodríguez N, Convit J, Araujo Z, De Waard JH: Tuberculosis en menores de 15 años en la población Warao de Venezuela. Invest Clin. 2002, 43: 35-48.PubMedGoogle Scholar
- Bayjanov JR, Molenaar D, Tzeneva V, Siezen RJ, Van Hijum SA: PhenoLink - a web-tool for linking phenotype to omics data for bacteria: application to gene-trait matching for Lactobacillus plantarum strains. BMC Genomics. 2012, 13: 170-PubMed CentralView ArticlePubMedGoogle Scholar
- Da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
- Diaz-Uriarte R: GeneSrF and varSelRF: a web-based tool and R package for gene selection and classification using random forest. BMC Bioinformatics. 2007, 8: 328-PubMed CentralView ArticlePubMedGoogle Scholar
- Coulson FR, Fryer AD: Muscarinic acetylcholine receptors and airway diseases. Pharmacol Ther. 2003, 98: 59-69.View ArticlePubMedGoogle Scholar
- Yamada H, Ohashi E, Abe T, Kusumi N, Li SA, Yoshida Y, Watanabe M, Tomizawa K, Kashiwakura Y, Kumon H, Matsui H, Takei K: Amphiphysin 1 is important for actin polymerization during phagocytosis. Mol Biol Cell. 2007, 18: 4669-4680.PubMed CentralView ArticlePubMedGoogle Scholar
- Lau C, Wang X, Song L, North M, Wiehler S, Proud D, Chow CW: Syk associates with clathrin and mediates phosphatidylinositol 3-kinase activation during human rhinovirus internalization. J Immunol. 2008, 180: 870-880.View ArticlePubMedGoogle Scholar
- Wang H, Jiang C: Influenza A virus H5N1 entry into host cells is through clathrin-dependent endocytosis. Sci China C Life Sci. 2009, 52: 464-469.View ArticlePubMedGoogle Scholar
- Wang Y, Yang Y, Liu X, Wang N, Cao H, Lu Y, Zhou H, Zheng J: Inhibition of clathrin/dynamin-dependent internalization interferes with LPS-mediated TRAM-TRIF-dependent signaling pathway. Cell Immunol. 2012, 274: 121-129.View ArticlePubMedGoogle Scholar
- Florian V, Schluter T, Bohnensack R: A new member of the sorting nexin family interacts with the C-terminus of P-selectin. Biochem Biophys Res Commun. 2001, 281: 1045-1050.View ArticlePubMedGoogle Scholar
- Mukae H, Ashitani J, Tokojima M, Ihi T, Kohno S, Matsukura S: Elevated levels of circulating adhesion molecules in patients with active pulmonary tuberculosis. Respirology. 2003, 8: 326-331.View ArticlePubMedGoogle Scholar
- Watanabe R, Inoue N, Westfall B, Taron CH, Orlean P, Takeda J, Kinoshita T: The first step of glycosylphosphatidylinositol biosynthesis is mediated by a complex of PIG-A, PIG-H, PIG-C and GPI1. EMBO J. 1998, 17: 877-885.PubMed CentralView ArticlePubMedGoogle Scholar
- Welin A, Winberg ME, Abdalla H, Sarndahl E, Rasmusson B, Stendahl O, Lerm M: Incorporation of Mycobacterium tuberculosis lipoarabinomannan into macrophage membrane rafts is a prerequisite for the phagosomal maturation block. Infect Immun. 2008, 76: 2882-2887.PubMed CentralView ArticlePubMedGoogle Scholar
- Austermann J, Nazmi AR, Muller-Tidow C, Gerke V: Characterization of the Ca2+ − regulated ezrin-S100P interaction and its role in tumor cell migration. J Biol Chem. 2008, 283: 29331-293340.PubMed CentralView ArticlePubMedGoogle Scholar
- Vergne I, Chua J, Singh SB, Deretic V: Cell biology of mycobacterium tuberculosis phagosome. Annu Rev Cell Dev Biol. 2004, 20: 367-394.View ArticlePubMedGoogle Scholar
- Shah AS, Ben-Shahar Y, Moninger TO, Kline JN, Welsh MJ: Motile cilia of human airway epithelia are chemosensory. Science. 2009, 325: 1131-1134.PubMed CentralView ArticlePubMedGoogle Scholar
- Deshpande DA, Wang WC, McIlmoyle EL, Robinett KS, Schillinger RM, An SS, Sham JS, Liggett SB: Bitter taste receptors on airway smooth muscle bronchodilate by localized calcium signaling and reverse obstruction. Nat Med. 2010, 16: 1299-1304.PubMed CentralView ArticlePubMedGoogle Scholar
- Lindquist ME, Lifland AW, Utley TJ, Santangelo PJ, Crowe JE: Respiratory syncytial virus induces host RNA stress granules to facilitate viral replication. J Virol. 2010, 84: 12274-12284.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang WC, Shyh-Chang N, Yang H, Rai A, Umashankar S, Ma S, Soh BS, Sun LL, Tai BC, Nga ME, Bhakoo KK, Jayapal SR, Nichane M, Yu Q, Ahmed DA, Tan C, Sing WP, Tam J, Thirugananam A, Noghabi MS, Pang YH, Ang HS, Mitchell W, Robson P, Kaldis P, Soo RA, Swarup S, Lim EH, Lim B: Glycine decarboxylase activity drives non-small cell lung cancer tumor-initiating cells and tumorigenesis. Cell. 2012, 148: 259-272.View ArticlePubMedGoogle Scholar
- Forwood JK, Thakur AS, Guncar G, Marfori M, Mouradov D, Meng W, Robinson J, Huber T, Kellie S, Martin JL, Hume DA, Kobe B: Structural basis for recruitment of tandem hotdog domains in acyl-CoA thioesterase 7 and its role in inflammation. Proc Natl Acad Sci USA. 2007, 104: 10382-10387.PubMed CentralView ArticlePubMedGoogle Scholar
- van der Meer W, Pickkers P, Scott CS, van der Hoeven JG, Gunnewiek JK: Hematological indices, inflammatory markers and neutrophil CD64 expression: comparative trends during experimental human endotoxemia. J Endotoxin Res. 2007, 13: 94-100.View ArticlePubMedGoogle Scholar
- Garcia-Sastre A, Biron CA: Type 1 interferons and the virus-host relationship: a lesson in detente. Science. 2006, 312: 879-882.View ArticlePubMedGoogle Scholar
- Parker D, Martin FJ, Soong G, Harfenist BS, Aguilar JL, Ratner AJ, Fitzgerald KA, Schindler C, Prince A: Streptococcus pneumoniae DNA initiates type I interferon signaling in the respiratory tract. MBio. 2011, 2: e00016-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Jacobsen M, Repsilber D, Gutschmidt A, Neher A, Feldmann K, Mollenkopf HJ, Ziegler A, Kaufmann SH: Candidate biomarkers for discrimination between infection and disease caused by Mycobacterium tuberculosis. J Mol Med (Berl). 2007, 85: 613-621.View ArticleGoogle Scholar
- Lu C, Wu J, Wang H, Wang S, Diao N, Wang F, Gao Y, Chen J, Shao L, Weng X, Zhang Y, Zhang W: Novel biomarkers distinguishing active tuberculosis from latent infection identified by gene expression profile of peripheral blood mononuclear cells. PLoS One. 2011, 6: e24290-PubMed CentralView ArticlePubMedGoogle Scholar
- Chan TY, Chan CH, Shek CC, Davies PD: Hypercalcemia in active pulmonary tuberculosis and its occurrence in relation to the radiographic extent of disease. Southeast Asian J Trop Med Public Health. 1992, 23: 702-704.PubMedGoogle Scholar
- Deniz O, Tozkoparan E, Yonem A, Ciftci F, Bozkanat E, Cakir E, Ozcan O, Narin Y, Bilgic H, Ekiz K, Demirci N: Low parathormone levels and hypercalcaemia in patients with pulmonary tuberculosis: relation to radiological extent of disease and tuberculin skin test. Int J Tuberc Lung Dis. 2005, 9: 317-321.PubMedGoogle Scholar
- Dosumu EA, Momoh JA: Hypercalcemia in patients with newly diagnosed tuberculosis in Abuja. Nigeria. Can Respir J. 2006, 13: 83-87.PubMedGoogle Scholar
- Payne HA, Menson E, Sharland M, Bryant PA: Symptomatic hypercalcaemia in paediatric tuberculosis. Eur Respir Rev. 2011, 20: 53-56.View ArticlePubMedGoogle Scholar
- Cadranel JL, Garabedian M, Milleron B, Guillozzo H, Valeyre D, Paillard F, Akoun G, Hance AJ: Vitamin D metabolism by alveolar immune cells in tuberculosis: correlation with calcium metabolism and clinical manifestations. Eur Respir J. 1994, 7: 1103-1110.PubMedGoogle Scholar
- Chang JM, Kuo MC, Kuo HT, Hwang SJ, Tsai JC, Chen HC, Lai YH: 1-alpha,25-Dihydroxyvitamin D3 regulates inducible nitric oxide synthase messenger RNA expression and nitric oxide release in macrophage-like RAW 264.7 cells. J Lab Clin Med. 2004, 143: 14-22.View ArticlePubMedGoogle Scholar
- Carlberg C, Campbell MJ: Vitamin D receptor signaling mechanisms: Integrated actions of a well-defined transcription factor. Steroids. 2012, online publication ahead of printGoogle Scholar
- Lewinsohn DA, Gennaro ML, Scholvinck L, Lewinsohn DM: Tuberculosis immunology in children: diagnostic and therapeutic challenges and opportunities. Int J Tuberc Lung Dis. 2004, 8: 658-674.PubMedGoogle Scholar
- Smith S, Jacobs RF, Wilson CB: Immunobiology of childhood tuberculosis: a window on the ontogeny of cellular immunity. J Pediatr. 1997, 131: 16-26.View ArticlePubMedGoogle Scholar
- World Health Organization (WHO): Global Tuberculosis Report 2012. 2012, Geneva, Switzerland: WHOGoogle Scholar
- DeGiorgio M, Jakobsson M, Rosenberg NA: Out of Africa: modern human origins special feature: explaining worldwide patterns of human genetic variation using a coalescent-based serial founder model of migration outward from Africa. Proc Natl Acad Sci USA. 2009, 106: 16057-16062.PubMed CentralView ArticlePubMedGoogle Scholar
- Jakobsson M, Scholz SW, Scheet P, Gibbs JR, VanLiere JM, Fung HC, Szpiech ZA, Degnan JH, Wang K, Guerreiro R, Bras JM, Schymick JC, Hernandez DG, Traynor BJ, Simon-Sanchez J, Matarin M, Britton A, van de Leemput J, Rafferty I, Bucan M, Cann HM, Hardy JA, Rosenberg NA, Singleton AB: Genotype, haplotype and copy-number variation in worldwide human populations. Nature. 2008, 451: 998-1003.View ArticlePubMedGoogle Scholar
- Wang C, Zollner S, Rosenberg NA: A Quantitative Comparison of the Similarity between Genes and Geography in Worldwide Human Populations. PLoS Genet. 2012, 8: e1002886-PubMed CentralView ArticlePubMedGoogle Scholar
- Marais BJ, Gie RP, Schaaf HS, Starke JR, Hesseling AC, Donald PR, Beyers N: A proposed radiological classification of childhood intra-thoracic tuberculosis. Pediatr Radiol. 2004, 34: 886-894.View ArticlePubMedGoogle Scholar
- Baldi P, Long AD: A Bayesian framework for the analysis of microarray expression data: regularized t -test and statistical inferences of gene changes. Bioinformatics. 2001, 17: 509-519.View ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445.PubMed CentralView ArticlePubMedGoogle Scholar
- Efron B, Tibshirani RJ: Improvements on cross-validation: the .632+ bootstrap method. J American Statistical Association. 1997, 92: 548-560.Google Scholar
- World Health Organization (WHO) Pneumonia Vaccine Trial Investigators Group: Standardization of interpretation of chest radiographs for the diagnosis of pneumonia in children. 2001, Geneva, Switzerland: WHOGoogle Scholar