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

MmpL3, Wag31, and PlrA are involved in coordinating polar growth with peptidoglycan metabolism and nutrient availability

Cell growth in mycobacteria involves cell wall expansion that is restricted to the cell poles. The DivIVA homolog Wag31 is required for this process, but the molecular mechanism and protein partners of Wag31 have not been described. In this study of Mycobacterium smegmatis, we identify a connection between wag31 and trehalose monomycolate (TMM) transporter mmpl3 in a suppressor screen and show that Wag31 and polar regulator PlrA are required for MmpL3's polar localization. In addition, the localization of PlrA and MmpL3 is responsive to nutrient and energy deprivation and inhibition of peptidoglycan metabolism. We show that inhibition of MmpL3 causes delocalized cell wall metabolism but does not delocalize MmpL3 itself. We found that cells with an MmpL3 C-terminal truncation, which is defective for localization, have only minor defects in polar growth but are impaired in their ability to downregulate cell wall metabolism under stress. Our work suggests that, in addition to its established function in TMM transport, MmpL3 has a second function in regulating global cell wall metabolism in response to stress. Our data are consistent with a model in which the presence of TMMs in the periplasm stimulates polar elongation and in which the connection between Wag31, PlrA, and the C-terminus of MmpL3 is involved in detecting and responding to stress in order to coordinate the synthesis of the different layers of the mycobacterial cell wall in changing conditions.

IMPORTANCE

This study is performed in Mycobacterium smegmatis, which is used as a model to understand the basic physiology of pathogenic mycobacteria such as Mycobacterium tuberculosis. In this work, we examine the function and regulation of three proteins involved in regulating cell wall elongation in mycobacterial cells, which occurs at the cell tips or poles. We find that Wag31, a regulator of polar elongation, works partly through the regulation of MmpL3, a transporter of cell wall constituents and an important drug target. Our work suggests that, beyond its transport function, MmpL3 has another function in controlling cell wall synthesis broadly in response to stress.

Diverse intrinsic properties shape transcript stability and stabilization in Mycolicibacterium smegmatis

Mycobacteria regulate transcript degradation to facilitate adaptation to environmental stress. However, the mechanisms underlying this regulation are unknown. Here we sought to gain understanding of the mechanisms controlling mRNA stability by investigating the transcript properties associated with variance in transcript stability and stress-induced transcript stabilization. We measured mRNA half-lives transcriptome-wide in Mycolicibacterium smegmatis in log phase growth and hypoxia-induced growth arrest. The transcriptome was globally stabilized in response to hypoxia, but transcripts of essential genes were generally stabilized more than those of non-essential genes. We then developed machine learning models that enabled us to identify the non-linear collective effect of a compendium of transcript properties on transcript stability and stabilization. We identified properties that were more predictive of half-life in log phase as well as properties that were more predictive in hypoxia, and many of these varied between leadered and leaderless transcripts. In summary, we found that transcript properties are differentially associated with transcript stability depending on both the transcript type and the growth condition. Our results reveal the complex interplay between transcript features and microenvironment that shapes transcript stability in mycobacteria.

MmpL3, Wag31 and PlrA are involved in coordinating polar growth with peptidoglycan metabolism and nutrient availability

Cell growth in mycobacteria involves cell wall expansion that is restricted to the cell poles. The DivIVA homolog Wag31 is required for this process, but the molecular mechanism and protein partners of Wag31 have not been described. In this study of Mycobacterium smegmatis, we identify a connection between wag31 and trehalose monomycolate (TMM) transporter mmpl3 in a suppressor screen, and show that Wag31 and polar regulator PlrA are required for MmpL3's polar localization. In addition, the localization of PlrA and MmpL3 are responsive to nutrient and energy deprivation and inhibition of peptidoglycan metabolism. We show that inhibition of MmpL3 causes delocalized cell wall metabolism, but does not delocalize MmpL3 itself. We found that cells with an MmpL3 C-terminal truncation, which is defective for localization, have only minor defects in polar growth, but are impaired in their ability to downregulate cell wall metabolism under stress. Our work suggests that, in addition to its established function in TMM transport, MmpL3 has a second function in regulating global cell wall metabolism in response to stress. Our data are consistent with a model in which the presence of TMMs in the periplasm stimulates polar elongation, and in which the connection between Wag31, PlrA and the C-terminus of MmpL3 is involved in detecting and responding to stress in order to coordinate synthesis of the different layers of the mycobacterial cell wall in changing conditions.

Prominent transcriptomic changes in Mycobacterium intracellulare under acidic and oxidative stress

BACKGROUND

Mycobacterium avium complex (MAC), including Mycobacterium intracellulare is a member of slow-growing mycobacteria and contributes to a substantial proportion of nontuberculous mycobacterial lung disease in humans affecting immunocompromised and elderly populations. Adaptation of pathogens in hostile environments is crucial in establishing infection and persistence within the host. However, the sophisticated cellular and molecular mechanisms of stress response in M. intracellulare still need to be fully explored. We aimed to elucidate the transcriptional response of M. intracellulare under acidic and oxidative stress conditions.

RESULTS

At the transcriptome level, 80 genes were shown [FC] ≥ 2.0 and p < 0.05 under oxidative stress with 10 mM hydrogen peroxide. Specifically, 77 genes were upregulated, while 3 genes were downregulated. In functional analysis, oxidative stress conditions activate DNA replication, nucleotide excision repair, mismatch repair, homologous recombination, and tuberculosis pathways. Additionally, our results demonstrate that DNA replication and repair system genes, such as dnaB, dinG, urvB, uvrD2, and recA, are indispensable for resistance to oxidative stress. On the contrary, 878 genes were shown [FC] ≥ 2.0 and p < 0.05 under acidic stress with pH 4.5. Among these genes, 339 were upregulated, while 539 were downregulated. Functional analysis highlighted nitrogen and sulfur metabolism pathways as the primary responses to acidic stress. Our findings provide evidence of the critical role played by nitrogen and sulfur metabolism genes in the response to acidic stress, including narGHIJ, nirBD, narU, narK3, cysND, cysC, cysH, ferredoxin 1 and 2, and formate dehydrogenase.

CONCLUSION

Our results suggest the activation of several pathways potentially critical for the survival of M. intracellulare under a hostile microenvironment within the host. This study indicates the importance of stress responses in M. intracellulare infection and identifies promising therapeutic targets.

The RsfSR two-component system regulates SigF function by monitoring the state of the respiratory electron transport chain in Mycobacterium smegmatis

