Dataset Information

Gene expression profiling of non-polyadenylated RNA-seq across species.

ABSTRACT: Transcriptomes are dynamic and unique, with each cell type/tissue, developmental stage and species expressing a different repertoire of RNA transcripts. Most mRNAs and well-characterized long noncoding RNAs are shaped with a 5' cap and 3' poly(A) tail, thus conventional transcriptome analyses typically start with the enrichment of poly(A)+ RNAs by oligo(dT) selection, followed by deep sequencing approaches. However, accumulated lines of evidence suggest that many RNA transcripts are processed by alternative mechanisms without 3' poly(A) tails and, therefore, fail to be enriched by oligo(dT) purification and are absent following deep sequencing analyses. We have described an enrichment strategy to purify non-polyadenylated (poly(A)-/ribo-) RNAs from human total RNAs by removal of both poly(A)+ RNA transcripts and ribosomal RNAs, which led to the identification of many novel RNA transcripts with non-canonical 3' ends in human. Here, we describe the application of non-polyadenylated RNA-sequencing in rhesus monkey and mouse cell lines/tissue, and further profile the transcription of non-polyadenylated RNAs across species, providing new resources for non-polyadenylated RNA identification and comparison across species.

SUBMITTER: Zhang XO

PROVIDER: S-EPMC4535946 | biostudies-literature | 2014 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Gene expression profiling of non-polyadenylated RNA-seq across species.

Zhang Xiao-Ou XO Yin Qing-Fei QF Chen Ling-Ling LL Yang Li L

Genomics data 20140803

Transcriptomes are dynamic and unique, with each cell type/tissue, developmental stage and species expressing a different repertoire of RNA transcripts. Most mRNAs and well-characterized long noncoding RNAs are shaped with a 5' cap and 3' poly(A) tail, thus conventional transcriptome analyses typically start with the enrichment of poly(A)+ RNAs by oligo(dT) selection, followed by deep sequencing approaches. However, accumulated lines of evidence suggest that many RNA transcripts are processed by ...[more]

PMID: 26484100

Similar Datasets

Project description:BackgroundThe advancement of the next-generation sequencing technology enables mapping gene expression at the single-cell level, capable of tracking cell heterogeneity and determination of cell subpopulations using single-cell RNA sequencing (scRNA-seq). Unlike the objectives of conventional RNA-seq where differential expression analysis is the integral component, the most important goal of scRNA-seq is to identify highly variable genes across a population of cells, to account for the discrete nature of single-cell gene expression and uniqueness of sequencing library preparation protocol for single-cell sequencing. However, there is lack of generic expression variation model for different scRNA-seq data sets. Hence, the objective of this study is to develop a gene expression variation model (GEVM), utilizing the relationship between coefficient of variation (CV) and average expression level to address the over-dispersion of single-cell data, and its corresponding statistical significance to quantify the variably expressed genes (VEGs).ResultsWe have built a simulation framework that generated scRNA-seq data with different number of cells, model parameters, and variation levels. We implemented our GEVM and demonstrated the robustness by using a set of simulated scRNA-seq data under different conditions. We evaluated the regression robustness using root-mean-square error (RMSE) and assessed the parameter estimation process by varying initial model parameters that deviated from homogeneous cell population. We also applied the GEVM on real scRNA-seq data to test the performance under distinct cases.ConclusionsIn this paper, we proposed a gene expression variation model that can be used to determine significant variably expressed genes. Applying the model to the simulated single-cell data, we observed robust parameter estimation under different conditions with minimal root mean square errors. We also examined the model on two distinct scRNA-seq data sets using different single-cell protocols and determined the VEGs. Obtaining VEGs allowed us to observe possible subpopulations, providing further evidences of cell heterogeneity. With the GEVM, we can easily find out significant variably expressed genes in different scRNA-seq data sets.

