RNIE: genome-wide prediction of bacterial intrinsic terminators

CryptKeeper: a negative design tool for reducing unintentional gene expression in bacteria

Foundational techniques in molecular biology-such as cloning genes, tagging biomolecules for purification or identification, and overexpressing recombinant proteins-rely on introducing non-native or synthetic DNA sequences into organisms. These sequences may be recognized by the transcription and translation machinery in their new context in unintended ways. The cryptic gene expression that sometimes results has been shown to produce genetic instability and mask experimental signals. Computational tools have been developed to predict individual types of gene expression elements, but it can be difficult for researchers to contextualize their collective output. Here, we introduce CryptKeeper, a software pipeline that visualizes predictions of Escherichia coli gene expression signals and estimates the translational burden possible from a DNA sequence. We investigate several published examples where cryptic gene expression in E. coli interfered with experiments. CryptKeeper accurately postdicts unwanted gene expression from both eukaryotic virus infectious clones and individual proteins that led to genetic instability. It also identifies off-target gene expression elements that resulted in truncations that confounded protein purification. Incorporating negative design using CryptKeeper into reverse genetics and synthetic biology workflows can help to mitigate cloning challenges and avoid unexplained failures and complications that arise from unintentional gene expression.

Exploring the transcription start sites and other genomic features facilitates the accurate identification and annotation of small RNAs across multiple stress conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs.

CryptKeeper: a negative design tool for reducing unintentional gene expression in bacteria

Foundational techniques in molecular biology-such as cloning genes, tagging biomolecules for purification or identification, and overexpressing recombinant proteins-rely on introducing non-native or synthetic DNA sequences into organisms. These sequences may be recognized by the transcription and translation machinery in their new context in unintended ways. The cryptic gene expression that sometimes results has been shown to produce genetic instability and mask experimental signals. Computational tools have been developed to predict individual types of gene expression elements, but it can be difficult for researchers to contextualize their collective output. Here, we introduce CryptKeeper, a software pipeline that visualizes predictions of bacterial gene expression signals and estimates the translational burden possible from a DNA sequence. We investigate several published examples where cryptic gene expression in E. coli interfered with experiments. CryptKeeper accurately postdicts unwanted gene expression from both eukaryotic virus infectious clones and individual proteins that led to genetic instability. It also identifies off-target gene expression elements that resulted in truncations that confounded protein purification. Incorporating negative design using CryptKeeper into reverse genetics and synthetic biology workflows can help to mitigate cloning challenges and avoid unexplained failures and complications that arise from unintentional gene expression.

CoCas9 is a compact nuclease from the human microbiome for efficient and precise genome editing

The expansion of the CRISPR-Cas toolbox is highly needed to accelerate the development of therapies for genetic diseases. Here, through the interrogation of a massively expanded repository of metagenome-assembled genomes, mostly from human microbiomes, we uncover a large variety (n = 17,173) of type II CRISPR-Cas loci. Among these we identify CoCas9, a strongly active and high-fidelity nuclease with reduced molecular size (1004 amino acids) isolated from an uncultivated Collinsella species. CoCas9 is efficiently co-delivered with its sgRNA through adeno associated viral (AAV) vectors, obtaining efficient in vivo editing in the mouse retina. With this study we uncover a collection of previously uncharacterized Cas9 nucleases, including CoCas9, which enriches the genome editing toolbox.

Detection of intrinsic transcription termination sites in bacteria: consensus from hairpin detection approaches

We compare the WebGeSTer and INtrinsic transcription TERmination hairPIN (INTERPIN) databases used for intrinsic transcription termination (ITT) site prediction in bacteria. The former deploys inverted nucleotide repeat detection for identification of RNA hairpin, while the latter a pair-potential function - the hairpin energy score evaluation being identical for both. We find INTERPIN more sensitive than WebGeSTer with about 6% and 51% additional predictions for ITTs in chromosomal and plasmid operons, respectively. INTERPIN hairpins are relatively shorter in length with ungapped stem, and even located in AT-rich segments, compared to GC-rich longer hairpins with a gapped stem in WebGeSTer. The GC%, length, and energy score from INTERPIN transcription units (TUs) are best inter-correlated while the lowest energy single hairpins from WebGeSTer, considered suitable for ITT, being the worst. Around 72% TUs from the two databases overlap, and ∼60% of all alternate ITT sites downstream of TUs overlap, of which 65% are cluster hairpins. This helps highlight hairpin features that can be used to identify termination sites in bacteria across different prediction methods. Overall, the pair-potential-function-based hairpins screened appear to be more consistent with the kinetic and thermodynamics processes of ITT known to date.Communicated by Ramaswamy H. Sarma.

Improved RNA stability estimation through Bayesian modeling reveals most Salmonella transcripts have subminute half-lives

RNA decay is a crucial mechanism for regulating gene expression in response to environmental stresses. In bacteria, RNA-binding proteins (RBPs) are known to be involved in posttranscriptional regulation, but their global impact on RNA half-lives has not been extensively studied. To shed light on the role of the major RBPs ProQ and CspC/E in maintaining RNA stability, we performed RNA sequencing of Salmonella enterica over a time course following treatment with the transcription initiation inhibitor rifampicin (RIF-seq) in the presence and absence of these RBPs. We developed a hierarchical Bayesian model that corrects for confounding factors in rifampicin RNA stability assays and enables us to identify differentially decaying transcripts transcriptome-wide. Our analysis revealed that the median RNA half-life in Salmonella in early stationary phase is less than 1 min, a third of previous estimates. We found that over half of the 500 most long-lived transcripts are bound by at least one major RBP, suggesting a general role for RBPs in shaping the transcriptome. Integrating differential stability estimates with cross-linking and immunoprecipitation followed by RNA sequencing (CLIP-seq) revealed that approximately 30% of transcripts with ProQ binding sites and more than 40% with CspC/E binding sites in coding or 3' untranslated regions decay differentially in the absence of the respective RBP. Analysis of differentially destabilized transcripts identified a role for ProQ in the oxidative stress response. Our findings provide insights into posttranscriptional regulation by ProQ and CspC/E, and the importance of RBPs in regulating gene expression.

Utilizing RNA-seq Data to Infer Bacterial Transcription Termination Sites and Validate Predictions

