Differential RNA-seq: the approach behind and the biological insight gained

Biological Insights from RNA-RNA Interactomes in Bacteria, as Revealed by RIL-seq

Bacteria reside in constantly changing environments and require rapid and precise adjustments of gene expression to ensure survival. Small regulatory RNAs (sRNAs) are a crucial element that bacteria utilize to achieve this. sRNAs are short RNA molecules that modulate gene expression usually through base-pairing interactions with target RNAs, primarily mRNAs. These interactions can lead to either negative outcomes such as mRNA degradation or translational repression or positive outcomes such as mRNA stabilization or translation enhancement. In recent years, high-throughput approaches such as RIL-seq (RNA interaction by ligation and sequencing) revolutionized the sRNA field by enabling the identification of sRNA targets on a global scale, unveiling intricate sRNA-RNA networks. In this review, we discuss the insights gained from investigating sRNA-RNA networks in well-studied bacterial species as well as in understudied bacterial species. Having a complete understanding of sRNA-mediated regulation is critical for the development of new strategies for controlling bacterial growth and combating bacterial infections.

Quantitative mapping of pseudouridines in bacteria RNA

RNA pseudouridylation is one of the most prevalent post-transcriptional modifications, occurring universally across all organisms. Although pseudouridines have been extensively studied in bacterial tRNAs and rRNAs, their presence and role in bacterial mRNA remain poorly characterized. Here, we used a bisulfite-based sequencing approach to provide a comprehensive and quantitative measurement of bacteria pseudouridines. As a proof of concept in E. coli, we identified 1,954 high-confidence sites in 1,331 transcripts, covering almost 30% of the transcriptome. Furthermore, pseudouridine mapping enabled the detection of differentially expressed genes associated with stress response that were unidentified using conventional RNA-seq approach. We also demonstrate that in addition to pseudouridine profiling, our approach can facilitate the discovery of previously unidentified transcripts. As an example, we identified a small RNA transcribed from the antisense strand of tRNA-Tyr which represses expression of distal genes. Finally, we mapped pseudouridines in oral microbiome samples of human subjects, demonstrating the broad applicability of our approach in complex microbiomes. Altogether, our work highlights the advantages of mapping bacterial pseudouridines and provides a tool to study posttranscription regulation in microbial communities.

Selective degradation of phage RNAs by the Csm6 ribonuclease provides robust type III CRISPR immunity in Streptococcus thermophilus

Type III CRISPR immune systems bind viral or plasmid RNA transcripts and activate Csm3/Cmr4 and Cas10 nucleases to uniquely cleave both invader RNA and DNA, respectively. Additionally, type III effector complexes generate cyclic oligoadenylate (cOA) signaling molecules to activate trans-acting, auxiliary Csm6/Csx1 ribonucleases, previously proposed to be non-specific in their in vivo RNA cleavage preference. Despite extensive in vitro studies, the nuclease requirements of type III systems in their native contexts remain poorly understood. Here we systematically investigated the in vivo roles for immunity of each of the three Streptococcus thermophilus (Sth) type III-A Cas nucleases and cOA signaling by challenging nuclease defective mutant strains with plasmid and phage infections. Our results reveal that RNA cleavage by Csm6 is both sufficient and essential for maintaining wild-type levels of immunity. Importantly, Csm6 RNase activity leads to immunity against even high levels of phage challenge without causing host cell dormancy or death. Transcriptomic analyses during phage infection indicated Csm6-mediated and crRNA-directed preferential cleavage of phage transcripts. Our findings highlight the critical role of Csm6 RNase activity in type III immunity and demonstrate specificity for invader RNA transcripts by Csm6 to ensure host cell survival upon phage infection.

De novo gene synthesis by an antiviral reverse transcriptase

Defense-associated reverse transcriptase (DRT) systems perform DNA synthesis to protect bacteria against viral infection, but the identities and functions of their DNA products remain largely unknown. We show that DRT2 systems encode an unprecedented immune pathway that involves de novo gene synthesis through rolling circle reverse transcription of a noncoding RNA (ncRNA). Programmed template jumping on the ncRNA generates a concatemeric cDNA, which becomes double-stranded upon viral infection. This DNA product constitutes a protein-coding, nearly endless open reading frame (neo) gene whose expression leads to potent cell growth arrest, restricting the viral infection. Our work highlights an elegant expansion of genome coding potential through RNA-templated gene creation and challenges conventional paradigms of genetic information encoded along the one-dimensional axis of genomic DNA.

GraphPro: An interpretable graph neural network-based model for identifying promoters in multiple species

Promoters are DNA sequences that bind with RNA polymerase to initiate transcription, regulating this process through interactions with transcription factors. Accurate identification of promoters is crucial for understanding gene expression regulation mechanisms and developing therapeutic approaches for various diseases. However, experimental techniques for promoter identification are often expensive, time-consuming, and inefficient, necessitating the development of accurate and efficient computational models for this task. Enhancing the model's ability to recognize promoters across multiple species and improving its interpretability pose significant challenges. In this study, we introduce a novel interpretable model based on graph neural networks, named GraphPro, for multi-species promoter identification. Initially, we encode the sequences using k-tuple nucleotide frequency pattern, dinucleotide physicochemical properties, and dna2vec. Subsequently, we construct two feature extraction modules based on convolutional neural networks and graph neural networks. These modules aim to extract specific motifs from the promoters, learn their dependencies, and capture the underlying structural features of the promoters, providing a more comprehensive representation. Finally, a fully connected neural network predicts whether the input sequence is a promoter. We conducted extensive experiments on promoter datasets from eight species, including Human, Mouse, and Escherichia coli. The experimental results show that the average Sn, Sp, Acc and MCC values of GraphPro are 0.9123, 0.9482, 0.8840 and 0.7984, respectively. Compared with previous promoter identification methods, GraphPro not only achieves better recognition accuracy on multiple species, but also outperforms all previous methods in cross-species prediction ability. Furthermore, by visualizing GraphPro's decision process and analyzing the sequences matching the transcription factor binding motifs captured by the model, we validate its significant advantages in biological interpretability. The source code for GraphPro is available at https://github.com/liuliwei1980/GraphPro.

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coli and Staphylococcus aureus and extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis.

Multimodal HLA-I genotype regulation by human cytomegalovirus US10 and resulting surface patterning

