Metabolomics,Unknown,Transcriptomics,Genomics,Proteomics

Dataset Information

High-resolution transcriptome analysis with long-read RNA sequencing

ABSTRACT: Ongoing improvements to next generation sequencing technologies are leading to longer sequencing read lengths, but a thorough understanding of the impact of longer reads on RNA sequencing analyses is lacking. To address this issue, we generated and compared two RNA sequencing datasets of differing read lengths -- 2x75 bp (L75) and 2x262 bp (L262) -- and investigated the impact of read length on various aspects of analysis, including the performance of currently available read-mapping tools, gene and transcript quantification, and detection of allele-specific expression patterns. Our results indicate that, while the scalability of read-mapping tools and the cost-effectiveness of long read protocol is an issue that requires further attention, longer reads enable more accurate quantification of diverse aspects of gene expression, including individual-specific patterns of allele-specific expression and alternative splicing. Two RNA-Seq datasets of differing read lengths (2x262 bp and 2x75 bp)

ORGANISM(S): Homo sapiens

SUBMITTER: Hyunghoon Cho

PROVIDER: E-GEOD-57862 | biostudies-arrayexpress |

REPOSITORIES: biostudies-arrayexpress

ACCESS DATA

Publications

High-resolution transcriptome analysis with long-read RNA sequencing.

Cho Hyunghoon H Davis Joe J Li Xin X Smith Kevin S KS Battle Alexis A Montgomery Stephen B SB

PloS one 20140924 9

RNA sequencing (RNA-seq) enables characterization and quantification of individual transcriptomes as well as detection of patterns of allelic expression and alternative splicing. Current RNA-seq protocols depend on high-throughput short-read sequencing of cDNA. However, as ongoing advances are rapidly yielding increasing read lengths, a technical hurdle remains in identifying the degree to which differences in read length influence various transcriptome analyses. In this study, we generated two ...[more]

PMID: 25251678

Similar Datasets

Project description:During interphase, the inactive X chromosome (Xi) adopts an unusual 3D configuration known as the Barr body and is largely transcriptionally silent. Despite the importance of X inactivation, little is known about the 3D configuration of Xi and its relationship to gene silencing. We recently showed that in humans, Xi exhibits two distinctive structural features. First, Xi is partitioned into two huge intervals, called superdomains, such that pairs of loci in each superdomain show an enhanced contact frequency with one another. The boundary between the two superdomains lies near DXZ4, a macrosatellite repeat spanning ~300kb, whose Xi allele extensively binds the protein CTCF. Second, Xi exhibits extremely large loops, up to 77Mb long, called superloops. DXZ4 lies at the anchor of several superloops. Here, we use 3D mapping to study the structure of Xi, focusing on the role of DXZ4. We show that superloops and superdomains are conserved across mammals. We develop a novel variant of our in situ Hi-C protocol, dubbed COLA (COncatemer Ligation Assay) to probe the higher order structures formed by the superloops. In COLA, in situ proximity ligation of multiple extremely short fragments produced by the enzyme CviJI is used to efficiently map simultaneous proximity among three or more loci. Using data from Hi-C and COLA, we demonstrate that DXZ4 and other superloop anchors tend to co-locate simultaneously within the same cells, a result that is confirmed by 3D-FISH. Finally, we examine the effects of deleting DXZ4 from Xi in human cells. Using in situ Hi-C, microscopy, and RNA-FISH, we show that superdomains and superloops disappear; that Xi frequently dissociates into multiple separate structures; and that transcriptional silencing on Xi is compromised. Deletion of DXZ4 from the active X chromosome (Xa) has no such effect. Thus, DXZ4 is essential for proper folding and silencing of Xi. Hi-C protocol was used on wildtype and DXZ4-deleted cells to examine the structure of Xi. A novel variant of our in situ Hi-C protocol, dubbed COLA (COncatemer Ligation Assay), was developed to probe the higher order structures formed by the superloops. This series also includes RNA-seq data on Retinal Pigmented Epithelial Cells (hTERT-RPE1). At the time of submission, processed data were available only for the RNA-seq samples. Submitter states that processed data files for HiC samples will be added to this series in the future.

Project description:In platelets, splicing and translation occur in the absence of a nucleus. However, the integrity and stability of mRNAs derived from megakaryocyte progenitor cells remain poorly quantified on a transcriptome-wide level. As circular RNAs (circRNAs) are resistant to degradation by exonucleases, their abundance relative to linear RNAs can be used as a surrogate marker for mRNA stability in the absence of transcription. Here we show that circRNAs are enriched in human platelets 17-188 fold relative to nucleated tissues, and 14-26 fold relative to samples digested with RNAseR to selectively remove linear RNA. We compare RNAseq read depths inside and outside circRNAs to provide in silico evidence of transcript circularity, show that exons within circRNAs are enriched ~13X in platelets relative to nucleated tissues, and identify 3162 genes significantly enriched for circRNAs including some where all RNAs appear to be derived from circular molecules. We also confirm that this is a feature of other anucleate cells through transcriptome sequencing of mature erythrocytes, demonstrate that circRNAs are not enriched in megakaryocytes, and that linear RNAs decay more rapidly than circRNAs in platelet preparations. Collectively, these results suggest that circulating platelets have lost on aveage over 90% of their progenitor mRNAs, and that translation in platelets occurrs against the backdrop of a highly degraded transcriptome. Finally, we find that transcripts classified as products of reverse transcriptase template switching are both enriched in platelets and resistant to decay, countering the recent suggestion that up to 50% of rearranged RNAs are artefacts. A single rRNA depleted total RNA sample was sequenced. This together with 25 publicly available rRNA depleted total RNA samples (including 3 from platelets) were analysed using PTESFinder v 1 (http://sourceforge.net/projects/ptesfinder-v1/) to identify back-splice junctions, characteristic of circRNA transcripts. The contribution of circRNA producing exons was analysed on a gene by gene basis as follows: All circRNA transcripts identified in any sample were first pooled to define exons which can contribute to circRNA generation using custom scripts (available on request). For each sample, expression estimates (RPKMI) across all circRNA producing exons were computed for each locus using the total size of exons (in bp) and the read counts mapping to them. Similarly, total size and exonic read counts for exons for which no circRNA were detected in any sample were used to compute expression estimates (RPKME) for non-circRNA producing exons for each locus. Abundance ratios (RPKMI/RPKME and RPKMI/RPKMI+RPKME) were calculated and compared between Platelets and human tissues using Wilcoxon signed-rank test. Please note that the '25sample_info_accn_no.txt' contains the accession numbers and tissue/cell type information for 25 samples analyzed together.

Dataset Information

High-resolution transcriptome analysis with long-read RNA sequencing

Publications

High-resolution transcriptome analysis with long-read RNA sequencing.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure