Concise gene signature for point-of-care classification of tuberculosis

Computational reassessment of RNA-seq data reveals key genes in active tuberculosis

BACKGROUND

Tuberculosis is a serious life-threatening disease among the top global health challenges and rapid and effective diagnostic biomarkers are vital for early diagnosis especially given the increasing prevalence of multidrug resistance.

METHODS

Two human whole blood microarray datasets, GSE42826 and GSE42830 were retrieved from publicly available gene expression omnibus (GEO) database. Deregulated genes (DEGs) were identified using GEO2R online tool and Gene Ontology (GO), protein-protein interaction (PPI) network analysis was performed using Metascape and STRING databases. Significant genes (n = 8) were identified using T-test/ANOVA and Molecular Complex Detection (MCODE) score ≥10, which was validated in GSE34608 dataset. The diagnostic potential of three biomarkers was assessed using Area Under Curve (AUC) of Receiver Operating Characteristic (ROC) plot. The transcriptional levels of these genes were also examined in a separate dataset GSE31348, to monitor the patterns of variation during tuberculosis treatment.

RESULTS

A total of 62 common DEGs (57 upregulated, 7 downregulated genes) were identified in two discovery datasets. GO functions and pathway enrichment analysis shed light on the functional roles of these DEGs in immune response and type-II interferon signaling. The genes in Module-1 (n = 18) were linked to innate immune response, interferon-gamma signaling. The common genes (n = 8) were validated in GSE34608 dataset, that corroborates the results obtained from discovery sets. The gene expression levels demonstrated responsiveness to Mtb infection during anti-TB therapy in GSE31348 dataset. In GSE34608 dataset, the expression levels of three specific genes, GBP5, IFITM3, and EPSTI1, emerged as potential diagnostic makers. In combination, these genes scored remarkable diagnostic performance with 100% sensitivity and 89% specificity, resulting in an impressive Area Under Curve (AUC) of 0.958. However, GBP5 alone showed the highest AUC of 0.986 with 100% sensitivity and 89% specificity.

CONCLUSIONS

The study presents valuable insights into the critical gene network perturbed during tuberculosis. These genes are determinants for assessing the effectiveness of an anti-TB response and distinguishing between active TB and healthy individuals. GBP5, IFITM3 and EPSTI1 emerged as candidate core genes in TB and holds potential as novel molecular targets for the development of interventions in the treatment of TB.

Comparison of gene set scoring methods for reproducible evaluation of tuberculosis gene signatures

BACKGROUND

Blood-based transcriptional gene signatures for tuberculosis (TB) have been developed with potential use to diagnose disease. However, an unresolved issue is whether gene set enrichment analysis of the signature transcripts alone is sufficient for prediction and differentiation or whether it is necessary to use the original model created when the signature was derived. Intra-method comparison is complicated by the unavailability of original training data and missing details about the original trained model. To facilitate the utilization of these signatures in TB research, comparisons between gene set scoring methods cross-data validation of original model implementations are needed.

METHODS

We compared the performance of 19 TB gene signatures across 24 transcriptomic datasets using both rrebuilt original models and gene set scoring methods. Existing gene set scoring methods, including ssGSEA, GSVA, PLAGE, Singscore, and Zscore, were used as alternative approaches to obtain the profile scores. The area under the ROC curve (AUC) value was computed to measure performance. Correlation analysis and Wilcoxon paired tests were used to compare the performance of enrichment methods with the original models.

RESULTS

For many signatures, the predictions from gene set scoring methods were highly correlated and statistically equivalent to the results given by the original models. In some cases, PLAGE outperformed the original models when considering signatures' weighted mean AUC values and the AUC results within individual studies.

CONCLUSION

Gene set enrichment scoring of existing gene sets can distinguish patients with active TB disease from other clinical conditions with equivalent or improved accuracy compared to the original methods and models. These data justify using gene set scoring methods of published TB gene signatures for predicting TB risk and treatment outcomes, especially when original models are difficult to apply or implement.

Age and sex influence antibody profiles associated with tuberculosis progression

Antibody features vary with tuberculosis (TB) disease state. Whether clinical variables, such as age or sex, influence associations between Mycobacterium tuberculosis-specific antibody responses and disease state is not well explored. Here we profiled Mycobacterium tuberculosis-specific antibody responses in 140 TB-exposed South African individuals from the Adolescent Cohort Study. We identified distinct response features in individuals progressing to active TB from non-progressing, matched controls. A multivariate antibody score differentially associated with progression (SeroScore) identified progressors up to 2 years before TB diagnosis, earlier than that achieved with the RISK6 transcriptional signature of progression. We validated these antibody response features in the Grand Challenges 6-74 cohort. Both the SeroScore and RISK6 correlated better with risk of TB progression in adolescents compared with adults, and in males compared with females. This suggests that age and sex are important, underappreciated modifiers of antibody responses associated with TB progression.

From immunology to artificial intelligence: revolutionizing latent tuberculosis infection diagnosis with machine learning

Latent tuberculosis infection (LTBI) has become a major source of active tuberculosis (ATB). Although the tuberculin skin test and interferon-gamma release assay can be used to diagnose LTBI, these methods can only differentiate infected individuals from healthy ones but cannot discriminate between LTBI and ATB. Thus, the diagnosis of LTBI faces many challenges, such as the lack of effective biomarkers from Mycobacterium tuberculosis (MTB) for distinguishing LTBI, the low diagnostic efficacy of biomarkers derived from the human host, and the absence of a gold standard to differentiate between LTBI and ATB. Sputum culture, as the gold standard for diagnosing tuberculosis, is time-consuming and cannot distinguish between ATB and LTBI. In this article, we review the pathogenesis of MTB and the immune mechanisms of the host in LTBI, including the innate and adaptive immune responses, multiple immune evasion mechanisms of MTB, and epigenetic regulation. Based on this knowledge, we summarize the current status and challenges in diagnosing LTBI and present the application of machine learning (ML) in LTBI diagnosis, as well as the advantages and limitations of ML in this context. Finally, we discuss the future development directions of ML applied to LTBI diagnosis.

Next-Generation TB Vaccines: Progress, Challenges, and Prospects

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), is a prevalent global infectious disease and a leading cause of mortality worldwide. Currently, the only available vaccine for TB prevention is Bacillus Calmette-Guérin (BCG). However, BCG demonstrates limited efficacy, particularly in adults. Efforts to develop effective TB vaccines have been ongoing for nearly a century. In this review, we have examined the current obstacles in TB vaccine research and emphasized the significance of understanding the interaction mechanism between MTB and hosts in order to provide new avenues for research and establish a solid foundation for the development of novel vaccines. We have also assessed various TB vaccine candidates, including inactivated vaccines, attenuated live vaccines, subunit vaccines, viral vector vaccines, DNA vaccines, and the emerging mRNA vaccines as well as virus-like particle (VLP)-based vaccines, which are currently in preclinical stages or clinical trials. Furthermore, we have discussed the challenges and opportunities associated with developing different types of TB vaccines and outlined future directions for TB vaccine research, aiming to expedite the development of effective vaccines. This comprehensive review offers a summary of the progress made in the field of novel TB vaccines.

Gene signature discovery and systematic validation across diverse clinical cohorts for TB prognosis and response to treatment

While blood gene signatures have shown promise in tuberculosis (TB) diagnosis and treatment monitoring, most signatures derived from a single cohort may be insufficient to capture TB heterogeneity in populations and individuals. Here we report a new generalized approach combining a network-based meta-analysis with machine-learning modeling to leverage the power of heterogeneity among studies. The transcriptome datasets from 57 studies (37 TB and 20 viral infections) across demographics and TB disease states were used for gene signature discovery and model training and validation. The network-based meta-analysis identified a common 45-gene signature specific to active TB disease across studies. Two optimized random forest regression models, using the full or partial 45-gene signature, were then established to model the continuum from Mycobacterium tuberculosis infection to disease and treatment response. In model validation, using pooled multi-cohort datasets to mimic the real-world setting, the model provides robust predictive performance for incipient to active TB risk over a 2.5-year period with an AUROC of 0.85, 74.2% sensitivity, and 78.3% specificity, which approximates the minimum criteria (>75% sensitivity and >75% specificity) within the WHO target product profile for prediction of progression to TB. Moreover, the model strongly discriminates active TB from viral infection (AUROC 0.93, 95% CI 0.91-0.94). For treatment monitoring, the TB scores generated by the model statistically correlate with treatment responses over time and were predictive, even before treatment initiation, of standard treatment clinical outcomes. We demonstrate an end-to-end gene signature model development scheme that considers heterogeneity for TB risk estimation and treatment monitoring.

Human GBP1 Is Involved in the Repair of Damaged Phagosomes/Endolysosomes

