Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi

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.

PIPETS: A statistically informed, gene-annotation agnostic analysis method to study bacterial termination using 3'-end sequencing

Background

Over the last decade the drop in short-read sequencing costs have allowed experimental techniques utilizing sequencing to address specific biological questions to proliferate, oftentimes outpacing standardized or effective analysis approaches for the data generated. There are growing amounts of bacterial 3'-end sequencing data, yet there is currently no commonly accepted analysis methodology for this datatype. Most data analysis approaches are somewhat ad hoc and, despite the presence of substantial signal within annotated genes, focus on genomic regions outside the annotated genes (e.g. 3' or 5' UTRs). Furthermore, lack of a systematic approach to analyzing such data makes it impossible to compare conclusions generated by different labs, using different organisms.

Results

We present PIPETS (Poisson Identification of PEaks from Term-Seq data), an R package available on Bioconductor that provides a novel analysis method for 3'-end sequencing data. PIPETS is a statistically informed, gene-annotation agnostic methodology. Across two different datasets PIPETS identified significant 3'-end termination signal across a wider range of annotated genomic contexts than existing analysis approaches, suggesting that existing approaches may miss biologically relevant signal. Furthermore, assessment of the previously called 3'-end positions not captured by PIPETS showed that they were uniformly very low coverage.

Conclusions

PIPETS provides a broadly applicable platform to compare 3'-end sequencing data sets across different organisms. It requires only the 3'-end sequencing data, and is broadly accessible to non-expert users.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Broadly conserved FlgV controls flagellar assembly and Borrelia burgdorferi dissemination in mice

Flagella propel pathogens through their environments yet are expensive to synthesize and are immunogenic. Thus, complex hierarchical regulatory networks control flagellar gene expression. Spirochetes are highly motile bacteria, but peculiarly in the Lyme spirochete Borrelia burgdorferi, the archetypal flagellar regulator σ28 is absent. We rediscovered gene bb0268 in B. burgdorferi as flgV, a broadly-conserved gene in the flagellar superoperon alongside σ28 in many Spirochaetes, Firmicutes and other phyla, with distant homologs in Epsilonproteobacteria. We found that B. burgdorferi FlgV is localized within flagellar motors. B. burgdorferi lacking flgV construct fewer and shorter flagellar filaments and are defective in cell division and motility. During the enzootic cycle, B. burgdorferi lacking flgV survive and replicate in Ixodes ticks but are attenuated for dissemination and infection in mice. Our work defines infection timepoints when spirochete motility is most crucial and implicates FlgV as a broadly distributed structural flagellar component that modulates flagellar assembly.

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.