The transcription termination process is an important part of the gene expression process in the cell. It has been studied extensively, but many aspects of the mechanism are not well understood. The widespread availability of experimental RNA-seq data from high-throughput experiments provides a unique opportunity to infer the end of the transcription units genome wide. This data is available for both Rho-dependent and Rho-independent termination pathways that drive transcription termination in bacteria. Our book chapter gives an overview of the current knowledge of Rho-independent transcription termination mechanisms and the prediction approaches currently deployed to infer the termination sites. Thereafter, we describe our method that uses cluster hairpins to detect Rho-independent transcription termination sites. These clusters are a group of hairpins that lies at <15 bp from each other and are together capable of enforcing the termination process. The idea of a group of hairpins being extensively used for transcription termination is new, and results show that at least 52% of the total cases are of this type, while in the remaining cases, a single strong hairpin is capable of driving transcription termination. The reads derived from the RNA-seq data for corresponding bacteria have been used to validate the predicted sites. The predictions that match these RNA-seq derived sites have higher confidence, and we find almost 98% of the predicted sites, including alternate termination sites, to match the RNA-seq data. We discuss the features of predicted hairpins in detail for a better understanding of the Rho-independent transcription termination mechanism in bacteria. We also explain how users can use the tools developed by us to do transcription terminator predictions and design their experiments through genome-level visualization of the transcription termination sites from the precomputed INTERPIN database.

Identification of over ten thousand candidate structured RNAs in viruses and phages

Structured RNAs play crucial roles in viruses, exerting influence over both viral and host gene expression. However, the extensive diversity of structured RNAs and their ability to act in cis or trans positions pose challenges for predicting and assigning their functions. While comparative genomics approaches have successfully predicted candidate structured RNAs in microbes on a large scale, similar efforts for viruses have been lacking. In this study, we screened over 5 million DNA and RNA viral sequences, resulting in the prediction of 10,006 novel candidate structured RNAs. These predictions are widely distributed across taxonomy and ecosystem. We found transcriptional evidence for 206 of these candidate structured RNAs in the human fecal microbiome. These candidate RNAs exhibited evidence of nucleotide covariation, indicative of selective pressure maintaining the predicted secondary structures. Our analysis revealed a diverse repertoire of candidate structured RNAs, encompassing a substantial number of putative tRNAs or tRNA-like structures, Rho-independent transcription terminators, and potentially cis-regulatory structures consistently positioned upstream of genes. In summary, our findings shed light on the extensive diversity of structured RNAs in viruses, offering a valuable resource for further investigations into their functional roles and implications in viral gene expression and pave the way for a deeper understanding of the intricate interplay between viruses and their hosts at the molecular level.

INTERPIN: A repository for intrinsic transcription termination hairpins in bacteria

The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for ∼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.

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.

Incomplete transcripts dominate the Mycobacterium tuberculosis transcriptome

Mycobacterium tuberculosis (Mtb) is a bacterial pathogen that causes tuberculosis, an infectious disease that inflicts major health and economic costs around the world 1 . Mtb encounters a diversity of environments during its lifecycle, and responds to these changing environments by reprogramming its transcriptional output 2 . However, the transcriptomic features of Mtb remain poorly characterized. In this work, we comprehensively profile the Mtb transcriptome using the SEnd-seq method that simultaneously captures the 5' and 3' ends of RNA 3 . Surprisingly, we find that the RNA coverage for most of the Mtb transcription units display a gradual drop-off within a 200-500 nucleotide window downstream of the transcription start site, yielding a massive number of incomplete transcripts with heterogeneous 3' ends. We further show that the accumulation of these short RNAs is mainly due to the intrinsically low processivity of the Mtb transcription machinery rather than trans-acting factors such as Rho. Finally, we demonstrate that transcription-translation coupling plays a critical role in generating full-length protein-coding transcripts in Mtb. In sum, our results depict a mycobacterial transcriptome that is dominated by incomplete RNA products, suggesting a distinctive set of transcriptional regulatory mechanisms that could be exploited for new therapeutics.

Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions

Ample evidence supports the importance of the microbiota on human health and disease. Recent studies suggest that extracellular vesicles are an important means of bacterial-host communication, in part via the transport of small RNAs (sRNAs). Bacterial sRNAs have been shown to co-precipitate with human and mouse RNA-induced silencing complex, hinting that some may regulate gene expression as eukaryotic microRNAs do. Bioinformatic tools, including those that can incorporate an sRNA's secondary structure, can be used to predict interactions between bacterial sRNAs and human messenger RNAs (mRNAs). Validation of these potential interactions using reproducible experimental methods is essential to move the field forward. This review will cover the evidence of interspecies communication via sRNAs, bioinformatic tools currently available to identify potential bacterial sRNA-host (specifically, human) mRNA interactions, and experimental methods to identify and validate those interactions.

Automated identification of sequence-tailored Cas9 proteins using massive metagenomic data

The identification of the protospacer adjacent motif (PAM) sequences of Cas9 nucleases is crucial for their exploitation in genome editing. Here we develop a computational pipeline that was used to interrogate a massively expanded dataset of metagenome and virome assemblies for accurate and comprehensive PAM predictions. This procedure allows the identification and isolation of sequence-tailored Cas9 nucleases by using the target sequence as bait. As proof of concept, starting from the disease-causing mutation P23H in the RHO gene, we find, isolate and experimentally validate a Cas9 which uses the mutated sequence as PAM. Our PAM prediction pipeline will be instrumental to generate a Cas9 nuclease repertoire responding to any PAM requirement.

Inverse folding based pre-training for the reliable identification of intrinsic transcription terminators

It is well-established that neural networks can predict or identify structural motifs of non-coding RNAs (ncRNAs). Yet, the neural network based identification of RNA structural motifs is limited by the availability of training data that are often insufficient for learning features of specific ncRNA families or structural motifs. Aiming to reliably identify intrinsic transcription terminators in bacteria, we introduce a novel pre-training approach that uses inverse folding to generate training data for predicting or identifying a specific family or structural motif of ncRNA. We assess the ability of neural networks to identify secondary structure by systematic in silico mutagenesis experiments. In a study to identify intrinsic transcription terminators as functionally well-understood RNA structural motifs, our inverse folding based pre-training approach significantly boosts the performance of neural network topologies, which outperform previous approaches to identify intrinsic transcription terminators. Inverse-folding based pre-training provides a simple, yet highly effective way to integrate the well-established thermodynamic energy model into deep neural networks for identifying ncRNA families or motifs. The pre-training technique is broadly applicable to a range of network topologies as well as different types of ncRNA families and motifs.

Intrinsic transcriptional terminators are essential regulators in determining the 3’ end of transcripts in bacteria. The underlying mechanism involves RNA secondary structure, where nucleotides fold into a specific hairpin motif. Identifying terminator sequences in bacterial genomes has conventionally been approached with well-established energy models for structural motifs. However, the folding mechanism of transcription terminators is understood only partially, limiting the success of energy-model based identification.

Neural networks have been proposed to overcome these limitations. However, their adoption for predicting and identifying RNA secondary structure has been a double edged sword: Neural networks promise to learn features that are not represented by the energy models, while they are black boxes that lack explicit modeling assumptions and may fail to account for features that are well understandable based on decades-old energy models.

