poly-A Selection Research Articles

Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion.

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Analysis Of Differentially Expressed Genes By RNA-Sequencing In 38 Cases of Rasmussen Encephalitis: New Insights on Non-Inflammatory Pathways (P12-9.002)

<h3>Objective:</h3> To analyze aberrant gene expression pathways involved in the pathophysiology of Rasmussen Encephalitis (RE) in a large cohort of brain specimens from affected children. <h3>Background:</h3> RE is a rare cause of pharmacoresistant epilepsy, characterized by progressive hemispheric inflammation resulting in seizures and, if left untreated, hemiplegia, hemianopia and cognitive decline. Data on gene expression in this condition are lacking. <h3>Design/Methods:</h3> We conducted RNA-sequencing of flash-frozen resected epileptogenic tissue from 38 children with RE. Control brain tissue included 10 brain specimens from perilesional resected areas in different patients with focal cortical dysplasia. RNA-sequencing was performed using standard protocol on Illumina HiSeq 2x150bp sequencing, single index platform with polyA selection library prep. After extraction of normalized gene hit counts, we utilized the <i>DESeq2</i> package with R version 4.2.1. The Wald test was used to generate p-values and log2 fold changes. Genes with an adjusted p-value &lt; 0.05 were called as differentially expressed genes (DEG) for each comparison. Gene enrichment analysis was then performed with the Reactome platform. <h3>Results:</h3> We identified 6184 upregulated and 3789 downregulated genes in RE vs. controls. Functional gene clustering analysis showed overrepresentation of several pathways involved in neuroinflammation including interferon alpha/beta (p=7.77E-11) and gamma (p=1.17E-12) signaling, IL-4/IL-13 (p=1.15E-2), IL-10 (p=2.18E-2) signaling, class I MHC-mediated antigen processing (p=4.22E-5), and non-inflammatory pathways such as nonsense-mediated RNA decay (p=2.04E-3), translation elongation cycle (p=5.42E-4). We observed downregulation of multiple neuronal pathways: neurotransmitter receptors and postsynaptic signal transmission (p=3.79E-7), activation of NMDA receptors (p=1.72E-5), GABA receptor activation (p=3.19E-2), glutamate binding (p=2.1E-2), presynaptic depolarization and Ca channel opening (p=1.49E-2). Extra-neuronal downregulated genes included: transcriptional regulation of MECP2 (p=4.49E-4), L1CAM interactions (p=1.51E-3), signaling by NTRKs (p=3.93E-2), calmodulin pathway (p=2.05E-1). <h3>Conclusions:</h3> With this study, we contribute to expanding the knowledge of pathophysiological inflammatory and non-inflammatory processes involved in Rasmussen’s encephalitis. <b>Disclosure:</b> Dr. Bartolini has received personal compensation in the range of $10,000-$49,999 for serving as an Editor, Associate Editor, or Editorial Advisory Board Member for American Academy of Neurology. The institution of Dr. Chan has received research support from TSC Alliance. Ms. Lob has nothing to disclose. Ms. Wu has nothing to disclose. Dr. Uzun has nothing to disclose. Dr. Oluigbo has nothing to disclose. The institution of Dr. Gaillard has received research support from all federal or foundation grants, NINDS, NIDCD, NICHD, NSF ,PERF. Dr. Chang has nothing to disclose. Gary W. Mathern has nothing to disclose. The institution of Dr. Liu has received research support from NIH/NINDS.

Abstract 1332: Extracted mRNA from frozen buffy coat samples shows promise for use in downstream sequencing analyses

Abstract Frozen blood samples are routinely collected and stored in medical biobanks across the world, often for periods of several decades, and later used for medical research such as biomarker studies. Some common biomarker targets have included DNA alterations, protein or antibody levels or metabolomics. However, despite the potential importance of gene-expression alterations in circulating immune cells, transcriptional studies on blood samples from biobanks have been limited due to preanalytical technical difficulties of extracting high-quality RNA. In the absence of an RNA preservative, RNA degrades rapidly during the thawing process. Our aim was, therefore, to investigate the possibility of using biobanked blood samples for transcriptomics analyses. First, we used a CryoXtract CXT 350 to establish a protocol capable of isolating high-quality RNA from frozen core aliquots of buffy coat without thawing the entire sample tube. Then, we performed RNA sequencing on RNA extracted from 20 pre-diagnostic buffy coat samples from colorectal cancer patients and matched controls included in the population-based, prospective Northern Sweden Health and Disease Study (NSHDS). For all samples, two different library preparation methods were used (Illumina TruSeq Stranded mRNA poly-A selection and Illumina Stranded Total RNA with Ribo-Zero Globin). We obtained RNA yields of more than 1 µg from a majority of the NSHDS samples (mean=2.57 µg), and more than 92% of samples had RIN values of 6 or higher, indicating suitability for downstream analyses. Preliminary random forest analyses showed expected discriminatory ability between men and women based on data generated from each of the two library preparation methods. In conclusion, our newly established RNA extraction method enables extraction of RNA from frozen, biobanked buffy coat samples, collected in standard tubes, that is of sufficient quantity and quality for use in downstream whole-genome transcriptomic analyses, such as RNA sequencing. Citation Format: Erik Bovinder Ylitalo, Linda Vidman, Sophia Harlid, Bethany Van Guelpen. Extracted mRNA from frozen buffy coat samples shows promise for use in downstream sequencing analyses [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1332.

Intra-Exon Motif Correlations as a Proxy Measure for Mean Per-Tile Sequence Quality Data in RNA-Seq.

Given the wide variability in the quality of next-generation sequencing data submitted to public repositories, it is essential to identify methods that can perform quality control on these data sets when additional quality control data, such as mean tile data, are missing from public repositories. In this study, we present evidence that correlating counts of reads corresponding to pairs of motifs separated over specific distances on individual exons can be used as a proxy mean tile data in the data sets we analyzed and hence could be used when mean tile data are not available. As test data sets we use the Homo sapiens in vitro transcribed (IVT) data set, and a Drosophila melanogaster data set comprising wild and mutant types. We find that a FastQC analysis of the available parts of these data sets demonstrates that the per-tile sequencing quality is good for all the data sets apart from the mutant-type data where the mutant-r3 data are worse than the mutant-r2 data. Correspondingly, intra-exon motif correlations are reasonably large for all data sets except this latter case where the mutant-r2 correlations are low and the mutant-r3 correlations close to zero. We propose that these extremely low correlations are indicative of bias of technical origin, such as flowcell errors. In addition to this, the intra-exon motif correlations as a function of both guanosine-cytosine (GC) content parameters are somewhat higher and less dependent on the GC content parameters in the IVT-Plasmids messenger RNA (mRNA) selection free RNA-Seq sample (control) than in the other RNA-Seq samples that did undergo mRNA selection: both ribosomal depletion (IVT-Only) and PolyA selection (IVT-PolyA, wild type, and mutant).

Abstract LB059: Subtype classification of pediatric high-grade glioma tumors by comparative transcriptomics

Abstract Pediatric high-grade glioma (pHGG) is a highly malignant and poorly understood cancer driven by diverse genetic and epigenetic mechanisms. Here, we use comparative RNA sequencing, outlier analysis, and spectral clustering approaches to analyze transcriptomic data of 1,543 pediatric brain tumor specimens from the UCSC Treehouse Childhood Cancer Initiative (Treehouse) and Open Pediatric Brain Tumor Atlas (OpenPBTA) to identify subpopulations of pHGG patients with characteristic gene expression profiles. We find that approximately half (45%) of pHGG tumors from OpenPBTA exist in three subgroups defined by high outlier-level expression either of: mitochondrially-encoded 12S and 16S rRNAs; genes enriched in the HSF1-mediated heat shock response and activation pathways; or six C/D box snoRNA (SNORD) genes originating from the paternally-expressed SNORD116 locus involved in Prader-Willi syndrome, a complex neurodevelopmental disorder. Interestingly, the same set of HSF1-dependent pathway genes is also significantly upregulated in a subset (~11%) of pHGG tumors from Treehouse, validating this finding in two independent compendia with different transcript isolation strategies (Treehouse, polyA selection; OpenPBTA, ribodepletion). Our work identifies distinct classes of tumors with outlier-level expression of genes with previously unknown roles in pHGG and provides a framework for subtyping tumors by comparative transcriptomics that is adaptable to any cancer type. We are currently investigating the molecular roles of HSF1-response genes and the imprinted SNORD116 gene cluster in pHGG. Our ongoing research into the biomolecular signatures and mechanisms of the three major tumor classes of pHGG as defined in our study will contribute to a greater understanding of pHGG disease manifestation and progression, and will inform strategies of tailored therapeutic interventions for children with this devastating disease. Citation Format: Gina D. Mawla, A. Geoffrey Lyle, Ellen T. Kephart, Katrina Learned, Holly C. Beale, Joshua E. Goldford, Olena M. Vaske. Subtype classification of pediatric high-grade glioma tumors by comparative transcriptomics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB059.

Skeletal muscle transcriptomics identifies common pathways in nerve crush injury and ageing

Motor unit remodelling involving repeated denervation and re-innervation occurs throughout life. The efficiency of this process declines with age contributing to neuromuscular deficits. This study investigated differentially expressed genes (DEG) in muscle following peroneal nerve crush to model motor unit remodelling in C57BL/6 J mice. Muscle RNA was isolated at 3 days post-crush, RNA libraries were generated using poly-A selection, sequenced and analysed using gene ontology and pathway tools. Three hundred thirty-four DEG were found in quiescent muscle from (26mnth) old compared with (4-6mnth) adult mice and these same DEG were present in muscle from adult mice following nerve crush. Peroneal crush induced 7133 DEG in muscles of adult and 699 DEG in muscles from old mice, although only one DEG (ZCCHC17) was found when directly comparing nerve-crushed muscles from old and adult mice. This analysis revealed key differences in muscle responses which may underlie the diminished ability of old mice to repair following nerve injury.

Abstract 2578: Calcium modulates apoptosis and differentiation in human normal colon 3D organoids irrespective of donor sex and biopsy site

Abstract Background: Observational studies indicate that calcium, a micronutrient involved in numerous cell-signaling cascades, may protect against colon cancer. Stratified analyses suggest the protective effect of calcium may depend on sex and anatomic subsite, but these hypotheses have been difficult to test experimentally. Here, we employ human-derived normal colon 3D organoids as a model of the healthy colon crypt to interrogate the effect of calcium exposure on gene expression and the effect of donor sex and biopsy anatomic subsite on the response to calcium. Methods: A total of 39 3D colon organoid lines derived from colon biopsies (22 proximal, 17 distal) obtained from unrelated subjects (21 female, 18 male) during screening colonoscopy were grown under low (1.7 mM) calcium conditions, and a matching set of the same lines were grown under high (5.0 mM) calcium conditions. After 72 hours, total RNA was extracted from organoid lines, libraries were prepared with poly-A selection, and bulk RNA-sequencing was used to measure gene expression. DNA was extracted from the same lines and genotyped using the OncoArray 500K beadchip. Missing genotypes were imputed to the TOPMed reference panel using the TOPMed Imputation Server. Polygenic risk scores for colorectal cancer were calculated for each organoid line using Plink and 135 of 141 genome-wide significant SNPs associated with colorectal cancer in recent genome-wide association studies. Relative expression between low and high calcium conditions was tested using generalized linear models as implemented in DESeq2. Functional enrichment of gene ontologies among differentially expressed genes was performed with ToppFun. The effects of donor sex, biopsy subsite, and organoid polygenic risk score on response to calcium were tested using single and multiple linear regression as implemented in the base R lm function. Results: Exposure to calcium induced differential expression of 767 genes with a false discovery rate of 5%. Differentially expressed genes were enriched for biological processes including extracellular structure organization (Bonferroni p=1.6E-3) and positive regulation of apoptosis (Bonferroni p=2.2E-2). Although stratified analyses demonstrated more differentially expressed genes in 3D colon organoids derived from male subjects, no statistically significant effect of sex, subsite, or polygenic risk score on response to calcium was observed. Conclusion: We tested the effect of calcium exposure on gene expression in human-derived normal colon 3D organoids and found that calcium modulated processes related to differentiation and apoptosis, but the response to calcium was not affected by sex, anatomic subsite, or polygenic risk score. Citation Format: Christopher H. Dampier, Matthew Devall, Jennifer Bryant, Stephen D. Eaton, Jeroen R. Huyghe, Andre E. Kim, Virginia Diez-Obrero, Franzel J. van Duijnhoven, Duncan C. Thomas, Steven M. Powell, W James Gauderman, Ulrike Peters, Graham Casey. Calcium modulates apoptosis and differentiation in human normal colon 3D organoids irrespective of donor sex and biopsy site [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2578.

A head-to-head comparison of ribodepletion and polyA selection approaches for Caenorhabditis elegans low input RNA-sequencing libraries.

A recent and powerful technique is to obtain transcriptomes from rare cell populations, such as single neurons in Caenorhabditis elegans, by enriching dissociated cells using fluorescent sorting. However, these cell samples often have low yields of RNA that present challenges in library preparation. This can lead to PCR duplicates, noisy gene expression for lowly expressed genes, and other issues that limit endpoint analysis. Furthermore, some common resources, such as sequence-specific kits for removing ribosomal RNA, are not optimized for nonmammalian samples. To advance library construction for such challenging samples, we compared two approaches for building RNAseq libraries from less than 10 nanograms of C. elegans RNA: SMARTSeq V4 (Takara), a widely used kit for selecting poly-adenylated transcripts; and SoLo Ovation (Tecan Genomics), a newly developed ribodepletion-based approach. For ribodepletion, we used a custom kit of 200 probes designed to match C. elegans rRNA gene sequences. We found that SoLo Ovation, in combination with our custom C. elegans probe set for rRNA depletion, detects an expanded set of noncoding RNAs, shows reduced noise in lowly expressed genes, and more accurately counts expression of long genes. The approach described here should be broadly useful for similar efforts to analyze transcriptomics when RNA is limiting.

Comparative evaluation for the globin gene depletion methods for mRNA sequencing using the whole blood-derived total RNAs

BackgroundThere are challenges in generating mRNA-Seq data from whole-blood derived RNA as globin gene and rRNA are frequent contaminants. Given the abundance of erythrocytes in whole blood, globin genes comprise some 80% or more of the total RNA. Therefore, depletion of globin gene RNA and rRNA are critical steps required to have adequate coverage of reads mapping to the reference transcripts and thus reduce the total cost of sequencing. In this study, we directly compared the performance of probe hybridization (GLOBINClear Kit and Globin-Zero Gold rRNA Removal Kit) and RNAse-H enzymatic depletion (NEBNext® Globin & rRNA Depletion Kit and Ribo-Zero Plus rRNA Depletion Kit) methods from 1 μg of whole blood-derived RNA on mRNA-Seq profiling. All RNA samples were treated with DNaseI for additional cleanup before the depletion step and were processed for poly-A selection for library generation.ResultsProbe hybridization revealed a better overall performance than the RNAse-H enzymatic depletion method, detecting a higher number of genes and transcripts without 3′ region bias. After depletion, samples treated with probe hybridization showed globin genes at 0.5% (±0.6%) of the total mapped reads; the RNAse-H enzymatic depletion had 3.2% (±3.8%). Probe hybridization showed more junction reads and transcripts compared with RNAse-H enzymatic depletion and also had a higher correlation (R > 0.9) than RNAse-H enzymatic depletion (R > 0.85).ConclusionIn this study, our results showed that 1 μg of high-quality RNA from whole blood could be routinely used for transcriptional profiling analysis studies with globin gene and rRNA depletion pre-processing. We also demonstrated that the probe hybridization depletion method is better suited to mRNA sequencing analysis with minimal effect on RNA quality during depletion procedures.

RNAseq Analysis of Rodent Spaceflight Experiments Is Confounded by Sample Collection Techniques.

SummaryTo understand the physiological changes that occur in response to spaceflight, mice are transported to the International Space Station (ISS) and housed for variable periods of time before euthanasia on-orbit or return to Earth. Sample collection under such difficult conditions introduces confounding factors that need to be identified and addressed. We found large changes in the transcriptome of mouse tissues dissected and preserved on-orbit compared with tissues from mice euthanized on-orbit, preserved, and dissected after return to Earth. Changes due to preservation method eclipsed those between flight and ground samples, making it difficult to identify spaceflight-specific changes. Follow-on experiments to interrogate the roles of euthanasia methods, tissue and carcass preservation protocols, and library preparation methods suggested that differences due to preservation protocols are exacerbated when coupled with polyA selection. This has important implications for the interpretation of existing datasets and the design of future experiments.

RNAseq Analysis of Rodent Spaceflight Experiments is Confounded by Sample Collection Techniques

To understand the physiological changes that occur in response to spaceflight, mice are transported to the International Space Station (ISS) and housed for variable periods of time before euthanasia on-orbit or return to Earth. Sample collection under such difficult conditions introduces confounding factors that need to be identified and addressed. We found large changes in the transcriptome of mouse tissues dissected and preserved on-orbit compared to tissues from mice euthanized on-orbit, preserved, and dissected after return to Earth. Changes due to preservation method eclipsed those between flight and ground samples making it difficult to identify spaceflight-specific changes. Follow-on experiments to interrogate the roles of euthanasia methods, tissue and carcass preservation protocols, and library preparation methods suggested that differences due to preservation protocols are exacerbated when coupled with polyA selection. This has important implications for the interpretation of existing datasets and the design of future experiments.

LITHOPHONE: Improving lncRNA Methylation Site Prediction Using an Ensemble Predictor.

N6-methyladenosine (m6A) is one of the most widely studied epigenetic modifications, which plays an important role in many biological processes, such as splicing, RNA localization, and degradation. Studies have shown that m6A on lncRNA has important functions, including regulating the expression and functions of lncRNA, regulating the synthesis of pre-mRNA, promoting the proliferation of cancer cells, and affecting cell differentiation and many others. Although a number of methods have been proposed to predict m6A RNA methylation sites, most of these methods aimed at general m6A sites prediction without noticing the uniqueness of the lncRNA methylation prediction problem. Since many lncRNAs do not have a polyA tail and cannot be captured in the polyA selection step of the most widely adopted RNA-seq library preparation protocol, lncRNA methylation sites cannot be effectively captured and are thus likely to be significantly underrepresented in existing experimental data affecting the accuracy of existing predictors. In this paper, we propose a new computational framework, LITHOPHONE, which stands for long noncoding RNA methylation sites prediction from sequence characteristics and genomic information with an ensemble predictor. We show that the methylation sites of lncRNA and mRNA have different patterns exhibited in the extracted features and should be differently handled when making predictions. Due to the used experiment protocols, the number of known lncRNA m6A sites is limited, and insufficient to train a reliable predictor; thus, the performance can be improved by combining both lncRNA and mRNA data using an ensemble predictor. We show that the newly developed LITHOPHONE approach achieved a reasonably good performance when tested on independent datasets (AUC: 0.966 and 0.835 under full transcript and mature mRNA modes, respectively), marking a substantial improvement compared with existing methods. Additionally, LITHOPHONE was applied to scan the entire human lncRNAome for all possible lncRNA m6A sites, and the results are freely accessible at: http://180.208.58.19/lith/.

Identification of Circular RNAs using RNA Sequencing.

Circular RNAs (circRNAs) are a class of non-coding RNAs involved in functions including micro-RNA (miRNA) regulation, mediation of protein-protein interactions, and regulation of parental gene transcription. In classical next generation RNA sequencing (RNA-seq), circRNAs are typically overlooked as a result of poly-A selection during construction of mRNA libraries, or are found at very low abundance, and are therefore difficult to isolate and detect. Here, a circRNA library construction protocol was optimized by comparing library preparation kits, pre-treatment options and various total RNA input amounts. Two commercially available whole transcriptome library preparation kits, with and without RNase R pre-treatment, and using variable amounts of total RNA input (1 to 4 µg), were tested. Lastly, multiple tissue types; including liver, lung, lymph node, and pancreas; as well as multiple brain regions; including the cerebellum, inferior parietal lobe, middle temporal gyrus, occipital cortex, and superior frontal gyrus; were compared to evaluate circRNA abundance across tissue types. Analysis of the generated RNA-seq data using six different circRNA detection tools (find_circ, CIRI, Mapsplice, KNIFE, DCC, and CIRCexplorer) revealed that a stranded total RNA library preparation kit with RNase R pre-treatment and 4 µg RNA input is the optimal method for identifying the highest relative number of circRNAs. Consistent with previous findings, the highest enrichment of circRNAs was observed in brain tissues compared to other tissue types.

SCAF4 and SCAF8, mRNA Anti-Terminator Proteins

SummaryAccurate regulation of mRNA termination is required for correct gene expression. Here, we describe a role for SCAF4 and SCAF8 as anti-terminators, suppressing the use of early, alternative polyadenylation (polyA) sites. The SCAF4/8 proteins bind the hyper-phosphorylated RNAPII C-terminal repeat domain (CTD) phosphorylated on both Ser2 and Ser5 and are detected at early, alternative polyA sites. Concomitant knockout of human SCAF4 and SCAF8 results in altered polyA selection and subsequent early termination, leading to expression of truncated mRNAs and proteins lacking functional domains and is cell lethal. While SCAF4 and SCAF8 work redundantly to suppress early mRNA termination, they also have independent, non-essential functions. SCAF8 is an RNAPII elongation factor, whereas SCAF4 is required for correct termination at canonical, distal transcription termination sites in the presence of SCAF8. Together, SCAF4 and SCAF8 coordinate the transition between elongation and termination, ensuring correct polyA site selection and RNAPII transcriptional termination in human cells.

Widespread and Functional RNA Circularization in Localized Prostate Cancer

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree ofcircRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Statistical inference of the rate of RNA polymerase II elongation by total RNA sequencing.

MotivationSequencing total RNA without poly-A selection enables us to obtain a transcriptomic profile of nascent RNAs undergoing transcription with co-transcriptional splicing. In general, the RNA-seq reads exhibit a sawtooth pattern in a gene, which is characterized by a monotonically decreasing gradient across introns in the 5’–3’ direction, and by substantially higher levels of RNA-seq reads present in exonic regions. Such patterns result from the process of underlying transcription elongation by RNA polymerase II, which traverses the DNA strand in a 5’–3’ direction as it performs a complex series of mRNA synthesis and processing. Therefore, data of sequenced total RNAs could be utilized to infer the rate of transcription elongation by solving the inverse problem.ResultsThough solving the inverse problem in total RNA-seq has the great potential, statistical methods have not yet been fully developed. We demonstrate what extent the newly developed method can be useful. The objective is to reconstruct the spatial distribution of transcription elongation rates in a gene from a given noisy, sawtooth-like profile. It is necessary to recover the signal source of the elongation rates separately from several types of nuisance factors, such as unobserved modes of co-transcriptionally occurring mRNA splicing, which exert significant influences on the sawtooth shape. The present method was tested using published total RNA-seq data derived from mouse embryonic stem cells. We investigated the spatial characteristics of the estimated elongation rates, focusing especially on the relation to promoter-proximal pausing of RNA polymerase II, nucleosome occupancy and histone modification patterns.Availability and implementationA C implementation of PolSter and sample data are available at https://github.com/yoshida-lab/PolSter.Supplementary information Supplementary data are available at Bioinformatics online.

Nutrient depletion and TOR inhibition induce 18S and 25S ribosomal RNAs resistant to a 5\u2032-phosphate-dependent exonuclease in Candida albicans and other yeasts

BackgroundMessenger RNA (mRNA) represents a small percentage of RNAs in a cell, with ribosomal RNA (rRNA) making up the bulk of it. To isolate mRNA from eukaryotes, typically poly-A selection is carried out. Recently, a 5´-phosphate-dependent, 5´→3´ processive exonuclease called Terminator has become available. It will digest only RNA that has a 5´-monophosphate end and therefore it is very useful to eliminate most of rRNAs in cell.ResultsWe have found that in the pathogenic yeast Candida albicans, while 18S and 25S components isolated from yeast in robust growth phase are easily eliminated by Terminator, those isolated from cells in the nutritionally diminished stationary phase, become resistant to digestion by this enzyme. Additional digestions with alkaline phosphatase, tobacco pyrophosphatase combined with Terminator point toward the 5′-prime end of 18S and 25S as the source of this resistance. Inhibition of TOR by rapamycin also induces resistance by these molecules. We also find that these molecules are incorporated into the ribosome and are not just produced incidentally. Finally, we show that three other yeasts show the same behavior.ConclusionsDigestion of RNA by Terminator has revealed 18S and 25S rRNA molecules different from the accepted processed ones seen in ribosome generation. The reason for these molecules and the underlying mechanism for their formation is unknown. The preservation of this behavior across these yeasts suggests a useful biological role for it, worthy of further inquiry.

Don't just dump your data and run: Authors should submit as much experimental information as possible when uploading sequence data.

If you have, in any way, been involved with genetic research over the past 10 years, then you have likely heard of the Sequence Read Archive (SRA), which is jointly housed at the National Center for Biotechnology (NCBI), the DNA Data Bank of Japan (DDBJ), and the European Bioinformatics Institute (EBI). And if you regularly work with genome or transcriptome sequence information, then you have probably extracted data from and/or deposited data into the SRA. For those who are unfamiliar with it, the SRA is an international public online archive for next‐generation sequencing (NGS) data, which was established about a decade ago under the guidance of the International Nucleotide Sequence Database Collaboration (INSDC) [1], [2]. Despite nearly being shut down in 2011 [3], it has grown at a staggering rate over the past 10 years. As of September 1, 2017, it housed over five quadrillion (1015) open‐access bases of NGS data, coming from thousands of different species and spanning the entire gamut of cellular and viral life. It contains DNA‐ and RNA‐sequencing (DNA‐seq and RNA‐seq) reads of every kind, from bisulfite‐seq to strand‐specific RNA‐seq to single‐cell DNA‐seq, and it accepts reads from every type of NGS platform, be it Illumina, Ion Torrent, or PacBio sequencing. In other words, the SRA is a crucial and central resource in the fast‐paced and increasingly important domain of contemporary genetic research. > … the SRA is a crucial and central resource in the fast‐paced and increasingly important domain of contemporary genetic research. ### The Sequence Read Archive The SRA can be easily accessed and searched via the NCBI (https://www.ncbi.nlm.nih.gov/sra), DDBJ (http://trace.ddbj.nig.ac.jp/dra/index_e.html), and EBI (http://www.ebi.ac.uk/ena/submit/read-submission) websites. Once there, you will find yourself at a sequencing‐read superstore. With a decent Wi‐Fi connection, a couple of keyword searches, and a few clicks of the trackpad, you can quickly …

Abstract POSTER-TECH-1109: Robust gene expression and mutation analyses from RNA-sequencing of formalin-fixed diagnostic tumor samples

Abstract Formalin-Fixed Paraffin Embedding (FFPE) has been the standard sample preparation for pathologists for decades. This created a mass of disease and normal tissue biospecimens with rich clinical annotation and patient follow-up data. If these FFPE archives proved useful for the next-generation sequence analyses, FFPE samples may be sequenced to investigate the complex genetic changes underlying tumor progression, therapy resistance and variability in disease outcome. Previous studies have successfully used FFPE-DNA for copy number analysis and mutation detection even though it is technically difficult to perform sequencing analyses on DNA isolated from FFPE samples. It is even more challenging to perform RNA sequence analysis because routine FFPE processing does not preserve high quality RNA. Aim of this study is to determine the extent to which FFPE RNA expression is correlated with matched frozen tumor samples. Six pairs of matched fresh frozen (FF) and FFPE ovarian tumor samples were processed for RNA sequence analysis. Multiple quality checks were performed before and after aligning the reads to the reference genome (hg19/GRCh37). Differential gene expression analysis was performed utilizing DESeq on aligned RNA-sequencing data, as well as data obtained by NanoString technology (250 selected genes). Both results were compared to each other. A pipeline called SNPiR was applied to detect mutation from the RNA sequencing data, and the mutation calls between the FF and FFPE samples were compared at the RNA and DNA level. For each analysis FF samples from the same cancer patients were used to verify the results from FFPE sample. The quality checks indicate that FFPE samples show uniform coverage of 5’ and 3’ of the transcripts whereas FF samples show 3’ bias in transcript coverage. This difference is likely attributable to differences in library preparation steps for FF and FFPE samples. The sequencing libraries for FF samples were done by poly-A selection while the libraries for FFPE samples were prepared using Ribo- Minus RNA sequencing kit. The differential gene expression analysis showed a good correlation between each matching FF and FFPE samples with one exception. Analyzing the NanoString gene counts show not only a remarkable correlation between the matching FF and FFPE sample but also a correlation between FF and FFPE samples between the two techniques (RNA-sequencing and NanoString). Unfortunately the mutation calling show variable results with 3 paired samples showing high rates of concordant mutation calls between FF and FFPE samples while the remaining 3 paired samples showing low rates of concordant mutation calls. Ongoing studies are focusing on biological and technical variability that may account for the observed differences in concordant calls. In summary, the results indicate that FFPE samples are a suitable replacement for FF samples in differential gene expression analysis. However, mutation analysis from FFPE samples is challenging and results are variable. Citation Format: Stefan Graw, Richard Meier, Kay Minn, Andrew K Godwin, Peter Beyerlein, Jeremy Chien. Robust gene expression and mutation analyses from RNA-sequencing of formalin-fixed diagnostic tumor samples [abstract]. In: Proceedings of the 10th Biennial Ovarian Cancer Research Symposium; Sep 8-9, 2014; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(16 Suppl):Abstract nr POSTER-TECH-1109.