Alternate splicing of transcripts upon Mycobacterium tuberculosis infection impacts the expression of functional protein domains

The Role of mRNA Alternative Splicing in Macrophages Infected with Mycobacterium tuberculosis: A Field Needing to Be Discovered

Mycobacterium tuberculosis (Mtb) is one of the major causes of human death. In its battle with humans, Mtb has fully adapted to its host and developed ways to evade the immune system. At the same time, the human immune system has developed ways to respond to Mtb. The immune system responds to viral and bacterial infections through a variety of mechanisms, one of which is alternative splicing. In this study, we summarized the overall changes in alternative splicing of the transcriptome after macrophages were infected with Mtb. We found that after infection with Mtb, cells undergo changes, including (1) directly reducing the expression of splicing factors, which affects the regulation of gene expression, (2) altering the original function of proteins through splicing, which can involve gene truncation or changes in protein domains, and (3) expressing unique isoforms that may contribute to the identification and development of tuberculosis biomarkers. Moreover, alternative splicing regulation of immune-related genes, such as IL-4, IL-7, IL-7R, and IL-12R, may be an important factor affecting the activation or dormancy state of Mtb. These will help to fully understand the immune response to Mtb infection, which is crucial for the development of tuberculosis biomarkers and new drug targets.

Large-scale analysis reveals splicing biomarkers for tuberculosis progression and prognosis

BACKGROUND

Emerging evidence suggests that aberrant alternative splicing (AS) may play an important role in tuberculosis (TB). However, current knowledge regarding the value of AS in TB progression and prognosis remains unclear.

METHOD

Public RNA-seq datasets related to TB progression and prognosis were searched and AS analyses were conducted based on SUPPA2. Percent spliced in (PSI) was used for quantifying AS events and multiple machine learning (ML) methods were employed to construct predictive models. Area under curve (AUC), sensitivity and specificity were calculated to evaluate the model performance.

RESULTS

A total of 1587 samples from 7 datasets were included. Among 923 TB-progression related differential AS events (DASEs), 3 events (GET1-skipping exon (SE), TPD52-alternative first exons (AF) and TIMM10-alternative 5' splice site (A5)) were selected as candidate biomarkers; however, their predictive performance was limited. For TB prognosis, 5 events (PHF23-AF, KIF1B-SE, MACROD2-alternative 3' splice site (A3), CD55-retained intron (RI) and GALNT11-AF) were selected as candidates from the 1282 DASEs. Six ML methods were used to integrate these 5 events and XGBoost outperformed than others. AUC, sensitivity and specificity of XGBoost model were 0.875, 81.1% and 83.5% in training set, while they were 0.805, 68.4% and 73.2% in test set.

CONCLUSION

GET1-SE, TPD52-AF and TIMM10-A5 showed limited role in predicting TB progression, while PHF23-AF, KIF1B-SE, MACROD2-A3, CD55-RI and GALNT11-AF could well predict TB prognosis and work as candidate biomarkers. This work preliminarily explored the value of AS in predicting TB progression and prognosis and offered potential targets for further research.

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.

PRMT5 epigenetically regulates the E3 ubiquitin ligase ITCH to influence lipid accumulation during mycobacterial infection

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), triggers enhanced accumulation of lipids to generate foamy macrophages (FMs). This process has been often attributed to the surge in the expression of lipid influx genes with a concomitant decrease in those involved in lipid efflux. Here, we define an Mtb-orchestrated modulation of the ubiquitination of lipid accumulation markers to enhance lipid accretion during infection. We find that Mtb infection represses the expression of the E3 ubiquitin ligase, ITCH, resulting in the sustenance of key lipid accrual molecules viz. ADRP and CD36, that are otherwise targeted by ITCH for proteasomal degradation. In line, overexpressing ITCH in Mtb-infected cells was found to suppress Mtb-induced lipid accumulation. Molecular analyses including loss-of-function and ChIP assays demonstrated a role for the concerted action of the transcription factor YY1 and the arginine methyl transferase PRMT5 in restricting the expression of Itch gene by conferring repressive symmetrical H4R3me2 marks on its promoter. Consequently, siRNA-mediated depletion of YY1 or PRMT5 rescued ITCH expression, thereby compromising the levels of Mtb-induced ADRP and CD36 and limiting FM formation during infection. Accumulation of lipids within the host has been implicated as a pro-mycobacterial process that aids in pathogen persistence and dormancy. In line, we found that perturbation of PRMT5 enzyme activity resulted in compromised lipid levels and reduced mycobacterial survival in mouse peritoneal macrophages (ex vivo) and in a therapeutic mouse model of TB infection (in vivo). These findings provide new insights into the role of PRMT5 and YY1 in augmenting mycobacterial pathogenesis. Thus, we posit that our observations could help design novel adjunct therapies and combinatorial drug regimen for effective anti-TB strategies.

