Search

Abstract Background Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts’ direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique. Results A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes. Conclusion The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs.

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cellular diversity with unprecedented resolution. However, many current methods are limited in capturing full-length transcripts and discerning strand orientation. We present RAG-seq, an innovative strand-specific total RNA sequencing technique that combines not-so-random (NSR) primers with Tn5 transposase-mediated tagmentation. RAG-seq overcomes previous limitations by delivering comprehensive transcript coverage and maintaining strand orientation, which is essential for accurate quantification of overlapping genes and detection of antisense transcripts. Through optimized reverse transcription with oligo dT primers, rRNA depletion via Depletion of Abundant Sequences by Hybridization (DASH), and linear amplification, RAG-seq enhances sensitivity and reproducibility, especially for low-input samples and single cells. Application to mouse oocytes and early embryos highlights RAG-seq's superior performance in identifying stage-specific antisense transcripts, shedding light on their regulatory roles during early development. This advancement represents a significant leap in transcriptome analysis within complex biological contexts., (© The Author(s) 2024. Published by Oxford University Press and Science Press on behalf of the Beijing Institute of Genomics, Chinese Academy of Sciences / China National Center for Bioinformation and Genetics Society of China.)

Background: Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts' direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique. Results: A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes. Conclusion: The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs. [ABSTRACT FROM AUTHOR]

R-loops are three-stranded nucleic acid structures with essential roles in many nuclear processes. However, their unchecked accumulation is associated with genome instability and is observed in neurodevelopmental diseases and cancers. Genome-wide profiling of R-loops in normal and diseased cells can help identify locations of pathogenic R-loops and advance efforts to attenuate them. We present an antibody-independent R-loop detection strategy, BisMapR, that combines nuclease-based R-loop isolation with non-denaturing bisulfite chemistry to produce genome-wide profiles that retain strand information. BisMapR achieves greater resolution and is faster than existing strand-specific R-loop profiling strategies. In mouse embryonic stem cells, we apply BisMapR to find that gene promoters form R-loops in both directions and uncover a subset of active enhancers that, despite being bidirectionally transcribed, form R-loops exclusively on one strand. BisMapR reveals a previously unnoticed feature of active enhancers and provides a tool to systematically examine their mechanisms in gene expression.

Abstract Background RNA-Sequencing (RNA-seq) is now commonly used to reveal quantitative spatiotemporal snapshots of the transcriptome, the structures of transcripts (splice variants and fusions) and landscapes of expressed mutations. However, standard approaches for library construction typically require relatively high amounts of input RNA, are labor intensive, and are time consuming. Methods Here, we report the outcome of a systematic effort to optimize and streamline steps in strand-specific RNA-seq library construction. Results This work has resulted in the identification of an optimized messenger RNA isolation protocol, a potent reverse transcriptase for cDNA synthesis, and an efficient chemistry and a simplified formulation of library construction reagents. We also present an optimization of bead-based purification and size selection designed to maximize the recovery of cDNA fragments. Conclusions These developments have allowed us to assemble a rapid high throughput pipeline that produces high quality data from amounts of total RNA as low as 25 ng. While the focus of this study is on RNA-seq sample preparation, some of these developments are also relevant to other next-generation sequencing library types.

Abstract Background RNA-Seq is now widely used as a research tool. Choices must be made whether to use paired-end (PE) or single-end (SE) sequencing, and whether to use strand-specific or non-specific (NS) library preparation kits. To date there has been no analysis of the effect of these choices on identifying differentially expressed genes (DEGs) between controls and treated samples and on downstream functional analysis. Results We undertook four mammalian transcriptomics experiments to compare the effect of SE and PE protocols on read mapping, feature counting, identification of DEGs and functional analysis. For three of these experiments we also compared a non-stranded (NS) and a strand-specific approach to mapping the paired-end data. SE mapping resulted in a reduced number of reads mapped to features, in all four experiments, and lower read count per gene. Up to 4.3% of genes in the SE data and up to 12.3% of genes in the NS data had read counts which were significantly different compared to the PE data. Comparison of DEGs showed the presence of false positives (average 5%, using voom) and false negatives (average 5%, using voom) using the SE reads. These increased further, by one or two percentage points, with the NS data. Gene ontology functional enrichment (GO) of the DEGs arising from SE or NS approaches, revealed striking differences in the top 20 GO terms, with as little as 40% concordance with PE results. Caution is therefore advised in the interpretation of such results. By comparison, there was overall consistency in gene set enrichment analysis results. Conclusions A strand-specific protocol should be used in library preparation to generate the most reliable and accurate profile of expression. Ideally PE reads are also recommended particularly for transcriptome assembly. Whilst SE reads produce a DEG list with around 5% of false positives and false negatives, this method can substantially reduce sequencing cost and this saving could be used to increase the number of biological replicates thereby increasing the power of the experiment. As SE reads, when used in association with gene set enrichment, can generate accurate biological results, this may be a desirable trade-off.

