Scalable and cost-effective ribonuclease-based rRNA depletion for transcriptomics

A practical introduction to holo-omics

Holo-omics refers to the joint study of non-targeted molecular data layers from host-microbiota systems or holobionts, which is increasingly employed to disentangle the complex interactions between the elements that compose them. We navigate through the generation, analysis, and integration of omics data, focusing on the commonalities and main differences to generate and analyze the various types of omics, with a special focus on optimizing data generation and integration. We advocate for careful generation and distillation of data, followed by independent exploration and analyses of the single omic layers to obtain a better understanding of the study system, before the integration of multiple omic layers in a final model is attempted. We highlight critical decision points to achieve this aim and flag the main challenges to address complex biological questions regarding the integrative study of host-microbiota relationships.

A scalable and cost-efficient rRNA depletion approach to enrich RNAs for molecular biology investigations

Transcriptomics analyses play pivotal roles in understanding the complex regulatory networks that govern cellular processes. The abundance of rRNAs, which account for 80%-90% of total RNA in eukaryotes, limits the detection and investigation of other transcripts. While mRNAs and long noncoding RNAs have poly(A) tails that are often used for positive selection, investigations of poly(A)- RNAs, such as circular RNAs, histone mRNAs, and small RNAs, typically require the removal of the abundant rRNAs for enrichment. Current approaches to deplete rRNAs for downstream molecular biology investigations are hampered by restrictive RNA input masses and high costs. To address these challenges, we developed rRNA Removal by RNaseH (rRRR), a method to efficiently deplete rRNAs from a wide range of human, mouse, and rat RNA inputs and of varying qualities at a cost 10- to 20-fold cheaper than other approaches. We used probe-based hybridization and enzymatic digestion to selectively target and remove rRNA molecules while preserving the integrity of non-rRNA transcripts. Comparison of rRRR to two commercially available approaches showed similar rRNA depletion efficiencies and comparable off-target effects. Our developed method provides researchers with a valuable tool for investigating gene expression and regulatory mechanisms across a wide range of biological systems at an affordable price that increases the accessibility for researchers to enter the field, ultimately advancing our understanding of cellular processes.

Identification of Determinants that Allow Maintenance of High-Level Fluoroquinolone Resistance in Acinetobacter baumannii

Acinetobacter baumannii is associated with multidrug resistant (MDR) infections in healthcare settings, with fluoroquinolones such as ciprofloxacin being currently ineffective. Clinical isolates largely harbor mutations in the GyrA and TopoIV fluoroquinolone targets, as well as mutations that increase expression of drug resistance-nodulation-division (RND) efflux pumps. Factors critical for maintaining fitness levels of pump overproducers are uncharacterized despite their prevalence in clinical isolates. We here identify proteins that contribute to the fitness of FQR strains overexpressing three known RND systems using high-density insertion mutagenesis. Overproduction of the AdeFGH efflux pump caused hypersensitization to defects in outer membrane homeostatic regulation, including lesions that reduced LOS biosynthesis and blocked production of the major A. baumannii porin. In contrast, AdeAB pump overproduction, which does not affect the outer membrane pump component, was relatively tolerant to loss of these functions, consistent with outer membrane protein overproduction being the primary disruptive component. Surprisingly, overproduction of proton-transporting efflux pumps had little impact on cytosolic pH, consistent with a compensatory response to pump activity. The most striking transcriptional changes were associated with AdeFGH pump overproduction, resulting in activation of the phenylacetate (PAA) degradation regulon. Disruption of the PAA pathway resulted in cytosolic acidification and defective expression of genes involved in protection from peroxide stress. These results indicate that the RND outer membrane protein overproduction is compensated by cytoplasmic buffering and maintenance of outer membrane integrity in A. baumannii to facilitate fitness of FQR isolates.

High-throughput transcriptomics of 409 bacteria-drug pairs reveals drivers of gut microbiota perturbation

Many drugs can perturb the gut microbiome, potentially leading to negative health consequences. However, mechanisms of most microorganism-drug responses have not been elucidated at the genetic level. Using high-throughput bacterial transcriptomics, we systematically characterized the gene expression profiles of prevalent human gut bacteria exposed to the most frequently prescribed orally administered pharmaceuticals. Across >400 drug-microorganism pairs, significant and reproducible transcriptional responses were observed, including pathways involved in multidrug resistance, metabolite transport, tartrate metabolism and riboflavin biosynthesis. Importantly, we discovered that statin-mediated upregulation of the AcrAB-TolC efflux pump in Bacteroidales species enhances microbial sensitivity to vitamin A and secondary bile acids. Moreover, gut bacteria carrying acrAB-tolC genes are depleted in patients taking simvastatin, suggesting that drug-efflux interactions generate collateral toxicity that depletes pump-containing microorganisms from patient microbiomes. This study provides a resource to further understand the drivers of drug-mediated microbiota shifts for better informed clinical interventions.