Mycobacterium tuberculosis infection leads to the formation of lipid-laden cells (foamy macrophages-FMs) that offer a favorable shelter for its persistence. During infection, we observe a significant reduction in the expression of the E3 ubiquitin ligase, ITCH. This repression allows the sustenance of key lipid accretion molecules (ADRP and CD36), by curbing their proteasomal degradation. Further, we show the repression of ITCH to be dependent on the concerted action of the bifunctional transcription factor, YY1 and the arginine methyl transferase, PRMT5. NOTCH signaling pathway was identified as a master-regulator of YY1 expression. In vitro and in vivo analyses revealed the significance of PRMT5 in regulating FM formation and consequently mycobacterial burden.

Past and Present Approaches to Diagnosis of Active Pulmonary Tuberculosis

Tuberculosis disease continues to contribute to the mortality burden globally. Due to the several shortcomings of the available diagnostic methods, tuberculosis disease continues to spread. The difficulty to obtain sputum among the very ill patients and the children also affects the quick diagnosis of tuberculosis disease. These challenges warrant investigating different sample types that can provide results in a short time. Highlighted in this review are the approved pulmonary tuberculosis diagnostic methods and ongoing research to improve its diagnosis. We used the PRISMA guidelines for systematic reviews to search for studies that met the selection criteria for this review. In this review we found out that enormous biosignature research is ongoing to identify host biomarkers that can be used as predictors of active PTB disease. On top of this, more research was also being done to improve already existing diagnostic tests. Host markers required more optimization for use in different settings given their varying sensitivity and specificity in PTB endemic and non-endemic settings.

Computational Methods to Study Human Transcript Variants in COVID-19 Infected Lung Cancer Cells

Microbes and viruses are known to alter host transcriptomes by means of infection. In light of recent challenges posed by the COVID-19 pandemic, a deeper understanding of the disease at the transcriptome level is needed. However, research about transcriptome reprogramming by post-transcriptional regulation is very limited. In this study, computational methods developed by our lab were applied to RNA-seq data to detect transcript variants (i.e., alternative splicing (AS) and alternative polyadenylation (APA) events). The RNA-seq data were obtained from a publicly available source, and they consist of mock-treated and SARS-CoV-2 infected (COVID-19) lung alveolar (A549) cells. Data analysis results show that more AS events are found in SARS-CoV-2 infected cells than in mock-treated cells, whereas fewer APA events are detected in SARS-CoV-2 infected cells. A combination of conventional differential gene expression analysis and transcript variants analysis revealed that most of the genes with transcript variants are not differentially expressed. This indicates that no strong correlation exists between differential gene expression and the AS/APA events in the mock-treated or SARS-CoV-2 infected samples. These genes with transcript variants can be applied as another layer of molecular signatures for COVID-19 studies. In addition, the transcript variants are enriched in important biological pathways that were not detected in the studies that only focused on differential gene expression analysis. Therefore, the pathways may lead to new molecular mechanisms of SARS-CoV-2 pathogenesis.

Interspecies RNA Interactome of Pathogen and Host in a Heritable Defensive Strategy