Human leucocyte antigen class I (HLA-I) molecules play a central role for both NK and T-cell responses that prevent serious human cytomegalovirus (HCMV) disease. To create opportunities for viral spread, several HCMV-encoded immunoevasins employ diverse strategies to target HLA-I. Among these, the glycoprotein US10 is so far insufficiently studied. While it was reported that US10 interferes with HLA-G expression, its ability to manipulate classical HLA-I antigen presentation remains unknown. In this study, we demonstrate that US10 recognizes and binds to all HLA-I (HLA-A, -B, -C, -E, -G) heavy chains. Additionally, impaired recruitment of HLA-I to the peptide loading complex was observed. Notably, the associated effects varied significantly dependending on HLA-I genotype and allotype: (i) HLA-A molecules evaded downregulation by US10, (ii) tapasin-dependent HLA-B molecules showed impaired maturation and cell surface expression, and (iii) β2m-assembled HLA-C, in particular HLA-C*05:01 and -C*12:03, and HLA-G were strongly retained in complex with US10 in the endoplasmic reticulum. These genotype-specific effects on HLA-I were confirmed through unbiased HLA-I ligandome analyses. Furthermore, in HCMV-infected fibroblasts inhibition of overlapping US10 and US11 transcription had little effect on HLA-A, but induced HLA-B antigen presentation. Thus, the US10-mediated impact on HLA-I results in multiple geno- and allotypic effects in a so far unparalleled and multimodal manner.

An expanded transcriptome atlas for Bacteroides thetaiotaomicron reveals a small RNA that modulates tetracycline sensitivity