Mouse guanylate-binding proteins (mGBPs) are recruited to various invasive pathogens, thereby conferring cell-autonomous immunity against these pathogens. However, whether and how human GBPs (hGBPs) target M. tuberculosis (Mtb) and L. monocytogenes (Lm) remains unclear. Here, we describe hGBPs association with intracellular Mtb and Lm, which was dependent on the ability of bacteria to induce disruption of phagosomal membranes. hGBP1 formed puncta structures which were recruited to ruptured endolysosomes. Furthermore, both GTP-binding and isoprenylation of hGBP1 were required for its puncta formation. hGBP1 was required for the recovery of endolysosomal integrity. In vitro lipid-binding assays demonstrated direct binding of hGBP1 to PI4P. Upon endolysosomal damage, hGBP1 was targeted to PI4P and PI(3,4)P2-positive endolysosomes in cells. Finally, live-cell imaging demonstrated that hGBP1 was recruited to damaged endolysosomes, and consequently mediated endolysosomal repair. In summary, we uncover a novel interferon-inducible mechanism in which hGBP1 contributes to the repair of damaged phagosomes/endolysosomes.

Advances in point-of-care genetic testing for personalized medicine applications

Breakthroughs within the fields of genomics and bioinformatics have enabled the identification of numerous genetic biomarkers that reflect an individual's disease susceptibility, disease progression, and therapy responsiveness. The personalized medicine paradigm capitalizes on these breakthroughs by utilizing an individual's genetic profile to guide treatment selection, dosing, and preventative care. However, integration of personalized medicine into routine clinical practice has been limited-in part-by a dearth of widely deployable, timely, and cost-effective genetic analysis tools. Fortunately, the last several decades have been characterized by tremendous progress with respect to the development of molecular point-of-care tests (POCTs). Advances in microfluidic technologies, accompanied by improvements and innovations in amplification methods, have opened new doors to health monitoring at the point-of-care. While many of these technologies were developed with rapid infectious disease diagnostics in mind, they are well-suited for deployment as genetic testing platforms for personalized medicine applications. In the coming years, we expect that these innovations in molecular POCT technology will play a critical role in enabling widespread adoption of personalized medicine methods. In this work, we review the current and emerging generations of point-of-care molecular testing platforms and assess their applicability toward accelerating the personalized medicine paradigm.

Latent tuberculosis and computational biology: A less-talked affair

Tuberculosis (TB) is a pervasive and devastating air-borne disease caused by the organisms belonging to the Mycobacterium tuberculosis (Mtb) complex. Currently, it is the global leader in infectious disease-related death in adults. The proclivity of TB to enter the latent state has become a significant impediment to the global effort to eradicate TB. Despite decades of research, latent tuberculosis (LTB) mechanisms remain poorly understood, making it difficult to develop efficient treatment methods. In this review, we seek to shed light on the current understanding of the mechanism of LTB, with an accentuation on the insights gained through computational biology. We have outlined various well-established computational biology components, such as omics, network-based techniques, mathematical modelling, artificial intelligence, and molecular docking, to disclose the crucial facets of LTB. Additionally, we highlighted important tools and software that may be used to conduct a variety of systems biology assessments. Finally, we conclude the article by addressing the possible future directions in this field, which might help a better understanding of LTB progression.

Systematic review of diagnostic and prognostic host blood transcriptomic signatures of tuberculosis disease in people living with HIV

Background HIV-associated tuberculosis (TB) has high mortality; however, current triage and prognostic tools offer poor sensitivity and specificity, respectively. We conducted a systematic review of diagnostic and prognostic host-blood transcriptomic signatures of TB in people living with HIV (PLHIV). Methods We systematically searched online databases for studies published in English between 1990-2020. Eligible studies included PLHIV of any age in test or validation cohorts, and used microbiological or composite reference standards for TB diagnosis. Inclusion was not restricted by setting or participant age. Study selection, quality appraisal using the QUADAS-2 tool, and data extraction were conducted independently by two reviewers. Thereafter, narrative synthesis of included studies, and comparison of signatures performance, was performed. Results We screened 1,580 records and included 12 studies evaluating 31 host-blood transcriptomic signatures in 10 test or validation cohorts of PLHIV that differentiated individuals with TB from those with HIV alone, latent Mycobacterium tuberculosis infection, or other diseases (OD). Two (2/10; 20%) cohorts were prospective (29 TB cases; 51 OD) and 8 (80%) case-control (353 TB cases; 606 controls) design. All cohorts (10/10) were recruited in Sub-Saharan Africa and 9/10 (90%) had a high risk of bias. Ten signatures (10/31; 32%) met minimum WHO Target Product Profile (TPP) criteria for TB triage tests. Only one study (1/12; 8%) evaluated prognostic performance of a transcriptomic signature for progression to TB in PLHIV, which did not meet the minimum WHO prognostic TPP. Conclusions Generalisability of reported findings is limited by few studies enrolling PLHIV, limited geographical diversity, and predominantly case-control design, which also introduces spectrum bias. New prospective cohort studies are needed that include PLHIV and are conducted in diverse settings. Further research exploring the effect of HIV clinical, virological, and immunological factors on diagnostic performance is necessary for development and implementation of TB transcriptomic signatures in PLHIV.

Comparison of gene set scoring methods for reproducible evaluation of multiple tuberculosis gene signatures

Rationale

Many blood-based transcriptional gene signatures for tuberculosis (TB) have been developed with potential use to diagnose disease, predict risk of progression from infection to disease, and monitor TB treatment outcomes. However, an unresolved issue is whether gene set enrichment analysis (GSEA) of the signature transcripts alone is sufficient for prediction and differentiation, or whether it is necessary to use the original statistical model created when the signature was derived. Intra-method comparison is complicated by the unavailability of original training data, missing details about the original trained model, and inadequate publicly-available software tools or source code implementing models. To facilitate these signatures' replicability and appropriate utilization in TB research, comprehensive comparisons between gene set scoring methods with cross-data validation of original model implementations are needed.

Objectives

We compared the performance of 19 TB gene signatures across 24 transcriptomic datasets using both re-rebuilt original models and gene set scoring methods to evaluate whether gene set scoring is a reasonable proxy to the performance of the original trained model. We have provided an open-access software implementation of the original models for all 19 signatures for future use.

Methods

We considered existing gene set scoring and machine learning methods, including ssGSEA, GSVA, PLAGE, Singscore, and Zscore, as alternative approaches to profile gene signature performance. The sample-size-weighted mean area under the curve (AUC) value was computed to measure each signature's performance across datasets. Correlation analysis and Wilcoxon paired tests were used to analyze the performance of enrichment methods with the original models.

Measurement and Main Results

For many signatures, the predictions from gene set scoring methods were highly correlated and statistically equivalent to the results given by the original diagnostic models. PLAGE outperformed all other gene scoring methods. In some cases, PLAGE outperformed the original models when considering signatures' weighted mean AUC values and the AUC results within individual studies.

Conclusion

Gene set enrichment scoring of existing blood-based biomarker gene sets can distinguish patients with active TB disease from latent TB infection and other clinical conditions with equivalent or improved accuracy compared to the original methods and models. These data justify using gene set scoring methods of published TB gene signatures for predicting TB risk and treatment outcomes, especially when original models are difficult to apply or implement.

Systematic review of diagnostic and prognostic host blood transcriptomic signatures of tuberculosis disease in people living with HIV

Background HIV-associated tuberculosis (TB) has high mortality; however, current triage and prognostic tools offer poor sensitivity and specificity, respectively. We conducted a systematic review of diagnostic and prognostic host-blood transcriptomic signatures of TB in people living with HIV (PLHIV). Methods We systematically searched online databases for studies published in English between 1990-2020. Eligible studies included PLHIV of any age in test or validation cohorts, and used microbiological or composite reference standards for TB diagnosis. Inclusion was not restricted by setting or participant age. Study selection, quality appraisal using the QUADAS-2 tool, and data extraction were conducted independently by two reviewers. Thereafter, narrative synthesis of included studies, and comparison of signatures performance, was performed. Results We screened 1,580 records and included 12 studies evaluating 31 host-blood transcriptomic signatures in 10 test or validation cohorts of PLHIV that differentiated individuals with TB from those with HIV alone, latent Mycobacterium tuberculosis infection, or other diseases (OD). Two (2/10; 20%) cohorts were prospective (29 TB cases; 51 OD) and 8 (80%) case-control (353 TB cases; 606 controls) design. All cohorts (10/10) were recruited in Sub-Saharan Africa and 9/10 (90%) had a high risk of bias. Ten signatures (10/31; 32%) met minimum WHO Target Product Profile (TPP) criteria for TB triage tests. Only one study (1/12; 8%) evaluated prognostic performance of a transcriptomic signature for progression to TB in PLHIV, which did not meet the minimum WHO prognostic TPP. Conclusions Generalisability of reported findings is limited by few studies enrolling PLHIV, limited geographical diversity, and predominantly case-control design, which also introduces spectrum bias. New prospective cohort studies are needed that include PLHIV and are conducted in diverse settings. Further research exploring the effect of HIV clinical, virological, and immunological factors on diagnostic performance is necessary for development and implementation of TB transcriptomic signatures in PLHIV.

Utility of a three-gene transcriptomic signature in the diagnosis of tuberculosis in a low-endemic hospital setting

BACKGROUND