In Mycobacterium smegmatis, the transcriptional activity of the alternative sigma factor SigF is posttranslationally regulated by the partner switching system consisting of SigF, the anti-SigF RsbW1, and three anti-SigF antagonists (RsfA, RsfB, and RsbW3). We previously demonstrated that expression of the SigF regulon is strongly induced in the Δaa3 mutant of M. smegmatis lacking the aa3 cytochrome c oxidase, the major terminal oxidase in the respiratory electron transport chain. Here, we identified and characterized the RsfSR two-component system involved in regulating the phosphorylation state of the major anti-SigF antagonist RsfB. RsfS (MSMEG_6130) is a histidine kinase with the cyclase/histidine kinase-associated sensing extracellular 3 domain at its N terminus, and RsfR (MSMEG_6131) is a receiver domain-containing protein phosphatase 2C-type phosphatase that can dephosphorylate phosphorylated RsfB. We demonstrated that phosphorylation of RsfR on Asp74 by RsfS reduces the phosphatase activity of RsfR toward phosphorylated RsfB and that the cellular abundance of the active unphosphorylated RsfB is increased in the Δaa3 mutant relative to the WT strain. We also demonstrated that the RsfSR two-component system is required for induction of the SigF regulon under respiration-inhibitory conditions such as inactivation of the cytochrome bcc1 complex and aa3 cytochrome c oxidase, as well as hypoxia, electron donor-limiting, high ionic strength, and low pH conditions. Collectively, our results reveal a key regulatory element involved in regulating the SigF signaling system by monitoring the state of the respiratory electron transport chain.

What the Hel: recent advances in understanding rifampicin resistance in bacteria

Rifampicin is a clinically important antibiotic that binds to, and blocks the DNA/RNA channel of bacterial RNA polymerase (RNAP). Stalled, nonfunctional RNAPs can be removed from DNA by HelD proteins; this is important for maintenance of genome integrity. Recently, it was reported that HelD proteins from high G+C Actinobacteria, called HelR, are able to dissociate rifampicin-stalled RNAPs from DNA and provide rifampicin resistance. This is achieved by the ability of HelR proteins to dissociate rifampicin from RNAP. The HelR-mediated mechanism of rifampicin resistance is discussed here, and the roles of HelD/HelR in the transcriptional cycle are outlined. Moreover, the possibility that the structurally similar HelD proteins from low G+C Firmicutes may be also involved in rifampicin resistance is explored. Finally, the discovery of the involvement of HelR in rifampicin resistance provides a blueprint for analogous studies to reveal novel mechanisms of bacterial antibiotic resistance.

Mycobacterial RNase E cleaves with a distinct sequence preference and controls the degradation rates of most Mycolicibacterium smegmatis mRNAs

The mechanisms and regulation of RNA degradation in mycobacteria have been subject to increased interest following the identification of interplay between RNA metabolism and drug resistance. Mycobacteria encode multiple ribonucleases predicted to participate in mRNA degradation and/or processing of stable RNAs. RNase E is hypothesized to play a major role in mRNA degradation because of its essentiality in mycobacteria and its role in mRNA degradation in gram-negative bacteria. Here, we defined the impact of RNase E on mRNA degradation rates transcriptome-wide in the nonpathogenic model Mycolicibacterium smegmatis. RNase E played a rate-limiting role in degradation of the transcripts encoded by at least 89% of protein-coding genes, with leadered transcripts often being more affected by RNase E repression than leaderless transcripts. There was an apparent global slowing of transcription in response to knockdown of RNase E, suggesting that M. smegmatis regulates transcription in responses to changes in mRNA degradation. This compensation was incomplete, as the abundance of most transcripts increased upon RNase E knockdown. We assessed the sequence preferences for cleavage by RNase E transcriptome-wide in M. smegmatis and Mycobacterium tuberculosis and found a consistent bias for cleavage in C-rich regions. Purified RNase E had a clear preference for cleavage immediately upstream of cytidines, distinct from the sequence preferences of RNase E in gram-negative bacteria. We furthermore report a high-resolution map of mRNA cleavage sites in M. tuberculosis, which occur primarily within the RNase E-preferred sequence context, confirming that RNase E has a broad impact on the M. tuberculosis transcriptome.

Gene recoding by synonymous mutations creates promiscuous intragenic transcription initiation in mycobacteria

IMPORTANCE

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, one of the deadliest infectious diseases worldwide. Previous studies have established that synonymous recoding to introduce rare codon pairings can attenuate viral pathogens. We hypothesized that non-optimal codon pairing could be an effective strategy for attenuating gene expression to create a live vaccine for Mtb. We instead discovered that these synonymous changes enabled the transcription of functional mRNA that initiated in the middle of the open reading frame and from which many smaller protein products were expressed. To our knowledge, this is one of the first reports that synonymous recoding of a gene in any organism can create or induce intragenic transcription start sites.

Transcription factor mce3R modulates antibiotics and disease persistence in Mycobacteriumtuberculosis

Transcription factors (TFs) of Mycobacterium tuberculosis (Mtb), an etiological agent of tuberculosis, regulate a network of pathways that help prolong the survival of Mtb inside the host. In this study, we have characterized a transcription repressor gene (mce3R) from the TetR family, that encodes for Mce3R protein in Mtb. We demonstrated that the mce3R gene is dispensable for the growth of Mtb on cholesterol. Gene expression analysis suggests that the transcription of genes belonging to the mce3R regulon is independent of the carbon source. We found that, in comparison to the wild type, the mce3R deleted strain (Δmce3R) generated more intracellular ROS and demonstrated reduced susceptibility to oxidative stress. Total lipid analysis suggests that mce3R regulon encoded proteins modulate the biosynthesis of cell wall lipids in Mtb. Interestingly, the absence of Mce3R increased the frequency of generation of antibiotic persisters in Mtb and imparted in-vivo growth advantage phenotype in guinea pigs. In conclusion, genes belonging to the mce3R regulon modulate the frequency of generation of persisters in Mtb. Hence, targeting mce3R regulon encoded proteins could potentiate the current regimen by eliminating persisters during Mtb infection.

K. pneumoniae and M. smegmatis infect epithelial cells via different strategies

Background

As the first line of defense, epithelial cells play a vital role in the initiation and control of both innate and adaptive immunity, which participate in the development of disease. Despite its therapeutic significance, little is understood about the specific interaction between pathogenic microorganisms and lung epithelial cells.

Methods

In this study, we performed a head-to-head comparison of the virulence and infection mechanisms of Klebsiella pneumoniae (K. pneumoniae) and Mycobacterium smegmatis (M. smegmatis), which represent Gram-negative/positive respiratory pathogens, respectively, in lung epithelial cell models for the first time.

Results

Through scanning electron microscopy combined with bacterial infection experiments, we confirmed the ability of K. pneumoniae and M. smegmatis strains to form biofilm and cord factor out of the cell wall. M. smegmatis has stronger adhesion and intracellular retention ability, while K. pneumoniae is more likely to induce acute infection. These pathogens could stay and proliferate in lung epithelial cells and stimulate the secretion of specific cytokines and chemokines through a gene transcription regulator. M. smegmatis infection can promote crosstalk among epithelial cells and other immune cells in the lung from a very early stage by prompting the secretion of pro-inflammatory cytokines. Meanwhile, there were significant correlations between K. pneumonia infection and higher levels of interleukin-15 (IL-15), interleukin-1Rα (IL-1Rα), fibroblast growth factor (FGF) basic, and granulocyte colony-stimulating factor (G-CSF). At the same time, K. pneumonia infection also led to changes in the expression of cytoskeletal proteins in epithelial cells.

