Efficient depletion of ribosomal RNA for RNA sequencing in planarians

RNaseH-based ribodepletion of total planarian RNA improves detection of longer and non-polyadenylated transcripts

The overwhelming majority of RNA species isolated from cells or tissues using organic extraction are ribosomal RNAs (rRNA), whereas a relatively small percentage are messenger RNAs (mRNA). For studies that seek to detect mRNA transcripts and measure changes in their expression, this lopsided ratio of desired transcripts to undesired transcripts creates a significant challenge to obtaining sensitive and reproducible results. One method for improving mRNA detection is to selectively amplify polyadenylated (polyA) mRNA molecules when generating RNA-seq libraries, a strategy that is generally very successful in many species. However, this strategy is less effective when starting with total RNA from some species e.g., the planarian species Schmidtea mediterranea (S.med), as it generates libraries that still contain significant and variable amounts of rRNA reads. Further, commercially available ribodepletion kits do not efficiently deplete rRNAs from these samples because their sequences are divergent from mammalian rRNAs. Here we report a customized, optimized, and economical ribodepletion strategy than allows the generation of comprehensive RNA-seq libraries with less than one percent rRNA contamination. We show that this method improves transcript detection, particularly for those without polyA tails (e.g., core histones) and those that are relatively long (e.g., microtubule motor proteins). Using this custom ribodepletion approach, we also detected many transcripts that are not represented in the most recent set of S.med gene annotations, including a subset that are likely expressed transposable elements (TEs). To facilitate future differential expression analyses of these newly identified loci, we created both an annotation file of the new loci we identified and a bioinformatic pipeline for generating additional annotations from future libraries. As significant recent research shows that TE activation is regulated and functionally important, the resources provided here will provide a starting point for investigating such mechanisms in planarians and other species with less conserved rRNA sequences.

A single workflow for multi-species blood transcriptomics

BACKGROUND

Blood transcriptomic analysis is widely used to provide a detailed picture of a physiological state with potential outcomes for applications in diagnostics and monitoring of the immune response to vaccines. However, multi-species transcriptomic analysis is still a challenge from a technological point of view and a standardized workflow is urgently needed to allow interspecies comparisons.

RESULTS

Here, we propose a single and complete total RNA-Seq workflow to generate reliable transcriptomic data from blood samples from humans and from animals typically used in preclinical models. Blood samples from a maximum of six individuals and four different species (rabbit, non-human primate, mouse and human) were extracted and sequenced in triplicates. The workflow was evaluated using different wet-lab and dry-lab criteria, including RNA quality and quantity, the library molarity, the number of raw sequencing reads, the Phred-score quality, the GC content, the performance of ribosomal-RNA and globin depletion, the presence of residual DNA, the strandness, the percentage of coding genes, the number of genes expressed, and the presence of saturation plateau in rarefaction curves. We identified key criteria and their associated thresholds to be achieved for validating the transcriptomic workflow. In this study, we also generated an automated analysis of the transcriptomic data that streamlines the validation of the dataset generated.

CONCLUSIONS

Our study has developed an end-to-end workflow that should improve the standardization and the inter-species comparison in blood transcriptomics studies. In the context of vaccines and drug development, RNA sequencing data from preclinical models can be directly compared with clinical data and used to identify potential biomarkers of value to monitor safety and efficacy.

CRISPR/Cas9-based depletion of 16S ribosomal RNA improves library complexity of single-cell RNA-sequencing in planarians

BACKGROUND

Single-cell RNA-sequencing (scRNA-seq) relies on PCR amplification to retrieve information from vanishingly small amounts of starting material. To selectively enrich mRNA from abundant non-polyadenylated transcripts, poly(A) selection is a key step during library preparation. However, some transcripts, such as mitochondrial genes, can escape this elimination and overwhelm libraries. Often, these transcripts are removed in silico, but whether physical depletion improves detection of rare transcripts in single cells is unclear.

RESULTS

We find that a single 16S ribosomal RNA is widely enriched in planarian scRNA-seq datasets, independent of the library preparation method. To deplete this transcript from scRNA-seq libraries, we design 30 single-guide RNAs spanning its length. To evaluate the effects of depletion, we perform a side-by-side comparison of the effects of eliminating the 16S transcript and find a substantial increase in the number of genes detected per cell, coupled with virtually complete loss of the 16S RNA. Moreover, we systematically determine that library complexity increases with a limited number of PCR cycles following CRISPR treatment. When compared to in silico depletion of 16S, physically removing it reduces dropout rates, retrieves more clusters, and reveals more differentially expressed genes.

