Non-coding RNA and its potential role in Mycobacterium tuberculosis pathogenesis

A virulence-associated small RNA MTS1338 activates an ABC transporter CydC for rifampicin efflux in Mycobacterium tuberculosis

The efficacy of the tuberculosis treatment is restricted by innate drug resistance of Mycobacterial tuberculosis and its ability to acquire resistance to all anti-tuberculosis drugs in clinical use. A profound understanding of bacterial ploys that decrease the effectiveness of drugs would identify new mechanisms for drug resistance, which would subsequently lead to the development of more potent TB therapies. In the current study, we identified a virulence-associated small RNA (sRNA) MTS1338-driven drug efflux mechanism in M. tuberculosis. The treatment of a frontline antitubercular drug rifampicin upregulated MTS1338 by >4-fold. Higher intrabacterial abundance of MTS1338 increased the growth rate of cells in rifampicin-treated conditions. This fact was attributed by the upregulation of an efflux protein CydC by MTS1338. Gel-shift assay identified a stable interaction of MTS1338 with the coding region of cydC mRNA thereby potentially stabilizing it at the posttranscriptional level. The drug efflux measurement assays revealed that cells with higher MTS1338 abundance accumulate less drug in the cells. This study identified a new regulatory mechanism of drug efflux controlled by an infection-induced sRNA in M. tuberculosis.

Unveiling the orchestration: mycobacterial small RNAs as key mediators in host-pathogen interactions

Small RNA (sRNA) molecules, a class of non-coding RNAs, have emerged as pivotal players in the regulation of gene expression and cellular processes. Mycobacterium tuberculosis and other pathogenic mycobacteria produce diverse small RNA species that modulate bacterial physiology and pathogenesis. Recent advances in RNA sequencing have enabled identification of novel small RNAs and characterization of their regulatory functions. This review discusses the multifaceted roles of bacterial small RNAs, covering their biogenesis, classification, and functional diversity. Small RNAs (sRNAs) play pivotal roles in orchestrating diverse cellular processes, ranging from gene silencing to epigenetic modifications, across a broad spectrum of organisms. While traditionally associated with eukaryotic systems, recent research has unveiled their presence and significance within bacterial domains as well. Unlike their eukaryotic counterparts, which primarily function within the context of RNA interference (RNAi) pathways, bacterial sRNAs predominantly act through base-pairing interactions with target mRNAs, leading to post-transcriptional regulation. This fundamental distinction underscores the necessity of elucidating the unique roles and regulatory mechanisms of bacterial sRNAs in bacterial adaptation and survival. By doing these myriad functions, they regulate bacterial growth, metabolism, virulence, and drug resistance. In Mycobacterium tuberculosis, apart from having various roles in the bacillus itself, small RNA molecules have emerged as key regulators of gene expression and mediators of host-pathogen interactions. Understanding sRNA regulatory networks in mycobacteria can drive our understanding of significant role they play in regulating virulence and adaptation to the host environment. Detailed functional characterization of Mtb sRNAs at the host-pathogen interface is required to fully elucidate the complex sRNA-mediated gene regulatory networks deployed by Mtb, to manipulate the host. A deeper understanding of this aspect could pave the development of novel diagnostic and therapeutic strategies for tuberculosis.

Transcriptional activation of the Mycobacterium tuberculosis virulence-associated small RNA MTS1338 by the response regulators DosR and PhoP

MTS1338, a distinctive small RNA in pathogenic mycobacteria, plays a crucial role in host-pathogen interactions during infection. Mycobacterial cells encounter heterogeneous stresses in macrophages, which highly upregulate MTS1338. A dormancy regulatory factor DosR regulates the intracellular abundance of MTS1338. Herein, we investigated the interplay of DosR and a low pH-inducible gene regulator PhoP binding to the MTS1338 promoter. We identified that DosR strongly binds to two regions upstream of the MTS1338 gene. The proximal region possesses a threefold higher affinity than the distal site, but the presence of both regions increased the affinity for DosR by > 10-fold. PhoP did not bind to the MTS1338 gene but binds to the DosR-bound MTS1338 gene, suggesting a concerted mechanism for MTS1338 expression.

Mycobacterium tuberculosis complex molecular networks and their regulation: Implications of strain heterogeneity on epigenetic diversity and transcriptome regulation

Tuberculosis has been a public health crisis since the 1900, which has caused the highest mortalities due to a single bacterial infection worldwide, that was recently further complicated by the Coronavirus disease 2019 pandemic. The causative agent of Tuberculosis, Mycobacterium tuberculosis, belongs to a genetically well-characterized family of strains known as the Mycobacterium tuberculosis complex, which has complicated progress made towards eradicating Tuberculosis due to pathogen-specific phenotypic differences in the members of this complex. Mycobacterium tuberculosis complex strains are genetically diverse human- and animal-adapted pathogens belonging to 7 lineages (Indo-Oceanic, East-Asian, East-African Indian, Euro-American, M. africanum West Africa 1, M. africanum West Africa 2 and Ethopia), respectively and the recently identified Lineage 8 and M. africanum Lineage 9. Genomic studies have revealed that Mycobacterium tuberculosis complex members are ∼99 % similar, however, due to selective pressure and adaptation to human host, they are prone to mutations that have resulted in development of drug resistance and phenotypic heterogeneity that impact strain virulence. Furthermore, members of the Mycobacterium tuberculosis complex have preferred geographic locations and possess unique phenotypic characteristics that is linked to their pathogenicity. Due to the recent advances in development next generation sequencing platforms, several studies have revealed epigenetic changes in genomic regions combined with "unique" gene regulatory mechanisms through non-coding RNAs that are responsible for strain-specific behaviour on in vitro and in vivo infection models. The current review provides up to date epigenetic patterns, gene regulation through non-coding RNAs, together with implications of these mechanisms in down-stream proteome and metabolome, which may be responsible for "unique" responses to infection by members of the Mycobacterium tuberculosis complex. Understanding lineage-specific molecular mechanisms during infection may provide novel drug targets and disease control measures towards World Health organization END-TB strategy.

Small RNA MTS1338 Configures a Stress Resistance Signature in Mycobacterium tuberculosis

In the course of evolution, Mycobacterium tuberculosis (Mtb), the etiological agent of tuberculosis, has developed sophisticated strategies to evade host immune response, including the synthesis of small non-coding RNAs (sRNAs), which regulate post-transcriptional pathways involved in the stress adaptation of mycobacteria. sRNA MTS1338 is upregulated in Mtb during its infection of cultured macrophages and in the model of chronic tuberculosis, suggesting involvement in host-pathogen interactions. Here, we analyzed the role of MTS1338 in the Mtb response to macrophage-like stresses in vitro. The Mtb strain overexpressing MTS1338 demonstrated enhanced survival ability under low pH, nitrosative, and oxidative stress conditions simulating the antimicrobial environment inside macrophages. Transcriptomic analysis revealed that in MTS1338-overexpressing Mtb, the stress factors led to the activation of a number of transcriptional regulators, toxin-antitoxin modules, and stress chaperones, about half of which coincided with the genes induced in Mtb phagocytosed by macrophages. We determined the MTS1338 "core regulon", consisting of 11 genes that were activated in all conditions under MTS1338 overexpression. Our findings indicate that MTS1338 is a stress-induced sRNA that promotes Mtb survival in macrophages by triggering adaptive transcriptional mechanisms in response to host antimicrobial defense reactions.

Using a whole genome co-expression network to inform the functional characterisation of predicted genomic elements from Mycobacterium tuberculosis transcriptomic data