Here, we introduce a pre-training approach for neural networks that uses energy-model based inverse folding of structural motifs. As we demonstrate, this approach “brings back the energy model” to identify transcriptional terminators and overcomes the limitations of previous energy-model based predictions. Our approach works for diverse types of neural networks, and is suitable for the identification of structural motifs of many other RNA molecules beyond transcriptional terminators.

Development and Application of Transcription Terminators for Polyhydroxylkanoates Production in Halophilic Halomonas bluephagenesis TD01

Halomonas bluephagenesis TD01 is one of the ideal chassis for low-cost industrial production based on “Next Generation Industrial Biotechnology,” yet the limited genetically regulatory parts such as transcriptional terminators, which are crucial for tuned regulations on gene expression, have hampered the engineering and applications of the strain. In this study, a series of intrinsic Rho-independent terminators were developed by either genome mining or rational design, and seven of them proved to exhibit higher efficiencies than the canonical strong T7 terminator, among which three terminators displayed high efficiencies over 90%. A preliminary modeling on the sequence-efficiency relationship of the terminators suggested that the poly U sequence regularity, the length and GC content of the stem, and the number and the size of hairpin loops remarkably affected the termination efficiency (TE). The rational and de novo designs of novel synthetic terminators based on the sequence-efficiency relationship and the “main contributor” engineering strategy proved to be effective, and fine-tuned polyhydroxylkanoates production was also achieved by the regulation of these native or synthetic terminators with different efficiencies. Furthermore, a perfectly positive correlation between the promoter activity and the TE was revealed in our study. The study enriches our knowledge of transcriptional termination via its sequence–strength relationship and enables the precise regulation of gene expression and PHA synthesis by intrinsic terminators, contributing to the extensive applications of H. bluephagenesis TD01 in the low-cost production of various chemicals.

Global RNA interactome of Salmonella discovers a 5' UTR sponge for the MicF small RNA that connects membrane permeability to transport capacity

The envelope of Gram-negative bacteria is a vital barrier that must balance protection and nutrient uptake. Small RNAs are crucial regulators of the envelope composition and function. Here, using RIL-seq to capture the Hfq-mediated RNA-RNA interactome in Salmonella enterica, we discover envelope-related riboregulators, including OppX. We show that OppX acts as an RNA sponge of MicF sRNA, a prototypical porin repressor. OppX originates from the 5' UTR of oppABCDF, encoding the major inner-membrane oligopeptide transporter, and sequesters MicF's seed region to derepress the synthesis of the porin OmpF. Intriguingly, OppX operates as a true sponge, storing MicF in an inactive complex without affecting its levels or stability. Conservation of the opp-OppX-MicF-ompF axis in related bacteria suggests that it serves an important mechanism, adjusting envelope porosity to specific transport capacity. These data also highlight the resource value of this Salmonella RNA interactome, which will aid in unraveling RNA-centric regulation in enteric pathogens.

Insights into 6S RNA in lactic acid bacteria (LAB)

BACKGROUND

6S RNA is a regulator of cellular transcription that tunes the metabolism of cells. This small non-coding RNA is found in nearly all bacteria and among the most abundant transcripts. Lactic acid bacteria (LAB) constitute a group of microorganisms with strong biotechnological relevance, often exploited as starter cultures for industrial products through fermentation. Some strains are used as probiotics while others represent potential pathogens. Occasional reports of 6S RNA within this group already indicate striking metabolic implications. A conceivable idea is that LAB with 6S RNA defects may metabolize nutrients faster, as inferred from studies of Echerichia coli. This may accelerate fermentation processes with the potential to reduce production costs. Similarly, elevated levels of secondary metabolites might be produced. Evidence for this possibility comes from preliminary findings regarding the production of surfactin in Bacillus subtilis, which has functions similar to those of bacteriocins. The prerequisite for its potential biotechnological utility is a general characterization of 6S RNA in LAB.

RESULTS

We provide a genomic annotation of 6S RNA throughout the Lactobacillales order. It laid the foundation for a bioinformatic characterization of common 6S RNA features. This covers secondary structures, synteny, phylogeny, and product RNA start sites. The canonical 6S RNA structure is formed by a central bulge flanked by helical arms and a template site for product RNA synthesis. 6S RNA exhibits strong syntenic conservation. It is usually flanked by the replication-associated recombination protein A and the universal stress protein A. A catabolite responsive element was identified in over a third of all 6S RNA genes. It is known to modulate gene expression based on the available carbon sources. The presence of antisense transcripts could not be verified as a general trait of LAB 6S RNAs.

CONCLUSIONS

Despite a large number of species and the heterogeneity of LAB, the stress regulator 6S RNA is well-conserved both from a structural as well as a syntenic perspective. This is the first approach to describe 6S RNAs and short 6S RNA-derived transcripts beyond a single species, spanning a large taxonomic group covering multiple families. It yields universal insights into this regulator and complements the findings derived from other bacterial model organisms.

Probing steps in DNA transcription using single-molecule methods

Transcriptional regulation is one of the key steps in determining gene expression. Diverse single-molecule techniques have been applied to characterize the stepwise progression of transcription, yielding complementary results. These techniques include, but are not limited to, fluorescence-based microscopy with single or multiple colors, force measuring and manipulating microscopy using magnetic field or light, and atomic force microscopy. Here, we summarize and evaluate these current methodologies in studying and resolving individual steps in the transcription reaction, which encompasses RNA polymerase binding, initiation, elongation, mRNA production, and termination. We also describe the advantages and disadvantages of each method for studying transcription.

Clusters of hairpins induce intrinsic transcription termination in bacteria

Intrinsic transcription termination (ITT) sites are currently identified by locating single and double-adjacent RNA hairpins downstream of the stop codon. ITTs for a limited number of genes/operons in only a few bacterial genomes are currently known. This lack of coverage is a lacuna in the existing ITT inference methods. We have studied the inter-operon regions of 13 genomes covering all major phyla in bacteria, for which good quality public RNA-seq data exist. We identify ITT sites in 87% of cases by predicting hairpin(s) and validate against 81% of cases for which the RNA-seq derived sites could be calculated. We identify 72% of these sites correctly, with 98% of them located ≤ 80 bases downstream of the stop codon. The predicted hairpins form a cluster (when present < 15 bases) in two-thirds of the cases, the remaining being single hairpins. The largest number of clusters is formed by two hairpins, and the occurrence decreases exponentially with an increasing number of hairpins in the cluster. Our study reveals that hairpins form an effective ITT unit when they act in concert in a cluster. Their pervasiveness along with single hairpin terminators corroborates a wider utilization of ITT mechanisms for transcription control across bacteria.

Predictable control of RNA lifetime using engineered degradation-tuning RNAs

