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

Quantitative and Longitudinal Assessment of Systemic Innate Immunity in Health and Disease Using a 2D Gene Model

Dysregulation of innate immunity is deeply involved in infectious and autoimmune diseases. For a better understanding of pathogenesis and improved management of these diseases, it is of vital importance to implement convenient monitoring of systemic innate immunity. Built upon our previous works on the host transcriptional response to infection in peripheral blood, we proposed a 2D gene model for the simultaneous assessment of two major components of systemic innate immunity, including VirSig as the signature of the host response to viral infection and BacSig as the signature of the host response to bacterial infection. The revelation of dysregulation in innate immunity by this 2D gene model was demonstrated with a wide variety of transcriptome datasets. In acute infection, distinctive patterns of VirSig and BacSig activation were observed in viral and bacterial infection. In comparison, both signatures were restricted to a defined range in the vast majority of healthy adults, regardless of age. In addition, BacSig showed significant elevation during pregnancy and an upward trend during development. In tuberculosis (TB), elevation of BacSig and VirSig was observed in a significant portion of active TB patients, and abnormal BacSig was also associated with a longer treatment course. In cystic fibrosis (CF), abnormal BacSig was observed in a subset of patients, and no overall change in BacSig abnormality was observed after the drug treatment. In systemic sclerosis-associated interstitial lung disease (SSc-ILD), significant elevation of VirSig and BacSig was observed in some patients, and treatment with a drug led to the further deviation of BacSig from the control level. In systemic lupus erythematosus (SLE), positivity for the anti-Ro autoantibody was associated with significant elevation of VirSig in SLE patients, and the additive effect of VirSig/BacSig activation was also observed in SLE patients during pregnancy. Overall, these data demonstrated that the 2D gene model can be used to assess systemic innate immunity in health and disease, with the potential clinical applications including patient stratification, prescription of antibiotics, understanding of pathogenesis, and longitudinal monitoring of treatment response.

Interferon-gamma release assay conversion after Mycobacterium tuberculosis exposure specifically associates with greater risk of progression to tuberculosis: A prospective cohort study in Leicester, UK

OBJECTIVES

We investigated whether quantifying the serial QuantiFERON-TB Gold (QFT) response improves tuberculosis (TB) risk stratification in pulmonary TB (PTB) contacts.

METHODS

A total of 297 untreated adult household PTB contacts, QFT tested at baseline and 3 months after index notification, were prospectively observed (median 1460 days). Normal variance of serial QFT responses was established in 46 extrapulmonary TB contacts. This informed categorisation of the response in QFT-positive PTB contacts as converters, persistently QFT-positive with significant increase (PPincrease), and without significant increase (PPno-increase).

RESULTS

In total, eight co-prevalent TB (disease ≤3 months after index notification) and 12 incident TB (>3 months after index notification) cases were diagnosed. Genetic linkage to the index strain was confirmed in all culture-positive progressors. The cumulative 2-year incident TB risk in QFT-positive contacts was 8.4% (95% confidence interval, 3.0-13.6%); stratifying by serial QFT response, significantly higher risk was observed in QFT converters (28%), compared with PPno-increase (4.8%) and PPincrease (3.7%). Converters were characterised by exposure to index cases with a shorter interval from symptom onset to diagnosis (median reduction 50.0 days, P = 0.013).

CONCLUSIONS

QFT conversion, rather than quantitative changes of a persistently positive serial QFT response, is associated with greater TB risk and exposure to rapidly progressive TB.

Diagnostic markers reflecting dysregulation of the host response in the transition to tuberculosis disease

Sustained control of Mycobacterium tuberculosis infection without evidence of disease is based on a finely tuned balance between pro- and anti-inflammatory responses. Loss of this balance leads to tuberculosis (TB) disease, in which exacerbated myeloid and neutrophil activation is common. Proteomic and transcriptomic assessment of the host response can detect increasing immune activation associated with TB disease progression several months before clinical disease. Future diagnostic methods based on measuring host response biomarkers that are able to detect this dysregulation could therefore be valuable in the early detection of TB disease progression.

Neutrophil-plasmacytoid dendritic cell interaction leads to production of type I IFN in response to Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) can cause a latent infection that sometimes progresses to clinically active tuberculosis (TB). Type I interferons (IFN-I) have been implicated in initiating the progression from latency to active TB, in part because IFN-I stimulated genes are the earliest genes to be upregulated in patients as they advance to active TB. Plasmacytoid dendritic cells (pDCs) are major producers of IFN-I during viral infections and in response to autoimmune-induced neutrophil extracellular traps. pDCs have also been suggested to be the major producers of IFN-I during Mtb infection of mice and nonhuman primates, but direct evidence has been lacking. Here, we found that Mtb did not stimulate isolated human pDCs to produce IFN-I, but human neutrophils infected with Mtb-activated co-cultured pDCs to do so. Mtb-infected neutrophils produced neutrophil extracellular traps, whose exposed DNA is a well-known mechanism to activate pDCs to secrete IFN-I. We conclude that pDCs contribute to the IFN-I response during Mtb infection by interacting with infected neutrophils which may then promote Mtb pathogenesis.