Host transcriptomic blood signatures have demonstrated diagnostic potential for tuberculosis (TB), requiring further validation across different geographical settings. Discriminating TB from other diseases with similar clinical manifestations is crucial for the development of an accurate immunodiagnostic tool. In this exploratory cohort study, we evaluated the performance of potential blood-based transcriptomic signatures in distinguishing TB disease from non-TB lower respiratory tract infections in hospitalised patients in a TB low-endemic country.

METHOD

Quantitative real-time polymerase chain reaction qPCR) was used to evaluate 26 previously published genes in blood from 31 patients (14 TB and 17 lower respiratory tract infection cases) admitted to Oslo University Hospital in Norway. The diagnostic accuracies of differentially expressed genes were determined by receiver operating characteristic curves.

RESULTS

A significant difference (p < .01) in the age distribution was observed between patients with TB (mean age, 40 ± 15 years) and lower respiratory tract infection (mean age 59 ± 12 years). Following adjustment for age, ETV7, GBP1, GBP5, P2RY14 and BLK were significantly differentially expressed between patients with TB and those with LRI. A general discriminant analysis generated a three-gene signature (BAFT2, ETV7 and CD1C), which diagnosed TB with an area under the receiver operating characteristic curve (AUC) of 0.86 (95% CI, 0.69 - 1.00), sensitivity of 69.23% (95% CI, 38.57%-90.91%) and specificity of 94.12% (95% CI, 71.31%-99.85%).

CONCLUSION

The three-genes signature may have potential to improve diagnosis of TB in a hospitalised low-burden setting. However, the influence of confounding variables or covariates such as age requires further evaluation in larger studies.

Functions of exosomal non-coding RNAs to the infection with Mycobacterium tuberculosis

Tuberculosis (TB) is a major infectious disease induced by Mycobacterium tuberculosis (M. tb) which causes the world's dominant fatal bacterial contagious disease. Increasing studies have indicated that exosomes may be a novel option for the diagnosis and treatment of TB. Exosomes are nanovesicles (30-150 nm) containing lipids, proteins and non-coding RNAs (ncRNAs) released from various cells, and can transfer their cargos and communicate between cells. Furthermore, exosomal ncRNAs exhibit diagnosis potential in bacterial infections, including TB. Additionally, differential exosomal ncRNAs regulate the physiological and pathological functions of M. tb-infected cells and act as diagnostic markers for TB. This current review explored the potential biological roles and the diagnostic application prospects of exosomal ncRNAs, and included recent information on their pathogenic and therapeutic functions in TB.

A protein signature associated with active tuberculosis identified by plasma profiling and network-based analysis

Annually, approximately 10 million people are diagnosed with active tuberculosis (TB), and 1.4 million die of the disease. If left untreated, each person with active TB will infect 10-15 new individuals. The lack of non-sputum-based diagnostic tests leads to delayed diagnoses of active pulmonary TB cases, contributing to continued disease transmission. In this exploratory study, we aimed to identify biomarkers associated with active TB. We assessed the plasma levels of 92 proteins associated with inflammation in individuals with active TB (n = 20), latent TB (n = 14), or healthy controls (n = 10). Using co-expression network analysis, we identified one module of proteins with strong association with active TB. We removed proteins from the module that had low abundance or were associated with non-TB diseases in published transcriptomic datasets, resulting in a 12-protein plasma signature that was highly enriched in individuals with pulmonary and extrapulmonary TB and was further associated with disease severity.

Development and Validation of a Parsimonious Tuberculosis Gene Signature Using the digital NanoString nCounter Platform

BACKGROUND

Blood-based biomarkers for diagnosing active tuberculosis (TB), monitoring treatment response, and predicting risk of progression to TB disease have been reported. However, validation of the biomarkers across multiple independent cohorts is scarce. A robust platform to validate TB biomarkers in different populations with clinical end points is essential to the development of a point-of-care clinical test. NanoString nCounter technology is an amplification-free digital detection platform that directly measures mRNA transcripts with high specificity. Here, we determined whether NanoString could serve as a platform for extensive validation of candidate TB biomarkers.

METHODS

The NanoString platform was used for performance evaluation of existing TB gene signatures in a cohort in which signatures were previously evaluated on an RNA-seq dataset. A NanoString codeset that probes 107 genes comprising 12 TB signatures and 6 housekeeping genes (NS-TB107) was developed and applied to total RNA derived from whole blood samples of TB patients and individuals with latent TB infection (LTBI) from South India. The TBSignatureProfiler tool was used to score samples for each signature. An ensemble of machine learning algorithms was used to derive a parsimonious biomarker.

RESULTS

Gene signatures present in NS-TB107 had statistically significant discriminative power for segregating TB from LTBI. Further analysis of the data yielded a NanoString 6-gene set (NANO6) that when tested on 10 published datasets was highly diagnostic for active TB.

CONCLUSIONS

The NanoString nCounter system provides a robust platform for validating existing TB biomarkers and deriving a parsimonious gene signature with enhanced diagnostic performance.

Isoniazid preventive therapy and tuberculosis transcriptional signatures in people with HIV

OBJECTIVES

To examine the association between isoniazid preventive therapy (IPT) or nontuberculous mycobacteria (NTM) sputum culture positivity and tuberculosis (TB) transcriptional signatures in people with HIV.

DESIGN

Cross-sectional study.

METHODS

We enrolled adults living with HIV who were IPT-naive or had completed IPT more than 6 months prior at HIV care clinics in western Kenya. We calculated TB signatures using gene expression data from qRT-PCR. We used multivariable linear regression to analyze the association between prior receipt of IPT or NTM sputum culture positivity with a transcriptional TB risk score, RISK6 (range 0-1). In secondary analyses, we explored the association between IPT or NTM positivity and four other TB transcriptional signatures.

RESULTS

Among 381 participants, 99.7% were receiving antiretroviral therapy and 86.6% had received IPT (completed median of 1.1 years prior). RISK6 scores were lower (mean difference 0.10; 95% confidence interval (CI): 0.06-0.15; P  < 0.001) among participants who received IPT than those who did not. In a model that adjusted for age, sex, duration of ART, and plasma HIV RNA, the RISK6 score was 52.8% lower in those with a history of IPT ( P  < 0.001). No significant association between year of IPT receipt and RISK6 scores was detected. There was no association between NTM sputum culture positivity and RISK6 scores.

CONCLUSION

In people with HIV, IPT was associated with significantly lower RISK6 scores compared with persons who did not receive IPT. These data support investigations of its performance as a TB preventive therapy response biomarker.

Transcriptional profiles predict treatment outcome in patients with tuberculosis and diabetes at diagnosis and at two weeks after initiation of anti-tuberculosis treatment

BACKGROUND

Globally, the tuberculosis (TB) treatment success rate is approximately 85%, with treatment failure, relapse and death occurring in a significant proportion of pulmonary TB patients. Treatment success is lower among people with diabetes mellitus (DM). Predicting treatment outcome early after diagnosis, especially in TB-DM patients, would allow early treatment adaptation for individuals and may improve global TB control.

METHODS

Samples were collected in a longitudinal cohort study of adult TB patients from South Africa (n  =  94) and Indonesia (n  =  81), who had concomitant DM (n  =  59), intermediate hyperglycaemia (n  =  79) or normal glycaemia/no DM (n  =  37). Treatment outcome was monitored, and patients were categorized as having a good (cured) or poor (failed, recurrence, died) outcome during treatment and 12 months follow-up. Whole blood transcriptional profiles before, during and at the end of TB treatment were characterized using unbiased RNA-Seq and targeted gene dcRT-MLPA.

FINDINGS

We report differences in whole blood transcriptome profiles, which were observed before initiation of treatment and throughout treatment, between patients with a good versus poor TB treatment outcome. An eight-gene and a 22-gene blood transcriptional signature distinguished patients with a good TB treatment outcome from patients with a poor TB treatment outcome at diagnosis (AUC = 0·815) or two weeks (AUC = 0·834) after initiation of TB treatment, respectively. High accuracy was obtained by cross-validating this signature in an external cohort (AUC = 0·749).

INTERPRETATION

These findings suggest that transcriptional profiles can be used as a prognostic biomarker for treatment failure and success, even in patients with concomitant DM.

FUNDING

The research leading to these results, as part of the TANDEM Consortium, received funding from the European Community's Seventh Framework Programme (FP7/2007-2013 Grant Agreement No. 305279) and the Netherlands Organization for Scientific Research (NWO-TOP Grant Agreement No. 91214038). The research leading to the results presented in the Indian validation cohort was supported by Research Council of Norway Global Health and Vaccination Research (GLOBVAC) projects: RCN 179342, 192534, and 248042, the University of Bergen (Norway).

Systematic review of host genomic biomarkers of invasive bacterial disease: Distinguishing bacterial from non-bacterial causes of acute febrile illness

Background

Infectious diseases play a significant role in the global burden of disease. The gold standard for the diagnosis of bacterial infection, bacterial culture, can lead to diagnostic delays and inappropriate antibiotic use. The advent of high- throughput technologies has led to the discovery of host-based genomic biomarkers of infection, capable of differentiating bacterial from other causes of infection, but few have achieved validation for use in a clinical setting.