The ability to tune RNA and gene expression dynamics is greatly needed for biotechnological applications. Native RNA stabilizers or engineered 5' stability hairpins have been used to regulate transcript half-life to control recombinant protein expression. However, these methods have been mostly ad hoc and hence lack predictability and modularity. Here, we report a library of RNA modules called degradation-tuning RNAs (dtRNAs) that can increase or decrease transcript stability in vivo and in vitro. dtRNAs enable modulation of transcript stability over a 40-fold dynamic range in Escherichia coli with minimal influence on translation initiation. We harness dtRNAs in messenger RNAs and noncoding RNAs to tune gene circuit dynamics and enhance CRISPR interference in vivo. Use of stabilizing dtRNAs in cell-free transcription-translation reactions also tunes gene and RNA aptamer production. Finally, we combine dtRNAs with toehold switch sensors to enhance the performance of paper-based norovirus diagnostics, illustrating the potential of dtRNAs for biotechnological applications.

Analysis of DnaK Expression from a Strain of Mycoplasma fermentans in Infected HCT116 Human Colon Carcinoma Cells

Several species of mycoplasmas, including Mycoplasma fermentans, are associated with certain human cancers. We previously isolated and characterized in our laboratory a strain of human mycoplasma M. fermentans subtype incognitus (MF-I1) able to induce lymphoma in a Severe Combined Immuno-Deficient (SCID) mouse model, and we demonstrated that its chaperone protein, DnaK, binds and reduces functions of human poly-ADP ribose polymerase-1 (PARP1) and ubiquitin carboxyl-terminal hydrolase protein-10 (USP10), which are required for efficient DNA repair and proper p53 activities, respectively. We also showed that other bacteria associated with human cancers (including Mycoplasmapneumoniae, Helicobacterpylori, Fusobacteriumnucleatum, Chlamydiathrachomatis, and Chlamydia pneumoniae) have closely related DnaK proteins, indicating a potential common mechanism of cellular transformation. Here, we quantify dnaK mRNA copy number by RT-qPCR analysis in different cellular compartments following intracellular MF-I1 infection of HCT116 human colon carcinoma cells. DnaK protein expression in infected cells was also detected and quantified by Western blot. The amount of viable intracellular mycoplasma reached a steady state after an initial phase of growth and was mostly localized in the cytoplasm of the invaded cells, while we detected a logarithmically increased number of viable extracellular bacteria. Our data indicate that, after invasion, MF-I1 is able to establish a chronic intracellular infection. Extracellular replication was more efficient while MF-I1 cultured in cell-free axenic medium showed a markedly reduced growth rate. We also identified modifications of important regulatory regions and heterogeneous lengths of dnaK mRNA transcripts isolated from intracellular and extracellular MF-I1. Both characteristics were less evident in dnaK mRNA transcripts isolated from MF-I1 grown in cell-free axenic media. Taken together, our data indicate that MF-I1, after establishing a chronic infection in eukaryotic cells, accumulates different forms of dnaK with efficient RNA turnover.

Evaluation and Comparison of the Efficiency of Transcription Terminators in Different Cyanobacterial Species

Cyanobacteria utilize sunlight to convert carbon dioxide into a wide variety of secondary metabolites and show great potential for green biotechnology applications. Although cyanobacterial synthetic biology is less mature than for other heterotrophic model organisms, there are now a range of molecular tools available to modulate and control gene expression. One area of gene regulation that still lags behind other model organisms is the modulation of gene transcription, particularly transcription termination. A vast number of intrinsic transcription terminators are now available in heterotrophs, but only a small number have been investigated in cyanobacteria. As artificial gene expression systems become larger and more complex, with short stretches of DNA harboring strong promoters and multiple gene expression cassettes, the need to stop transcription efficiently and insulate downstream regions from unwanted interference is becoming more important. In this study, we adapted a dual reporter tool for use with the CyanoGate MoClo Assembly system that can quantify and compare the efficiency of terminator sequences within and between different species. We characterized 34 intrinsic terminators in Escherichia coli, Synechocystis sp. PCC 6803, and Synechococcus elongatus UTEX 2973 and observed significant differences in termination efficiencies. However, we also identified five terminators with termination efficiencies of >96% in all three species, indicating that some terminators can behave consistently in both heterotrophic species and cyanobacteria.

Discovery and characterization of a fourth class of guanidine riboswitches

Riboswitches are RNAs that specifically sense a small molecule and regulate genes accordingly. The recent discovery of guanidine-binding riboswitches revealed the biological significance of this compound, and uncovered genes related to its biology. For example, certain sugE genes encode guanidine exporters and are activated by the riboswitches to reduce toxic levels of guanidine in the cell. In order to study guanidine biology and riboswitches, we applied a bioinformatics strategy for discovering additional guanidine riboswitches by searching for new candidate motifs associated with sugE genes. Based on in vitro and in vivo experiments, we determined that one of our six best candidates is a new structural class of guanidine riboswitches. The expression of a genetic reporter was induced 80-fold in response to addition of 5 mM guanidine in Staphylococcus aureus. This new class, called the guanidine-IV riboswitch, reveals additional guanidine-associated protein domains that are extremely rarely or never associated with previously established guanidine riboswitches. Among these protein domains are two transporter families that are structurally distinct from SugE, and could represent novel types of guanidine exporters. These results establish a new metabolite-binding RNA, further validate a bioinformatics method for finding riboswitches and suggest substrate specificities for as-yet uncharacterized transporter proteins.

Transcriptional Organization of the Stability Module of Broad-Host-Range Plasmid RA3, from the IncU Group

The broad-host-range (BHR) conjugative plasmids have developed diverse adaptive mechanisms defining the range of their promiscuity. The BHR conjugative RA3 plasmid, the archetype of the IncU group, can transfer between, replicate in, and be maintained in representatives of Alpha-, Beta-, and Gammaproteobacteria Its stability module encompasses ten open reading frames (ORFs) apparently organized into five operons, all transcribed in the same direction from several strong promoters that are tightly regulated either by autorepressors or by global plasmid-encoded regulators. In this paper, we demonstrate that owing to an efficient RNA polymerase (RNAP) read-through, the transcription from the first promoter, orf02p, may continue through the whole module. Moreover, an analysis of mRNA produced from the wild-type (WT) stability module and its deletion variants deprived of particular internal transcription initiation sites reveals that in fact each operon may be transcribed from any upstream promoter, giving rise to multicistronic transcripts of variable length and creating an additional level of gene expression control by transcript dosage adjustment. The gene expression patterns differ among various hosts, indicating that promoter recognition, regulation, and the RNAP read-through mechanisms are modulated in a species-specific manner.IMPORTANCE The efficiently disseminating conjugative or mobilizable BHR plasmids play key roles in the horizontal spread of genetic information between closely related and phylogenetically distant species, which can be harmful from the medical, veterinary, or industrial point of view. Understanding the mechanisms determining the plasmid's ability to function in diverse hosts is essential to help limit the spread of undesirable plasmid-encoded traits, e.g., antibiotic resistance. The range of a plasmid's promiscuity depends on the adaptations of its transfer, replication, and stability functions to the various hosts. IncU plasmids, with the archetype plasmid RA3, are considered to constitute a reservoir of antibiotic resistance genes in aquatic environments; however, the molecular mechanisms determining their adaptability to a broad range of hosts are rather poorly characterized. Here, we present the transcriptional organization of the stability module and show that the gene transcript dosage effect is an important determinant of the stable maintenance of RA3 in different hosts.

Precision design of stable genetic circuits carried in highly-insulated E. coli genomic landing pads

Genetic circuits have many applications, from guiding living therapeutics to ordering process in a bioreactor, but to be useful they have to be genetically stable and not hinder the host. Encoding circuits in the genome reduces burden, but this decreases performance and can interfere with native transcription. We have designed genomic landing pads in Escherichia coli at high-expression sites, flanked by ultrastrong double terminators. DNA payloads >8 kb are targeted to the landing pads using phage integrases. One landing pad is dedicated to carrying a sensor array, and two are used to carry genetic circuits. NOT/NOR gates based on repressors are optimized for the genome and characterized in the landing pads. These data are used, in conjunction with design automation software (Cello 2.0), to design circuits that perform quantitatively as predicted. These circuits require fourfold less RNA polymerase than when carried on a plasmid and are stable for weeks in a recA+ strain without selection. This approach enables the design of synthetic regulatory networks to guide cells in environments or for applications where plasmid use is infeasible.

The minimal meningococcal ProQ protein has an intrinsic capacity for structure-based global RNA recognition

FinO-domain proteins are a widespread family of bacterial RNA-binding proteins with regulatory functions. Their target spectrum ranges from a single RNA pair, in the case of plasmid-encoded FinO, to global RNA regulons, as with enterobacterial ProQ. To assess whether the FinO domain itself is intrinsically selective or promiscuous, we determine in vivo targets of Neisseria meningitidis, which consists of solely a FinO domain. UV-CLIP-seq identifies associations with 16 small non-coding sRNAs and 166 mRNAs. Meningococcal ProQ predominantly binds to highly structured regions and generally acts to stabilize its RNA targets. Loss of ProQ alters transcript levels of >250 genes, demonstrating that this minimal ProQ protein impacts gene expression globally. Phenotypic analyses indicate that ProQ promotes oxidative stress resistance and DNA damage repair. We conclude that FinO domain proteins recognize some abundant type of RNA shape and evolve RNA binding selectivity through acquisition of additional regions that constrain target recognition.

FinO-domain proteins are bacterial RNA-binding proteins with a wide range of target specificities. Here, the authors employ UV CLIP-seq and show that minimal ProQ protein of Neisseria meningitidis binds to various small non-coding RNAs and mRNAs involved in virulence.

Navigation through the twists and turns of RNA sequencing technologies: Application to bacterial regulatory RNAs

Discovered in the 1980s, small regulatory RNAs (sRNAs) are now considered key actors in virtually all aspects of bacterial physiology and virulence. Together with transcriptional and translational regulatory proteins, they integrate and often are hubs of complex regulatory networks, responsible for bacterial response/adaptation to various perceived stimuli. The recent development of powerful RNA sequencing technologies has facilitated the identification and characterization of sRNAs (length, structure and expression conditions) and their RNA targets in several bacteria. Nevertheless, it could be very difficult for non-experts to understand the advantages and drawbacks related to each offered option and, consequently, to make an informed choice. Therefore, the main goal of this review is to provide a guide to navigate through the twists and turns of high-throughput RNA sequencing technologies, with a specific focus on those applied to the study of sRNAs. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Small RNA Mcr11 requires the transcription factor AbmR for stable expression and regulates genes involved in the central metabolism of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the etiologic agent of tuberculosis, must adapt to host-associated environments during infection by modulating gene expression. Small regulatory RNAs (sRNAs) are key regulators of bacterial gene expression, but their roles in Mtb are not well understood. Here, we address the expression and function of the Mtb sRNA Mcr11, which is associated with slow bacterial growth and chronic infections in mice. We found that stable expression of Mcr11 requires multiple factors specific to TB-complex bacteria, including the AbmR transcription factor. Bioinformatic analyses used to predict regulatory targets of Mcr11 identified 7-11 nucleotide regions with potential for direct base-pairing with Mcr11 immediately upstream of Rv3282, fadA3, and lipB. mcr11-dependent regulation of these genes was demonstrated using qRT-PCR and found to be responsive to the presence of fatty acids. Mutation of the putative Mcr11 base-pairing site upstream of lipB in a promoter reporter strain resulted in significant de-repression of lipB expression, similar to that observed in mcr11-deleted Mtb. These studies establish Mcr11's roles in regulating growth and central metabolism in Mtb. Our finding that multiple TB-complex-specific factors are required for production of stable Mcr11 also emphasizes the need to better understand mechanisms of sRNA expression and stability in TB.

Revisiting intrinsic transcription termination in mycobacteria: U-tract downstream of secondary structure is dispensable for termination

In Escherichia coli, G/C-rich hairpin structure followed by a U-tract in the 3' region of the nascent RNA are crucial determinants for intrinsic or factor independent transcription termination. In mycobacteria, there is a scarcity of such intrinsic terminators. However, secondary structures having G/C-rich stem devoid of any U's or with suboptimal U-tracts were identified earlier as terminators and found to be functional both in vitro and in vivo. Two different observations - that a mycobacterial RNA polymerase (RNAP) does not function at intrinsic terminators devoid of U-tracts and the identification of an altogether new motif for termination in mycobacteria necessitated re-examining a number of putative terminators for their function as terminators. When these in silico identified non-canonical terminators were subjected to experimental validation, they were found to dissociate RNA from the elongating RNAP. Termination is observed when the U-tracts were reduced, or totally absent both in vitro and in vivo. Our results, thus indicate that the presence of U-tract following the G/C-rich stem in an intrinsic terminator may not be an essential determinant for transcription termination in mycobacteria.

A multivariate prediction model for Rho-dependent termination of transcription

Bacterial transcription termination proceeds via two main mechanisms triggered either by simple, well-conserved (intrinsic) nucleic acid motifs or by the motor protein Rho. Although bacterial genomes can harbor hundreds of termination signals of either type, only intrinsic terminators are reliably predicted. Computational tools to detect the more complex and diversiform Rho-dependent terminators are lacking. To tackle this issue, we devised a prediction method based on Orthogonal Projections to Latent Structures Discriminant Analysis [OPLS-DA] of a large set of in vitro termination data. Using previously uncharacterized genomic sequences for biochemical evaluation and OPLS-DA, we identified new Rho-dependent signals and quantitative sequence descriptors with significant predictive value. Most relevant descriptors specify features of transcript C>G skewness, secondary structure, and richness in regularly-spaced 5'CC/UC dinucleotides that are consistent with known principles for Rho-RNA interaction. Descriptors collectively warrant OPLS-DA predictions of Rho-dependent termination with a ∼85% success rate. Scanning of the Escherichia coli genome with the OPLS-DA model identifies significantly more termination-competent regions than anticipated from transcriptomics and predicts that regions intrinsically refractory to Rho are primarily located in open reading frames. Altogether, this work delineates features important for Rho activity and describes the first method able to predict Rho-dependent terminators in bacterial genomes.

High-resolution RNA 3'-ends mapping of bacterial Rho-dependent transcripts

Transcription termination in bacteria can occur either via Rho-dependent or independent (intrinsic) mechanisms. Intrinsic terminators are composed of a stem-loop RNA structure followed by a uridine stretch and are known to terminate in a precise manner. In contrast, Rho-dependent terminators have more loosely defined characteristics and are thought to terminate in a diffuse manner. While transcripts ending in an intrinsic terminator are protected from 3′-5′ exonuclease digestion due to the stem-loop structure of the terminator, it remains unclear what protects Rho-dependent transcripts from being degraded. In this study, we mapped the exact steady-state RNA 3′ ends of hundreds of Escherichia coli genes terminated either by Rho-dependent or independent mechanisms. We found that transcripts generated from Rho-dependent termination have precise 3′-ends at steady state. These termini were localized immediately downstream of energetically stable stem-loop structures, which were not followed by uridine rich sequences. We provide evidence that these structures protect Rho-dependent transcripts from 3′-5′ exonucleases such as PNPase and RNase II, and present data localizing the Rho-utilization (rut) sites immediately downstream of these protective structures. This study represents the first extensive in-vivo map of exact RNA 3′-ends of Rho-dependent transcripts in E. coli.

Prediction and diversity of tracrRNAs from type II CRISPR-Cas systems

Type II CRISPR-Cas9 systems require a small RNA called the trans-activating CRISPR RNA (tracrRNA) in order to function. The prediction of these non-coding RNAs in prokaryotic genomes is challenging because they have dissimilar structures, having short stems (3-6 bp) and non-canonical base-pairs e.g. G-A. Much of the tracrRNA is involved in base-pairing interactions with the CRISPR RNA, or itself, or in RNA-protein interactions with Cas9. Here we develop a new bioinformatic tool to predict tracrRNAs. On an experimentally verified test set the algorithm achieved a high sensitivity and specificity, and a low false discovery rate (FDR) on genome analysis. Analysis of representative RefSeq genomes (5462) detected 275 tracrRNAs from 165 genera. These tracrRNAs could be grouped into 15 clusters which were used to build covariance models. These clusters included Streptococci and Staphylococci tracrRNAs from the CRISPR-Cas9 systems which are currently used for gene editing. Compensating base changes observed in the models were consistent with the experimental structures of single guide RNAs (sgRNAs). Other clusters, for which there are not yet structures available, were predicted to form novel tracrRNA folds. These clusters included a large and divergent tracrRNA set from Bacteroidetes. These computational models contribute to the understanding of CRISPR-Cas biology, and will assist in the design of further engineered CRISPR-Cas9 systems. The tracrRNA prediction software is available through a galaxy web server.

RhoTermPredict: an algorithm for predicting Rho-dependent transcription terminators based on Escherichia coli, Bacillus subtilis and Salmonella enterica databases

BACKGROUND

In bacterial genomes, there are two mechanisms to terminate the DNA transcription: the "intrinsic" or Rho-independent termination and the Rho-dependent termination. Intrinsic terminators are characterized by a RNA hairpin followed by a run of 6-8 U residues relatively easy to identify using one of the numerous available prediction programs. In contrast, Rho-dependent termination is mediated by the Rho protein factor that, firstly, binds to ribosome-free mRNA in a site characterized by a C > G content and then reaches the RNA polymerase to induce its release. Conversely on intrinsic terminators, the computational prediction of Rho-dependent terminators in prokaryotes is a very difficult problem because the sequence features required for the function of Rho are complex and poorly defined. This is the reason why it still does not exist an exhaustive Rho-dependent terminators prediction program.

RESULTS

In this study we introduce RhoTermPredict, the first published algorithm for an exhaustive Rho-dependent terminators prediction in bacterial genomes. RhoTermPredict identifies these elements based on a previously proposed consensus motif common to all Rho-dependent transcription terminators. It essentially searches for a 78 nt long RUT site characterized by a C > G content and with regularly spaced C residues, followed by a putative pause site for the RNA polymerase. We tested RhoTermPredict performances by using available genomic and transcriptomic data of the microorganism Escherichia coli K-12, both in limited-length sequences and in the whole-genome, and available genomic sequences from Bacillus subtilis 168 and Salmonella enterica LT2 genomes. We also estimated the overlap between the predictions of RhoTermPredict and those obtained by the predictor of intrinsic terminators ARNold webtool. Our results demonstrated that RhoTermPredict is a very performing algorithm both for limited-length sequences (F₁-score obtained about 0.7) and for a genome-wide analysis. Furthermore the degree of overlap with ARNold predictions was very low.

CONCLUSIONS

Our analysis shows that RhoTermPredict is a powerful tool for Rho-dependent terminators search in the three analyzed genomes and could fill this gap in computational genomics. We conclude that RhoTermPredict could be used in combination with an intrinsic terminators predictor in order to predict all the transcription terminators in bacterial genomes.

A Fluorogenic Assay To Monitor Rho-Dependent Termination of Transcription

Transcription termination mediated by the ring-shaped, ATP-dependent Rho motor is a multipurpose regulatory mechanism specific to bacteria and constitutes an interesting target for the development of new antibiotics. Although Rho-dependent termination can punctuate gene expression or contribute to the protection of the genome at hundreds of sites within a given bacterium, its exact perimeter and site- or species-specific features remain insufficiently characterized. New advanced approaches are required to explore thoroughly the diversity of Rho-dependent terminators and the complexity of associated mechanisms. Current in vitro analyses of Rho-dependent termination rely on radiolabeling, gel electrophoresis, and phosphorimaging of transcription reaction products and are thus hazardous, inconvenient, and low-throughput. To address these limitations, we have developed the first in vitro assay using a fluorescence detection modality to study Rho-dependent transcription termination. This powerful experimental tool accurately estimates terminator strengths in a matter of minutes and is optimized for a microplate reader format allowing multiplexed characterization of putative terminator sequences and mechanisms or high-throughput screening of new drugs targeting Rho-dependent termination.

Computational operon prediction in whole-genomes and metagenomes