Conclusions

Our results emphasized the immunoprotection and immunomodulation of lung epithelial cells against exogenous pathogenic microorganisms, indicating that different pathogens damaged the host through different strategies and induced varying innate immune responses. At the same time, they provided important clues and key immune factors for dealing with complicated pulmonary infections.

The small non-coding RNA B11 regulates multiple facets of Mycobacterium abscessus virulence

Mycobacterium abscessus causes severe disease in patients with cystic fibrosis. Little is known in M. abscessus about the roles of small regulatory RNAs (sRNA) in gene regulation. We show that the sRNA B11 controls gene expression and virulence-associated phenotypes in this pathogen. B11 deletion from the smooth strain ATCC_19977 produced a rough strain, increased pro-inflammatory signaling and virulence in multiple infection models, and increased resistance to antibiotics. Examination of clinical isolate cohorts identified isolates with B11 mutations or reduced expression. We used RNAseq and proteomics to investigate the effects of B11 on gene expression and test the impact of mutations found in clinical isolates. Over 200 genes were differentially expressed in the deletion mutant. Strains with the clinical B11 mutations showed expression trends similar to the deletion mutant, suggesting partial loss of function. Among genes upregulated in the B11 mutant, there was a strong enrichment for genes with B11-complementary sequences in their predicted ribosome binding sites (RBS), consistent with B11 functioning as a negative regulator that represses translation via base-pairing to RBSs. Comparing the proteomes similarly revealed that upregulated proteins were strongly enriched for B11-complementary sequences. Intriguingly, genes upregulated in the absence of B11 included components of the ESX-4 secretion system, critical for M. abscessus virulence. Many of these genes had B11-complementary sequences at their RBSs, which we show is sufficient to mediate repression by B11 through direct binding. Altogether, our data show that B11 acts as a direct negative regulator and mediates (likely indirect) positive regulation with pleiotropic effects on gene expression and clinically important phenotypes in M. abscessus. The presence of hypomorphic B11 mutations in clinical strains is consistent with the idea that lower B11 activity may be advantageous for M. abscessus in some clinical contexts. This is the first report on an sRNA role in M. abscessus.

Mycobacterium tuberculosis Rv0494 Protein Contributes to Mycobacterial Persistence

Purpose

Fatty acid metabolism plays an important role in the survival and pathogenesis of Mycobacterium tuberculosis. During dormancy, lipids are considered to be the main source of energy. A previous study found that Rv0494 is a starvation-inducible, lipid-responsive transcriptional regulator. However, the role of Rv0494 in bacterial persister survival has not been studied.

Methods

We constructed a Rv0494 deletion mutant strain of Mycobacterium tuberculosis H37Rv and evaluated the susceptibility of the mutant strain to antibiotics using a persistence test.

Results

We found that mutations in Rv0494 lead to survival defects of persisters, which reflected in increased sensitivity to isoniazid.

Conclusion

We conclude that Rv0494 is important for persister survival and may serve as a good target for developing new antibiotics that kill persister bacteria for improved treatment of persistent bacterial infections.

Transcription regulation by CarD in mycobacteria is guided by basal promoter kinetics

Bacterial pathogens like Mycobacterium tuberculosis (Mtb) employ transcription factors to adapt their physiology to the diverse environments within their host. CarD is a conserved bacterial transcription factor that is essential for viability in Mtb. Unlike classical transcription factors that recognize promoters by binding to specific DNA sequence motifs, CarD binds directly to the RNA polymerase to stabilize the open complex intermediate (RPo) during transcription initiation. We previously showed using RNA-sequencing that CarD is capable of both activating and repressing transcription in vivo. However, it is unknown how CarD achieves promoter-specific regulatory outcomes in Mtb despite binding indiscriminate of DNA sequence. We propose a model where CarD's regulatory outcome depends on the promoter's basal RPo stability and test this model using in vitro transcription from a panel of promoters with varying levels of RPo stability. We show that CarD directly activates full-length transcript production from the Mtb ribosomal RNA promoter rrnAP3 (AP3) and that the degree of transcription activation by CarD is negatively correlated with RPo stability. Using targeted mutations in the extended -10 and discriminator region of AP3, we show that CarD directly represses transcription from promoters that form relatively stable RPo. DNA supercoiling also influenced RPo stability and affected the direction of CarD regulation, indicating that the outcome of CarD activity can be regulated by factors beyond promoter sequence. Our results provide experimental evidence for how RNA polymerase-binding transcription factors like CarD can exert specific regulatory outcomes based on the kinetic properties of a promoter.

Structural and functional basis of the universal transcription factor NusG pro-pausing activity in Mycobacterium tuberculosis

Transcriptional pauses mediate regulation of RNA biogenesis. DNA-encoded pause signals trigger pausing by stabilizing RNA polymerase (RNAP) swiveling and inhibiting DNA translocation. The N-terminal domain (NGN) of the only universal transcription factor, NusG/Spt5, modulates pausing through contacts to RNAP and DNA. Pro-pausing NusGs enhance pauses, whereas anti-pausing NusGs suppress pauses. Little is known about pausing and NusG in the human pathogen Mycobacterium tuberculosis (Mtb). We report that MtbNusG is pro-pausing. MtbNusG captures paused, swiveled RNAP by contacts to the RNAP protrusion and nontemplate-DNA wedged between the NGN and RNAP gate loop. In contrast, anti-pausing Escherichia coli (Eco) NGN contacts the MtbRNAP gate loop, inhibiting swiveling and pausing. Using CRISPR-mediated genetics, we show that pro-pausing NGN is required for mycobacterial fitness. Our results define an essential function of mycobacterial NusG and the structural basis of pro- versus anti-pausing NusG activity, with broad implications for the function of all NusG orthologs.

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.

Gene recoding by synonymous mutations creates promiscuous intragenic transcription initiation in mycobacteria

Each genome encodes some codons more frequently than their synonyms (codon usage bias), but codons are also arranged more frequently into specific pairs (codon pair bias). Recoding viral genomes and yeast or bacterial genes with non-optimal codon pairs has been shown to decrease gene expression. Gene expression is thus importantly regulated not only by the use of particular codons but by their proper juxtaposition. We therefore hypothesized that non-optimal codon pairing could likewise attenuate Mtb genes. We explored the role of codon pair bias by recoding Mtb genes ( rpoB, mmpL3, ndh ) and assessing their expression in the closely related and tractable model organism M. smegmatis . To our surprise, recoding caused the expression of multiple smaller protein isoforms from all three genes. We confirmed that these smaller proteins were not due to protein degradation, but instead issued from new transcription initiation sites positioned within the open reading frame. New transcripts gave rise to intragenic translation initiation sites, which in turn led to the expression of smaller proteins. We next identified the nucleotide changes associated with these new sites of transcription and translation. Our results demonstrated that apparently benign, synonymous changes can drastically alter gene expression in mycobacteria. More generally, our work expands our understanding of the codon-level parameters that control translation and transcription initiation.

