Relaxed selection drives a noisy noncoding transcriptome in members of the Mycobacterium tuberculosis complex

WhiB-like proteins: Diversity of structure, function and mechanism

The WhiB-Like (Wbl) proteins are a large family of iron-sulfur (Fe-S) cluster-containing transcription factors exclusively found in the phylum Actinobacteria, including the notable genera like Mycobacteria, Streptomycetes and Corynebacteria. These proteins play pivotal roles in diverse biological processes, such as cell development, redox stress response and antibiotic resistance. Members of the Wbl family exhibit remarkable diversity in their sequences, structures and functions, attracting great attention since their first discovery. This review highlights the most recent breakthroughs in understanding the structural and mechanistic aspects of Wbl-dependent transcriptional regulation.

Discovery of a novel type IIb RelBE toxin-antitoxin system in Mycobacterium tuberculosis defined by co-regulation with an antisense RNA

Toxin‐antitoxin loci regulate adaptive responses to stresses associated with the host environment and drug exposure. Phylogenomic studies have shown that Mycobacterium tuberculosis encodes a naturally expanded type II toxin‐antitoxin system, including ParDE/RelBE superfamily members. Type II toxins are presumably regulated exclusively through protein–protein interactions with type II antitoxins. However, experimental observations in M. tuberculosis indicated that additional control mechanisms regulate RelBE2 type II loci under host‐associated stress conditions. Herein, we describe for the first time a novel antisense RNA, termed asRelE2, that co‐regulates RelE2 production via targeted processing by the Mtb RNase III, Rnc. We find that convergent expression of this coding‐antisense hybrid TA locus, relBE2‐asrelE2, is controlled in a cAMP‐dependent manner by the essential cAMP receptor protein transcription factor, Crp, in response to the host‐associated stresses of low pH and nutrient limitation. Ex vivo survival studies with relE2 and asrelE2 knockout strains showed that RelE2 contributes to Mtb survival in activated macrophages and low pH to nutrient limitation. To our knowledge, this is the first report of a novel tripartite type IIb TA loci and antisense post‐transcriptional regulation of a type II TA loci.

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.

uORF-mediated riboregulation controls transcription of whiB7/wblC antibiotic resistance gene

WhiB7/WblC is a transcriptional factor of actinomycetes conferring intrinsic resistance to multiple translation-inhibitory antibiotics. It positively autoregulates its own transcription in response to the same antibiotics. The presence of a uORF and a potential Rho-independent transcription terminator in the 5' leader region has suggested a possibility that the whiB7/wblC gene is regulated via a uORF-mediated transcription attenuation. However, experimental evidence for the molecular mechanism to explain how antibiotic stress suppresses the attenuator, if any, and induces transcription of the whiB7/wblC gene has been lacking. Here we report that the 5' leader sequences of the whiB7/wblC genes in sub-clades of actinomycetes include conserved antiterminator RNA structures. We confirmed that the putative antiterminator in the whiB7/wblC leader sequences of both Streptomyces and Mycobacterium indeed suppresses Rho-independent transcription terminator and facilitates transcription readthrough, which is required for WhiB7/WblC-mediated antibiotic resistance. The antibiotic-mediated suppression of the attenuator can be recapitulated by amino acid starvation, indicating that translational inhibition of uORF by multiple signals is a key to induce whiB7/wblC expression. Our findings of a mechanism leading to intrinsic antibiotic resistance could provide an alternative to treat drug-resistant mycobacteria.

Genome-scale analyses of transcriptional start sites in Mycobacterium marinum under normoxic and hypoxic conditions

Background

Hypoxic stress plays a critical role in the persistence of Mycobacterium tuberculosis (Mtb) infection, but the mechanisms underlying this adaptive response remain ill defined.

Material and methods

In this study, using M. marinum as a surrogate, we analyzed hypoxic responses at the transcriptional level by Cappable-seq and regular RNA-seq analyses.

Results

A total of 6808 transcriptional start sites (TSSs) were identified under normoxic and hypoxic conditions. Among these TSSs, 1112 were upregulated and 1265 were downregulated in response to hypoxic stress. Using SigE-recognized consensus sequence, we identified 59 SigE-dependent promoters and all were upregulated under hypoxic stress, suggesting an important role for SigE in this process. We also compared the performance of Cappable-seq and regular RNA-seq using the same RNA samples collected from normoxic and hypoxic conditions, and confirmed that Cappable-seq is a valuable approach for global transcriptional regulation analyses.

Conclusions

Our results provide insights and information for further characterization of responses to hypoxia in mycobacteria, and prove that Cappable-seq is a valuable approach for global transcriptional studies in mycobacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07572-8.

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.

Discovery of ANTAR-RNAs and their Mechanism of Action in Mycobacteria

