The Contribution of RNA Decay Quantitative Trait Loci to Inter-Individual Variation in Steady-State Gene Expression Levels

Recent gene expression QTL (eQTL) mapping studies have provided considerable insight into the genetic basis for inter-individual regulatory variation. However, a limitation of all eQTL studies to date, which have used measurements of steady-state gene expression levels, is the inability to directly distinguish between variation in transcription and decay rates. To address this gap, we performed a genome-wide study of variation in gene-specific mRNA decay rates across individuals. Using a time-course study design, we estimated mRNA decay rates for over 16,000 genes in 70 Yoruban HapMap lymphoblastoid cell lines (LCLs), for which extensive genotyping data are available. Considering mRNA decay rates across genes, we found that: (i) as expected, highly expressed genes are generally associated with lower mRNA decay rates, (ii) genes with rapid mRNA decay rates are enriched with putative binding sites for miRNA and RNA binding proteins, and (iii) genes with similar functional roles tend to exhibit correlated rates of mRNA decay. Focusing on variation in mRNA decay across individuals, we estimate that steady-state expression levels are significantly correlated with variation in decay rates in 10% of genes. Somewhat counter-intuitively, for about half of these genes, higher expression is associated with faster decay rates, possibly due to a coupling of mRNA decay with transcriptional processes in genes involved in rapid cellular responses. Finally, we used these data to map genetic variation that is specifically associated with variation in mRNA decay rates across individuals. We found 195 such loci, which we named RNA decay quantitative trait loci (“rdQTLs”). All the observed rdQTLs are located near the regulated genes and therefore are assumed to act in cis. By analyzing our data within the context of known steady-state eQTLs, we estimate that a substantial fraction of eQTLs are associated with inter-individual variation in mRNA decay rates.

[1]  B. Weir,et al.  The quantitative genetics of transcription. , 2005, Trends in genetics : TIG.

[2]  L. Liang,et al.  A genome-wide association study of global gene expression , 2007, Nature Genetics.

[3]  Pan Du,et al.  lumi: a pipeline for processing Illumina microarray , 2008, Bioinform..

[4]  Joseph K. Pickrell,et al.  DNaseI sensitivity QTLs are a major determinant of human expression variation , 2011, Nature.

[5]  Eran Segal,et al.  Transient transcriptional responses to stress are generated by opposing effects of mRNA production and degradation , 2008, Molecular systems biology.

[6]  C. Reilly,et al.  Genome-wide analysis of mRNA decay in resting and activated primary human T lymphocytes. , 2002, Nucleic acids research.

[7]  Jacek Majewski,et al.  The study of eQTL variations by RNA-seq: from SNPs to phenotypes. , 2011, Trends in genetics : TIG.

[8]  John D. Storey,et al.  Genetic interactions between polymorphisms that affect gene expression in yeast , 2005, Nature.

[9]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[10]  M. Choder,et al.  Promoter Elements Regulate Cytoplasmic mRNA Decay , 2011, Cell.

[11]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[12]  Sophie Janssens,et al.  Ubiquitin: tool and target for intracellular NF-kappaB inhibitors. , 2006, Trends in immunology.

[13]  D. Price,et al.  Flavopiridol Inactivates P-TEFb and Blocks Most RNA Polymerase II Transcription in Vivo * , 2001, The Journal of Biological Chemistry.

[14]  Toshihiro Tanaka.  The International HapMap Project , 2003, Nature.

[15]  N. Friedman,et al.  Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells , 2011, Nature Biotechnology.

[16]  Jörg Hackermüller,et al.  mRNA Openers and Closers: Modulating AU‐Rich Element‐Controlled mRNA Stability by a Molecular Switch in mRNA Secondary Structure , 2004, Chembiochem : a European journal of chemical biology.

[17]  Yoav Gilad,et al.  The effects of EBV transformation on gene expression levels and methylation profiles. , 2011, Human molecular genetics.

[18]  M. Culbertson,et al.  Looking at mRNA decay pathways through the window of molecular evolution. , 2003, Current opinion in genetics & development.

[19]  O. Bensaude,et al.  The Transcriptional Inhibitors, Actinomycin D and α-Amanitin, Activate the HIV-1 Promoter and Favor Phosphorylation of the RNA Polymerase II C-terminal Domain* , 1999, The Journal of Biological Chemistry.

[20]  Wenqian Hu,et al.  Co-translational mRNA decay in Saccharomyces cerevisiae , 2009, Nature.

[21]  M. Stephens,et al.  High-Resolution Mapping of Expression-QTLs Yields Insight into Human Gene Regulation , 2008, PLoS genetics.

[22]  M. Gerstein,et al.  Variation in Transcription Factor Binding Among Humans , 2010, Science.