IMPORTANCE

Mycobacterium tuberculosis ( Mtb ) is the causative agent of tuberculosis, one of the deadliest infectious diseases worldwide. Previous studies have established that synonymous recoding to introduce rare codon pairings can attenuate viral pathogens. We hypothesized that non-optimal codon pairing could be an effective strategy for attenuating gene expression to create a live vaccine for Mtb . We instead discovered that these synonymous changes enabled the transcription of functional mRNA that initiated in the middle of the open reading frame and from which many smaller protein products were expressed. To our knowledge, this is the first report that synonymous recoding of a gene in any organism can create or induce intragenic transcription start sites.

Mycobacterial Regulatory Systems Involved in the Regulation of Gene Expression Under Respiration-Inhibitory Conditions

Mycobacterium tuberculosis is the causative agent of tuberculosis. M. tuberculosis can survive in a dormant state within the granuloma, avoiding the host-mounting immune attack. M. tuberculosis bacilli in this state show increased tolerance to antibiotics and stress conditions, and thus the transition of M. tuberculosis to the nonreplicating dormant state acts as an obstacle to tuberculosis treatment. M. tuberculosis in the granuloma encounters hostile environments such as hypoxia, nitric oxide, reactive oxygen species, low pH, and nutrient deprivation, etc., which are expected to inhibit respiration of M. tuberculosis. To adapt to and survive in respiration-inhibitory conditions, it is required for M. tuberculosis to reprogram its metabolism and physiology. In order to get clues to the mechanism underlying the entry of M. tuberculosis to the dormant state, it is important to understand the mycobacterial regulatory systems that are involved in the regulation of gene expression in response to respiration inhibition. In this review, we briefly summarize the information regarding the regulatory systems implicated in upregulation of gene expression in mycobacteria exposed to respiration-inhibitory conditions. The regulatory systems covered in this review encompass the DosSR (DevSR) two-component system, SigF partner switching system, MprBA-SigE-SigB signaling pathway, cAMP receptor protein, and stringent response.

Artemisia afra and Artemisia annua Extracts Have Bactericidal Activity against Mycobacterium tuberculosis in Physiologically Relevant Carbon Sources and Hypoxia

Mycobacterium tuberculosis (Mtb) is a deadly pathogen and causative agent of human tuberculosis, causing ~1.5 million deaths every year. The increasing drug resistance of this pathogen necessitates novel and improved treatment strategies. A crucial aspect of the host-pathogen interaction is bacterial nutrition. In this study, Artemisia annua and Artemisia afra dichloromethane extracts were tested for bactericidal activity against Mtb strain mc26230 under hypoxia and various infection-associated carbon sources (glycerol, glucose, and cholesterol). Both extracts showed significant bactericidal activity against Mtb, regardless of carbon source. Based on killing curves, A. afra showed the most consistent bactericidal activity against Mtb for all tested carbon sources, whereas A. annua showed the highest bactericidal activity in 7H9 minimal media with glycerol. Both extracts retained their bactericidal activity against Mtb under hypoxic conditions. Further investigations are required to determine the mechanism of action of these extracts and identify their active constituent compounds.

Loss of RNase J leads to multi-drug tolerance and accumulation of highly structured mRNA fragments in Mycobacterium tuberculosis

Despite the existence of well-characterized, canonical mutations that confer high-level drug resistance to Mycobacterium tuberculosis (Mtb), there is evidence that drug resistance mechanisms are more complex than simple acquisition of such mutations. Recent studies have shown that Mtb can acquire non-canonical resistance-associated mutations that confer survival advantages in the presence of certain drugs, likely acting as stepping-stones for acquisition of high-level resistance. Rv2752c/rnj, encoding RNase J, is disproportionately mutated in drug-resistant clinical Mtb isolates. Here we show that deletion of rnj confers increased tolerance to lethal concentrations of several drugs. RNAseq revealed that RNase J affects expression of a subset of genes enriched for PE/PPE genes and stable RNAs and is key for proper 23S rRNA maturation. Gene expression differences implicated two sRNAs and ppe50-ppe51 as important contributors to the drug tolerance phenotype. In addition, we found that in the absence of RNase J, many short RNA fragments accumulate because they are degraded at slower rates. We show that the accumulated transcript fragments are targets of RNase J and are characterized by strong secondary structure and high G+C content, indicating that RNase J has a rate-limiting role in degradation of highly structured RNAs. Taken together, our results demonstrate that RNase J indirectly affects drug tolerance, as well as reveal the endogenous roles of RNase J in mycobacterial RNA metabolism.

Genome-Wide Transcription Start Sites Mapping in Methylorubrum Grown with Dichloromethane and Methanol

Dichloromethane (DCM, methylene chloride) is a toxic halogenated volatile organic compound massively used for industrial applications, and consequently often detected in the environment as a major pollutant. DCM biotransformation suggests a sustainable decontamination strategy of polluted sites. Among methylotrophic bacteria able to use DCM as a sole source of carbon and energy for growth, Methylorubrum extorquens DM4 is a longstanding reference strain. Here, the primary 5′-ends of transcripts were obtained using a differential RNA-seq (dRNA-seq) approach to provide the first transcription start site (TSS) genome-wide landscape of a methylotroph using DCM or methanol. In total, 7231 putative TSSs were annotated and classified with respect to their localization to coding sequences (CDSs). TSSs on the opposite strand of CDS (antisense TSS) account for 31% of all identified TSSs. One-third of the detected TSSs were located at a distance to the start codon inferior to 250 nt (average of 84 nt) with 7% of leaderless mRNA. Taken together, the global TSS map for bacterial growth using DCM or methanol will facilitate future studies in which transcriptional regulation is crucial, and efficient DCM removal at polluted sites is limited by regulatory processes.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

An essential periplasmic protein coordinates lipid trafficking and is required for asymmetric polar growth in mycobacteria

Mycobacteria, including the human pathogen Mycobacterium tuberculosis, grow by inserting new cell wall material at their poles. This process and that of division are asymmetric, producing a phenotypically heterogeneous population of cells that respond non-uniformly to stress (Aldridge et al., 2012; Rego et al., 2017). Surprisingly, deletion of a single gene - lamA - leads to more symmetry, and to a population of cells that is more uniformly killed by antibiotics (Rego et al., 2017). How does LamA create asymmetry? Here, using a combination of quantitative time-lapse imaging, bacterial genetics, and lipid profiling, we find that LamA recruits essential proteins involved in cell wall synthesis to one side of the cell - the old pole. One of these proteins, MSMEG_0317, here renamed PgfA, was of unknown function. We show that PgfA is a periplasmic protein that interacts with MmpL3, an essential transporter that flips mycolic acids in the form of trehalose monomycolate (TMM), across the plasma membrane. PgfA interacts with a TMM analog suggesting a direct role in TMM transport. Yet our data point to a broader function as well, as cells with altered PgfA levels have differences in the abundance of other lipids and are differentially reliant on those lipids for survival. Overexpression of PgfA, but not MmpL3, restores growth at the old poles in cells missing lamA. Together, our results suggest that PgfA is a key determinant of polar growth and cell envelope composition in mycobacteria, and that the LamA-mediated recruitment of this protein to one side of the cell is a required step in the establishment of cellular asymmetry.

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.

Promotech: a general tool for bacterial promoter recognition

Promoters are genomic regions where the transcription machinery binds to initiate the transcription of specific genes. Computational tools for identifying bacterial promoters have been around for decades. However, most of these tools were designed to recognize promoters in one or few bacterial species. Here, we present Promotech, a machine-learning-based method for promoter recognition in a wide range of bacterial species. We compare Promotech's performance with the performance of five other promoter prediction methods. Promotech outperforms these other programs in terms of area under the precision-recall curve (AUPRC) or precision at the same level of recall. Promotech is available at https://github.com/BioinformaticsLabAtMUN/PromoTech .

DNA binding and gene regulatory functions of MSMEG_2295, a repressor encoded by the dinB2 operon of Mycobacterium smegmatis

MSMEG_2295 is a TetR family protein encoded by the first gene of a Mycobacterium smegmatis (Msm) operon that expresses the gene for DinB2 (MSMEG_2294), a translesion DNA repair enzyme. We have carried out investigations to understand its function by performing DNA binding studies and gene knockout experiments. We found that the protein binds to a conserved inverted repeat sequence located upstream of the dinB2 operon and several other genes. Using a knockout of MSMEG_2295, we show that MSMEG_2295 controls the expression of at least five genes, the products of which could potentially influence carbohydrate and fatty acid metabolism as well as antibiotic and oxidative stress resistance. We have demonstrated that MSMEG_2295 is a repressor by performing complementation analysis. Knocking out of MSMEG_2295 had a significant impact on pyruvate metabolism. Pyruvate dehydrogenase activity was virtually undetectable in ΔMSMEG_2295, although in the complemented strain, it was high. We also show that knocking out of MSMEG_2295 causes resistance to H₂O₂, reversed in the complemented strain. We have further found that the mycobacterial growth inhibitor plumbagin, a compound of plant origin, acts as an inducer of MSMEG_2295 regulated genes. We, therefore, establish that MSMEG_2295 functions by exerting its role as a repressor of multiple Msm genes and that by doing so, it plays a vital role in controlling pyruvate metabolism and response to oxidative stress.

RNA Sequencing for Transcript 5'-End Mapping in Mycobacteria

Next-generation sequencing technologies facilitate the analysis of multiple important properties of transcriptomes in addition to gene expression levels. Here, we describe a method for mapping RNA 5' ends in Mycobacterium tuberculosis and Mycobacterium smegmatis, which allows the determination of transcription start sites (TSSs), comparative analysis of promoter usage under different conditions, and mapping of endoribonucleolytic cleavage sites. We describe in detail the procedures for constructing RNA sequencing libraries appropriate for RNA 5' end mapping using an Illumina sequencing platform, as well as bioinformatic procedures for data analysis.

Lsr2 and Its Novel Paralogue Mediate the Adjustment of Mycobacterium smegmatis to Unfavorable Environmental Conditions

Nucleoid-associated proteins (NAPs) are the most abundant proteins involved in bacterial chromosome organization and global transcription regulation. The mycobacterial NAP family includes many diverse proteins; some are unique to actinobacteria, and many are crucial for survival under stress (e.g., HupB and Lsr2) and/or optimal growth conditions (e.g., mycobacterial integration host factor [mIHF]).

Lsr2 is a nucleoid-associated protein (NAP) that has been found strictly in actinobacteria, including mycobacteria. It is a functional homolog of histone-like nucleoid-structuring protein (H-NS); it acts as a DNA-bridging protein that plays a role in chromosomal organization and transcriptional regulation. To date, the studies on Lsr2 have focused mainly on Mycobacterium tuberculosis. In this study, we analyze the role of Lsr2 as a transcription factor in Mycobacterium smegmatis, a saprophytic bacterium whose natural habitat (soil and water) substantially differs from those of the obligatory mycobacterial pathogens. Our chromatin immunoprecipitation-sequencing (ChIP-seq) data revealed that Lsr2 binds preferentially to AT-rich regions of the M. smegmatis chromosome. We found that Lsr2 acts mainly as a repressor, controlling gene expression either directly by binding promoter regions or indirectly through DNA loop formation and DNA coating. One of the Lsr2-repressed genes encodes polyketide synthase (MSMEG_4727), which is involved in the synthesis of lipooligosaccharides (LOSs). An M. smegmatis strain deprived of Lsr2 produces more LOSs, which is mirrored by changes in the smoothness of cells and their susceptibilities to antibiotics. Unlike M. tuberculosis, M. smegmatis additionally encodes a paralogue of Lsr2, MSMEG_1060, which is a novel member of the mycobacterial NAP family. The Lsr2 and MSMEG_1060 proteins exhibit different DNA binding specificities and chromosomal localizations. Our results suggest that these proteins help M. smegmatis cells cope with stress conditions, including hypoxia and exposure to antibiotics. Thus, the present work provides novel insight into the role of Lsr2 paralogues in the ability of a saprophytic mycobacterial species to adjust to environmental changes.

IMPORTANCE Nucleoid-associated proteins (NAPs) are the most abundant proteins involved in bacterial chromosome organization and global transcription regulation. The mycobacterial NAP family includes many diverse proteins; some are unique to actinobacteria, and many are crucial for survival under stress (e.g., HupB and Lsr2) and/or optimal growth conditions (e.g., mycobacterial integration host factor [mIHF]). Here, we present a comprehensive study concerning two functional homologues of mycobacterial H-NS: Lsr2 and its paralogue from M. smegmatis, MSMEG_1060. We found that Lsr2 plays a role in transcriptional regulation, mainly by repressing gene expression via DNA loop formation and/or DNA-coating mechanisms. Intriguingly, the number of Lsr2-mediated genes was found to increase under hypoxia. Compared to Lsr2, MSMEG_1060 exhibits a different DNA binding specificity and chromosomal localization. Since tuberculosis remains a serious worldwide health problem, studies on stress response-mediating agents, such as Lsr2, may contribute to the development of novel antituberculosis drugs.

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.

Insights into the evolutionary history of the virulent factor HBHA of Mycobacterium tuberculosis