Abstract To distinguish between three types of Siniperca chuatsi rhabdovirus (SCRV) viral RNA (vRNA, cRNA, and mRNA) and investigate SCRV transcription and replication dynamics in Chinese perch brain CPB cells, a novel, strand-specific, reverse transcriptase quantitative real-time PCR (RT-qPCR) assay was established. The method is based on strand-specific reverse transcription, using tagged primers to add a 'tag' sequence at the 5′ end. We used the 'tag' sequence as the forward primer and a strand-specific reverse primer to quantify the three types of RNA. Three types of synthetic viral RNA were used as reference standards for validation and quantification. These assays were optimized to produce a standard curve from 102 to 107 copies/μL, with an efficiency of 91–101% and an R2 value of 0.9949–0.9999. The coefficients of variation for repeatability and reproducibility were less than 2.85% and 5.52%, respectively. Using this method, specific target RNA was detected at a 3500–70,000 fold higher level than other types of RNA. This method was also used to evaluate the dynamics of vRNA, cRNA and mRNA synthesis in CPB cells infected with SCRV. The results indicate that the intracellular dynamics of vRNA, cRNA and mRNA are different. In the earliest phase of SCRV infection, all three types of viral RNA increased very slowly. The copy number of vRNA and mRNA increased exponentially from 4 h post infection, while cRNA increased from 6 h post infection. The amount of cRNA was lower than vRNA and mRNA throughout the infection. The novel, strand-specific RT-qPCR method developed in this study provides critical data to aid the understanding of transcription and replication during SCRV infection. Highlights • A RT-qPCR for distinguish between three types of SCRV viral RNA was developed. • The intracellular dynamics of vRNA, cRNA and mRNA are different. • The amount of cRNA was lower than vRNA and mRNA throughout the SCRV infection. [ABSTRACT FROM AUTHOR]

Background: RNA-Seq is now widely used as a research tool. Choices must be made whether to use paired-end (PE) or single-end (SE) sequencing, and whether to use strand-specific or non-specific (NS) library preparation kits. To date there has been no analysis of the effect of these choices on identifying differentially expressed genes (DEGs) between controls and treated samples and on downstream functional analysis. Results: We undertook four mammalian transcriptomics experiments to compare the effect of SE and PE protocols on read mapping, feature counting, identification of DEGs and functional analysis. For three of these experiments we also compared a non-stranded (NS) and a strand-specific approach to mapping the paired-end data. SE mapping resulted in a reduced number of reads mapped to features, in all four experiments, and lower read count per gene. Up to 4.3% of genes in the SE data and up to 12.3% of genes in the NS data had read counts which were significantly different compared to the PE data. Comparison of DEGs showed the presence of false positives (average 5%, using voom) and false negatives (average 5%, using voom) using the SE reads. These increased further, by one or two percentage points, with the NS data. Gene ontology functional enrichment (GO) of the DEGs arising from SE or NS approaches, revealed striking differences in the top 20 GO terms, with as little as 40% concordance with PE results. Caution is therefore advised in the interpretation of such results. By comparison, there was overall consistency in gene set enrichment analysis results. Conclusions: A strand-specific protocol should be used in library preparation to generate the most reliable and accurate profile of expression. Ideally PE reads are also recommended particularly for transcriptome assembly. Whilst SE reads produce a DEG list with around 5% of false positives and false negatives, this method can substantially reduce sequencing cost and this saving could be used to increase the number of biological replicates thereby increasing the power of the experiment. As SE reads, when used in association with gene set enrichment, can generate accurate biological results, this may be a desirable trade-off. [ABSTRACT FROM AUTHOR]