A whole genome co-expression network was created using Mycobacterium tuberculosis transcriptomic data from publicly available RNA-sequencing experiments covering a wide variety of experimental conditions. The network includes expressed regions with no formal annotation, including putative short RNAs and untranslated regions of expressed transcripts, along with the protein-coding genes. These unannotated expressed transcripts were among the best-connected members of the module sub-networks, making up more than half of the 'hub' elements in modules that include protein-coding genes known to be part of regulatory systems involved in stress response and host adaptation. This data set provides a valuable resource for investigating the role of non-coding RNA, and conserved hypothetical proteins, in transcriptomic remodelling. Based on their connections to genes with known functional groupings and correlations with replicated host conditions, predicted expressed transcripts can be screened as suitable candidates for further experimental validation.

Identification of RNU44 as an Endogenous Reference Gene for Normalizing Cell-Free RNA in Tuberculosis

Background

Normalization of cell-free RNA (cf-RNA) concentration can be affected by variable experimental conditions and thus impact the performance of their diagnostic potential. Our study aimed to identify appropriate endogenous reference genes for cf-RNA biomarker evaluation in the diagnosis of tuberculosis (TB).

Methods

Subjects consisting of patients with TB with and without malignancy, patients with pneumonia, and healthy controls were recruited. Candidate reference genes were screened and identified by literature reviewing and RNA-Seq analysis. Expression levels of the candidate genes were determined by reverse-transcription real-time quantitative polymerase chain reaction in plasma from patients with TB and healthy controls. The stability of gene expression was assessed by geNorm, NormFinder, BestKeeper, the Comparative Delta Ct method, and RefFinder. Differential expression of 2 small RNAs (sRNAs) encoding by genome of Mycobacterium tuberculosis in plasma of patients with TB were determined by both absolute quantification and relative quantification with candidate reference genes.

Results

According to the stability ranking analyzed with the 5 computational programs, the top 4 candidates—miR-93, RNU44, miR-16, and glyceraldehyde 3-phosphate dehydrogenase—were used to normalize the transcript levels of 2 mycobacterial sRNAs, MTS2823 and MTS1338, which were observed to have higher copy numbers in the plasma of patients with TB. Normalization with RNU44 displayed significantly higher levels of the 2 M tuberculosis sRNAs in the patients’ plasma than those of healthy controls.

Conclusions

RNU44 was demonstrated as a proper endogenous gene for cf-RNA normalization in TB diagnosis.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

Immune regulation and emerging roles of noncoding RNAs in Mycobacterium tuberculosis infection

Tuberculosis, caused by Mycobacterium tuberculosis, engenders an onerous burden on public hygiene. Congenital and adaptive immunity in the human body act as robust defenses against the pathogens. However, in coevolution with humans, this microbe has gained multiple lines of mechanisms to circumvent the immune response to sustain its intracellular persistence and long-term survival inside a host. Moreover, emerging evidence has revealed that this stealthy bacterium can alter the expression of demic noncoding RNAs (ncRNAs), leading to dysregulated biological processes subsequently, which may be the rationale behind the pathogenesis of tuberculosis. Meanwhile, the differential accumulation in clinical samples endows them with the capacity to be indicators in the time of tuberculosis suffering. In this article, we reviewed the nearest insights into the impact of ncRNAs during Mycobacterium tuberculosis infection as realized via immune response modulation and their potential as biomarkers for the diagnosis, drug resistance identification, treatment evaluation, and adverse drug reaction prediction of tuberculosis, aiming to inspire novel and precise therapy development to combat this pathogen in the future.

Transcriptional regulation and drug resistance in Mycobacterium tuberculosis

Bacterial drug resistance is one of the major challenges to present and future human health, as the continuous selection of multidrug resistant bacteria poses at serious risk the possibility to treat infectious diseases in the near future. One of the infection at higher risk to become incurable is tuberculosis, due to the few drugs available in the market against Mycobacterium tuberculosis. Drug resistance in this species is usually due to point mutations in the drug target or in proteins required to activate prodrugs. However, another interesting and underexplored aspect of bacterial physiology with important impact on drug susceptibility is represented by the changes in transcriptional regulation following drug exposure. The main regulators involved in this phenomenon in M. tuberculosis are the sigma factors, and regulators belonging to the WhiB, GntR, XRE, Mar and TetR families. Better understanding the impact of these regulators in survival to drug treatment might contribute to identify new drug targets and/or to design new strategies of intervention.

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.

Non-Coding RNAs in Tuberculosis Epidemiology: Platforms and Approaches for Investigating the Genome's Dark Matter

A growing amount of information about the different types, functions, and roles played by non-coding RNAs (ncRNAs) is becoming available, as more and more research is done. ncRNAs have been identified as potential therapeutic targets in the treatment of tuberculosis (TB), because they may be essential regulators of the gene network. ncRNA profiling and sequencing has recently revealed significant dysregulation in tuberculosis, primarily due to aberrant processes of ncRNA synthesis, including amplification, deletion, improper epigenetic regulation, or abnormal transcription. Despite the fact that ncRNAs may have a role in TB characteristics, the detailed mechanisms behind these occurrences are still unknown. The dark matter of the genome can only be explored through the development of cutting-edge bioinformatics and molecular technologies. In this review, ncRNAs' synthesis and functions are discussed in detail, with an emphasis on the potential role of ncRNAs in tuberculosis. We also focus on current platforms, experimental strategies, and computational analyses to explore ncRNAs in TB. Finally, a viewpoint is presented on the key challenges and novel techniques for the future and for a wide-ranging therapeutic application of ncRNAs.

Machine Learning of All Mycobacterium tuberculosis H37Rv RNA-seq Data Reveals a Structured Interplay between Metabolism, Stress Response, and Infection

Mycobacterium tuberculosis is one of the most consequential human bacterial pathogens, posing a serious challenge to 21st century medicine. A key feature of its pathogenicity is its ability to adapt its transcriptional response to environmental stresses through its transcriptional regulatory network (TRN). While many studies have sought to characterize specific portions of the M. tuberculosis TRN, and some studies have performed system-level analysis, few have been able to provide a network-based model of the TRN that also provides the relative shifts in transcriptional regulator activity triggered by changing environments. Here, we compiled a compendium of nearly 650 publicly available, high quality M. tuberculosis RNA-sequencing data sets and applied an unsupervised machine learning method to obtain a quantitative, top-down TRN. It consists of 80 independently modulated gene sets known as “iModulons,” 41 of which correspond to known regulons. These iModulons explain 61% of the variance in the organism’s transcriptional response. We show that iModulons (i) reveal the function of poorly characterized regulons, (ii) describe the transcriptional shifts that occur during environmental changes such as shifting carbon sources, oxidative stress, and infection events, and (iii) identify intrinsic clusters of regulons that link several important metabolic systems, including lipid, cholesterol, and sulfur metabolism. This transcriptome-wide analysis of the M. tuberculosis TRN informs future research on effective ways to study and manipulate its transcriptional regulation and presents a knowledge-enhanced database of all published high-quality RNA-seq data for this organism to date.

IMPORTANCEMycobacterium tuberculosis H37Rv is one of the world's most impactful pathogens, and a large part of the success of the organism relies on the differential expression of its genes to adapt to its environment. The expression of the organism's genes is driven primarily by its transcriptional regulatory network, and most research on the TRN focuses on identifying and quantifying clusters of coregulated genes known as regulons. While previous studies have relied on molecular measurements, in the manuscript we utilized an alternative technique that performs machine learning to a large data set of transcriptomic data. This approach is less reliant on hypotheses about the role of specific regulatory systems and allows for the discovery of new biological findings for already collected data. A better understanding of the structure of the M. tuberculosis TRN will have important implications in the design of improved therapeutic approaches.

