Nanopore Bacterial Direct RNA Sequencing
Bacterial transcriptomes pose unique challenges that standard RNA sequencing methods were not designed to solve. Bacterial mRNAs lack poly(A) tails, carry abundant ribosomal RNA contamination (95–98% of total RNA), and are organized into polycistronic operons whose boundaries and internal regulatory features are masked by short-read fragmentation. CD Genomics offers Nanopore Bacterial Direct RNA Sequencing (DRS) — an Oxford Nanopore Technologies (ONT)-based approach optimized for bacterial RNA that directly sequences native RNA molecules without reverse transcription or PCR amplification, capturing full-length bacterial transcripts, operon structures, and single-base RNA modifications in a single experiment.
Unlike conventional bacterial RNA-seq that relies on fragmentation and computational transcript reconstruction, our bacterial DRS workflow employs bacterial-specific rRNA depletion followed by polyadenylation and 3′ adapter ligation, enabling direct detection of complete transcripts from 5′ to 3′ end. The platform simultaneously delivers transcriptome sequence, operon architecture, and base-resolution RNA modification profiling (m6A, m5C, pseudouridine Ψ, inosine, and 2′-O-methylation) — providing an integrated view of bacterial gene expression and epitranscriptomic regulation without the biases introduced by amplification or cDNA conversion.
Why Choose Nanopore Bacterial Direct RNA Sequencing?
- Full-length bacterial transcript resolution — Reads span complete transcripts including 5′-UTR, CDS, intergenic regions, and 3′-UTR without fragmentation or assembly
- Direct RNA modification detection — Identifies m6A, m5C, Ψ, inosine, and 2′-O-methylation at single-base resolution without antibodies or chemical treatment
- Operon structure mapping — Long reads capture polycistronic transcript architectures, revealing internal promoters, terminators, and co-transcription patterns
- Bacterial-optimized workflow — rRNA depletion and polyadenylation protocols specifically designed for bacterial RNA (no poly-A tail requirement)
Introduction: The Bacterial Transcriptome Challenge
Bacterial gene expression is fundamentally different from eukaryotic transcription. Genes are organized into polycistronic operons, transcripts lack the poly(A) tails that eukaryotic mRNAs carry, and regulatory complexity is encoded at the level of transcript boundaries, internal promoters, and RNA modifications that modulate stability and translation efficiency. Short-read RNA sequencing — the current standard for bacterial transcriptomics — fragments RNA into 50–300 bp segments, destroying the connectivity between proximal genes within an operon and forcing computational reconstruction that frequently misassigns or entirely misses complex transcript architectures.
The biological consequences are significant. Operon boundaries defined by short-read data often miss internal transcription start sites, attenuator structures, and condition-dependent alternative transcripts. RNA modifications — particularly m6A and pseudouridine — are increasingly recognized as key regulators of bacterial RNA stability, virulence gene expression, and stress adaptation, but detecting these modifications requires either multiple separate assays (MeRIP-seq, bisulfite sequencing, etc.) or a direct RNA sequencing approach that preserves the native chemical state of every base. Our Nanopore Bacterial Direct RNA Sequencing service is designed specifically to address these requirements using a bacterial-optimized ONT direct RNA workflow.
What is Bacterial Direct RNA Sequencing?
Bacterial Direct RNA Sequencing is an Oxford Nanopore Technology-based method that sequences native bacterial RNA molecules directly — without chemical conversion, antibody enrichment, or PCR amplification. The core adaptation for bacterial samples lies in the library preparation: total bacterial RNA undergoes species-specific rRNA depletion to remove the abundant ribosomal RNA fraction (which can constitute >95% of bacterial total RNA), followed by polyadenylation to add synthetic poly(A) tails that enable adapter ligation (since native bacterial mRNAs lack poly(A) tails), 3′ adapter ligation, reverse transcription, and motor protein loading for nanopore sequencing.
This is fundamentally different from standard bacterial RNA-seq approaches. Short-read methods (Illumina) fragment the RNA and rely on computational assembly to reconstruct transcripts, losing single-molecule connectivity. cDNA-based long-read methods (PacBio Iso-Seq, ONT cDNA-seq) require reverse transcription and PCR amplification, which lose the native RNA modification signal. Bacterial DRS preserves the original RNA molecule and captures both sequence and modification information simultaneously by measuring the ionic current changes as each nucleotide passes through the nanopore — modified bases produce characteristic signal deviations from their unmodified counterparts, enabling dual-layer readout without additional experimental steps.
The resulting data supports operon structure mapping (including polycistronic transcript boundaries, internal promoters, and terminator sites), full-length transcript isoform identification, quantitative expression profiling, and single-base RNA modification calling — all from a single library preparation and sequencing run.
Key Advantages of Bacterial Direct RNA Sequencing
Scientific Advantages
- Single-base RNA modification profiling on native bacterial RNA
Direct nanopore sequencing reads the chemical state of every base as it exists in the bacterial sample. Modifications including m6A, m5C, pseudouridine (Ψ), inosine, and 2′-O-methylation produce characteristic ionic current signatures detectable at single-nucleotide resolution without antibody enrichment, bisulfite conversion, or chemical labeling. This is particularly valuable for bacterial epitranscriptomics, where modifications function as rapid regulatory switches in stress adaptation and virulence gene expression — dynamics invisible to cDNA-based sequencing methods.
- Operon-resolved transcript discovery and architecture mapping
Long native RNA reads spanning full operon structures eliminate the need for computational transcript assembly. Each read represents a complete transcript molecule, enabling direct detection of polycistronic operon boundaries, internal transcription start sites, attenuator structures, and condition-dependent alternative transcripts. Published data from E. coli bacterial DRS studies have extended the boundaries of over 225 operons beyond their current annotations (Tan et al., Nucleic Acids Research, 2024), demonstrating the resolution advantage over short-read approaches.
Business & Project Advantages
- Multiple data types from a single bacterial RNA sequencing run
Our bacterial DRS generates transcriptome sequence, operon architecture maps, and epitranscriptome modification data concurrently. Researchers receive expression quantifications, full-length transcript annotations, operon boundary calls, and base-resolution modification profiles from one library preparation — reducing both project cost and turnaround complexity compared to running separate RNA-seq, MeRIP-seq, and operon mapping experiments. For researchers seeking transcript-level expression data, our Full-Length Transcriptome Profiling service provides complementary analysis options.
- Customized bioinformatics with bacterial-specific pipelines
We provide end-to-end bioinformatics support aligned with bacterial transcriptome analysis needs. Standard deliverables include basecalled FASTQ files, genome-aligned BAM files with modification probability tracks, transcript-level quantification matrices, and operon structure annotation reports. Advanced analysis — including differential modification analysis, condition-specific operon boundary detection, and multi-strain comparative transcriptomics — is available through our Long-Read Sequencing Data Analysis pipeline.
Applications of Nanopore Bacterial Direct RNA Sequencing
Operon Structure and Transcript Architecture Mapping
- Full-length resolution of polycistronic operon boundaries, including leader regions, internal promoters, and terminator sites
- Identification of condition-dependent alternative transcripts and intra-operonic transcription start sites
- Cross-strain and cross-species comparative analysis of operon conservation and divergence
Bacterial Epitranscriptomics and RNA Modification Biology
- Single-base m6A mapping across the bacterial transcriptome to identify modification-regulated transcripts in stress response, virulence, and metabolic regulation
- Co-occurrence analysis of multiple modification types on individual bacterial RNA molecules
- Dynamic RNA modification profiling across growth phases, stress conditions (heat shock, antibiotic exposure, biofilm formation), or host infection time courses
Antibiotic Resistance and Pathogenesis Research
- Detection of modification-mediated regulation of antibiotic resistance gene expression at the operon level
- Identification of resistance-associated transcript isoforms and fusion transcripts in clinical isolates
- Integration with Long-Read Sequencing of RNA Methylation and Antibiotic Resistance Gene (ARG) Analysis for multi-faceted epitranscriptomic characterization of resistant strains
Microbial Community and Host-Microbe Interaction Transcriptomics
- Direct RNA sequencing of bacterial pathogens during host infection — capturing in situ transcript architectures and modification states
- Transcript-level resolution in bacterial co-culture and defined community models
- Multi-kingdom RNA analysis for host-pathogen interaction studies using Microbial Genomics with Long-Read Sequencing
Technology Overview — How Bacterial Direct RNA Sequencing Works
1. Total Bacterial RNA Extraction and Quality Assessment
High-integrity total RNA is extracted from bacterial cultures, clinical isolates, or enriched microbial samples using optimized lysis and RNA purification protocols. RNA integrity and purity are assessed by microfluidic electrophoresis. DNase treatment is performed to eliminate genomic DNA contamination, which is particularly critical for bacterial samples due to the high gDNA-to-RNA ratio.
2. Bacterial-Specific rRNA Depletion
Ribosomal RNA (16S, 23S, 5S) is removed using species-matched or closely related probe-based depletion systems. This step is essential because bacterial rRNA constitutes 95–98% of total RNA — far higher than eukaryotic total RNA. Effective rRNA depletion is the single most important factor determining bacterial DRS data quality and mRNA-mapping read yield.
3. Polyadenylation and Adapter Ligation
Since native bacterial mRNAs lack poly(A) tails, a controlled polyadenylation step adds synthetic poly(A) tails using E. coli Poly(A) Polymerase. This is followed by 3′ adapter ligation, which attaches the ONT-compatible adapter carrying the motor protein docking site to the polyadenylated RNA 3′ end. This two-step process ensures that all RNA molecules — coding and non-coding — are captured for sequencing.
Figure 1. Bacterial DRS library preparation workflow: starting from total bacterial RNA through rRNA depletion, polyadenylation, adapter ligation, reverse transcription, and nanopore sequencing with real-time basecalling and modification detection.
4. Reverse Transcription and Motor Protein Loading
Reverse transcription is performed at optimized temperatures to generate an RNA-cDNA hybrid while minimizing secondary structure stalling common in structured bacterial RNAs (rRNA remnants, tRNA, ncRNA). The motor protein is then loaded onto the hybrid, preparing the library for nanopore sequencing.
5. Nanopore Direct RNA Sequencing and Signal Analysis
The RNA-cDNA hybrid is loaded onto an ONT flow cell. The motor protein unwinds the duplex and processively feeds the RNA strand through the nanopore. Real-time ionic current measurements are basecalled using ONT’s Dorado neural-network basecallers with modified-base-aware models. Signal-level analysis (Nanocompore, m6Anet, Tombo) compares per-read ionic current profiles against reference signals for single-base modification probability estimation, while full-length reads are simultaneously aligned to the bacterial reference genome for operon mapping and transcript quantification.
Bioinformatics Analysis
| Analysis Feature | Basic | Advanced |
| Dorado super-accuracy basecalling | ✓ | ✓ |
| Read QC, filtering, adapter trimming | ✓ | ✓ |
| Bacterial genome alignment (minimap2) | ✓ | ✓ |
| Full-length transcript detection and quantification | ✓ | ✓ |
| Operon boundary and polycistronic structure annotation | ✓ | ✓ |
| Single-base RNA modification calling (m6A, m5C, Ψ) | — | ✓ |
| Condition-dependent operon structure comparison | — | ✓ |
| Differential modification analysis across conditions | — | ✓ |
| Custom visualization and publication-ready figures | — | ✓ |
For detailed bioinformatics support options, see our Long-Read Sequencing Data Analysis Services.
Choosing the Right Bacterial RNA Sequencing Approach
The table below compares our three ONT-based RNA sequencing services specifically for bacterial transcriptome analysis to guide researchers in selecting the best approach for their experimental goals.
| Feature | Bacterial Direct RNA Sequencing | Standard Direct RNA Sequencing | Full-Length cDNA Sequencing |
| RNA capture method | rRNA depletion + polyadenylation + 3′ adapter | Poly(A) capture (eukaryotic) | rRNA depletion + random priming |
| Optimized for bacterial RNA | ✔ Yes — bacterial-specific rRNA depletion + polyadenylation | ✘ No — requires poly(A) tails | ✔ Partial — rRNA depletion available |
| Native RNA modifications | ✔ Direct detection preserved | ✔ Direct detection preserved | ✘ Lost during RT-PCR |
| Operon structure resolution | ✔ Full-length polycistronic reads | ✘ Bacterial mRNAs lack poly(A) tails | ✔ Full-length reads with cDNA bias |
| Throughput | Moderate | Low–Moderate | High |
| Best for | Comprehensive bacterial transcriptome + epitranscriptome + operon analysis | Eukaryotic transcriptome modification analysis | High-throughput bacterial isoform discovery |
Sample Requirements for Bacterial Direct RNA Sequencing
| Category | Requirement | Notes |
| Sample type | Total RNA from bacterial culture, clinical isolate, or enriched sample | Fresh or RNAlater-stabilized; DNase-treated recommended |
| Minimum input | ≥2 μg (standard); ≥1 μg (high-yield optimized workflow) | Bacterial RNA yield varies by species and growth conditions |
| RNA integrity | RIN ≥ 7 (standard); degraded RNA assessed case-by-case | Bacterial RNA degrades rapidly — flash-frozen pellets recommended over extracted RNA for long-term storage |
| rRNA depletion | Species-matched probes required | Standard probes: E. coli, S. aureus, P. aeruginosa, B. subtilis, M. tuberculosis; custom probe design available for non-standard species |
| Recommended depth | 5–15 million reads per sample | Higher depth recommended for low-abundance transcript detection and modification calling |
| Growth condition documentation | Detailed growth phase, medium, and treatment conditions recommended | Operon structure and modification profiles are condition-dependent |
Please refer to our Sample Submission Guidelines for detailed instructions on bacterial RNA preparation and shipping.
Why Choose CD Genomics for Bacterial Direct RNA Sequencing
Bacterial transcriptomics expertise across diverse species
CD Genomics has extensive experience processing bacterial RNA from a wide range of species — from Gram-negative model organisms (E. coli, P. aeruginosa) to Gram-positive pathogens (S. aureus, B. subtilis), mycobacteria, and fastidious clinical isolates. This experience informs protocol optimization for each species’ specific rRNA depletion requirements and RNA extraction characteristics.
Optimized bacterial DRS workflow with RNA004 chemistry
Our bacterial DRS protocol has been refined through iterative optimization of rRNA depletion conditions, polyadenylation parameters, and adapter ligation efficiency specifically for bacterial RNA. We deploy the latest ONT RNA004 kit chemistry, which has demonstrated significantly improved sequencing yield for bacterial direct RNA applications compared to earlier chemistries.
End-to-end project support from experimental design to publication
We manage every stage of your bacterial DRS project: initial feasibility assessment (including difficult-to-culture or low-biomass samples), bacterial RNA extraction and QC, rRNA depletion optimization, ONT sequencing on PromethION instruments, and a comprehensive bioinformatics pipeline that delivers operon annotations, modification probability tracks, and expression matrices ready for downstream analysis.
Integrated operon mapping and epitranscriptomics in a single workflow
Unlike providers who separate transcript sequencing from modification analysis, our bacterial DRS service delivers both data types from the same library preparation. This integration is particularly valuable for bacterial systems, where operon structure and RNA modifications are functionally interconnected — modification status can influence operon transcript stability, and operon context determines which transcripts carry specific modifications.
Case Study: Analysis of Bacterial Transcriptome and Epitranscriptome Using Nanopore Direct RNA Sequencing
Tan L, Guo Z, Shao Y, Ye L, Wang M, Deng X, Chen S, Li R. Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing. Nucleic Acids Res. 2024;52(15):8746–8762. DOI: 10.1093/nar/gkae601. (CC BY 4.0)
1. Background
Short-read RNA sequencing of bacterial transcriptomes cannot resolve operon boundaries or identify full-length transcript architectures, and standard eukaryotic DRS protocols are not optimized for bacterial RNA (no poly-A tails, >95% rRNA content). The authors developed a bacterial-specific DRS pre-processing pipeline and evaluated its performance in E. coli and S. aureus, with the goal of establishing a single-assay method for simultaneous transcriptome reconstruction and epitranscriptome analysis.
2. Methods
Total RNA from E. coli K-12 and S. aureus was rRNA-depleted using species-specific probes, polyadenylated with E. coli Poly(A) Polymerase, and processed for ONT direct RNA sequencing on R9.4.1 flowcells. In vitro transcribed (IVT) unmodified RNA libraries were generated from the same bacterial strains as negative controls for modification detection. Data were basecalled with Guppy, aligned with minimap2, and used for operon boundary detection with a custom pipeline. Modification calling was performed using Nanocompore and Tombo, comparing native DRS signals to IVT controls.
3. Results
Figure 2. Operon architecture reconstruction and RNA modification detection in bacterial transcriptomes. From Tan L, Guo Z, Shao Y, et al. Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing. Nucleic Acids Res. 2024;52(15):8746-8762. (CC BY 4.0)
Key Findings
- 225 E. coli operon boundaries extended beyond existing annotations through full-length DRS reads
- 89 S. aureus operon boundaries extended — demonstrating cross-species applicability
- 75 high-confidence m6A candidates identified in E. coli protein-coding transcripts
- Method enables single-assay integration of transcriptome mapping and modification detection
- Bacterial-specific pre-processing (rRNA depletion + polyadenylation) is essential for DRS success in bacterial samples
4. Conclusions
This study established that Nanopore direct RNA sequencing, with appropriate bacterial-specific pre-processing modifications, can effectively characterize both the transcriptome and epitranscriptome of bacterial species in a single experiment. The data demonstrate that bacterial DRS provides operon-level resolution unattainable by short-read methods while simultaneously detecting RNA modifications that would require separate, dedicated assays.
FAQs
Q: Can Bacterial Direct RNA Sequencing work on any bacterial species?
Yes, with appropriate optimization. The key requirement is species-matched or closely related rRNA depletion probes. We maintain a panel of standard probes for common laboratory and pathogenic species (E. coli, S. aureus, P. aeruginosa, B. subtilis, M. tuberculosis) and can design custom probes for non-standard species. RNA extraction protocols may also require optimization for bacteria with thick cell walls or high RNase content.
Q: How does the polyadenylation step affect the data?
The polyadenylation step uses E. coli Poly(A) Polymerase to add synthetic poly(A) tails to all RNA molecules, enabling uniform adapter ligation. This controlled addition does not interfere with downstream sequence alignment or transcript quantification. Importantly, it does not affect modification detection — the native modifications on the original RNA sequence are preserved and detected through the nanopore signal analysis.
Q: What RNA modifications can be detected in bacterial samples?
Our standard analysis detects m6A across the bacterial transcriptome using Nanocompore comparative analysis against unmodified IVT controls (recommended for bacterial samples to distinguish true modifications from sequence-specific signal variation). Extended analysis can additionally profile m5C, pseudouridine (Ψ), inosine, and 2′-O-methylation signal patterns upon request. For bacterial samples, we strongly recommend including matched unmodified controls (IVT libraries) for reliable modification calling.
Q: What is the minimum RNA input required?
The standard protocol requires ≥2 μg of high-quality total bacterial RNA (RIN ≥ 7). A high-yield optimized workflow is available for samples with ≥1 μg. For low-biomass samples (clinical isolates, difficult-to-culture species) or degraded RNA, we recommend consulting our project scientists for feasibility assessment.
Q: What bioinformatics deliverables are provided for operon analysis?
Standard deliverables include full-length transcript annotations with operon boundary calls, transcript-level quantification tables, and genome-aligned BAM files for visualization in IGV. Advanced deliverables add condition-specific operon structure comparisons, internal promoter and terminator site predictions, differential modification analysis, and custom visualization. We support common bacterial reference genomes and can work with user-provided genome assemblies.
Demo Data and Deliverable Examples
Deliverable 1 — Full-length bacterial transcript and operon annotation
Genome browser tracks showing complete bacterial transcript structures identified from DRS reads, including polycistronic operon organization with annotated CDS features, intergenic regions, and 5′/3′-UTR boundaries. Each read represents a full native transcript spanning from transcription start site to terminator.
Deliverable 2 — Single-base m6A modification probability profiles
Per-nucleotide modification probability tracks along bacterial transcripts, displayed as bigWig signal tracks. Modification probabilities are derived from Nanocompore comparative analysis against matched IVT unmodified controls, with confidence scores and genomic coordinates for each candidate modification site.
Deliverable 3 — Condition-dependent operon structure comparison
Comparative visualization of operon architecture under control vs. treatment conditions, highlighting condition-specific transcript boundary shifts, internal promoter activation, and differential transcript abundance within operons.
Deliverable 4 — Transcript-level expression quantification matrix
Count and TPM matrix across all bacterial samples, with transcript annotations including operon membership, gene locus tags, product descriptions, and COG functional categories.
References
- Tan L, Guo Z, Shao Y, Ye L, Wang M, Deng X, Chen S, Li R. Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing. Nucleic Acids Res. 2024;52(15):8746-8762. DOI: 10.1093/nar/gkae601. (CC BY 4.0)
- Riquelme-Barrios S, Vásquez-Camus L, Cusack SA, Burdack K, Petrov DP, Yeşiltaç-Tosun GN, Kaiser S, Giehr P, Jung K. Direct RNA sequencing of the Escherichia coli epitranscriptome uncovers alterations under heat stress. Nucleic Acids Res. 2025;53(6):gkaf175. DOI: 10.1093/nar/gkaf175. (CC BY 4.0)
- Guo Z, Shao Y, Tan L, Lu B, Deng X, Chen S, Li R. Enhanced detection of RNA modifications in Escherichia coli utilizing direct RNA sequencing. Cell Rep Methods. 2025;5(9):101168. DOI: 10.1016/j.crmeth.2025.101168. (CC BY 4.0)