Plasticity in gene expression allows bacteria to adapt to diverse environments. This is particularly relevant in the dynamic niche of the human intestinal tract; however, transcriptional networks remain largely unknown for gut-resident bacteria. Here we apply differential RNA sequencing (RNA-seq) and conventional RNA-seq to the model gut bacterium Bacteroides thetaiotaomicron to map transcriptional units and profile their expression levels across 15 in vivo-relevant growth conditions. We infer stress- and carbon source-specific transcriptional regulons and expand the annotation of small RNAs (sRNAs). Integrating this expression atlas with published transposon mutant fitness data, we predict conditionally important sRNAs. These include MasB, which downregulates tetracycline tolerance. Using MS2 affinity purification and RNA-seq, we identify a putative MasB target and assess its role in the context of the MasB-associated phenotype. These data-publicly available through the Theta-Base web browser ( http://micromix.helmholtz-hiri.de/bacteroides/ )-constitute a valuable resource for the microbiome community.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

Genome-Scale Analysis of the Structure and Function of RNA Pathways and Networks in Pseudomonas aeruginosa

In recent years, several genome-wide approaches based on RNA sequencing (RNA-seq) have been developed. These methods allow a comprehensive and dynamic view of the structure and function of the multi-layered RNA pathways and networks. Many of these approaches, including the promising one of single-cell transcriptome analysis, have been successfully applied to Pseudomonas aeruginosa. However, we are only at the beginning because only a few surrounding conditions have been considered. Here, we aim to illustrate the different types of approaches based on RNA-seq that will lead us in the future to a better understanding of the dynamics of RNA biology in P. aeruginosa.

Discovering Novel Bacterial Small RNA by RNA-seq Analysis Toolkit ANNOgesic

ANNOgesic is an RNA-seq analysis pipeline that can detect sRNAs and many other genomic features in bacteria and archaea. In addition to listing sRNA candidates, ANNOgesic also generates various formats of data files for visual examination and downstream experimental design. Based on validations from previous studies, the sRNA predictions are accurate and reliable. In this chapter, we outline the sRNA detection algorithm, important parameters used, step-by-step execution, and data interpretation with a B. pertussis study as an example. Following those procedures, novel sRNA can be revealed by ANNOgesic.

Genome-wide transcriptional response to silver stress in extremely halophilic archaeon Haloferax alexandrinus DSM 27206 T

BACKGROUND

The extremely halophilic archaeon Haloferax (Hfx.) alexandrinus DSM 27206 T was previously documented for the ability to biosynthesize silver nanoparticles while mechanisms underlying its silver tolerance were overlooked. In the current study, we aimed to assess the transcriptional response of this haloarchaeon to varying concentrations of silver, seeking a comprehensive understanding of the molecular determinants underpinning its heavy metal tolerance.

RESULTS

The growth curves confirmed the capacity of Hfx. alexandrinus to surmount silver stress, while the SEM-EDS analysis illustrated the presence of silver nanoparticles in cultures exposed to 0.5 mM silver nitrate. The RNA-Seq based transcriptomic analysis of Hfx. alexandrinus cells exposed to 0.1, 0.25, and 0.5 mM silver nitrate revealed the differential expression of multiple sets of genes potentially employed in heavy-metal stress response, genes mostly related to metal transporters, basic metabolism, oxidative stress response and cellular motility. The RT-qPCR analysis of selected transcripts was conducted to verify and validate the generated RNA-Seq data.

CONCLUSIONS

Our results indicated that copA, encoding the copper ATPase, is essential for the survival of Hfx. alexandrinus cells in silver-containing saline media. The silver-exposed cultures underwent several metabolic adjustments that enabled the activation of enzymes involved in the oxidative stress response and impairment of the cellular movement capacity. To our knowledge, this study represents the first comprehensive analysis of gene expression in halophillic archaea facing increased levels of heavy metals.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

Identification of NAD-RNA species and ADPR-RNA decapping in Archaea

NAD is a coenzyme central to metabolism that also serves as a 5'-terminal cap for bacterial and eukaryotic transcripts. Thermal degradation of NAD can generate nicotinamide and ADP-ribose (ADPR). Here, we use LC-MS/MS and NAD captureSeq to detect and identify NAD-RNAs in the thermophilic model archaeon Sulfolobus acidocaldarius and in the halophilic mesophile Haloferax volcanii. None of the four Nudix proteins of S. acidocaldarius catalyze NAD-RNA decapping in vitro, but one of the proteins (Saci_NudT5) promotes ADPR-RNA decapping. NAD-RNAs are converted into ADPR-RNAs, which we detect in S. acidocaldarius total RNA. Deletion of the gene encoding the 5'-3' exonuclease Saci-aCPSF2 leads to a 4.5-fold increase in NAD-RNA levels. We propose that the incorporation of NAD into RNA acts as a degradation marker for Saci-aCPSF2. In contrast, ADPR-RNA is processed by Saci_NudT5 into 5'-p-RNAs, providing another layer of regulation for RNA turnover in archaeal cells.

The environmentally-regulated interplay between local three-dimensional chromatin organisation and transcription of proVWX in E. coli

Nucleoid associated proteins (NAPs) maintain the architecture of bacterial chromosomes and regulate gene expression. Thus, their role as transcription factors may involve three-dimensional chromosome re-organisation. While this model is supported by in vitro studies, direct in vivo evidence is lacking. Here, we use RT-qPCR and 3C-qPCR to study the transcriptional and architectural profiles of the H-NS (histone-like nucleoid structuring protein)-regulated, osmoresponsive proVWX operon of Escherichia coli at different osmolarities and provide in vivo evidence for transcription regulation by NAP-mediated chromosome re-modelling in bacteria. By consolidating our in vivo investigations with earlier in vitro and in silico studies that provide mechanistic details of how H-NS re-models DNA in response to osmolarity, we report that activation of proVWX in response to a hyperosmotic shock involves the destabilization of H-NS-mediated bridges anchored between the proVWX downstream and upstream regulatory elements (DRE and URE), and between the DRE and ygaY that lies immediately downstream of proVWX. The re-establishment of these bridges upon adaptation to hyperosmolarity represses the operon. Our results also reveal additional structural features associated with changes in proVWX transcript levels such as the decompaction of local chromatin upstream of the operon, highlighting that further complexity underlies the regulation of this model operon. H-NS and H-NS-like proteins are wide-spread amongst bacteria, suggesting that chromosome re-modelling may be a typical feature of transcriptional control in bacteria.

Methodologies for bacterial ribonuclease characterization using RNA-seq

Bacteria adjust gene expression at the post-transcriptional level through an intricate network of small regulatory RNAs and RNA-binding proteins, including ribonucleases (RNases). RNases play an essential role in RNA metabolism, regulating RNA stability, decay, and activation. These enzymes exhibit species-specific effects on gene expression, bacterial physiology, and different strategies of target recognition. Recent advances in high-throughput RNA sequencing (RNA-seq) approaches have provided a better understanding of the roles and modes of action of bacterial RNases. Global studies aiming to identify direct targets of RNases have highlighted the diversity of RNase activity and RNA-based mechanisms of gene expression regulation. Here, we review recent RNA-seq approaches used to study bacterial RNases, with a focus on the methods for identifying direct RNase targets.

Low Molecular Weight and High Deacetylation Degree Chitosan Batch Alleviates Pathogenesis, Toxin Accumulation, and Fusarium Gene Regulation in Barley Leaf Pathosystem

Fusarium graminearum is a cosmopolitan fungal pathogen that destroys cereal production, in terms of loss of yield and grain contamination with mycotoxins, worldwide. Chitosan is a natural biopolymer abundant in the environment with proven antifungal properties that also acts as a plant immunity elicitor. Despite a number of articles, there is a lack of systematic comparison of antifungal activity of diverse batches of chitosan. The current study aimed to test the inhibitory effects of a collection of diverse chitosan samples on the growth and production of F. graminearum toxins, validated by changes in the Fusarium transcriptome. Experiments included testing antifungal activity of different chitosan samples, the application of the best performing one in vitro to investigate the impact on F. graminearum growth, followed by analyzing its effect on Fusarium toxins accumulation, and Fusarium transcriptomics in the barley leaf pathosystem. Confirmatory antifungal assays revealed that CS_10, a specific batch of chitosan, retarded Fusarium growth with an application concentration of 200 ppm, significantly reducing toxin synthesis and disease symptoms in Fusarium-inoculated barley leaves. RNA-Seq analysis of F. graminearum in barley leaf pathosystem exposed to CS_10 showed a list of differentially expressed genes involved in redox balance, cell respiration, nutrient transport, cell wall degradation enzymes, ergosterol biosynthesis, and trichothecenes production. The genes functioning in these essential pathways are discussed and assigned as critical checkpoints to control Fusarium infections. The results suggest some important molecular targets in F. graminearum that may be suitable in gene-specific targeting or transgene-free methods, such as spray-induced gene silencing during host-pathogen interactions.

Characterization of Bacterial Transcriptional Regulatory Networks in Escherichia coli through Genome-Wide In Vitro Run-Off Transcription/RNA-seq (ROSE)

The global characterization of transcriptional regulatory networks almost exclusively uses in vivo conditions, thereby providing a snapshot on multiple regulatory interactions at the same time. To complement these approaches, we developed and applied a method for characterizing bacterial promoters genome-wide by in vitro transcription coupled to transcriptome sequencing specific for native 5'-ends of transcripts. This method, called ROSE (run-off transcription/RNA-sequencing), only requires chromosomal DNA, ribonucleotides, RNA polymerase (RNAP) core enzyme, and a specific sigma factor, recognizing the corresponding promoters, which have to be analyzed. ROSE was performed on E. coli K-12 MG1655 genomic DNA using Escherichia coli RNAP holoenzyme (including σ70) and yielded 3226 transcription start sites, 2167 of which were also identified in in vivo studies, and 598 were new. Many new promoters not yet identified by in vivo experiments might be repressed under the tested conditions. Complementary in vivo experiments with E. coli K-12 strain BW25113 and isogenic transcription factor gene knockout mutants of fis, fur, and hns were used to test this hypothesis. Comparative transcriptome analysis demonstrated that ROSE could identify bona fide promoters that were apparently repressed in vivo. In this sense, ROSE is well-suited as a bottom-up approach for characterizing transcriptional networks in bacteria and ideally complementary to top-down in vivo transcriptome studies.

Decoding murine cytomegalovirus

The genomes of both human cytomegalovirus (HCMV) and murine cytomegalovirus (MCMV) were first sequenced over 20 years ago. Similar to HCMV, the MCMV genome had initially been proposed to harbor ≈170 open reading frames (ORFs). More recently, omics approaches revealed HCMV gene expression to be substantially more complex comprising several hundred viral ORFs. Here, we provide a state-of-the art reannotation of lytic MCMV gene expression based on integrative analysis of a large set of omics data. Our data reveal 365 viral transcription start sites (TiSS) that give rise to 380 and 454 viral transcripts and ORFs, respectively. The latter include >200 small ORFs, some of which represented the most highly expressed viral gene products. By combining TiSS profiling with metabolic RNA labelling and chemical nucleotide conversion sequencing (dSLAM-seq), we provide a detailed picture of the expression kinetics of viral transcription. This not only resulted in the identification of a novel MCMV immediate early transcript encoding the m166.5 ORF, which we termed ie4, but also revealed a group of well-expressed viral transcripts that are induced later than canonical true late genes and contain an initiator element (Inr) but no TATA- or TATT-box in their core promoters. We show that viral upstream ORFs (uORFs) tune gene expression of longer viral ORFs expressed in cis at translational level. Finally, we identify a truncated isoform of the viral NK-cell immune evasin m145 arising from a viral TiSS downstream of the canonical m145 mRNA. Despite being ≈5-fold more abundantly expressed than the canonical m145 protein it was not required for downregulating the NK cell ligand, MULT-I. In summary, our work will pave the way for future mechanistic studies on previously unknown cytomegalovirus gene products in an important virus animal model.

Antisense transcription and its roles in adaption to environmental stress in E. coli

It has been reported that a highly varying proportion (1% ∼ 93%) of genes in various prokaryotes have antisense RNA (asRNA) transcription. However, the extent of the pervasiveness of asRNA transcription in the well-studied E. coli K12 strain has thus far been an issue of debate. Furthermore, very little is known about the expression patterns and functions of asRNAs under various conditions. To fill these gaps, we determined the transcriptomes and proteomes of E. coli K12 at multiple time points in five culture conditions using strand-specific RNA-seq, differential RNA-seq, and quantitative mass spectrometry methods. To reduce artifacts of possible transcriptional noise, we identified asRNA using stringent criteria with biological replicate verification and transcription start sites (TSSs) information included. We identified a total of 660 asRNAs, which were generally short and largely condition-dependently transcribed. We found that the proportions of the genes which had asRNA transcription highly depended on the culture conditions and time points. We classified the transcriptional activities of the genes in six transcriptional modes according to their relative levels of asRNA to mRNA. Many genes changed their transcriptional modes at different time points of the culture conditions, and such transitions can be described in a well-defined manner. Intriguingly, the protein levels and mRNA levels of genes in the sense-only/sense-dominant mode were moderately correlated, but the same was not true for genes in the balanced/antisense-dominant mode, in which asRNAs were at a comparable or higher level to mRNAs. These observations were further validated by western blot on candidate genes, where an increase in asRNA transcription diminished gene expression in one case and enhanced it in another. These results suggest that asRNAs may directly or indirectly regulate translation by forming duplexes with cognate mRNAs. Thus, asRNAs may play an important role in the bacterium's responses to environmental changes during growth and adaption to different environments.

IMPORTANCE

The cis -antisense RNA (asRNA) is a type of understudied RNA molecules in prokaryotes, which is believed to be important in regulating gene expression. Our current understanding of asRNA is constrained by inconsistent reports about its identification and properties. These discrepancies are partially caused by a lack of sufficient samples, biological replicates, and culture conditions. This study aimed to overcome these disadvantages and identified 660 putative asRNAs using integrated information from strand-specific RNA-seq, differential RNA-seq, and mass spectrometry methods. In addition, we explored the relative expression between asRNAs and sense RNAs and investigated asRNA regulated transcriptional activity changes over different culture conditions and time points. Our work strongly suggests that asRNAs may play a crucial role in bacterium's responses to environmental changes during growth and adaption to different environments.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

RNA Preparation and RNA-Seq Bioinformatics for Comparative Transcriptomics

The principal transcriptome analysis is the determination of differentially expressed genes across experimental conditions. For this, the next-generation sequencing of RNA (RNA-seq) has several advantages over other techniques, such as the capability of detecting all the transcripts in one assay over RT-qPCR, such as its higher accuracy and broader dynamic range over microarrays or the ability to detect novel transcripts, including non-coding RNA molecules, at nucleotide-level resolution over both techniques. Despite these advantages, many microbiology laboratories have not yet applied RNA-seq analyses to their investigations. The high cost of the equipment for next-generation sequencing is no longer an issue since this intermediate part of the analysis can be provided by commercial or central services. Here, we detail a protocol for the first part of the analysis, the RNA extraction and an introductory protocol to the bioinformatics analysis of the sequencing data to generate the differential expression results.

Escherichia coli transcriptome assembly from a compendium of RNA-seq data sets

Owing to the complexities of bacterial RNA biology, the transcriptomes of even the best studied bacteria are not fully understood. To help elucidate the transcriptional landscape of E. coli, we compiled a compendium of 3,376 RNA-seq data sets composed of more than 7 trillion sequenced bases, which we evaluate with a transcript assembly pipeline. We report expression profiles for all annotated E. coli genes as well as 5,071 other transcripts. Additionally, we observe hundreds of instances of co-transcribed genes that are novel with respect to existing operon databases. By integrating data from a large number of sequencing experiments corresponding to a wide range of conditions, we are able to obtain a comprehensive view of the E. coli transcriptome.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Bacterial Transcription Factors Bind to Coding Regions and Regulate Internal Cryptic Promoters

Transcription factors (TFs) regulate transcription by binding to the specific sequences at the promoter region. However, the mechanisms and functions of TFs binding within the coding sequences (CDS) remain largely elusive in prokaryotes. To this end, we collected 409 data sets for bacterial TFs, including 104 chromatin immunoprecipitation sequencing (ChIP-seq) assays and 305 data sets from the systematic evolution of ligands by exponential enrichment (SELEX) in seven model bacteria. Interestingly, these TFs displayed the same binding capabilities for both coding and intergenic regions. Subsequent biochemical and genetic experiments demonstrated that several TFs bound to the coding regions and regulated the transcription of the binding or adjacent genes. Strand-specific RNA sequencing revealed that these CDS-binding TFs regulated the activity of the cryptic promoters, resulting in the altered transcription of the corresponding antisense RNA. TF RhpR hindered the transcriptional elongation of a subgenic transcript within a CDS. A ChIP-seq and Ribo-seq coanalysis revealed that RhpR influenced the translational efficiency of binding genes. Taken together, the present study reveals three regulatory mechanisms of CDS-bound TFs within individual genes, operons, and antisense RNAs, which demonstrate the variability of the regulatory mechanisms of TFs and expand upon the complexity of bacterial transcriptomes.

SAPPHIRE.CNN: Implementation of dRNA-seq-driven, species-specific promoter prediction using convolutional neural networks

Data availability is a consistent bottleneck for the development of bacterial species-specific promoter prediction software. In this work we leverage genome-wide promoter datasets generated with dRNA-seq in the Gram-negative bacteria Pseudomonas aeruginosa and Salmonella enterica for promoter prediction. Convolutional neural networks are presented as an optimal architecture for model training and are further modified and tailored for promoter prediction. The resulting predictors reach high binary accuracies (95% and 94.9%) on test sets and outperform each other when predicting promoters in their associated species. SAPPHIRE.CNN is available online and can also be downloaded to run locally. Our results indicate a dependency of binary promoter classification on an organism’s GC content and a decreased performance of our classifiers on genera they were not trained for, further supporting the need for dedicated, species-specific promoter classification tools.

Key players in regulatory RNA realm of bacteria

Precise regulation of gene expression is crucial for living cells to adapt for survival in diverse environmental conditions. Among the common cellular regulatory mechanisms, RNA-based regulators play a key role in all domains of life. Discovery of regulatory RNAs have made a paradigm shift in molecular biology as many regulatory functions of RNA have been identified beyond its canonical roles as messenger, ribosomal and transfer RNA. In the complex regulatory RNA network, riboswitches, small RNAs, and RNA thermometers can be identified as some of the key players. Herein, we review the discovery, mechanism, and potential therapeutic use of these classes of regulatory RNAs mainly found in bacteria. Being highly adaptive organisms that inhabit a broad range of ecological niches, bacteria have adopted tight and rapid-responding gene regulation mechanisms. This review aims to highlight how bacteria utilize versatile RNA structures and sequences to build a sophisticated gene regulation network.

•

The three major classes of prokaryotic ncRNAs and their characterized mechanisms of operation in gene regulation.

•

sRNAs emerging as major players in global gene regulatory networks.

•

Riboswitch mediated gene control mechanisms through on/off switches in response to ligand binding.

•

RNA thermo sensors for temperature-dependent gene expression.

•

Therapeutic importance of ncRNAs and computational approaches involved in the discovery of ncRNAs.

Mapping the Complex Transcriptional Landscape of the Phytopathogenic Bacterium Dickeya dadantii

Dickeya dadantii is a phytopathogenic bacterium that causes soft rot in a wide range of plant hosts worldwide and a model organism for studying virulence gene regulation. The present study provides a comprehensive and annotated transcriptomic map of D. dadantii obtained by a computational method combining five independent transcriptomic data sets: (i) paired-end RNA sequencing (RNA-seq) data for a precise reconstruction of the RNA landscape; (ii) DNA microarray data providing transcriptional responses to a broad variety of environmental conditions; (iii) long-read Nanopore native RNA-seq data for isoform-level transcriptome validation and determination of transcription termination sites; (iv) differential RNA sequencing (dRNA-seq) data for the precise mapping of transcription start sites; (v) in planta DNA microarray data for a comparison of gene expression profiles between in vitro experiments and the early stages of plant infection. Our results show that transcription units sometimes coincide with predicted operons but are generally longer, most of them comprising internal promoters and terminators that generate alternative transcripts of variable gene composition. We characterize the occurrence of transcriptional read-through at terminators, which might play a basal regulation role and explain the extent of transcription beyond the scale of operons. We finally highlight the presence of noncontiguous operons and excludons in the D. dadantii genome, novel genomic arrangements that might contribute to the basal coordination of transcription. The highlighted transcriptional organization may allow D. dadantii to finely adjust its gene expression program for a rapid adaptation to fast-changing environments. IMPORTANCE This is the first transcriptomic map of a Dickeya species. It may therefore significantly contribute to further progress in the field of phytopathogenicity. It is also one of the first reported applications of long-read Nanopore native RNA-seq in prokaryotes. Our findings yield insights into basal rules of coordination of transcription that might be valid for other bacteria and may raise interest in the field of microbiology in general. In particular, we demonstrate that gene expression is coordinated at the scale of transcription units rather than operons, which are larger functional genomic units capable of generating transcripts with variable gene composition for a fine-tuning of gene expression in response to environmental changes. In line with recent studies, our findings indicate that the canonical operon model is insufficient to explain the complexity of bacterial transcriptomes.

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

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

RNase III-CLASH of multi-drug resistant Staphylococcus aureus reveals a regulatory mRNA 3'UTR required for intermediate vancomycin resistance

Treatment of methicillin-resistant Staphylococcus aureus infections is dependent on the efficacy of last-line antibiotics including vancomycin. Treatment failure is commonly linked to isolates with intermediate vancomycin resistance (termed VISA). These isolates have accumulated point mutations that collectively reduce vancomycin sensitivity, often by thickening the cell wall. Changes in regulatory small RNA expression have been correlated with antibiotic stress in VISA isolates however the functions of most RNA regulators is unknown. Here we capture RNA-RNA interactions associated with RNase III using CLASH. RNase III-CLASH uncovers hundreds of novel RNA-RNA interactions in vivo allowing functional characterisation of many sRNAs for the first time. Surprisingly, many mRNA-mRNA interactions are recovered and we find that an mRNA encoding a long 3' untranslated region (UTR) (termed vigR 3'UTR) functions as a regulatory 'hub' within the RNA-RNA interaction network. We demonstrate that the vigR 3'UTR promotes expression of folD and the cell wall lytic transglycosylase isaA through direct mRNA-mRNA base-pairing. Deletion of the vigR 3'UTR re-sensitised VISA to glycopeptide treatment and both isaA and vigR 3'UTR deletions impact cell wall thickness. Our results demonstrate the utility of RNase III-CLASH and indicate that S. aureus uses mRNA-mRNA interactions to co-ordinate gene expression more widely than previously appreciated.