Methods

A systematic review was performed. PubMed/Ovid Medline, Ovid Embase and Scopus databases were searched for relevant studies from inception up to 30/03/2022 with forward and backward citation searching of key references. Studies assessing the diagnostic performance of human host genomic biomarkers of bacterial infection were included. Study selection and assessment of quality were conducted by two independent reviewers. A meta-analysis was undertaken using a diagnostic random-effects model. The review was registered with PROSPERO (ID: CRD42021208462).

Findings

Seventy-two studies evaluating the performance of 116 biomarkers in 16,216 patients were included. Forty-six studies examined TB-specific biomarker performance and twenty-four studies assessed biomarker performance in a paediatric population. The results of pooled sensitivity, specificity, negative and positive likelihood ratio, and diagnostic odds ratio of genomic biomarkers of bacterial infection were 0.80 (95% CI 0.78 to 0.82), 0.86 (95% CI 0.84 to 0.88), 0.18 (95% CI 0.16 to 0.21), 5.5 (95% CI 4.9 to 6.3), 30.1 (95% CI 24 to 37), respectively. Significant between-study heterogeneity (I² 77%) was present.

Interpretation

Host derived genomic biomarkers show significant potential for clinical use as diagnostic tests of bacterial infection however, further validation and attention to test platform is warranted before clinical implementation can be achieved.

Funding

No funding received.

Increasing reproducibility, robustness, and generalizability of biomarker selection from meta-analysis using Bayesian methodology

A major limitation of gene expression biomarker studies is that they are not reproducible as they simply do not generalize to larger, real-world, heterogeneous populations. Frequentist multi-cohort gene expression meta-analysis has been frequently used as a solution to this problem to identify biomarkers that are truly differentially expressed. However, the frequentist meta-analysis framework has its limitations-it needs at least 4-5 datasets with hundreds of samples, is prone to confounding from outliers and relies on multiple-hypothesis corrected p-values. To address these shortcomings, we have created a Bayesian meta-analysis framework for the analysis of gene expression data. Using real-world data from three different diseases, we show that the Bayesian method is more robust to outliers, creates more informative estimates of between-study heterogeneity, reduces the number of false positive and false negative biomarkers and selects more generalizable biomarkers with less data. We have compared the Bayesian framework to a previously published frequentist framework and have developed a publicly available R package for use.

Diagnosis of pulmonary tuberculosis via identification of core genes and pathways utilizing blood transcriptional signatures: a multicohort analysis

Background

Blood transcriptomics can be used for confirmation of tuberculosis diagnosis or sputumless triage, and a comparison of their practical diagnostic accuracy is needed to assess their usefulness. In this study, we investigated potential biomarkers to improve our understanding of the pathogenesis of active pulmonary tuberculosis (PTB) using bioinformatics methods.

Methods

Differentially expressed genes (DEGs) were analyzed between PTB and healthy controls (HCs) based on two microarray datasets. Pathways and functional annotation of DEGs were identified and ten hub genes were selected. They were further analyzed and selected, then verified with an independent sample set. Finally, their diagnostic power was further evaluated between PTB and HCs or other diseases.

Results

62 DEGs mostly related to type I IFN pathway, IFN-γ-mediated pathway, etc. in GO term and immune process, and especially RIG-I-like receptor pathway were acquired. Among them, OAS1, IFIT1 and IFIT3 were upregulated and were the main risk factors for predicting PTB, with adjusted risk ratios of 1.36, 3.10, and 1.32, respectively. These results further verified that peripheral blood mRNA expression levels of OAS1, IFIT1 and IFIT3 were significantly higher in PTB patients than HCs (all P < 0.01). The performance of a combination of these three genes (three-gene set) had exceeded that of all pairwise combinations of them in discriminating TB from HCs, with mean AUC reaching as high as 0.975 with a sensitivity of 94.4% and a specificity of 100%. The good discernibility capacity was evaluated d via 7 independent datasets with an AUC of 0.902, as well as mean sensitivity of 87.9% and mean specificity of 90.2%. In regards to discriminating PTB from other diseases (i.e., initially considered to be possible TB, but rejected in differential diagnosis), the three-gene set equally exhibited an overall strong ability to separate PTB from other diseases with an AUC of 0.999 (sensitivity: 99.0%; specificity: 100%) in the training set, and 0.974 with a sensitivity of 96.4% and a specificity of 98.6% in the test set.

Conclusion

The described commonalities and unique signatures in the blood profiles of PTB and the other control samples have considerable implications for PTB biosignature design and future diagnosis, and provide insights into the biological processes underlying PTB.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-022-02035-4.

Prospective multicentre head-to-head validation of host blood transcriptomic biomarkers for pulmonary tuberculosis by real-time PCR

Background

Sensitive point-of-care screening tests are urgently needed to identify individuals at highest risk of tuberculosis. We prospectively tested performance of host-blood transcriptomic tuberculosis signatures.

Methods

Adults without suspicion of tuberculosis were recruited from five endemic South African communities. Eight parsimonious host-blood transcriptomic tuberculosis signatures were measured by microfluidic RT-qPCR at enrolment. Upper respiratory swab specimens were tested with a multiplex bacterial-viral RT-qPCR panel in a subset of participants. Diagnostic and prognostic performance for microbiologically confirmed prevalent and incident pulmonary tuberculosis was tested in all participants at baseline and during active surveillance through 15 months follow-up, respectively.

Results

Among 20,207 HIV-uninfected and 963 HIV-infected adults screened; 2923 and 861 were enroled. There were 61 HIV-uninfected (weighted prevalence 1.1%) and 10 HIV-infected (prevalence 1.2%) tuberculosis cases at baseline. Parsimonious signature diagnostic performance was superior among symptomatic (AUCs 0.85-0.98) as compared to asymptomatic (AUCs 0.61-0.78) HIV-uninfected participants. Thereafter, 24 HIV-uninfected and 9 HIV-infected participants progressed to incident tuberculosis (1.1 and 1.0 per 100 person-years, respectively). Among HIV-uninfected individuals, prognostic performance for incident tuberculosis occurring within 6-12 months was higher relative to 15 months. 1000 HIV-uninfected participants were tested for respiratory microorganisms and 413 HIV-infected for HIV plasma viral load; 7/8 signature scores were higher (p < 0.05) in participants with viral respiratory infections or detectable HIV viraemia than those without.

Conclusions

Several parsimonious tuberculosis transcriptomic signatures met triage test targets among symptomatic participants, and incipient test targets within 6 months. However, the signatures were upregulated with viral infection and offered poor specificity for diagnosing sub-clinical tuberculosis.

Single-cell transcriptomic landscape identifies the expansion of peripheral blood monocytes as an indicator of HIV-1-TB co-infection

Certain circulating cell subsets are involved in immune dysregulation in human immunodeficiency virus type 1 (HIV-1) and tuberculosis (TB) co-infection; however, the characteristics and role of these subclusters are unknown. Peripheral blood mononuclear cells (PBMCs) of patients with HIV-1 infection alone (HIV-pre) and those with HIV-1-TB co-infection without anti-TB treatment (HIV-pre & TB-pre) and with anti-TB treatment for 2 weeks (HIV-pre & TB-pos) were subjected to single-cell RNA sequencing (scRNA-seq) to characterize the transcriptome of different immune cell subclusters. We obtained > 60,000 cells and identified 32 cell subclusters based on gene expression. The proportion of immune-cell subclusters was altered in HIV-1-TB co-infected individuals compared with that in HIV-pre-group, indicating immune dysregulation corresponding to different disease states. The proportion of an inflammatory CD14+CD16+ monocyte subset was higher in the HIV-pre & TB-pre group than in the HIV-pre group; this was validated in an additional cohort (n = 80) via a blood cell differential test, which also demonstrated a good discriminative performance (area under the curve, 0.8046). These findings depicted the atlas of immune PBMC subclusters in HIV-1-TB co-infection and demonstrate that monocyte subsets in peripheral blood might serve as a discriminating biomarker for diagnosis of HIV-1-TB co-infection.

Blood transcriptomics reveal the evolution and resolution of the immune response in tuberculosis

Blood transcriptomics have revealed major characteristics of the immune response in active TB, but the signature early after infection is unknown. In a unique clinically and temporally well-defined cohort of household contacts of active TB patients that progressed to TB, we define minimal changes in gene expression in incipient TB increasing in subclinical and clinical TB. While increasing with time, changes in gene expression were highest at 30 d before diagnosis, with heterogeneity in the response in household TB contacts and in a published cohort of TB progressors as they progressed to TB, at a bulk cohort level and in individual progressors. Blood signatures from patients before and during anti-TB treatment robustly monitored the treatment response distinguishing early and late responders. Blood transcriptomics thus reveal the evolution and resolution of the immune response in TB, which may help in clinical management of the disease.

Prediction of anti-tuberculosis treatment duration based on a 22-gene transcriptomic model

BACKGROUND