Emerging applications of single-cell profiling in precision medicine of atherosclerosis

Atherosclerosis is a chronic, progressive, inflammatory disease that occurs in the arterial wall. Despite recent advancements in treatment aimed at improving efficacy and prolonging survival, atherosclerosis remains largely incurable. In this review, we discuss emerging single-cell sequencing techniques and their novel insights into atherosclerosis. We provide examples of single-cell profiling studies that reveal phenotypic characteristics of atherosclerosis plaques, blood, liver, and the intestinal tract. Additionally, we highlight the potential clinical applications of single-cell analysis and propose that combining this approach with other techniques can facilitate early diagnosis and treatment, leading to more accurate medical interventions.

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

IMPORTANCE

Microbes present one of the most diverse sources of biochemistry in nature, and mRNA sequencing provides a comprehensive view of this biological activity by quantitatively measuring microbial transcriptomes. However, efficient mRNA capture for sequencing presents significant challenges in prokaryotes as mRNAs are not poly-adenylated and typically make up less than 5% of total RNA compared with rRNAs that exceed 80%. Recently developed methods for sequencing bacterial mRNA typically rely on depleting rRNA by tiling large probe sets against rRNAs; however, such approaches are expensive, time-consuming, and challenging to scale to varied bacterial species and complex microbial communities. Therefore, we developed EMBR-seq+, a method that requires fewer than 10 short oligonucleotides per rRNA to achieve up to 99% rRNA depletion in diverse bacterial species. Finally, EMBR-seq+ resulted in a deeper view of the transcriptome, enabling systematic quantification of how microbial interactions result in altering the transcriptional state of bacteria within co-cultures.

Rational probe design for efficient rRNA depletion and improved metatranscriptomic analysis of human microbiomes

The microbiota that colonize the human gut and other tissues are dynamic, varying both in composition and functional state between individuals and over time. Gene expression measurements can provide insights into microbiome composition and function. However, efficient and unbiased removal of microbial ribosomal RNA (rRNA) presents a barrier to acquiring metatranscriptomic data. Here we describe a probe set that achieves efficient enzymatic rRNA removal of complex human-associated microbial communities. We demonstrate that the custom probe set can be further refined through an iterative design process to efficiently deplete rRNA from a range of human microbiome samples. Using synthetic nucleic acid spike-ins, we show that the rRNA depletion process does not introduce substantial quantitative error in gene expression profiles. Successful rRNA depletion allows for efficient characterization of taxonomic and functional profiles, including during the development of the human gut microbiome. The pan-human microbiome enzymatic rRNA depletion probes described here provide a powerful tool for studying the transcriptional dynamics and function of the human microbiome.

Single-cell massively-parallel multiplexed microbial sequencing (M3-seq) identifies rare bacterial populations and profiles phage infection

Bacterial populations are highly adaptive. They can respond to stress and survive in shifting environments. How the behaviours of individual bacteria vary during stress, however, is poorly understood. To identify and characterize rare bacterial subpopulations, technologies for single-cell transcriptional profiling have been developed. Existing approaches show some degree of limitation, for example, in terms of number of cells or transcripts that can be profiled. Due in part to these limitations, few conditions have been studied with these tools. Here we develop massively-parallel, multiplexed, microbial sequencing (M3-seq)-a single-cell RNA-sequencing platform for bacteria that pairs combinatorial cell indexing with post hoc rRNA depletion. We show that M3-seq can profile bacterial cells from different species under a range of conditions in single experiments. We then apply M3-seq to hundreds of thousands of cells, revealing rare populations and insights into bet-hedging associated with stress responses and characterizing phage infection.

Current concepts, advances, and challenges in deciphering the human microbiota with metatranscriptomics

Metatranscriptomics refers to the analysis of the collective microbial transcriptome of a sample. Its increased utilization for the characterization of human-associated microbial communities has enabled the discovery of many disease-state related microbial activities. Here, we review the principles of metatranscriptomics-based analysis of human-associated microbial samples. We describe strengths and weaknesses of popular sample preparation, sequencing, and bioinformatics approaches and summarize strategies for their use. We then discuss how human-associated microbial communities have recently been examined and how their characterization may change. We conclude that metatranscriptomics insights into human microbiotas under health and disease have not only expanded our knowledge on human health, but also opened avenues for rational antimicrobial drug use and disease management.

Droplet-based high-throughput single microbe RNA sequencing by smRandom-seq

Bacteria colonize almost all parts of the human body and can differ significantly. However, the population level transcriptomics measurements can only describe the average bacteria population behaviors, ignoring the heterogeneity among bacteria. Here, we report a droplet-based high-throughput single-microbe RNA-seq assay (smRandom-seq), using random primers for in situ cDNA generation, droplets for single-microbe barcoding, and CRISPR-based rRNA depletion for mRNA enrichment. smRandom-seq showed a high species specificity (99%), a minor doublet rate (1.6%), a reduced rRNA percentage (32%), and a sensitive gene detection (a median of ~1000 genes per single E. coli). Furthermore, smRandom-seq successfully captured transcriptome changes of thousands of individual E. coli and discovered a few antibiotic resistant subpopulations displaying distinct gene expression patterns of SOS response and metabolic pathways in E. coli population upon antibiotic stress. smRandom-seq provides a high-throughput single-microbe transcriptome profiling tool that will facilitate future discoveries in microbial resistance, persistence, microbe-host interaction, and microbiome research.

Single-cell transcriptomics and data analyses for prokaryotes-Past, present and future concepts

Transcriptomics, or more specifically mRNA sequencing, is a powerful tool to study gene expression at the single-cell level (scRNA-seq) which enables new insights into a plethora of biological processes. While methods for single-cell RNA-seq in eukaryotes are well established, application to prokaryotes is still challenging. Reasons for that are rigid and diverse cell wall structures hampering lysis, the lack of polyadenylated transcripts impeding mRNA enrichment, and minute amounts of RNA requiring amplification steps before sequencing. Despite those obstacles, several promising scRNA-seq approaches for bacteria have been published recently, albeit difficulties in the experimental workflow and data processing and analysis remain. In particular, bias is often introduced by amplification which makes it difficult to distinguish between technical noise and biological variation. Future optimization of experimental procedures and data analysis algorithms are needed for the improvement of scRNA-seq but also to aid in the emergence of prokaryotic single-cell multi-omics. to help address 21st century challenges in the biotechnology and health sector.

Small RNA MTS1338 Configures a Stress Resistance Signature in Mycobacterium tuberculosis

In the course of evolution, Mycobacterium tuberculosis (Mtb), the etiological agent of tuberculosis, has developed sophisticated strategies to evade host immune response, including the synthesis of small non-coding RNAs (sRNAs), which regulate post-transcriptional pathways involved in the stress adaptation of mycobacteria. sRNA MTS1338 is upregulated in Mtb during its infection of cultured macrophages and in the model of chronic tuberculosis, suggesting involvement in host-pathogen interactions. Here, we analyzed the role of MTS1338 in the Mtb response to macrophage-like stresses in vitro. The Mtb strain overexpressing MTS1338 demonstrated enhanced survival ability under low pH, nitrosative, and oxidative stress conditions simulating the antimicrobial environment inside macrophages. Transcriptomic analysis revealed that in MTS1338-overexpressing Mtb, the stress factors led to the activation of a number of transcriptional regulators, toxin-antitoxin modules, and stress chaperones, about half of which coincided with the genes induced in Mtb phagocytosed by macrophages. We determined the MTS1338 "core regulon", consisting of 11 genes that were activated in all conditions under MTS1338 overexpression. Our findings indicate that MTS1338 is a stress-induced sRNA that promotes Mtb survival in macrophages by triggering adaptive transcriptional mechanisms in response to host antimicrobial defense reactions.

Blocking Abundant RNA Transcripts by High-Affinity Oligonucleotides during Transcriptome Library Preparation

BACKGROUND

RNA sequencing has become the gold standard for transcriptome analysis but has an inherent limitation of challenging quantification of low-abundant transcripts. In contrast to microarray technology, RNA sequencing reads are proportionally divided in function of transcript abundance. Therefore, low-abundant RNAs compete against highly abundant - and sometimes non-informative - RNA species.

RESULTS

We developed an easy-to-use strategy based on high-affinity RNA-binding oligonucleotides to block reverse transcription and PCR amplification of specific RNA transcripts, thereby substantially reducing their abundance in the final sequencing library. To demonstrate the broad application potential of our method, we applied it to different transcripts and library preparation strategies, including YRNAs in small RNA sequencing of human blood plasma, mitochondrial rRNAs in both 3' end sequencing and long-read sequencing, and MALAT1 in single-cell 3' end sequencing. We demonstrate that the blocking strategy is highly efficient, reproducible, specific, and generally results in better transcriptome coverage and complexity.

CONCLUSION

Our method does not require modifications of the library preparation procedure apart from simply adding blocking oligonucleotides to the RT reaction and can thus be easily integrated into virtually any RNA sequencing library preparation protocol.

Profiling of bacterial transcriptome from ultra-low input with MiniBac-seq

There is a lack of appropriate methods for preparing bacterial RNA-seq library with ultra-low amount of RNA. To address this issue, we developed miniBac-seq, a strand-specific method for high-quality library construction from sub-nanogram of total RNA, which is 100-fold lower than the current benchmark kit and dramatically reduces preparation cost ($28 + $15 × samples). We further demonstrated the high sensitivity of miniBac-seq via detecting more than 500 genes from amount of total RNA equivalent to that of a single bacterial cell. Finally, we profiled the transcriptome of growth-arrested bacteria in isogenic culture of Escherichia coli. This subpopulation of bacteria is generally low in abundance but is a potent reservoir of antibiotic persistence, and their gene expression has been largely unknown due to technical limitations. Using miniBac-seq, we identified potential molecular driver towards arrested growth as well as antibiotic tolerance. Our method thus expands the capacity to quantify bacterial transcriptome in situ, which is useful to the understanding of bacterial physiology and regulation in their native contexts.

Efficient and specific DNA oligonucleotide rRNA probe-based rRNA removal in Talaromyces marneffei

Emerging evidence showed that lncRNAs play important roles in a wide range of biological processes of fungi such as Saccharomyces cerevisiae. However, systemic identification of lncRNAs in non-model fungi is a challenging task as the efficiency of rRNA removal has been proved to be affected by mismatches of universal rRNA-targeting probes of commercial kits, which forces deeper sequencing depth and increases costs. Here, we developed a low-cost and simple rRNA depletion method (rProbe) that could efficiently remove more than 99% rRNA in both yeast and mycelium samples of Talaromyces marneffei. The efficiency and robustness of rProbe were demonstrated to outperform the Illumina Ribo-Zero kit. Using rProbe RNA-seq, we identified 115 differentially expressed lncRNAs and constructed lncRNA-mRNA co-expression network related to dimorphic switch of T. marneffei. Our rRNA removal method has the potential to be a useful tool to explore non-coding transcriptomes of non-model fungi by adjusting rRNA probe sequences species specifically.

Therapeutic efficacy of antimalarial drugs targeting DosRS signaling in Mycobacterium abscessus

A search for alternative Mycobacterium abscessus treatments led to our interest in the two-component regulator DosRS, which, in Mycobacterium tuberculosis, is required for the bacterium to establish a state of nonreplicating, drug-tolerant persistence in response to a variety of host stresses. We show here that the genetic disruption of dosRS impairs the adaptation of M. abscessus to hypoxia, resulting in decreased bacterial survival after oxygen depletion, reduced tolerance to a number of antibiotics in vitro and in vivo, and the inhibition of biofilm formation. We determined that three antimalarial drugs or drug candidates, artemisinin, OZ277, and OZ439, can target DosS-mediated hypoxic signaling in M. abscessus and recapitulate the phenotypic effects of genetically disrupting dosS. OZ439 displayed bactericidal activity comparable to standard-of-care antibiotics in chronically infected mice, in addition to potentiating the activity of antibiotics used in combination. The identification of antimalarial drugs as potent inhibitors and adjunct inhibitors of M. abscessus in vivo offers repurposing opportunities that could have an immediate impact in the clinic.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Unique Features of Mycobacterium abscessus Biofilms Formed in Synthetic Cystic Fibrosis Medium

Characterizing Mycobacterium abscessus complex (MABSC) biofilms under host-relevant conditions is essential to the design of informed therapeutic strategies targeted to this persistent, drug-tolerant, population of extracellular bacilli. Using synthetic cystic fibrosis medium (SCFM) which we previously reported to closely mimic the conditions encountered by MABSC in actual cystic fibrosis (CF) sputum and a new model of biofilm formation, we show that MABSC biofilms formed under these conditions are substantially different from previously reported biofilms grown in standard laboratory media in terms of their composition, gene expression profile and stress response. Extracellular DNA (eDNA), mannose-and glucose-containing glycans and phospholipids, rather than proteins and mycolic acids, were revealed as key extracellular matrix (ECM) constituents holding clusters of bacilli together. None of the environmental cues previously reported to impact biofilm development had any significant effect on SCFM-grown biofilms, most likely reflecting the fact that SCFM is a nutrient-rich environment in which MABSC finds a variety of ways of coping with stresses. Finally, molecular determinants were identified that may represent attractive new targets for the development of adjunct therapeutics targeting MABSC biofilms in persons with CF.

A Novel Family of RNA-Binding Proteins Regulate Polysaccharide Metabolism in Bacteroides thetaiotaomicron

Human gut microbiome composition is constantly changing, and diet is a major driver of these changes. Gut microbial species that persist in mammalian hosts for long periods of time must possess mechanisms for sensing and adapting to nutrient shifts to avoid being outcompeted. Global regulatory mechanisms mediated by RNA-binding proteins (RBPs) that govern responses to nutrient shifts have been characterized in Proteobacteria and Firmicutes but remain undiscovered in the Bacteroidetes. Here, we report the identification of RBPs that are broadly distributed across the Bacteroidetes, with many genomes encoding multiple copies. Genes encoding these RBPs are highly expressed in many Bacteroides species. A purified RBP, RbpB, from Bacteroides thetaiotaomicron binds to single-stranded RNA in vitro with an affinity similar to other characterized regulatory RBPs. B. thetaiotaomicron mutants lacking RBPs show dramatic shifts in expression of polysaccharide utilization and capsular polysaccharide loci, suggesting that these RBPs may act as global regulators of polysaccharide metabolism. A B. thetaiotaomicron ΔrbpB mutant shows a growth defect on dietary sugars belonging to the raffinose family of oligosaccharides (RFOs). The ΔrbpB mutant had reduced expression of BT1871, encoding a predicted RFO-degrading melibiase, compared to the wild-type strain. Mutation of BT1871 confirmed that the enzyme it encodes is essential for growth on melibiose and promotes growth on the RFOs raffinose and stachyose. Our data reveal that RbpB is required for optimal expression of BT1871 and other polysaccharide-related genes, suggesting that we have identified an important new family of global regulatory proteins in the Bacteroidetes. IMPORTANCE The human colon houses hundreds of bacterial species, including many belonging to the genus Bacteroides, that aid in breaking down our food to keep us healthy. Bacteroides have many genes responsible for breaking down different dietary carbohydrates, and complex regulatory mechanisms ensure that specific genes are only expressed when the right carbohydrates are available. In this study, we discovered that Bacteroides use a family of RNA-binding proteins as global regulators to coordinate expression of carbohydrate utilization genes. The ability to turn different carbohydrate utilization genes on and off in response to changing nutrient conditions is critical for Bacteroides to live successfully in the gut, and thus the new regulators we have identified may be important for life in the host.

RiboRid: A low cost, advanced, and ultra-efficient method to remove ribosomal RNA for bacterial transcriptomics

RNA sequencing techniques have enabled the systematic elucidation of gene expression (RNA-Seq), transcription start sites (differential RNA-Seq), transcript 3′ ends (Term-Seq), and post-transcriptional processes (ribosome profiling). The main challenge of transcriptomic studies is to remove ribosomal RNAs (rRNAs), which comprise more than 90% of the total RNA in a cell. Here, we report a low-cost and robust bacterial rRNA depletion method, RiboRid, based on the enzymatic degradation of rRNA by thermostable RNase H. This method implemented experimental considerations to minimize nonspecific degradation of mRNA and is capable of depleting pre-rRNAs that often comprise a large portion of RNA, even after rRNA depletion. We demonstrated the highly efficient removal of rRNA up to a removal efficiency of 99.99% for various transcriptome studies, including RNA-Seq, Term-Seq, and ribosome profiling, with a cost of approximately $10 per sample. This method is expected to be a robust method for large-scale high-throughput bacterial transcriptomic studies.

Removal of ribosomal RNAs, a major constituent (over 90%) of cellular RNA is a critical experimental step for transcriptomic studies that deal with messenger RNAs. In this manuscript, we describe a robust method to subtract ribosomal RNA from various RNA samples. The method is based on the enzymatic degradation of target RNA by short complementary DNA and RNA:DNA duplex specific nuclease. The method comprises carefully designed experimental procedures to minimize experimental bias and unwanted removal of messenger RNAs. We validate the method on various types of transcriptomic studies for seven diverse bacterial species. This method successfully removed ribosomal RNA with over 99% of efficiency and it was comparable to commercial systems even for degraded RNA samples at a fraction of a cost.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Optimized design of antisense oligomers for targeted rRNA depletion

RNA sequencing (RNA-seq) is extensively used to quantify gene expression transcriptome-wide. Although often paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA), many applications require alternate approaches to counteract the high proportion of ribosomal RNA (rRNA) in total RNA. Recently, digestion using RNaseH and antisense DNA oligomers tiling target rRNAs has emerged as an alternative to commercial rRNA depletion kits. Here, we present a streamlined, more economical RNaseH-mediated rRNA depletion with substantially lower up-front costs, using shorter antisense oligos only sparsely tiled along the target RNA in a 5-min digestion reaction. We introduce a novel Web tool, Oligo-ASST, that simplifies oligo design to target regions with optimal thermodynamic properties, and additionally can generate compact, common oligo pools that simultaneously target divergent RNAs, e.g. across different species. We demonstrate the efficacy of these strategies by generating rRNA-depletion oligos for Xenopus laevis and for zebrafish, which expresses two distinct versions of rRNAs during embryogenesis. The resulting RNA-seq libraries reduce rRNA to <5% of aligned reads, on par with poly(A) selection, and also reveal expression of many non-adenylated RNA species. Oligo-ASST is freely available at https://mtleelab.pitt.edu/oligo to design antisense oligos for any taxon or to target any abundant RNA for depletion.

Methods for isolation of messenger RNA from biological samples

RNA molecules contain many chemical modifications that can regulate a variety of biological processes. Messenger RNA (mRNA) molecules are critical components in the central dogma of molecular biology. The discovery of reversible chemical modifications in eukaryotic mRNA brings forward a new research field in RNA modification-mediated regulation of gene expression. The modifications in mRNA generally exist in low abundance. The use of highly pure mRNA is critical for the confident identification of new modifications as well as for the accurate quantification of existing modifications in mRNA. In addition, isolation of highly pure mRNA is the first step in many biological research studies. Therefore, the methods for isolating highly pure mRNA are important for mRNA-based downstream studies. A variety of methods for isolating mRNA have been developed in the past few decades and new methods continuously emerge. This review focuses on the methodologies and protocols for isolating mRNA populations. In addition, we discuss the advantages and limitations of these methods. We hope this paper will provide a general view of mRNA isolation strategies and facilitate studies that involve mRNA modifications and functions.

Coupled Transcriptomics for Differential Expression Analysis and Determination of Transcription Start Sites: Design and Bioinformatics

The term coupled transcriptomics is coined to describe a design of an RNA-seq experiment intended for both differential expression analysis and genome-wide determination of the transcription start sites (TSS). The minimal requirements for the first analysis are two experimental conditions with at least two biological replicates enabling statistical tests. The second analysis involves the bioinformatics comparison of the data generated from a control RNA-seq library with another library enriched in primary transcripts using Terminator^™ 5'-phosphate-dependent exonuclease, in an experiment denominated differential RNA-seq (dRNA-seq). Usually, dRNA-seq is carried out with specific protocols for library construction, different of those used for common differential expression analysis. Our experimental design allows to use the same data for both analyses, reducing the number of libraries to be generated and sequenced. This is a guide for designing a coupled transcriptomics experiment and for the subsequent bioinformatics procedures. The proposed methods can be applied to the detection and study of small RNA genes.

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

RNA-Seq enables the identification and quantification of RNA molecules, often with the aim of detecting differentially expressed genes (DEGs). Although RNA-Seq evolved into a standard technique, there is no universal gold standard for these data's computational analysis. On top of that, previous studies proved the irreproducibility of RNA-Seq studies. Here, we present a portable, scalable, and parallelizable Nextflow RNA-Seq pipeline to detect DEGs, which assures a high level of reproducibility. The pipeline automatically takes care of common pitfalls, such as ribosomal RNA removal and low abundance gene filtering. Apart from various visualizations for the DEG results, we incorporated downstream pathway analysis for common species as Homo sapiens and Mus musculus. We evaluated the DEG detection functionality while using qRT-PCR data serving as a reference and observed a very high correlation of the logarithmized gene expression fold changes.

