Reverse complementary matches simultaneously promote both back-splicing and exon-skipping

Background Circular RNAs (circRNAs) play diverse roles in different biological and physiological environments and are always expressed in a tissue-specific manner. Especially, circRNAs are enriched in the brain tissues of almost all investigated species, including humans, mice, Drosophila, etc. Although circRNAs were found in C. elegans, the neuron-specific circRNA data is not available yet. Exon-skipping is found to be correlated to circRNA formation, but the mechanisms that link them together are not clear. Results Here, through large-scale neuron isolation from the first larval (L1) stage of C. elegans followed by RNA sequencing with ribosomal RNA depletion, the neuronal circRNA data in C. elegans were obtained. Hundreds of novel circRNAs were annotated with high accuracy. circRNAs were highly expressed in the neurons of C. elegans and were positively correlated to the levels of their cognate linear mRNAs. Disruption of reverse complementary match (RCM) sequences in circRNA flanking introns effectively abolished circRNA formation. In the zip-2 gene, deletion of either upstream or downstream RCMs almost eliminated the production of both the circular and the skipped transcript. Interestingly, the 13-nt RCM in zip-2 is highly conserved across five nematode ortholog genes, which show conserved exon-skipping patterns. Finally, through in vivo one-by-one mutagenesis of all the splicing sites and branch points required for exon-skipping and back-splicing in the zip-2 gene, I showed that back-splicing still happened without exon-skipping, and vice versa. Conclusions Through protocol optimization, total RNA obtained from sorted neurons is increased to hundreds of nanograms. circRNAs highly expressed in the neurons of C. elegans are more likely to be derived from genes also highly expressed in the neurons. RCMs are abundant in circRNA flanking introns, and RCM-deletion is an efficient way to knockout circRNAs. More importantly, these RCMs are not only required for back-splicing but also promote the skipping of exon(s) to be circularized. Finally, RCMs in circRNA flanking introns can directly promote both exon-skipping and back-splicing, providing a new explanation for the correlation between them.

Dual RNA-seq analysis of in vitro infection multiplicity and RNA depletion methods in Chlamydia-infected epithelial cells

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host–pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

SPAAC-NAD-seq, a sensitive and accurate method to profile NAD+-capped transcripts

Nicotinamide adenine diphosphate (NAD+) is a novel messenger RNA 5′ cap in Escherichia coli, yeast, mammals, and Arabidopsis. Transcriptome-wide identification of NAD+-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide–alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m7G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m7G-capped RNAs (m7G-RNAs). In addition, the Cu2+ ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide–alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m7G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m7G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis. Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes.

Effective ribosomal RNA depletion for single-cell total RNA-seq by scDASH

A decade since its invention, single-cell RNA sequencing (scRNA-seq) has become a mainstay technology for profiling transcriptional heterogeneity in individual cells. Yet, most existing scRNA-seq methods capture only polyadenylated mRNA to avoid the cost of sequencing non-messenger transcripts, such as ribosomal RNA (rRNA), that are usually not of-interest. Hence, there are not very many protocols that enable single-cell analysis of total RNA. We adapted a method called DASH (Depletion of Abundant Sequences by Hybridisation) to make it suitable for depleting rRNA sequences from single-cell total RNA-seq libraries. Our analyses show that our single-cell DASH (scDASH) method can effectively deplete rRNAs from sequencing libraries with minimal off-target non-specificity. Importantly, as a result of depleting the rRNA, the rest of the transcriptome is significantly enriched for detection.

Paired rRNA-depleted and polyA-selected RNA sequencing data and supporting multi-omics data from human T cells

Both poly(A) enrichment and ribosomal RNA depletion are commonly used for RNA sequencing. Either has its advantages and disadvantages that may lead to biases in the downstream analyses. To better access these effects, we carried out both ribosomal RNA-depleted and poly(A)-selected RNA-seq for CD4+ T naive cells isolated from 40 healthy individuals from the Blueprint Project. For these 40 individuals, the genomic and epigenetic data were also available. This dataset offers a unique opportunity to understand how library construction influences differential gene expression, alternative splicing and molecular QTL (quantitative loci) analyses for human primary cells.

Dual RNA-Seq analysis of in vitro infection multiplicity in Chlamydia-infected epithelial cells

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 hours), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

Architectural RNA in chromatin organization

RNA plays a well-established architectural role in the formation of membraneless interchromatin nuclear bodies. However, a less well-known role of RNA is in organizing chromatin, whereby specific RNAs have been found to recruit chromatin modifier proteins. Whether or not RNA can act as an architectural molecule for chromatin remains unclear, partly because dissecting the architectural role of RNA from its regulatory role remains challenging. Studies that have addressed RNA's architectural role in chromatin organization rely on in situ RNA depletion using Ribonuclease A (RNase A) and suggest that RNA plays a major direct architectural role in chromatin organization. In this review, we will discuss these findings, candidate chromatin architectural long non-coding RNAs and possible mechanisms by which RNA, along with RNA binding proteins might be mediating chromatin organization.