Pervasive translation in Mycobacterium tuberculosis

Integrated sequence and -omic features reveal novel small proteome of Mycobacterium tuberculosis

Bioinformatic studies on small proteins are under-represented due to difficulties in annotation posed by their small size. However, recent discoveries emphasize the functional significance of small proteins in cellular processes including cell signaling, metabolism, and adaptation to stress. In this study, we utilized a Random Forest classifier trained on sequence features, RNA-Seq, and Ribo-Seq data to uncover small proteins (smORFs) in M. tuberculosis. Independent predictions for the exponential and starvation conditions resulted in 695 potential smORFs. We examined the functional implications of these smORFs using homology searches, LC-MS/MS, and ChIP-seq data, testing their expression in diverse growth conditions, and identifying protein domains. We provide evidence that some of these smORFs could be part of operons, or exist as upstream ORFs. This expanded data resource for the proteins of M. tuberculosis would aid in fine-tuning the existing protein and gene regulatory networks, thereby improving system-wide studies. The primary goal of this study was to uncover and characterize smORFs in M. tuberculosis through bioinformatic analysis, shedding light on their functional roles and genomic organization. Further investigation of these potential smORFs would provide valuable insights into the genome organization and functional diversity of the M. tuberculosis proteome.

A novel regulatory interplay between atypical B12 riboswitches and uORF translation in Mycobacterium tuberculosis

Vitamin B12 is an essential cofactor in all domains of life and B12-sensing riboswitches are some of the most widely distributed riboswitches. Mycobacterium tuberculosis, the causative agent of tuberculosis, harbours two B12-sensing riboswitches. One controls expression of metE, encoding a B12-independent methionine synthase, the other controls expression of ppe2 of uncertain function. Here, we analysed ligand sensing, secondary structure and gene expression control of the metE and ppe2 riboswitches. Our results provide the first evidence of B12 binding by these riboswitches and show that they exhibit different preferences for individual isoforms of B12, use distinct regulatory and structural elements and act as translational OFF switches. Based on our results, we propose that the ppe2 switch represents a new variant of Class IIb B12-sensing riboswitches. Moreover, we have identified short translated open reading frames (uORFs) upstream of metE and ppe2, which modulate the expression of their downstream genes. Translation of the metE uORF suppresses MetE expression, while translation of the ppe2 uORF is essential for PPE2 expression. Our findings reveal an unexpected regulatory interplay between B12-sensing riboswitches and the translational machinery, highlighting a new level of cis-regulatory complexity in M. tuberculosis. Attention to such mechanisms will be critical in designing next-level intervention strategies.

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis (M. tb), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

Promoter capture drives the emergence of proto-genes in Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli Long-Term Evolution Experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time-span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, thereby serving as raw material for new gene emergence. Most proto-genes result either from insertion element activity or chromosomal translocations that fused pre-existing regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, persist stably, and can serve as potential substrates for new gene formation.

Exposing the small protein load of bacterial life

The ever-growing repertoire of genomic techniques continues to expand our understanding of the true diversity and richness of prokaryotic genomes. Riboproteogenomics laid the foundation for dynamic studies of previously overlooked genomic elements. Most strikingly, bacterial genomes were revealed to harbor robust repertoires of small open reading frames (sORFs) encoding a diverse and broadly expressed range of small proteins, or sORF-encoded polypeptides (SEPs). In recent years, continuous efforts led to great improvements in the annotation and characterization of such proteins, yet many challenges remain to fully comprehend the pervasive nature of small proteins and their impact on bacterial biology. In this work, we review the recent developments in the dynamic field of bacterial genome reannotation, catalog the important biological roles carried out by small proteins and identify challenges obstructing the way to full understanding of these elusive proteins.

Small proteins in Gram-positive bacteria

Small proteins comprising less than 100 amino acids have been often ignored in bacterial genome annotations. About 10 years ago, focused efforts started to investigate whole peptidomes, which resulted in the discovery of a multitude of small proteins, but only a number of them have been characterized in detail. Generally, small proteins can be either membrane or cytosolic proteins. The latter interact with larger proteins, RNA or even metal ions. Here, we summarize our current knowledge on small proteins from Gram-positive bacteria with a special emphasis on the model organism Bacillus subtilis. Our examples include membrane-bound toxins of type I toxin-antitoxin systems, proteins that block the assembly of higher order structures, regulate sporulation or modulate the RNA degradosome. We do not consider antimicrobial peptides. Furthermore, we present methods for the identification and investigation of small proteins.

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.

Alternative Reading Frames are an Underappreciated Source of Protein Sequence Novelty

Protein-coding DNA sequences can be translated into completely different amino acid sequences if the nucleotide triplets used are shifted by a non-triplet amount on the same DNA strand or by translating codons from the opposite strand. Such "alternative reading frames" of protein-coding genes are a major contributor to the evolution of novel protein products. Recent studies demonstrating this include examples across the three domains of cellular life and in viruses. These sequences increase the number of trials potentially available for the evolutionary invention of new genes and also have unusual properties which may facilitate gene origin. There is evidence that the structure of the standard genetic code contributes to the features and gene-likeness of some alternative frame sequences. These findings have important implications across diverse areas of molecular biology, including for genome annotation, structural biology, and evolutionary genomics.