The World Health Organization recommends standardised treatment durations for patients with tuberculosis (TB). We identified and validated a host-RNA signature as a biomarker for individualised therapy durations for patients with drug-susceptible (DS)- and multidrug-resistant (MDR)-TB.

METHODS

Adult patients with pulmonary TB were prospectively enrolled into five independent cohorts in Germany and Romania. Clinical and microbiological data and whole blood for RNA transcriptomic analysis were collected at pre-defined time points throughout therapy. Treatment outcomes were ascertained by TBnet criteria (6-month culture status/1-year follow-up). A whole-blood RNA therapy-end model was developed in a multistep process involving a machine-learning algorithm to identify hypothetical individual end-of-treatment time points.

RESULTS

50 patients with DS-TB and 30 patients with MDR-TB were recruited in the German identification cohorts (DS-GIC and MDR-GIC, respectively); 28 patients with DS-TB and 32 patients with MDR-TB in the German validation cohorts (DS-GVC and MDR-GVC, respectively); and 52 patients with MDR-TB in the Romanian validation cohort (MDR-RVC). A 22-gene RNA model (TB22) that defined cure-associated end-of-therapy time points was derived from the DS- and MDR-GIC data. The TB22 model was superior to other published signatures to accurately predict clinical outcomes for patients in the DS-GVC (area under the curve 0.94, 95% CI 0.9-0.98) and suggests that cure may be achieved with shorter treatment durations for TB patients in the MDR-GIC (mean reduction 218.0 days, 34.2%; p<0.001), the MDR-GVC (mean reduction 211.0 days, 32.9%; p<0.001) and the MDR-RVC (mean reduction of 161.0 days, 23.4%; p=0.001).

CONCLUSION

Biomarker-guided management may substantially shorten the duration of therapy for many patients with MDR-TB.

Gene Set Enrichment Analysis Reveals Individual Variability in Host Responses in Tuberculosis Patients

Group-aggregated responses to tuberculosis (TB) have been well characterized on a molecular level. However, human beings differ and individual responses to infection vary. We have combined a novel approach to individual gene set analysis (GSA) with the clustering of transcriptomic profiles of TB patients from seven datasets in order to identify individual molecular endotypes of transcriptomic responses to TB. We found that TB patients differ with respect to the intensity of their hallmark interferon (IFN) responses, but they also show variability in their complement system, metabolic responses and multiple other pathways. This variability cannot be sufficiently explained with covariates such as gender or age, and the molecular endotypes are found across studies and populations. Using datasets from a Cynomolgus macaque model of TB, we revealed that transcriptional signatures of different molecular TB endotypes did not depend on TB progression post-infection. Moreover, we provide evidence that patients with molecular endotypes characterized by high levels of IFN responses (IFN-rich), suffered from more severe lung pathology than those with lower levels of IFN responses (IFN-low). Harnessing machine learning (ML) models, we derived gene signatures classifying IFN-rich and IFN-low TB endotypes and revealed that the IFN-low signature allowed slightly more reliable overall classification of TB patients from non-TB patients than the IFN-rich one. Using the paradigm of molecular endotypes and the ML-based predictions allows more precisely tailored treatment regimens, predicting treatment-outcome with higher accuracy and therefore bridging the gap between conventional treatment and precision medicine.

Multi-country evaluation of RISK6, a 6-gene blood transcriptomic signature, for tuberculosis diagnosis and treatment monitoring

There is a crucial need for non-sputum-based TB tests. Here, we evaluate the performance of RISK6, a human-blood transcriptomic signature, for TB screening, triage and treatment monitoring. RISK6 performance was also compared to that of two IGRAs: one based on RD1 antigens (QuantiFERON-TB Gold Plus, QFT-P, Qiagen) and one on recombinant M. tuberculosis HBHA expressed in Mycobacterium smegmatis (IGRA-rmsHBHA). In this multicenter prospective nested case-control study conducted in Bangladesh, Georgia, Lebanon and Madagascar, adult non-immunocompromised patients with bacteriologically confirmed active pulmonary TB (ATB), latent TB infection (LTBI) and healthy donors (HD) were enrolled. ATB patients were followed-up during and after treatment. Blood RISK6 scores were assessed using quantitative real-time PCR and evaluated by area under the receiver-operating characteristic curve (ROC AUC). RISK6 performance to discriminate ATB from HD reached an AUC of 0.94 (95% CI 0.89-0.99), with 90.9% sensitivity and 87.8% specificity, thus achieving the minimal WHO target product profile for a non-sputum-based TB screening test. Besides, RISK6 yielded an AUC of 0.93 (95% CI 0.85-1) with 90.9% sensitivity and 88.5% specificity for discriminating ATB from LTBI. Moreover, RISK6 showed higher performance (AUC 0.90, 95% CI 0.85-0.94) than IGRA-rmsHBHA (AUC 0.75, 95% CI 0.69-0.82) to differentiate TB infection stages. Finally, RISK6 signature scores significantly decreased after 2 months of TB treatment and continued to decrease gradually until the end of treatment reaching scores obtained in HD. We confirmed the performance of RISK6 signature as a triage TB test and its utility for treatment monitoring.

The Evaluation and Validation of Blood-Derived Novel Biomarkers for Precise and Rapid Diagnosis of Tuberculosis in Areas With High-TB Burden

Tuberculosis (TB) remains a highly contagious public health threat. Precise and prompt diagnosis and monitoring of treatment responses are urgently needed for clinics. To pursue novel and satisfied host blood-derived biomarkers, we streamlined a bioinformatic pipeline by integrating differentially expressed genes, a gene co-expression network, and short time-series analysis to mine the published transcriptomes derived from whole blood of TB patients in the GEO database, followed by validating the diagnostic performance of biomarkers in both independent datasets and blood samples of Chinese patients using quantitative real-time PCR (qRT-PCR). We found that four genes, namely UBE2L6 (Ubiquitin/ISG15-conjugating enzyme E2 L6), BATF2 (Basic leucine zipper transcriptional factor ATF-like), SERPING1 (Plasma protease C1 inhibitor), and VAMP5 (Vesicle-associated membrane protein 5), had high diagnostic value for active TB. The transcription levels of these four gene combinations can reach up to 88% sensitivity and 78% specificity (average) for the diagnosis of active TB; the highest sensitivity can achieve 100% by parallel of BATF2 and VAMP5, and the highest specificity can reach 89.5% through a combination of SERPIG1, UBE2L6, and VAMP5, which were significantly higher than 75.3% sensitivity and 69.1% specificity by T-SPOT.TB in the same patients. Quite unexpectedly, the gene set can assess the efficacy of anti-TB response and differentiate active TB from Latent TB infection. The data demonstrated these four biomarkers might have great potency and advantage over IGRAs in the diagnosis of TB.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Mycobacterium tuberculosis infection up-regulates MFN2 expression to promote NLRP3 inflammasome formation

Tuberculosis (TB), caused by the infection of Mycobacterium tuberculosis (MTB), is one of the leading causes of death worldwide, especially in children. However, the mechanisms by which MTB infects its cellular host, activates an immune response, and triggers inflammation remain unknown. Mitochondria play important roles in the initiation and activation of the nucleotide-binding oligomerization domain-like receptor with a pyrin domain 3 (NLRP3) inflammasome, where mitochondria-associated endoplasmic reticulum membranes (MAMs) may serve as the platform for inflammasome assembly and activation. Additionally, mitofusin 2 (MFN2) is implicated in the formation of MAMs, but, the roles of mitochondria and MFN2 in MTB infection have not been elucidated. Using mircroarry profiling of TB patients and in vitro MTB stimulation of macrophages, we observed an up-regulation of MFN2 in the peripheral blood mononuclear cells of active TB patients. Furthermore, we found that MTB stimulation by MTB-specific antigen ESAT-6 or lysate of MTB promoted MFN2 interaction with NLRP3 inflammasomes, resulting in the assembly and activation of the inflammasome and, subsequently, IL-1β secretion. These findings suggest that MFN2 and mitochondria play important role in the pathogen-host interaction during MTB infection.

Validation of Differentially Expressed Immune Biomarkers in Latent and Active Tuberculosis by Real-Time PCR

