Expression evolution in gene groups: Inferring non-neutral regulatory change in pathways from transcriptional profiling data
Ontology highlight
ABSTRACT:
ORGANISM(S): Saccharomyces cerevisiae x Saccharomyces mikatae
SUBMITTER: Yulia Mostovoy
PROVIDER: E-GEOD-38875 | biostudies-arrayexpress |
REPOSITORIES: biostudies-arrayexpress
Publications
Inferring evolutionary histories of pathway regulation from transcriptional profiling data.
Schraiber Joshua G JG Mostovoy Yulia Y Hsu Tiffany Y TY Brem Rachel B RB
PLoS computational biology 20131010 10
One of the outstanding challenges in comparative genomics is to interpret the evolutionary importance of regulatory variation between species. Rigorous molecular evolution-based methods to infer evidence for natural selection from expression data are at a premium in the field, and to date, phylogenetic approaches have not been well-suited to address the question in the small sets of taxa profiled in standard surveys of gene expression. We have developed a strategy to infer evolutionary histories ...[more]
Similar Datasets
Project description:We profiled the strand-specific transcriptome of Saccharomyces cerevisiae lacking SWR1 and HTZ1 and compared it to the wild-type transcriptome 3'-end RNA-seq was performed on RNA from Saccharomyces cerevisiae wild-type and swr1 htz1 strains (two biological replicates each)
2016-06-01 | E-GEOD-46535 | biostudies-arrayexpress
Project description:Alternative cleavage and polyadenylation (APA) generates diverse mRNA isoforms. We developed 3' region extraction and deep sequencing (3'READS) to address mispriming issues that commonly plague poly(A) site (pA) identification, and we used the method to comprehensively map pAs in the mouse genome. Thorough annotation of gene 3' ends revealed over 5,000 previously overlooked pAs (~8% of total) flanked by A-rich sequences, underscoring the necessity of using an accurate tool for pA mapping. About 79% of mRNA genes and 66% of long noncoding RNA genes undergo APA, but these two gene types have distinct usage patterns for pAs in introns and upstream exons. Quantitative analysis of APA isoforms by 3'READS indicated that promoter-distal pAs, regardless of intron or exon locations, become more abundant during embryonic development and cell differentiation and that upregulated isoforms have stronger pAs, suggesting global modulation of the 3' endM-bM-^@M-^Sprocessing activity in development and differentiation. 3'READS to map pAs in mouse genome
2012-12-17 | E-GEOD-42398 | biostudies-arrayexpress
Project description:Histone modifications are important markers of function and chromatin state, yet the DNA elements that direct them to specific locations in the genome are poorly understood. Here we use the genetic variation in Yoruba lymphoblastoid cell lines as a natural experiment to identify genetic differences that affect histone marks and to better understand their relationship with transcriptional regulation. Across the genome, we identified hundreds of quantitative trait loci that impact histone modification or RNA polymerase (PolII) occupancy. In many cases the same variant is associated with quantitative changes in multiple histone marks and PolII, as well as in DNaseI sensitivity and nucleosome positioning, indicating that these molecular phenotypes often share a single underlying genetic cause. Variants that impact chromatin at distal regulatory sites frequently also direct changes in chromatin and gene expression at associated promoters; while most of these distal regulators enhance promoter activity, some act as distal chromatin silencers. Finally, we find that polymorphisms in transcription factor binding sites are often causally responsible for variation in local histone modification. In summary, the class of variants identified here generate coordinated changes in chromatin both locally and sometimes at distant locations, frequently drive changes in gene expression, and likely play an important role in the genetics of complex traits. ChIP-seq of RNA Polymerase II and 4 histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K27me3) in 10 unrelated Yoruba HapMap lymphoblastoid cell lines
2013-10-01 | E-GEOD-47991 | biostudies-arrayexpress
Project description:This project aims to investigate the metabolic pathways expressed by the active microbial community occurring at the deep continental subsurface. Subsurface chemoLithoautotrophic Microbial Ecosystems (SLiMEs) under oligotrophic conditions are supported by H2; however, the overall ecological trophic structures of these communities are poorly understood. Some deep, fluid-filled fractures in the Witwatersrand Basin, South Africa appear to support inverted trophic pyramids wherein methanogens contributing <5% of the total DNA apparently produce CH4 that supports the rest of the community. Here we show the active metabolic relationships of one such trophic structure by combining metatranscriptomic assemblies, metaproteomic and stable isotopic data, and thermodynamic modeling. Four autotrophic β-proteobacteria genera that are capable of oxidizing sulfur by denitrification dominate. They co-occur with sulfate reducers, anaerobic methane oxidizers and methanogens, which each comprises <5% of the total community. Defining trophic levels of microbial chemolithoautotrophs by the number of transfers from the initial abiotic H2-driven CO2 fixation, we propose a top-down cascade influence of the metabolic consumers that enhances the fitness of the metabolic producers to explain the inverted biomass pyramid of a multitrophic SLiME. Symbiotic partnerships are pivotal in the deep biosphere on and potentially beyond the Earth.
2018-10-27 | PXD004634 | Pride
Project description:The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3' untranslated regions (3'UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-Seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-Seq reads substantially increased and improved existing 3'UTR annotations, resulting in confidently identified 3'UTRs for more than 78.79% of the annotated protein-coding genes in zebrafish. Most zebrafish genes undergo alternative CPA with more than a thousand genes using different dominant 3'UTRs at different stages. 3'UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3'UTRs are highly expressed in the ovary yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3'UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At two hours post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3'UTRs provide a resource for studying gene regulation during vertebrate development. 3P-Seq was used to map the 3' ends of protein-coding genes in the zebrafish genome
2012-06-17 | E-GEOD-37453 | 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 liver X receptors (LXRs) are nuclear receptors that form permissive heterodimers with retinoid X receptor (RXR) and are important regulators of lipid metabolism in the liver. We have recently shown that RXR agonist-induced hypertriglyceridemia and hepatic steatosis in mice is dependent on LXR and correlates with an LXR-dependent hepatic induction of lipogenic genes. To further investigate the role of RXR and LXR in the regulation of hepatic gene expression, we have mapped the ligand-regulated genome-wide binding of these factors in mouse liver. We find that the RXR agonist bexarotene primarily increases the genomic binding of RXR, whereas the LXR agonist T0901317 greatly increases both LXR and RXR binding. Functional annotation of putative direct LXR target genes revealed a significant association with classical LXR-regulated pathways as well as PPAR signaling pathways, and subsequent ChIP-seq mapping of PPARM-NM-1 binding demonstrated binding of PPARM-NM-1 to 71-88% of the identified LXR:RXR binding sites. Sequence analysis of shared binding regions combined with sequential ChIP on selected sites indicate that LXR:RXR and PPARM-NM-1:RXR bind to degenerate response elements in a mutually exclusive manner. Together our findings suggest extensive and unexpected cross-talk between hepatic LXR and PPARM-NM-1 at the level of binding to shared genomic sites LXR, RXR, PPARalpha and RNA Polymerase II ChIP-seq on livers from female C57BL/6 wild-type and/or LXRM-NM-1/M-NM-2-deficient mice (13 weeks of age, n=1) treated by oral gavage once daily for 14 days with the RXR agonist bexarotene (100 mg/kg body weight [mpk], in 1% carboxymethylcellulose), the LXR agonist T0901317 (T09, 30 mpk) or vehicle alone.
2012-01-31 | E-GEOD-35262 | biostudies-arrayexpress
Project description:Sequencing of the 3â end of poly(A)+ RNA identifies cleavage and polyadenylation sites (pAs) and measures transcript expression. We previously developed a method, 3â region extraction and deep sequencing (3âREADS), to address mispriming issues that often plague 3â end sequencing. Here we report a new version, named 3âREADS+, which has vastly improved accuracy and sensitivity. Using a special locked nucleic acid oligo to capture poly(A)+ RNA and to remove bulk of the poly(A) tail, 3âREADS+ generates RNA fragments with an optimal number of terminal As that balance data quality and detection of genuine pAs. With improved RNA ligation steps for efficiency, the method shows much higher sensitivity (over two orders of magnitude) compared to the previous version. Using 3âREADS+, we have uncovered a sizable fraction of previously overlooked pAs located next to or within a stretch of adenylate residues in human genes, and more accurately assessed the frequency of alternative cleavage and polyadenylation (APA) in HeLa cells (~50%). 3âREADS+ will be a useful tool to accurately study APA and to analyze gene expression by 3â end counting, especially when the amount of input total RNA is limited. Nine 3'READS+ libraries were made with different amounts (100 ng, 200 ng, 400 ng, 1000 ng, 5000 ng, 15000 ng) of input Hela RNA.
2016-08-14 | E-GEOD-84170 | biostudies-arrayexpress
Project description:Analysis of bacterial fraction collected on GF/F filters post pre-filtration on 1um filter. 15L were filtered from Bering Strait (BSt) surface water and Chukchi Sea (station 2) bottom waters.
2019-01-25 | PXD006472 | Pride
Project description:Alternative polyadenylation is an important cellular mechanism that enables generation of mRNA isoforms that differ in their 3' untranslated regions (3' UTRs) and consequently in their susceptibility to miRNA and RNA binding protein mediated regulation. A dramatic change in polyadenylation site usage, leading to the systematic expression of short 3’ UTR isoforms is known to occur upon induction of proliferation in resting cells. To understand the functional consequences of short 3’ UTR isoform expression we used 3' end sequencing and quantitative mass spectroscopy to determine polyadenylation site use, mRNA and protein levels in murine and human naive and activated T cells. We found that while the process and its impact on the susceptibility to miRNA and RNA binding protein mediated regulation are evolutionarily conserved, the conservation is poor at the level of individual orthologous genes. Contrary to the common belief, we did not find that transcriptome-wide 3' UTR shortening leads to a matched increase in mRNA and protein levels of genes with tandem polyadenylation sites. 3' ends of transcripts were profiled by high-throughput sequencing in murine and human naive and activated T cells.
2014-09-25 | E-GEOD-54950 | biostudies-arrayexpress