Development of ONT-cappable-seq to unravel the transcriptional landscape of Pseudomonas phages

RNA sequencing has become the method of choice to study the transcriptional landscape of phage-infected bacteria. However, short-read RNA sequencing approaches generally fail to capture the primary 5' and 3' boundaries of transcripts, confounding the discovery of key transcription initiation and termination events as well as operon architectures. Yet, the elucidation of these elements is crucial for the understanding of the strategy of transcription regulation during the infection process, which is currently lacking beyond a handful of model phages. We developed ONT-cappable-seq, a specialized long-read RNA sequencing technique that allows end-to-end sequencing of primary prokaryotic transcripts using the Nanopore sequencing platform. We applied ONT-cappable-seq to study transcription of Pseudomonas aeruginosa phage LUZ7, obtaining a comprehensive genome-wide map of viral transcription start sites, terminators, and complex operon structures that fine-regulate gene expression. Our work provides new insights in the RNA biology of a non-model phage, unveiling distinct promoter architectures, putative small non-coding viral RNAs, and the prominent regulatory role of terminators during infection. The robust workflow presented here offers a framework to obtain a global, yet fine-grained view of phage transcription and paves the way for standardized, in-depth transcription studies for microbial viruses or bacteria in general.

Selective inhibition of miRNA processing by a herpesvirus-encoded miRNA