Project description:Leishmania is a protozoan parasite that alternates its life cycle between the sand-fly vector and the mammalian host. This alternation involves environmental changes and leads the parasite to dynamic modifications in morphology, metabolism, cellular signaling and regulation of gene expression to allow for a rapid adaptation to new conditions. The L-arginine pathway in L. amazonensis is important during the parasite life cycle and interferes in the establishment and maintenance of the infection in mammalian macrophages. Host arginase is an immune-regulatory enzyme that can reduce the production of nitric oxide by activated macrophages, directing the availability of L-arginine to the polyamine pathway, resulting in parasite replication. In this work, we performed transcriptional profiling to identify differentially expressed genes in L. amazonensis wild-type (La-WT) versus L. amazonensis arginase knockout (La-arg-) promastigotes and axenic amastigotes.A total of 8253 transcripts were identified in La-WT and La-arg- promastigotes and axenic amastigotes, about 60% of them codifying hypothetical proteins and 443 novel transcripts, which did not match any previously annotated genes. Our RNA-seq data revealed that 85% of genes were constitutively expressed. The comparison of transcriptome and metabolome data showed lower levels of arginase and higher levels of glutamate-5-kinase in La-WT axenic amastigotes compared to promastigotes. The absence of arginase activity in promastigotes increased the levels of pyrroline 5-carboxylate reductase, but decreased the levels of arginosuccinate synthase, pyrroline 5-carboxylate dehydrogenase, acetylornithine deacetylase and spermidine synthase transcripts levels. These observations can explain previous metabolomic data pointing to the increase of L-arginine, citrulline and L-glutamate and reduction of aspartate, proline, ornithine and putrescine. Altogether, these results indicate that arginase activity is important in Leishmania gene expression modulation during differentiation and adaptation to environmental changes. Here, we confirmed this hypothesis with the identification of differential gene expression of the enzymes involved in biosynthesis of amino acids, arginine and proline metabolism and arginine biosynthesis.All data provided information about the transcriptomic profiling and the expression levels of La-WT and La-arg- promastigotes and axenic amastigotes. These findings revealed the importance of arginase in parasite survival and differentiation, and indicated the existence of a coordinated response in the absence of arginase activity related to arginine and polyamine pathways.

Project description:Deep sequencing was used to investigate the subchronic effects of 1 ppm aflatoxin B1 (AFB1), a potent hepatocarcinogen, on the male rat liver transcriptome prior to onset of histopathological lesions or tumors. We hypothesized RNA-Seq would reveal more differentially expressed genes (DEG) than microarray analysis, including low copy and novel transcripts related to AFB1's carcinogenic activity compared to feed controls (CTRL). Paired-end reads were mapped to the rat genome (Rn4) with TopHat and further analyzed by DESeq and Cufflinks-Cuffdiff pipelines to identify differentially expressed transcripts, new exons and unannotated transcripts. PCA and cluster analysis of DEGs showed clear separation between AFB1 and CTRL treatments and concordance among group replicates. qPCR of eight high and medium DEGs and three low DEGs showed good comparability among RNA-Seq and microarray transcripts. DESeq analysis identified 1,026 differentially expressed transcripts at greater than two-fold change (p<0.005) compared to 626 transcripts by microarray due to base pair resolution of transcripts by RNA-Seq, probe placement within transcripts or an absence of probes to detect novel transcripts, splice variants and exons. Pathway analysis among DEGs revealed signaling of Ahr, Nrf2, GSH, xenobiotic, cell cycle, extracellular matrix, and cell differentiation networks consistent with pathways leading to AFB1 carcinogenesis, including almost 200 upregulated transcripts controlled by E2f1-related pathways related to kinetochore structure, mitotic spindle assembly and tissue remodeling. We report 49 novel, differentially-expressed transcripts including confirmation by PCR-cloning of two unique, unannotated, hepatic AFB1-responsive transcripts (HAfT's) on chromosomes 1.q55 and 15.q11, overexpressed by 10 to 25-fold. Several potentially novel exons were found and exon refinements were made including AFB1 exon-specific induction of homologous family members, Ugt1a6 and Ugt1a7c. We find the rat transcriptome contains many previously unidentified, AFB1-responsive exons and transcripts supporting RNA-Seq's capabilities to provide new insights into AFB1-mediated gene expression leading to hepatocellular carcinoma.