Tuberculosis (TB) remains a major global threat and diagnosis of active TB ((ATB) both extra-pulmonary (EPTB), pulmonary (PTB)) and latent TB (LTBI) infection remains challenging, particularly in high-burden countries which still rely heavily on conventional methods. Although molecular diagnostic methods are available, e.g., Cepheid GeneXpert, they are not universally available in all high TB burden countries. There is intense focus on immune biomarkers for use in TB diagnosis, which could provide alternative low-cost, rapid diagnostic solutions. In our previous gene expression studies, we identified peripheral blood leukocyte (PBL) mRNA biomarkers in a non-human primate TB aerosol-challenge model. Here, we describe a study to further validate select mRNA biomarkers from this prior study in new cohorts of patients and controls, as a prerequisite for further development. Whole blood mRNA was purified from ATB patients recruited in the UK and India, LTBI and two groups of controls from the UK (i) a low TB incidence region (CNTRLA) and (ii) individuals variably-domiciled in the UK and Asia ((CNTRLB), the latter TB high incidence regions). Seventy-two mRNA biomarker gene targets were analyzed by qPCR using the Roche Lightcycler 480 qPCR platform and data analyzed using GeneSpring™ 14.9 bioinformatics software. Differential expression of fifty-three biomarkers was confirmed between MTB infected, LTBI groups and controls, seventeen of which were significant using analysis of variance (ANOVA): CALCOCO2, CD52, GBP1, GBP2, GBP5, HLA-B, IFIT3, IFITM3, IRF1, LOC400759 (GBP1P1), NCF1C, PF4V1, SAMD9L, S100A11, TAF10, TAPBP, and TRIM25. These were analyzed using receiver operating characteristic (ROC) curve analysis. Single biomarkers and biomarker combinations were further assessed using simple arithmetic algorithms. Minimal combination biomarker panels were delineated for primary diagnosis of ATB (both PTB and EPTB), LTBI and identifying LTBI individuals at high risk of progression which showed good performance characteristics. These were assessed for suitability for progression against the standards for new TB diagnostic tests delineated in the published World Health Organization (WHO) technology product profiles (TPPs).

Peripheral Blood Mucosal-Associated Invariant T Cells in Tuberculosis Patients and Healthy Mycobacterium tuberculosis-Exposed Controls

BACKGROUND

In human blood, mucosal-associated invariant T (MAIT) cells are abundant T cells that recognize antigens presented on non-polymorphic major histocompatibility complex-related 1 (MR1) molecules. The MAIT cells are activated by mycobacteria, and prior human studies indicate that blood frequencies of MAIT cells, defined by cell surface markers, decline during tuberculosis (TB) disease, consistent with redistribution to the lungs.

METHODS

We tested whether frequencies of blood MAIT cells were altered in patients with TB disease relative to healthy Mycobacterium tuberculosis-exposed controls from Peru and South Africa. We quantified their frequencies using MR1 tetramers loaded with 5-(2-oxopropylideneamino)-6-D-ribitylaminouracil.

RESULTS

Unlike findings from prior studies, frequencies of blood MAIT cells were similar among patients with TB disease and latent and uninfected controls. In both cohorts, frequencies of MAIT cells defined by MR1-tetramer staining and coexpression of CD161 and the T-cell receptor alpha variable gene TRAV1-2 were strongly correlated. Disease severity captured by body mass index or TB disease transcriptional signatures did not correlate with MAIT cell frequencies in patients with TB.

CONCLUSIONS

Major histocompatibility complex (MHC)-related 1-restrictied MAIT cells are detected at similar levels with tetramers or surface markers. Unlike MHC-restricted T cells, blood frequencies of MAIT cells are poor correlates of TB disease but may play a role in pathophysiology.

Identification of Reduced Host Transcriptomic Signatures for Tuberculosis Disease and Digital PCR-Based Validation and Quantification

Recently, host whole blood gene expression signatures have been identified for diagnosis of tuberculosis (TB). Absolute quantification of the concentrations of signature transcripts in blood have not been reported, but would facilitate diagnostic test development. To identify minimal transcript signatures, we applied a transcript selection procedure to microarray data from African adults comprising 536 patients with TB, other diseases (OD) and latent TB (LTBI), divided into training and test sets. Signatures were further investigated using reverse transcriptase (RT)-digital PCR (dPCR). A four-transcript signature (GBP6, TMCC1, PRDM1, and ARG1) measured using RT-dPCR distinguished TB patients from those with OD (area under the curve (AUC) 93.8% (CI_95% 82.2-100%). A three-transcript signature (FCGR1A, ZNF296, and C1QB) differentiated TB from LTBI (AUC 97.3%, CI_95%: 93.3-100%), regardless of HIV. These signatures have been validated across platforms and across samples offering strong, quantitative support for their use as diagnostic biomarkers for TB.

Human tuberculosis and Mycobacterium tuberculosis complex: A review on genetic diversity, pathogenesis and omics approaches in host biomarkers discovery

Mycobacterium tuberculosis complex (MTBC) refers to a group of mycobacteria encompassing nine members of closely related species that causes tuberculosis in animals and humans. Among the nine members, Mycobacterium tuberculosis (M. tuberculosis) remains the main causative agent for human tuberculosis that results in high mortality and morbidity globally. In general, MTBC species are low in diversity but exhibit distinctive biological differences and phenotypes among different MTBC lineages. MTBC species are likely to have evolved from a common ancestor through insertions/deletions processes resulting in species speciation with different degrees of pathogenicity. The pathogenesis of human tuberculosis is complex and remains poorly understood. It involves multi-interactions or evolutionary co-options between host factors and bacterial determinants for survival of the MTBC. Granuloma formation as a protection or survival mechanism in hosts by MTBC remains controversial. Additionally, MTBC species are capable of modulating host immune response and have adopted several mechanisms to evade from host immune attack in order to survive in humans. On the other hand, current diagnostic tools for human tuberculosis are inadequate and have several shortcomings. Numerous studies have suggested the potential of host biomarkers in early diagnosis of tuberculosis, in disease differentiation and in treatment monitoring. "Multi-omics" approaches provide holistic views to dissect the association of MTBC species with humans and offer great advantages in host biomarkers discovery. Thus, in this review, we seek to understand how the genetic variations in MTBC lead to species speciation with different pathogenicity. Furthermore, we also discuss how the host and bacterial players contribute to the pathogenesis of human tuberculosis. Lastly, we provide an overview of the journey of "omics" approaches in host biomarkers discovery in human tuberculosis and provide some interesting insights on the challenges and directions of "omics" approaches in host biomarkers innovation and clinical implementation.

Transcriptomic signatures for diagnosing tuberculosis in clinical practice: a prospective, multicentre cohort study

Background

Blood transcriptomic signatures for diagnosis of tuberculosis have shown promise in case-control studies, but none have been prospectively designed or validated in adults presenting with the full clinical spectrum of suspected tuberculosis, including extrapulmonary tuberculosis and common differential diagnoses that clinically resemble tuberculosis. We aimed to evaluate the diagnostic accuracy of transcriptomic signatures in patients presenting with clinically suspected tuberculosis in routine practice.

Methods

The Validation of New Technologies for Diagnostic Evaluation of Tuberculosis (VANTDET) study was nested within a prospective, multicentre cohort study in secondary care in England (IDEA 11/H0722/8). Patients (aged ≥16 years) suspected of having tuberculosis in the routine clinical inpatient and outpatient setting were recruited at ten National Health Service hospitals in England for IDEA and were included in VANTDET if they provided consent for genomic analysis. Patients had whole blood taken for microarray analysis to measure abundance of transcripts and were followed up for 6–12 months to determine final diagnoses on the basis of predefined diagnostic criteria. The diagnostic accuracy of six signatures derived from the cohort and three previously published transcriptomic signatures with potentially high diagnostic performance were assessed by calculating area under the receiver-operating characteristic curves (AUC-ROCs), sensitivities, and specificities.

Findings

Between Nov 25, 2011, and Dec 31, 2013, 1162 participants were enrolled. 628 participants (aged ≥16 years) were included in the analysis, of whom 212 (34%) had culture-confirmed tuberculosis, 89 (14%) had highly probable tuberculosis, and 327 (52%) had tuberculosis excluded. The novel signature with highest performance for identifying all active tuberculosis gave an AUC-ROC of 0·87 (95% CI 0·81–0·92), sensitivity of 77% (66–87), and specificity of 84% (74–91). The best-performing published signature gave an AUC-ROC of 0·83 (0·80–0·86), sensitivity of 78% (73–83), and specificity of 76% (70–80). For detecting highly probable tuberculosis, the best novel signature yielded results of 0·86 (0·71–0·95), 77% (56–94%), and 77% (57–95%). None of the relevant cohort-derived or previously published signatures achieved the WHO-defined targets of paired sensitivity and specificity for a non-sputum-based diagnostic test.

Interpretation

In a clinically representative cohort in routine practice in a low-incidence setting, transcriptomic signatures did not have adequate accuracy for diagnosis of tuberculosis, including in patients with highly probable tuberculosis where the unmet need is greatest. These findings suggest that transcriptomic signatures have little clinical utility for diagnostic assessment of suspected tuberculosis.

Funding

National Institute for Health Research.

Pandemics Throughout History

The emergence and spread of infectious diseases with pandemic potential occurred regularly throughout history. Major pandemics and epidemics such as plague, cholera, flu, severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) have already afflicted humanity. The world is now facing the new coronavirus disease 2019 (COVID-19) pandemic. Many infectious diseases leading to pandemics are caused by zoonotic pathogens that were transmitted to humans due to increased contacts with animals through breeding, hunting and global trade activities. The understanding of the mechanisms of transmission of pathogens to humans allowed the establishment of methods to prevent and control infections. During centuries, implementation of public health measures such as isolation, quarantine and border control helped to contain the spread of infectious diseases and maintain the structure of the society. In the absence of pharmaceutical interventions, these containment methods have still been used nowadays to control COVID-19 pandemic. Global surveillance programs of water-borne pathogens, vector-borne diseases and zoonotic spillovers at the animal-human interface are of prime importance to rapidly detect the emergence of infectious threats. Novel technologies for rapid diagnostic testing, contact tracing, drug repurposing, biomarkers of disease severity as well as new platforms for the development and production of vaccines are needed for an effective response in case of pandemics.