Herpesviruses have mastered host cell modulation and immune evasion to augment productive infection, life-long latency and reactivation^1,2. A long appreciated, yet undefined relationship exists between the lytic-latent switch and viral non-coding RNAs^3,4. Here we identify viral microRNA (miRNA)-mediated inhibition of host miRNA processing as a cellular mechanism that human herpesvirus 6A (HHV-6A) exploits to disrupt mitochondrial architecture, evade intrinsic host defences and drive the switch from latent to lytic virus infection. We demonstrate that virus-encoded miR-aU14 selectively inhibits the processing of multiple miR-30 family members by direct interaction with the respective primary (pri)-miRNA hairpin loops. Subsequent loss of miR-30 and activation of the miR-30-p53-DRP1 axis triggers a profound disruption of mitochondrial architecture. This impairs induction of type I interferons and is necessary for both productive infection and virus reactivation. Ectopic expression of miR-aU14 triggered virus reactivation from latency, identifying viral miR-aU14 as a readily druggable master regulator of the herpesvirus lytic-latent switch. Our results show that miRNA-mediated inhibition of miRNA processing represents a generalized cellular mechanism that can be exploited to selectively target individual members of miRNA families. We anticipate that targeting miR-aU14 will provide new therapeutic options for preventing herpesvirus reactivations in HHV-6-associated disorders.