Project description:BackgroundWhile RNA-sequencing (RNA-seq) is becoming a powerful technology in transcriptome profiling, one significant shortcoming of the first-generation RNA-seq protocol is that it does not retain the strand specificity of origin for each transcript. Without strand information it is difficult and sometimes impossible to accurately quantify gene expression levels for genes with overlapping genomic loci that are transcribed from opposite strands. It has recently become possible to retain the strand information by modifying the RNA-seq protocol, known as strand-specific or stranded RNA-seq. Here, we evaluated the advantages of stranded RNA-seq in transcriptome profiling of whole blood RNA samples compared with non-stranded RNA-seq, and investigated the influence of gene overlaps on gene expression profiling results based on practical RNA-seq datasets and also from a theoretical perspective.ResultsOur results demonstrated a substantial impact of stranded RNA-seq on transcriptome profiling and gene expression measurements. As many as 1751 genes in Gencode Release 19 were identified to be differentially expressed when comparing stranded and non-stranded RNA-seq whole blood samples. Antisense and pseudogenes were significantly enriched in differential expression analyses. Because stranded RNA-seq retains strand information of a read, we can resolve read ambiguity in overlapping genes transcribed from opposite strands, which provides a more accurate quantification of gene expression levels compared with traditional non-stranded RNA-seq. In the human genome, it is not uncommon to find genomic loci where both strands encode distinct genes. Among the over 57,800 annotated genes in Gencode release 19, there are an estimated 19 % (about 11,000) of overlapping genes transcribed from the opposite strands. Based on our whole blood mRNA-seq datasets, the fraction of overlapping nucleotide bases on the same and opposite strands were estimated at 2.94 % and 3.1 %, respectively. The corresponding theoretical estimations are 3 % and 3.6 %, well in agreement with our own findings.ConclusionsStranded RNA-seq provides a more accurate estimate of transcript expression compared with non-stranded RNA-seq, and is therefore the recommended RNA-seq approach for future mRNA-seq studies.

Project description:Graphical abstract The cGAS-STING pathway, orchestrating complicated transcriptome-wide immune responses, is essential for host antiviral defense but can also drive immunopathology in severe COVID-19. Here, we performed time-course RNA-Seq experiments to dissect the transcriptome expression dynamics at the gene-isoform level after cGAS-STING pathway activation. The in-depth time-course transcriptome after cGAS-STING pathway activation within 12 h enabled quantification of 48,685 gene isoforms. By employing regression models, we obtained 13,232 gene isoforms with expression patterns significantly associated with the process of cGAS-STING pathway activation, which were named activation-associated isoforms. The combination of hierarchical and k-means clustering algorithms revealed four major expression patterns of activation-associated isoforms, including two clusters with increased expression patterns enriched in cell cycle, autophagy, antiviral innate-immune functions, and COVID-19 coronavirus disease pathway, and two clusters showing decreased expression pattern that mainly involved in ncRNA metabolism, translation process, and mRNA processing. Importantly, by merging four clusters of activation-associated isoforms, we identified three types of genes that underwent isoform usage alteration during the cGAS-STING pathway activation. We further found that genes exhibiting protein-coding and non-protein-coding gene isoform usage alteration were strongly enriched for the factors involved in innate immunity and RNA splicing. Notably, overexpression of an enriched splicing factor, EFTUD2, shifted transcriptome towards the cGAS-STING pathway activated status and promoted protein-coding isoform abundance of several key regulators of the cGAS-STING pathway. Taken together, our results revealed the isoform-level gene expression dynamics of the cGAS-STING pathway and uncovered novel roles of splicing factors in regulating cGAS-STING pathway mediated immune responses.

Dataset Information

Gene expression profiling of non-polyadenylated RNA-seq across species.

Publications

Gene expression profiling of non-polyadenylated RNA-seq across species.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets