Sporadic Autism Exomes Reveal a Highly Interconnected Protein Network of De Novo Mutations

Total Page:16

File Type:pdf, Size:1020Kb

Sporadic Autism Exomes Reveal a Highly Interconnected Protein Network of De Novo Mutations LETTER doi:10.1038/nature10989 Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations Brian J. O’Roak1,LauraVives1, Santhosh Girirajan1,EmreKarakoc1, Niklas Krumm1,BradleyP.Coe1,RoieLevy1,ArthurKo1,CholiLee1, Joshua D. Smith1, Emily H. Turner1, Ian B. Stanaway1, Benjamin Vernot1, Maika Malig1, Carl Baker1, Beau Reilly2,JoshuaM.Akey1, Elhanan Borenstein1,3,4,MarkJ.Rieder1, Deborah A. Nickerson1, Raphael Bernier2, Jay Shendure1 &EvanE.Eichler1,5 It is well established that autism spectrum disorders (ASD) have a per generation, in close agreement with our previous observations4, strong genetic component; however, for at least 70% of cases, the yet in general, higher than previous studies, indicating increased underlying genetic cause is unknown1. Under the hypothesis that sensitivity (Supplementary Table 2 and Supplementary Table 4)7. de novo mutations underlie a substantial fraction of the risk for We also observed complex classes of de novo mutation including: five developing ASD in families with no previous history of ASD or cases of multiple mutations in close proximity; two events consistent related phenotypes—so-called sporadic or simplex families2,3—we with paternal germline mosaicism (that is, where both siblings con- sequenced all coding regions of the genome (the exome) for tained a de novo event observed in neither parent); and nine events parent–child trios exhibiting sporadic ASD, including 189 new showing a weak minor allele profile consistent with somatic mosaicism trios and 20 that were previously reported4. Additionally, we also (Supplementary Table 3 and Supplementary Figs 2 and 3). sequenced the exomes of 50 unaffected siblings corresponding to Of the severe de novo events, 28% (33 of 120) are predicted to these new (n 5 31) and previously reported trios (n 5 19)4, for a truncate the protein. The distribution of synonymous, missense and total of 677 individual exomes from 209 families. Here we show nonsense changes corresponds well with a random mutation model7 that de novo point mutations are overwhelmingly paternal in (Supplementary Fig. 4 and Supplementary Table 2). However, the origin (4:1 bias) and positively correlated with paternal age, con- difference in nonsense rates between de novo and rare singleton events sistent with the modest increased risk for children of older fathers (not present in 1,779 other exomes) is striking (4:1) and suggests to develop ASD5. Moreover, 39% (49 of 126) of the most severe or strong selection against new nonsense events (Fisher’s exact test, disruptive de novo mutations map to a highly interconnected P , 0.0001). In contrast with a recent report8, we find no significant b-catenin/chromatin remodelling protein network ranked signifi- difference in mutation rate between affected and unaffected indivi- cantly for autism candidate genes. In proband exomes, recurrent duals; however, we do observe a trend towards increased non- protein-altering mutations were observed in two genes: CHD8 and synonymous rates in probands, consistent with the findings of ref. 9 NTNG1. Mutation screening of six candidate genes in 1,703 ASD (Supplementary Tables 1 and 2). probands identified additional de novo, protein-altering muta- Given the association of ASD with increased paternal age5 and our tions in GRIN2B, LAMC3 and SCN1A. Combined with copy previous observations4, we used molecular cloning, read-pair informa- number variant (CNV) data, these results indicate extreme locus tion, and obligate carrier status to identify informative markers linked heterogeneity but also provide a target for future discovery, to 51 de novo events and observed a marked paternal bias (41:10; diagnostics and therapeutics. binomial P , 1.4 3 1025; Fig. 1a and Supplementary Tables 3 and 5). We selected 189 autism trios from the Simons Simplex Collection This provides strong direct evidence that the germline mutation rate in (SSC)6, which included males significantly impaired with autism and protein-coding regions is, on average, substantially higher in males. A intellectual disability (n 5 47), a female sample set (n 5 56) of which similar finding was recently reported for de novo CNVs10. In addition, 26 were cognitively impaired, and samples chosen at random from the we observe that the number of de novo events is positively correlated remaining males in the collection (n 5 86) (Supplementary Table 1 with increasing paternal age (Spearman’s rank correlation 5 0.19; and Supplementary Fig. 1). In general, we excluded samples known to P , 0.008; Fig. 1b). Together, these observations are consistent with carry large de novo CNVs2. Exome sequencing was performed as the hypothesis that the modest increased risk for children of older described previously4, but with an expanded target definition (see fathers to develop ASD5 is the result of an increased mutation rate. Methods). We achieved sufficient coverage for both parents and child Using sequence read-depth methods in 122 of the 189 families, we to call genotypes for, on average, 29.5 megabases (Mb) of haploid scanned ASD probands for either de novo CNVs or rare (,1% of exome coding sequence (Supplementary Table 1). In addition, we controls), inherited CNVs. Individual events were validated by either performed copy number analysis on 122 of these families, using a array CGH or genotyping array (see Methods). We identified 76 events combination of the exome data, array comparative genomic hybrid- in 53 individuals, including six de novo (median size 467 kilobases ization (CGH), and genotyping arrays, thereby providing a more com- (kb)) and 70 inherited (median size 155 kb) CNVs (Supplementary prehensive view of rare variation. Table 6). These include disruptions of EHMT1 (Kleefstra’s syndrome, In the 189 new probands, we validated 248 de novo events, 225 single Online Mendelian Inheritance in Man (OMIM) accession 610253), nucleotide variants (SNVs), 17 small insertions/deletions (indels), and CNTNAP4 (reported in children with developmental delay and aut- six CNVs (Supplementary Table 2). These included 181 non- ism11) and the 16p11.2 duplication (OMIM 611913) associated with synonymous changes, of which 120 were classified as severe based developmental delay, bipolar disorder and schizophrenia. on sequence conservation and/or biochemical properties (Methods We performed a multivariate analysis on non-verbal IQ (NVIQ), and Supplementary Table 3). The observed point mutation rate in verbal IQ (VIQ) and the load of ‘extreme’ de novo mutations—where coding sequence was ,1.3 events per trio or 2.17 3 1028 per base extreme is defined as point mutations that truncate proteins, intersect 1Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington 98195, USA. 2Department of Psychiatry and Behavioral Sciences, University of Washington, Seattle, Washington 98195, USA. 3Department of Computer Science and Engineering, University of Washington, Seattle, Washington 98195, USA. 4Santa Fe Institute, Santa Fe, New Mexico 87501, USA. 5Howard Hughes Medical Institute, Seattle, Washington 98195, USA. 246 | NATURE | VOL 485 | 10 MAY 2012 ©2012 Macmillan Publishers Limited. All rights reserved LETTER RESEARCH a b c 140 120 A A T T 80 100 T T 41 10 Non-verbal IQ A G paternal maternal 40 60 Paternal age (months) A G events events 250 350 450 550 T 012 3+ 20 012+ Number of de novo coding mutations Number of extreme de novo mutations d Chr18: 40000000 40500000 41000000 41500000 Cases Controls 18q12.3 SETBP1 SLC14A2 EPG5 SLC14A1 SIGLEC15 Figure 1 | De novo mutation events in autism spectrum disorder. mutation events (0, n 5 138; 1, n 5 41; 21, n 5 10), both with and without a, Haplotype phasing using informative markers shows a strong parent-of- CNVs (Supplementary Discussion). d, Browser images showing CNVs origin bias with 41 of 51 de novo events occurring on the paternally inherited identified in the del(18)(q12.2q21.1) syndrome region. The truncating point haplotype. Arrows represent sequence reads from paternal (blue) or maternal mutation in SETBP1 occurs within the critical region, identifying the likely (red) haplotypes. b, c, Box and whisker plots for 189 SSC probands. b,The causative locus. Each red (deletion) and green (duplication) line represents an paternal estimated age at conception versus the number of observed de novo identified CNV in cases (solid lines) versus controls (dashed lines), with point mutations (0, n 5 53; 1, n 5 65; 2, n 5 44; 31, n 5 27). c, Decreased non- arrowheads showing point mutation. verbal IQ is significantly associated with an increasing number of extreme Mendelian or ASD loci (n 5 57), or de novo CNVs that intersect genes The de novo mutations included truncating events in syndromic (n 5 5) (Fig. 1c and Supplementary Discussion). NVIQ, but not VIQ, intellectual disability genes (MBD5 (mental retardation, autosomal decreased significantly (P , 0.01) with increased number of events. dominant 1, OMIM 156200), RPS6KA3 (Coffin–Lowry syndrome, Covariant analysis of the samples with CNV data showed that this OMIM 303600) and DYRK1A (the Down’s syndrome candidate finding was strengthened, but not exclusively driven, by the presence gene, OMIM 600855)), and missense variants in loci associated with of either de novo or rare CNVs (Supplementary Fig. 5). syndromic ASD, including CHD7, PTEN (macrocephaly/autism Among the de novo events, we identified 62 top ASD risk con- syndrome, OMIM 605309) and TSC2 (tuberous sclerosis complex, tributing mutations based on the deleteriousness of the
Recommended publications
  • Supp Material.Pdf
    Legends for Supplemental Figures and Tables Figure S1. Expression of Tlx during retinogenesis. (A) Staged embryos were stained for β- galactosidase knocked into the Tlx locus to indicate Tlx expression. Tlx was expressed in the neural blast layer in the early phase of neural retina development (blue signal). (B) Expression of Tlx in neural retina was quantified using Q-PCR at multiple developmental stages. Figure S2. Expression of p27kip1 and cyclin D1 (Ccnd1) at various developmental stages in wild-type or Tlx-/- retinas. (A) Q-PCR analysis of p27kip1 mRNA expression. (B) Western blotting analysis of p27kip1 protein expression. (C) Q-PCR analysis of cyclin D1 mRNA expression. Figure S3. Q-PCR analysis of mRNA expression of Sf1 (A), Lrh1 (B), and Atn1 (C) in wild-type mouse retinas. RNAs from testis and liver were used as controls. Table S1. List of genes dysregulated both at E15.5 and P0 Tlx-/- retinas. Gene E15.5 P0 Cluste Gene Title Fold Fold r Name p-value p-value Change Change nuclear receptor subfamily 0, group B, Nr0b1 1.65 0.0024 2.99 0.0035 member 1 1 Pou4f3 1.91 0.0162 2.39 0.0031 POU domain, class 4, transcription factor 3 1 Tcfap2d 2.18 0.0000 2.37 0.0001 transcription factor AP-2, delta 1 Zic5 1.66 0.0002 2.02 0.0218 zinc finger protein of the cerebellum 5 1 Zfpm1 1.85 0.0030 1.88 0.0025 zinc finger protein, multitype 1 1 Pten 1.60 0.0155 1.82 0.0131 phospatase and tensin homolog 2 Itgb5 -1.85 0.0063 -1.85 0.0007 integrin beta 5 2 Gpr49 6.86 0.0001 15.16 0.0001 G protein-coupled receptor 49 3 Cmkor1 2.60 0.0007 2.72 0.0013
    [Show full text]
  • Seq2pathway Vignette
    seq2pathway Vignette Bin Wang, Xinan Holly Yang, Arjun Kinstlick May 19, 2021 Contents 1 Abstract 1 2 Package Installation 2 3 runseq2pathway 2 4 Two main functions 3 4.1 seq2gene . .3 4.1.1 seq2gene flowchart . .3 4.1.2 runseq2gene inputs/parameters . .5 4.1.3 runseq2gene outputs . .8 4.2 gene2pathway . 10 4.2.1 gene2pathway flowchart . 11 4.2.2 gene2pathway test inputs/parameters . 11 4.2.3 gene2pathway test outputs . 12 5 Examples 13 5.1 ChIP-seq data analysis . 13 5.1.1 Map ChIP-seq enriched peaks to genes using runseq2gene .................... 13 5.1.2 Discover enriched GO terms using gene2pathway_test with gene scores . 15 5.1.3 Discover enriched GO terms using Fisher's Exact test without gene scores . 17 5.1.4 Add description for genes . 20 5.2 RNA-seq data analysis . 20 6 R environment session 23 1 Abstract Seq2pathway is a novel computational tool to analyze functional gene-sets (including signaling pathways) using variable next-generation sequencing data[1]. Integral to this tool are the \seq2gene" and \gene2pathway" components in series that infer a quantitative pathway-level profile for each sample. The seq2gene function assigns phenotype-associated significance of genomic regions to gene-level scores, where the significance could be p-values of SNPs or point mutations, protein-binding affinity, or transcriptional expression level. The seq2gene function has the feasibility to assign non-exon regions to a range of neighboring genes besides the nearest one, thus facilitating the study of functional non-coding elements[2]. Then the gene2pathway summarizes gene-level measurements to pathway-level scores, comparing the quantity of significance for gene members within a pathway with those outside a pathway.
    [Show full text]
  • A Clinicopathological and Molecular Genetic Analysis of Low-Grade Glioma in Adults
    A CLINICOPATHOLOGICAL AND MOLECULAR GENETIC ANALYSIS OF LOW-GRADE GLIOMA IN ADULTS Presented by ANUSHREE SINGH MSc A thesis submitted in partial fulfilment of the requirements of the University of Wolverhampton for the degree of Doctor of Philosophy Brain Tumour Research Centre Research Institute in Healthcare Sciences Faculty of Science and Engineering University of Wolverhampton November 2014 i DECLARATION This work or any part thereof has not previously been presented in any form to the University or to any other body whether for the purposes of assessment, publication or for any other purpose (unless otherwise indicated). Save for any express acknowledgments, references and/or bibliographies cited in the work, I confirm that the intellectual content of the work is the result of my own efforts and of no other person. The right of Anushree Singh to be identified as author of this work is asserted in accordance with ss.77 and 78 of the Copyright, Designs and Patents Act 1988. At this date copyright is owned by the author. Signature: Anushree Date: 30th November 2014 ii ABSTRACT The aim of the study was to identify molecular markers that can determine progression of low grade glioma. This was done using various approaches such as IDH1 and IDH2 mutation analysis, MGMT methylation analysis, copy number analysis using array comparative genomic hybridisation and identification of differentially expressed miRNAs using miRNA microarray analysis. IDH1 mutation was present at a frequency of 71% in low grade glioma and was identified as an independent marker for improved OS in a multivariate analysis, which confirms the previous findings in low grade glioma studies.
    [Show full text]
  • Prox1regulates the Subtype-Specific Development of Caudal Ganglionic
    The Journal of Neuroscience, September 16, 2015 • 35(37):12869–12889 • 12869 Development/Plasticity/Repair Prox1 Regulates the Subtype-Specific Development of Caudal Ganglionic Eminence-Derived GABAergic Cortical Interneurons X Goichi Miyoshi,1 Allison Young,1 Timothy Petros,1 Theofanis Karayannis,1 Melissa McKenzie Chang,1 Alfonso Lavado,2 Tomohiko Iwano,3 Miho Nakajima,4 Hiroki Taniguchi,5 Z. Josh Huang,5 XNathaniel Heintz,4 Guillermo Oliver,2 Fumio Matsuzaki,3 Robert P. Machold,1 and Gord Fishell1 1Department of Neuroscience and Physiology, NYU Neuroscience Institute, Smilow Research Center, New York University School of Medicine, New York, New York 10016, 2Department of Genetics & Tumor Cell Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee 38105, 3Laboratory for Cell Asymmetry, RIKEN Center for Developmental Biology, Kobe 650-0047, Japan, 4Laboratory of Molecular Biology, Howard Hughes Medical Institute, GENSAT Project, The Rockefeller University, New York, New York 10065, and 5Cold Spring Harbor Laboratory, Cold Spring Harbor, New York 11724 Neurogliaform (RELNϩ) and bipolar (VIPϩ) GABAergic interneurons of the mammalian cerebral cortex provide critical inhibition locally within the superficial layers. While these subtypes are known to originate from the embryonic caudal ganglionic eminence (CGE), the specific genetic programs that direct their positioning, maturation, and integration into the cortical network have not been eluci- dated. Here, we report that in mice expression of the transcription factor Prox1 is selectively maintained in postmitotic CGE-derived cortical interneuron precursors and that loss of Prox1 impairs the integration of these cells into superficial layers. Moreover, Prox1 differentially regulates the postnatal maturation of each specific subtype originating from the CGE (RELN, Calb2/VIP, and VIP).
    [Show full text]
  • Systematic Detection of Divergent Brain Proteins in Human Evolution and Their Roles in Cognition
    bioRxiv preprint doi: https://doi.org/10.1101/658658; this version posted June 3, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Systematic detection of divergent brain proteins in human evolution and their roles in cognition Guillaume Dumas1,*, Simon Malesys1 and Thomas Bourgeron1 Affiliations: 1 Human Genetics and Cognitive Functions, Institut Pasteur, UMR3571 CNRS, Université de Paris, Paris, (75015) France * Corresponding author: [email protected] 1 bioRxiv preprint doi: https://doi.org/10.1101/658658; this version posted June 3, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Abstract The human brain differs from that of other primates, but the underlying genetic mechanisms remain unclear. Here we measured the evolutionary pressures acting on all human protein- coding genes (N=17,808) based on their divergence from early hominins such as Neanderthal, and non-human primates. We confirm that genes encoding brain-related proteins are among the most conserved of the human proteome. Conversely, several of the most divergent proteins in humans compared to other primates are associated with brain-associated diseases such as micro/macrocephaly, dyslexia, and autism. We identified specific eXpression profiles of a set of divergent genes in ciliated cells of the cerebellum, that might have contributed to the emergence of fine motor skills and social cognition in humans.
    [Show full text]
  • Supplementary Table 1: Adhesion Genes Data Set
    Supplementary Table 1: Adhesion genes data set PROBE Entrez Gene ID Celera Gene ID Gene_Symbol Gene_Name 160832 1 hCG201364.3 A1BG alpha-1-B glycoprotein 223658 1 hCG201364.3 A1BG alpha-1-B glycoprotein 212988 102 hCG40040.3 ADAM10 ADAM metallopeptidase domain 10 133411 4185 hCG28232.2 ADAM11 ADAM metallopeptidase domain 11 110695 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 195222 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 165344 8751 hCG20021.3 ADAM15 ADAM metallopeptidase domain 15 (metargidin) 189065 6868 null ADAM17 ADAM metallopeptidase domain 17 (tumor necrosis factor, alpha, converting enzyme) 108119 8728 hCG15398.4 ADAM19 ADAM metallopeptidase domain 19 (meltrin beta) 117763 8748 hCG20675.3 ADAM20 ADAM metallopeptidase domain 20 126448 8747 hCG1785634.2 ADAM21 ADAM metallopeptidase domain 21 208981 8747 hCG1785634.2|hCG2042897 ADAM21 ADAM metallopeptidase domain 21 180903 53616 hCG17212.4 ADAM22 ADAM metallopeptidase domain 22 177272 8745 hCG1811623.1 ADAM23 ADAM metallopeptidase domain 23 102384 10863 hCG1818505.1 ADAM28 ADAM metallopeptidase domain 28 119968 11086 hCG1786734.2 ADAM29 ADAM metallopeptidase domain 29 205542 11085 hCG1997196.1 ADAM30 ADAM metallopeptidase domain 30 148417 80332 hCG39255.4 ADAM33 ADAM metallopeptidase domain 33 140492 8756 hCG1789002.2 ADAM7 ADAM metallopeptidase domain 7 122603 101 hCG1816947.1 ADAM8 ADAM metallopeptidase domain 8 183965 8754 hCG1996391 ADAM9 ADAM metallopeptidase domain 9 (meltrin gamma) 129974 27299 hCG15447.3 ADAMDEC1 ADAM-like,
    [Show full text]
  • Transcriptome Sequencing and Genome-Wide Association Analyses Reveal Lysosomal Function and Actin Cytoskeleton Remodeling in Schizophrenia and Bipolar Disorder
    Molecular Psychiatry (2015) 20, 563–572 © 2015 Macmillan Publishers Limited All rights reserved 1359-4184/15 www.nature.com/mp ORIGINAL ARTICLE Transcriptome sequencing and genome-wide association analyses reveal lysosomal function and actin cytoskeleton remodeling in schizophrenia and bipolar disorder Z Zhao1,6,JXu2,6, J Chen3,6, S Kim4, M Reimers3, S-A Bacanu3,HYu1, C Liu5, J Sun1, Q Wang1, P Jia1,FXu2, Y Zhang2, KS Kendler3, Z Peng2 and X Chen3 Schizophrenia (SCZ) and bipolar disorder (BPD) are severe mental disorders with high heritability. Clinicians have long noticed the similarities of clinic symptoms between these disorders. In recent years, accumulating evidence indicates some shared genetic liabilities. However, what is shared remains elusive. In this study, we conducted whole transcriptome analysis of post-mortem brain tissues (cingulate cortex) from SCZ, BPD and control subjects, and identified differentially expressed genes in these disorders. We found 105 and 153 genes differentially expressed in SCZ and BPD, respectively. By comparing the t-test scores, we found that many of the genes differentially expressed in SCZ and BPD are concordant in their expression level (q ⩽ 0.01, 53 genes; q ⩽ 0.05, 213 genes; q ⩽ 0.1, 885 genes). Using genome-wide association data from the Psychiatric Genomics Consortium, we found that these differentially and concordantly expressed genes were enriched in association signals for both SCZ (Po10 − 7) and BPD (P = 0.029). To our knowledge, this is the first time that a substantially large number of genes show concordant expression and association for both SCZ and BPD. Pathway analyses of these genes indicated that they are involved in the lysosome, Fc gamma receptor-mediated phagocytosis, regulation of actin cytoskeleton pathways, along with several cancer pathways.
    [Show full text]
  • Supplemental Information
    Supplemental information Dissection of the genomic structure of the miR-183/96/182 gene. Previously, we showed that the miR-183/96/182 cluster is an intergenic miRNA cluster, located in a ~60-kb interval between the genes encoding nuclear respiratory factor-1 (Nrf1) and ubiquitin-conjugating enzyme E2H (Ube2h) on mouse chr6qA3.3 (1). To start to uncover the genomic structure of the miR- 183/96/182 gene, we first studied genomic features around miR-183/96/182 in the UCSC genome browser (http://genome.UCSC.edu/), and identified two CpG islands 3.4-6.5 kb 5’ of pre-miR-183, the most 5’ miRNA of the cluster (Fig. 1A; Fig. S1 and Seq. S1). A cDNA clone, AK044220, located at 3.2-4.6 kb 5’ to pre-miR-183, encompasses the second CpG island (Fig. 1A; Fig. S1). We hypothesized that this cDNA clone was derived from 5’ exon(s) of the primary transcript of the miR-183/96/182 gene, as CpG islands are often associated with promoters (2). Supporting this hypothesis, multiple expressed sequences detected by gene-trap clones, including clone D016D06 (3, 4), were co-localized with the cDNA clone AK044220 (Fig. 1A; Fig. S1). Clone D016D06, deposited by the German GeneTrap Consortium (GGTC) (http://tikus.gsf.de) (3, 4), was derived from insertion of a retroviral construct, rFlpROSAβgeo in 129S2 ES cells (Fig. 1A and C). The rFlpROSAβgeo construct carries a promoterless reporter gene, the β−geo cassette - an in-frame fusion of the β-galactosidase and neomycin resistance (Neor) gene (5), with a splicing acceptor (SA) immediately upstream, and a polyA signal downstream of the β−geo cassette (Fig.
    [Show full text]
  • NICU Gene List Generator.Xlsx
    Neonatal Crisis Sequencing Panel Gene List Genes: A2ML1 - B3GLCT A2ML1 ADAMTS9 ALG1 ARHGEF15 AAAS ADAMTSL2 ALG11 ARHGEF9 AARS1 ADAR ALG12 ARID1A AARS2 ADARB1 ALG13 ARID1B ABAT ADCY6 ALG14 ARID2 ABCA12 ADD3 ALG2 ARL13B ABCA3 ADGRG1 ALG3 ARL6 ABCA4 ADGRV1 ALG6 ARMC9 ABCB11 ADK ALG8 ARPC1B ABCB4 ADNP ALG9 ARSA ABCC6 ADPRS ALK ARSL ABCC8 ADSL ALMS1 ARX ABCC9 AEBP1 ALOX12B ASAH1 ABCD1 AFF3 ALOXE3 ASCC1 ABCD3 AFF4 ALPK3 ASH1L ABCD4 AFG3L2 ALPL ASL ABHD5 AGA ALS2 ASNS ACAD8 AGK ALX3 ASPA ACAD9 AGL ALX4 ASPM ACADM AGPS AMELX ASS1 ACADS AGRN AMER1 ASXL1 ACADSB AGT AMH ASXL3 ACADVL AGTPBP1 AMHR2 ATAD1 ACAN AGTR1 AMN ATL1 ACAT1 AGXT AMPD2 ATM ACE AHCY AMT ATP1A1 ACO2 AHDC1 ANK1 ATP1A2 ACOX1 AHI1 ANK2 ATP1A3 ACP5 AIFM1 ANKH ATP2A1 ACSF3 AIMP1 ANKLE2 ATP5F1A ACTA1 AIMP2 ANKRD11 ATP5F1D ACTA2 AIRE ANKRD26 ATP5F1E ACTB AKAP9 ANTXR2 ATP6V0A2 ACTC1 AKR1D1 AP1S2 ATP6V1B1 ACTG1 AKT2 AP2S1 ATP7A ACTG2 AKT3 AP3B1 ATP8A2 ACTL6B ALAS2 AP3B2 ATP8B1 ACTN1 ALB AP4B1 ATPAF2 ACTN2 ALDH18A1 AP4M1 ATR ACTN4 ALDH1A3 AP4S1 ATRX ACVR1 ALDH3A2 APC AUH ACVRL1 ALDH4A1 APTX AVPR2 ACY1 ALDH5A1 AR B3GALNT2 ADA ALDH6A1 ARFGEF2 B3GALT6 ADAMTS13 ALDH7A1 ARG1 B3GAT3 ADAMTS2 ALDOB ARHGAP31 B3GLCT Updated: 03/15/2021; v.3.6 1 Neonatal Crisis Sequencing Panel Gene List Genes: B4GALT1 - COL11A2 B4GALT1 C1QBP CD3G CHKB B4GALT7 C3 CD40LG CHMP1A B4GAT1 CA2 CD59 CHRNA1 B9D1 CA5A CD70 CHRNB1 B9D2 CACNA1A CD96 CHRND BAAT CACNA1C CDAN1 CHRNE BBIP1 CACNA1D CDC42 CHRNG BBS1 CACNA1E CDH1 CHST14 BBS10 CACNA1F CDH2 CHST3 BBS12 CACNA1G CDK10 CHUK BBS2 CACNA2D2 CDK13 CILK1 BBS4 CACNB2 CDK5RAP2
    [Show full text]
  • Identification of Gene Modules Associated with Warfarin Dosage by a Genome-Wide DNA Methylation Study
    ORIGINAL ARTICLES Department of Clinical Pharmacology1, Xiangya Hospital, Institute of Clinical Pharmacology2, Central South University, Hunan Key Laboratory of Pharmacogenetics; Changsha; Department of Cardio-Thoracic Surgery3, the Second Xiangya Hospital of Central South University, Changsha; Key Laboratory of Bio-resources and Eco-environment4 (Ministry of Education), College of Life Science, Sichuan University, Chengdu, China Identification of gene modules associated with warfarin dosage by a genome-wide DNA methylation study ZHIYING LUO1, 2, RONG LIU1, 2, BAO SUN1, 2, XINMING ZHOU3, ZHAOQIAN LIU1, 2, HONGHAO ZHOU1, 2, HENG XU4, XI LI1,2,*,#, WEI ZHANG1,2,*,# Received December 12, 2017, accepted January 10, 2018 *Corresponding authors: Prof. Wei Zhang, Xi Li, Department of Clinical Pharmacology, Xiangya Hospital, Central South University, 110 Xiangya Rode, Kaifu district, Changsha, Hunan 410008, China [email protected]; [email protected] #Wei Zhang and Xi Li contributed equally to this study. Pharmazie 73: 288–293 (2018) doi: 10.1691/ph.2018.7319 Objective: To identify warfarin dose-associated DNA methylation changes, we conducted the first genomewide DNA methylation association study. Method: A total of 22 patients who required an extreme warfarin dosage from VKORC1 -1639AA & CYP2C9*1*1 genotype group were enrolled in this study. The Illumina Infinium Human- Methylation450 platform was used to perform genome-scale DNA methylation profiling, identifying differentially methylated CpG sites by a nonparametric test. WGCNA was used to analyze the association between gene modules and extreme warfarin dosage. Results: For a total of 378,313 CpG sites that passed the quality control processes, we identified eight differentially methylated CpG probes (p<0.05) showing altered DNA methylation level (>20%) between two extreme dose groups.
    [Show full text]
  • Effects of Chronic Stress on Prefrontal Cortex Transcriptome in Mice Displaying Different Genetic Backgrounds
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Springer - Publisher Connector J Mol Neurosci (2013) 50:33–57 DOI 10.1007/s12031-012-9850-1 Effects of Chronic Stress on Prefrontal Cortex Transcriptome in Mice Displaying Different Genetic Backgrounds Pawel Lisowski & Marek Wieczorek & Joanna Goscik & Grzegorz R. Juszczak & Adrian M. Stankiewicz & Lech Zwierzchowski & Artur H. Swiergiel Received: 14 May 2012 /Accepted: 25 June 2012 /Published online: 27 July 2012 # The Author(s) 2012. This article is published with open access at Springerlink.com Abstract There is increasing evidence that depression signaling pathway (Clic6, Drd1a,andPpp1r1b). LA derives from the impact of environmental pressure on transcriptome affected by CMS was associated with genetically susceptible individuals. We analyzed the genes involved in behavioral response to stimulus effects of chronic mild stress (CMS) on prefrontal cor- (Fcer1g, Rasd2, S100a8, S100a9, Crhr1, Grm5,and tex transcriptome of two strains of mice bred for high Prkcc), immune effector processes (Fcer1g, Mpo,and (HA)and low (LA) swim stress-induced analgesia that Igh-VJ558), diacylglycerol binding (Rasgrp1, Dgke, differ in basal transcriptomic profiles and depression- Dgkg,andPrkcc), and long-term depression (Crhr1, like behaviors. We found that CMS affected 96 and 92 Grm5,andPrkcc) and/or coding elements of dendrites genes in HA and LA mice, respectively. Among genes (Crmp1, Cntnap4,andPrkcc) and myelin proteins with the same expression pattern in both strains after (Gpm6a, Mal,andMog). The results indicate significant CMS, we observed robust upregulation of Ttr gene contribution of genetic background to differences in coding transthyretin involved in amyloidosis, seizures, stress response gene expression in the mouse prefrontal stroke-like episodes, or dementia.
    [Show full text]
  • A Master Autoantigen-Ome Links Alternative Splicing, Female Predilection, and COVID-19 to Autoimmune Diseases
    bioRxiv preprint doi: https://doi.org/10.1101/2021.07.30.454526; this version posted August 4, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. A Master Autoantigen-ome Links Alternative Splicing, Female Predilection, and COVID-19 to Autoimmune Diseases Julia Y. Wang1*, Michael W. Roehrl1, Victor B. Roehrl1, and Michael H. Roehrl2* 1 Curandis, New York, USA 2 Department of Pathology, Memorial Sloan Kettering Cancer Center, New York, USA * Correspondence: [email protected] or [email protected] 1 bioRxiv preprint doi: https://doi.org/10.1101/2021.07.30.454526; this version posted August 4, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Abstract Chronic and debilitating autoimmune sequelae pose a grave concern for the post-COVID-19 pandemic era. Based on our discovery that the glycosaminoglycan dermatan sulfate (DS) displays peculiar affinity to apoptotic cells and autoantigens (autoAgs) and that DS-autoAg complexes cooperatively stimulate autoreactive B1 cell responses, we compiled a database of 751 candidate autoAgs from six human cell types. At least 657 of these have been found to be affected by SARS-CoV-2 infection based on currently available multi-omic COVID data, and at least 400 are confirmed targets of autoantibodies in a wide array of autoimmune diseases and cancer.
    [Show full text]