Microbial diversity in unique environmental settings enables abrupt responses catalysed by altering the gene regulation and formation of gene clusters called operons. Operons increases bacterial adaptability, which in turn increases their survival. This review article presents the emergence of computational operon prediction methods for whole microbial genomes and metagenomes, and discusses their strengths and limitations. Most of the whole-genome operon prediction methods struggle to generalize on unrelated genomes. The applicability of universal whole-genome operon prediction methods to metagenomic data is an interesting yet less investigated question. We have evaluated the potential of various operon prediction features for genomic and metagenomic data. Most of operon prediction methods with high accuracy have been compiled into databases. Despite of the high predictive performance, the data among many databases are not completely consistent for similar species. We performed a correlation analysis between the computationally predicted operon databases and experimentally validated data for Escherichia coli, Bacillus subtilis and Mycobacterium tuberculosis. Operon prediction for most of the less characterized microbes cannot be verified due to absence of experimentally validated operons. The generation of validated information for other microbes would test the authenticity of operon databases for other less annotated microbes as well. Advances in sequencing technologies and development of better analysis methods will help researchers to overcome the technological hurdles (such as long sequencing reads and improved contig size) and further improve operon predictions and better utilize operonic information.

Part by Part: Synthetic Biology Parts Used in Solventogenic Clostridia

The solventogenic Clostridia are of interest to the chemical industry because of their natural ability to produce chemicals such as butanol, acetone and ethanol from diverse feedstocks. Their use as whole cell factories presents multiple metabolic engineering targets that could lead to improved sustainability and profitability of Clostridium industrial processes. However, engineering efforts have been held back by the scarcity of genetic and synthetic biology tools. Over the past decade, genetic tools to enable transformation and chromosomal modifications have been developed, but the lack of a broad palette of synthetic biology parts remains one of the last obstacles to the rapid engineered improvement of these species for bioproduction. We have systematically reviewed existing parts that have been used in the modification of solventogenic Clostridia, revealing a narrow range of empirically chosen and nonengineered parts that are in current use. The analysis uncovers elements, such as promoters, transcriptional terminators and ribosome binding sites where increased fundamental knowledge is needed for their reliable use in different applications. Together, the review provides the most comprehensive list of parts used and also presents areas where an improved toolbox is needed for full exploitation of these industrially important bacteria.

Detection of 224 candidate structured RNAs by comparative analysis of specific subsets of intergenic regions

The discovery of structured non-coding RNAs (ncRNAs) in bacteria can reveal new facets of biology and biochemistry. Comparative genomics analyses executed by powerful computer algorithms have successfully been used to uncover many novel bacterial ncRNA classes in recent years. However, this general search strategy favors the discovery of more common ncRNA classes, whereas progressively rarer classes are correspondingly more difficult to identify. In the current study, we confront this problem by devising several methods to select subsets of intergenic regions that can concentrate these rare RNA classes, thereby increasing the probability that comparative sequence analysis approaches will reveal their existence. By implementing these methods, we discovered 224 novel ncRNA classes, which include ROOL RNA, an RNA class averaging 581 nt and present in multiple phyla, several highly conserved and widespread ncRNA classes with properties that suggest sophisticated biochemical functions and a multitude of putative cis-regulatory RNA classes involved in a variety of biological processes. We expect that further research on these newly found RNA classes will reveal additional aspects of novel biology, and allow for greater insights into the biochemistry performed by ncRNAs.

The primary transcriptome of Neisseria meningitidis and its interaction with the RNA chaperone Hfq

Neisseria meningitidis is a human commensal that can also cause life-threatening meningitis and septicemia. Despite growing evidence for RNA-based regulation in meningococci, their transcriptome structure and output of regulatory small RNAs (sRNAs) are incompletely understood. Using dRNA-seq, we have mapped at single-nucleotide resolution the primary transcriptome of N. meningitidis strain 8013. Annotation of 1625 transcriptional start sites defines transcription units for most protein-coding genes but also reveals a paucity of classical σ70-type promoters, suggesting the existence of activators that compensate for the lack of −35 consensus sequences in N. meningitidis. The transcriptome maps also reveal 65 candidate sRNAs, a third of which were validated by northern blot analysis. Immunoprecipitation with the RNA chaperone Hfq drafts an unexpectedly large post-transcriptional regulatory network in this organism, comprising 23 sRNAs and hundreds of potential mRNA targets. Based on this data, using a newly developed gfp reporter system we validate an Hfq-dependent mRNA repression of the putative colonization factor PrpB by the two trans-acting sRNAs RcoF1/2. Our genome-wide RNA compendium will allow for a better understanding of meningococcal transcriptome organization and riboregulation with implications for colonization of the human nasopharynx.

Depleting Mycobacterium tuberculosis of the transcription termination factor Rho causes pervasive transcription and rapid death

Rifampicin, which inhibits bacterial RNA polymerase, provides one of the most effective treatments for tuberculosis. Inhibition of the transcription termination factor Rho is used to treat some bacterial infections, but its importance varies across bacteria. Here we show that Rho of Mycobacterium tuberculosis functions to both define the 3' ends of mRNAs and silence substantial fragments of the genome. Brief inactivation of Rho affects over 500 transcripts enriched for genes of foreign DNA elements and bacterial virulence factors. Prolonged inactivation of Rho causes extensive pervasive transcription, a genome-wide increase in antisense transcripts, and a rapid loss of viability of replicating and non-replicating M. tuberculosis in vitro and during acute and chronic infection in mice. Collectively, these data suggest that inhibition of Rho may provide an alternative strategy to treat tuberculosis with an efficacy similar to inhibition of RNA polymerase.

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

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

IMPORTANCE

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

Studying RNA Homology and Conservation with Infernal: From Single Sequences to RNA Families

Emerging high-throughput technologies have led to a deluge of putative non-coding RNA (ncRNA) sequences identified in a wide variety of organisms. Systematic characterization of these transcripts will be a tremendous challenge. Homology detection is critical to making maximal use of functional information gathered about ncRNAs: identifying homologous sequence allows us to transfer information gathered in one organism to another quickly and with a high degree of confidence. ncRNA presents a challenge for homology detection, as the primary sequence is often poorly conserved and de novo secondary structure prediction and search remain difficult. This unit introduces methods developed by the Rfam database for identifying "families" of homologous ncRNAs starting from single "seed" sequences, using manually curated sequence alignments to build powerful statistical models of sequence and structure conservation known as covariance models (CMs), implemented in the Infernal software package. We provide a step-by-step iterative protocol for identifying ncRNA homologs and then constructing an alignment and corresponding CM. We also work through an example for the bacterial small RNA MicA, discovering a previously unreported family of divergent MicA homologs in genus Xenorhabdus in the process. © 2016 by John Wiley & Sons, Inc.

Regulation of mRNA Decay in Bacteria

