Whole transcriptome RNA-Seq allelic expression in human brain

BackgroundMeasuring allelic RNA expression ratios is a powerful approach for detecting cis-acting regulatory variants, RNA editing, loss of heterozygosity in cancer, copy number variation, and allele-specific epigenetic gene silencing. Whole transcriptome RNA sequencing (RNA-Seq) has emerged as a genome-wide tool for identifying allelic expression imbalance (AEI), but numerous factors bias allelic RNA ratio measurements. Here, we compare RNA-Seq allelic ratios measured in nine different human brain regions with a highly sensitive and accurate SNaPshot measure of allelic RNA ratios, identifying factors affecting reliable allelic ratio measurement. Accounting for these factors, we subsequently surveyed the variability of RNA editing across brain regions and across individuals.ResultsWe find that RNA-Seq allelic ratios from standard alignment methods correlate poorly with SNaPshot, but applying alternative alignment strategies and correcting for observed biases significantly improves correlations. Deploying these methods on a transcriptome-wide basis in nine brain regions from a single individual, we identified genes with AEI across all regions (SLC1A3, NHP2L1) and many others with region-specific AEI. In dorsolateral prefrontal cortex (DLPFC) tissues from 14 individuals, we found evidence for frequent regulatory variants affecting RNA expression in tens to hundreds of genes, depending on stringency for assigning AEI. Further, we find that the extent and variability of RNA editing is similar across brain regions and across individuals.ConclusionsThese results identify critical factors affecting allelic ratios measured by RNA-Seq and provide a foundation for using this technology to screen allelic RNA expression on a transcriptome-wide basis. Using this technology as a screening tool reveals tens to hundreds of genes harboring frequent functional variants affecting RNA expression in the human brain. With respect to RNA editing, the similarities within and between individuals leads us to conclude that this post-transcriptional process is under heavy regulatory influence to maintain an optimal degree of editing for normal biological function.

[1]  W. Sadee,et al.  Searching for polymorphisms that affect gene expression and mRNA processing: Example ABCB1 (MDR1) , 2006, The AAPS Journal.

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

[3]  J. Leek,et al.  Temporal dynamics and genetic control of transcription in the human prefrontal cortex , 2011, Nature.

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

[5]  Ruijie Liu,et al.  Allele-specific expression analysis methods for high-density SNP microarray data , 2012, Bioinform..

[6]  D. Stephan,et al.  A survey of genetic human cortical gene expression , 2007, Nature Genetics.

[7]  Wolfgang Viechtbauer,et al.  Conducting Meta-Analyses in R with the metafor Package , 2010 .

[8]  Jacek Majewski,et al.  Genome-wide analysis of transcript isoform variation in humans , 2008, Nature Genetics.

[9]  David I. Smith,et al.  Tumor Transcriptome Sequencing Reveals Allelic Expression Imbalances Associated with Copy Number Alterations , 2010, PloS one.

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

[11]  T. Mikkelsen,et al.  Altered adenosine-to-inosine RNA editing in human cancer. , 2007, Genome research.

[12]  Stanley F. Nelson,et al.  Identification of allele-specific alternative mRNA processing via transcriptome sequencing , 2012, Nucleic acids research.

[13]  L. Saksida,et al.  Assessment of Glutamate Transporter GLAST (EAAT1)-Deficient Mice for Phenotypes Relevant to the Negative and Executive/Cognitive Symptoms of Schizophrenia , 2009, Neuropsychopharmacology.

[14]  A. Malhotra,et al.  Dopamine Transporter Gene Variant Affecting Expression in Human Brain is Associated with Bipolar Disorder , 2011, Neuropsychopharmacology.

[15]  Daniel A. Skelly,et al.  A powerful and flexible statistical framework for testing hypotheses of allele-specific gene expression from RNA-seq data. , 2011, Genome research.

[16]  Mats Ensterö,et al.  Large-scale mRNA sequencing determines global regulation of RNA editing during brain development. , 2009, Genome research.

[17]  E. Levanon,et al.  Consistent levels of A-to-I RNA editing across individuals in coding sequences and non-conserved Alu repeats , 2010, BMC Genomics.

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

[19]  Mee Hyang Ko,et al.  Bioimaging of the unbalanced expression of microRNA9 and microRNA9* during the neuronal differentiation of P19 cells , 2008, The FEBS journal.

[20]  W. Sadee,et al.  Allelic expression of serotonin transporter (SERT) mRNA in human pons: lack of correlation with the polymorphism SERTLPR , 2006, Molecular Psychiatry.

[21]  Andrew D. Johnson,et al.  Allelic Expression Imbalance of Human mu Opioid Receptor (OPRM1) Caused by Variant A118G* , 2005, Journal of Biological Chemistry.

[22]  Shirley A. Miller,et al.  A simple salting out procedure for extracting DNA from human nucleated cells. , 1988, Nucleic acids research.

[23]  R. Baloh,et al.  Mutation in the glutamate transporter EAAT1 causes episodic ataxia, hemiplegia, and seizures , 2005, Neurology.