The tracrRNA in CRISPR Biology and Technologies

CRISPR-Cas adaptive immune systems in bacteria and archaea utilize short CRISPR RNAs (crRNAs) to guide sequence-specific recognition and clearance of foreign genetic material. Multiple crRNAs are stored together in a compact format called a CRISPR array that is transcribed and processed into the individual crRNAs. While the exact processing mechanisms vary widely, some CRISPR-Cas systems, including those encoding the Cas9 nuclease, rely on a trans-activating crRNA (tracrRNA). The tracrRNA was discovered in 2011 and was quickly co-opted to create single-guide RNAs as core components of CRISPR-Cas9 technologies. Since then, further studies have uncovered processes extending beyond the traditional role of tracrRNA in crRNA biogenesis, revealed Cas nucleases besides Cas9 that are dependent on tracrRNAs, and established new applications based on tracrRNA engineering. In this review, we describe the biology of the tracrRNA and how its ongoing characterization has garnered new insights into prokaryotic immune defense and enabled key technological advances.

Promotech: a general tool for bacterial promoter recognition

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

Comparative genomics provides structural and functional insights into Bacteroides RNA biology

Bacteria employ noncoding RNA molecules for a wide range of biological processes, including scaffolding large molecular complexes, catalyzing chemical reactions, defending against phages, and controlling gene expression. Secondary structures, binding partners, and molecular mechanisms have been determined for numerous small noncoding RNAs (sRNAs) in model aerobic bacteria. However, technical hurdles have largely prevented analogous analyses in the anaerobic gut microbiota. While experimental techniques are being developed to investigate the sRNAs of gut commensals, computational tools and comparative genomics can provide immediate functional insight. Here, using Bacteroides thetaiotaomicron as a representative microbiota member, we illustrate how comparative genomics improves our understanding of the RNA biology in an understudied gut bacterium. We investigate putative RNA-binding proteins and predict a Bacteroides cold-shock protein homologue to have an RNA-related function. We apply an in-silico protocol incorporating both sequence and structural analysis to determine the consensus structures and conservation of nine Bacteroides noncoding RNA families. Using structure probing, we validate and refine these predictions, and deposit them in the Rfam database. Through synteny analyses, we illustrate how genomic co-conservation can serve as a predictor of sRNA function. Altogether, this work showcases the power of RNA informatics for investigating the RNA biology of anaerobic microbiota members.

RNA landscape of the emerging cancer-associated microbe Fusobacterium nucleatum

Fusobacterium nucleatum, long known as a constituent of the oral microflora, has recently garnered renewed attention for its association with several different human cancers. The growing interest in this emerging cancer-associated bacterium contrasts with a paucity of knowledge about its basic gene expression features and physiological responses. As fusobacteria lack all established small RNA-associated proteins, post-transcriptional networks in these bacteria are also unknown. In the present study, using differential RNA-sequencing, we generate high-resolution global RNA maps for five clinically relevant fusobacterial strains-F. nucleatum subspecies nucleatum, animalis, polymorphum and vincentii, as well as F. periodonticum-for early, mid-exponential growth and early stationary phase. These data are made available in an online browser, and we use these to uncover fundamental aspects of fusobacterial gene expression architecture and a suite of non-coding RNAs. Developing a vector for functional analysis of fusobacterial genes, we discover a conserved fusobacterial oxygen-induced small RNA, FoxI, which serves as a post-transcriptional repressor of the major outer membrane porin FomA. Our findings provide a crucial step towards delineating the regulatory networks enabling F. nucleatum adaptation to different environments, which may elucidate how these bacteria colonize different compartments of the human body.

Halobacterium salinarum and Haloferax volcanii Comparative Transcriptomics Reveals Conserved Transcriptional Processing Sites

Post-transcriptional processing of messenger RNA is an important regulatory strategy that allows relatively fast responses to changes in environmental conditions. In halophile systems biology, the protein perspective of this problem (i.e., ribonucleases which implement the cleavages) is generally more studied than the RNA perspective (i.e., processing sites). In the present in silico work, we mapped genome-wide transcriptional processing sites (TPS) in two halophilic model organisms, Halobacterium salinarum NRC-1 and Haloferax volcanii DS2. TPS were established by reanalysis of publicly available differential RNA-seq (dRNA-seq) data, searching for non-primary (monophosphorylated RNAs) enrichment. We found 2093 TPS in 43% of H. salinarum genes and 3515 TPS in 49% of H. volcanii chromosomal genes. Of the 244 conserved TPS sites found, the majority were located around start and stop codons of orthologous genes. Specific genes are highlighted when discussing antisense, ribosome and insertion sequence associated TPS. Examples include the cell division gene ftsZ2, whose differential processing signal along growth was detected and correlated with post-transcriptional regulation, and biogenesis of sense overlapping transcripts associated with IS200/IS605. We hereby present the comparative, transcriptomics-based processing site maps with a companion browsing interface.

A Computational Framework for Identifying Promoter Sequences in Nonmodel Organisms Using RNA-seq Data Sets