Background: The development of multicellular organisms is coordinated by various gene regulatory mechanisms that ensure correct spatio-temporal patterns of gene expression. Recently, the role of antisense transcription in gene regulation has moved into focus of research. To characterize genome-wide patterns of antisense transcription and to study their evolutionary conservation, we sequenced a strand-specific RNA-seq library of the nematode Pristionchus pacificus. Results: We identified 1112 antisense configurations of which the largest group represents 465 antisense transcripts (ASTs) that are fully embedded in introns of their host genes. We find that most ASTs show homology to protein-coding genes and are overrepresented in proteomic data. Together with the finding, that expression levels of ASTs and host genes are uncorrelated, this indicates that most ASTs in P. pacificus do not represent non-coding RNAs and do not exhibit regulatory functions on their host genes. We studied the evolution of antisense gene pairs across 20 nematode genomes, showing that the majority of pairs is lineage-specific and even the highly conserved vps-4, ddx-27, and sel-2 loci show abundant structural changes including duplications, deletions, intron gains and loss of antisense transcription. In contrast, host genes in general, are remarkably conserved and encode exceptionally long introns leading to unusually large blocks of conserved synteny. Conclusions: Our study has shown that in P. pacificus antisense transcription as such does not define non-coding RNAs but is rather a feature of highly conserved genes with long introns. We hypothesize that the presence of regulatory elements imposes evolutionary constraint on the intron length, but simultaneously, their large size makes them a likely target for translocation of genomic elements including protein-coding genes that eventually end up as ASTs. [ABSTRACT FROM AUTHOR]

Reverse-transcription quantitative polymerase chain reaction (RT-qPCR) is widely used to quantify viral RNA genomes for diagnostics and research, yet conventional RT-qPCR protocols are unable to accurately distinguish between the different viral RNA species that exist during infection. Here we show that false-priming and self-priming occur during reverse transcription with several published Zika virus (ZIKV) primer sets. We developed a RT-qPCR assay using tagged primers and thermostable reverse transcriptase, which greatly reduced the occurrence of nonspecific cDNA products. Furthermore, we optimized the assay for use in multiplex qPCR which allows for simultaneous quantitative detection of positive-strand and negative-strand ZIKV RNA along with an internal control from both human and mosquito cells. Importantly, this assay is sensitive enough to study early stages of virus infection in vitro. Strikingly, using this assay, we detected ZIKV negative-strand RNA as early as 3 h post-infection in mammalian cell culture, at a time point prior to the onset of positive-strand RNA synthesis. Overall, the strand-specific RT-qPCR assay developed herein is a valuable tool to quantify ZIKV RNA and to study viral replication dynamics during infection. The application of these findings has the potential to increase accuracy of RNA detection methods for a variety of viral pathogens. [Display omitted] • Self-primed cDNA is amplified by widely-used ZIKV qPCR primer sets. • Use of tagged primers and thermostable RT increases strand-specificity for RT-qPCR. • Multiplexed qPCR allows simultaneous quantitation of (+) and (-) strand viral RNAs. • Strand-specific RT-qPCR can detect fewer than one copy of viral RNA per cell. [ABSTRACT FROM AUTHOR]