In vitro activity of new combinations of β-lactam and β-lactamase inhibitors against the Mycobacterium tuberculosis complex

As meropenem-clavulanic acid is recommended for the treatment of drug-resistant tuberculosis, the repurposing of new carbapenem combinations may provide new treatment options, including oral alternatives. Therefore, we studied the in vitro activities of meropenem-vaborbactam, meropenem-clavulanic acid, and tebipenem-clavulanic acid. One hundred nine Mycobacterium tuberculosis complex (MTBC) clinical isolates were tested, of which 69 were pan-susceptible and the remaining pyrazinamide- or multidrug-resistant. Broth microdilution MICs were determined using the EUCAST reference method. Meropenem and tebipenem were tested individually and in combination with vaborbactam 8 mg/L and clavulanic-acid 2 and 4 mg/L, respectively. Whole-genome sequencing was performed to explore resistance mechanisms. Clavulanic acid lowered the modal tebipenem MIC approximately 16-fold (from 16 to 1 mg/L). The modal meropenem MIC was reduced twofold by vaborbactam compared with an approximately eightfold decrease by clavulanic acid. The only previously described high-confidence carbapenem resistance mutation, crfA T62A, was shared by a subgroup of lineage 4.3.4.1 isolates and did not correlate with elevated MICs. The presence of a β-lactamase inhibitor reduced the MTBC MICs of tebipenem and meropenem. The resulting MIC distribution was lowest for the orally available drugs tebipenem-clavulanic acid. Whether this in vitro activity translates to similar or greater clinical efficacy of tebipenem-clavulanic acid compared with the currently WHO-endorsed meropenem-clavulanic acid requires clinical studies. IMPORTANCE Repurposing of already approved antibiotics, such as β-lactams in combination with β-lactamase inhibitors, may provide new treatment alternatives for drug-resistant tuberculosis. Meropenem-clavulanic acid was more active in vitro compared to meropenem-vaborbactam. Notably, tebipenem-clavulanic acid showed even better activity, raising the potential of an all-oral treatment option. Clinical data are needed to investigate whether the better in vitro activity of tebipenem-clavulanic acid correlates with greater clinical efficacy compared with the currently WHO-endorsed meropenem-clavulanic acid.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis ( M. tb ), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Importance

Understanding mechanisms of resistance and tolerance in Mycobacterium tuberculosis ( M. tb ) can help us develop new treatments that capitalize on M. tb 's vulnerabilities. Here we used transcriptomics to study both the regulation of biofilm formation in clinical isolates as well as how those regulatory systems adapt to new environments. We find that closely related clinical populations have diverse strategies for growth under biofilm conditions, and that genetic background plays a large role in determining the trajectory of evolution. These results have implications for future treatment strategies that may be informed by our knowledge of the evolutionary constraints on strain(s) from an individual infection. This work provides new information about the mechanisms of biofilm formation in M. tb and outlines a framework for population level approaches for studying bacterial adaptation.

Between noise and function: Toward a taxonomy of the non-canonical translatome

Eukaryotic genomes are pervasively translated, but the properties of translated sequences outside of canonical genes are poorly understood. A new study in Cell Systems reveals a large translatome that is not under significant evolutionary constraint but is still an active part of diverse cellular systems.

Premature termination of transcription is shaped by Rho and translated uORFS in Mycobacterium tuberculosis

Little is known about the decisions behind transcription elongation versus termination in the human pathogen Mycobacterium tuberculosis (M.TB). By applying Term-seq to M.TB we found that the majority of transcription termination is premature and associated with translated regions, i.e., within previously annotated or newly identified open reading frames. Computational predictions and Term-seq analysis, upon depletion of termination factor Rho, suggests that Rho-dependent transcription termination dominates all transcription termination sites (TTS), including those associated with regulatory 5' leaders. Moreover, our results suggest that tightly coupled translation, in the form of overlapping stop and start codons, may suppress Rho-dependent termination. This study provides detailed insights into novel M.TB cis-regulatory elements, where Rho-dependent, conditional termination of transcription and translational coupling together play major roles in gene expression control. Our findings contribute to a deeper understanding of the fundamental regulatory mechanisms that enable M.TB adaptation to the host environment offering novel potential points of intervention.

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.

The small mycobacterial ribosomal protein, bS22, modulates aminoglycoside accessibility to its 16S rRNA helix-44 binding site