Whole blood RNA signatures in tuberculosis patients receiving H56:IC31 vaccine as adjunctive therapy

Introduction

Therapeutic vaccination in tuberculosis (TB) represents a Host Directed Therapy strategy which enhances immune responses in order to improve clinical outcomes and shorten TB treatment. Previously, we have shown that the subunit H56:IC31 vaccine induced both humoral and cellular immune responses when administered to TB patients adjunctive to standard TB treatment (TBCOX2 study, NCT02503839). Here we present the longitudinal whole blood gene expression patterns in H56:IC31 vaccinated TB patients compared to controls receiving standard TB treatment only.

Methods

The H56:IC31 group (N=11) and Control group (N=7) underwent first-line TB treatment for 182 days. The H56:IC31 group received 5 micrograms of the H56:IC31 vaccine (Statens Serum Institut; SSI, Valneva Austria GmbH) intramuscularly at day 84 and day 140. Total RNA was extracted from whole blood samples collected in PAXgene tubes on days 0, 84, 98, 140, 154, 182 and 238. The expression level of 183 immune-related genes was measured by high-throughput microfluidic qPCR (Biomark HD system, Standard BioTools).

Results

The targeted gene expression profiling unveiled the upregulation of modules such as interferon (IFN) signalling genes, pattern recognition receptors and small nucleotide guanosine triphosphate (GTP)-ases in the vaccinated group compared to controls two weeks after administration of the first H56:IC31 vaccine. Additionally, the longitudinal analysis of the Adolescent Cohort Study-Correlation of Risk (ACS-COR) signature showed a progressive downregulation in both study arms towards the end of TB treatment, in congruence with reported treatment responses and clinical improvements. Still, two months after the end of TB treatment, vaccinated patients, and especially those developing both cellular and humoral vaccine responses, showed a lower expression of the ACS-COR genes compared to controls.

Discussion

Our data report gene expression patterns following H56:IC31 vaccination which might be interpreted as a lower risk of relapse in therapeutically vaccinated patients. Further studies are needed to conclude if these gene expression patterns could be used as prognostic biosignatures for therapeutic TB vaccine responses.

PET-CT-guided characterisation of progressive, preclinical tuberculosis infection and its association with low-level circulating Mycobacterium tuberculosis DNA in household contacts in Leicester, UK: a prospective cohort study

BACKGROUND

Incipient tuberculosis, a progressive state of Mycobacterium tuberculosis infection with an increased risk of developing into tuberculosis disease, remains poorly characterised. Animal models suggest an association of progressive infection with bacteraemia. Circulating M tuberculosis DNA has previously been detected in pulmonary tuberculosis by use of Actiphage, a bacteriophage-based real-time PCR assay. We aimed to investigate whether serial [18F]fluorodeoxyglucose ([18F]FDG)-PET-CT could be used to characterise the state and progressive trajectory of incipient tuberculosis, and examine whether these PET-CT findings are associated with Actiphage-based detection of circulating M tuberculosis DNA.

METHODS

We did a prospective 12-month cohort study in healthy, asymptomatic adults (aged ≥16 years) who were household contacts of patients with pulmonary tuberculosis, and who had a clinical phenotype of latent tuberculosis infection, in Leicester, UK. Actiphage testing of participants' blood samples was done at baseline, and [18F]FDG PET-CT at baseline and after 3 months. Baseline PET-CT features were classified as positive, indeterminate, or negative, on the basis of the quantitation (maximum standardised uptake value [SUVmax]) and distribution of [18F]FDG uptake. Microbiological sampling was done at amenable sites of [18F]FDG uptake. Changes in [18F]FDG uptake after 3 months were quantitatively categorised as progressive, stable, or resolving. Participants received treatment if features of incipient tuberculosis, defined as microbiological detection of M tuberculosis or progressive PET-CT change, were identified.

FINDINGS

20 contacts were recruited between Aug 5 and Nov 5, 2020; 16 of these participants had a positive result on IFNγ release assay (QuantiFERON-TB Gold Plus [QFT]) indicating tuberculosis infection. Baseline PET-CT scans were positive in ten contacts (all QFT positive), indeterminate in six contacts (three QFT positive), and negative in four contacts (three QFT positive). Four of eight PET-CT-positive contacts sampled had M tuberculosis identified (three through culture, one through Xpert MTB/RIF Ultra test) from intrathoracic lymph nodes or bronchial wash and received full antituberculosis treatment. Two further unsampled PET-CT-positive contacts were also treated: one with [18F]FDG uptake in the lung (SUVmax 9·4) received empirical antituberculosis treatment and one who showed progressive [18F]FDG uptake received preventive treatment. The ten untreated contacts with [18F]FDG uptake at baseline (seven QFT positive) had stable or resolving changes at follow-up and remained free of tuberculosis disease after 12 months. A positive baseline Actiphage test was associated with the presence of features of incipient tuberculosis requiring treatment (p=0·018).

INTERPRETATION

