Protocol Citation: Brandon Robin, Sebastien Bontemps-Gallo 2025. Workflow to shed light on gram-negative bacteria transcriptomes by using direct cDNA Oxford Nanopore sequencing: from sample to data collections . protocols.io https://dx.doi.org/10.17504/protocols.io.8epv5rkw5g1b/v1
Manuscript citation:
Robin, B., Baillez, A., Le Guillouzer, S., Lecoeur, C., Sebbane, F., & Bontemps-Gallo, S. (2025). Exploring temperature-dependent transcriptomic adaptations in Yersinia pestis using direct cDNA sequencing by Oxford Nanopore Technologies. Scientific Report. DOI : 10.1038/s41598-025-05662-1
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: August 28, 2024
Last Modified: June 10, 2025
Protocol Integer ID: 106553
Keywords: RNA-Seq, Yersinia pestis, temperature adaptation, Oxford Nanopore Technology, bacterial transcriptomic, negative bacteria transcriptome, effective method for bacterial transcriptomic, new insights into bacterial transcriptional adaptation, oxford nanopore minion sequencer, bacterial transcriptional adaptation, detailed workflow for whole transcriptome rna, using direct cdna oxford nanopore, direct cdna oxford nanopore, whole transcriptome, whole transcriptome rna, other pathogen, pathogen, other pathogens of public health importance, regulated gene, rna, latest ont r10 chemistry, gene expression, temperature, gene, leveraging direct cdna
Funders Acknowledgements:
French National Research Agency - ANR
Grant ID: ANR-21-CE15-0047 RESISTANT
Disclaimer
DISCLAIMER – FOR INFORMATIONAL PURPOSES ONLY; USE AT YOUR OWN RISK
The protocol content here is for informational purposes only and does not constitute legal, medical, clinical, or safety advice, or otherwise; content added to protocols.io is not peer reviewed and may not have undergone a formal approval of any kind. Information presented in this protocol should not substitute for independent professional judgment, advice, diagnosis, or treatment. Any action you take or refrain from taking using or relying upon the information presented here is strictly at your own risk. You agree that neither the Company nor any of the authors, contributors, administrators, or anyone else associated with protocols.io, can be held responsible for your use of the information contained in or linked to this protocol or any of our Sites/Apps and Services.
Abstract
We developed a detailed workflow for whole transcriptome RNA-Seq analysis using the Oxford Nanopore MinION sequencer to profile the pathogen’s gene expression under conditions mimicking its infection cycle. By leveraging direct cDNA sequencing with the latest ONT R10 chemistry, we identified temperature-regulated genes, revealing new insights into bacterial transcriptional adaptations. Our workflow provides a reproducible and cost-effective method for bacterial transcriptomics, contributing to the study of temperature-dependent regulatory mechanisms in Y. pestis and other pathogens of public health importance.
Protocol materials
Maxima H Minus First Strand cDNA Synthesis KitThermo FisherCatalog #K1652
Reduce viscosity and clear the lysate by filtration through NucleoSpin‱ Filter (violet ring) in a collection tube centrifuging 11000 x g, Room temperature, 00:01:00.
Discard the NucleoSpin‱ Filter and add350 µL Ethanol 70% to the homogenized lysate and mix by pipetting up and down (5 times).
Pipet lysate up and down 2–3 times and load the lysate to the NucleoSpin‱ RNA Column (light blue ring) placed in a collection tube.
11000 x g, Room temperature, 00:00:30 and place the column in a new Collection Tube (2 mL).
Add350 µL MDB (Membrane Desalting Buffer) and 11000 x g, Room temperature, 00:01:00. Discard flow through and place the column back into the collection tube.
Prepare DNase reaction mixture in a sterile 1.5 mL microcentrifuge tube: For each isolation, add 10 μL reconstituted rDNase to 90 μL Reaction Buffer for rDNase. Mix by flicking the tube.
Apply95 µL DNase reaction mixture directly onto the center of the silica membrane of the column.
Incubate00:15:00 at Room temperature.
Add200 µL Buffer RAW2 and 11000 x g, Room temperature, 00:00:30.
Place the column into a new Collection Tube (2 mL).
Add600 µL Buffer RA3 and 11000 x g, Room temperature, 00:00:30.
Discard flow through and place the column back into the collection tube.
Add250 µL Buffer RA3 and 11000 x g, Room temperature, 00:00:30.
Discard flow through and place the column back into the collection tube.
11000 x g, Room temperature, 00:02:00 to dry the membrane completely.
Place the column into a 1.5 mL nuclease-free collection tube provided with the kit.
Place the 1.5 mL nuclease-free collection tube and the column with the cap opened at 50 °C for 00:02:00 to dry the membrane completely.
Elute the RNA in 40 µL RNase-free H2O (prewarmed at 50 °C) and incubate at 50 °C for 00:02:00.
11000 x g, Room temperature, 00:01:00.
Recoverand deposit the eluate onto the center of the silica membrane of the column for a second elution and 11000 x g, Room temperature, 00:01:00.
Aliquot3 µL eluate to evaluate the quantity and quality of total RNA and store total RNA at -80 °C.
Analyze1 µL to evaluate the RNA concentration of total RNA .
MICROBExpress‱ is designed to rapidly enrich bacterial mRNA from purified total RNA by removing the 16S and 23S ribosomal RNAs. MICROBExpress is designed so that small RNAs (including tRNA and 5S rRNA) remain in the enriched mRNA population if they were present in the total RNA treated in the procedure.
Pipet200 µL Binding Buffer into a 1.5 mL tube provided with the kit.
Addtotal RNA (2-10 µg in a maximum volume of 15 μL) to the Binding Buffer and vortex gently to mix.
Note
For the proof of concept, we used the maximal quantity supported by the kit, 10 µg
Add4 µL Capture Oligo Mix to the total RNA in the Binding Buffer.
Vortex gently to mix, and microfuge briefly to get the mixture to the bottom of the tube.
Heat denaturate RNA to 70 °C for 00:10:00
Anneal at 37 °C for 00:15:00 minimum to 01:00:00 maximum (it will result in only a very modest increase in rRNA removal).
During annealing incubation, prepare the Oligo MagBeads.
Withdraw 50 µL Oligo MagBeadsper sample to a 1.5 mL tube (max 500 µL per tube, i.e. max for 10 samples per tube).
Note
Record the volume of Oligo MagBeads withdrawn.
Capture the Oligo MagBeads by placing the tube on a magnetic stand. Leave the tube on the stand until all of the Oligo MagBeads are arranged inside the tube near the magnet.
Note
This will take 3 min or more depending the magnet stand used.
Leaving the tube on the magnet stand, carefully remove the supernatant by aspiration and discard the supernatant.
Wash the Oligo MagBeads with an equal volume of Nuclease-free H2O than the recorded volume in step 32.1 by brief and gentle vortexing.
Leaving the tube on the magnet stand, carefully remove the supernatant by aspiration and discard the supernatant.
Equilibrate the Oligo MagBeads with an equal volume of Binding Buffer (recorded volume in step 32.1) by brief and gentle vortexing.
Leaving the tube on the magnet stand, carefully remove the supernatant by aspiration and discard the supernatant.
Resuspend the Oligo MagBeads in an equal volume of Binding Buffer (recorded volume in step 31.1) by gently tapping the tube.
Prewarm the Oligo MagBead slurry at 37 °C as well as the Wash Solution before proceeding.
Add50 µL prewarmed prepared Oligo MagBeads to the RNA/Capture Oligo Mix (incubated at step 32).
Incubate at 37 °C for 00:15:00.
Capture the Oligo MagBeads by placing the tube on the magnet stand.
Carefully aspirate the supernatant to avoid dislodging the Oligo MagBeads.
Transferit to aCollection TubeOn ice
Add100 µL prewarmed Wash Solution to the captured Oligo MagBeads. Resuspend the beads by brief and gentle vortexing.
Recapture the Oligo MagBeads on the magnet stand and carefully recover the supernatant.
Pool this supernatant with the RNA already in the Collection Tube On ice.
Add35 µL 3 M Sodium Acetate, 7 µL 5 mg/mL Glycogen,1200 µL 100% Ethanol to the Collection Tube.
Precipitate at least at-20 °COvernight.
Note
Prepare 70% Ethanol and store at -20 °C for the following washes.
15000 x g, 4°C, 00:30:00
Carefully discard the supernatant by inverting the tube. The pellet could be invisible.
Wash by adding 750 µL cold 70% Ethanol, vortex briefly.
15000 x g, 4°C, 00:15:00 and carefully decant and discard the supernatant by inverting the tube.
Repeat the wash steps 42-43
Short-spin the tube and carefully remove any remaining supernatant with a pipettor avoiding to dislodge the pellet.
Air dry the pellet for 00:05:00 maximum
Resuspend the mRNA-enriched pellet in 18 µL Nuclease-free Water and incubate00:15:00 atRoom temperature. If necessary, vortex the sample vigorously to resuspend the RNA.
Note
If the RNA solution is brown, there is probably a small amount of Oligo Magbeads remaining in the sample. To remove them, put the tube on the magnet stand for 3 min and move the enriched mRNA solution to a new RNase-free tube.
Aliquot2 µL eluate to evaluate the quantity and quality of mRNA-enriched sample and store it at -80 °C.
Quantification obtained with the Qubit 4 Fluorometer: 50-200 ng/µL
RIN obtained with the Bioanalyzer: N/A
The RIN is calculated with the 16S and 23S rRNA peaks. However, the smaller the 16S and 23S rRNA peaks, the more efficient the depletion. The efficiency depends on the bacteria of interest and the starting total RNA quantity.
In a 1.5 ml microcentrifuge tube, set up 3’ polyadenylation reaction as follows:
A
B
Reagent
Volume
Non-polyadenylated mRNA-enriched
X µl
10X E. coli poly(A) polymerase buffer
2 µl
ATP (10 mM)
2 µl
Nuclease-free water
15-X µl
E. coli poly(A) polymerase (5 U/µl)
1 µl
TOTAL
20 µl
Note
The quantity of starting mRNA-enriched depends on the quantity obtained after rRNA depletion. We recommend using more than 300 ng and less than 10 µg. Here, we used 1,5 µg.
Incubate at 37 °C for 00:00:45. The recommendations from Oxford Nanopore are between 0.5 and 1.5 minutes maximum.
Stop the reaction by adding5 µL 50 mM EDTA.
sparQ PureMag BeadsQuantabioCatalog #95196-060
Equilibrate at Room temperature for at least 00:30:00before use and vigorously vortex to resuspend the beads.
Add100 µL SparQ PureMag Beads (0.8X for RNA fragments >150 bp, decreased the ratio until 0.5X to select long fragments or increase up to 0.9X to keep the short fragments).
Incubate on a Hula Mixer 100 rpm, 00:05:00 at Room temperature.
MICROTUBE 0.5 ML PP NEUTRALFisher ScientificCatalog #12655512
Expected result
The expected concentration of 3'-polyadenylated mRNA-enriched as a function of input material: 15-100 ng/µL
Store the samples at -80 °C or proceed immediately with the library preparation, keeping your sample On ice.
5 - Reverse transcription and strand-switching
2h 30m
Note
From this point, we followed the SQK-LSK114 ONT protocol (version DCS_9187_v114_revG_19Apr2023) coupled to the SQK-NBD114.24 protocol (version NBA_9168_v114_revL_15Sep2022) with few adjustments.
It is important to check the newly available protocol update on the Oxford Nanopore Community to take notes about the potential new major recommendations.
Transfer100 ng of polyA mRNA into 0.2 mL PCR tube and adjust the volume to 7.5 µL with nuclease-free H2O.
Mix by flicking the tube to avoid unwanted shearing and spin down briefly.
Add2.5 µL 2 µM VN primer and 1 µL 10 mM dNTPs to the PCR tube.
MICROTUBE 0.5 ML PP NEUTRALFisher ScientificCatalog #12655512
Expected result
The expected concentration of cDNA: 0.5-2 ng/µL
To calculate the cDNA molarity, analyze 1 μl of the strand-switched cDNA length with a Bioanalyzer.
Note
If the mean cDNA length is known, the molarity could be calculated with the NEBioCalculator.
Store the samples at -20 °C for long-term storage or at 4 °C for short-term storage.
7 - End-prep
1h 30m
ONT Native barcoding sequencing kit v14 (24)Oxford Nanopore TechnologiesCatalog #SQk-NBD114.24
The ONT Native barcoding sequencing kit v14 (24) contains 24 barcodes. In other words, the maximum number of samples to be simultaneously run in compatible R10.4.1 flow cells is 24.
Dilute the DNA Control Sample (DCS) by adding 105 µL Elution Buffer (EB)directly to one DCS tube. The final concentration is 10 ng/µL.
Mix gently by pipetting and spin down.
In 0.2 ml PCR tubes, aliquot the same highest quantity possible of cDNA per sample up to 11.5 µL.
Note
Expect a minimum of 70 fmol for all samples summed. To date, we have not go below 5 fmol of cDNA per sample with 24 samples. Here, we used 14 fmol of cDNA per sample.
Make up each sample to 11.5 μl using Nuclease-free H2O. Mix gently by pipetting and spin down.
Calculate the DNA Control Sample volume (X µL) to put in each cDNA sample. To know, 1 µL of DCS should be used for 200 fmol (130 ng for 1 kb cDNA, the maximum quantity per sample) of cDNA sample.
Note
Here, for 14 fmol of cDNA, we used 0.07 µL DCS per sample.
Prepare the following end-prep mix:
A
B
Reagent
Volume per sample
DNA Control Sample
X µL
Ultra II End-prep Reaction Buffer
1.75 µL
Ultra II End-prep Enzyme Mix
0.75 µL
Nuclease-free H2O
1 - X µL
Mix gently by flicking the tube and spin down.
NEBNext Ultra II End Repair/dA-Tailing Module - 96 rxnsNew England BiolabsCatalog #E7546L
Add3.5 µL end-prep mix to the 11.5 µL of cDNA sample. Mix by pipetting and spin down.
Incubate in a thermal cycler 20 °C00:05:00 and 65 °C00:05:00.
Meanwhile, by vortexing, resuspend the room-temperature equilibrated AMPure XP beads (AXP) provided in the ONT kit and add15 µL AXP into as many 1.5 ml Eppendorf DNA LoBind tubes as samples.
Transfer each cDNA sample into the 1.5 ml Eppendorf DNA LoBind tubes containing AXP.
Mix by flicking the tube.
Incubate on a Hula mixer for 00:05:00 at Room temperature.
Spin down the sample and pellet on a magnet stand until the eluate is clear and colorless. Pipet off the supernatant.
Wash the beads with 200 µL fresh 80 % Ethanol without disturbing the pellet.
Discard the ethanol using a pipette.
Repeat the previous step.
Spin down and place the tube back on the magnet stand. Pipet off any residual ethanol.
Allow to dry for 00:00:30, but do not dry the pellet to the point of cracking.
Remove the tube from the magnet stand and resuspend the pellet in 9 µL Nuclease-free H2O.
Incubate on a Hula mixer 100 rpm, 00:02:00 at Room temperature.
Briefly spin down the tube and pellet the beads on the magnet stand until the eluate is clear and colorless, for at least 00:01:00.
Remove and retain9 µL eluate into a 1.5 mL Eppendorf LoBind tube.
Analyze1 µL end-prepped cDNA with a Qubit fluorometer using the Qubit dsDNA HS assay kit.
MICROTUBE 0.5 ML PP NEUTRALFisher ScientificCatalog #12655512
Expected result
The expected concentration of end-prepped cDNA: 1-4 ng/µL1
Take forward the highest equimolar mass possible of each sample to be barcoded up to 7.5 µL into the native barcode ligation step. However, you may store the samples at 4°C overnight.
Note
To date, we fruitfully tested the range 7-30 ng cDNA per sample.
Here, we pooled 15 ng (16 fmol) of each end-prepped cDNA sample.
8 - Native barcode ligation
2h
ONT Native barcoding sequencing kit v14 (24)Oxford Nanopore TechnologiesCatalog #SQk-NBD114.24
Select a unique barcode for each sample to be run together on the same flow cell.
Up to 24 samples can be barcoded and combined in one experiment.
In 0.2 ml PCR tubes, add the reagents in the following order:
A
B
Reagent
Volume
End-prepped DNA
7.5 µL
Native Barcode (NB01-24)
2.5 µL
Blunt/TA Ligase Master Mix
10 µL
Thoroughly mix the reaction by gently pipetting and briefly spinning down
Blunt/TA Ligase Master Mix - 50 rxnsNew England BiolabsCatalog #M0367S
Incubate for 00:20:00 at Room temperature.
The following volume of EDTA depends on the cap color EDTA tube in the ONT kit.
Add2 µL clear cap EDTA or 4 µL blue cap EDTA to each tube, mix thoroughly by pipetting and spin down briefly.
Pool all the barcoded samples in a 1.5 mL Eppendorf DNA LoBind tube.
By vortexing, resuspend the room-temperature equilibrated AMPure XP beads (AXP) provided.
Add AMPure XP Beads (AXP) to the pooled reaction, and mix by pipetting for a 0.4X clean.
In other words: 14.67 µL AXP per sample used if you added the clear cap EDTA
OR 16 µL AXP per sample used if you added the blue cap EDTA
Note
Here, we added blue cap EDTA in 20 barcoded samples, thus we performed a 320 µL AXP clean.
Incubate on a Hula mixer 100 rpm, 00:10:00 at Room temperature.
Spin down the sample and pellet on a magnet stand for 00:05:00 until the eluate is clear and colorless. Pipet off the supernatant.
Keeping the tube on the magnet stand, wash the beads with 700 µL fresh 80% Ethanol without disturbing the pellet. Discard the ethanol using a pipette.
Repeat the previous step.
Spin down and place the tube back on the magnet stand. Pipet off any residual ethanol.
Allow to dry for 00:00:30, but do not dry the pellet to the point of cracking.
Remove the tube from the magnet stand and resuspend the pellet in 33 µL Nuclease-free H2O by gently flicking the tube.
Incubate for 00:10:00 at 37 °C. Agitate the sample for 00:00:10 by gently flicking the tube every 00:02:00 to encourage DNA elution.
Pellet the beads on a magnet stand until the eluate is clear and colorless.
Remove and retain33 µL eluate into a clean 1.5 mL Eppendorf DNA LoBind tube.
Analyze1 µL barcoded cDNA with a Qubit fluorometer using the Qubit dsDNA HS assay kit.
The expected concentrations of adapted DNA library: 3-10 ng/µL13-20 nanomolar (nM)
Make up the library to12 µL containing 35-50 fmol (for 1-10kb library length) or 100 fmol (for <1kb library length).
The prepared library is ready for loading onto the R10.4.1 flow cell. Store the library on ice until loading.
Note
Here, we loaded around 100 fmol.
10 - Priming and loading the SpotON flow cell
45m
On the computer, plug in the MinION, check for the operating system updates and prohibit it during the run. Get the last MinKNOW software version and start a hardware check.
ONT will replace any flow cell with fewer than 800 nanopores when the result is reported within two days of performing the flow cell check, and when the storage recommendations have been followed. The warranty on this product is 12 weeks from receipt by the customer.
Meanwhile, at Room temperature, prepare the flow cell priming mix as directed below:
A
B
Reagent
Volume
Flow Cell Flush (FCF)
1170 µL
50 mg/ml Bovine Serum Albumin (BSA)
5 µL
Flow Cell Tether (FCT)
30 µL
Mix by pipetting.
ONT Native barcoding sequencing kit v14 (24)Oxford Nanopore TechnologiesCatalog #SQk-NBD114.24
Slide the flow cell priming port cover clockwise to open the priming port and draw back a small volume to
remove any bubbles in the following way:
1) Set a P1000 pipette to 200 μl
2) Insert the tip into the priming port
3) Turn the wheel until you can see a small volume of yellow buffer entering the pipette tip
Note
Visually check that there is continuous buffer from the priming port to the beginning of the waste channel after the sensor array to be sure that the array of pores is covered by buffer at all times.
Avoiding the introduction of air bubbles, load800 µL priming mix into the flow cell via the openedpriming port (form a drop at the end of the pipette tip before inserting into the priming port and gently load the priming mix). You can see the yellow buffer enter the waste channel.
Wait for 00:05:00.
During this time, in a 1.5 ml Eppendorf DNA LoBind tube, prepare the library as follows:
A
B
Reagent
Volume
Sequencing Buffer (SB)
37.5 μL
(Pipet mix immediately before use) Library Beads (LIB)
25.5 μL
adapted cDNA library
12 μL
Gently lift the SpotON sample port cover (to make the SpotON sample port accessible and enable the creation of a drop before sample loading).
To complete the flow cell priming, load200 µL priming mix into the flow cell via the openedpriming port (not the SpotON sample port), avoiding the introduction of air bubbles as already mentioned.
Mix the prepared library gently by pipetting and immediately load 75 µL prepared library to the flow cell via the SpotON sample port in a dropwise fashion. Ensure each drop flows into the port before adding the next.
Gently replace the SpotON sample port cover, making sure the bung enters the SpotON port and close the priming port.
Close the device lid and set up a sequencing run on MinKNOW.
Fill in the fields step-by-step (Experiment ID, Sample ID) and select the kit SQK-NBD114.24.
In sequencing analysis settings, basecalling and barcoding should be ON, select Trim barcodes since each barcode data is separately stored, turn the alignment ON and check the output settings. Make sure there is enough free disk space (at least 500 GB) and use POD5 file extension to rebasecall. The MinQ Score should be 8. To determine the sequencing duration, choose 72 h or less, or a number of reads depending on the objectives.
Note
We recommend downloading the sequence FASTA file of your strain to benefit from the real-time analysis in terms of sequence coverage. In this way, you can stop the run when the required sequencing depth is obtained (10X minimum for transcriptomes).