Abstract:ABSTRACT RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6–1.7 kbp, while the median size of the 5′- and 3′-untranslated regions (UTRs) were 30–90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100–200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community. IMPORTANCE Our understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Deciphering Bacterial and Archaeal Transcriptional Dark Matter and Its Architectural Complexity

Deciphering transcript architectural complexity in bacteria and archaea

Discovering putative peptides encoded from non-coding RNAs in ribosome profiling data of Arabidopsis thaliana.

Bacterial transcriptomics: what is beyond the RNA horiz-ome?

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

A Novel Computational Framework for Genome-Scale Alternative Transcription Units Prediction.

Robust identification of noncoding RNA from transcriptomes requires phylogenetically-informed sampling

Deep Learning Prediction of Ribosome Profiling with Translatomer Reveals Translational Regulation and Interprets Disease Variants

Ab initio construction of a eukaryotic transcriptome by massively parallel mRNA sequencing

Srnascanner: a Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes.

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

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

Illuminating the dark side of the human transcriptome with long read transcript sequencing

Landscape of the Dark Transcriptome Revealed Through Re-mining Massive RNA-Seq Data.

Condition-Dependent Transcriptome Reveals High-Level Regulatory Architecture in Bacillus subtilis

An Automated Approach for Global Identification of Srna-Encoding Regions in RNA-Seq Data from Mycobacterium Tuberculosis.

Systematic Analysis of Transcribed Loci in ENCODE Regions Using RACE Sequencing Reveals Extensive Transcription in the Human Genome

Comparative Analysis of the Transcriptome Across Distant Species.

Landscape of transcription in human cells

De novo computational prediction of non-coding RNA genes in prokaryotic genomes.