Microbiological and inflammatory features of incipient tuberculosis can be visualised on PET-CT and are associated with M tuberculosis detection in the blood, supporting the development of pathogen-directed blood biomarkers of tuberculosis risk.

FUNDING

MRC Confidence in Concept.

Mechanopathology of biofilm-like Mycobacterium tuberculosis cords

Mycobacterium tuberculosis (Mtb) cultured axenically without detergent forms biofilm-like cords, a clinical identifier of virulence. In lung-on-chip (LoC) and mouse models, cords in alveolar cells contribute to suppression of innate immune signaling via nuclear compression. Thereafter, extracellular cords cause contact-dependent phagocyte death but grow intercellularly between epithelial cells. The absence of these mechanopathological mechanisms explains the greater proportion of alveolar lesions with increased immune infiltration and dissemination defects in cording-deficient Mtb infections. Compression of Mtb lipid monolayers induces a phase transition that enables mechanical energy storage. Agent-based simulations demonstrate that the increased energy storage capacity is sufficient for the formation of cords that maintain structural integrity despite mechanical perturbation. Bacteria in cords remain translationally active despite antibiotic exposure and regrow rapidly upon cessation of treatment. This study provides a conceptual framework for the biophysics and function in tuberculosis infection and therapy of cord architectures independent of mechanisms ascribed to single bacteria.

Granulocytes subsets and their divergent functions in host resistance to Mycobacterium tuberculosis - a 'tipping-point' model of disease exacerbation

Granulocytes are innate immune effector cells with essential functions in host resistance to bacterial infections. I will discuss emerging evidence that during Mycobacterium tuberculosis infection, counter-intuitively, eosinophils are host-protective while neutrophils are host detrimental. Additionally, I will propose a 'tipping-point' model in which neutrophils are an integral part of a feedforward loop driving tuberculosis disease exacerbation.

Exposure of Mycobacterium tuberculosis to human alveolar lining fluid shows temporal and strain-specific adaptation to the lung environment

Upon infection, Mycobacterium tuberculosis ( M.tb ) reaches the alveolar space and comes in close contact with human alveolar lining fluid (ALF) for an uncertain period of time prior to its encounter with alveolar cells. We showed that homeostatic ALF hydrolytic enzymes modify the M.tb cell envelope, driving M.tb -host cell interactions. Still, the contribution of ALF during M.tb infection is poorly understood. Here, we exposed 4 M.tb strains with different levels of virulence, transmissibility, and drug resistance (DR) to physiological concentrations of human ALF for 15-min and 12-h, and performed RNA sequencing. Gene expression analysis showed a temporal and strain-specific adaptation to human ALF. Differential expression (DE) of ALF-exposed vs. unexposed M.tb revealed a total of 397 DE genes associated with lipid metabolism, cell envelope and processes, intermediary metabolism and respiration, and regulatory proteins, among others. Most DE genes were detected at 12-h post-ALF exposure, with DR- M.tb strain W-7642 having the highest number of DE genes. Interestingly, genes from the KstR2 regulon, which controls the degradation of cholesterol C and D rings, were significantly upregulated in all strains post-ALF exposure. These results indicate that M.tb -ALF contact drives initial metabolic and physiologic changes in M.tb , with potential implications in infection outcome.

IMPORTANCE

Tuberculosis, caused by airborne pathogen Mycobacterium tuberculosis ( M.tb ), is one of the leading causes of mortality worldwide. Upon infection, M.tb reaches the alveoli and gets in contact with human alveolar lining fluid (ALF), where ALF hydrolases modify the M.tb cell envelope driving subsequent M.tb -host cell interactions. Still, the contributions of ALF during infection are poorly understood. We exposed 4 M.tb strains to ALF for 15-min and 12-h and performed RNA sequencing, demonstrating a temporal and strain-specific adaptation of M.tb to ALF. Interestingly, genes associated with cholesterol degradation were highly upregulated in all strains. This study shows for the first time that ALF drives global metabolic changes in M.tb during the initial stages of the infection, with potential implications in disease outcome. Biologically relevant networks and common and strain-specific bacterial determinants derived from this study could be further investigated as potential therapeutic candidates.

A half-century of research on tuberculosis: Successes and challenges

Great progress has been made over the past half-century, but TB remains a formidable global health problem, particularly in low- and middle-income countries. Understanding the mechanisms of pathogenesis and necessary and sufficient conditions for protection are critical. The need for inexpensive and sensitive point-of-care diagnostic tests for earlier detection of infection and disease, shorter and less-toxic drug regimens for drug-sensitive and -resistant TB, and a more effective vaccine than BCG is immense. New and better tools, greater support for international research, collaborations, and training will be required to dramatically reduce the burden of this devastating disease which still kills 1.6 million people annually.

A contrast set mining based approach for cancer subtype analysis