Gram-negative and gram-positive bacteria use a variety of enzymatic pathways to degrade mRNAs. Although several recent reviews have outlined these pathways, much less attention has been paid to the regulation of mRNA decay. The functional half-life of a particular mRNA, which affects how much protein is synthesized from it, is determined by a combination of multiple factors. These include, but are not necessarily limited to, (a) stability elements at either the 5' or the 3' terminus, (b) posttranscriptional modifications, (c) ribosome density on individual mRNAs, (d) small regulatory RNA (sRNA) interactions with mRNAs, (e) regulatory proteins that alter ribonuclease binding affinities, (f) the presence or absence of endonucleolytic cleavage sites, (g) control of intracellular ribonuclease levels, and (h) physical location within the cell. Changes in physiological conditions associated with environmental alterations can significantly alter the impact of these factors in the decay of a particular mRNA.

Novel Utilization of Terminators in the Design of Biologically Adjustable Synthetic Filters

Terminators, which signal the end of transcription processes, are typically placed behind the last coding sequence of an operon to prevent interference between transcript units in most biologically synthetic systems. Here, we seek to extend the usability of terminators in genetic system design by using terminators as regulatory genetic parts. Terminators with different impacts on their upstream and downstream genes are characterized in detail via dynamic modeling to predict the behavior of the overall genetic system. Some nonlinear effects of terminators observed in our terminator measurements potentially facilitate regulation of gene expression. Through dynamic modeling in silico, we find that such genetic systems may behave like genetic filters. In agreement with the simulations, we successfully implement genetic high-pass and bandpass filters in vivo, demonstrating the potential of using terminators as regulatory parts. The genetic bandpass filter in this work is implemented through the interdependence between genetic parts, in which the termination efficiency of a terminator varies with the strength of the upstream promoter. This design strategy for a bandpass filter requires fewer base pairs than the conventional strategy of concatenating high-pass and low-pass filters. Our results show that this novel utilization of terminators as regulatory parts may provide a new perspective for efficient design of genetic circuits. We believe that further exploration of the complicated dynamics of terminators is important in the development of synthetic biology.

Global RNA recognition patterns of post-transcriptional regulators Hfq and CsrA revealed by UV crosslinking in vivo

The molecular roles of many RNA‐binding proteins in bacterial post‐transcriptional gene regulation are not well understood. Approaches combining in vivo UV crosslinking with RNA deep sequencing (CLIP‐seq) have begun to revolutionize the transcriptome‐wide mapping of eukaryotic RNA‐binding protein target sites. We have applied CLIP‐seq to chart the target landscape of two major bacterial post‐transcriptional regulators, Hfq and CsrA, in the model pathogen Salmonella Typhimurium. By detecting binding sites at single‐nucleotide resolution, we identify RNA preferences and structural constraints of Hfq and CsrA during their interactions with hundreds of cellular transcripts. This reveals 3′‐located Rho‐independent terminators as a universal motif involved in Hfq–RNA interactions. Additionally, Hfq preferentially binds 5′ to sRNA‐target sites in mRNAs, and 3′ to seed sequences in sRNAs, reflecting a simple logic in how Hfq facilitates sRNA–mRNA interactions. Importantly, global knowledge of Hfq sites significantly improves sRNA‐target predictions. CsrA binds AUGGA sequences in apical loops and targets many Salmonella virulence mRNAs. Overall, our generic CLIP‐seq approach will bring new insights into post‐transcriptional gene regulation by RNA‐binding proteins in diverse bacterial species.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

New classes of self-cleaving ribozymes revealed by comparative genomics analysis

Enzymes made of RNA catalyze reactions that are essential for protein synthesis and RNA processing. However, such natural ribozymes are exceedingly rare, as evident by the fact that the discovery rate for new classes has dropped to one per decade from about one per year during the 1980s. Indeed, only 11 distinct ribozyme classes have been experimentally validated to date. Recently, we recognized that self-cleaving ribozymes frequently associate with certain types of genes from bacteria. Herein this synteny was exploited to identify divergent architectures for two previously known ribozyme classes and to discover additional noncoding RNA motifs that are self-cleaving RNA candidates. Three new self-cleaving classes, named twister sister, pistol and hatchet, have been identified from this collection, suggesting that even more ribozymes remain hidden in modern cells.

An empirical strategy to detect bacterial transcript structure from directional RNA-seq transcriptome data

Background

As sequencing costs are being lowered continuously, RNA-seq has gradually been adopted as the first choice for comparative transcriptome studies with bacteria. Unlike microarrays, RNA-seq can directly detect cDNA derived from mRNA transcripts at a single nucleotide resolution. Not only does this allow researchers to determine the absolute expression level of genes, but it also conveys information about transcript structure. Few automatic software tools have yet been established to investigate large-scale RNA-seq data for bacterial transcript structure analysis.

Results

In this study, 54 directional RNA-seq libraries from Salmonella serovar Typhimurium (S. Typhimurium) 14028s were examined for potential relationships between read mapping patterns and transcript structure. We developed an empirical method, combined with statistical tests, to automatically detect key transcript features, including transcriptional start sites (TSSs), transcriptional termination sites (TTSs) and operon organization. Using our method, we obtained 2,764 TSSs and 1,467 TTSs for 1331 and 844 different genes, respectively. Identification of TSSs facilitated further discrimination of 215 putative sigma 38 regulons and 863 potential sigma 70 regulons. Combining the TSSs and TTSs with intergenic distance and co-expression information, we comprehensively annotated the operon organization in S. Typhimurium 14028s.

Conclusions

Our results show that directional RNA-seq can be used to detect transcriptional borders at an acceptable resolution of ±10-20 nucleotides. Technical limitations of the RNA-seq procedure may prevent single nucleotide resolution. The automatic transcript border detection methods, statistical models and operon organization pipeline that we have described could be widely applied to RNA-seq studies in other bacteria. Furthermore, the TSSs, TTSs, operons, promoters and unstranslated regions that we have defined for S. Typhimurium 14028s may constitute valuable resources that can be used for comparative analyses with other Salmonella serotypes.

Electronic supplementary material

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

Annotating RNA motifs in sequences and alignments

RNA performs a diverse array of important functions across all cellular life. These functions include important roles in translation, building translational machinery and maturing messenger RNA. More recent discoveries include the miRNAs and bacterial sRNAs that regulate gene expression, the thermosensors, riboswitches and other cis-regulatory elements that help prokaryotes sense their environment and eukaryotic piRNAs that suppress transposition. However, there can be a long period between the initial discovery of a RNA and determining its function. We present a bioinformatic approach to characterize RNA motifs, which are critical components of many RNA structure-function relationships. These motifs can, in some instances, provide researchers with functional hypotheses for uncharacterized RNAs. Moreover, we introduce a new profile-based database of RNA motifs--RMfam--and illustrate some applications for investigating the evolution and functional characterization of RNA. All the data and scripts associated with this work are available from: https://github.com/ppgardne/RMfam.