Engineering microorganisms into biological factories that convert renewable feedstocks into valuable materials is a major goal of synthetic biology; however, for many nonmodel organisms, we do not yet have the genetic tools, such as suites of strong promoters, necessary to effectively engineer them. In this work, we developed a computational framework that can leverage standard RNA-seq data sets to identify sets of constitutive, strongly expressed genes and predict strong promoter signals within their upstream regions. The framework was applied to a diverse collection of RNA-seq data measured for the methanotroph Methylotuvimicrobium buryatense 5GB1 and identified 25 genes that were constitutively, strongly expressed across 12 experimental conditions. For each gene, the framework predicted short (27-30 nucleotide) sequences as candidate promoters and derived -35 and -10 consensus promoter motifs (TTGACA and TATAAT, respectively) for strong expression in M. buryatense. This consensus closely matches the canonical E. coli sigma-70 motif and was found to be enriched in promoter regions of the genome. A subset of promoter predictions was experimentally validated in a XylE reporter assay, including the consensus promoter, which showed high expression. The pmoC, pqqA, and ssrA promoter predictions were additionally screened in an experiment that scrambled the -35 and -10 signal sequences, confirming that transcription initiation was disrupted when these specific regions of the predicted sequence were altered. These results indicate that the computational framework can make biologically meaningful promoter predictions and identify key pieces of regulatory systems that can serve as foundational tools for engineering diverse microorganisms for biomolecule production.

An RNA-centric global view of Clostridioides difficile reveals broad activity of Hfq in a clinically important gram-positive bacterium

The gram-positive human pathogen Clostridioides difficile has emerged as the leading cause of antibiotic-associated diarrhea. However, little is known about the bacterium's transcriptome architecture and mechanisms of posttranscriptional control. Here, we have applied transcription start site and termination mapping to generate a single-nucleotide-resolution RNA map of C. difficile 5' and 3' untranslated regions, operon structures, and noncoding regulators, including 42 sRNAs. Our results indicate functionality of many conserved riboswitches and predict cis-regulatory RNA elements upstream of multidrug resistance (MDR)-type ATP-binding cassette (ABC) transporters and transcriptional regulators. Despite growing evidence for a role of Hfq in RNA-based gene regulation in C. difficile, the functions of Hfq-based posttranscriptional regulatory networks in gram-positive pathogens remain controversial. Using Hfq immunoprecipitation followed by sequencing of bound RNA species (RIP-seq), we identify a large cohort of transcripts bound by Hfq and show that absence of Hfq affects transcript stabilities and steady-state levels. We demonstrate sRNA expression during intestinal colonization by C. difficile and identify infection-related signals impacting its expression. As a proof of concept, we show that the utilization of the abundant intestinal metabolite ethanolamine is regulated by the Hfq-dependent sRNA CDIF630nc_085. Overall, our study lays the foundation for understanding clostridial riboregulation with implications for the infection process and provides evidence for a global role of Hfq in posttranscriptional regulation in a gram-positive bacterium.

Shaping the Bacterial Epitranscriptome-5'-Terminal and Internal RNA Modifications

All domains of life utilize a diverse set of modified ribonucleotides that can impact the sequence, structure, function, stability, and the fate of RNAs, as well as their interactions with other molecules. Today, more than 160 different RNA modifications are known that decorate the RNA at the 5'-terminus or internal RNA positions. The boost of next-generation sequencing technologies sets the foundation to identify and study the functional role of RNA modifications. The recent advances in the field of RNA modifications reveal a novel regulatory layer between RNA modifications and proteins, which is central to developing a novel concept called "epitranscriptomics." The majority of RNA modifications studies focus on the eukaryotic epitranscriptome. In contrast, RNA modifications in prokaryotes are poorly characterized. This review outlines the current knowledge of the prokaryotic epitranscriptome focusing on mRNA modifications. Here, it is described that several internal and 5'-terminal RNA modifications either present or likely present in prokaryotic mRNA. Thereby, the individual techniques to identify these epitranscriptomic modifications, their writers, readers and erasers, and their proposed functions are explored. Besides that, still unanswered questions in the field of prokaryotic epitranscriptomics are pointed out, and its future perspectives in the dawn of next-generation sequencing technologies are outlined.

PPD: A Manually Curated Database for Experimentally Verified Prokaryotic Promoters

As a key region, promoter plays a key role in transcription regulation. A eukaryotic promoter database called EPD has been constructed to store eukaryotic POL II promoters. Although there are some promoter databases for specific prokaryotic species or specific promoter type, such as RegulonDB for Escherichia coli K-12, DBTBS for Bacillus subtilis and Pro54DB for sigma 54 promoter, because of the diversity of prokaryotes and the development of sequencing technology, huge amounts of prokaryotic promoters are scattered in numerous published articles, which is inconvenient for researchers to explore the process of gene regulation in prokaryotes. In this study, we constructed a Prokaryotic Promoter Database (PPD), which records the experimentally validated promoters in prokaryotes, from published articles. Up to now, PPD has stored 129,148 promoters across 63 prokaryotic species manually extracted from published papers. We provided a friendly interface for users to browse, search, blast, visualize, submit and download data. The PPD will provide relatively comprehensive resources of prokaryotic promoter for the study of prokaryotic gene transcription. The PPD is freely available and easy accessed at http://lin-group.cn/database/ppd/.

Coupled Transcriptomics for Differential Expression Analysis and Determination of Transcription Start Sites: Design and Bioinformatics

The term coupled transcriptomics is coined to describe a design of an RNA-seq experiment intended for both differential expression analysis and genome-wide determination of the transcription start sites (TSS). The minimal requirements for the first analysis are two experimental conditions with at least two biological replicates enabling statistical tests. The second analysis involves the bioinformatics comparison of the data generated from a control RNA-seq library with another library enriched in primary transcripts using Terminator^™ 5'-phosphate-dependent exonuclease, in an experiment denominated differential RNA-seq (dRNA-seq). Usually, dRNA-seq is carried out with specific protocols for library construction, different of those used for common differential expression analysis. Our experimental design allows to use the same data for both analyses, reducing the number of libraries to be generated and sequenced. This is a guide for designing a coupled transcriptomics experiment and for the subsequent bioinformatics procedures. The proposed methods can be applied to the detection and study of small RNA genes.

Tetrachloroethene respiration in Sulfurospirillum species is regulated by a two-component system as unraveled by comparative genomics, transcriptomics, and regulator binding studies

Energy conservation via organohalide respiration (OHR) in dehalogenating Sulfurospirillum species is an inducible process. However, the gene products involved in tetrachloroethene (PCE) sensing and signal transduction have not been unambiguously identified. Here, genome sequencing of Sulfurospirillum strains defective in PCE respiration and comparative genomics, which included the PCE-respiring representatives of the genus, uncovered the genetic inactivation of a two-component system (TCS) in the OHR gene region of the natural mutants. The assumption that the TCS gene products serve as a PCE sensor that initiates gene transcription was supported by the constitutive low-level expression of the TCS operon in fumarate-adapted cells of Sulfurospirillum multivorans. Via RNA sequencing, eight transcriptional units were identified in the OHR gene region, which includes the TCS operon, the PCE reductive dehalogenase operon, the gene cluster for norcobamide biosynthesis, and putative accessory genes with unknown functions. The OmpR-family response regulator (RR) encoded in the TCS operon was functionally characterized by promoter-binding assays. The RR bound a cis-regulatory element that contained a consensus sequence of a direct repeat (CTATW) separated by 17 bp. Its location either overlapping the -35 box or 50 bp further upstream indicated different regulatory mechanisms. Sequence variations in the regulator binding sites identified in the OHR gene region were in accordance with differences in the transcript levels of the respective gene clusters forming the PCE regulon. The results indicate the presence of a fine-tuned regulatory network controlling PCE metabolism in dehalogenating Sulfurospirillum species, a group of metabolically versatile organohalide-respiring bacteria.

Single-Nucleotide RNA Maps for the Two Major Nosocomial Pathogens Enterococcus faecalis and Enterococcus faecium

Enterococcus faecalis and faecium are two major representative clinical strains of the Enterococcus genus and are sadly notorious to be part of the top agents responsible for nosocomial infections. Despite their critical implication in worldwide public healthcare, essential and available resources such as deep transcriptome annotations remain poor, which also limits our understanding of post-transcriptional control small regulatory RNA (sRNA) functions in these bacteria. Here, using the dRNA-seq technique in combination with ANNOgesic analysis, we successfully mapped and annotated transcription start sites (TSS) of both E. faecalis V583 and E. faecium AUS0004 at single nucleotide resolution. Analyzing bacteria in late exponential phase, we capture ~40% (E. faecalis) and 43% (E. faecium) of the annotated protein-coding genes, determine 5′ and 3′ UTR (untranslated region) length, and detect instances of leaderless mRNAs. The transcriptome maps revealed sRNA candidates in both bacteria, some found in previous studies and new ones. Expression of candidate sRNAs is being confirmed under biologically relevant environmental conditions. This comprehensive global TSS mapping atlas provides a valuable resource for RNA biology and gene expression analysis in the Enterococci. It can be accessed online at www.helmholtz-hiri.de/en/datasets/enterococcus through an instance of the genomic viewer JBrowse.

Introducing differential RNA-seq mapping to track the early infection phase for Pseudomonas phage ɸKZ

As part of the ongoing renaissance of phage biology, more phage genomes are becoming available through DNA sequencing. However, our understanding of the transcriptome architecture that allows these genomes to be expressed during host infection is generally poor. Transcription start sites (TSSs) and operons have been mapped for very few phages, and an annotated global RNA map of a phage – alone or together with its infected host – is not available at all. Here, we applied differential RNA-seq (dRNA-seq) to study the early, host takeover phase of infection by assessing the transcriptome structure of Pseudomonas aeruginosa jumbo phage ɸKZ, a model phage for viral genetics and structural research. This map substantially expands the number of early expressed viral genes, defining TSSs that are active ten minutes after ɸKZ infection. Simultaneously, we record gene expression changes in the host transcriptome during this critical metabolism conversion. In addition to previously reported upregulation of genes associated with amino acid metabolism, we observe strong activation of genes with functions in biofilm formation (cdrAB) and iron storage (bfrB), as well as an activation of the antitoxin ParD. Conversely, ɸKZ infection rapidly down-regulates complexes IV and V of oxidative phosphorylation (atpCDGHF and cyoABCDE). Taken together, our data provide new insights into the transcriptional organization and infection process of the giant bacteriophage ɸKZ and adds a framework for the genome-wide transcriptomic analysis of phage–host interactions.

Genetic circuit characterization by inferring RNA polymerase movement and ribosome usage

To perform their computational function, genetic circuits change states through a symphony of genetic parts that turn regulator expression on and off. Debugging is frustrated by an inability to characterize parts in the context of the circuit and identify the origins of failures. Here, we take snapshots of a large genetic circuit in different states: RNA-seq is used to visualize circuit function as a changing pattern of RNA polymerase (RNAP) flux along the DNA. Together with ribosome profiling, all 54 genetic parts (promoters, ribozymes, RBSs, terminators) are parameterized and used to inform a mathematical model that can predict circuit performance, dynamics, and robustness. The circuit behaves as designed; however, it is riddled with genetic errors, including cryptic sense/antisense promoters and translation, attenuation, incorrect start codons, and a failed gate. While not impacting the expected Boolean logic, they reduce the prediction accuracy and could lead to failures when the parts are used in other designs. Finally, the cellular power (RNAP and ribosome usage) required to maintain a circuit state is calculated. This work demonstrates the use of a small number of measurements to fully parameterize a regulatory circuit and quantify its impact on host.

An RNA-centric view on gut Bacteroidetes

Bacteria employ noncoding RNAs to maintain cellular physiology, adapt global gene expression to fluctuating environments, sense nutrients, coordinate their interaction with companion microbes and host cells, and protect themselves against bacteriophages. While bacterial RNA research has made fundamental contributions to biomedicine and biotechnology, the bulk of our knowledge of RNA biology stems from the study of a handful of aerobic model species. In comparison, RNA research is lagging in many medically relevant obligate anaerobic species, in particular the numerous commensal bacteria comprising our gut microbiota. This review presents a guide to RNA-based regulatory mechanisms in the phylum Bacteroidetes, focusing on the most abundant bacterial genus in the human gut, Bacteroides spp. This includes recent case reports on riboswitches, an mRNA leader, cis- and trans-encoded small RNAs (sRNAs) in Bacteroides spp., and a survey of CRISPR-Cas systems across Bacteroidetes. Recent work from our laboratory now suggests the existence of hundreds of noncoding RNA candidates in Bacteroides thetaiotaomicron, the emerging model organism for functional microbiota research. Based on these collective observations, we predict mechanistic and functional commonalities and differences between Bacteroides sRNAs and those of other model bacteria, and outline open questions and tools needed to boost Bacteroidetes RNA research.