CONCLUSIONS

Our results show that abundant transcripts reduce the retrieval of informative transcripts in scRNA-seq and distort the analysis. Physical removal of these contaminants enables the detection of rare transcripts at lower sequencing depth, and also outperforms in silico depletion. Importantly, this method can be easily customized to deplete any abundant transcript from scRNA-seq libraries.

CRISPR/Cas9-based depletion of 16S ribosomal RNA improves library complexity of single-cell RNA-sequencing

Background

Single-cell RNA-sequencing (scRNA-seq) relies on PCR amplification to retrieve information from vanishingly small amounts of starting material. To selectively enrich mRNA from abundant non-polyadenylated transcripts, poly(A) selection is a key step during library preparation. However, some transcripts, such as mitochondrial genes, can escape this elimination and overwhelm libraries. Often, these transcripts are removed in silico, but whether physical depletion improves detection of rare transcripts in single cells is unclear.

Results

We find that a single 16S ribosomal RNA is widely enriched in planarian scRNA-seq datasets, independent of the library preparation method. To deplete this transcript from scRNA-seq libraries, we design 30 single-guide RNAs spanning its length. To evaluate the effects of depletion, we perform a side-by-side comparison of the effects of eliminating the 16S transcript and find a substantial increase in the number of genes detected per cell, coupled with virtually complete loss of the 16S RNA. Moreover, we systematically determine that library complexity increases with a limited number of PCR cycles following CRISPR treatment. When compared to in silico depletion of 16S, physically removing it reduces dropout rates, retrieves more clusters, and reveals more differentially-expressed genes.

Conclusions

Our results show that abundant transcripts reduce the retrieval of informative transcripts in scRNA-seq and distort the analysis. Physical removal of these contaminants enables the detection of rare transcripts at lower sequencing depth, and also outperforms in silico depletion. Importantly, this method can be easily customized to deplete any abundant transcript from scRNA-seq libraries.

Nectary development in Cleome violacea

Nectaries are a promising frontier for plant evo-devo research, and are particularly fascinating given their diversity in form, position, and secretion methods across angiosperms. Emerging model systems permit investigations of the molecular basis for nectary development and nectar secretion across a range of taxa, which addresses fundamental questions about underlying parallelisms and convergence. Herein, we explore nectary development and nectar secretion in the emerging model taxa, Cleome violacea (Cleomaceae), which exhibits a prominent adaxial nectary. First, we characterized nectary anatomy and quantified nectar secretion to establish a foundation for quantitative and functional gene experiments. Next, we leveraged RNA-seq to establish gene expression profiles of nectaries across three key stages of development: pre-anthesis, anthesis, and post-fertilization. We then performed functional studies on five genes that were putatively involved in nectary and nectar formation: CvCRABSCLAW (CvCRC), CvAGAMOUS (CvAG), CvSHATTERPROOF (CvSHP), CvSWEET9, and a highly expressed but uncharacterized transcript. These experiments revealed a high degree of functional convergence to homologues from other core Eudicots, especially Arabidopsis. CvCRC, redundantly with CvAG and CvSHP, are required for nectary initiation. Concordantly, CvSWEET9 is essential for nectar formation and secretion, which indicates that the process is eccrine based in C. violacea. While demonstration of conservation is informative to our understanding of nectary evolution, questions remain. For example, it is unknown which genes are downstream of the developmental initiators CvCRC, CvAG, and CvSHP, or what role the TCP gene family plays in nectary initiation in this family. Further to this, we have initiated a characterization of associations between nectaries, yeast, and bacteria, but more research is required beyond establishing their presence. Cleome violacea is an excellent model for continued research into nectary development because of its conspicuous nectaries, short generation time, and close taxonomic distance to Arabidopsis.

Inhibition of EZH2 Causes Retrotransposon Derepression and Immune Activation in Porcine Lung Alveolar Macrophages