Efficient and cost-effective bacterial mRNA sequencing from low input samples through ribosomal RNA depletion

BACKGROUND

RNA sequencing is a powerful approach to quantify the genome-wide distribution of mRNA molecules in a population to gain deeper understanding of cellular functions and phenotypes. However, unlike eukaryotic cells, mRNA sequencing of bacterial samples is more challenging due to the absence of a poly-A tail that typically enables efficient capture and enrichment of mRNA from the abundant rRNA molecules in a cell. Moreover, bacterial cells frequently contain 100-fold lower quantities of RNA compared to mammalian cells, which further complicates mRNA sequencing from non-cultivable and non-model bacterial species. To overcome these limitations, we report EMBR-seq (Enrichment of mRNA by Blocked rRNA), a method that efficiently depletes 5S, 16S and 23S rRNA using blocking primers to prevent their amplification.

RESULTS

EMBR-seq results in 90% of the sequenced RNA molecules from an E. coli culture deriving from mRNA. We demonstrate that this increased efficiency provides a deeper view of the transcriptome without introducing technical amplification-induced biases. Moreover, compared to recent methods that employ a large array of oligonucleotides to deplete rRNA, EMBR-seq uses a single or a few oligonucleotides per rRNA, thereby making this new technology significantly more cost-effective, especially when applied to varied bacterial species. Finally, compared to existing commercial kits for bacterial rRNA depletion, we show that EMBR-seq can be used to successfully quantify the transcriptome from more than 500-fold lower starting total RNA.

CONCLUSIONS

EMBR-seq provides an efficient and cost-effective approach to quantify global gene expression profiles from low input bacterial samples.

Prokaryotic single-cell RNA sequencing by in situ combinatorial indexing

Despite longstanding appreciation of gene expression heterogeneity in isogenic bacterial populations, affordable and scalable technologies for studying single bacterial cells have been limited. Although single-cell RNA sequencing (scRNA-seq) has revolutionized studies of transcriptional heterogeneity in diverse eukaryotic systems1-13, the application of scRNA-seq to prokaryotes has been hindered by their extremely low mRNA abundance14-16, lack of mRNA polyadenylation and thick cell walls17. Here, we present prokaryotic expression profiling by tagging RNA in situ and sequencing (PETRI-seq)-a low-cost, high-throughput prokaryotic scRNA-seq pipeline that overcomes these technical obstacles. PETRI-seq uses in situ combinatorial indexing11,12,18 to barcode transcripts from tens of thousands of cells in a single experiment. PETRI-seq captures single-cell transcriptomes of Gram-negative and Gram-positive bacteria with high purity and low bias, with median capture rates of more than 200 mRNAs per cell for exponentially growing Escherichia coli. These characteristics enable robust discrimination of cell states corresponding to different phases of growth. When applied to wild-type Staphylococcus aureus, PETRI-seq revealed a rare subpopulation of cells undergoing prophage induction. We anticipate that PETRI-seq will have broad utility in defining single-cell states and their dynamics in complex microbial communities.

Construction of prokaryotic strand-specific primary-transcripts saturated RNASeq library by controlled heat magnesium-dependent mRNA degradation

The main limiting factors for RNA-Seq analysis are quality and quantity of the isolated mRNA. In prokaryotes, the proportion of messenger RNA to total RNA is rather low. Therefore, the main strategy of library preparation for sequencing is mRNA enrichment. Ribosomal and transfer RNAs, both monophosphorylated at the 5'-ends, are the major fractions of total RNA, while the bulk of primary transcripts is triphosphorylated at the 5'-teminus. Due to its low molecular weight, transfer RNA could be easily removed by a quick precipitation in LiCl solution. Ribosomal RNA may be degraded enzymatically by 5'-end terminal exonuclease XRN-1. These steps allow enriching samples in mRNA during the first stages of RNA-Seq library preparation. The desired level of fragmentation of enriched mRNA necessary for the 2^nd generation sequencing can be controlled by the duration of incubation at elevated temperatures in the presence of Mg²⁺-ions. Here, we describe a simple protocol for construction of the primary prokaryotic mRNA-saturated library without long depletion procedures.

Improved bacterial RNA-seq by Cas9-based depletion of ribosomal RNA reads