Communication with bacteria deeply impacts the life history traits of their hosts. Through specific molecules and metabolites, bacteria can promote short- and long-term phenotypic and behavioral changes in the nematode Caenorhabditis elegans. The chronic exposure of C. elegans to pathogens promotes the adaptive behavior in the host’s progeny called pathogen-induced diapause formation (PIDF). PIDF is a pathogen avoidance strategy induced in the second generation of animals infected and can be recalled transgenerationally. This behavior requires the RNA interference machinery and specific nematode and bacteria small RNAs (sRNAs). In this work, we assume that RNAs from both species co-exist and can interact with each other. Under this principle, we explore the potential interspecies RNA interactions during PIDF-triggering conditions, using transcriptomic data from the holobiont. We study two transcriptomics datasets: first, the dual sRNA expression of Pseudomonas aeruginosa PAO1 and C. elegans in a transgenerational paradigm for six generations and second, the simultaneous expression of sRNAs and mRNA in intergenerational PIDF. We focus on those bacterial sRNAs that are systematically overexpressed in the intestines of animals compared with sRNAs expressed in host-naïve bacteria. We selected diverse in silico methods that represent putative mechanisms of RNA-mediated interspecies interaction. These interactions are as follows: heterologous perfect and incomplete pairing between bacterial RNA and host mRNA; sRNAs of similar sequence expressed in both species that could mimic each other; and known or predicted eukaryotic motifs present in bacterial transcripts. We conclude that a broad spectrum of tools can be applied for the identification of potential sRNA and mRNA targets of the interspecies RNA interaction that can be subsequently tested experimentally.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case–control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Is targeting dysregulation in apoptosis splice variants in Mycobacterium tuberculosis (MTB) host interactions and splicing factors resulting in immune evasion by MTB strategies a possibility?

Mycobacterium tuberculosis (Mtb), is one of the foremost organisms causing mortality in humans, and has been for most of human history. When faced with an infection the human immune system is ordinarily very competent in killing both extracellular and intracellular bacilli. However, Mtb is able to evade the host immune system and is even able to establish a persistent infectious reservoir by "hiding" in the immune cells of the host. While the mechanisms by which the bacteria accomplishes this are not fully understood, it is known that the bacterium can subvert cellular processes in cells such as macrophages that prevent the lysis of the bacteria or the cell undergoing apoptosis. They are also able to interfere with immune cell signalling. One of the greatest effects that Mtb has is too alter the transcriptome of the macrophage. An easy way for the bacterium to accomplish this is to alter the alternative splicing patterns of the host. This can lead to a large change in the population of different protein isoforms, some of which have very different functions when compared to the original protein. At the same time the long history of Mtb infecting humans have led to specific immune reactions that occur in the host immune system in order to fight the infection. Many of these specific reactions involve new isoforms of host defence proteins. In this way the human host can use alternate splicing to create new isoforms of immune- related proteins that are more effective in defending against Mtb.

From Gene to Protein-How Bacterial Virulence Factors Manipulate Host Gene Expression During Infection

Bacteria evolved many strategies to survive and persist within host cells. Secretion of bacterial effectors enables bacteria not only to enter the host cell but also to manipulate host gene expression to circumvent clearance by the host immune response. Some effectors were also shown to evade the nucleus to manipulate epigenetic processes as well as transcription and mRNA procession and are therefore classified as nucleomodulins. Others were shown to interfere downstream with gene expression at the level of mRNA stability, favoring either mRNA stabilization or mRNA degradation, translation or protein stability, including mechanisms of protein activation and degradation. Finally, manipulation of innate immune signaling and nutrient supply creates a replicative niche that enables bacterial intracellular persistence and survival. In this review, we want to highlight the divergent strategies applied by intracellular bacteria to evade host immune responses through subversion of host gene expression via bacterial effectors. Since these virulence proteins mimic host cell enzymes or own novel enzymatic functions, characterizing their properties could help to understand the complex interactions between host and pathogen during infections. Additionally, these insights could propose potential targets for medical therapy.

The E-Cadherin Cleavage Associated to Pathogenic Bacteria Infections Can Favor Bacterial Invasion and Transmigration, Dysregulation of the Immune Response and Cancer Induction in Humans

Once bound to the epithelium, pathogenic bacteria have to cross epithelial barriers to invade their human host. In order to achieve this goal, they have to destroy the adherens junctions insured by cell adhesion molecules (CAM), such as E-cadherin (E-cad). The invasive bacteria use more or less sophisticated mechanisms aimed to deregulate CAM genes expression or to modulate the cell-surface expression of CAM proteins, which are otherwise rigorously regulated by a molecular crosstalk essential for homeostasis. Apart from the repression of CAM genes, a drastic decrease in adhesion molecules on human epithelial cells can be obtained by induction of eukaryotic endoproteases named sheddases or through synthesis of their own (prokaryotic) sheddases. Cleavage of CAM by sheddases results in the release of soluble forms of CAM. The overexpression of soluble CAM in body fluids can trigger inflammation and pro-carcinogenic programming leading to tumor induction and metastasis. In addition, the reduction of the surface expression of E-cad on epithelia could be accompanied by an alteration of the anti-bacterial and anti-tumoral immune responses. This immune response dysfunction is likely to occur through the deregulation of immune cells homing, which is controlled at the level of E-cad interaction by surface molecules α_E integrin (CD103) and lectin receptor KLRG1. In this review, we highlight the central role of CAM cell-surface expression during pathogenic microbial invasion, with a particular focus on bacterial-induced cleavage of E-cad. We revisit herein the rapidly growing body of evidence indicating that high levels of soluble E-cad (sE-cad) in patients’ sera could serve as biomarker of bacterial-induced diseases.

Gut Microbiota Has a Widespread and Modifiable Effect on Host Gene Regulation

The composition of the gut microbiome has been associated with various aspects of human health, but the mechanism of this interaction is still unclear. We utilized a cellular system to characterize the effect of the microbiome on human gene expression. We showed that some of these changes in expression may be mediated by changes in chromatin accessibility. Furthermore, we validate the role of a specific microbe and show that changes in its abundance can modify the host gene expression response. These results show an important role of gut microbiota in regulating host gene expression and suggest that manipulation of microbiome composition could be useful in future therapies.

Variation in gut microbiome is associated with wellness and disease in humans, and yet the molecular mechanisms by which this variation affects the host are not well understood. A likely mechanism is that of changing gene regulation in interfacing host epithelial cells. Here, we treated colonic epithelial cells with live microbiota from five healthy individuals and quantified induced changes in transcriptional regulation and chromatin accessibility in host cells. We identified over 5,000 host genes that change expression, including 588 distinct associations between specific taxa and host genes. The taxa with the strongest influence on gene expression alter the response of genes associated with complex traits. Using ATAC-seq, we showed that a subset of these changes in gene expression are associated with changes in host chromatin accessibility and transcription factor binding induced by exposure to gut microbiota. We then created a manipulated microbial community with titrated doses of Collinsella, demonstrating that manipulation of the composition of the microbiome under both natural and controlled conditions leads to distinct and predictable gene expression profiles in host cells. Taken together, our results suggest that specific microbes play an important role in regulating expression of individual host genes involved in human complex traits. The ability to fine-tune the expression of host genes by manipulating the microbiome suggests future therapeutic routes.

IMPORTANCE The composition of the gut microbiome has been associated with various aspects of human health, but the mechanism of this interaction is still unclear. We utilized a cellular system to characterize the effect of the microbiome on human gene expression. We showed that some of these changes in expression may be mediated by changes in chromatin accessibility. Furthermore, we validate the role of a specific microbe and show that changes in its abundance can modify the host gene expression response. These results show an important role of gut microbiota in regulating host gene expression and suggest that manipulation of microbiome composition could be useful in future therapies.

RNA Splicing: A New Paradigm in Host-Pathogen Interactions

RNA splicing brings diversity to the eukaryotic proteome. Different spliced variants of a gene may differ in their structure, function, localization, and stability influencing protein stoichiometry and physiological outcomes. Alternate spliced variants of different genes are known to associate with various chronic pathologies including cancer. Emerging evidence suggests precise regulation of splicing as fundamental to normal well-being. In this context, infection-induced alternative splicing has emerged as a new pivot of host function, which pathogenic microbes can alter-directly or indirectly-to tweak the host immune responses against the pathogen. The implications of these findings are vast, and although not explored much in the case of pathogenic infections, we present here examples from splicing mediated regulation of immune responses across a variety of conditions and explore how this fascinating finding brings a new paradigm to host-pathogen interactions.