In Mycobacterium tuberculosis, heparin-binding hemagglutinin (HBHAMT) has a relevant role in infection. It is also present in non-virulent mycobacteria and ancient actinobacteria, such as Rhodococcus opacus. To have a better understanding of the underlying mechanisms that shaped the evolutionary divergence of these proteins, we performed a comprehensive phylogenetic analysis of the regulatory sequences that drive the expression of hbha in saprophytic and pathogenic mycobacterial species. The alignment of the hbha loci showed the appearance of intergenic sequences containing regulatory elements upstream the hbha gene; this sequence arrangement is present only in slow-growing pathogenic mycobacteria. The heterologous expression of HBHAMT in oleaginous R. opacus PD630 results in protein binding to lipid droplets, as it happens with HBHA proteins from saprophytic mycobacteria. We hypothesize that mycobacterial hbha gene cluster underwent functional divergence during the evolutionary differentiation of slow-growing pathogenic mycobacteria. We propose here an evolutionary scenario to explain the structural and functional divergence of HBHA in fast and slow-growing mycobacteria.

A snapshot of translation in Mycobacterium tuberculosis during exponential growth and nutrient starvation revealed by ribosome profiling

Mycobacterium tuberculosis, which causes tuberculosis, can undergo prolonged periods of non-replicating persistence in the host. The mechanisms underlying this are not fully understood, but translational regulation is thought to play a role. A large proportion of mRNA transcripts expressed in M. tuberculosis lack canonical bacterial translation initiation signals, but little is known about the implications of this for fine-tuning of translation. Here, we perform ribosome profiling to characterize the translational landscape of M. tuberculosis under conditions of exponential growth and nutrient starvation. Our data reveal robust, widespread translation of non-canonical transcripts and point toward different translation initiation mechanisms compared to canonical Shine-Dalgarno transcripts. During nutrient starvation, patterns of ribosome recruitment vary, suggesting that regulation of translation in this pathogen is more complex than originally thought. Our data represent a rich resource for others seeking to understand translational regulation in bacterial pathogens.

Re-sensitization of Mycobacterium smegmatis to Rifampicin Using CRISPR Interference Demonstrates Its Utility for the Study of Non-essential Drug Resistance Traits

A greater understanding of the genes involved in antibiotic resistance in Mycobacterium tuberculosis (Mtb) is necessary for the design of improved therapies. Clustered regularly interspaced short palindromic repeat interference (CRISPRi) has been previously utilized in mycobacteria to identify novel drug targets by the demonstration of gene essentiality. The work presented here shows that it can also be usefully applied to the study of non-essential genes involved in antibiotic resistance. The expression of an ADP-ribosyltransferase (Arr) involved in rifampicin resistance in Mycobacterium smegmatis was silenced using CRISPRi and the impact on rifampicin susceptibility was measured. Gene silencing resulted in a decrease in the minimum inhibitory concentration (MIC) similar to that previously reported in an arr deletion mutant. There is contradictory evidence for the toxicity of Streptococcus pyogenes dCas9 (dCas9Spy) in the literature. In this study the expression of dCas9Spy in M. smegmatis showed no impact on viability. Silencing was achieved with concentrations of the aTc inducer lower than previously described and with shorter induction times. Finally, designing small guide RNAs (sgRNAs) that target transcription initiation, or the early stages of elongation had the most impact on rifampicin susceptibility. This study demonstrates that CRISPRi based gene silencing can be as impactful as gene deletion for the study of non-essential genes and further contributes to the knowledge on the design and induction of sgRNAs for CRISPRi. This approach can be applied to other non-essential antimicrobial resistance genes such as drug efflux pumps.

Induction of the cydAB Operon Encoding the bd Quinol Oxidase Under Respiration-Inhibitory Conditions by the Major cAMP Receptor Protein MSMEG_6189 in Mycobacterium smegmatis

The respiratory electron transport chain (ETC) of Mycobacterium smegmatis is terminated with two terminal oxidases, the aa ₃ cytochrome c oxidase and the cytochrome bd quinol oxidase. The bd quinol oxidase with a higher binding affinity for O₂ than the aa ₃ oxidase is known to play an important role in aerobic respiration under oxygen-limiting conditions. Using relevant crp1 (MSMEG_6189) and crp2 (MSMEG_0539) mutant strains of M. smegmatis, we demonstrated that Crp1 plays a predominant role in induction of the cydAB operon under ETC-inhibitory conditions. Two Crp-binding sequences were identified upstream of the cydA gene, both of which are necessary for induction of cydAB expression under ETC-inhibitory conditions. The intracellular level of cAMP in M. smegmatis was found to be increased under ETC-inhibitory conditions. The crp2 gene was found to be negatively regulated by Crp1 and Crp2, which appears to lead to significantly low cellular abundance of Crp2 relative to Crp1 in M. smegmatis. Our RNA sequencing analyses suggest that in addition to the SigF partner switching system, Crp1 is involved in induction of gene expression in M. smegmatis exposed to ETC-inhibitory conditions.

The diversity of Shine-Dalgarno sequences sheds light on the evolution of translation initiation

Shine-Dalgarno (SD) sequences, the core element of prokaryotic ribosome-binding sites, facilitate mRNA translation by base-pair interaction with the anti-SD (aSD) sequence of 16S rRNA. In contrast to this paradigm, an inspection of thousands of prokaryotic species unravels tremendous SD sequence diversity both within and between genomes, whereas aSD sequences remain largely static. The pattern has led many to suggest unidentified mechanisms for translation initiation. Here we review known translation-initiation pathways in prokaryotes. Moreover, we seek to understand the cause and consequence of SD diversity through surveying recent advances in biochemistry, genomics, and high-throughput genetics. These findings collectively show: (1) SD:aSD base pairing is beneficial but nonessential to translation initiation. (2) The 5' untranslated region of mRNA evolves dynamically and correlates with organismal phylogeny and ecological niches. (3) Ribosomes have evolved distinct usage of translation-initiation pathways in different species. We propose a model portraying the SD diversity shaped by optimization of gene expression, adaptation to environments and growth demands, and the species-specific prerequisite of ribosomes to initiate translation. The model highlights the coevolution of ribosomes and mRNA features, leading to functional customization of the translation apparatus in each organism.

The Partner Switching System of the SigF Sigma Factor in Mycobacterium smegmatis and Induction of the SigF Regulon Under Respiration-Inhibitory Conditions

The partner switching system (PSS) of the SigF regulatory pathway in Mycobacterium smegmatis has been previously demonstrated to include the anti-sigma factor RsbW (MSMEG_1803) and two anti-sigma factor antagonists RsfA and RsfB. In this study, we further characterized two additional RsbW homologs and revealed the distinct roles of three RsbW homologs [RsbW1 (MSMEG_1803), RsbW2 (MSMEG_6129), and RsbW3 (MSMEG_1787)] in the SigF PSS. RsbW1 and RsbW2 serve as the anti-sigma factor of SigF and the protein kinase phosphorylating RsfB, respectively, while RsbW3 functions as an anti-SigF antagonist through its protein interaction with RsbW1. Using relevant mutant strains, RsfB was demonstrated to be the major anti-SigF antagonist in M. smegmatis. The phosphorylation state of Ser-63 was shown to determine the functionality of RsfB as an anti-SigF antagonist. RsbW2 was demonstrated to be the only protein kinase that phosphorylates RsfB in M. smegmatis. Phosphorylation of Ser-63 inactivates RsfB to render it unable to interact with RsbW1. Our comparative RNA sequencing analysis of the wild-type strain of M. smegmatis and its isogenic Δaa₃ mutant strain lacking the aa₃ cytochrome c oxidase of the respiratory electron transport chain revealed that expression of the SigF regulon is strongly induced under respiration-inhibitory conditions in an RsfB-dependent way.