Abstract: Murine norovirus (MNV), currently the only norovirus that efficiently replicates in cell culture, is often used as a model system to understand the molecular mechanisms of norovirus replication. MNV is a single stranded positive sense RNA virus of the Caliciviridae family. Replication of MNV involves the synthesis of both full length genomic and sub-genomic RNAs. The replication of these RNAs involves the synthesis of negative strand intermediates. To understand the molecular mechanism of RNA replication and the role of viral and host factors in virus replication, it is necessary to quantify accurately both positive and negative sense RNA molecules of the viral RNA during replication. Increasingly, strand specific reverse transcription-quantitative PCR (RT-qPCR) is becoming the method of choice for this kind of quantitation. Many strategies have been developed to avoid the false priming property of reverse transcriptase and to amplify specifically one strand in the presence of excess opposite strand. In this report, a SYBR based, real time RT-qPCR assay was developed to detect and quantify specifically the negative and the positive sense RNAs of MNV genomic RNA. This assay is based on using a tagged RT primer containing a non-viral sequence at the 5′ end of the viral strand specific sequence. This non-viral sequence is then used to amplify selectively the strand specific cDNA at the PCR stage. This assay can be used for a range of MNV strains including MNV-1 and 3, as these are now widely accepted for use in molecular studies. The specificity of this assay was determined by its ability to quantify one strand in the presence of up to 106 copies of competitor opposite sense RNA. Using this assay, the production of both strands of MNV-1 RNA was monitored during viral single step growth curve. [Copyright &y& Elsevier]

Abstract: The detection of replicative intermediate RNAs as markers of active replication of RNA viruses is an essential tool to investigate pathogenesis in acute viral infections, as well as in their long-term sequelae. In this regard, strand-specific PCR has been used widely to distinguish (−) and (+) enteroviral RNAs in pathogenesis studies of diseases such as dilated cardiomyopathy. It has been generally assumed that oligonucleotide-primed reverse transcription of a given RNA generates only the corresponding specific cDNA, thus assuring the specificity of a PCR product amplified from it. Nevertheless, such assumed strand-specificity is a fallacy, because falsely primed cDNAs can be produced by RNA reverse transcription in the absence of exogenously added primers, (cDNAprimer(−)), and such falsely primed cDNAs are amplifiable by PCR in the same way as the correctly primed cDNAs. Using as a prototype the coxsackievirus B5 (CVB5), a (+) strand RNA virus, it was shown that cDNAprimer(−) renders the differential detection of viral (−) and (+) RNAs by conventional PCR virtually impossible, due to gross non-specificity. Using in vitro transcribed CVB5 RNAs (+) and (−), it was shown that cDNAprimer(−) could be removed effectively by magnetic physical separation of correctly primed biotinylated cDNA. Such strategy enabled truly strand-specific detection of RNA (−) and (+), not only for CVB5, but also for other non-polio enteroviruses. These findings indicate that previous conclusions supporting a role for the persistence of actively replicating enterovirus in the pathogenesis of chronic myocarditis should be regarded with strong skepticism and purification of correctly primed cDNA should be used for strand-specific PCR of viral RNA in order to obtain reliable information on this important subject. [Copyright &y& Elsevier]

Abstract: The fish rhabdovirus, Infectious hematopoietic necrosis virus (IHNV), is an important pathogen of salmonids. Cell culture assays have traditionally been used to quantify levels of IHNV in samples; however, real-time or quantitative RT-PCR assays have been proposed as a rapid alternative. For viruses having a single-stranded, negative-sense RNA genome, standard qRT-PCR assays do not distinguish between the negative-sense genome and positive-sense RNA species including mRNA and anti-genome. Thus, these methods do not determine viral genome copy number. This study reports development of strand-specific, qRT-PCR assays that use tagged primers for enhancing strand specificity during cDNA synthesis and quantitative PCR. Protocols were developed for positive-strand specific (pss-qRT-PCR) and negative-strand specific (nss-qRT-PCR) assays for IHNV glycoprotein (G) gene sequences. Validation with synthetic RNA transcripts demonstrated the assays could discriminate the correct strand with greater than 1000-fold fidelity. The number of genome copies in livers of IHNV-infected fish determined by nss-qRT-PCR was, on average, 8000-fold greater than the number of infectious units as determined by plaque assay. We also compared the number of genome copies with the quantity of positive-sense RNA and determined that the ratio of positive-sense molecules to negative-sense genome copies was, on average, 2.7:1. Potential future applications of these IHNV strand-specific qRT-PCR assays are discussed. [Copyright &y& Elsevier]