Treatment of tuberculosis continues to be challenging due to the widespread latent form of the disease and the emergence of antibiotic-resistant strains of the pathogen, Mycobacterium tuberculosis. Bacterial ribosomes are a common and effective target for antibiotics. Several second line anti-tuberculosis drugs, e.g. kanamycin, amikacin, and capreomycin, target ribosomal RNA to inhibit protein synthesis. However, M. tuberculosis can acquire resistance to these drugs, emphasizing the need to identify new drug targets. Previous cryo-EM structures of the M. tuberculosis and M. smegmatis ribosomes identified two novel ribosomal proteins, bS22 and bL37, in the vicinity of two crucial drug-binding sites: the mRNA-decoding center on the small (30S), and the peptidyl-transferase center on the large (50S) ribosomal subunits, respectively. The functional significance of these two small proteins is unknown. In this study, we observe that an M. smegmatis strain lacking the bs22 gene shows enhanced susceptibility to kanamycin compared to the wild-type strain. Cryo-EM structures of the ribosomes lacking bS22 in the presence and absence of kanamycin suggest a direct role of bS22 in modulating the 16S rRNA kanamycin-binding site. Our structures suggest that amino-acid residue Lys-16 of bS22 interacts directly with the phosphate backbone of helix 44 of 16S rRNA to influence the micro-configuration of the kanamycin-binding pocket. Our analysis shows that similar interactions occur between eukaryotic homologues of bS22, and their corresponding rRNAs, pointing to a common mechanism of aminoglycoside resistance in higher organisms.

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.

Miniprep assisted proteomics (MAP) for rapid proteomics sample preparation

Complete enzymatic digestion of proteins for bottom-up proteomics is substantially improved by use of detergents for denaturation and solubilization. Detergents however, are incompatible with many proteases and highly detrimental to LC-MS/MS. Recently; filter-based methods have seen wide use due to their capacity to remove detergents and harmful reagents prior to digestion and mass spectrometric analysis. We hypothesized that non-specific protein binding to negatively charged silica-based filters would be enhanced by addition of lyotropic salts, similar to DNA purification. We sought to exploit these interactions and investigate if low-cost DNA purification spin-filters, 'Minipreps,' efficiently and reproducibly bind proteins for digestion and LC-MS/MS analysis. We propose a new method, Miniprep Assisted Proteomics (MAP), for sample preparation. We demonstrate binding capacity, performance, recovery and identification rates for proteins and whole-cell lysates using MAP. MAP recovered equivalent or greater protein yields from 0.5-50 μg analyses benchmarked against commercial trapping preparations. Nano UHPLC-MS/MS proteome profiling of lysates of Escherichia coli had 99.3% overlap vs. existing approaches and reproducibility of replicate minipreps was 98.8% at the 1% FDR protein level. Label Free Quantitative proteomics was performed and 91.2% of quantified proteins had a %CV <20% (2044/2241). Miniprep Assisted Proteomics can be performed in minutes, shows low variability, high recovery and proteome depth. This suggests a significant role for adventitious binding in developing new proteomics sample preparation techniques. MAP represents an efficient, ultra-low-cost alternative for sample preparation in a commercially obtainable device that costs ∼$0.50 (USD) per miniprep.

Mycobacterium smegmatis: The Vanguard of Mycobacterial Research

The genus Mycobacterium contains several slow-growing human pathogens, including Mycobacterium tuberculosis, Mycobacterium leprae, and Mycobacterium avium. Mycobacterium smegmatis is a nonpathogenic and fast growing species within this genus. In 1990, a mutant of M. smegmatis, designated mc²155, that could be transformed with episomal plasmids was isolated, elevating M. smegmatis to model status as the ideal surrogate for mycobacterial research. Classical bacterial models, such as Escherichia coli, were inadequate for mycobacteria research because they have low genetic conservation, different physiology, and lack the novel envelope structure that distinguishes the Mycobacterium genus. By contrast, M. smegmatis encodes thousands of conserved mycobacterial gene orthologs and has the same cell architecture and physiology. Dissection and characterization of conserved genes, structures, and processes in genetically tractable M. smegmatis mc²155 have since provided previously unattainable insights on these same features in its slow-growing relatives. Notably, tuberculosis (TB) drugs, including the first-line drugs isoniazid and ethambutol, are active against M. smegmatis, but not against E. coli, allowing the identification of their physiological targets. Furthermore, Bedaquiline, the first new TB drug in 40 years, was discovered through an M. smegmatis screen. M. smegmatis has become a model bacterium, not only for M. tuberculosis, but for all other Mycobacterium species and related genera. With a repertoire of bioinformatic and physical resources, including the recently established Mycobacterial Systems Resource, M. smegmatis will continue to accelerate mycobacterial research and advance the field of microbiology.

Ribosome profiling enhances understanding of mycobacterial translation

A recent addition to the -omics toolkit, ribosome profiling, enables researchers to gain insight into the process and regulation of translation by mapping fragments of mRNA protected from nuclease digestion by ribosome binding. In this review, we discuss how ribosome profiling applied to mycobacteria has led to discoveries about translational regulation. Using case studies, we show that the traditional view of “canonical” translation mechanisms needs expanding to encompass features of mycobacterial translation that are more widespread than previously recognized. We also discuss the limitations of the method and potential future developments that could yield further insight into the fundamental biology of this important human pathogen.