Arrayed CRISPRi and quantitative imaging describe the morphotypic landscape of essential mycobacterial genes

Mycobacterium tuberculosis possesses a large number of genes of unknown or predicted function, undermining fundamental understanding of pathogenicity and drug susceptibility. To address this challenge, we developed a high-throughput functional genomics approach combining inducible CRISPR-interference and image-based analyses of morphological features and sub-cellular chromosomal localizations in the related non-pathogen, M. smegmatis. Applying automated imaging and analysis to 263 essential gene knockdown mutants in an arrayed library, we derive robust, quantitative descriptions of bacillary morphologies consequent on gene silencing. Leveraging statistical-learning, we demonstrate that functionally related genes cluster by morphotypic similarity and that this information can be used to inform investigations of gene function. Exploiting this observation, we infer the existence of a mycobacterial restriction-modification system, and identify filamentation as a defining mycobacterial response to histidine starvation. Our results support the application of large-scale image-based analyses for mycobacterial functional genomics, simultaneously establishing the utility of this approach for drug mechanism-of-action studies.

Low-Cost Quantitative Photothermal Genetic Detection of Pathogens on a Paper Hybrid Device Using a Thermometer

Tuberculosis (TB), one of the deadliest infectious diseases, is caused by Mycobacterium tuberculosis (MTB) and remains a public health problem nowadays. Conventional MTB DNA detection methods require sophisticated infrastructure and well-trained personnel, which leads to increasing complexity and high cost for diagnostics and limits their wide accessibility in low-resource settings. To address these issues, we have developed a low-cost photothermal biosensing method for the quantitative genetic detection of pathogens such as MTB DNA on a paper hybrid device using a thermometer. First, DNA capture probes were simply immobilized on paper through a one-step surface modification process. After DNA sandwich hybridization, oligonucleotide-functionalized gold nanoparticles (AuNPs) were introduced on paper and then catalyzed the oxidation reaction of 3,3',5,5'-tetramethylbenzidine (TMB). The produced oxidized TMB, acting as a strong photothermal agent, was used for the photothermal biosensing of MTB DNA under 808 nm laser irradiation. Under optimal conditions, the on-chip quantitative detection of the target DNA was readily achieved using an inexpensive thermometer as a signal recorder. This method does not require any expensive analytical instrumentation but can achieve higher sensitivity and there are no color interference issues, compared to conventional colorimetric methods. The method was further validated by detecting genomic DNA with high specificity. To the best of our knowledge, this is the first photothermal biosensing strategy for quantitative nucleic acid analysis on microfluidics using a thermometer, which brings fresh inspirations on the development of simple, low-cost, and miniaturized photothermal diagnostic platforms for quantitative detection of a variety of diseases at the point of care.

MamA essentiality in Mycobacterium smegmatis is explained by the presence of an apparent cognate restriction endonuclease

OBJECTIVE

Restriction-Modification (R-M) systems are ubiquitous in bacteria and were considered for years as rudimentary immune systems that protect bacterial cells from foreign DNA. Currently, these R-M systems are recognized as important players in global gene expression and other cellular processes such us virulence and evolution of genomes. Here, we report the role of the unique DNA methyltransferase in Mycobacterium smegmatis, which shows a moderate degree of sequence similarity to MamA, a previously characterized methyltransferase that affects gene expression in Mycobacterium tuberculosis and is important for survival under hypoxic conditions.

RESULTS

We found that depletion of mamA levels impairs growth and produces elongated cell bodies. Microscopy revealed irregular septation and unevenly distributed DNA, with large areas devoid of DNA and small DNA-free cells. Deletion of MSMEG_3214, a predicted endonuclease-encoding gene co-transcribed with mamA, restored the WT growth phenotype in a mamA-depleted background. Our results suggest that the mamA-depletion phenotype can be explained by DNA cleavage by the apparent cognate restriction endonuclease MSMEG_3214. In addition, in silico analysis predicts that both MamA methyltransferase and MSMEG_3214 endonuclease recognize the same palindromic DNA sequence. We propose that MamA and MSMEG_3214 constitute a previously undescribed R-M system in M. smegmatis.

Regulation of mRNA Stability During Bacterial Stress Responses

Bacteria have a remarkable ability to sense environmental changes, swiftly regulating their transcriptional and posttranscriptional machinery as a response. Under conditions that cause growth to slow or stop, bacteria typically stabilize their transcriptomes in what has been shown to be a conserved stress response. In recent years, diverse studies have elucidated many of the mechanisms underlying mRNA degradation, yet an understanding of the regulation of mRNA degradation under stress conditions remains elusive. In this review we discuss the diverse mechanisms that have been shown to affect mRNA stability in bacteria. While many of these mechanisms are transcript-specific, they provide insight into possible mechanisms of global mRNA stabilization. To that end, we have compiled information on how mRNA fate is affected by RNA secondary structures; interaction with ribosomes, RNA binding proteins, and small RNAs; RNA base modifications; the chemical nature of 5' ends; activity and concentration of RNases and other degradation proteins; mRNA and RNase localization; and the stringent response. We also provide an analysis of reported relationships between mRNA abundance and mRNA stability, and discuss the importance of stress-associated mRNA stabilization as a potential target for therapeutic development.

Polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression in mycobacteria

Genome-wide transcriptomic analyses have revealed abundant expressed short open reading frames (ORFs) in bacteria. Whether these short ORFs, or the small proteins they encode, are functional remains an open question. One quarter of mycobacterial mRNAs are leaderless, beginning with a 5'-AUG or GUG initiation codon. Leaderless mRNAs often encode unannotated short ORFs as the first gene of a polycistronic transcript. Here, we show that polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression. Detailed mutational analysis shows that one polycysteine short ORF controls expression of the downstream genes. Our data indicate that ribosomes stalled in the polycysteine tract block mRNA structures that otherwise sequester the ribosome-binding site of the 3'gene. We assessed endogenous proteomic responses to cysteine limitation in Mycobacterium smegmatis using mass spectrometry. Six cysteine metabolic loci having unannotated polycysteine-encoding leaderless short ORF architectures responded to cysteine limitation, revealing widespread cysteine-responsive attenuation in mycobacteria. Individual leaderless short ORFs confer independent operon-level control, while their shared dependence on cysteine ensures a collective response mediated by ribosome pausing. We propose the term ribulon to classify ribosome-directed regulons. Regulon-level coordination by ribosomes on sensory short ORFs illustrates one utility of the many unannotated short ORFs expressed in bacterial genomes.

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.

Vitamin B and Vitamin C Affect DNA Methylation and Amino Acid Metabolism in Mycobacterium bovis BCG