Modulatory Impact of the sRNA Mcr11 in Two Clinical Isolates of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is a successful pathogen causing tuberculosis (TB) disease in humans. It has been shown, that some circulating strains of Mtb in TB endemic populations, are more virulent and more transmissible than others, which may be related to their evolved adaptations to modulate the host immune responses. Underlying these adaptations to the stressful conditions, different genetic regulatory networks involved sRNAs that are mostly unknown for Mtb. We have previously shown that Mcr11 is one of the main sRNAs that determine transcriptomic differences among the Colombian clinical isolates UT127 and UT205 compared to the laboratory strain H37Rv. We found that the knock-down of mcr11 using CRISPRi has a major impact on phenotypic traits, especially in the clinical isolate UT205. Through the analysis of RNA-seq during the knock-down of mcr11 in UT205, we found a downregulation of genes mainly involved in lipid synthesis, lipid metabolism, ribosomal proteins, transport systems, respiratory and energy systems, membrane and cell wall components, intermediary metabolism, lipoproteins and virulence genes. One of the most interesting genes showing transcriptomic changes is OprA (encoded by the gene rv0516c), which has been involved in the K⁺ regulation. Overall, our data may suggest that one of the prominent roles of the sRNA Mcr11 is to regulate genes that control Mtb growth and osmoregulation.

Challenges in defining the functional, non-coding, expressed genome of members of the Mycobacterium tuberculosis complex

A definitive transcriptome atlas for the non-coding expressed elements of the members of the Mycobacterium tuberculosis complex (MTBC) does not exist. Incomplete lists of non-coding transcripts can be obtained for some of the reference genomes (e.g., M. tuberculosis H37Rv) but to what extent these transcripts have homologues in closely related species or even strains is not clear. This has implications for the analysis of transcriptomic data; non-coding parts of the transcriptome are often ignored in the absence of formal, reliable annotation. Here, we review the state of our knowledge of non-coding RNAs in pathogenic mycobacteria, emphasizing the disparities in the information included in commonly used databases. We then proceed to review ways of combining computational solutions for predicting the non-coding transcriptome with experiments that can help refine and confirm these predictions.

MTS1338 in Mycobacterium tuberculosis promotes detoxification of reactive oxygen species under oxidative stress

Diverse mechanisms exist in Mycobacterium tuberculosis for adaptation to stresses leading to its persistence in the hostile environment of macrophages. Small RNA (sRNA)-mediated regulatory mechanisms have been scarcely explored in M. tuberculosis. MTS1338, a sRNA present only in pathogenic mycobacteria, was discovered to be highly abundant during infection and significantly contributes to host-pathogen interaction. A variety of stresses have been implicated for its accumulation. Herein, we showed that MTS1338 is an oxidative stress induced sRNA, which promotes the detoxification of reactive oxygen species (ROS) under oxidative stress. Current study identified a new role of MTS1338 in M. tuberculosis under oxidative stress.

Effect of Sequence on the Interactions of Divalent Cations with M-Box Riboswitches from Mycobacterium tuberculosis and Bacillus subtilis

RNA is highly negatively charged and often acquires complex structures that require the presence of divalent cations. Subtle changes in conformation resulting from changes in sequence can affect the way ions associate with RNA. Riboswitches are RNA molecules that are involved in the control of gene expression in bacteria and are excellent systems for testing the effects of sequence variations on the conformation of RNA because they contain a highly conserved binding pocket but present sequence variability among different organisms. In this work, we have compared the aptamer domain of a proposed M-box riboswitch from Mycobacterium tuberculosis with the aptamer domain of a validated M-box riboswitch from Bacillus subtilis. We have in vitro transcribed and purified wild-type (WT) M-box riboswitches from M. tuberculosis and B. subtilis as well as a variety of mutated aptamers in which regions from one riboswitch have been replaced with regions from the other riboswitch. We have used ultraviolet unfolding experiments and circular dichroism to characterize the interactions of WT and related M-box riboswitches with divalent cations. Our results show that M-box from M. tuberculosis associates with Mg²⁺ and Sr²⁺ in a similar fashion while M-box from B. subtilis discriminates between these two ions and appears to associate better with Mg²⁺. Our overall results show that M-box from M. tuberculosis interacts differently with cations than M-box from B. subtilis and suggest conformational differences between these two riboswitches.

Mycobacterium tuberculosis protein PPE2 binds to DNA region containing promoter activity

ppe2 gene is predicted to be present in operon with non pe/ppe genes, cobq1 and cobu as ppe2-cobq1-cobu. Thus, we were interested to investigate the role of ppe2 in operonic organization. We performed microscale thermophoresis (MST) experiment which revealed that PPE2 protein could bind to upstream DNA segments of ppe2-cobq1-cobu operon. Upstream region of ppe2 had shown promoter activity in β-gal assay. In this study, for the first time, a physical interaction between PPE2 protein and DNA fragment was reported, suggesting that PPE2 protein plays a role in the regulation of the putative ppe2-cobq1-cobu operon, via unknown mechanism.

The Roles of Host Noncoding RNAs in Mycobacterium tuberculosis Infection

Tuberculosis remains a major health problem. Mycobacterium tuberculosis, the causative agent of tuberculosis, can replicate and persist in host cells. Noncoding RNAs (ncRNAs) widely participate in various biological processes, including Mycobacterium tuberculosis infection, and play critical roles in gene regulation. In this review, we summarize the latest reports on ncRNAs (microRNAs, piRNAs, circRNAs and lncRNAs) that regulate the host response against Mycobacterium tuberculosis infection. In the context of host-Mycobacterium tuberculosis interactions, a broad and in-depth understanding of host ncRNA regulatory mechanisms may lead to potential clinical prospects for tuberculosis diagnosis and the development of new anti-tuberculosis therapies.

Applications of Transcriptomics and Proteomics for Understanding Dormancy and Resuscitation in Mycobacterium tuberculosis

Mycobacterium tuberculosis can survive within its host for extended periods of time without any clinical symptoms of disease and reactivate when the immune system is weakened. A detailed understanding of how M. tuberculosis enters into and exits out of dormancy, is necessary in order to develop new strategies for tackling tuberculosis. Omics methodologies are unsupervised and unbiased to any hypothesis, making them useful tools for the discovery of new drug targets. This review summarizes the findings of transcriptomic and proteomic approaches toward understanding dormancy and reactivation of M. tuberculosis. Within the granuloma of latently infected individuals, the bacteria are dormant, with a marked slowdown of growth, division and metabolism. In vitro models have attempted to simulate these features by subjecting the bacterium to hypoxia, nutrient starvation, potassium depletion, growth in the presence of vitamin C, or growth in the presence of long-chain fatty acids. The striking feature of all the models is the upregulation of the DosR regulon, which includes the transcriptional regulator Rv0081, one of the central hubs of dormancy. Also upregulated are chaperone proteins, fatty acid and cholesterol degrading enzymes, the sigma factors SigE and SigB, enzymes of the glyoxylate and the methylcitrate cycle, the Clp proteases and the transcriptional regulator ClgR. Further, there is increased expression of genes involved in mycobactin synthesis, fatty acid degradation, the glyoxylate shunt and gluconeogenesis, in granulomas formed in vitro from peripheral blood mononuclear cells from latently infected individuals compared to naïve individuals. Genes linked to aerobic respiration, replication, transcription, translation and cell division, are downregulated during dormancy in vitro, but upregulated during reactivation. Resuscitation in vitro is associated with upregulation of genes linked to the synthesis of mycolic acids, phthiocerol mycocerosate (PDIM) and sulfolipids; ribosome biosynthesis, replication, transcription and translation, cell division, and genes encoding the five resuscitation promoting factors (Rpfs). The expression of proteases, transposases and insertion sequences, suggests genome reorganization during reactivation.