A major challenge for RNA-seq analysis of gene expression is to achieve sufficient coverage of informative nonribosomal transcripts. In eukaryotic samples, this is typically achieved by selective oligo(dT)-priming of messenger RNAs to exclude ribosomal RNA (rRNA) during cDNA synthesis. However, this strategy is not compatible with prokaryotes in which functional transcripts are generally not polyadenylated. To overcome this, we adopted DASH (depletion of abundant sequences by hybridization), initially developed for eukaryotic cells, to improve both the sensitivity and depth of bacterial RNA-seq. DASH uses the Cas9 nuclease to remove unwanted cDNA sequences prior to library amplification. We report the design, evaluation, and optimization of DASH experiments for standard bacterial short-read sequencing approaches, including software for automated guide RNA (gRNA) design for Cas9-mediated cleavage in bacterial rDNA sequences. Using these gRNA pools, we effectively removed rRNA reads (56%-86%) in RNA-seq libraries from two different model bacteria, the Gram-negative pathogen Salmonella enterica and the anaerobic gut commensal Bacteroides thetaiotaomicron DASH works robustly, even with subnanogram amounts of input RNA. Its efficiency, high sensitivity, ease of implementation, and low cost (∼$5 per sample) render DASH an attractive alternative to rRNA removal protocols, in particular for material-constrained studies where conventional ribodepletion techniques fail.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

A Simple, Cost-Effective, and Robust Method for rRNA Depletion in RNA-Sequencing Studies

The profiling of gene expression by RNA sequencing (RNA-seq) has enabled powerful studies of global transcriptional patterns in all organisms, including bacteria. Because the vast majority of RNA in bacteria is rRNA, it is standard practice to deplete the rRNA from a total RNA sample such that the reads in an RNA-seq experiment derive predominantly from mRNA. One of the most commonly used commercial kits for rRNA depletion, the Ribo-Zero kit from Illumina, was recently discontinued abruptly and for an extended period of time. Here, we report the development of a simple, cost-effective, and robust method for depleting rRNA that can be easily implemented by any lab or facility. We first developed an algorithm for designing biotinylated oligonucleotides that will hybridize tightly and specifically to the 23S, 16S, and 5S rRNAs from any species of interest. Precipitation of these oligonucleotides bound to rRNA by magnetic streptavidin-coated beads then depletes rRNA from a complex, total RNA sample such that ∼75 to 80% of reads in a typical RNA-seq experiment derive from mRNA. Importantly, we demonstrate a high correlation of RNA abundance or fold change measurements in RNA-seq experiments between our method and the Ribo-Zero kit. Complete details on the methodology are provided, including open-source software for designing oligonucleotides optimized for any bacterial species or community of interest.IMPORTANCE The ability to examine global patterns of gene expression in microbes through RNA sequencing has fundamentally transformed microbiology. However, RNA-seq depends critically on the removal of rRNA from total RNA samples. Otherwise, rRNA would comprise upward of 90% of the reads in a typical RNA-seq experiment, limiting the reads coming from mRNA or requiring high total read depth. A commonly used kit for rRNA subtraction from Illumina was recently unavailable for an extended period of time, disrupting routine rRNA depletion. Here, we report the development of a "do-it-yourself" kit for rapid, cost-effective, and robust depletion of rRNA from total RNA. We present an algorithm for designing biotinylated oligonucleotides that will hybridize to the rRNAs from a target set of species. We then demonstrate that the designed oligonucleotides enable sufficient rRNA depletion to produce RNA-seq data with 75 to 80% of reads coming from mRNA. The methodology presented should enable RNA-seq studies on any species or metagenomic sample of interest.

Meeting report of the third annual Tri-Service Microbiome Consortium symposium

The Tri-Service Microbiome Consortium (TSMC) was founded to enhance collaboration, coordination, and communication of microbiome research among U.S. Department of Defense (DoD) organizations and to facilitate resource, material and information sharing among consortium members. The 2019 annual symposium was held 22-24 October 2019 at Wright-Patterson Air Force Base in Dayton, OH. Presentations and discussions centered on microbiome-related topics within five broad thematic areas: 1) human microbiomes; 2) transitioning products into Warfighter solutions; 3) environmental microbiomes; 4) engineering microbiomes; and 5) microbiome simulation and characterization. Collectively, the symposium provided an update on the scope of current DoD microbiome research efforts, highlighted innovative research being done in academia and industry that can be leveraged by the DoD, and fostered collaborative opportunities. This report summarizes the presentations and outcomes of the 3rd annual TSMC symposium.