Vitamins are essential nutrients and key cofactors of enzymes that regulate cellular metabolism, and also activate the immune system. Recent studies have shown that vitamin B1 (V_{B 1}) and vitamin C (Vc) can inhibit Mycobacterium tuberculosis growth, but the precise mechanism is still not well understood. In the present study, we have used RNA-sequencing (RNA-seq), liquid chromatography coupled to mass spectrometry (LC-MS) and single-molecule real-time (SMRT) sequencing to analyze the transcriptional, metabolic and methylation profiles of Mycobacterium bovis BCG when treated with V_{B 1} and Vc. Our results show that, after vitamin treatment, variant metabolites were mainly clustered in pathways related to amino acid metabolism. Treatment with both vitamins significantly up-regulated the gene encoding cysteine synthase A. Additionally, only BCG that was treated with V_C showed m4c modifications. Genes harboring this methylation were up-regulated, suggesting that m4c methylation can promote gene transcription to some extent. Overall, this study contributes to the understanding of the effects of V_{B 1} and V_C, and suggests that these vitamins constitute potential anti-tuberculosis drugs.

The Impact of Leadered and Leaderless Gene Structures on Translation Efficiency, Transcript Stability, and Predicted Transcription Rates in Mycobacterium smegmatis

Regulation of gene expression is critical for Mycobacterium tuberculosis to tolerate stressors encountered during infection and for nonpathogenic mycobacteria such as Mycobacterium smegmatis to survive environmental stressors. Unlike better-studied models, mycobacteria express ∼14% of their genes as leaderless transcripts. However, the impacts of leaderless transcript structures on mRNA half-life and translation efficiency in mycobacteria have not been directly tested. For leadered transcripts, the contributions of 5' untranslated regions (UTRs) to mRNA half-life and translation efficiency are similarly unknown. In M. tuberculosis and M. smegmatis, the essential sigma factor, SigA, is encoded by a transcript with a relatively short half-life. We hypothesized that the long 5' UTR of sigA causes this instability. To test this, we constructed fluorescence reporters and measured protein abundance, mRNA abundance, and mRNA half-life and calculated relative transcript production rates. The sigA 5' UTR conferred an increased transcript production rate, shorter mRNA half-life, and decreased apparent translation rate compared to a synthetic 5' UTR commonly used in mycobacterial expression plasmids. Leaderless transcripts appeared to be translated with similar efficiency as those with the sigA 5' UTR but had lower predicted transcript production rates. A global comparison of M. tuberculosis mRNA and protein abundances failed to reveal systematic differences in protein/mRNA ratios for leadered and leaderless transcripts, suggesting that variability in translation efficiency is largely driven by factors other than leader status. Our data are also discussed in light of an alternative model that leads to different conclusions and suggests leaderless transcripts may indeed be translated less efficiently.IMPORTANCE Tuberculosis, caused by Mycobacterium tuberculosis, is a major public health problem killing 1.5 million people globally each year. During infection, M. tuberculosis must alter its gene expression patterns to adapt to the stress conditions it encounters. Understanding how M. tuberculosis regulates gene expression may provide clues for ways to interfere with the bacterium's survival. Gene expression encompasses transcription, mRNA degradation, and translation. Here, we used Mycobacterium smegmatis as a model organism to study how 5' untranslated regions affect these three facets of gene expression in multiple ways. We furthermore provide insight into the expression of leaderless mRNAs, which lack 5' untranslated regions and are unusually prevalent in mycobacteria.

Elucidating the functional role of Mycobacterium smegmatis recX in stress response

The RecX protein has attracted considerable interest because the recX mutants exhibit multiple phenotypes associated with RecA functions. To further our understanding of the functional relationship between recA and recX, the effect of different stress treatments on their expression profiles, cell yield and viability were investigated. A significant correlation was found between the expression of Mycobacterium smegmatis recA and recX genes at different stages of growth, and in response to different stress treatments albeit recX exhibiting lower transcript and protein abundance at the mid-log and stationary phases of the bacterial growth cycle. To ascertain their roles in vivo, a targeted deletion of the recX and recArecX was performed in M. smegmatis. The growth kinetics of these mutant strains and their sensitivity patterns to different stress treatments were assessed relative to the wild-type strain. The deletion of recA affected normal cell growth and survival, while recX deletion showed no significant effect. Interestingly, deletion of both recX and recA genes results in a phenotype that is intermediate between the phenotypes of the ΔrecA mutant and the wild-type strain. Collectively, these results reveal a previously unrecognized role for M. smegmatis recX and support the notion that it may regulate a subset of the yet unknown genes involved in normal cell growth and DNA-damage repair.

mRNA Degradation Rates Are Coupled to Metabolic Status in Mycobacterium smegmatis

The logistics of tuberculosis therapy are difficult, requiring multiple drugs for many months. Mycobacterium tuberculosis survives in part by entering nongrowing states in which it is metabolically less active and thus less susceptible to antibiotics. Basic knowledge on how M. tuberculosis survives during these low-metabolism states is incomplete, and we hypothesize that optimized energy resource management is important. Here, we report that slowed mRNA turnover is a common feature of mycobacteria under energy stress but is not dependent on the mechanisms that have generally been postulated in the literature. Finally, we found that mRNA stability and growth status can be decoupled by a drug that causes growth arrest but increases metabolic activity, indicating that mRNA stability responds to metabolic status rather than to growth rate per se. Our findings suggest a need to reorient studies of global mRNA stabilization to identify novel mechanisms that are presumably responsible.

The success of Mycobacterium tuberculosis as a human pathogen is due in part to its ability to survive stress conditions, such as hypoxia or nutrient deprivation, by entering nongrowing states. In these low-metabolism states, M. tuberculosis can tolerate antibiotics and develop genetically encoded antibiotic resistance, making its metabolic adaptation to stress crucial for survival. Numerous bacteria, including M. tuberculosis, have been shown to reduce their rates of mRNA degradation under growth limitation and stress. While the existence of this response appears to be conserved across species, the underlying bacterial mRNA stabilization mechanisms remain unknown. To better understand the biology of nongrowing mycobacteria, we sought to identify the mechanistic basis of mRNA stabilization in the nonpathogenic model Mycobacterium smegmatis. We found that mRNA half-life was responsive to energy stress, with carbon starvation and hypoxia causing global mRNA stabilization. This global stabilization was rapidly reversed when hypoxia-adapted cultures were reexposed to oxygen, even in the absence of new transcription. The stringent response and RNase levels did not explain mRNA stabilization, nor did transcript abundance. This led us to hypothesize that metabolic changes during growth cessation impact the activities of degradation proteins, increasing mRNA stability. Indeed, bedaquiline and isoniazid, two drugs with opposing effects on cellular energy status, had opposite effects on mRNA half-lives in growth-arrested cells. Taken together, our results indicate that mRNA stability in mycobacteria is not directly regulated by growth status but rather is dependent on the status of energy metabolism.