Alveolar macrophages (AMs) form the first defense line against various respiratory pathogens, and their immune response has a profound impact on the outcome of respiratory infection. Enhancer of zeste homolog 2 (EZH2), which catalyzes the trimethylation of H3K27 for epigenetic repression, has gained increasing attention for its immune regulation function, yet its exact function in AMs remains largely obscure. Using porcine 3D4/21 AM cells as a model, we characterized the transcriptomic and epigenomic alterations after the inhibition of EZH2. We found that the inhibition of EZH2 causes transcriptional activation of numerous immune genes and inhibits the subsequent infection by influenza A virus. Interestingly, specific families of transposable elements, particularly endogenous retrovirus elements (ERVs) and LINEs which belong to retrotransposons, also become derepressed. While some of the derepressed ERV families are pig-specific, a few ancestral families are known to be under EZH2-mediated repression in humans. Given that derepression of ERVs can promote innate immune activation through "viral mimicry", we speculate that ERVs may also contribute to the coinciding immune activation in AMs after the inhibition of EZH2. Overall, this study improves the understanding of the EZH2-related immune regulation in AMs and provides novel insights into the epigenetic regulation of retrotransposons in pigs.

Isolation and Library Preparation of Planarian piRNAs

In planarian flatworms, piRNAs and SMEDWI (Schmidtea mediterranea PIWI) proteins are both essential for the animals' impressive regenerative ability and for their survival. A knockdown of SMEDWI proteins disrupts the specification of the planarian germline and impairs stem cell differentiation, resulting in lethal phenotypes. As the molecular targets of PIWI proteins and thus their biological function are determined by PIWI-bound small RNAs, termed piRNAs (for PIWI-interacting RNAs), it is imperative to study the wealth of PIWI-bound piRNAs using next-generation sequencing-based techniques. Prior to sequencing, piRNAs bound to individual SMEDWI proteins must be isolated. To that end, we established an immunoprecipitation protocol that can be applied to all planarian SMEDWI proteins. Co-immunoprecipitated piRNAs are visualized by using qualitative radioactive 5'-end labeling, which detects even trace amounts of small RNAs. Next, isolated piRNAs are subjected to a library preparation protocol that has been optimized for the efficient capture of piRNAs, whose 3'-ends carry a 2'-O-methyl modification. Successfully prepared piRNA libraries are subjected to Illumina-based next-generation sequencing. Obtained data are analyzed as presented in the accompanying manuscript.

Ubiquitination of stalled ribosomes enables mRNA decay via HBS-1 and NONU-1 in vivo

As ribosomes translate the genetic code, they can encounter a variety of obstacles that hinder their progress. If ribosomes stall for prolonged times, cells suffer due to the loss of translating ribosomes and the accumulation of aberrant protein products. Thus to protect cells, stalled ribosomes experience a series of reactions to relieve the stall and degrade the offending mRNA, a process known as No-Go mRNA Decay (NGD). While much of the machinery for NGD is known, the precise ordering of events and factors along this pathway has not been tested. Here, we deploy C. elegans to unravel the coordinated events comprising NGD. Utilizing a novel reporter and forward and reverse genetics, we identify the machinery required for NGD. Our subsequent molecular analyses define a functional requirement for ubiquitination on at least two ribosomal proteins (eS10 and uS10), and we show that ribosomes lacking ubiquitination sites on eS10 and uS10 fail to perform NGD in vivo. We show that the nuclease NONU-1 acts after the ubiquitin ligase ZNF-598, and discover a novel requirement for the ribosome rescue factors HBS-1/PELO-1 in mRNA decay via NONU-1. Taken together, our work demonstrates mechanisms by which ribosomes signal to effectors of mRNA repression, and we delineate links between repressive factors working toward a well-defined NGD pathway.

Genome-Wide Analysis of Planarian piRNAs

In planarian flatworms, the piRNA pathway is operated by three PIWI proteins, termed SMEDWI-1, SMEDWI-2, and SMEDWI-3 (SMEDWI = Schmidtea mediterranea PIWI). The interplay between these three PIWI proteins and their associated small noncoding RNAs, termed piRNAs, fuels the outstanding regenerative abilities of planarians, enables tissue homeostasis, and, ultimately, ensures animal survival. As the molecular targets of PIWI proteins are determined by the sequences of their co-bound piRNAs, it is imperative to identify these sequences by next-generation sequencing applications. Following sequencing, the genomic targets and the regulatory potential of the isolated piRNA populations need to be uncovered. To that end, here we present a bioinformatics analysis pipeline for processing and systematic characterization of planarian piRNAs. The pipeline includes steps for the removal of PCR duplicates based on unique molecular identifier (UMI) sequences, and it accounts for piRNA multimapping to different loci in the genome. Importantly, our protocol also includes a fully automated pipeline that is freely available at GitHub. Together with the piRNA isolation and library preparation protocol (see accompanying chapter), the presented computational pipeline enables researchers to explore the functional role of the piRNA pathway in flatworm biology.