The task of detecting common and unique characteristics among different cancer subtypes is an important focus of research that aims to improve personalized therapies. Unlike current approaches mainly based on predictive techniques, our study aims to improve the knowledge about the molecular mechanisms that descriptively led to cancer, thus not requiring previous knowledge to be validated. Here, we propose an approach based on contrast set mining to capture high-order relationships in cancer transcriptomic data. In this way, we were able to extract valuable insights from several cancer subtypes in the form of highly specific genetic relationships related to functional pathways affected by the disease. To this end, we have divided several cancer gene expression databases by the subtype associated with each sample to detect which gene groups are related to each cancer subtype. To demonstrate the potential and usefulness of the proposed approach we have extensively analysed RNA-Seq gene expression data from breast, kidney, and colon cancer subtypes. The possible role of the obtained genetic relationships was further evaluated through extensive literature research, while its prognosis was assessed via survival analysis, finding gene expression patterns related to survival in various cancer subtypes. Some gene associations were described in the literature as potential cancer biomarkers while other results have been not described yet and could be a starting point for future research.

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.

Host-pathogen relationship in retreated tuberculosis with major rifampicin resistance-conferring mutations

Introduction

It is assumed that host defense systems eliminating the pathogen and regulating tissue damage make a strong impact on the outcome of tuberculosis (TB) disease and that these processes are affected by rifampicin (RIF) resistance-conferring mutations of Mycobacterium tuberculosis (Mtb). However, the host responses to the pathogen harboring different mutations have not been studied comprehensively in clinical settings. We analyzed clinico-epidemiological factors and blood transcriptomic signatures associated with major rpoB mutations conferring RIF resistance in a cohort study.

Methods

Demographic data were collected from 295 active pulmonary TB patients with treatment history in Hanoi, Vietnam. When recruited, drug resistance-conferring mutations and lineage-specific variations were identified using whole-genome sequencing of clinical Mtb isolates. Before starting retreatment, total RNA was extracted from the whole blood of HIV-negative patients infected with Mtb that carried either the rpoB H445Y or rpoB S450L mutation, and the total RNA was subjected to RNA sequencing after age-gender matching. The individual RNA expression levels in the blood sample set were also measured using real-time RT-PCR. Logistic and linear regression models were used to assess possible associations.

Results

In our cohort, rpoB S450L and rpoB H445Y were major RIF resistance-conferring mutations [32/87 (36.8%) and 15/87 (17.2%), respectively]. H445Y was enriched in the ancient Beijing genotype and was associated with nonsynonymous mutations of Rv1830 that has been reported to regulate antibiotic resilience. H445Y was also more frequently observed in genetically clustered strains and in samples from patients who had received more than one TB treatment episode. According to the RNA sequencing, gene sets involved in the interferon-γ and-α pathways were downregulated in H445Y compared with S450L. The qRT-PCR analysis also confirmed the low expression levels of interferon-inducible genes, including BATF2 and SERPING1, in the H445Y group, particularly in patients with extensive lesions on chest X-ray.

Discussion

Our study results showed that rpoB mutations as well as Mtb sublineage with additional genetic variants may have significant effects on host response. These findings strengthen the rationale for investigation of host-pathogen interactions to develop countermeasures against epidemics of drug-resistant TB.

Downregulation of monocyte miRNAs: implications for immune dysfunction and disease severity in drug-resistant tuberculosis

Background

Monocyte miRNAs govern both protective and pathological responses during tuberculosis (TB) through their differential expression and emerged as potent targets for biomarker discovery and host-directed therapeutics. Thus, this study examined the miRNA profile of sorted monocytes across the TB disease spectrum [drug-resistant TB (DR-TB), drug-sensitive TB (DS-TB), and latent TB] and in healthy individuals (HC) to understand the underlying pathophysiology and their regulatory mechanism.

Methods

We sorted total monocytes including three subsets (HLA-DR⁺CD14⁺, HLA-DR⁺CD14⁺CD16⁺, and HLA-DR⁺CD16⁺cells) from peripheral blood mononuclear cells (PBMCs) of healthy and TB-infected individuals through flow cytometry and subjected them to NanoString-based miRNA profiling.

Results

The outcome was the differential expression of 107 miRNAs particularly the downregulation of miRNAs in the active TB groups (both drug-resistant and drug-sensitive). The miRNA profile revealed differential expression signatures: i) decline of miR-548m in DR-TB alone, ii) decline of miR-486-3p in active TB but significant elevation only in LTB iii) elevation of miR-132-3p only in active TB (DR-TB and DS-TB) and iv) elevation of miR-150-5p in DR-TB alone. The directionality of functions mediated by monocyte miRNAs from Gene Set Enrichment Analysis (GSEA) facilitated two phenomenal findings: i) a bidirectional response between active disease (activation profile in DR-TB and DS-TB compared to LTB and HC) and latent infection (suppression profile in LTB vs HC) and ii) hyper immune activation in the DR-TB group compared to DS-TB.

Conclusion

Thus, monocyte miRNA signatures provide pathological clues for altered monocyte function, drug resistance, and disease severity. Further studies on monocyte miRNAs may shed light on the immune regulatory mechanism for tuberculosis.

Identification and control for the effects of bioinformatic globin depletion on human RNA-seq differential expression analysis

