Digital gene expression and global mapping of polyadenylation sites with PolyA-Seq
Ontology highlight
ABSTRACT:
ORGANISM(S): Homo sapiens
SUBMITTER: Adnan Derti
PROVIDER: E-GEOD-30198 | biostudies-arrayexpress |
REPOSITORIES: biostudies-arrayexpress
Publications
A quantitative atlas of polyadenylation in five mammals.
Derti Adnan A Garrett-Engele Philip P Macisaac Kenzie D KD Stevens Richard C RC Sriram Shreedharan S Chen Ronghua R Rohl Carol A CA Johnson Jason M JM Babak Tomas T
Genome research 20120327 6
We developed PolyA-seq, a strand-specific and quantitative method for high-throughput sequencing of 3' ends of polyadenylated transcripts, and used it to globally map polyadenylation (polyA) sites in 24 matched tissues in human, rhesus, dog, mouse, and rat. We show that PolyA-seq is as accurate as existing RNA sequencing (RNA-seq) approaches for digital gene expression (DGE), enabling simultaneous mapping of polyA sites and quantitative measurement of their usage. In human, we confirmed 158,533 ...[more]
Similar Datasets
Project description:For some neurological disorders, disease is primarily RNA-mediated due to expression of non-coding microsatellite expansion RNAs (RNAexp). Toxicity is thought to result from enhanced binding of proteins to these expansions and depletion from their normal cellular targets. However, experimental evidence for this sequestration model is lacking. Here, we use HITS-CLIP and pre-mRNA processing analysis of human control versus myotonic dystrophy (DM) brains to provide compelling evidence for this RNA toxicity model. MBNL2 binds directly to DM repeat expansions in the brain resulting in depletion from its normal RNA targets with downstream effects on alternative splicing and polyadenylation. Similar RNA processing defects were detected in Mbnl compound knockout mice, highlighted by dysregulation of Mapt splicing and fetal tau isoform expression in adults. These results demonstrate that MBNL proteins are directly sequestered by RNAexp in the DM brain and introduce a powerful experimental tool to evaluate RNA-mediated toxicity in other expansion diseases. HITS-CLIP analysis was performed to identify RNA binding sites of MBNL2 in control, DM type 1 (DM1), and DM type 2 (DM2) autopsy-derived brain (n=3). Two regions of the brain selected for the study included the frontal cortex and hippocampus. Libraries were sequenced and wiggle files were generated for each biological replicate as well as three pooled biological replicates (3BRs) per group (control, DM1, and DM2). In addition, differential CLIP analysis (dCLIP) was performed to normalize binding data between groups and identify statistically significant changes in binding. The dCLIP analysis generated bedgraph files representing normalized binding profiles of MBNL2 in each group (control, DM1, and DM2) for visualization and comparative analysis. PolyA-seq was performed on control, DM1, and DM2 autopsy-derived brain samples (hippocampus and frontal cortex, n=3) as well as wild-type (WT) and Mbnl1; Mbnl2 conditional double knockout (Mbnl1-/-; Mbnl2c/c; Nestin-Cre+/- or DKO) brain (n=3). Libraries were sequenced and the resultant files were processed and aligned to the reference genomes (hg19 and mm10). Further computational processing was performed to remove internal oligo(dT) mis-priming events, identify valid polyA sites, and trim to the exact polyA sites. BedGraph files were generated for each group in human (control, DM1, and DM2) and mouse (WT and Mbnl DKO) for comparative visualization.
2015-07-30 | E-GEOD-68890 | biostudies-arrayexpress
Project description:We sought to more precisely characterize the different alpha-synuclein (aSyn) 3M-bM-^@M-^YUTR mRNA species in normal and PD human brain. High-throughput, whole-transcriptome sequencing of the 3M-bM-^@M-^YUTR ends of polyadenylated mRNA transcripts (termed pA-RNAseq; see Methods) was performed on a cohort of 17 unaffected and 17 PD cerebral cortical tissue samples. This revealed 5 aSyn 3M-bM-^@M-^YUTR isoforms, with lengths of 290, 480, 560, 1070 and 2520 nt. Of these, the 560 nt and 2520 nt forms were predominant. The existence and relative preponderance of these species was further confirmed by Northern Blot. We next hypothesized, that aSyn 3M-bM-^@M-^YUTR selection might be altered in PD. Comparison of pA-RNAseq profiles from PD and unaffected cerebral cortex samples revealed an increase in the preponderance of the long 3M-bM-^@M-^YUTR species (>560 nt) relative to shorter species (<560 nt). Such a relative increase in aSynL was confirmed by Quantitative real-time RT-PCR (rt-qPCR) and appeared specific for PD, as the increase was also observed by comparison to RNA from amyotrophic lateral sclerosis patient samples. We note that the modified aSyn 3M-bM-^@M-^YUTR selection associated with PD patient tissue was detected in cerebral cortex tissue, which typically harbors pathological evidence of the disease process without frank cell loss; thus, this phenotype is unlikely to be a secondary consequence of neurodegeneration. Comparison of 3'UTR ends of alpha-synuclein in PD and unaffected brain cortex
2012-10-04 | E-GEOD-40710 | biostudies-arrayexpress
Project description:The tumorigenesis of small intestinal neuroendocrine tumors (NETs) is poorly understood. Recent studies have associated alternative polyadenylation with proliferation, cell transformation and cancer. Polyadenylation is the process in which the pre-mRNA is cleaved at a polyA site and a polyA tail is added. Genes with two or more polyA sites can undergo alternative polyadenylation. This produces two or more distinct mRNA isoforms with different 3M-bM-^@M-^Y untranslated regions. Additionally, alternative polyadenylation can also produce mRNAs containing different 3M-bM-^@M-^Y-terminal coding regions. Therefore, alternative polyadenylation alters both the repertoire and the expression level of proteins. Here we used high-throughput sequencing data to map polyA sites and characterize polyadenylation genome-wide in three small intestinal neuroendocrine tumors and a reference sample. In the tumors sixteen genes showed significant changes of alternative polyadenylation pattern, which lead to either the 3M-bM-^@M-^Y truncation of mRNA coding regions or 3M-bM-^@M-^Y untranslated regions. Among these, 11 genes had been previously associated with cancer, with 4 genes being known tumor suppressors: DCC, PDZD2, MAGI1 and DACT2. We validated the alternative polyadenylation in 3 out of 3 cases with Q-RT-PCR. Our findings suggest that changes of alternative polyadenylation pattern in these 16 genes could be involved in the tumorigenesis of small intestinal neuroendocrine tumors. Furthermore, they also point to alternative polyadenylation as a new target for both diagnostic and treatment of small intestinal neuroendocrine tumors. The identified genes with alternative polyadenylation specific to the small intestinal neuroendocrine tumors could be further tested as diagnostic markers and drug targets for disease prevention and treatment. PolyA-seq profiling of 3 human neuroendocrine tumors compared and pituitary using Direct RNA Sequencing from Helicos Biosciences Technology
2014-04-10 | E-GEOD-56657 | biostudies-arrayexpress
Project description:The 3' ends of most Drosophila melanogaster genes are poorly annotated or are determined by only a single EST or cDNA clone. To enhance the annotation of poly(A) site use in Drosophila, we performed deep sequencing on RNA isolated from 29 dissected tissues using an approach designed to enrich for poly(A) spanning reads. From these experiments, we identified 1.4 million poly(A) spanning reads leading to the identification of many new poly(A) sites and the identification of many tissue-specific poly(A) sites. For data usage terms and conditions, please refer to http://www.genome.gov/27528022 and http://www.genome.gov/Pages/Research/ENCODE/ENCODEDataReleasePolicyFinal2008.pdf RNA from 29 dissected Drosophila melanogaster tissues (in duplicate) were used to prepare polyA enriched RNA-Seq libraries. Briefly, total RNA was poly(A) selected, fragmented, and ligated to 5' and 3' RNA linkers. These libraries were amplified using Illumina paired-end primers, and subsequently reamplified using a 3' primer complementary to the 3' adapter but containing 6 Ts at the 3' end. The libraries were also multiplexed and up to 12 samples mixed per lane and sequenced on an Illumina GAIIx using paired-end 76 bp reads, or an illumina HiSeq 2000 using paired-end 100 bp reads. All reads were mapped to the Drosophila melanogaster genome to identify unmapped reads. Unmapped reads containing at least 10 A residues at the 3' end were identified, the terminal A residues trimmed, realigned to the genome to identify uniquely mapped reads. Such reads were identified as polyA spanning reads
2011-11-23 | E-GEOD-33905 | biostudies-arrayexpress
Project description:Our RNA-seq experiment (E-MTAB-5383) and subsequent in vitro analysis suggested that snoRD50a is a trans-acting RNA that can function in mRNA 3' processing. To test this hypothesis at transcriptomic level, we performed PAS-seq experiments to compare the polyA+ transcripts profile between negative control ASO (nc ASO) and snoRD50a ASO treated Hela cells, ASO (antisense oligonucletode) technique is a commonly used technique to deplete nuclear RNAs. snoRD50a ASO is used to deplete endogenous snoRD50a. See related experiments: E-MTAB-5383; E-MTAB-5384.
2017-05-31 | E-MTAB-5389 | biostudies-arrayexpress
Project description:Prostate adenocarcinoma and matched adjacent normal samples were profiled by deep transcriptional sequencing to analyze transcription-induced chimeras and gene fusions. Reference samples from the MAQC and brain and universal reference libraries were also sequenced. Two-condition experiment plus reference samples: Prostate adenocarcinoma versus matched normal from three separate patients, plus brain and universal reference samples from the MAQC project.
2010-09-23 | E-GEOD-24283 | biostudies-arrayexpress
Project description:A Parallel Analysis of six mouse tissues to investigate cleavage products and degradome Survey of polyA transcripts having 5' monophosphates (a characteristic of RNA cleavage events) in six mouse tissues.
2010-07-30 | E-GEOD-22627 | biostudies-arrayexpress
Project description:Long non-coding RNAs (lncRNAs) comprise a diverse class of transcripts that can regulate molecular and cellular processes in brain development and diseasee. LncRNAs exhibit cell type- and tissue-specific expression, but little is known about the expression and function of lncRNAs in the developing human brain. Here, we deeply profiled lncRNAs from polyadenylated and total RNA obtained from human neocortex at different stages of development and integrated this resource to analyze the transcriptomes of single cells. While lncRNAs were generally detected at low levels in whole tissues, single cell transcriptomics revealed that many lncRNAs are abundantly expressed in individual cells and are cell type-specific. Furthermore, we used CRISRPi to show that LOC646329, a lncRNA enriched in radial glia but detected at low abundance in tissues, regulates cell proliferation. The discrete and abundant expression of lncRNAs among individual cells has important implications for both their biological function and utility for distinguishing neural cell types. 16 Bulk Tissue Samples from GW13-23; 226 Single Cells from GW19.5-23.5 ------------------ bulk_tpm.polya.txt: bulk RNA-seq expression; using polyA full reference scell_ncounts.genes.thresh.txt: single cell RNA-seq expression; using polyA stringent reference; includes 50 GW16, GW21, GW21p3 cells previously analyzed (Pollen et. al. 2014) polya_RNA_stringent_ref.gtf: bulk RNA-seq polyA stringent transcriptome reference polya_RNA_full_ref.gtf: bulk RNA-seq polyA full transcriptome reference total_RNA_stringent_ref.gtf: bulk RNA-seq total stringent transcriptome reference total_RNA_full_ref.gtf: bulk RNA-seq total full transcriptome reference GW13_1_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW13_1_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW13_1_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW13_1_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW14.5_1_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW14.5_1_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW14.5_1_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW14.5_1_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW16_1_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW16_1_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW16_1_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW16_1_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW16_2_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW16_2_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW16_2_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW16_2_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW21_1_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW21_1_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW21_1_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW21_1_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW21_2_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW21_2_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW21_2_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW21_2_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW23_1_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW23_1_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW23_1_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW23_1_total_plus.bw: strand-specific bulk RNA-seq alignment signal GW23_2_polya_minus.bw: strand-specific bulk RNA-seq alignment signal GW23_2_polya_plus.bw: strand-specific bulk RNA-seq alignment signal GW23_2_total_minus.bw: strand-specific bulk RNA-seq alignment signal GW23_2_total_plus.bw: strand-specific bulk RNA-seq alignment signal
2016-03-28 | E-GEOD-71315 | biostudies-arrayexpress
Project description:RNA transcripts are distributed non-uniformly in the oocytes of many animals, such that newly-divided embryo cells (blastomeres) inherit distinct transcriptomes following fertilization. In animals such as the frog, Xenopus laevis, programmed transcript regionalization directs early embryonic axis formation and is essential for normal development. However, it is unknown whether such transcriptome asymmetry directs embryogenesis in mammals, or indeed, whether it occurs at all. We here address this by transcript profiling of matching sub-cellular structures and single cells in mouse oocytes and early embryos and analytical strategies exploiting the paired data structure. Spindle samples contained a set of transcripts that was distinguishable from that of the unfertilized metaphase II (mII) oocytes from which each spindle was microsurgically dissected. Immediately following fertilization, cytokinesis produces a 1-cell embryo (zygote) and associated spindle-enriched second polar body (Pb2) whose transcript profiles also differed one from the other, partially reflecting the enrichment of spindle-associated transcripts. Non-uniform transcript distribution within zygotes did not lead to programmed transcriptome asymmetry between the blastomeres of nascent two-cell embryos, or between the second mitotic products of 3-cell embryos. These findings suggest that mammalian oocytes and zygotes exhibit transcript regionalization without subsequent transcriptome asymmetries between respective early mitotic products. This contrasts the situation in Xenopus and places constraints on the ability of maternal transcriptomic prepatterning to prescribe early mammalian development. 16 Zygote vs. Polar body pairs, 10 spindle vs. oocyte pairs, 22 2-cell embryos and 8 3-cell embryos were analysed. Although the raw data are two channel, only the Cy5 signal of each file was analyzed. The Cy5 channel for each gene is normalized to the average Cy5 intensity of the gene across all samples.
2011-04-21 | E-GEOD-27396 | biostudies-arrayexpress
Project description:Soil Samples Collected from Along the DC-Metro Area of the Potomac River
2023-10-06 | E-MTAB-13365 | biostudies-arrayexpress