(A) Schematic of classification of Term-seq 3´ ends and Rho termination sites relative to an annotated ORF. 3´ ends and termination sites were defined as: primary (purple colored, located on the …
Schematic of total RNA-seq (modified from the RNAtag-seq methodology, Shishkin et al., 2015), Term-seq (modified from Dar et al., 2016) and DirectRNA-seq (Ozsolak and Milos, 2011). RNA 3´ end …
(A) Principal component analysis plot to show correlation among total RNA-seq and Term-seq replicates. Read counts at annotated ORFs for total RNA-seq or the union of 3´ ends characterized from all …
(A) RNA-seq screenshot of the sugE (gdx) locus displaying sequencing reads from LB 0.4 total RNA-seq, LB 0.4 Term-seq and DirectRNA-seq ±BCM treatment. Total and Term-seq tracks represent an overlay …
(A) RNA-seq screenshot of the thiM locus displaying sequencing reads from LB 0.4 total RNA-seq, LB 0.4 Term-seq and DirectRNA-seq ± BCM treatment. Total and Term-seq tracks represent an overlay of …
(A) Sequence of the mdtJI 5´ UTR. The transcription start site (green shaded nucleotide) determined by dRNA-seq (Thomason et al., 2015) and 3´ end (red shaded nucleotide) determined by Term-seq …
Unannotated uORF (mdtU) sequences were selected by searching for short (20–40 amino acid) ORFs with <300 nt distance from the start of an mdtJ ortholog in 1742 sequenced gammaproteobacterial …
(A) Sequence of documented region of sRNA-mRNA pairing. 3´ end determined by Term-seq is highlighted in red. Start codon of the corresponding ORF is indicated with green text. (B) Northern analysis …
Transcription start sites (nucleotide shaded green) determined by dRNA-seq (Thomason et al., 2015) and 3´ ends (nucleotide shaded red) determined by Term-seq (current study) are indicated. Start …
(A) RNA-seq screenshot of the ChiZ and IspZ loci displaying sequencing reads from the LB 0.4 growth condition from dRNA-seq (Thomason et al., 2015, HS2 samples), total RNA-seq and Term-seq. Total …
Secondary structures of ChiX (A) and IspZ (B) predicted by sfold (http://sfold.wadsworth.org/cgi-bin/srna.pl; Ding et al., 2004). The base-pairing regions are highlighted in yellow, and predicted …
(A) Northern analysis for ChiZ in WT (GSO989) and rhoR66S mutant (GSO990) cells. Cells were grown to OD600 ~0.4 or 2.0 after a dilution of the overnight culture. Total RNA was extracted, separated …
(A) RNA-seq screenshots of the ftsI, aceK, rlmD, mglC, and ampG mRNAs containing putative internal (int) sRNAs. Sequencing reads from the LB 0.4 dRNA-seq (Thomason et al., 2015, HS2 samples), total …
(A) RIL-seq screenshot showing RybB chimeras at the ftsO locus. Data are from Hfq-FLAG LB RIL-seq performed 150 min after a dilution of the overnight culture, (Melamed et al., 2020, RIL-seq …
(A) Northern analysis showing reciprocal effect of RybB on FtsO. RNA was extracted from WT (GSO982) cells harboring either the pBR or pBR-RybB at 150 and 360 min after dilution of the overnight …
(A) Histogram of percent nucleotide sequence identity at wobble positions across the ftsI gene. Clustal (Madeira et al., 2019) was used to align the following species against the E. coli K12 MG1655 …
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Strain, strain background (Escherichia coli) | MG1655 (WT) | this study | N/A | see Supplementary file 4 for derivatives |
Antibody | mouse monoclonal anti-FLAG-M2-HRP | Sigma-Aldrich | Cat#A8592 | WB (1:2000) |
Antibody | rabbit polyclonal anti-OmpC | Biorbyt | Cat#orb6940 | WB (1:500) |
Antibody | donkey polyclonal peroxidase labeled anti-rabbit | GE Healthcare | Cat#NIF824 | WB (1:1000) |
Recombinant DNA reagent | pBRplac | Guillier and Gottesman, 2006 | N/A | see Supplementary file 4 for derivatives |
Recombinant DNA reagent | pNM46 (pBRplac-lacI) | this study | N/A | see Supplementary file 4 for derivatives |
Recombinant DNA reagent | pKD46 | Datsenko and Wanner, 2000 | N/A | see Supplementary file 4 for derivatives |
Recombinant DNA reagent | pMM1 | Stringer et al., 2014 | N/A | see Supplementary file 4 for derivatives |
Sequence-based reagent | northern probes and primers | this study | N/A | see Supplementary file 4 |
Chemical compound, drug | spermidine | Sigma-Aldrich | Cat#S2626-1G | |
Other | TRIzol reagent | Thermo Fisher Scientific | Cat#15596018 | RNA extractions |
Other | RNA-sequencing reagents | Melamed et al., 2018 | N/A | |
Other | ureagel-8 | National Diagnostics | Cat#EC-838 | acrylamide northern solution |
Other | ureagel complete | National Diagnostics | Cat#EC-841 | acrylamide northern solution |
Other | NuSieve 3:1 agarose | Lonza | Cat#50090 | agarose for northern blotting |
Other | 37% formaldehyde | Fisher Scientific | Cat#BP531-500 | |
Other | RiboRuler high range RNA ladder | Thermo Fisher Scientific | Cat#SM1821 | |
Other | RiboRuler low range RNA ladder | Thermo Fisher Scientific | Cat#SM1831 | |
Other | Zeta-Probe blotting membrane | Bio-Rad | Cat#1620159 | northern membrane |
Other | ULTRAhyb-oligo hybridization buffer | Thermo Fisher Scientific | Cat#AM8663 | |
Other | γ-32P ATP | PerkinElmer | Cat#NEG035C010MC | |
Other | T4 polynucleotide kinase | New England Biolabs | Cat#M0201L | |
Other | Illustra MicroSpin G-50 columns | GE Healthcare | Cat#27533001 | |
Other | mini-PROTEAN TGX gels | Bio-Rad | Cat#456–1086 | |
Other | EZ-RUN pre-stained Rec protein ladder | Fisher Scientific | Cat#BP3603-500 | |
Other | nitrocellulose membrane | Thermo Fisher Scientific | Cat#LC2000 | western membrane |
Other | SuperSignal West Pico PLUS chemiluminescent substrate | Thermo Fisher Scientific | Cat#34580 | |
Software, algorithm | lcdb-wf | ‘lcdb-wf’. Dale et al., GitHub Repository | v1.5.3 | https://github.com/lcdb/lcdb-wf |
Software, algorithm | sra-tools | SRA Toolkit Development Team | v2.9.1_1 | http://ncbi.github.io/sra-tools/ |
Software, algorithm | cutadapt | Martin, 2011 | v2.3 | https://github.com/marcelm/cutadapt |
Software, algorithm | fastqc | Wingett and Andrews, 2018 | v0.11.8 | https://qubeshub.org/resources/fastqc |
Software, algorithm | bwa | Li and Durbin, 2010 | v0.7.17 | http://bio-bwa.sourceforge.net/ |
Software, algorithm | bowtie2 | Langmead and Salzberg, 2012 | v2.3.5 | http://bowtie-bio.sourceforge.net/bowtie2 |
Software, algorithm | samtools | Li et al., 2009 | v1.9 | http://www.htslib.org/ |
Software, algorithm | subread | Liao et al., 2013 | v1.6.4 | http://subread.sourceforge.net/ |
Software, algorithm | multiqc | Ewels et al., 2016 | v1.7 | https://github.com/ewels/MultiQC |
Software, algorithm | picard | ‘Picard Toolkit.’ 2019. Broad Institute, GitHub Repository | v2.20.0 | https://github.com/broadinstitute/picard |
Software, algorithm | deeptools | Ramírez et al., 2016 | v3.2.1 | https://deeptools.readthedocs.io/ |
Software, algorithm | termseq-peaks | ‘termseq-peaks.’ 2020. NICHD-BSPC, GitHub Repository | N/A | https://github.com/NICHD-BSPC/termseq-peaks |
Software, algorithm | bedtools | Quinlan and Hall, 2010 | v2.27.1 | https://bedtools.readthedocs.io/ |
Software, algorithm | pybedtools | Dale et al., 2011 | v0.8.0 | https://github.com/daler/pybedtools |
Software, algorithm | ucsc-toolkit | Kent et al., 2010 | v357 | https://doi.org/10.1093/bioinformatics/btq351 |
Software, algorithm | biopython SeqIO | Cock et al., 2009 | v1.73 | biopython.org |
Software, algorithm | rhoterm-peaks | ‘rhoterm-peaks’ 2020. GitHub repository | N/A | https://github.com/gbaniulyte/rhoterm-peaks |
Software, algorithm | CLG Genomics Workbench | Qiagen | v8.5.1 | alignment of DirectRNA-seq reads |
Software, algorithm | CLUSTAL | Madeira et al., 2019 | v2.1 | https://www.ebi.ac.uk/Tools/msa/clustalo/ |
Software, algorithm | BOXSHADE | v3.21 | https://embnet.vital-it.ch/software/BOX_form.html | |
Software, algorithm | Sfold | Ding et al., 2004 | v2.2 | http://sfold.wadsworth.org/cgi-bin/srna.pl |
3´ ends identified by Term-seq.
Curated 3´ ends for E. coli MG1655 (WT) grown to OD600 ~0.4 and 2.0 in LB, and to OD600 ~0.4 in minimal (M63) glucose medium. The genomic coordinate of the 3´ end (3´ end position), DNA strand, average RNA-seq read count of the 3´ end of both biological replicates, 3´ end classification (see Materials and methods), the gene annotation for the classification (details), and the sequence surrounding the 3´ end (40 bp upstream and 10 bp downstream, 3´ end nucleotide red and bolded) make up the columns of the table. The data for each growth condition are displayed on a separate tab.
Identification of Rho termination regions using DirectRNA-seq ±BCM.
All identified Rho termination regions are represented, defined as regions with at least one genomic coordinate with a significance score <1e−4. Rho scores were calculated for each genomic position by comparing DirectRNA-seq coverage in windows 800 nt upstream and downstream in the treated (+BCM) and untreated (-BCM) samples (see materials and methods). The 3´ genomic coordinate with the highest Rho score, DNA strand, 3´ end classification (see materials and methods), the gene annotation for the classification (details), the read coverage in the 800 nt windows upstream and downstream the 3´ end position ±BCM, the Rho score, and the p-value from the Fisher’s exact test (significance score) make up the columns of the table. An ‘undefined’ Rho score indicates one that could not be calculated due to zero reads in the -BCM downstream region. A significance score of ‘N/A’ indicates that the significance score was too low to accurately report.
Analysis of 3´ ends in 5´ UTRs and within coding sequences.
Term-seq identified 3´ ends for E. coli MG1655 (WT) that were between 200 nt upstream of an ORF and the corresponding stop codon. The genomic coordinate of the 3´ end (3´ end position), DNA strand, average RNA-seq read count of the 3´ end of both biological replicates, 3´ end classification (see materials and methods), the gene annotation for the classification (details), the location of the 3´ end relative to an ORF – upstream ORF or internal (ORF classification), and the gene annotation of the associated ORF (upstream ORF/internal details) make up the columns of the table. The data for each growth condition are displayed on a separate tab. For 3´ ends identified in WT E. coli grown to OD600 ~0.4 in LB, each 3´ end was also given a Rho score, an assessment of whether it is in a Rho termination region, and an intrinsic terminator score (see materials and methods). The read coverage in 800 nt windows upstream and downstream of each 3´ end position ±BCM from the DirectRNA-seq data, used to calculate Rho scores, is presented. An ‘undefined’ Rho score indicates one that could not be calculated due to zero reads in the ±BCM downstream region. Intrinsic terminator scores > 3.0 are suggestive of intrinsic termination (Chen et al., 2013). An ‘undefined’ intrinsic terminator score indicates one that could not be calculated because the sequence could not be folded into a recognizable secondary structure. Characterized regulatory elements in the LB 0.4 condition are also noted.
List of strains together with plasmids (tab 1) and oligonucleotides (tab 2) used in this study.