When profiling blood samples by RNA-sequencing (RNA-seq), RNA from haemoglobin (Hgb) can account for up to 70% of the transcriptome. Due to considerations of sequencing depth and power to detect biological variation, Hgb RNA is typically depleted prior to sequencing by hybridisation-based methods; an alternative approach is to deplete reads arising from Hgb RNA bioinformatically. In the present study, we compared the impact of these two approaches on the outcome of differential gene expression analysis performed using RNA-seq data from 58 human tuberculosis (TB) patient or contact whole blood samples-29 globin kit-depleted and 29 matched non-depleted-a subset of which were taken at TB diagnosis and at six months post-TB treatment from the same patient. Bioinformatic depletion of Hgb genes from the non-depleted samples (bioinformatic-depleted) substantially reduced library sizes (median = 57.24%) and fewer long non-coding, micro, small nuclear and small nucleolar RNAs were captured in these libraries. Profiling published TB gene signatures across all samples revealed inferior correlation between kit-depleted and bioinformatic-depleted pairs when the proportion of reads mapping to Hgb genes was higher in the non-depleted sample, particularly at the TB diagnosis time point. A set of putative "globin-fingerprint" genes were identified by directly comparing kit-depleted and bioinformatic-depleted samples at each timepoint. Two TB treatment response signatures were also shown to have decreased differential performance when comparing samples at TB diagnosis to six months post-TB treatment when profiled on the bioinformatic-depleted samples compared with their kit-depleted counterparts. These results demonstrate that failure to deplete Hgb RNA prior to sequencing has a negative impact on the sensitivity to detect disease-relevant gene expression changes even when bioinformatic removal is performed.

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.

The effect of infection risk on female blood transcriptomics

Defenses against pathogens can take on many forms. For instance, behavioral avoidance of diseased conspecifics is widely documented. Interactions with these infectious conspecifics can also, however, lead to physiological changes in uninfected animals, an effect that is much less well understood. These changes in behavior and physiology are particularly important to study in a reproductive context, where they can impact reproductive decisions and offspring quality. Here, we studied how an acute (3 h) exposure to an immune-challenged male affected female blood transcriptomics and behavior. We predicted that females paired with immune-challenged males would reduce eating and drinking behaviors (as avoidance behaviors) and that their blood would show activation of immune and stress responses. We used female Japanese quail as a study system because they have been shown to respond to male traits, in terms of their own physiology and egg investment. Only two genes showed significant differential expression due to treatment, including an increase in the threonine dehydrogenase (TDH) transcript, an enzyme important for threonine breakdown. However, hundreds of genes in pathways related to activation of immune responses showed coordinated up-regulation in females exposed to immune-challenged males. Suppressed pathways revealed potential changes to metabolism and reduced responsiveness to glucocorticoids. Contrary to our prediction, we found that females paired with immune-challenged males increased food consumption. Water consumption was not changed by treatment. These findings suggest that even short exposure to diseased conspecifics can trigger both behavioral and physiological responses in healthy animals.

Classification of tuberculosis-related programmed cell death-related patient subgroups and associated immune cell profiling

Background

Tuberculosis (TB) is the deadliest communicable disease in the world with the exception of the ongoing COVID-19 pandemic. Programmed cell death (PCD) patterns play key roles in the development and progression of many disease states such that they may offer value as effective biomarkers or therapeutic targets that can aid in identifying and treating TB patients.

Materials and methods

The Gene Expression Omnibus (GEO) was used to gather TB-related datasets after which immune cell profiles in these data were analyzed to examine the potential TB-related loss of immune homeostasis. Profiling of differentially expressed PCD-related genes was performed, after which candidate hub PCD-associated genes were selected via a machine learning approach. TB patients were then stratified into two subsets based on the expression of PCD-related genes via consensus clustering. The potential roles of these PCD-associated genes in other TB-related diseases were further examined.

Results

In total, 14 PCD-related differentially expressed genes (DEGs) were identified and highly expressed in TB patient samples and significantly correlated with the abundance of many immune cell types. Machine learning algorithms enabled the selection of seven hub PCD-related genes that were used to establish PCD-associated patient subgroups, followed by the validation of these subgroups in independent datasets. These findings, together with GSVA results, indicated that immune-related pathways were significantly enriched in TB patients exhibiting high levels of PCD-related gene expression, whereas metabolic pathways were significantly enriched in the other patient group. Single cell RNA-seq (scRNA-seq) further highlighted significant differences in the immune status of these different TB patient samples. Furthermore, we used CMap to predict five potential drugs for TB-related diseases.

Conclusion

These results highlight clear enrichment of PCD-related gene expression in TB patients and suggest that this PCD activity is closely associated with immune cell abundance. This thus indicates that PCD may play a role in TB progression through the induction or dysregulation of an immune response. These findings provide a foundation for further research aimed at clarifying the molecular drivers of TB, the selection of appropriate diagnostic biomarkers, and the design of novel therapeutic interventions aimed at treating this deadly infectious disease.

Progress on diagnosis and treatment of latent tuberculosis infection