A cap 0-dependent mRNA capture method to analyze the yeast transcriptome

Analysis of the protein coding transcriptome by the RNA sequencing requires either enrichment of the desired fraction of coding transcripts or depletion of the abundant non-coding fraction consisting mainly of rRNA. We propose an alternative mRNA enrichment strategy based on the RNA-binding properties of the human IFIT1, an antiviral protein recognizing cap 0 RNA. Here, we compare for Saccharomyces cerevisiae an IFIT1-based mRNA pull-down with yeast targeted rRNA depletion by the RiboMinus method. IFIT1-based RNA capture depletes rRNA more effectively, producing high quality RNA-seq data with an excellent coverage of the protein coding transcriptome, while depleting cap-less transcripts such as mitochondrial or some non-coding RNAs. We propose IFIT1 as a cost effective and versatile tool to prepare mRNA libraries for a variety of organisms with cap 0 mRNA ends, including diverse plants, fungi and eukaryotic microbes.

Ribosomal RNA-Depletion Provides an Efficient Method for Successful Dual RNA-Seq Expression Profiling of a Marine Sponge Holobiont

Investigations of host-symbiont interactions can benefit enormously from a complete and reliable holobiont gene expression profiling. The most efficient way to acquire holobiont transcriptomes is to perform RNA-Seq on both host and symbionts simultaneously. However, optimal methods for capturing both host and symbiont mRNAs are still under development, particularly when the host is a eukaryote and the symbionts are bacteria or archaea. Traditionally, poly(A)-enriched libraries have been used to capture eukaryotic mRNA, but the ability of this method to adequately capture bacterial mRNAs is unclear because of the short half-life of the bacterial transcripts. Here, we address this gap in knowledge with the aim of helping others to choose an appropriate RNA-Seq approach for analysis of animal host-bacterial symbiont transcriptomes. Specifically, we compared transcriptome bias, depth and coverage achieved by two different mRNA capture and sequencing strategies applied to the marine demosponge Amphimedon queenslandica holobiont. Annotated genomes of the sponge host and the three most abundant bacterial symbionts, which can comprise up to 95% of the adult microbiome, are available. Importantly, this allows for transcriptomes to be accurately mapped to these genomes, and thus quantitatively assessed and compared. The two strategies that we compare here are (i) poly(A) captured mRNA-Seq (Poly(A)-RNA-Seq) and (ii) ribosomal RNA depleted RNA-Seq (rRNA-depleted-RNA-Seq). For the host sponge, we find no significant difference in transcriptomes generated by the two different mRNA capture methods. However, for the symbiont transcriptomes, we confirm the expectation that the rRNA-depleted-RNA-Seq performs much better than the Poly(A)-RNA-Seq. This comparison demonstrates that RNA-Seq by ribosomal RNA depletion is an effective and reliable method to simultaneously capture gene expression in host and symbionts and thus to analyse holobiont transcriptomes.

Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

An Easy, Cost-Effective, and Scalable Method to Deplete Human Ribosomal RNA for RNA-seq

RNA sequencing (RNA-seq) is a powerful and increasingly prevalent method to characterize and quantify the transcriptome. Ribosomes are extremely abundant, however, and approximately 80% of total RNA is ribosomal RNA (rRNA). Therefore, to detect and quantify less abundant yet biologically important transcripts such as messenger RNA (mRNA) and long noncoding RNAs (lncRNA), it is essential to minimize the rRNA being sequenced. Although commercial methods exist to deplete rRNA from total RNA samples before sequencing, they are expensive and require specific amounts of input RNA, and the most commonly used kit is no longer available as a stand-alone product. Here, we present an optimized rRNA depletion protocol using RNase H and DNA oligonucleotides complementary to human rRNA transcripts. This protocol includes guidelines for DNA oligo preparation, RNA:DNA hybridization, RNase H cleavage and RNA cleanup, and benchmarking of rRNA depletion. The method is flexible because the user can include additional complementary DNA oligos directed against any abundant transcript in their particular system. Furthermore, the performance of this rRNA depletion approach is comparable to or better than that of commercial kits, at a fraction of the cost and across a wide range of input RNA amounts. © 2021 Wiley Periodicals LLC. Basic Protocol: Specific depletion of rRNA transcripts from human total RNA Support Protocol: Preparation of the rRNA depletion DNA oligonucleotide pool.