Tuberculosis comorbidity with rheumatoid arthritis: Gene signatures, associated biomarkers, and screening

Rheumatoid arthritis (RA) is known to be related to an elevated risk of infections because of its pathobiology and the use of immunosuppressive therapies. Reactivation of latent tuberculosis (TB) infection is a serious issue in patients with RA, especially after receiving anti-TNFs therapy. TNF blocking reinforces the TB granuloma formation and maintenance and the growth of Mycobacterium tuberculosis (Mtb). After intercurrent of TB infection, the standard recommendation is that the treatment with TNF inhibitors to be withheld despite its impressive effect on suppression of inflammation until the infection has resolved. Knowing pathways and mechanisms that are common between two diseases might help to find the mechanistic basis of this comorbidity, as well as provide us a new approach to apply them as therapeutic targets or diagnostic biomarkers. Also, screening for latent TB before initiation of an anti-TNF therapy can minimize complications. This review summarizes the shared gene signature between TB and RA and discusses the biomarkers for early detection of this infection, and screening procedures as well.

Tuberculosis-Associated MicroRNAs: From Pathogenesis to Disease Biomarkers

Tuberculosis (TB) caused by Mycobacterium tuberculosis is one of the most lethal infectious diseases with estimates of approximately 1.4 million human deaths in 2018. M. tuberculosis has a well-established ability to circumvent the host immune system to ensure its intracellular survival and persistence in the host. Mechanisms include subversion of expression of key microRNAs (miRNAs) involved in the regulation of host innate and adaptive immune response against M. tuberculosis. Several studies have reported differential expression of miRNAs during active TB and latent tuberculosis infection (LTBI), suggesting their potential use as biomarkers of disease progression and response to anti-TB therapy. This review focused on the miRNAs involved in TB pathogenesis and on the mechanism through which miRNAs induced during TB modulate cell antimicrobial responses. An attentive study of the recent literature identifies a group of miRNAs, which are differentially expressed in active TB vs. LTBI or vs. treated TB and can be proposed as candidate biomarkers.

Factors Associated with Mortality in Critically Ill Patients Diagnosed with Hospital Acquired Infections

Objective

Evaluate host and pathogen factors associated with mortality in those with hospital acquired infections (HAI) in a tertiary intensive care unit in Brazil.

Methods

Observational and analytical cohort single center study in a general intensive care unit (ICU) in Northeastern Brazil between January 2016 and August 2018, including those over 18 years of age admitted to the ICU found to have a HAI.

Results

A total of 165 patients were included, with a mean age of 72 years and male predominance (53.3%) and observed mortality of 46%. Mortality in those with HAI was significantly associated with older age, increased ICU length of stay and readmission to the ICU in univariate analysis. Multivariate analysis revealed that development of septic shock and obtundation during ICU admission was significantly associated with an increased risk of death (OR: 6.94, 95% CI 1.23–39.27, OR: 2.48, 95% CI 1.17–5.29, respectively). A trend towards mortality risk was noted in those with increased age and prior cardiovascular disease. Surprisingly, mortality risk was independent of site of infection, type of pathogen and antibiotic resistance. Furthermore, having more than one HAI over the course of the ICU admission did not impact mortality.

Conclusion

Risk of death in those with HAI is associated with obtundation and septic shock, in addition to vasopressor use. Host factors, rather than pathogen-specific characteristics or infecting site, impact risk of death related to HAI in the ICU.

RISK6, a 6-gene transcriptomic signature of TB disease risk, diagnosis and treatment response

Improved tuberculosis diagnostics and tools for monitoring treatment response are urgently needed. We developed a robust and simple, PCR-based host-blood transcriptomic signature, RISK6, for multiple applications: identifying individuals at risk of incident disease, as a screening test for subclinical or clinical tuberculosis, and for monitoring tuberculosis treatment. RISK6 utility was validated by blind prediction using quantitative real-time (qRT) PCR in seven independent cohorts. Prognostic performance significantly exceeded that of previous signatures discovered in the same cohort. Performance for diagnosing subclinical and clinical disease in HIV-uninfected and HIV-infected persons, assessed by area under the receiver-operating characteristic curve, exceeded 85%. As a screening test for tuberculosis, the sensitivity at 90% specificity met or approached the benchmarks set out in World Health Organization target product profiles for non-sputum-based tests. RISK6 scores correlated with lung immunopathology activity, measured by positron emission tomography, and tracked treatment response, demonstrating utility as treatment response biomarker, while predicting treatment failure prior to treatment initiation. Performance of the test in capillary blood samples collected by finger-prick was noninferior to venous blood collected in PAXgene tubes. These results support incorporation of RISK6 into rapid, capillary blood-based point-of-care PCR devices for prospective assessment in field studies.

A transcriptional blood signature distinguishes early tuberculosis disease from latent tuberculosis infection and uninfected individuals in a Vietnamese cohort

OBJECTIVES

Global tuberculosis (TB) control is restricted by the failure to detect an estimated 3.3 million TB cases annually. In the majority of TB endemic settings, sputum smear microscopy is used to diagnose TB, but this test is insensitive for TB in its early stages. The objective of this study is to establish a concise gene signature that discriminates between individuals with early TB disease, latent TB infection (LTBI) and those without infection.

METHODS

This is a case control study nested within a cluster-randomised trial of population screening for active TB using Xpert MTB/RIF. Whole blood samples from 303 participants with active TB (97), LTBI (92) and uninfected individuals (114) were subject to transcriptomic analysis of selected target genes based on a systematic review of previous studies.

RESULTS

Analysis of 82 genes identified a pattern of differentially expressed genes in TB disease. A seven gene signature was identified that distinguished between TB disease and no TB disease with an AUC of 0.86 (95% CI: 0.80-0.91), and between TB disease from LTBI with an AUC of 0.88 (95% CI: 0.82-0.93).

CONCLUSION

This gene signature accurately distinguishes early TB disease from those without TB disease or infection, in the context of community-wide TB screening. It could be used as a non-sputum based screening tool or triage test to detect prevalent cases of TB in the community.

Biomarker Signatures of Quality for Engineering Nasal Chondrocyte-Derived Cartilage

The definition of quality controls for cell therapy and engineered product manufacturing processes is critical for safe, effective, and standardized clinical implementation. Using the example context of cartilage grafts engineered from autologous nasal chondrocytes, currently used for articular cartilage repair in a phase II clinical trial, we outlined how gene expression patterns and generalized linear models can be introduced to define molecular signatures of identity, purity, and potency. We first verified that cells from the biopsied nasal cartilage can be contaminated by cells from a neighboring tissue, namely perichondrial cells, and discovered that they cannot deposit cartilaginous matrix. Differential analysis of gene expression enabled the definition of identity markers for the two cell populations, which were predictive of purity in mixed cultures. Specific patterns of expression of the same genes were significantly correlated with cell potency, defined as the capacity to generate tissues with histological and biochemical features of hyaline cartilage. The outlined approach can now be considered for implementation in a good manufacturing practice setting, and offers a paradigm for other regenerative cellular therapies.

Blood transcriptional biomarkers for active pulmonary tuberculosis in a high-burden setting: a prospective, observational, diagnostic accuracy study

BACKGROUND

Blood transcriptional signatures are candidates for non-sputum triage or confirmatory tests of tuberculosis. Prospective head-to-head comparisons of their diagnostic accuracy in real-world settings are necessary to assess their clinical use. We aimed to compare the diagnostic accuracy of candidate transcriptional signatures identified by systematic review, in a setting with a high burden of tuberculosis and HIV.

METHODS

We did a prospective observational study nested within a diagnostic accuracy study of sputum Xpert MTB/RIF (Xpert) and Xpert MTB/RIF Ultra (Ultra) tests for pulmonary tuberculosis. We recruited consecutive symptomatic adults aged 18 years or older self-presenting to a tuberculosis clinic in Cape Town, South Africa. Participants provided blood for RNA sequencing, and sputum samples for liquid culture and molecular testing using Xpert and Ultra. We assessed the diagnostic accuracy of candidate blood transcriptional signatures for active tuberculosis (including those intended to distinguish active tuberculosis from other diseases) identified by systematic review, compared with culture or Xpert MTB/RIF positivity as the standard reference. In our primary analysis, patients with tuberculosis were defined as those with either a positive liquid culture or Xpert result. Patients with missing blood RNA or sputum results were excluded. Our primary objective was to benchmark the diagnostic accuracy of candidate transcriptional signatures against the WHO target product profile (TPP) for a tuberculosis triage test.

FINDINGS

