Profiling the immune response to Mycobacterium tuberculosis Beijing family infection: a perspective from the transcriptome

Effect of IL-17A on the immune response to pulmonary tuberculosis induced by high- and low-virulence strains of Mycobacterium bovis

Tuberculosis (TB) is an infectious, chronic, and progressive disease occurring globally. Human TB is caused mainly by Mycobacterium tuberculosis (M. tuberculosis), while the main causative agent of bovine TB is Mycobacterium bovis (M. bovis). The latter is one of the most important cattle pathogens and is considered the main cause of zoonotic TB worldwide. The mechanisms responsible for tissue damage (necrosis) during post-primary TB remain elusive. Recently, IL-17A was reported to be important for protection against M. tuberculosis infection, but it is also related to the production of an intense inflammatory response associated with necrosis. We used two M. bovis isolates with different levels of virulence and high IL-17A production to study this important cytokine's contrasting functions in a BALB/c mouse model of pulmonary TB. In the first part of the study, the gene expression kinetics and cellular sources of IL-17A were determined by real time PCR and immunohistochemistry respectively. Non-infected lungs showed low production of IL-17A, particularly by the bronchial epithelium, while lungs infected with the low-virulence 534 strain showed high IL-17A expression on Day 3 post-infection, followed by a decrease in expression in the early stage of the infection and another increase during late infection, on Day 60, when very low bacillary burdens were found. In contrast, infection with the highly virulent strain 04-303 induced a peak of IL-17A expression on Day 14 of infection, 1 week before extensive pulmonary necrosis was seen, being lymphocytes and macrophages the most important sources. In the second part of the study, the contribution of IL-17A to immune protection and pulmonary necrosis was evaluated by suppressing IL-17A via the administration of specific blocking antibodies. Infection with M. bovis strain 534 and treatment with IL-17A neutralizing antibodies did not affect mouse survival but produced a significant increase in bacillary load and a non-significant decrease in inflammatory infiltrate and granuloma area. In contrast, mice infected with the highly virulent 04-303 strain and treated with IL-17A blocking antibodies showed a significant decrease in survival, an increase in bacillary loads on Day 24 post-infection, and significantly more and earlier necrosis. Our results suggest that high expression of IL-17A is more related to protection than necrosis in a mouse model of pulmonary TB induced by M. bovis strains.

Backtranslation of human RNA biosignatures of tuberculosis disease risk into the preclinical pipeline is condition dependent

It is not clear whether human progression to active tuberculosis disease (TB) risk signatures are viable endpoint criteria for evaluations of treatments in clinical or preclinical development. TB is the deadliest infectious disease globally and more efficacious vaccines are needed to reduce this mortality. However, the immune correlates of protection for either preventing infection with Mycobacterium tuberculosis or preventing TB disease have yet to be completely defined, making the advancement of candidate vaccines through the pipeline slow, costly, and fraught with risk. Human-derived correlate of risk (COR) gene signatures, which identify an individual's risk to progressing to active TB disease, provide an opportunity for evaluating new therapies for TB with clear and defined endpoints. Though prospective clinical trials with longitudinal sampling are prohibitively expensive, characterization of COR gene signatures is practical with preclinical models. Using a 3Rs (Replacement, Reduction and Refinement) approach we reanalyzed heterogeneous publicly available transcriptional datasets to determine whether a specific set of COR signatures are viable endpoints in the preclinical pipeline. We selected RISK6, Sweeney3 and BATF2 human-derived blood-based RNA biosignatures because they require relatively few genes to assign a score and have been carefully evaluated across several clinical cohorts. Excitingly, these data provide proof-of-concept that human COR signatures seem to have high fidelity across several tissue types in the preclinical TB model pipeline and show best performance when the model most closely reflected human infection or disease conditions. Human-derived COR signatures offer an opportunity for high-throughput preclinical endpoint criteria of vaccine and drug therapy evaluations.

One Sentence Summary

Human-derived biosignatures of tuberculosis disease progression were evaluated for their predictive fidelity across preclinical species and derived tissues using available public data sets.

Bifurcation analysis of a tuberculosis progression model for drug target identification