Sticky problems: extraction of nucleic acids from molluscs

Traditional molecular methods and omics-techniques across molluscan taxonomy increasingly inform biology of Mollusca. Recovery of DNA and RNA for such studies is challenged by common biological properties of the highly diverse molluscs. Molluscan biomineralization, adhesive structures and mucus involve polyphenolic proteins and mucopolysaccharides that hinder DNA extraction or copurify to inhibit enzyme-catalysed molecular procedures. DNA extraction methods that employ the detergent hexadecyltrimethylammoniumbromide (CTAB) to remove these contaminants importantly facilitate molecular-level study of molluscs. Molluscan pigments may stain DNA samples and interfere with spectrophotometry, necessitating gel electrophoresis or fluorometry for accurate quantification. RNA can reliably be extracted but the 'hidden break' in 28S rRNA of molluscs (like most protostomes) causes 18S and 28S rRNA fragments to co-migrate electrophoretically. This challenges the standard quality control based on the ratio of 18S and 28S rRNA, developed for deuterostome animals. High-AT content in molluscan rRNA prevents the effective purification of polyadenylated mRNA. Awareness of these matters aids the continuous expansion of molecular malacology, enabling work also with museum specimens and next-generation sequencing, with the latter imposing unprecedented demands on DNA quality. Alternative methods to extract nucleic acids from molluscs are available from literature and, importantly, from communications with others who study the molecular biology of molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

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.

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.

Overexpression of the vitamin D receptor (VDR) induces skeletal muscle hypertrophy

OBJECTIVE

The Vitamin D receptor (VDR) has been positively associated with skeletal muscle mass, function and regeneration. Mechanistic studies have focused on the loss of the receptor, with in vivo whole-body knockout models demonstrating reduced myofibre size and function and impaired muscle development. To understand the mechanistic role upregulation of the VDR elicits in muscle mass/health, we studied the impact of VDR over-expression (OE) in vivo before exploring the importance of VDR expression upon muscle hypertrophy in humans.

METHODS

Wistar rats underwent in vivo electrotransfer (IVE) to overexpress the VDR in the Tibialis anterior (TA) muscle for 10 days, before comprehensive physiological and metabolic profiling to characterise the influence of VDR-OE on muscle protein synthesis (MPS), anabolic signalling and satellite cell activity. Stable isotope tracer (D2O) techniques were used to assess sub-fraction protein synthesis, alongside RNA-Seq analysis. Finally, human participants underwent 20 wks of resistance exercise training, with body composition and transcriptomic analysis.

RESULTS