One fourth of the global population has been infected with Mycobacterium tuberculosis, and about 5%-10% of the infected individuals with latent tuberculosis infection (LTBI) will convert to active tuberculosis (ATB). Correct diagnosis and treatment of LTBI are important in ending the tuberculosis epidemic. Current methods for diagnosing LTBI, such as tuberculin skin test (TST) and interferon-γ release assay (IGRA), have limitations. Some novel biomarkers, such as transcriptome derived host genes in peripheral blood cells, will help to distinguish LTBI from ATB. More emphasis should be placed on surveillance in high-risk groups, including patients with HIV infection, those using biological agents, organ transplant recipients and those in close contact with ATB patients. For those with LTBI, treatment should be based on the risk of progression to ATB and the potential benefit. Prophylactic LTBI regimens include isoniazid monotherapy for 6 or 9 months, rifampicin monotherapy for 4 months, weekly rifapentine plus isoniazid for 3 months (3HP regimen) and daily rifampicin plus isoniazid for 3 months (3HR regimen). The success of the one month rifapentine plus isoniazid daily regimen (1HP regimen) suggests the feasibility of an ultra-short treatment strategy although its efficacy needs further assessment. Prophylactic treatment of LTBI in close contact with MDR-TB patients is another challenge, and the regimens include new anti-tuberculosis drugs such as bedaquiline, delamanid, fluoroquinolone and their combinations, which should be carefully evaluated. This article summarizes the current status of diagnosis and treatment of LTBI and its future development direction.

A Persistent Tuberculosis Outbreak in the UK Is Characterized by Hydrophobic fadB4-Deficient Mycobacterium tuberculosis That Replicates Rapidly in Macrophages

The genetic diversity of Mycobacterium tuberculosis can influence disease severity and transmissibility. To better understand how this diversity influences individuals and communities, we phenotyped M. tuberculosis that was causing a persistent outbreak in the East Midlands, United Kingdom. Compared to nonoutbreak isolates, bacilli had higher lipid contents and more hydrophobic cell surfaces. In macrophage infection models, the bacteria increased more rapidly, provoked the enhanced accumulation of macrophage lipid droplets and enhanced the secretion of IL-1β. Natural deletions in fadB4, nrdB, and plcC distinguished the outbreak isolates from other lineage 3 isolates in the region. fadB4 is annotated with a putative role in cell envelope biosynthesis, so the loss of this gene has the potential to alter the interactions of bacteria with immune cells. Reintroduction of fadB4 to the outbreak strain led to a phenotype that more closely resembled those of nonoutbreak strains. The improved understanding of the microbiological characteristics and the corresponding genetic polymorphisms that associate with outbreaks have the potential to inform tuberculosis control. IMPORTANCE Tuberculosis (TB) killed 1.5 million people in 2020 and affects every country. The extent to which the natural genetic diversity of Mycobacterium tuberculosis influences disease manifestation at both the individual and epidemiological levels remains poorly understood. Insights into how pathogen polymorphisms affect patterns of TB have the potential to translate into clinical and public health practice. Two distinct lineage 3 strains isolated from local TB outbreaks, one of which (CH) was rapidly terminated and the other of which (Lro) persistently transmitted for over a decade, provided us with an opportunity to study these issues. We compared genome sequences, microbiological characteristics, and early immune responses that were evoked upon infection. Our results indicate that the natural lack of fadB4 in the Lro strain contributes to its unique features.

Highlights of the ERS Lung Science Conference 2022

This article presents the highlights of the ERS Lung Science Conference 2022, including a session organised by the Early Career Member Committee (ECMC) dedicated to career development https://bit.ly/3tarCXc.

Gene expression profiling identifies candidate biomarkers for new latent tuberculosis infections. A cohort study

Objective

To determine the gene expression profile in individuals with new latent tuberculosis infection (LTBI), and to compare them with people with active tuberculosis (TB) and those exposed to TB but not infected.

Design

A prospective cohort study. Recruitment and follow-up were conducted between September 2016 to December 2018. Gene expression and data processing and analysis from April 2019 to April 2021.

Setting

Two male Colombian prisons.

Participants

15 new tuberculin skin test (TST) converters (negative TST at baseline that became positive during follow-up), 11 people that continued with a negative TST after two years of follow-up, and 10 people with pulmonary ATB.

Main outcome measures

Gene expression profile using RNA sequencing from PBMC samples. The differential expression was assessed using the DESeq2 package in Bioconductor. Genes with |logFC| >1.0 and an adjusted p-value < 0.1 were differentially expressed. We analyzed the differences in the enrichment of KEGG pathways in each group using InterMiner.

Results

The gene expression was affected by the time of incarceration. We identified group-specific differentially expressed genes between the groups: 289 genes in people with a new LTBI and short incarceration (less than three months of incarceration), 117 in those with LTBI and long incarceration (one or more years of incarceration), 26 in ATB, and 276 in the exposed but non-infected individuals. Four pathways encompassed the largest number of down and up-regulated genes among individuals with LTBI and short incarceration: cytokine signaling, signal transduction, neutrophil degranulation, and innate immune system. In individuals with LTBI and long incarceration, the only enriched pathway within up-regulated genes was Emi1 phosphorylation.

Conclusions

