Hidden Genetic Variation in LCA9-Associated Congenital Blindness Explained by 5’UTR Mutations and Copy-Number Variations of NMNAT1 Frauke Coppieters, Anne Laure Todeschini, Takuro Fujimaki, Annelot Baert, Marieke de Bruyne, Caroline van Cauwenbergh, Hannah Verdin, Miriam Bauwens, Maté Ongenaert, Mineo Kondo, et al.
To cite this version:
Frauke Coppieters, Anne Laure Todeschini, Takuro Fujimaki, Annelot Baert, Marieke de Bruyne, et al.. Hidden Genetic Variation in LCA9-Associated Congenital Blindness Explained by 5’UTR Muta- tions and Copy-Number Variations of NMNAT1. Human Mutation, Wiley, 2015, 36 (12), pp.1188- 1196. 10.1002/humu.22899. hal-02119195
HAL Id: hal-02119195 https://hal.archives-ouvertes.fr/hal-02119195 Submitted on 3 May 2019
HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. RESEARCH ARTICLE
OFFICIAL JOURNAL Hidden Genetic Variation in LCA9-Associated Congenital Blindness Explained by 5UTR Mutations and Copy-Number
Variations of NMNAT1 www.hgvs.org
Frauke Coppieters,1 † Anne Laure Todeschini,2 † Takuro Fujimaki,3 Annelot Baert,1 Marieke De Bruyne,1 Caroline Van Cauwenbergh,1 Hannah Verdin,1 Miriam Bauwens,1 Mate´ Ongenaert,1 Mineo Kondo,4 Franc¸oiseMeire,5 Akira Murakami,3 Reiner A. Veitia,2 Bart P. Leroy,1,6,7 and Elfride De Baere1∗ 1Center for Medical Genetics Ghent, Ghent University, Ghent, Belgium; 2Institut Jacques Monod, UMR 7592 CNRS-Universite´ Paris Diderot, Paris, France; 3Department of Ophthalmology, Juntendo University Graduate School of Medicine, Tokyo, Japan; 4Department of Ophthalmology, Mie University Graduate School of Medicine, Mie, Japan; 5Department of Ophthalmology, Queen Fabiola Children’s University Hospital, Brussels, Belgium; 6Department of Ophthalmology, Ghent University Hospital, Ghent, Belgium; 7Division of Ophthalmology, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania Communicated by Daniel Schorderet Received 15 April 2015; accepted revised manuscript 19 August 2015. Published online 28 August 2015 in Wiley Online Library (www.wiley.com/humanmutation). DOI: 10.1002/humu.22899
Introduction ABSTRACT: Leber congenital amaurosis (LCA) is a se- Leber congenital amaurosis (LCA; MIM #204000) is the earli- vere autosomal-recessive retinal dystrophy leading to congen- est and most severe inherited retinal dystrophy, causing profound ital blindness. A recently identified LCA gene is NMNAT1, visual deficiency, nystagmus, and an undetectable or severely re- located in the LCA9 locus. Although most mutations in blind- ness genes are coding variations, there is accumulating evidence duced electroretinogram (ERG) in the 1st year of life [Leber, 1869; for hidden noncoding defects or structural variations (SVs). Franceschetti and Dieterle, 1954]. The disease may develop as a The starting point of this study was an LCA9-associated con- cone-rod (type I) or rod-cone dystrophy (type II), with associated sanguineous family in which no coding mutations were found features further impairing visual function [Hanein et al., 2004]. So in the LCA9 region. Exploring the untranslated regions of far, 19 disease genes have been associated with autosomal-recessive NMNAT1 revealed a novel homozygous 5UTR variant, c.- LCA (RetNet). The encoded proteins function in retina-specific as 70A>T. Moreover, an adjacent 5UTR variant, c.-69C>T, well as general processes [den Hollander et al., 2008]. In recent was identified in a second consanguineous family displaying a years, genetic causes have been found in up to 70% of cases, fa- similar phenotype. Both 5 UTR variants resulted in decreased cilitated by massively parallel sequencing (MPS) technologies such NMNAT1 mRNA abundance in patients’ lymphocytes, and as targeted resequencing of linkage intervals or whole-exome se- caused decreased luciferase activity in human retinal pigment quencing (WES) [Neveling et al., 2013]. Although most mutations epithelial RPE-1 cells. Second, we unraveled pseudohomozy- underlying LCA are located in coding regions, there is accumulat- gosity of a coding NMNAT1 mutation in two unrelated LCA ing evidence for genetic defects in noncoding regions such as deep patients by the identification of two distinct heterozygous par- intronic mutations or regulatory mutations in enhancers, promot- tial NMNAT1 deletions. Molecular characterization of the ers, or untranslated regions (UTRs). A striking example is a deep breakpoint junctions revealed a complex Alu-rich genomic ar- chitecture. Our study uncovered hidden genetic variation in intronic mutation of CEP290, accounting for 15% of congenital NMNAT1-associated LCA and emphasized a shift from coding blindness in Europe [Coppieters et al., 2010]. to noncoding regulatory mutations and repeat-mediated SVs in In 2012, four groups independently identified the gene encod- the molecular pathogenesis of heterogeneous recessive disor- ing nicotinamide nucleotide adenylyltransferase 1 (NMNAT1;MIM ders such as hereditary blindness. #608700) as a new LCA disease gene through WES [Chiang et al., Hum Mutat 00:1–9, 2015. Published 2015 Wiley Periodicals, 2012; Falk et al., 2012; Koenekoop et al., 2012; Perrault et al., 2012]. Inc.∗ Interestingly, all studies reported severe phenotypes with macular coloboma or atrophic macular lesions in almost all cases, repre- KEY WORDS: Leber congenital amaurosis; LCA9; senting one of the strongest genotype–phenotype correlations ever NMNAT1; 5 UTR variants; Alu-mediated deletions observed in LCA. Subsequent sequencing of the NMNAT1 coding re- gion in 740 additional prescreened LCA cases revealed homozygous or compound heterozygous mutations in 45 probands (6.1%) from different origins [Chiang et al., 2012; Falk et al., 2012; Koenekoop Additional Supporting Information may be found in the online version of this article. et al., 2012; Perrault et al., 2012]. The majority of NMNAT1 muta- †These authors contributed equally to this work. tions causing LCA are missense changes, predicted to alter protein ∗Correspondence to: Elfride DeBaere, Center for Medical Genetics Ghent, structure and/or abolish function based on the NMNAT1 crystal Ghent University Hospital, De Pintelaan 185, B-9000 Ghent, Belgium. E-mail: El- structure [Zhou et al., 2002], and in vivo as well as in vitro enzymatic [email protected] activity assays [Chiang et al., 2012; Falk et al., 2012; Koenekoop et al., 2012; Perrault et al., 2012; Siemiatkowska et al., 2014]. Until now, Contract grant sponsors: Research Foundation Flanders (FWO) (FWO 3G079711, only five nonsense, three frameshift, and one splice-site mutation FWO/KAN/1520913N); Ghent University Special Research Fund (BOF15/GOA/011); Bel- have been identified. In addition, the segregation of a read-through spo IAP project P7/43 (Belgian Medical Genomics Initiative: BeMGI).
∗∗ C 2015 The Authors. Human Mutation published by Wiley Periodicals, Inc. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distributioninanymedium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made. mutation in the original LCA9 family sustained the elucidation of the with PLINK [Purcell et al., 2007]. A healthy sibling (IV:3) was born last remaining LCA locus [Keen et al., 2003; Koenekoop et al., 2012]. later and underwent SNP genotyping later in the course of this study. Thus far, only coding NMNAT1 mutations have been reported. NMNAT1, which is localized to the nucleus, catalyzes NAD+ MPS in F1 synthesis in the last step of a salvage synthesis pathway that recycles nicotinamide back to NAD+. Both a fusion protein The four largest IBD regions common between both oldest af- (Ube4b/Nmnat1) from a spontaneous mouse model [Coleman and fected siblings of F1 (IV:1 and IV:2) were selected for MPS using Freeman, 2010] and NMNAT1 overexpression protect axons from the NimbleGen Sequence Capture 385K array. The array design degeneration [Zhu et al., 2013; Fang et al., 2014], opening thera- was based on the RefSeq (NCBI), AceView (NCBI), and RNAgene peutic perspectives for retinal neurodegeneration. (UCSC) resources and included all coding exons, 50 bp upstream Here, we aimed to identify the hidden genetic defect in a LCA9- intronic sequence, 20 bp downstream intronic sequence, UTRs (Ref- associated family. The absence of a NMNAT1 coding mutation Seq), and 500 bp upstream genomic region (putative promoter, Ref- prompted us to explore the UTRs, revealing a novel 5 UTR vari- Seq) (NCBI36/hg18). Two samples enriched with this method, one ant, c.-70A>T. Subsequently, a second adjacent 5 UTR variant, of which was the proband of F1, were sequenced on a Roche, Vil- c.-69C>T, was identified in an unrelated family with LCA and mac- voorde, Belgium GS FLX Titanium run. The CLC Genomics Work- ular involvement. Functional studies of both 5 UTRvariantssuggest bench (version 6.5; CLC bio, Aarhus, Denmark) was employed for a loss-of-function effect. Second, we aimed at unraveling the pseu- read mapping against the human genome reference sequence (NCBI, dohomozygosity of a coding NMNAT1 mutation in two unrelated GRCh37.p5), postmapping duplicate read removal, coverage anal- patients with a similar form of LCA, leading to the identification ysis, and quality-based variant detection. of two distinct heterozygous partial deletions of NMNAT1.Our study uncovers, for the first time, hidden genetic variation and a NMNAT1 Sequencing on gDNA and cDNA Level in F1 and F2 complex molecular pathogenesis of NMNAT1-associated LCA and is exemplary for congenital blindness and for other heterogeneous All probands underwent gDNA sequencing of exon 1, all coding autosomal-recessive diseases in general. regions and intron–exon boundaries of NMNAT1. cDNA sequenc- ing was carried out using intron-spanning amplicons in all individ- Materials and Methods uals of F1 and F2. PCR primers were designed with Primer3Plus [Untergasser et al., 2007]. Sanger sequencing was performed using Ethics Statement the BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosys- tems, ThermoFisher Scientific, Waltham, MA) on a 3730xl DNA This study was conducted following the tenets of the Declara- Analyzer (Applied Biosystems). tion of Helsinki and ethical approval was given by the local ethics committee (Ghent University Hospital). Evaluation of Sequence Changes Patients Mutation nomenclature uses numbering with the A of the ini- tiation codon ATG as +1(NM022787.3) (http://www.hgvs.org/ This study included 101 consenting subjects initially diagnosed mutnomen). All variants and individuals have been submitted to the with LCA (74) or early-onset retinal dystrophy (EORD) (27) and re- following locus-specific database: http://www.lovd.nl/NMNAT1. In ferred for molecular testing by an ophthalmologist and/or geneticist. silico evaluation of novel NMNAT1 variants was done with the All affected individuals were subjected to standard ophthalmologic Alamut mutation interpretation software v.2.3 (Interactive Biosoft- evaluation, among which visual acuity measurement, fundoscopy, ware, Rouen, France). Segregation analysis of disease alleles was per- ERG, infrared, visual field assessment, and autofluorescence imag- formed if possible, whereas gDNA obtained from unrelated healthy ing. individuals of Belgian origin (>173) was used as a control panel. Genomic DNA (gDNA) was extracted from leukocytes using theGentraPuregeneBloodKit(QIAGENBeneluxBV,Antwerp, Belgium), whereas RNA was extracted from leukocytes using the NMNAT1 and LZIC mRNA Quantification RNeasy Mini kit and the RNase-Free DNase set (Qiagen), followed Intron-spanning primer pairs for NMNAT1 and LZIC mRNA by cDNA synthesis (iScript cDNA Synthesis kit; Bio-Rad Laborato- quantification were designed using PrimerXL (http://www. ries NV, Temse, Belgium). A selection of lymphocyte cultures were primerxl.org). qPCR-based mRNA quantification was performed treated with puromycin to block nonsense-mediated decay. If avail- as previously described [Derveaux et al., 2010]. At least two ref- able, parents and/or siblings were also genotyped. Most probands erence genes were used for data normalization. Data analysis was previously underwent at least LCA chip analysis, testing 344 to performedwithqbase+ software (Biogazelle, Ghent, Belgium). P- 780 mutations in eight to 15 LCA and EORD genes (LCA chip values were calculated in qbase+ using a two-sided two-sample t-test Versions 2004–2012; Asper Ophthalmics, Tartu, Estonia, Estonia) with unequal variances. Six healthy controls were included in which (http://www.asperbio.com). no mutations were identified following Sanger sequencing of exon 1, all coding regions, and intron–exon boundaries of NMNAT1. Identity-by-Descent Mapping in F1
At the onset of this study, gDNA was available of two affected cases Bisulfite Sequencing (IV:1 and IV:2) of F1, originating from Niger. These individuals un- derwent genome-wide identity-by-descent (IBD) mapping using Bisulfite conversion of DNA was performed using the EZ DNA Affymetrix GeneChip Human Mapping 250K Nsp arrays accord- Methylation kit and the Human Methylated and Non-Methylated ing to manufacturer protocols (DNAVision, Charleroi, Belgium). DNA Set (Zymo Research, Freiburg, Germany), followed by Sanger Single-nucleotide polymorphism (SNP) genotypes were analyzed sequencing and analysis in SeqScape v2.5 (Applied Biosystems).
2 HUMAN MUTATION, Vol. 00, No. 0, 1–9, 2015 Luciferase Assays severe LCA with progressive macular involvement. Genome-wide IBD mapping using SNP arrays in the two oldest affected individuals Three different luciferase constructs were used: one for the NM- of F1, IV:1 and IV:2, revealed four common regions larger than NAT1 promoter and two for the NMNAT1 5 UTR. For the pro- 3 Mb (Supp. Table S1). Two of these could be partially excluded moter construct, we used the commercially available NMNAT1 later when the healthy sibling (IV:3) was born, as the latter child was pLightSwitch Promreportervector,whichwasmutatedtogen- homozygous for the same haplotype. Three regions remained in erate the c.-70A>T (F1) and c.-69C>T(F2)mutantconstructs total, together harboring 22.9 Mb and 321 genes. Interestingly, the (SwitchGear Genomics, La Hulpe, Belgium). For the 5 UTR con- largest region contained the LCA9 locus (Supp. Fig. S1) [Keen et al., structs, wild-type and mutant DNA oligos (Integrated DNA Tech- 2003]. nologies, BVBA, Leuven, Belgium) were cloned into the pGL3 re- MPS of the IBD regions was performed at the time when only IV:1 porter vectors pGL3-prom and pGL3-basic (Promega Benelux BV, and IV:2 were available, and included capturing and resequencing Leiden, The Netherlands), the latter containing a minimal CMV of all coding exons, intron–exon boundaries, UTRs, and putative promoter sequence. promoters of the four common IBD regions (Supp. Table S1). Of RPE-1 (human retinal pigment epithelial) cells were used to per- the 1.6 Mb target region, 87.5% was covered by the NimbleGen form functional analysis. They were grown in DMEM-F12 medium, Sequence Capture 385K platform design considering an offset of supplemented with 10% FBS and 1% penicillin/streptomycin (Ther- 100 bp. MPS generated 539,559 reads, of which 463,208 (86%) moFisher Scientific, Waltham, MA Gibco, ThermoFisher Scientific, unique reads were aligned against the human genome reference Waltham, MA). RPE-1 cells were seeded 16 hr prior to transfection sequence (NCBI, GRCh37.p5), resulting in a minimal coverage of 4 2 at a density of 4.10 cells/cm . Each experiment was performed in 10x and 20x for 93% and 85%, respectively, of the target regions. six replicates in 96-well culture plates. Cells were transfected with In total, 6,153 variants were detected with a coverage and variant 100 ng of DNA per well using fuGENE-HD (SwitchGear Genomics). allele frequency equal to or above 5x and 70%, respectively. Variants Forty-eight hours after transfection, cells were washed with PBS be- were prioritized based on their predicted effect on the protein (non- fore lysis and firefly luciferase and renilla luciferase activities were sense/frameshift > missense > silent) or splicing, and subsequently assessed by the Dual-Luciferase Reporter Assay System (Promega). subjected to segregation analysis in the family and control samples. Luciferase results are reported as relative light units (RLU, i.e., the However, this approach did not allow the identification of the causal ratio of the Renilla luciferase and firefly luciferase activities). Each defect. value is the mean of six independent experiments, and standard Following the identification of NMNAT1 as the causal disease error bar represent the standard error of the mean. gene for the LCA9 locus [Koenekoop et al., 2012], MPS data were revisited and subjected to an in-depth analysis of the five exons of Copy-Number Variation Analysis NMNAT1, including the UTRs. Supp. Table S2 lists all NMNAT1 variants identified in the proband of F1 through MPS. Of the nine For each of the five exons, a qPCR amplicon was designed variants, four were located in an intron, one upstream of the gene, (http://www.primerxl.org). Identified deletions were further delin- one in the 5 UTR and three in the 3 UTR. Two variants were novel, eated at the 5 and 3 end using additional qPCR amplicons. qPCR- with only one predicted to affect splicing, namely, c.-70A>T(Supp. based copy-number variation (CNV) analysis was performed as Fig. S2). This variant is located in exon 1, which is part of the 5 UTR, previously described, using two reference genes for normalization and segregates with the disease in the family (Fig. 1). of the relative quantities [D’Haene et al., 2010]. Data analysis was performedwithqbase+ software (Biogazelle). The junction prod- ucts were amplified with the iProof High-Fidelity DNA Polymerase (Bio-Rad), followed by Sanger sequencing using internal sequencing Second Adjacent NMNAT1 5UTR Mutation in Unrelated primers. Family Bioinformatic Analysis of the Deletion Junctions Sanger sequencing of the 5 UTR, all coding regions and intron– exon boundaries of NMNAT1 in 98 additional LCA (71) and EORD Bioinformatic analysis of the deletions junctions was performed (27) probands revealed potential pathogenic coding variants in five as previously described [Verdin et al., 2013]. LCA families (Supp. Table S3). Biallelic coding mutations were In short, known repetitive elements intersecting the breakpoints found in F5, F6, and F7. In the probands of F8 and F9, however, only were identified using the RepeatMasker track in the UCSC genome a single heterozygous variant was identified. Unfortunately, no phe- browser, followed by BLAST2 (NCBI) analysis to determine the notype data are available for these two patients. In both individuals, percentage of sequence identity between the different elements. The Sanger sequencing of the 5 UTR did not reveal a second mutation. presence of DNA sequences leading to non-B DNA conformations For F8, no DNA was available anymore to perform CNV analysis. In in the breakpoints was examined using the following tools: Re- F9, qPCR analysis did not reveal a CNV, pointing to other mutation peatAround [Goios et al., 2006], QGRS Mapper [Kikin et al., 2006], mechanisms in NMNAT1 or to the involvement of another retinal and nBMST [Cer et al., 2012]. Finally, the presence of 40 previously disease gene. described sequence motifs was analyzed with Fuzznuc [Rice et al., Interestingly, a homozygous 5 UTR variant, c.-69C>T, was 2000; Abeysinghe et al., 2003]. found in F2, an unrelated consanguineous family from Morocco (Table 1). This variant is located one nucleotide downstream of Results c.-70A>T; however, unlike c.-70A>T, it is not predicted to affect splicing (data not shown). The c.-69C>T variant was the only NM- First NMNAT1 5UTR Mutation in LCA9-Linked Family NAT1 variant identified in the proband of F2 and segregates with disease. Similar to F1, the affected individuals of F2 display LCA The starting point of this study was a consanguineous family with progressive macular involvement, which is in agreement with from Niger (F1) in which all affected individuals suffered from the typical NMNAT1-associated phenotype (Fig. 1).
HUMAN MUTATION, Vol. 00, No. 0, 1–9, 2015 3 Figure 1. Segregation analysis, associated phenotype, and mRNA quantification of the NMNAT1 5 UTR mutations. A: Segregation analysis of the c.-70A>T mutation in F1 and the c.-69C>T mutation in F2. B: Fundus picture of right eye of female patient IV:1 (F1) at age 7 years; note area of atrophy of outer retinal layers with some intraretinal pigment migration and highlighted luteal pigment in macula; mottled aspect of peripheral retinal pigment epithelium illustrates disease in outer retinal layers alternating with anatomically better preserved spots. C: Fundus picture of right eye of male patient IV:2 (F2) at age 21 years; note large area of total retinal and choroidal atrophy in macular area; this area was smaller with atrophy limited only to outer retinal layers but not choroid at age 7 years (data not shown); also note grayish area of outer retinal atrophy in retinal periphery, with limited intraretinal pigment migration. D: qPCR quantification of NMNAT1 mRNA abundance on lymphocyte cDNA of all available family members of F1 and F2, and six healthy controls. The abundance of NMNAT1 mRNA was significantly lower in the affected individuals of F1 in comparison with the controls (P = 5,473E-5). No significant difference was observed however in the affected individuals of F2 in comparison with controls (P = 0.072) (unpaired t-test). Abbreviations used: C, control; F, family; M, mutant allele, +, wild-type allele; UTR, untranslated region.
Functional Assessment of c.-70A>T and c.-69C>T c.-70A>T was predicted to affect splicing, lymphocyte cultures of F1 were treated with and without puromycine, a protein synthe- NMNAT1 cDNA sequencing and quantification sis inhibitor used to prevent nonsense-mediated mRNA decay. No splicing defects were identified. However, haplotype analysis at the Subsequent Sanger sequencing was performed on lymphocyte- cDNA level of three NMNAT1 variants previously found at the gDNA derived cDNA from all available family members of F1 and F2. As level revealed allele-specific expression in the parents and the healthy
4 HUMAN MUTATION, Vol. 00, No. 0, 1–9, 2015 Table 1. NMNAT1 5UTR Mutations and Copy-Number Variations Explaining Hidden Genetic Variation in This Study
Allele 1 Allele 2 Family Exon cDNA Protein Reference Exon/intron cDNA Protein Reference