[24]  W. Sadee,et al.  Tryptophan hydroxylase 2 (TPH2) haplotypes predict levels of TPH2 mRNA expression in human pons , 2007, Molecular Psychiatry.

[25]  Mathieu Blanchette,et al.  Global patterns of cis variation in human cells revealed by high-density allelic expression analysis , 2009, Nature Genetics.

[26]  D. Clayton,et al.  Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing , 2009, Human molecular genetics.

[27]  Joaquín Goñi,et al.  Allele-Specific Gene Expression Is Widespread Across the Genome and Biological Processes , 2009, PloS one.

[28]  Ryan M. Smith,et al.  Nicotinic α5 receptor subunit mRNA expression is associated with distant 5′ upstream polymorphisms , 2011, European Journal of Human Genetics.

[29]  Identifying breast cancer risk loci by global differential allele-specific expression (DASE) analysis in mammary epithelial transcriptome , 2012, BMC Genomics.

[30]  Wenwei Zhang,et al.  Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome , 2012, Nature Biotechnology.

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

[32]  Jian Wang,et al.  Next-generation DNA sequencing-based assay for measuring allelic expression imbalance (AEI) of candidate neuropsychiatric disorder genes in human brain , 2011, BMC Genomics.

[33]  Leonardo Fazio,et al.  Polymorphisms in human dopamine D2 receptor gene affect gene expression, splicing, and neuronal activity during working memory , 2007, Proceedings of the National Academy of Sciences.

[34]  N. Alon,et al.  Monoallelic expression determines oncogenic progression and outcome in benign and malignant brain tumors. , 2012, Cancer research.

[35]  Michael Krawczak,et al.  Statistical inference of allelic imbalance from transcriptome data , 2011, Human mutation.

[36]  M. Gerstein,et al.  AlleleSeq: analysis of allele-specific expression and binding in a network framework , 2011, Molecular systems biology.

[37]  Kasia Wilczek,et al.  Role of miRNAs in Neuronal Differentiation from Human Embryonic Stem Cell—Derived Neural Stem Cells , 2012, Stem Cell Reviews and Reports.

[38]  Ryan M. Smith,et al.  Multiple Regulatory Variants Modulate Expression of 5-Hydroxytryptamine 2A Receptors in Human Cortex , 2013, Biological Psychiatry.

[39]  W. Sadee,et al.  Allelic mRNA expression of X-linked monoamine oxidase a (MAOA) in human brain: dissection of epigenetic and genetic factors. , 2006, Human molecular genetics.

[40]  Andrew D. Johnson,et al.  Polymorphisms affecting gene transcription and mRNA processing in pharmacogenetic candidate genes: detection through allelic expression imbalance in human target tissues , 2008, Pharmacogenetics and genomics.

[41]  Andrew D. Johnson,et al.  RNA Structures Affected By Single Nucleotide Polymorphisms In Transcribed Regions Of The Human Genome , 2011 .

[42]  Liang Chen.  Characterization and comparison of human nuclear and cytosolic editomes , 2013, Proceedings of the National Academy of Sciences.

[43]  John C. Marioni,et al.  Effect of read-mapping biases on detecting allele-specific expression from RNA-sequencing data , 2009, Bioinform..

[44]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[45]  Doron Lancet,et al.  Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification , 2005, Bioinform..

[46]  J. Rinn,et al.  Modular regulatory principles of large non-coding RNAs , 2012, Nature.

[47]  K. Gunderson,et al.  Genome-wide assessment of imprinted expression in human cells , 2011, Genome Biology.

[48]  P. Sullivan,et al.  Quantitative Analysis of Focused A-To-I RNA Editing Sites by Ultra-High-Throughput Sequencing in Psychiatric Disorders , 2012, PloS one.

[49]  Luigi Ferrucci,et al.  Abundant Quantitative Trait Loci Exist for DNA Methylation and Gene Expression in Human Brain , 2010, PLoS genetics.

[50]  Pavel V. Baranov,et al.  Darned in 2013: inclusion of model organisms and linking with Wikipedia , 2012, Nucleic Acids Res..

[51]  J. Shendure,et al.  Global survey of escape from X inactivation by RNA-sequencing in mouse. , 2010, Genome research.

[52]  C. Molony,et al.  Genetic analysis of genome-wide variation in human gene expression , 2004, Nature.

[53]  J. Reifman,et al.  A new strategy to reduce allelic bias in RNA-Seq readmapping , 2012, Nucleic acids research.

[54]  W. Sadee,et al.  Human N-acetyltransferase 1 *10 and *11 alleles increase protein expression through distinct mechanisms and associate with sulfamethoxazole-induced hypersensitivity , 2011, Pharmacogenetics and genomics.

[55]  P. Hevezi,et al.  Gene expression analyses reveal molecular relationships among 20 regions of the human CNS , 2006, Neurogenetics.

[56]  Jae-Hyung Lee,et al.  Accurate identification of A-to-I RNA editing in human by transcriptome sequencing. , 2012, Genome research.