Mycobacterium tuberculosis Small RNA MTS1338 Confers Pathogenic Properties to Non-Pathogenic Mycobacterium smegmatis

Small non-coding RNAs play a key role in bacterial adaptation to various stresses. Mycobacterium tuberculosis small RNA MTS1338 is upregulated during mycobacteria infection of macrophages, suggesting its involvement in the interaction of the pathogen with the host. In this study, we explored the functional effects of MTS1338 by expressing it in non-pathogenic Mycobacterium smegmatis that lacks the MTS1338 gene. The results indicated that MTS1338 slowed the growth of the recombinant mycobacteria in culture and increased their survival in RAW 264.7 macrophages, where the MTS1338-expressing strain significantly (p < 0.05) reduced the number of mature phagolysosomes and changed the production of cytokines IL-1β, IL-6, IL-10, IL-12, TGF-β, and TNF-α compared to those of the control strain. Proteomic and secretomic profiling of recombinant and control strains revealed differential expression of proteins involved in the synthesis of main cell wall components and in the regulation of iron metabolism (ESX-3 secretion system) and response to hypoxia (furA, whiB4, phoP). These effects of MTS1338 expression are characteristic for M. tuberculosis during infection, suggesting that in pathogenic mycobacteria MTS1338 plays the role of a virulence factor supporting the residence of M. tuberculosis in the host.

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

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

Systematic evaluation, verification and comparison of tuberculosis-related non-coding RNA diagnostic panels

We systematically summarized tuberculosis (TB)‐related non‐coding RNA (ncRNA) diagnostic panels, validated and compared panel performance. We searched TB‐related ncRNA panels in PubMed, OVID and Web of Science up to 28 February 2020, and available datasets in GEO, SRA and EBI ArrayExpress up to 1 March 2020. We rebuilt models and synthesized the results of each model in validation sets by bivariate mixed models. Specificity at 90% sensitivity, area under curve (AUC) and inconsistence index (I²) were calculated. NcRNA biofunctions were analysed. Nineteen models based on 18 ncRNA panels (miRNA, lncRNA, circRNA and snoRNA panels) and 18 datasets were included. Limited available datasets only allowed to evaluate miRNA panels further. Cui 2017 and Latorre 2015 exhibited specificity >70% at 90% sensitivity and AUC >80% in all validation sets. Cui 2017 showed higher specificity at 90% sensitivity (92%) and AUC (95%) and lower heterogeneity (I² = 0%) in ethological‐confirmation validation sets. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis indicated that most ncRNAs in panels involved in immune cell activation, oxidative stress, and Wnt and MAPK signalling pathway. Cui 2017 outperformed other models in both all available and aetiological‐confirmed validation sets, meeting the criteria of target product profile of WHO. This work provided a basis for clinical choice of TB‐related ncRNA diagnostic panels to a certain extent.

Development and Optimization of Chromosomally-Integrated Fluorescent Mycobacterium tuberculosis Reporter Constructs

Mycobacterium tuberculosis resides in the lungs in various lesion types with unique microenvironmental conditions. This diversity is in line with heterogeneous disease progression and divergent drug efficiency. Fluorescent reporter strains can be used to decipher the micromilieu and to guide future treatment regimens. Current reporters using replicating plasmids, however, are not suitable for long-term mouse infections or studies in non-human primates. Using a combination of recombinant DNA and protein optimization techniques, we have developed reporter strains based on integrative plasmids, which exhibit stimulus-response characteristics and fluorescence intensities comparable to those based on replicating plasmids. We successfully applied the concepts by constructing a multi-color reporter strain able to detect simultaneous changes in environmental pH, Mg²⁺ concentrations, and protein expression levels.

sncRNA-1 Is a Small Noncoding RNA Produced by Mycobacterium tuberculosis in Infected Cells That Positively Regulates Genes Coupled to Oleic Acid Biosynthesis

Nearly one third of the world’s population is infected with Mycobacterium tuberculosis (Mtb). While much work has focused on the role of different Mtb encoded proteins in pathogenesis, recent studies have revealed that Mtb also transcribes many noncoding RNAs whose functions remain poorly characterized. We performed RNA sequencing and identified a subset of Mtb H37Rv-encoded small RNAs (<30 nts in length) that were produced in infected macrophages. Designated as smaller noncoding RNAs (sncRNAs), three of these predominated the read counts. Each of the three, sncRNA-1, sncRNA-6, and sncRNA-8 had surrounding sequences with predicted stable secondary RNA stem loops. Site-directed mutagenesis of the precursor sequences suggest the existence of a hairpin loop dependent RNA processing mechanism. A functional assessment of sncRNA-1 suggested that it positively regulated two mycobacterial transcripts involved in oleic acid biosynthesis. Complementary loss- and gain- of-function approaches revealed that sncRNA-1 positively supports Mtb growth and survival in nutrient-depleted cultures as well as in infected macrophages. Overall, the findings reveal that Mtb produces sncRNAs in infected cells, with sncRNA-1 modulating mycobacterial gene expression including genes coupled to oleic acid biogenesis.

Does Mycobacterium bovis persist in cattle in a non-replicative latent state as Mycobacterium tuberculosis in human beings?

Members of the Mycobacterium tuberculosis complex (MTBC) are responsible for tuberculosis in several mammals. In this complex, Mycobacterium tuberculosis and Mycobacterium bovis, which are closely related, show host preference for humans and cattle, respectively. Although human and bovine tuberculosis are clinically similar, M. tuberculosis mostly causes latent infection in humans, whereas M. bovis frequently leads to an acute infection in cattle. This review attempts to connect the pathology in experimental animal models as well as the cellular responses to M. bovis and M. tuberculosis regarding the differences in protein expression and regulatory mechanisms of both pathogens that could explain their apparent divergent latency behaviour. The occurrence of latent bovine tuberculosis (bTB) would represent a serious complication for the eradication of the disease in cattle, with the risk of onward transmission to humans. Thus, understanding the physiological events that may lead to the state of latency in bTB could assist in the development of appropriate prevention and control tools.

MTS1338, A Small Mycobacterium tuberculosis RNA, Regulates Transcriptional Shifts Consistent With Bacterial Adaptation for Entering Into Dormancy and Survival Within Host Macrophages

Small non-coding RNAs play a significant role in bacterial adaptation to changing environmental conditions. We investigated the dynamics of expression of MTS1338, a small non-coding RNA of Mycobacterium tuberculosis, in the mouse model in vivo, regulation of its expression in the infected macrophages, and the consequences of its overexpression in bacterial cultures. Here we demonstrate that MTS1338 significantly contributes to host-pathogen interactions. Activation of the host immune system triggered NO-inducible up-regulation of MTS1338 in macrophage-engulfed mycobacteria. Constitutive overexpression of MTS1338 in cultured mycobacteria improved their survival in vitro under low pH conditions. MTS1338 up-regulation launched a spectrum of shifts in the transcriptome profile similar to those reported for M. tuberculosis adaptation to hostile intra-macrophage environment. Using the RNA-seq approach, we demonstrate that gene expression changes accompanying MTS1338 overexpression indicate reduction in translational activity and bacterial growth. These changes indicate mycobacteria entering the dormant state. Taken together, our results suggest a direct involvement of this sRNA in the interplay between mycobacteria and the host immune system during infectious process.