Muscle VDR-OE yielded total protein and RNA accretion, manifesting in increased myofibre area, i.e., hypertrophy. The observed increases in MPS were associated with enhanced anabolic signalling, reflecting translational efficiency (e.g., mammalian target of rapamycin (mTOR-signalling), with no effects upon protein breakdown markers being observed. Additionally, RNA-Seq illustrated marked extracellular matrix (ECM) remodelling, while satellite cell content, markers of proliferation and associated cell-cycled related gene-sets were upregulated. Finally, induction of VDR mRNA correlated with muscle hypertrophy in humans following long-term resistance exercise type training.

CONCLUSION

VDR-OE stimulates muscle hypertrophy ostensibly via heightened protein synthesis, translational efficiency, ribosomal expansion and upregulation of ECM remodelling-related gene-sets. Furthermore, VDR expression is a robust marker of the hypertrophic response to resistance exercise in humans. The VDR is a viable target of muscle maintenance through testable Vitamin D molecules, as active molecules and analogues.

Organism-specific depletion of highly abundant RNA species from bacterial total RNA

High-throughput sequencing has become a standard tool for transcriptome analysis. The depletion of overrepresented RNA species from sequencing libraries plays a key role in establishing potent and cost-efficient RNA-seq routines. Commercially available kits are known to obtain good results for the reduction of ribosomal RNA (rRNA). However, we found that the transfer-messenger RNA (tmRNA) was frequently highly abundant in rRNA-depleted samples of Pseudomonas aeruginosa, consuming up to 25 % of the obtained reads. The tmRNA fraction was particularly high in samples taken from stationary cultures. This suggests that overrepresentation of this RNA species reduces the mRNA fraction when cells are grown under challenging conditions. Here, we present an RNase-H-based depletion protocol that targets the tmRNA in addition to ribosomal RNAs. We were able to increase the mRNA fraction to 93–99% and therefore outperform not only the commercially Ribo-off kit (Vazyme) operating by the same principle but also the formerly widely used Ribo-Zero kit (Illumina). Maximizing the read share of scientifically interesting RNA species enhances the discriminatory potential of next-generation RNA-seq experiments and, therefore, can contribute to a better understanding of the transcriptomic landscape of bacterial pathogens and their used mechanisms in host infection.

Ribo-Pop: simple, cost-effective, and widely applicable ribosomal RNA depletion

The measurement of RNA abundance derived from massively parallel sequencing experiments is an essential technique. Methods that reduce ribosomal RNA levels are usually required prior to sequencing library construction because ribosomal RNA typically comprises the vast majority of a total RNA sample. For some experiments, ribosomal RNA depletion is favored over poly(A) selection because it offers a more inclusive representation of the transcriptome. However, methods to deplete ribosomal RNA are generally proprietary, complex, inefficient, applicable to only specific species, or compatible with only a narrow range of RNA input levels. Here, we describe Ribo-Pop (ribosomal RNA depletion for popular use), a simple workflow and antisense oligo design strategy that we demonstrate works over a wide input range and can be easily adapted to any organism with a sequenced genome. We provide a computational pipeline for probe selection, a streamlined 20-min protocol, and ready-to-use oligo sequences for several organisms. We anticipate that our simple and generalizable "open source" design strategy would enable virtually any laboratory to pursue full transcriptome sequencing in their organism of interest with minimal time and resource investment.

The piRNA pathway in planarian flatworms: new model, new insights

PIWI-interacting RNAs (piRNAs) are small regulatory RNAs that associate with members of the PIWI clade of the Argonaute superfamily of proteins. piRNAs are predominantly found in animal gonads. There they silence transposable elements (TEs), regulate gene expression and participate in DNA methylation, thus orchestrating proper germline development. Furthermore, PIWI proteins are also indispensable for the maintenance and differentiation capabilities of pluripotent stem cells in free-living invertebrate species with regenerative potential. Thus, PIWI proteins and piRNAs seem to constitute an essential molecular feature of somatic pluripotent stem cells and the germline. In keeping with this hypothesis, both PIWI proteins and piRNAs are enriched in neoblasts, the adult stem cells of planarian flatworms, and their presence is a prerequisite for the proper regeneration and perpetual tissue homeostasis of these animals. The piRNA pathway is required to maintain the unique biology of planarians because, in analogy to the animal germline, planarian piRNAs silence TEs and ensure stable genome inheritance. Moreover, planarian piRNAs also contribute to the degradation of numerous protein-coding transcripts, a function that may be critical for neoblast differentiation. This review gives an overview of the planarian piRNA pathway and of its crucial function in neoblast biology.

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.

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

Bacterial RNA sequencing (RNA-seq) is a powerful approach for quantitatively delineating the global transcriptional profiles of microbes in order to gain deeper understanding of their physiology and function. Cost-effective bacterial RNA-seq requires efficient physical removal of ribosomal RNA (rRNA), which otherwise dominates transcriptomic reads. However, current methods to effectively deplete rRNA of diverse non-model bacterial species are lacking. Here, we describe a probe and ribonuclease based strategy for bacterial rRNA removal. We implemented the method using either chemically synthesized oligonucleotides or amplicon-based single-stranded DNA probes and validated the technique on three novel gut microbiota isolates from three distinct phyla. We further showed that different probe sets can be used on closely related species. We provide a detailed methods protocol, probe sets for >5000 common microbes from RefSeq, and an online tool to generate custom probe libraries. This approach lays the groundwork for large-scale and cost-effective bacterial transcriptomics studies.