R-loops are three-stranded, DNA:RNA hybrid-containing structures that form naturally throughout the genome as a consequence of transcription. Accurately determining the genomic locations and strand of origin of R-loops is critical to understanding their roles in gene regulation and disease. Here, we describe a nuclease-based protocol for genome-wide and strand-specific R-loop detection, which we term MapR. This method targets native R-loops for cleavage and release using a modified RNase H enzyme, followed by deep sequencing. An extension of the protocol, BisMapR, can additionally introduce strand specificity via non-denaturing bisulfite conversion of the R-loop's single-stranded DNA component. MapR and BisMapR identify R-loops with high resolution and low background, can be performed with low cell input, and require short experimental time., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

CLIP-Seq (Deep Sequencing after in vivo Crosslinking and Immunoprecipitation, HITS-CLIP) has emerged as a key method for the study of RNA-binding proteins (RBPs), as it can scrutinize the RNAs bound by an RBP in vivo, with minimum manipulation of biological samples. CLIP-Seq is best used to reveal changes of the RNA cargo of an RBP and differences on binding patterns of the bound RNAs in living cells in different genetic backgrounds or after experimental treatment, rather than simply identifying RNA species. It is therefore crucial that a reference of the steady state levels of the RNAs present in the samples used for the CLIP-Seq experiment is included in the bioinformatic analysis. A simple directional RNA-Seq method was developed that uses the same oligonucleotides and the same PCR amplification steps as our CLIP-Seq method, which therefore can be analyzed using the same bioinformatic pipeline as the CLIP-Seq data. This greatly simplifies and streamlines the analysis process, and at the same time reduces the chances of protocol-specific artifacts and biases interfering with data interpretation. Some considerations on ways to integrate CLIP-Seq and RNA-Seq analyses are also provided herein., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

R-loops are three-stranded nucleic acid structures with essential roles in many nuclear processes. However, their unchecked accumulation is associated with genome instability and is observed in neurodevelopmental diseases and cancers. Genome-wide profiling of R-loops in normal and diseased cells can help identify locations of pathogenic R-loops and advance efforts to attenuate them. We present an antibody-independent R-loop detection strategy, BisMapR, that combines nuclease-based R-loop isolation with non-denaturing bisulfite chemistry to produce genome-wide profiles that retain strand information. BisMapR achieves greater resolution and is faster than existing strand-specific R-loop profiling strategies. In mouse embryonic stem cells, we apply BisMapR to find that gene promoters form R-loops in both directions and uncover a subset of active enhancers that, despite being bidirectionally transcribed, form R-loops exclusively on one strand. BisMapR reveals a previously unnoticed feature of active enhancers and provides a tool to systematically examine their mechanisms in gene expression., Competing Interests: PW No competing interests declared, KS Reviewing editor, eLife, (© 2021, Wulfridge and Sarma.)

Background: RNA-Sequencing (RNA-seq) is now commonly used to reveal quantitative spatiotemporal snapshots of the transcriptome, the structures of transcripts (splice variants and fusions) and landscapes of expressed mutations. However, standard approaches for library construction typically require relatively high amounts of input RNA, are labor intensive, and are time consuming. METHODS: Here, we report the outcome of a systematic effort to optimize and streamline steps in strand-specific RNA-seq library construction. RESULTS: This work has resulted in the identification of an optimized messenger RNA isolation protocol, a potent reverse transcriptase for cDNA synthesis, and an efficient chemistry and a simplified formulation of library construction reagents. We also present an optimization of bead-based purification and size selection designed to maximize the recovery of cDNA fragments., Conclusions: These developments have allowed us to assemble a rapid high throughput pipeline that produces high quality data from amounts of total RNA as low as 25 ng. While the focus of this study is on RNA-seq sample preparation, some of these developments are also relevant to other next-generation sequencing library types.

This unit presents protocols for construction of next-generation sequencing (NGS) directional RNA sequencing libraries for the Illumina HiSeq and MiSeq from a wide variety of input RNA sources. The protocols are based on the New England Biolabs (NEB) small RNA library preparation set for Illumina, although similar kits exist from different vendors. The protocol preserves the orientation of the original RNA in the final sequencing library, enabling strand-specific analysis of the resulting data. These libraries have been used for differential gene expression analysis and small RNA discovery and are currently being tested for de novo transcriptome assembly. The protocol is robust and applicable to a broad range of RNA input types and RNA quality, making it ideal for high-throughput laboratories., (Copyright © 2014 John Wiley & Sons, Inc.)