Mycobacterium tuberculosis 6C sRNA binds multiple mRNA targets via C-rich loops independent of RNA chaperones

Bacterial small regulatory RNAs (sRNAs) are the most abundant class of post-transcriptional regulators and have been well studied in Gram-negative bacteria. Little is known about the functions and mechanisms of sRNAs in high GC Gram-positive bacteria including Mycobacterium and Streptomyces. Here, we performed an in-depth study of 6C sRNA of Mycobacterium tuberculosis, which is conserved among high GC Gram-positive bacteria. Forty-seven genes were identified as possible direct targets of 6C sRNA and 15 of them were validated using an in vivo translational lacZ fusion system. We found that 6C sRNA plays a pleotropic role and regulates genes involved in various cellular processes, including DNA replication and protein secretion. Mapping the interactions of 6C sRNA with mRNA targets panD and dnaB revealed that the C-rich loops of 6C sRNA act as direct binding sites to mRNA targets. Unlike in Gram-negative bacteria where RNA binding proteins Hfq and ProQ are required, the interactions of 6C sRNA with mRNAs appear to be independent of RNA chaperones. Our findings suggest that the multiple G–C pairings between single stranded regions are sufficient to establish stable interactions between 6C sRNA and mRNA targets, providing a mechanism for sRNAs in high GC Gram-positive bacteria.

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

On a stake-out: Mycobacterial small RNA identification and regulation

Persistence of mycobacteria in the hostile environment of human macrophage is pivotal for its successful pathogenesis. Rapid adaptation to diverse stresses is the key aspect for their survival in the host cells. A range of heterogeneous mechanisms operate in bacteria to retaliate stress conditions. Small RNAs (sRNA) have been implicated in many of those mechanisms in either a single or multiple regulatory networks to post-transcriptionally modulate bacterial gene expression. Although small RNA profiling in mycobacteria by advanced technologies like deep sequencing, tilling microarray etc. have identified hundreds of sRNA, however, a handful of those small RNAs have been unearthed with precise regulatory mechanism. Extensive investigations on sRNA-mediated gene regulations in eubacteria like Escherichia coli revealed the existence of a plethora of distinctive sRNA mechanisms e.g. base pairing, protein sequestration, RNA decoy etc. Increasing studies on mycobacterial sRNA also discovered several eccentric mechanisms where sRNAs act at the posttranscriptional stage to either activate or repress target gene expression that lead to promote mycobacterial survival in stresses. Several intrinsic features like high GC content, absence of any homologue of abundant RNA chaperones, Hfq and ProQ, isolate sRNA mechanisms of mycobacteria from that of other bacteria. An insightful approach has been taken in this review to describe sRNA identification and its regulations in mycobacterial species especially in Mycobacterium tuberculosis.

Resuscitation of Dormant "Non-culturable" Mycobacterium tuberculosis Is Characterized by Immediate Transcriptional Burst

Under unfavorable conditions such as host immune responses and environmental stresses, human pathogen Mycobacterium tuberculosis may acquire the dormancy phenotype characterized by "non-culturability" and a substantial decrease of metabolic activity and global transcription rates. Here, we found that the transition of M. tuberculosis from the dormant "non-culturable" (NC) cells to fully replicating population in vitro occurred not earlier than 7 days after the start of the resuscitation process, with predominant resuscitation over this time interval evidenced by shortening apparent generation time up to 2.8 h at the beginning of resuscitation. The early resuscitation phase was characterized by constant, albeit low, incorporation of radioactive uracil, indicating de novo transcription immediately after the removal of the stress factor, which resulted in significant changes of the M. tuberculosis transcriptional profile already after the first 24 h of resuscitation. This early response included transcriptional upregulation of genes encoding enzymes of fatty acid synthase system type I (FASI) and type II (FASII) responsible for fatty acid/mycolic acid biosynthesis, and regulatory genes, including whiB6 encoding a redox-sensing transcription factor. The second resuscitation phase took place 4 days after the resuscitation onset, i.e., still before the start of active cell division, and included activation of central metabolism genes encoding NADH dehydrogenases, ATP-synthases, and ribosomal proteins. Our results demonstrate, for the first time, that the resuscitation of dormant NC M. tuberculosis is characterized by immediate activation of de novo transcription followed by the upregulation of genes controlling key metabolic pathways and then, cell multiplication.

Variation in pre-therapy levels of selected Mycobacterium tuberculosis transcripts in sputum and their relationship with 2-month culture conversion

Background:The abundance of transcripts arising fromMycobacterium tuberculosis(MTB) in sputum pre-chemotherapy may enhance our understanding of factors influencing treatment response. We hypothesized that differences in the prevalence of pre-existing slowly metabolizing MTB in sputum may be partially responsible for differences in the rate of sputum clearance during treatment.Methods:Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to characterize a selected limited transcription profile of MTB in sputum pre-chemotherapy and assess inter-individual variation. The difference in cycle threshold (Ct) per gene, normalized to 16S, between exponential/stationary phase culture and sputum was calculated and stratified by 2-month culture converter status.Results:HIV-1 uninfected patients with rifampicin-susceptible tuberculosis provided sputum pre-chemotherapy; 11 patients were negative for MTB culture after two months of therapy and 8 remained culture-positive.Increasedicl1andprpDandrpsN2:rpsN1in sputum relative to culture suggested cholesterol utilization and a low-zinc environment respectively. IncreasedhspXand decreasedatpAandnuoGrelative to exponential culture suggested a slowly metabolizing subpopulation of MTB. While the thehspXhiatpAlonuoGlosignal varied, we did not observe statistically significant enrichment of this phenotype in the non-converter population nor an association with MTB-lineage.Conclusion:Differential abundance of selected informative transcripts suggested a metabolically less-active subpopulation with a prevalence that varied between individual untreated patients.

Hypoxia Is Not a Main Stress When Mycobacterium tuberculosis Is in a Dormancy-Like Long-Chain Fatty Acid Environment

The capacity of Mycobacterium tuberculosis (Mtb) to sense, respond and adapt to a variable and hostile environment within the host makes it one of the most successful human pathogens. During different stages of infection, Mtb is surrounded by a plethora of lipid molecules and current evidence points out the relevance of fatty acids during the infectious process. In this study, we have compared the transcriptional response of Mtb to hypoxia in cultures supplemented with a mix of even long-chain fatty acids or dextrose as main carbon sources. Using RNA sequencing, we have identified differential expressed genes in early and late hypoxia, defined according to the in vitro Wayne and Hayes model, and compared the results with the exponential phase of growth in both carbon sources. We show that the number of genes over-expressed in the lipid medium was quite low in both, early and late hypoxia, relative to conditions including dextrose, with the exception of transcripts of stable and non-coding RNAs, which were more expressed in the fatty acid medium. We found that sigB and sigE were over-expressed in the early phase of hypoxia, confirming their pivotal role in early adaptation to low oxygen concentration independently of the carbon source. A drastic contrast was found with the transcriptional regulatory factors at early hypoxia. Only 2 transcriptional factors were over-expressed in early hypoxia in the lipid medium compared to 37 that were over-expressed in the dextrose medium. Instead of Rv0081, known to be the central regulator of hypoxia in dextrose, Rv2745c (ClgR), seems to play a main role in hypoxia in the fatty acid medium. The low level of genes associated to the stress-response during their adaptation to hypoxia in fatty acids, suggests that this lipid environment makes hypoxia a less stressful condition for the tubercle bacilli. Taken all together, these results indicate that the presence of lipid molecules shapes the metabolic response of Mtb to an adaptive state for different stresses within the host, including hypoxia. This fact could explain the success of Mtb to establish long-term survival during latent infection.

