Integrated Bioinformatics Analysis of Aberrantly-Methylated

Total Page:16

File Type:pdf, Size:1020Kb

Integrated Bioinformatics Analysis of Aberrantly-Methylated Shen et al. BMC Ophthalmology (2020) 20:119 https://doi.org/10.1186/s12886-020-01392-2 RESEARCH ARTICLE Open Access Integrated bioinformatics analysis of aberrantly-methylated differentially- expressed genes and pathways in age- related macular degeneration Yinchen Shen1,2†,MoLi3†, Kun Liu1,2, Xiaoyin Xu1,2, Shaopin Zhu1,2, Ning Wang1,2, Wenke Guo4, Qianqian Zhao4, Ping Lu4, Fudong Yu4 and Xun Xu1,2* Abstract Background: Age-related macular degeneration (AMD) represents the leading cause of visual impairment in the aging population. The goal of this study was to identify aberrantly-methylated, differentially-expressed genes (MDEGs) in AMD and explore the involved pathways via integrated bioinformatics analysis. Methods: Data from expression profile GSE29801 and methylation profile GSE102952 were obtained from the Gene Expression Omnibus database. We analyzed differentially-methylated genes and differentially-expressed genes using R software. Functional enrichment and protein–protein interaction (PPI) network analysis were performed using the R package and Search Tool for the Retrieval of Interacting Genes online database. Hub genes were identified using Cytoscape. Results: In total, 827 and 592 genes showed high and low expression, respectively, in GSE29801; 4117 hyper-methylated genes and 511 hypo-methylated genes were detected in GSE102952. Based on overlap, we categorized 153 genes as hyper-methylated, low-expression genes (Hyper-LGs) and 24 genes as hypo-methylated, high-expression genes (Hypo-HGs). Four Hyper-LGs (CKB, PPP3CA, TGFB2, SOCS2) overlapped with AMD risk genes in the Public Health Genomics and Precision Health Knowledge Base. KEGG pathway enrichment analysis indicated that Hypo-HGs were enriched in the calcium signaling pathway, whereas Hyper-LGs were enriched in sphingolipid metabolism. In GO analysis, Hypo-HGs were enriched in fibroblast migration, membrane raft, and coenzyme binding, among others. Hyper-LGs were enriched in mRNA transport, nuclear speck, and DNA binding, among others. In PPI network analysis, 23 nodes and two edges were established from Hypo-HGs, and 151 nodes and 73 edges were established from Hyper-LGs. Hub genes (DHX9, MAPT, PAX6) showed the greatest overlap. (Continued on next page) * Correspondence: [email protected] †Yinchen Shen and Mo Li are co-first authors. 1Department of Ophthalmology, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai 200080, People’s Republic of China 2National Clinical Research Center for Eye Diseases, Shanghai Key Laboratory of Ocular Fundus Diseases, Shanghai Engineering Center for Visual Science and Photomedicine, Shanghai Engineering Center for Precise Diagnosis and Treatment of Eye Diseases, Shanghai, China Full list of author information is available at the end of the article © The Author(s). 2020 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Shen et al. BMC Ophthalmology (2020) 20:119 Page 2 of 12 (Continued from previous page) Conclusion: This study revealed potentially aberrantly MDEGs and pathways in AMD, which might improve the understanding of this disease. Keywords: Methylation, Gene expression, Age-related macular degeneration, Choroidal neovascularization, Bioinformatics analysis Background gene deletions, gene mutations, and gene insertions, Age-related macular degeneration (AMD) is the leading among other changes. Microarrays based on high- cause of adult blindness in developed countries, particu- throughput platforms are useful and efficient tools to larly in those over 55 years of age. Worldwide, this con- search for meaningful genes and epigenetic alterations dition accounts for 7% of cases of blindness [1, 2] and is for the identification of diagnostic or prognostic bio- expected to affect 288 million people by the year 2040 markers [22]. To better explore the molecular mechan- [3]. AMD first appears as drusen (dry AMD) and ism underlying AMD, it is necessary to conduct full advances to late AMD (wet AMD) characterized by transcriptome analysis at the tissue level; however, as choroidal neovascularization (CNV) [4]. Drusen is com- mentioned previously herein, obtaining ocular tissue posed of lipoproteinaceous deposits and acellular debris samples is difficult. Compared to retina or choroid tis- [5]. CNV involves the growth of new abnormal blood sue, blood samples are relatively easy to obtain from pa- vessels originating from the choroid through a break in tients with AMD. Therefore, it would be helpful to the Bruch’s membrane, which then invade the subretinal evaluate gene expression in the ocular tissue of patients pigment epithelium or sub-retinal space, resulting in se- with AMD as biomarkers in the blood. In the present vere vision loss [6, 7]. AMD affects central fine vision, study, data from gene expression profiling microarrays significantly impairs a patient’s ability to drive, read, and of human retinal and choroidal samples from the Iowa recognize faces, and greatly affects quality of life [8]. As and Oregon cohort AMD and control donors (GSE29801: for wet-AMD (i.e. advanced stages), anti-vascular endo- https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE2 thelial growth factor (anti-VEGF) therapy was shown to 9801)[23], as well as gene methylation profiling microar- be effective and has become the first choice for the treat- rays of the peripheral blood of subjects with AMD ment of CNV [9]. However, anti-VEGF therapy requires (GSE102952: https://www.ncbi.nlm.nih.gov/geo/query/acc. repeated intra-vitreal injections, which are associated cgi?acc=GSE102952), were integrated and analyzed using a with a risk of infection and treatment burden for both series of bioinformatics tools. More precise screening re- the patients and the ophthalmologists [10]. Moreover, sults were obtained by overlapping these two AMD data some patients have poor response to the drugs after a sets. Few studies have attempted to combine gene expres- long-term treatment [11]. Apart from anti-VEGF drugs, sion profile microarrays and gene methylation profile some other therapies have also impacted macular disease microarrays to understand the development of AMD. The treatment and showed their effectiveness, such as dexa- introduction of DNA methylation characteristics in the methasone implant [12–14]. blood is useful to understand the characteristics of AMD Previous studies indicated that many environmental disease at the tissue level. factors are associated with an increased risk of AMD, In the present study, data from gene expression profiling such as age, race, smoking, obesity, and hypertension microarrays and gene methylation profiling microarrays [15, 16]. Additionally, genetic factors are regarded as were integrated and analyzed. Functional enrichment and important for the initiation and progression of AMD protein–protein interaction (PPI) network analyses of [17, 18]. Comparing to tumor tissue samples, ocular fun- screened genes were performed using the R package “clus- dus tissue samples (i.e. retina and choroid) of patients terProfiler” and Search Tool for the Retrieval of Interact- with AMD are quite difficult to obtain in real world. ing Genes (STRING) online database. We identified The difficulties to obtain human fundus tissues re- methylated genes in the peripheral blood, which might be stricted our understanding of this blinding disease. Over useful as biomarkers for the precise diagnosis and treat- the past few years, most genetic studies on AMD were ment of AMD. case-control genome-wide association studies (GWASs) of single-nucleotide polymorphisms in the patients’ per- ipheral blood [19–21], which were valuable but provided Methods limited information. The need for ethics approval was waived by the ethics With the rapid development of gene assay technology, committee of Shanghai General Hospital, Shanghai Jiao studies of disease pathogenesis are no longer limited to Tong University School of Medicine, Shanghai, China (the document is attached as an additional file). Shen et al. BMC Ophthalmology (2020) 20:119 Page 3 of 12 Microarray data information (GSE102952: https://www.ncbi.nlm.nih.gov/geo/query/ We identified methylated, differentially-expressed genes acc.cgi?acc=GSE102952). The methylation profile data of (MDEGs) between AMD and control samples by analyz- all nine patients with AMD and nine normal controls ing mRNA microarray and methylation profiling data- were included in this analysis. sets. In this study, a gene expression profiling dataset (GSE29801: https://www.ncbi.nlm.nih.gov/geo/query/
Recommended publications
  • A Genomic Screen for Activators of the Antioxidant Response Element
    A genomic screen for activators of the antioxidant response element Yanxia Liu†, Jonathan T. Kern‡, John R. Walker§, Jeffrey A. Johnson‡, Peter G. Schultz§¶ʈ, and Hendrik Luesch†ʈ †Department of Medicinal Chemistry, University of Florida, 1600 Southwest Archer Road, Gainesville, FL 32610; ‡Department of Pharmaceutical Sciences, School of Pharmacy, University of Wisconsin, 777 Highland Avenue, Madison, WI 53705; §Genomics Institute of the Novartis Research Foundation, 10675 John Jay Hopkins Drive, San Diego, CA 92121; and ¶Department of Chemistry, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037 Contributed by Peter G. Schultz, January 31, 2007 (sent for review December 13, 2006) The antioxidant response element (ARE) is a cis-acting regulatory activation, translocates into the nucleus and transcriptionally acti- enhancer element found in the 5؅ flanking region of many phase II vates ARE-dependent genes after recruiting Maf proteins (2). The detoxification enzymes. Up-regulation of ARE-dependent target upstream regulatory mechanisms by which ARE-activating signals genes is known to have neuroprotective effects; yet, the mechanism are linked to Nrf2 remain to be fully elucidated. It has been of activation is largely unknown. By screening an arrayed collection demonstrated that reactive sulfhydryl groups of Keap1 are sensors of Ϸ15,000 full-length expression cDNAs in the human neuroblas- for induction of phase II genes (13), leading to the proposal that the toma cell line IMR-32 with an ARE-luciferase reporter, we have Nrf2/Keap1 interaction represents a cytoplasmic sensor for oxida- identified several cDNAs not previously associated with ARE activa- tive stress. However, 1-phosphatidylinositol 3-kinase (PI3K), tion.
    [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]
  • Noelia Díaz Blanco
    Effects of environmental factors on the gonadal transcriptome of European sea bass (Dicentrarchus labrax), juvenile growth and sex ratios Noelia Díaz Blanco Ph.D. thesis 2014 Submitted in partial fulfillment of the requirements for the Ph.D. degree from the Universitat Pompeu Fabra (UPF). This work has been carried out at the Group of Biology of Reproduction (GBR), at the Department of Renewable Marine Resources of the Institute of Marine Sciences (ICM-CSIC). Thesis supervisor: Dr. Francesc Piferrer Professor d’Investigació Institut de Ciències del Mar (ICM-CSIC) i ii A mis padres A Xavi iii iv Acknowledgements This thesis has been made possible by the support of many people who in one way or another, many times unknowingly, gave me the strength to overcome this "long and winding road". First of all, I would like to thank my supervisor, Dr. Francesc Piferrer, for his patience, guidance and wise advice throughout all this Ph.D. experience. But above all, for the trust he placed on me almost seven years ago when he offered me the opportunity to be part of his team. Thanks also for teaching me how to question always everything, for sharing with me your enthusiasm for science and for giving me the opportunity of learning from you by participating in many projects, collaborations and scientific meetings. I am also thankful to my colleagues (former and present Group of Biology of Reproduction members) for your support and encouragement throughout this journey. To the “exGBRs”, thanks for helping me with my first steps into this world. Working as an undergrad with you Dr.
    [Show full text]
  • Molecular Evolutionary Analysis of Plastid Genomes in Nonphotosynthetic Angiosperms and Cancer Cell Lines
    The Pennsylvania State University The Graduate School Department or Biology MOLECULAR EVOLUTIONARY ANALYSIS OF PLASTID GENOMES IN NONPHOTOSYNTHETIC ANGIOSPERMS AND CANCER CELL LINES A Dissertation in Biology by Yan Zhang 2012 Yan Zhang Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy Dec 2012 The Dissertation of Yan Zhang was reviewed and approved* by the following: Schaeffer, Stephen W. Professor of Biology Chair of Committee Ma, Hong Professor of Biology Altman, Naomi Professor of Statistics dePamphilis, Claude W Professor of Biology Dissertation Adviser Douglas Cavener Professor of Biology Head of Department of Biology *Signatures are on file in the Graduate School iii ABSTRACT This thesis explores the application of evolutionary theory and methods in understanding the plastid genome of nonphotosynthetic parasitic plants and role of mutations in tumor proliferations. We explore plastid genome evolution in parasitic angiosperms lineages that have given up the primary function of plastid genome – photosynthesis. Genome structure, gene contents, and evolutionary dynamics were analyzed and compared in both independent and related parasitic plant lineages. Our studies revealed striking similarities in changes of gene content and evolutionary dynamics with the loss of photosynthetic ability in independent nonphotosynthetic plant lineages. Evolutionary analysis suggests accelerated evolution in the plastid genome of the nonphotosynthetic plants. This thesis also explores the application of phylogenetic and evolutionary analysis in cancer biology. Although cancer has often been likened to Darwinian process, very little application of molecular evolutionary analysis has been seen in cancer biology research. In our study, phylogenetic approaches were used to explore the relationship of several hundred established cancer cell lines based on multiple sequence alignments constructed with variant codons and residues across 494 and 523 genes.
    [Show full text]
  • Proteomic Analysis of Ubiquitin Ligase KEAP1 Reveals Associated Proteins That Inhibit NRF2 Ubiquitination
    Published OnlineFirst February 4, 2013; DOI: 10.1158/0008-5472.CAN-12-4400 Cancer Molecular and Cellular Pathobiology Research Proteomic Analysis of Ubiquitin Ligase KEAP1 Reveals Associated Proteins That Inhibit NRF2 Ubiquitination Bridgid E. Hast1, Dennis Goldfarb2, Kathleen M. Mulvaney1, Michael A. Hast4, Priscila F. Siesser1, Feng Yan1, D. Neil Hayes3, and Michael B. Major1,2 Abstract Somatic mutations in the KEAP1 ubiquitin ligase or its substrate NRF2 (NFE2L2) commonly occur in human cancer, resulting in constitutive NRF2-mediated transcription of cytoprotective genes. However, many tumors display high NRF2 activity in the absence of mutation, supporting the hypothesis that alternative mechanisms of pathway activation exist. Previously, we and others discovered that via a competitive binding mechanism, the proteins WTX (AMER1), PALB2, and SQSTM1 bind KEAP1 to activate NRF2. Proteomic analysis of the KEAP1 protein interaction network revealed a significant enrichment of associated proteins containing an ETGE amino acid motif, which matches the KEAP1 interaction motif found in NRF2. Like WTX, PALB2, and SQSTM1, we found that the dipeptidyl peptidase 3 (DPP3) protein binds KEAP1 via an "ETGE" motif to displace NRF2, thus inhibiting NRF2 ubiquitination and driving NRF2-dependent transcription. Comparing the spectrum of KEAP1-interacting proteins with the genomic profile of 178 squamous cell lung carcinomas characterized by The Cancer Genome Atlas revealed amplification and mRNA overexpression of the DPP3 gene in tumors with high NRF2 activity but lacking NRF2 stabilizing mutations. We further show that tumor-derived mutations in KEAP1 are hypomorphic with respect to NRF2 inhibition and that DPP3 overexpression in the presence of these mutants further promotes NRF2 activation.
    [Show full text]
  • A Peripheral Blood Gene Expression Signature to Diagnose Subclinical Acute Rejection
    CLINICAL RESEARCH www.jasn.org A Peripheral Blood Gene Expression Signature to Diagnose Subclinical Acute Rejection Weijia Zhang,1 Zhengzi Yi,1 Karen L. Keung,2 Huimin Shang,3 Chengguo Wei,1 Paolo Cravedi,1 Zeguo Sun,1 Caixia Xi,1 Christopher Woytovich,1 Samira Farouk,1 Weiqing Huang,1 Khadija Banu,1 Lorenzo Gallon,4 Ciara N. Magee,5 Nader Najafian,5 Milagros Samaniego,6 Arjang Djamali ,7 Stephen I. Alexander,2 Ivy A. Rosales,8 Rex Neal Smith,8 Jenny Xiang,3 Evelyne Lerut,9 Dirk Kuypers,10,11 Maarten Naesens ,10,11 Philip J. O’Connell,2 Robert Colvin,8 Madhav C. Menon,1 and Barbara Murphy1 Due to the number of contributing authors, the affiliations are listed at the end of this article. ABSTRACT Background In kidney transplant recipients, surveillance biopsies can reveal, despite stable graft function, histologic features of acute rejection and borderline changes that are associated with undesirable graft outcomes. Noninvasive biomarkers of subclinical acute rejection are needed to avoid the risks and costs associated with repeated biopsies. Methods We examined subclinical histologic and functional changes in kidney transplant recipients from the prospective Genomics of Chronic Allograft Rejection (GoCAR) study who underwent surveillance biopsies over 2 years, identifying those with subclinical or borderline acute cellular rejection (ACR) at 3 months (ACR-3) post-transplant. We performed RNA sequencing on whole blood collected from 88 indi- viduals at the time of 3-month surveillance biopsy to identify transcripts associated with ACR-3, developed a novel sequencing-based targeted expression assay, and validated this gene signature in an independent cohort.
    [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]
  • Population-Haplotype Models for Mapping and Tagging Structural Variation Using Whole Genome Sequencing
    Population-haplotype models for mapping and tagging structural variation using whole genome sequencing Eleni Loizidou Submitted in part fulfilment of the requirements for the degree of Doctor of Philosophy Section of Genomics of Common Disease Department of Medicine Imperial College London, 2018 1 Declaration of originality I hereby declare that the thesis submitted for a Doctor of Philosophy degree is based on my own work. Proper referencing is given to the organisations/cohorts I collaborated with during the project. 2 Copyright Declaration The copyright of this thesis rests with the author and is made available under a Creative Commons Attribution Non-Commercial No Derivatives licence. Researchers are free to copy, distribute or transmit the thesis on the condition that they attribute it, that they do not use it for commercial purposes and that they do not alter, transform or build upon it. For any reuse or redistribution, researchers must make clear to others the licence terms of this work 3 Abstract The scientific interest in copy number variation (CNV) is rapidly increasing, mainly due to the evidence of phenotypic effects and its contribution to disease susceptibility. Single nucleotide polymorphisms (SNPs) which are abundant in the human genome have been widely investigated in genome-wide association studies (GWAS). Despite the notable genomic effects both CNVs and SNPs have, the correlation between them has been relatively understudied. In the past decade, next generation sequencing (NGS) has been the leading high-throughput technology for investigating CNVs and offers mapping at a high-quality resolution. We created a map of NGS-defined CNVs tagged by SNPs using the 1000 Genomes Project phase 3 (1000G) sequencing data to examine patterns between the two types of variation in protein-coding genes.
    [Show full text]
  • A Genetic Mapping System in Caenorhabditis Elegans Based On
    z A Genetic Mapping System in Caenorhabditis elegansBased on Polymorphic Sequence-Tagged Sites Benjamin D. Williams,* Bertold Schrank,* Chau Huynh,* Ratna Shownkeent and Robert H. Waterston* *Department of Genetics, Washington University School ofMedicine, St. Louis, Missouri631 IO, and ?Medical Research Council Laboratory of Molecular Biology, Cambridge CB2 2QH, United Kingdom Manuscript received December 26, 1991 Accepted for publication March 9, 1992 ABSTRACT We devised an efficient genetic mapping system in the nematode Caenorhabditis elegans which is based upon the differences in number and location of the transposable element Tcl between the Bristol and Bergerac strains. Using the nearly completed physical map of theC. elegans genome, we selected 40 widely distributed sites which contain a Tcl element in the Bergerac strain, but not in the Bristol strain. For each site a polymerase chain reaction assay was designed that can distinguish between the Bergerac Tcl-containing site and the Bristol “empty” site. By combining appropriate assays in a single reaction, one can score multiple sites within single worms. This permits a mutation to be rapidly mapped,first to a linkage group and thento a chromosomal subregion, through analysis of only a small number of progeny from a single interstrain cross. EQUENCE-tagged sites (STSs), which havere- Previously, Tcls have been used in efforts to mo- S cently been proposed as common landmarks for lecularly clone specific genes, both by direct transpo- genomic physical mapping (OLSONet al. 1989), are son tagging(MOERMAN, BENIAN and WATERSTON short, unique genomic regions that areeasily detected 1986) andby genetic mappingof the Bristol/Bergerac by polymerase chain reaction(PCR) amplification Tcl polymorphisms nearest a gene to identify molec- (SAIKIet al.
    [Show full text]
  • Mouse Fam76b Knockout Project (CRISPR/Cas9)
    https://www.alphaknockout.com Mouse Fam76b Knockout Project (CRISPR/Cas9) Objective: To create a Fam76b knockout Mouse model (C57BL/6J) by CRISPR/Cas-mediated genome engineering. Strategy summary: The Fam76b gene (NCBI Reference Sequence: NM_176836 ; Ensembl: ENSMUSG00000037808 ) is located on Mouse chromosome 9. 10 exons are identified, with the ATG start codon in exon 1 and the TGA stop codon in exon 10 (Transcript: ENSMUST00000059579). Exon 3~8 will be selected as target site. Cas9 and gRNA will be co-injected into fertilized eggs for KO Mouse production. The pups will be genotyped by PCR followed by sequencing analysis. Note: Exon 3 starts from about 15.04% of the coding region. Exon 3~8 covers 66.47% of the coding region. The size of effective KO region: ~7248 bp. The KO region does not have any other known gene. Page 1 of 9 https://www.alphaknockout.com Overview of the Targeting Strategy Wildtype allele 5' gRNA region gRNA region 3' 1 3 4 5 6 7 8 10 Legends Exon of mouse Fam76b Knockout region Page 2 of 9 https://www.alphaknockout.com Overview of the Dot Plot (up) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 900 bp section upstream of Exon 3 is aligned with itself to determine if there are tandem repeats. No significant tandem repeat is found in the dot plot matrix. So this region is suitable for PCR screening or sequencing analysis. Overview of the Dot Plot (down) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 2000 bp section downstream of Exon 8 is aligned with itself to determine if there are tandem repeats.
    [Show full text]
  • Dynamics of Meiotic Sex Chromosome Inactivation And
    bioRxiv preprint doi: https://doi.org/10.1101/665372; this version posted July 8, 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-ND 4.0 International license. 1 Dynamics of Meiotic Sex Chromosome Inactivation 2 and Pachytene Activation in Mice Spermatogenesis 3 4 Ábel Vértesy1,2; Javier Frias-Aldeguer1,4; Zeliha Sahin1,3; Nicolas Rivron1,4; Alexander van 5 Oudenaarden1,2 and Niels Geijsen1,5 6 7 1. Hubrecht Institute-KNAW (Royal Netherlands Academy of Arts and Sciences) and 8 University Medical Center, 3584 CT Utrecht, The Netherlands 9 2. Department of Genetics, Center for Molecular Medicine, Cancer Genomics Netherlands, 10 University Medical Center Utrecht, The Netherlands 11 3. Amsterdam UMC, University of Amsterdam, Clinical Genetics, Amsterdam Medical 12 Research, Meibergdreef 9, Amsterdam, Netherlands 13 4. MERLN Institute for Technology-Inspired Regenerative Medicine, Maastricht University, 14 The Netherlands 15 5. Faculty of Veterinary Medicine, Clinical Sciences of Companion Animals, Utrecht 16 University, The Netherlands 17 18 Abstract 19 During germ cell development, cells undergo a drastic switch from mitosis to meiosis to 20 form haploid germ cells. Sequencing and computational technologies now allow studying 21 development at the single-cell level. Here we developed a multiplexed trajectory 22 reconstruction to create a high-resolution developmental map of spermatogonia and 23 prophase-I spermatocytes from testes of a Dazl-GFP reporter mouse. We identified three 24 main transitions in the meiotic prophase-I: meiotic entry, the meiotic sex chromosome 25 inactivation (MSCI), and concomitant pachytene activation.
    [Show full text]
  • Supplementary Information Contents
    Supplementary Information Contents Supplementary Methods: Additional methods descriptions Supplementary Results: Biology of suicidality-associated loci Supplementary Figures Supplementary Figure 1: Flow chart of UK Biobank participants available for primary analyses (Ordinal GWAS and PRS analysis) Supplementary Figure 2: Flow chart of UK Biobank participants available for secondary analyses. The flow chart of participants is the same as Supplementary Figure 1 up to the highlighted box. Relatedness exclusions were applied for A) the DSH GWAS considering the categories Controls, Contemplated self-harm and Actual self-ham and B) the SIA GWAS considering the categories Controls, Suicidal ideation and attempted suicide. Supplementary Figure 3: Manhattan plot of GWAS of ordinal DSH in UK Biobank (N=100 234). Dashed red line = genome wide significance threshold (p<5x10-5). Inset: QQ plot for genome-wide association with DSH. Red line = theoretical distribution under the null hypothesis of no association. Supplementary Figure 4: Manhattan plot of GWAS of ordinal SIA in UK Biobank (N=108 090). Dashed red line = genome wide significance threshold (p<5x10-5). Inset: QQ plot for genome-wide association with SIA. Red line = theoretical distribution under the null hypothesis of no association. Supplementary Figure 5: Manhattan plot of gene-based GWAS of ordinal suicide in UK Biobank (N=122 935). Dashed red line = genome wide significance threshold (p<5x10-5). Inset: QQ plot for genome-wide association with suicidality in UK Biobank. Red line = theoretical distribution under the null hypothesis of no association. Supplementary Figure 6: Manhattan plot of gene-based GWAS of ordinal DSH in UK Biobank (N=100 234).
    [Show full text]