Non-coding RNAs play pivotal roles in bacterial signaling. However, RNAs from certain phyla (specially high-GC actinobacteria) still remain elusive. Here, by re-engineering the existing genome-wide search approach, we discover a family of structurally conserved RNAs that are present ubiquitously across actinobacteria, including mycobacteria. In vitro analysis shows that RNAs belonging to this family bind response-regulator proteins that contain the widely prevalent ANTAR domain. The Mycobacterium tuberculosis ANTAR protein gets phosphorylated by a histidine kinase and interacts with RNA only in its phosphorylated state. These newly identified RNAs reside only in certain transcripts and typically overlap with the ribosome-binding site, regulating translation of these transcripts. In this way, the RNAs directly link signaling pathways to translational control, thus expanding the mechanistic tool kit available for ANTAR-based control of gene expression. In mycobacteria, we find that RNAs targeted by ANTAR proteins majorly encode enzymes of lipid metabolism and associated redox pathways. This now allows us to identify the key genes that mediate ANTAR-dependent control of lipid metabolism. Our study establishes the identity and wide prevalence of ANTAR-target RNAs in mycobacteria, bringing RNA-mediated regulation in these bacteria to the center stage.

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.

Defining the Transcriptional and Post-transcriptional Landscapes of Mycobacterium smegmatis in Aerobic Growth and Hypoxia

The ability of Mycobacterium tuberculosis to infect, proliferate, and survive during long periods in the human lungs largely depends on the rigorous control of gene expression. Transcriptome-wide analyses are key to understanding gene regulation on a global scale. Here, we combine 5′-end-directed libraries with RNAseq expression libraries to gain insight into the transcriptome organization and post-transcriptional mRNA cleavage landscape in mycobacteria during log phase growth and under hypoxia, a physiologically relevant stress condition. Using the model organism Mycobacterium smegmatis, we identified 6,090 transcription start sites (TSSs) with high confidence during log phase growth, of which 67% were categorized as primary TSSs for annotated genes, and the remaining were classified as internal, antisense, or orphan, according to their genomic context. Interestingly, over 25% of the RNA transcripts lack a leader sequence, and of the coding sequences that do have leaders, 53% lack a strong consensus Shine-Dalgarno site. This indicates that like M. tuberculosis, M. smegmatis can initiate translation through multiple mechanisms. Our approach also allowed us to identify over 3,000 RNA cleavage sites, which occur at a novel sequence motif. To our knowledge, this represents the first report of a transcriptome-wide RNA cleavage site map in mycobacteria. The cleavage sites show a positional bias toward mRNA regulatory regions, highlighting the importance of post-transcriptional regulation in gene expression. We show that in low oxygen, a condition associated with the host environment during infection, mycobacteria change their transcriptomic profiles and endonucleolytic RNA cleavage is markedly reduced, suggesting a mechanistic explanation for previous reports of increased mRNA half-lives in response to stress. In addition, a number of TSSs were triggered in hypoxia, 56 of which contain the binding motif for the sigma factor SigF in their promoter regions. This suggests that SigF makes direct contributions to transcriptomic remodeling in hypoxia-challenged mycobacteria. Taken together, our data provide a foundation for further study of both transcriptional and posttranscriptional regulation in mycobacteria.

New insights into the transposition mechanisms of IS6110 and its dynamic distribution between Mycobacterium tuberculosis Complex lineages

The insertion Sequence IS6110, only present in the pathogens of the Mycobacterium tuberculosis Complex (MTBC), has been the gold-standard epidemiological marker for TB for more than 25 years, but biological implications of IS6110 transposition during MTBC adaptation to humans remain elusive. By studying 2,236 clinical isolates typed by IS6110-RFLP and covering the MTBC, we remarked a lineage-specific content of IS6110 being higher in modern globally distributed strains. Once observed the IS6110 distribution in the MTBC, we selected representative isolates and found a correlation between the normalized expression of IS6110 and its abundance in MTBC chromosomes. We also studied the molecular regulation of IS6110 transposition and we found a synergistic action of two post-transcriptional mechanisms: a -1 ribosomal frameshift and a RNA pseudoknot which interferes translation. The construction of a transcriptionally active transposase resulted in 20-fold increase of the transposition frequency. Finally, we examined transposition in M. bovis and M. tuberculosis during laboratory starvation and in a mouse infection model of TB. Our results shown a higher transposition in M. tuberculosis, that preferably happens during TB infection in mice and after one year of laboratory culture, suggesting that IS6110 transposition is dynamically adapted to the host and to adverse growth conditions.

Since the pioneering discovery of transposition by Barbara McClintock in eukaryotes and later in prokaryotes by Robert W. Hedges and Alan E. Jacob, it has become clear the key role of mobile genetics elements in chromosome remodelling, microbial evolution and host adaptation. The insertion sequence IS6110 is widely recognized for its utility in TB diagnosis and epidemiology because it is only present in the M. tuberculosis Complex (MTBC) and its transposition provides an excellent chromosomal polymorphic variability allowing the study of recent TB transmission. This inherent feature of IS6110 leads us to hypothesize that IS6110 plays a crucial role during the TB infectious cycle. However, the biological significance of IS6110 has been hindered by its almost exclusive use as an epidemiological marker. Here, we study the regulatory mechanisms and the distribution of IS6110 in the different MTBC lineages. We discuss the potential biological implications of IS6110, that is much more than an excellent TB epidemiological tool. Since IS6110 could play an important role in the adaptation of MTBC to the host, this study opens new avenues to decipher the biological roles of IS6110 in TB pathogenesis.