Ms1 RNA increases the amount of RNA polymerase in Mycobacterium smegmatis

Ms1 is a sRNA recently found in mycobacteria and several other actinobacterial species. Ms1 interacts with the RNA polymerase (RNAP) core devoid of sigma factors, which differs from 6S RNA that binds to RNAP holoenzymes containing the primary sigma factor. Here we show that Ms1 is the most abundant non-rRNA transcript in stationary phase in Mycobacterium smegmatis. The accumulation of Ms1 stems from its high-level synthesis combined with decreased degradation. We identify the Ms1 promoter, P_Ms1 , and cis-acting elements important for its activity. Furthermore, we demonstrate that PNPase (an RNase) contributes to the differential accumulation of Ms1 during growth. Then, by comparing the transcriptomes of wt and ΔMs1 strains from stationary phase, we reveal that Ms1 affects the intracellular level of RNAP. The absence of Ms1 results in decreased levels of the mRNAs encoding β and β' subunits of RNAP, which is also reflected at the protein level. Thus, the ΔMs1 strain has a smaller pool of RNAPs available when the transcriptional demand increases. This contributes to the inability of the ΔMs1 strain to rapidly react to environmental changes during outgrowth from stationary phase.

Recent Trends in System-Scale Integrative Approaches for Discovering Protective Antigens Against Mycobacterial Pathogens

Mycobacterial infections are one of the deadliest infectious diseases still posing a major health burden worldwide. The battle against these pathogens needs to focus on novel approaches and key interventions. In recent times, availability of genome scale data has revolutionized the fields of computational biology and immunoproteomics. Here, we summarize the cutting-edge ‘omics’ technologies and innovative system scale strategies exploited to mine the available data. These may be targeted using high-throughput technologies to expedite the identification of novel antigenic candidates for the rational next generation vaccines and serodiagnostic development against mycobacterial pathogens for which traditional methods have been failing.

Regulatory RNA in Mycobacterium tuberculosis, back to basics

Since the turn of the millenium, RNA-based control of gene expression has added an extra dimension to the central dogma of molecular biology. Still, the roles of Mycobacterium tuberculosis regulatory RNAs and the proteins that facilitate their functions remain elusive, although there can be no doubt that RNA biology plays a central role in the baterium's adaptation to its many host environments. In this review, we have presented examples from model organisms and from M. tuberculosis to showcase the abundance and versatility of regulatory RNA, in order to emphasise the importance of these 'fine-tuners' of gene expression.

Small RNAs in Bacterial Virulence and Communication

Bacterial pathogens must endure or adapt to different environments and stresses during transmission and infection. Posttranscriptional gene expression control by regulatory RNAs, such as small RNAs and riboswitches, is now considered central to adaptation in many bacteria, including pathogens. The study of RNA-based regulation (riboregulation) in pathogenic species has provided novel insight into how these bacteria regulate virulence gene expression. It has also uncovered diverse mechanisms by which bacterial small RNAs, in general, globally control gene expression. Riboregulators as well as their targets may also prove to be alternative targets or provide new strategies for antimicrobials. In this article, we present an overview of the general mechanisms that bacteria use to regulate with RNA, focusing on examples from pathogens. In addition, we also briefly review how deep sequencing approaches have aided in opening new perspectives in small RNA identification and the study of their functions. Finally, we discuss examples of riboregulators in two model pathogens that control virulence factor expression or survival-associated phenotypes, such as stress tolerance, biofilm formation, or cell-cell communication, to illustrate how riboregulation factors into regulatory networks in bacterial pathogens.

Small regulatory bacterial RNAs regulating the envelope stress response

Most bacteria encode a large repertoire of RNA-based regulatory mechanisms. Recent discoveries have revealed that the expression of many genes is controlled by a plethora of base-pairing noncoding small regulatory RNAs (sRNAs), regulatory RNA-binding proteins and RNA-degrading enzymes. Some of these RNA-based regulated processes respond to stress conditions and are involved in the maintenance of cellular homeostasis. They achieve it by either direct posttranscriptional repression of several mRNAs, including blocking access to ribosome and/or directing them to RNA degradation when the synthesis of their cognate proteins is unwanted, or by enhanced translation of some key stress-regulated transcriptional factors. Noncoding RNAs that regulate the gene expression by binding to regulatory proteins/transcriptional factors often act negatively by sequestration, preventing target recognition. Expression of many sRNAs is positively regulated by stress-responsive sigma factors like RpoE and RpoS, and two-component systems like PhoP/Q, Cpx and Rcs. Some of these regulatory RNAs act via a feedback mechanism on their own regulators, which is best reflected by recent discoveries, concerning the regulation of cell membrane composition by sRNAs in Escherichia coli and Salmonella, which are highlighted here.

Expression, maturation and turnover of DrrS, an unusually stable, DosR regulated small RNA in Mycobacterium tuberculosis

Mycobacterium tuberculosis depends on the ability to adjust to stresses encountered in a range of host environments, adjustments that require significant changes in gene expression. Small RNAs (sRNAs) play an important role as post-transcriptional regulators of prokaryotic gene expression, where they are associated with stress responses and, in the case of pathogens, adaptation to the host environment. In spite of this, the understanding of M. tuberculosis RNA biology remains limited. Here we have used a DosR-associated sRNA as an example to investigate multiple aspects of mycobacterial RNA biology that are likely to apply to other M. tuberculosis sRNAs and mRNAs. We have found that accumulation of this particular sRNA is slow but robust as cells enter stationary phase. Using reporter gene assays, we find that the sRNA core promoter is activated by DosR, and we have renamed the sRNA DrrS for DosR Regulated sRNA. Moreover, we show that DrrS is transcribed as a longer precursor, DrrS+, which is rapidly processed to the mature and highly stable DrrS. We characterise, for the first time in mycobacteria, an RNA structural determinant involved in this extraordinary stability and we show how the addition of a few nucleotides can lead to acute destabilisation. Finally, we show how this RNA element can enhance expression of a heterologous gene. Thus, the element, as well as its destabilising derivatives may be employed to post-transcriptionally regulate gene expression in mycobacteria in combination with different promoter variants. Moreover, our findings will facilitate further investigations into the severely understudied topic of mycobacterial RNA biology and into the role that regulatory RNA plays in M. tuberculosis pathogenesis.

Comprehensive Essentiality Analysis of the Mycobacterium tuberculosis Genome via Saturating Transposon Mutagenesis

For decades, identifying the regions of a bacterial chromosome that are necessary for viability has relied on mapping integration sites in libraries of random transposon mutants to find loci that are unable to sustain insertion. To date, these studies have analyzed subsaturated libraries, necessitating the application of statistical methods to estimate the likelihood that a gap in transposon coverage is the result of biological selection and not the stochasticity of insertion. As a result, the essentiality of many genomic features, particularly small ones, could not be reliably assessed. We sought to overcome this limitation by creating a completely saturated transposon library in Mycobacterium tuberculosis In assessing the composition of this highly saturated library by deep sequencing, we discovered that a previously unknown sequence bias of the Himar1 element rendered approximately 9% of potential TA dinucleotide insertion sites less permissible for insertion. We used a hidden Markov model of essentiality that accounted for this unanticipated bias, allowing us to confidently evaluate the essentiality of features that contained as few as 2 TA sites, including open reading frames (ORF), experimentally identified noncoding RNAs, methylation sites, and promoters. In addition, several essential regions that did not correspond to known features were identified, suggesting uncharacterized functions that are necessary for growth. This work provides an authoritative catalog of essential regions of the M. tuberculosis genome and a statistical framework for applying saturating mutagenesis to other bacteria.