[23]  A. Millar,et al.  Genome-Wide Analysis of mRNA Decay Rates and Their Determinants in Arabidopsis thaliana[W] , 2007, The Plant Cell Online.

[24]  D. Galson,et al.  Resveratrol triggers the pro-apoptotic endoplasmic reticulum stress response and represses pro-survival XBP1 signaling in human multiple myeloma cells. , 2011, Experimental hematology.

[25]  J. Cherry,et al.  Identification of unstable transcripts in Arabidopsis by cDNA microarray analysis: Rapid decay is associated with a group of touch- and specific clock-controlled genes , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[26]  F. Bruggeman,et al.  Introduction to systems biology. , 2007, EXS.

[27]  Daniel Rios,et al.  Ensembl 2011 , 2010, Nucleic Acids Res..

[28]  L. Lim,et al.  MicroRNA targeting specificity in mammals: determinants beyond seed pairing. , 2007, Molecular cell.

[29]  Joshua T. Burdick,et al.  Mapping determinants of human gene expression by regional and genome-wide association , 2005, Nature.

[30]  Y. Pilpel,et al.  Transcriptome Kinetics Is Governed by a Genome-Wide Coupling of mRNA Production and Degradation: A Role for RNA Pol II , 2011, PLoS genetics.

[31]  T. Pastinen.  Genome-wide allele-specific analysis: insights into regulatory variation , 2010, Nature Reviews Genetics.

[32]  K. Hansen,et al.  Genome-Wide Identification of Alternative Splice Forms Down-Regulated by Nonsense-Mediated mRNA Decay in Drosophila , 2009, PLoS genetics.

[33]  John D. Storey,et al.  Precision and functional specificity in mRNA decay , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Arkady B. Khodursky,et al.  Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[35]  Gonçalo R. Abecasis,et al.  Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma , 2007, Nature.

[36]  M. Loda,et al.  Fatty acid synthase as a potential therapeutic target in cancer. , 2010, Future oncology.

[37]  Joseph K. Pickrell,et al.  Understanding mechanisms underlying human gene expression variation with RNA sequencing , 2010, Nature.

[38]  Jeffrey Wilusz,et al.  The highways and byways of mRNA decay , 2007, Nature Reviews Molecular Cell Biology.

[39]  C. Burge,et al.  Most mammalian mRNAs are conserved targets of microRNAs. , 2008, Genome research.

[40]  M. Choder,et al.  RNA Polymerase II Subunits Link Transcription and mRNA Decay to Translation , 2010, Cell.

[41]  C. Burge,et al.  Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets , 2005, Cell.

[42]  J. Wilusz,et al.  Bringing the role of mRNA decay in the control of gene expression into focus. , 2004, Trends in genetics : TIG.

[43]  S. Nanko,et al.  Impaired feedback regulation of XBP1 as a genetic risk factor for bipolar disorder , 2003, Nature Genetics.

[44]  Yongtao Guan,et al.  Practical Issues in Imputation-Based Association Mapping , 2008, PLoS genetics.

[45]  R. Guigó,et al.  Transcriptome genetics using second generation sequencing in a Caucasian population , 2010, Nature.

[46]  T. Hudson,et al.  Influence of human genome polymorphism on gene expression. , 2006, Human molecular genetics.

[47]  Daniel R. Larson,et al.  Single-Molecule mRNA Decay Measurements Reveal Promoter- Regulated mRNA Stability in Yeast , 2011, Cell.

[48]  John D. Storey,et al.  Statistical significance for genomewide studies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[49]  A. Gartel,et al.  A novel transcriptional inhibitor induces apoptosis in tumor cells and exhibits antiangiogenic activity. , 2006, Cancer research.

[50]  Sophie Janssens,et al.  Ubiquitin: tool and target for intracellular NF-kappaB inhibitors. , 2006, Trends in immunology.

[51]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[52]  M. W. Foster,et al.  Integrating ethics and science in the International HapMap Project , 2004, Nature Reviews Genetics.

[53]  Scott A. Rifkin,et al.  Revealing the architecture of gene regulation: the promise of eQTL studies. , 2008, Trends in genetics : TIG.

[54]  Matthew Stephens,et al.  Dissecting the regulatory architecture of gene expression QTLs , 2012, Genome Biology.

[55]  Zhaohui S. Qin,et al.  A second generation human haplotype map of over 3.1 million SNPs , 2007, Nature.

[56]  A. Khodursky,et al.  Life after transcription--revisiting the fate of messenger RNA. , 2003, Trends in genetics : TIG.

[57]  G. Churchill,et al.  Variation in gene expression within and among natural populations , 2002, Nature Genetics.

[58]  H. M. Sobell.  Actinomycin and DNA transcription. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[59]  S. Peltz,et al.  The cap-to-tail guide to mRNA turnover , 2001, Nature Reviews Molecular Cell Biology.