Recent infection with MTB is associated with an identifiable RNA pattern related to innate immune system pathways that can be used to prioritize LTBI treatment for those at greatest risk for developing active TB.

Whole transcriptome sequencing reveals neutrophils’ transcriptional landscape associated with active tuberculosis

Neutrophils have been recognized to play an important role in the pathogenesis of tuberculosis in recent years. Interferon-induced blood transcriptional signatures in ATB are predominantly driven by neutrophils. In this study, we performed global RNA-seq on peripheral blood neutrophils from active tuberculosis patients (ATB, n=15); latent tuberculosis infections (LTBI, n=22); and healthy controls (HC, n=21). The results showed that greater perturbations of gene expression patterns happened in neutrophils from ATB individuals than HC or those with LTBI, and a total of 344 differentially expressed genes (DEGs) were observed. Functional enrichment analysis showed that besides the interferon signaling pathway, multiple pattern recognition receptor pathways were significantly activated in ATB, such as NOD-like receptors and Toll-like receptors. Meanwhile, we also observed that the expression of genes related to endocytosis, secretory granules, and neutrophils degranulation were downregulated. Our data also showed that the NF-κB signaling pathway might be inhibited in patients with ATB, which could increase Mycobacterium tuberculosis survival and lead to active tuberculosis status. Furthermore, we validated the accuracy of some differentially expressed genes in an independent cohort using quantitative PCR, and obtained three novel genes (RBM3, CSRNP1, SRSF5) with the ability to discriminate active tuberculosis from LTBI and HC.

Contribution and Future of High-Throughput Transcriptomics in Battling Tuberculosis

While Tuberculosis (TB) infection remains a serious challenge worldwide, big data and “omic” approaches have greatly contributed to the understanding of the disease. Transcriptomics have been used to tackle a wide variety of queries including diagnosis, treatment evolution, latency and reactivation, novel target discovery, vaccine response or biomarkers of protection. Although a powerful tool, the elevated cost and difficulties in data interpretation may hinder transcriptomics complete potential. Technology evolution and collaborative efforts among multidisciplinary groups might be key in its exploitation. Here, we discuss the main fields explored in TB using transcriptomics, and identify the challenges that need to be addressed for a real implementation in TB diagnosis, prevention and therapy.

It Takes a Village: The Multifaceted Immune Response to Mycobacterium tuberculosis Infection and Vaccine-Induced Immunity

Despite co-evolving with humans for centuries and being intensely studied for decades, the immune correlates of protection against Mycobacterium tuberculosis (Mtb) have yet to be fully defined. This lapse in understanding is a major lag in the pipeline for evaluating and advancing efficacious vaccine candidates. While CD4+ T helper 1 (TH1) pro-inflammatory responses have a significant role in controlling Mtb infection, the historically narrow focus on this cell population may have eclipsed the characterization of other requisite arms of the immune system. Over the last decade, the tuberculosis (TB) research community has intentionally and intensely increased the breadth of investigation of other immune players. Here, we review mechanistic preclinical studies as well as clinical anecdotes that suggest the degree to which different cell types, such as NK cells, CD8+ T cells, γ δ T cells, and B cells, influence infection or disease prevention. Additionally, we categorically outline the observed role each major cell type plays in vaccine-induced immunity, including Mycobacterium bovis bacillus Calmette-Guérin (BCG). Novel vaccine candidates advancing through either the preclinical or clinical pipeline leverage different platforms (e.g., protein + adjuvant, vector-based, nucleic acid-based) to purposefully elicit complex immune responses, and we review those design rationales and results to date. The better we as a community understand the essential composition, magnitude, timing, and trafficking of immune responses against Mtb, the closer we are to reducing the severe disease burden and toll on human health inflicted by TB globally.

Clinical predictors of pulmonary tuberculosis among South African adults with HIV

BACKGROUND

Tuberculosis (TB) clinical prediction rules rely on presence of symptoms, however many undiagnosed cases in the community are asymptomatic. This study aimed to explore the utility of clinical factors in predicting TB among people with HIV not seeking care.

METHODS

Baseline data were analysed from an observational cohort of ambulant adults with HIV in South Africa. Participants were tested for Mycobacterium tuberculosis (Mtb) sensitisation (interferon-γ release assay, IGRA) and microbiologically-confirmed prevalent pulmonary TB disease at baseline, and actively surveilled for incident TB through 15 months. Multivariable LASSO regression with post-selection inference was used to test associations with Mtb sensitisation and TB disease.

FINDINGS

Between March 22, 2017, and May 15, 2018, 861 participants were enrolled; Among 851 participants included in the analysis, 94·5% were asymptomatic and 45·9% sensitised to Mtb. TB prevalence was 2·0% at baseline and incidence 2·3/100 person-years through 15 months follow-up. Study site was associated with baseline Mtb sensitisation (p < 0·001), prevalent (p < 0·001), and incident TB disease (p = 0·037). Independent of site, higher CD4 counts (per 50 cells/mm³, aOR 1·48, 95%CI 1·12-1·77, p = 0·006) were associated with increased IGRA positivity, and participants without TB disease (aOR 0·80, 95%CI 0·69-0·94, p = 0·006) had reduced IGRA positivity; no variables were independently associated with prevalent TB. Mixed ancestry (aHR 1·49, 95%CI 1·30->1000, p = 0·005) and antiretroviral initiation (aHR 1·48, 95%CI 1·01-929·93, p = 0·023) were independently associated with incident TB. Models incorporating clinical features alone performed poorly in diagnosing prevalent (AUC 0·65, 95%CI 0·44-0·85) or predicting progression to incident (0·67, 0·46-0·88) TB.