IMPORTANCE

Sequencing of transposon-insertion mutant libraries has become a widely used tool for probing the functions of genes under various conditions. The Himar1 transposon is generally believed to insert with equal probabilities at all TA dinucleotides, and therefore its absence in a mutant library is taken to indicate biological selection against the corresponding mutant. Through sequencing of a saturated Himar1 library, we found evidence that TA dinucleotides are not equally permissive for insertion. The insertion bias was observed in multiple prokaryotes and influences the statistical interpretation of transposon insertion (TnSeq) data and characterization of essential genomic regions. Using these insights, we analyzed a fully saturated TnSeq library for M. tuberculosis, enabling us to generate a comprehensive catalog of in vitro essentiality, including ORFs smaller than those found in any previous study, small (noncoding) RNAs (sRNAs), promoters, and other genomic features.

Trehalose-6-Phosphate-Mediated Toxicity Determines Essentiality of OtsB2 in Mycobacterium tuberculosis In Vitro and in Mice

Trehalose biosynthesis is considered an attractive target for the development of antimicrobials against fungal, helminthic and bacterial pathogens including Mycobacterium tuberculosis. The most common biosynthetic route involves trehalose-6-phosphate (T6P) synthase OtsA and T6P phosphatase OtsB that generate trehalose from ADP/UDP-glucose and glucose-6-phosphate. In order to assess the drug target potential of T6P phosphatase, we generated a conditional mutant of M. tuberculosis allowing the regulated gene silencing of the T6P phosphatase gene otsB2. We found that otsB2 is essential for growth of M. tuberculosis in vitro as well as for the acute infection phase in mice following aerosol infection. By contrast, otsB2 is not essential for the chronic infection phase in mice, highlighting the substantial remodelling of trehalose metabolism during infection by M. tuberculosis. Blocking OtsB2 resulted in the accumulation of its substrate T6P, which appears to be toxic, leading to the self-poisoning of cells. Accordingly, blocking T6P production in a ΔotsA mutant abrogated otsB2 essentiality. T6P accumulation elicited a global upregulation of more than 800 genes, which might result from an increase in RNA stability implied by the enhanced neutralization of toxins exhibiting ribonuclease activity. Surprisingly, overlap with the stress response caused by the accumulation of another toxic sugar phosphate molecule, maltose-1-phosphate, was minimal. A genome-wide screen for synthetic lethal interactions with otsA identified numerous genes, revealing additional potential drug targets synergistic with OtsB2 suitable for combination therapies that would minimize the emergence of resistance to OtsB2 inhibitors.

Trehalose biosynthesis is considered an attractive target for the development of new drugs against various microbial pathogens including Mycobacterium tuberculosis. In this human pathogen, two partially redundant pathways mediate trehalose biosynthesis. The OtsA-OtsB2 pathway, which dominates in culture, involves trehalose-6-phosphate (T6P) synthase OtsA and T6P phosphatase OtsB2. While OtsA is dispensable, OtsB2 is strictly essential for growth of M. tuberculosis. Using conditional gene silencing, we here show that essentiality of OtsB2 is linked to accumulation of its substrate T6P, which exhibits direct or indirect toxic effects. Regulated gene expression in vivo revealed that OtsB2 is required to establish an acute infection of M. tuberculosis in a mouse infection model, but is surprisingly fully dispensable during the chronic infection phase. This highlights that trehalose metabolism of M. tuberculosis is substantially remodelled during infection.

Role of Small Noncoding RNAs in Bacterial Metabolism

The study of prokaryotic small RNAs is one of the most important directions in modern molecular biology. In the last decade, multiple short regulatory transcripts have been found in prokaryotes, and for some of them functional roles have been elucidated. Bacterial small RNAs are implicated in the regulation of transcription and translation, and they affect mRNA stability and gene expression via different mechanisms, including changes in mRNA conformation and interaction with proteins. Most small RNAs are expressed in response to external factors, and they help bacteria to adapt to changing environmental conditions. Bacterial infections of various origins remain a serious medical problem, despite significant progress in fighting them. Discovery of mechanisms that bacteria employ to survive in infected organisms and ways to block these mechanisms is promising for finding new treatments for bacterial infections. Regulation of pathogenesis with small RNAs is an attractive example of such mechanisms. This review considers the role of bacterial small RNAs in adaptation to stress conditions. We pay special attention to the role of small RNAs in Mycobacterium tuberculosis infection, in particular during establishment and maintenance of latent infection.

The Mycobacterium tuberculosis transcriptional landscape under genotoxic stress

Background

As an intracellular human pathogen, Mycobacterium tuberculosis (Mtb) is facing multiple stressful stimuli inside the macrophage and the granuloma. Understanding Mtb responses to stress is essential to identify new virulence factors and pathways that play a role in the survival of the tubercle bacillus. The main goal of this study was to map the regulatory networks of differentially expressed (DE) transcripts in Mtb upon various forms of genotoxic stress. We exposed Mtb cells to oxidative (H₂O₂ or paraquat), nitrosative (DETA/NO), or alkylation (MNNG) stress or mitomycin C, inducing double-strand breaks in the DNA. Total RNA was isolated from treated and untreated cells and subjected to high-throughput deep sequencing. The data generated was analysed to identify DE genes encoding mRNAs, non-coding RNAs (ncRNAs), and the genes potentially targeted by ncRNAs.

Results

The most significant transcriptomic alteration with more than 700 DE genes was seen under nitrosative stress. In addition to genes that belong to the replication, recombination and repair (3R) group, mainly found under mitomycin C stress, we identified DE genes important for bacterial virulence and survival, such as genes of the type VII secretion system (T7SS) and the proline-glutamic acid/proline-proline-glutamic acid (PE/PPE) family. By predicting the structures of hypothetical proteins (HPs) encoded by DE genes, we found that some of these HPs might be involved in mycobacterial genome maintenance. We also applied a state-of-the-art method to predict potential target genes of the identified ncRNAs and found that some of these could regulate several genes that might be directly involved in the response to genotoxic stress.

Conclusions

Our study reflects the complexity of the response of Mtb in handling genotoxic stress. In addition to genes involved in genome maintenance, other potential key players, such as the members of the T7SS and PE/PPE gene family, were identified. This plethora of responses is detected not only at the level of DE genes encoding mRNAs but also at the level of ncRNAs and their potential targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3132-1) contains supplementary material, which is available to authorized users.

Long non-coding RNA UCA1 promotes the tumorigenesis in pancreatic cancer

The contribution of long non-coding RNAs (lncRNAs) to tumorigenesis and metastasis of pancreatic cancer (PC) remains largely unknown. Urothelial cancer-associated 1 (UCA1), which is an originally identified lncRNA in bladder cancer, has be proved to play a pivotal role in bladder cancer progression and embryonic development. In this study, we detected the mRNA expression of UCA1 in 128 PC patients by qRT-PCR, and found that UCA1 expression was significantly, up-regulated in tumor tissues than that in matched adjacent non-tumor tissues (p<0.05). Clinicopathological analysis demonstrated that UCA1 expression in PC significantly correlated with malignant potential factors such as tumor size (p=0.021), depth of invasion (p=0.033), CA19-9 level (p=0.034) and tumor stage (p=0.013). Cox proportional hazards regression analysis also confirmed that high UCA1 expression was an independent prognostic biomarker of PC (p=0.046), which led to an obviously shorter 5-year overall survival (OS) compared to those patients with low UCA1 expressions (p=0.018). Furthermore, we effectively down-regulated UCA1 mRNA expression by transfecting RNA interfere fragments into SW-1990 cells, and our results in vitro indicated that down-regulation of UCA1 could effectively inhibit the cell proliferative activities, induce apoptotic rate and cause cell cycle arrest in PC cells (p<0.05). Meanwhile, UCA1 expression negative-correlated with p27 in PC tissues (r²=0.46, p<0.01), and knockdown of p27 partly abrogated the cell proliferative activities caused by UCA1 (p<0.05). Our results raised the possibility of using UCA1 as a potential prognostic biomarker and therapy target of PC, and down-regulation of UCA1 might be considered to be a novel molecular treatment strategy for patients with PC.