Between Feb 12, 2016, and July 18, 2017, we obtained paired sputum and RNA sequencing data from 181 participants, 54 (30%) of whom had confirmed pulmonary tuberculosis. Of 27 eligible signatures identified by systematic review, four achieved the highest diagnostic accuracy with similar area under the receiver operating characteristic curves (Sweeney3: 90·6% [95% CI 85·6-95·6]; Kaforou25: 86·9% [80·9-92·9]; Roe3: 86·9% [80·3-93·5]; and BATF2: 86·8% [80·6-93·1]), independent of age, sex, HIV status, previous tuberculosis, or sputum smear result. At test thresholds that gave 70% specificity (the minimum WHO TPP specificity for a triage test), these four signatures achieved sensitivities between 83·3% (95% CI 71·3-91·0) and 90·7% (80·1-96·0). No signature met the optimum criteria, of 95% sensitivity and 80% specificity proposed by WHO for a triage test, or the minimum criteria (of 65% sensitivity and 98% specificity) for a confirmatory test, but all four correctly identified Ultra-positive, culture-negative patients.

INTERPRETATION

Selected blood transcriptional signatures met the minimum WHO benchmarks for a tuberculosis triage test but not for a confirmatory test. Further development of the signatures is warranted to investigate their possible effects on clinical and health economic outcomes as part of a triage strategy, or when used as add-on confirmatory test in conjunction with the highly sensitive Ultra test for Mycobacterium tuberculosis DNA.

FUNDING

Royal Society Newton Advanced Fellowship, Wellcome Trust, National Institute of Health Research, and UK Medical Research Council.

Concise whole blood transcriptional signatures for incipient tuberculosis: a systematic review and patient-level pooled meta-analysis

BACKGROUND

Multiple blood transcriptional signatures have been proposed for identification of active and incipient tuberculosis. We aimed to compare the performance of systematically identified candidate signatures for incipient tuberculosis and to benchmark these against WHO targets.

METHODS

We did a systematic review and individual participant data meta-analysis. We searched Medline and Embase for candidate whole blood mRNA signatures discovered with the primary objective of diagnosis of active or incipient tuberculosis, compared with controls who were healthy or had latent tuberculosis infection. We tested the performance of eligible signatures in whole blood transcriptomic datasets, in which sampling before tuberculosis diagnosis was done and time to disease was available. Culture-confirmed and clinically or radiologically diagnosed pulmonary or extrapulmonary tuberculosis cases were included. Non-progressor (individuals who remained tuberculosis-free during follow-up) samples with less than 6 months of follow-up from the date of sample collection were excluded, as were participants with prevalent tuberculosis and those who received preventive therapy. Scores were calculated for candidate signatures for each participant in the pooled dataset. Receiver operating characteristic curves, sensitivities, and specificities were examined using prespecified intervals to tuberculosis (<3 months, <6 months, <1 year, and <2 years) from sample collection. This study is registered with PROSPERO, number CRD42019135618.

RESULTS

We tested 17 candidate mRNA signatures in a pooled dataset from four eligible studies comprising 1126 samples. This dataset included 183 samples from 127 incipient tuberculosis cases in South Africa, Ethiopia, The Gambia, and the UK. Eight signatures (comprising 1-25 transcripts) that predominantly reflect interferon and tumour necrosis factor-inducible gene expression, had equivalent diagnostic accuracy for incipient tuberculosis over a 2-year period with areas under the receiver operating characteristic curves ranging from 0·70 (95% CI 0·64-0·76) to 0·77 (0·71-0·82). The sensitivity of all eight signatures declined with increasing disease-free time interval. Using a threshold derived from two SDs above the mean of uninfected controls to prioritise specificity and positive-predictive value, the eight signatures achieved sensitivities of 24·7-39·9% over 24 months and of 47·1-81·0% over 3 months, with corresponding specificities of more than 90%. Based on pre-test probability of 2%, the eight signatures achieved positive-predictive values ranging from 6·8-9·4% over 24 months and 11·2-14·4% over 3 months. When using biomarker thresholds maximising sensitivity and specificity with equal weighting to both, no signature met the minimum WHO target product profile parameters for incipient tuberculosis biomarkers over a 2-year period.

INTERPRETATION

Blood transcriptional biomarkers reflect short-term risk of tuberculosis and only exceed WHO benchmarks if applied to 3-6-month intervals. Serial testing among carefully selected target groups might be required for optimal implementation of these biomarkers.

FUNDING

Wellcome Trust and National Institute for Health Research.

Expression and production of the SERPING1-encoded endogenous complement regulator C1-inhibitor in multiple cohorts of tuberculosis patients

BACKGROUND

To facilitate better discrimination between patients with active tuberculosis (TB) and latent TB infection (LTBI), whole blood transcriptomic studies have been performed to identify novel candidate host biomarkers. SERPING1, which encodes C1-inhibitor (C1-INH), the natural inhibitor of the C1-complex has emerged as candidate biomarker. Here we collated and analysed SERPING1 expression data and subsequently determined C1-INH protein levels in four cohorts of patients with TB.

METHODS

SERPING1 expression data were extracted from online deposited datasets. C1-INH protein levels were determined by ELISA in sera from individuals with active TB, LTBI as well as other disease controls in geographically diverse cohorts.

FINDINGS

SERPING1 expression was increased in patients with active TB compared to healthy controls (8/11 cohorts), LTBI (13/14 cohorts) and patients with other (non-TB) lung-diseases (7/7 cohorts). Serum levels of C1-INH were significantly increased in The Gambia and Italy in patients with active TB relative to the endemic controls but not in South Africa or Korea. In the largest cohort (n = 50), with samples collected longitudinally, normalization of C1-INH levels following successful TB treatment was observed. This cohort, also showed the most abundant increase in C1-INH, and a positive correlation between C1q and C1-INH levels. Combined presence of increased levels of both C1q and C1-INH had high specificity for active TB (96 %) but only very modest sensitivity 38 % compared to the endemic controls.

INTERPRETATION

SERPING1 transcript expression is increased in TB patients, while serum protein levels of C1-INH were increased in half of the cohorts analysed.

Single-cell transcriptomics of blood reveals a natural killer cell subset depletion in tuberculosis

Background

Tuberculosis (TB) continues to be a critical global health problem, which killed millions of lives each year. Certain circulating cell subsets are thought to differentially modulate the host immune response towards Mycobacterium tuberculosis (Mtb) infection, but the nature and function of these subsets is unclear.

Methods

Peripheral blood mononuclear cells (PBMC) were isolated from healthy controls (HC), latent tuberculosis infection (LTBI) and active tuberculosis (TB) and then subjected to single-cell RNA sequencing (scRNA-seq) using 10 × Genomics platform. Unsupervised clustering of the cells based on the gene expression profiles using the Seurat package and passed to tSNE for clustering visualization. Flow cytometry was used to validate the subsets identified by scRNA-Seq.

Findings

Cluster analysis based on differential gene expression revealed both known and novel markers for all main PBMC cell types and delineated 29 cell subsets. By comparing the scRNA-seq datasets from HC, LTBI and TB, we found that infection changes the frequency of immune-cell subsets in TB. Specifically, we observed gradual depletion of a natural killer (NK) cell subset (CD3-CD7+GZMB+) from HC, to LTBI and TB. We further verified that the depletion of CD3-CD7+GZMB+ subset in TB and found an increase in this subset frequency after anti-TB treatment. Finally, we confirmed that changes in this subset frequency can distinguish patients with TB from LTBI and HC.

Interpretation

We propose that the frequency of CD3-CD7+GZMB+ in peripheral blood could be used as a novel biomarker for distinguishing TB from LTBI and HC.

Fund

The study was supported by Natural Science Foundation of China (81770013, 81525016, 81772145, 81871255 and 91942315), National Science and Technology Major Project (2017ZX10201301), Science and Technology Project of Shenzhen (JCYJ20170412101048337) and Guangdong Provincial Key Laboratory of Regional Immunity and Diseases (2019B030301009). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Use of DosR Dormancy Antigens from Mycobacterium tuberculosis for Serodiagnosis of Active and Latent Tuberculosis

As more than two billion people possibly have a latent tuberculosis (LTB) infection, early LTB diagnosis is crucial for the efficient control and elimination of tuberculosis (TB). The aim of this study is to detect the serum antibody responses to dormancy-related DosR regulon antigens of Mycobacterium tuberculosis for the diagnosis of active and latent TB infections. A membrane array with 25 latency antigens detected by silver-enhanced gold nanoparticles was used to determine the corresponding cognate antibody levels in clinical serum samples from healthy controls, TB patients, and individuals with LTB. The array is sized to fit into a 24-well ELISA plate and follows an ELISA-like experimental procedure without expensive instrumentation. Linear discriminant analysis (LDA) of the resulting antibody profiling data set identified a panel of nine DosR antigens with significant discriminatory capability among different subjects with ≥90% sensitivity, specificity, and overall accuracy. Furthermore, the high predictive performance validated by an independent test sample set reflects the robustness and reliability of the LDA classification model. Our current data demonstrate that the nine DosR antigen combination associated with the proposed membrane array platform is a clinically feasible approach for distinguishing different TB infection statuses. The proposed methodology in this study could be further developed for multiple disease serodiagnoses with high sensitivity and specificity.