[60]  Doron Betel,et al.  The microRNA.org resource: targets and expression , 2007, Nucleic Acids Res..

[61]  N. Slonim,et al.  A universal framework for regulatory element discovery across all genomes and data types. , 2007, Molecular cell.

[62]  Paul Scheet,et al.  A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. , 2006, American journal of human genetics.

[63]  C. Friedel,et al.  Genome-Wide Assessment of AU-Rich Elements by the AREScore Algorithm , 2012, PLoS genetics.

[64]  T. Kawamoto,et al.  Expression of the gene for Dec2, a basic helix-loop-helix transcription factor, is regulated by a molecular clock system. , 2004, The Biochemical journal.

[65]  Life Technologies,et al.  A map of human genome variation from population-scale sequencing , 2011 .

[66]  Matthew S. Sachs,et al.  Early nonsense: mRNA decay solves a translational problem , 2006, Nature Reviews Molecular Cell Biology.

[67]  Y. Pilpel,et al.  Regulatory mechanisms and networks couple the different phases of gene expression. , 2011, Trends in genetics : TIG.

[68]  M. Magnasco,et al.  Decay rates of human mRNAs: correlation with functional characteristics and sequence attributes. , 2003, Genome research.

[69]  John D. Storey,et al.  Mapping the Genetic Architecture of Gene Expression in Human Liver , 2008, PLoS biology.

[70]  John C. Marioni,et al.  Gene Expression Differences Among Primates Are Associated With Changes in a Histone Epigenetic Modification , 2011, Genetics.

[71]  Uwe Ohler,et al.  Integrative regulatory mapping indicates that the RNA-binding protein HuR couples pre-mRNA processing and mRNA stability. , 2011, Molecular cell.

[72]  J. Ross,et al.  mRNA stability in mammalian cells. , 1995, Microbiological reviews.

[73]  Hong Sun,et al.  Stra13 expression is associated with growth arrest and represses transcription through histone deacetylase (HDAC)-dependent and HDAC-independent mechanisms. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[74]  D. Hartl,et al.  Dominance and the evolutionary accumulation of cis- and trans-effects on gene expression , 2008, Proceedings of the National Academy of Sciences.

[75]  I. Tirosh,et al.  Coupled Evolution of Transcription and mRNA Degradation , 2011, PLoS biology.

[76]  K. Gunsalus,et al.  Combinatorial microRNA target predictions , 2005, Nature Genetics.

[77]  D. Altshuler,et al.  A map of human genome variation from population-scale sequencing , 2010, Nature.

[78]  H. Stefánsson,et al.  Genetics of gene expression and its effect on disease , 2008, Nature.

[79]  M. Gorospe,et al.  Stability Regulation of mRNA and the Control of Gene Expression , 2005, Annals of the New York Academy of Sciences.

[80]  J. LaMarre,et al.  AU-rich elements and the control of gene expression through regulated mRNA stability , 2004, Animal Health Research Reviews.

[81]  J. Pérez-Ortín.  Genomics of mRNA turnover. , 2008, Briefings in functional genomics & proteomics.

[82]  L. Maquat,et al.  Regulation of cytoplasmic mRNA decay , 2012, Nature Reviews Genetics.

[83]  Andreas Rosenwald,et al.  Genomic-scale measurement of mRNA turnover and the mechanisms of action of the anti-cancer drug flavopiridol , 2001, Genome Biology.

[84]  Gregory D. Schuler,et al.  Database resources of the National Center for Biotechnology Information: update , 2004, Nucleic acids research.

[85]  A. Bertolotti,et al.  Exploiting the selectivity of protein phosphatase 1 for pharmacological intervention , 2013, The FEBS journal.

[86]  Colin N. Dewey,et al.  A Genome-Wide Map of Conserved MicroRNA Targets in C. elegans , 2006, Current Biology.

[87]  Alexei A. Sharov,et al.  Database for mRNA Half-Life of 19 977 Genes Obtained by DNA Microarray Analysis of Pluripotent and Differentiating Mouse Embryonic Stem Cells , 2008, DNA research : an international journal for rapid publication of reports on genes and genomes.

[88]  Christina Backes,et al.  GeneTrail—advanced gene set enrichment analysis , 2007, Nucleic Acids Res..

[89]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[90]  R. Stoughton,et al.  Genetics of gene expression surveyed in maize, mouse and man , 2003, Nature.

[91]  Joel S. Parker,et al.  Novel mRNA Targets for Tristetraprolin (TTP) Identified by Global Analysis of Stabilized Transcripts in TTP-Deficient Fibroblasts , 2006, Molecular and Cellular Biology.

[92]  Emilie Lalonde,et al.  RNA sequencing reveals the role of splicing polymorphisms in regulating human gene expression. , 2011, Genome research.