Tuberculosis (TB) is a major cause of morbidity and mortality worldwide. The emergence and rapid spread of drug-resistant M. tuberculosis strains urge us to develop novel treatments. Experimental trials are constrained by laboratory capacity, insufficient funds, low number of laboratory animals and obsolete technology. Systems-level approaches to quantitatively study TB can overcome these limitations. Previously, we proposed a mathematical model describing the key regulatory mechanisms underlying the pathological progression of TB. Here, we systematically explore the effect of parameter variations on disease outcome. We find five bifurcation parameters that steer the clinical outcome of TB: number of bacteria phagocytosed per macrophage, macrophages death, macrophage killing by bacteria, macrophage recruitment, and phagocytosis of bacteria. The corresponding bifurcation diagrams show all-or-nothing dose-response curves with parameter regions mapping onto bacterial clearance, persistent infection, or history-dependent clearance or infection. Importantly, the pathogenic stage strongly affects the sensitivity of the host to these parameter variations. We identify parameter values corresponding to a latent-infection model of TB, where disease progression occurs significantly slower than in progressive TB. Two-dimensional bifurcation analyses uncovered synergistic parameter pairs that could act as efficient compound therapeutic approaches. Through bifurcation analysis, we reveal how modulation of specific regulatory mechanisms could steer the clinical outcome of TB.

The Detection of Mutations and Genotyping of Drug-Resistant Mycobacterium tuberculosis Strains Isolated from Patients in the Rural Eastern Cape Province

Drug-resistant tuberculosis (DR-TB) is still a major public health concern in South Africa. Mutations in M. tuberculosis can cause varying levels of phenotypic resistance to anti-TB medications. There have been no prior studies on gene mutations and the genotyping of DR-TB in the rural Eastern Cape Province; hence, we aimed to identify DR-TB mutations, genetic diversity, and allocated lineages among patients in this area. Using Xpert^® MTB/RIF, we assessed the rifampin resistance of sputum samples collected from 1157 patients suspected of having tuberculosis. GenoType MTBDR plus VER 2.0 was used for the detection of mutations causing resistance to anti-TB medications. The next step was to spoligotype 441 isolates. The most prevalent rifampin resistance-conferring mutations were in rpoB codon S531L in INH-resistant strains; the katG gene at codon S315TB and the inhA gene at codon C-15TB had the most mutations; 54.5% and 24.7%, respectively. In addition, 24.6% of strains showed mutations in both the rpoB and inhA genes, while 69.9% of strains showed mutations in both the katG and rpoB genes. Heteroresistance was seen in 17.9% of all cases in the study. According to spoligotyping analysis, Beijing families predominated. Investigation of the evolutionary lineages of M. tuberculosis isolates can be carried out using the information provided by the study's diversity of mutations. In locations wherein these mutations have been discovered, decision-making regarding the standardization of treatment regimens or individualized treatment may be aided by the detection frequency of rpoB, katG, and inhA mutations in various study areas.

Clinical isolates of Mycobacterium tuberculosis with different genotypes exhibit distinct host macrophage responses in vitro

Introduction. Mycobacterium tuberculosis (M.tb), the causative agent of tuberculosis, can survive as an intracellular parasite after entering macrophages via phagocytosis. M.tb strains are genotypically distinct and engage in diverse pathogen-host interactions, with different host immune responses triggered by different M.tb strains. Importantly, differences in intracellular accumulation and triggering of host macrophage responses during early infection stages are key determinants that shape the final outcomes of host innate immune responses to different M.tb strains.Hypothesis/Gap Statement. Clinical M.tb strains with different genotypes elicit different host innate immune responses in vitro.Aim. This work aimed to compare host innate immune responses elicited by genotypically diverse, clinically derived M.tb strains in vitro.Methodology. RAW264.7 cells were infected with three lineage 2 and lineage 4 clinically derived M.tb strains and strain H37Rv. Strains were evaluated for differences in intracellular growth, induction of macrophage apoptosis, and induction of expression of proinflammatory cytokines and associated pattern recognition receptors.Results. Highly variable cytokine profiles were observed subsequent to RAW264.7 cell infection with the different strains. The Beijing genotype strain, a modern Beijing strain belonging to lineage 2, induced milder host proinflammatory responses and less apoptosis and exhibited greater intracellular growth as compared to the other strains. Moreover, mRNA expression levels of iNOS in Beijing and MANU2 genotype strains exceeded corresponding levels obtained for the T1 genotype strain. Meanwhile, mRNA expression levels of toll-like receptor (TLR)-encoding genes TLR2 and TLR7 in macrophages infected with the Beijing genotype strain were higher than corresponding levels observed in MANU2 genotype strain-infected macrophages.Conclusion. The higher intracellular survival rate and lower level of host cell apoptosis associated with macrophage infection with the Beijing genotype strain indicated greater virulence of this strain relative to that of the other strains. Furthermore, in vitro immune responses induced by the Beijing genotype strain were unique in that this strain induced a weaker inflammatory response than was induced by T1 or MANU2 genotype strains. Nevertheless, additional evidence is needed to confirm that Beijing genotype strains possess greater virulence than strains with other genotypes.

Do B-1 cells play a role in response to Mycobacterium tuberculosis Beijing lineages?

We highlight the need to include an analysis of the B-1 B cell subset to complement the characterization of the cell-mediated immune response to the Mycobacterium tuberculosis Beijing lineage. The literature describes the B-1 cell repertoire's involvement in the cell-mediated response within granulomas, which is different from the classic antibody response B cells are generally associated with. Specifically, the B-1 B cell subset migrates from other compartments along with other cells to the infection site. We provide details to complement the reported results from Cerezo-Cortes et al.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

OBJECTIVE

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

METHODS

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

RESULTS

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

CONCLUSION

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Multidrug-Resistant Tuberculosis—Diagnostic Procedures and Treatment of Two Beijing-like TB Cases

The Beijing/W genotype is one of the major molecular families of Mycobacterium tuberculosis complex (MTBC), responsible for approximately 50% of tuberculosis (TB) cases in Far East Asia and at least 25% of TB cases globally. Studies have revealed that the Beijing genotype family is associated with a more severe clinical course of TB, increased ability to spread compared to other genotypes, and an unpredictable response to treatment. Based on the profile of spacers 35–43 in the Direct Repeat (DR) locus of the MTBC genome determined by spoligotyping, classical (typical) and modern (Beijing-like) clones can be identified within the Beijing family. While the modern and ancient Beijing strains appear to be closely related at the genetic level, there are marked differences in their drug resistance, as well as their ability to spread and cause disease. This paper presents two cases of drug-resistant tuberculosis caused by rare mycobacteria from the Beijing family: the Beijing 265 and Beijing 541 subtypes. The genotypes of isolated strains were linked with the clinical course of TB, and an attempt was made to initially assess whether the Beijing subtype can determine treatment outcomes in patients.

Mycobacterium tuberculosis infection drives a type I IFN signature in lung lymphocytes

Mycobacterium tuberculosis (Mtb) infects 25% of the world's population and causes tuberculosis (TB), which is a leading cause of death globally. A clear understanding of the dynamics of immune response at the cellular level is crucial to design better strategies to control TB. We use the single-cell RNA sequencing approach on lung lymphocytes derived from healthy and Mtb-infected mice. Our results show the enrichment of the type I IFN signature among the lymphoid cell clusters, as well as heat shock responses in natural killer (NK) cells from Mtb-infected mice lungs. We identify Ly6A as a lymphoid cell activation marker and validate its upregulation in activated lymphoid cells following infection. The cross-analysis of the type I IFN signature in human TB-infected peripheral blood samples further validates our results. These findings contribute toward understanding and characterizing the transcriptional parameters at a single-cell depth in a highly relevant and reproducible mouse model of TB.

Close Related Drug-Resistance Beijing Isolates of Mycobacterium tuberculosis Reveal a Different Transcriptomic Signature in a Murine Disease Progression Model

Mycobacterium tuberculosis (MTB) lineage 2/Beijing is associated with high virulence and drug resistance worldwide. In Colombia, the Beijing genotype has circulated since 1997, predominantly on the pacific coast, with the Beijing-Like SIT-190 being more prevalent. This genotype conforms to a drug-resistant cluster and shows a fatal outcome in patients. To better understand virulence determinants, we performed a transcriptomic analysis with a Beijing-Like SIT-190 isolate (BL-323), and Beijing-Classic SIT-1 isolate (BC-391) in progressive tuberculosis (TB) murine model. Bacterial RNA was extracted from mice lungs on days 3, 14, 28, and 60. On average, 0.6% of the total reads mapped against MTB genomes and of those, 90% against coding genes. The strains were independently associated as determined by hierarchical cluster and multidimensional scaling analysis. Gene ontology showed that in strain BL-323 enriched functions were related to host immune response and hypoxia, while proteolysis and protein folding were enriched in the BC-391 strain. Altogether, our results suggested a differential bacterial transcriptional program when evaluating these two closely related strains. The data presented here could potentially impact the control of this emerging, highly virulent, and drug-resistant genotype.