Gene expression models based on a reference laboratory strain are poor predictors of Mycobacterium tuberculosis complex transcriptional diversity

Every year, species of the Mycobacterium tuberculosis complex (MTBC) kill more people than any other infectious disease caused by a single agent. As a consequence of its global distribution and parallel evolution with the human host the bacteria is not genetically homogeneous. The observed genetic heterogeneity has relevance at different phenotypic levels, from gene expression to epidemiological dynamics. However, current systems biology datasets have focused on the laboratory reference strain H37Rv. By using large expression datasets testing the role of almost two hundred transcription factors, we have constructed computational models to grab the expression dynamics of Mycobacterium tuberculosis H37Rv genes. However, we have found that many of those transcription factors are deleted or likely dysfunctional across strains of the MTBC. As a result, we failed to predict expression changes in strains with a different genetic background when compared with experimental data. These results highlight the importance of designing systems biology approaches that take into account the genetic diversity of tubercle bacilli, or any other pathogen, if we want to identify universal targets for vaccines, diagnostics and treatments.

Functional characterization of the Mycobacterium abscessus genome coupled with condition specific transcriptomics reveals conserved molecular strategies for host adaptation and persistence

Background

Mycobacterium abscessus subsp. abscessus (MAB) is a highly drug resistant mycobacterium and the most common respiratory pathogen among the rapidly growing non-tuberculous mycobacteria. MAB is also one of the most deadly of the emerging cystic fibrosis (CF) pathogens requiring prolonged treatment with multiple antibiotics. In addition to its “mycobacterial” virulence genes, the genome of MAB harbours a large accessory genome, presumably acquired via lateral gene transfer including homologs shared with the CF pathogens Pseudomonas aeruginosa and Burkholderia cepacia. While multiple genome sequences are available there is little functional genomics data available for this important pathogen.

Results

We report here the first multi-omics approach to characterize the primary transcriptome, coding potential and potential regulatory regions of the MAB genome utilizing differential RNA sequencing (dRNA-seq), RNA-seq, Ribosome profiling and LC-MS proteomics. In addition we attempt to address the genomes contribution to the molecular systems that underlie MAB’s adaptation and persistence in the human host through an examination of MABs transcriptional response to a number of clinically relevant conditions. These include hypoxia, exposure to sub-inhibitory concentrations of antibiotics and growth in an artificial sputum designed to mimic the conditions within the cystic fibrosis lung.

Conclusions

Our integrated map provides the first comprehensive view of the primary transcriptome of MAB and evidence for the translation of over one hundred new short open reading frames (sORFs). Our map will act as a resource for ongoing functional genomics characterization of MAB and our transcriptome data from clinically relevant stresses informs our understanding of MAB’s adaptation to life in the CF lung. MAB’s adaptation to growth in artificial CF sputum highlights shared metabolic strategies with other CF pathogens including P. aeruginosa and mirrors the transcriptional responses that lead to persistence in mycobacteria. These strategies include an increased reliance on amino acid metabolism, and fatty acid catabolism and highlights the relevance of the glyoxylate shunt to growth in the CF lung. Our data suggests that, similar to what is seen in chronically persisting P. aeruginosa, progression towards a biofilm mode of growth would play a more prominent role in a longer-term MAB infection. Finally, MAB’s transcriptional response to antibiotics highlights the role of antibiotic modifications enzymes, active transport and the evolutionarily conserved WhiB7 driven antibiotic resistance regulon.

Electronic supplementary material

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

Induction of a stable sigma factor SigR by translation-inhibiting antibiotics confers resistance to antibiotics

Antibiotic-producing streptomycetes are rich sources of resistance mechanisms against endogenous and exogenous antibiotics. An ECF sigma factor σ^R (SigR) is known to govern the thiol-oxidative stress response in Streptomyces coelicolor. Amplification of this response is achieved by producing an unstable isoform of σ^R called σ^R′. In this work, we present evidence that antibiotics induce the SigR regulon via a redox-independent pathway, leading to antibiotic resistance. The translation-inhibiting antibiotics enhanced the synthesis of stable σ^R, eliciting a prolonged response. WblC/WhiB7, a WhiB-like DNA-binding protein, is responsible for inducing sigRp1 transcripts encoding the stable σ^R. The amount of WblC protein and its binding to the sigRp1 promoter in vivo increased upon antibiotic treatment. A similar phenomenon appears to exist in Mycobacterium tuberculosis as well. These findings reveal a novel antibiotic-induced resistance mechanism conserved among actinomycetes, and also give an explicit example of overlap in cellular damage and defense mechanisms between thiol-oxidative and anti- translational stresses.