RNA-Seq analysis uncovers non-coding small RNA system of Mycobacterium neoaurum in the metabolism of sterols to accumulate steroid intermediates

BACKGROUND

Understanding the metabolic mechanism of sterols to produce valuable steroid intermediates in mycobacterium by a noncoding small RNA (sRNA) view is still limited. In the work, RNA-seq was implemented to investigate the noncoding transcriptome of Mycobacterium neoaurum (Mn) in the transformation process of sterols to valuable steroid intermediates, including 9α-hydroxy-4-androstene-3,17-dione (9OHAD), 1,4-androstadiene-3,17-dione (ADD), and 22-hydroxy-23, 24-bisnorchola-1,4-dien-3-one (1,4-BNA).

RESULTS

A total of 263 sRNA candidates were predicted from the intergenic regions in Mn. Differential expression of sRNA candidates was explored in the wide type Mn with vs without sterol addition, and the steroid intermediate producing Mn strains vs wide type Mn with sterol addition, respectively. Generally, sRNA candidates were differentially expressed in various strains, but there were still some shared candidates with outstandingly upregulated or downregulated expression in these steroid producing strains. Accordingly, four regulatory networks were constructed to reveal the direct and/or indirect interactions between sRNA candidates and their target genes in four groups, including wide type Mn with vs without sterol addition, 9OHAD, ADD, and BNA producing strains vs wide type Mn with sterol addition, respectively. Based on these constructed networks, several highly focused sRNA candidates were discovered to be prevalent in the networks, which showed comprehensive regulatory roles in various cellular processes, including lipid transport and metabolism, amino acid transport and metabolism, signal transduction, cell envelope biosynthesis and ATP synthesis. To explore the functional role of sRNA candidates in Mn cells, we manipulated the overexpression of candidates 131 and 138 in strain Mn-9OHAD, which led to enhanced production of 9OHAD from 1.5- to 2.3-fold during 6 d' fermentation and a slight effect on growth rate.

CONCLUSIONS

This study revealed the complex and important regulatory roles of noncoding small RNAs in the metabolism of sterols to produce steroid intermediates in Mn, further analysis of which will promote the better understanding about the molecular metabolism of these sRNA candidates and open a broad range of opportunities in the field.

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b

Background

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models.

Results

The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses.

Conclusion

The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2528-2) contains supplementary material, which is available to authorized users.

Discovery and profiling of small RNAs responsive to stress conditions in the plant pathogen Pectobacterium atrosepticum

Background

Small RNAs (sRNAs) have emerged as important regulatory molecules and have been studied in several bacteria. However, to date, there have been no whole-transcriptome studies on sRNAs in any of the Soft Rot Enterobacteriaceae (SRE) group of pathogens. Although the main ecological niches for these pathogens are plants, a significant part of their life cycle is undertaken outside their host within adverse soil environment. However, the mechanisms of SRE adaptation to this harsh nutrient-deficient environment are poorly understood.

Results

In the study reported herein, by using strand-specific RNA-seq analysis and in silico sRNA predictions, we describe the sRNA pool of Pectobacterium atrosepticum and reveal numerous sRNA candidates, including those that are induced during starvation-activated stress responses. Consequently, strand-specific RNA-seq enabled detection of 137 sRNAs and sRNA candidates under starvation conditions; 25 of these sRNAs were predicted for this bacterium in silico. Functional annotations were computationally assigned to 68 sRNAs. The expression of sRNAs in P. atrosepticum was compared under growth-promoting and starvation conditions: 68 sRNAs were differentially expressed with 47 sRNAs up-regulated under nutrient-deficient conditions. Conservation analysis using BLAST showed that most of the identified sRNAs are conserved within the SRE. Subsequently, we identified 9 novel sRNAs within the P. atrosepticum genome.

Conclusions

Since many of the identified sRNAs are starvation-induced, the results of our study suggests that sRNAs play key roles in bacterial adaptive response. Finally, this work provides a basis for future experimental characterization and validation of sRNAs in plant pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2376-0) contains supplementary material, which is available to authorized users.

Dormant non-culturable Mycobacterium tuberculosis retains stable low-abundant mRNA

Background

Dormant Mycobacterium tuberculosis bacilli are believed to play an important role in latent tuberculosis infection. Previously, we have demonstrated that cultivation of M. tuberculosis in K⁺-deficient medium resulted in generation of dormant cells. These bacilli were non-culturable on solid media (a key feature of dormant M. tuberculosis in vivo) and characterized by low metabolism and tolerance to anti-tuberculosis drugs. The dormant bacteria demonstrated a high potential to reactivation after K⁺ reintroduction even after prolonged persistence under rifampicin. In this work, we studied the transcriptome and stability of transcripts in persisting dormant bacilli under arrest of mRNA de novo synthesis.

Results

RNA-seq-based analysis of the dormant non-culturable population obtained under rifampicin exposure revealed a 30–50-fold decrease of the total mRNA level, indicating global transcriptional repression. However, the analysis of persisting transcripts displayed a cohort of mRNA molecules coding for biosynthetic enzymes, proteins involved in adaptation and repair processes, detoxification, and control of transcription initiation. This ‘dormant transcriptome’ demonstrated considerable stability during M. tuberculosis persistence and mRNA de novo synthesis arrest. On the contrary, several small non-coding RNAs showed increased abundance on dormancy. Interestingly, M. tuberculosis entry into dormancy was accompanied by the cleavage of 23S ribosomal RNA at a specific point located outside the ribosome catalytic center.

Conclusions

Dormant non-culturable M. tuberculosis bacilli are characterized by a global transcriptional repression. At the same time, the dormant bacilli retain low-abundant mRNAs, which are considerably stable during in vitro persistence, reflecting their readiness for translation upon early resuscitation steps. Increased abundance of non-coding RNAs on dormancy may indicate their role in the entry into and maintenance of M. tuberculosis dormant non-culturable state.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2197-6) contains supplementary material, which is available to authorized users.

Dissemination of 6S RNA among bacteria

6S RNA is a highly abundant small non-coding RNA widely spread among diverse bacterial groups. By competing with DNA promoters for binding to RNA polymerase (RNAP), the RNA regulates transcription on a global scale. RNAP produces small product RNAs derived from 6S RNA as template, which rearranges the 6S RNA structure leading to dissociation of 6S RNA:RNAP complexes. Although 6S RNA has been experimentally analysed in detail for some species, such as Escherichia coli and Bacillus subtilis, and was computationally predicted in many diverse bacteria, a complete and up-to-date overview of the distribution among all bacteria is missing. In this study we searched with new methods for 6S RNA genes in all currently available bacterial genomes. We ended up with a set of 1,750 6S RNA genes, of which 1,367 are novel and bona fide, distributed among 1,610 bacteria, and had a few tentative candidates among the remaining 510 assembled bacterial genomes accessible. We were able to confirm two tentative candidates by Northern blot analysis. We extended 6S RNA genes of the Flavobacteriia significantly in length compared to the present Rfam entry. We describe multiple homologs of 6S RNAs (including split 6S RNA genes) and performed a detailed synteny analysis.