INTERPRETATION

CD4 count and antiretroviral initiation, proxies for immune status and HIV stage, were associated with Mtb sensitisation and TB disease. Inadequate performance of clinical prediction models may reflect predominantly subclinical disease diagnosed in this setting and unmeasured local site factors affecting transmission and progression.

FUNDING

The CORTIS-HR study was funded by the Bill & Melinda Gates Foundation (OPP1151915) and by the Strategic Health Innovation Partnerships Unit of the South African Medical Research Council with funds received from the South African Department of Science and Technology. The regulatory sponsor was the University of Cape Town.

Longitudinal Variations of M. tuberculosis-Induced IFN-γ Responses in HIV-Negative Pregnant Women Exposed to Tuberculosis

Introduction

Pregnancy triggers an alteration of the immune functions and increases the risk of developing the active tuberculosis (TB) symptoms in exposed women. The effect of pregnancy on the Mycobacterium tuberculosis-specific immune responses used for most of the TB immunodiagnostic assays is not well documented. Here we investigated the changes in the M. tuberculosis-specific IFN-γ production in age-matched pregnant and non-pregnant women according to their TB exposition status.

Methods

We conducted a prospective cohort study on HIV-seronegative pregnant and non-pregnant women with compatible pulmonary TB symptoms addressed to TB healthcare facilities in Antananarivo, Madagascar. Active pulmonary TB was bacteriologically assessed with culture from sputum samples. Clinical data and blood samples were collected at inclusion and after 6 months of follow-up for each individual included. Whole blood samples were stimulated with QuantiFERON TB-Gold Plus (QFT-P) assay antigens. Plasma IFN-γ concentrations were then assessed by ELISA.

Results

A total of 284 women were investigated for the study including 209 pregnant women without confirmed TB (pNTB), 24 pregnant women with bacteriologically confirmed active TB (pATB), 16 non-pregnant women with active TB (ATB), and 35 non-pregnant healthy donors (HC). At inclusion, IFN-γ responses are lower in the pregnant women compared to their age-matched non-pregnant counterparts and independently of their TB status. Among the pregnant women, higher concentrations of M. tuberculosis-specific IFN-γ were observed in those exposed to TB, but with a lower magnitude in the active TB compared to the latently infected pregnant women (p < 0.05 with TB1 and p < 0.01 with TB2). After 6 months of follow-up, the M. tuberculosis-specific IFN-γ responses return to their baseline concentrations except for the pregnant women treated for TB for which none of the QFT-P positive reversed to negative (0%, 0/10) at the end of their TB treatment.

Conclusion

These results support the concept of specific immune priorities characterized by a concomitant reduction in inflammatory immunity during pregnancy and corroborate the important role of activating the M. tuberculosis-specific immune responses to control the infection when the pregnant women are exposed to the pathogen.

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.

Single-cell profiling reveals distinct subsets of CD14+ monocytes drive blood immune signatures of active tuberculosis

Introduction

Previous studies suggest that monocytes are an important contributor to tuberculosis (TB)-specific immune signatures in blood.

Methods

Here, we carried out comprehensive single-cell profiling of monocytes in paired blood samples of active TB (ATB) patients at diagnosis and mid-treatment, and healthy controls.

Results

At diagnosis, ATB patients displayed increased monocyte-to-lymphocyte ratio, increased frequency of CD14+CD16- and intermediate CD14+CD16+ monocytes, and upregulation of interferon signaling genes that significantly overlapped with previously reported blood TB signatures in both CD14+ subsets. In this cohort, we identified additional transcriptomic and functional changes in intermediate CD14+CD16+ monocytes, such as the upregulation of inflammatory and MHC-II genes, and increased capacity to activate T cells, reflecting overall increased activation in this population. Single-cell transcriptomics revealed that distinct subsets of intermediate CD14+CD16+ monocytes were responsible for each gene signature, indicating significant functional heterogeneity within this population. Finally, we observed that changes in CD14+ monocytes were transient, as they were no longer observed in the same ATB patients mid-treatment, suggesting they are associated with disease resolution.

Discussion

Together, our study demonstrates for the first time that both intermediate and classical monocytes individually contribute to blood immune signatures of ATB and identifies novel subsets and associated gene signatures that may hold disease relevance.

Signature required: The transcriptional response to tuberculosis

The majority of humans infected with Mycobacterium tuberculosis never experience clinical symptoms or signs, but predicting those who will remains out of reach. Here, we discuss recent studies that reveal patterns and pathways that determine who is at highest risk for progression.