CODAS Syndrome Is Associated with Mutations of LONP1, Encoding Mitochondrial AAA+ Lon Protease ARTICLE CODAS Syndrome Is Associated with Mutations of LONP1, Encoding Mitochondrial AAAþ Lon Protease Kevin A. Strauss,1,2,3,14,* Robert N. Jinks,3,14 Erik G. Puffenberger,1,3,14 Sundararajan Venkatesh,4 Kamalendra Singh,4,6 Iteen Cheng,5 Natalie Mikita,5 Jayapalraja Thilagavathi,4 Jae Lee,4 Stefan Sarafianos,6 Abigail Benkert,1,3 Alanna Koehler,3 Anni Zhu,3 Victoria Trovillion,3 Madeleine McGlincy,3 Thierry Morlet,7 Matthew Deardorff,8,9 A. Micheil Innes,10 Chitra Prasad,11 Albert E. Chudley,12,13 Irene Nga Wing Lee,5 and Carolyn K. Suzuki4,14 CODAS syndrome is a multi-system developmental disorder characterized by cerebral, ocular, dental, auricular, and skeletal anomalies. Using whole-exome and Sanger sequencing, we identified four LONP1 mutations inherited as homozygous or compound-heterozygous combinations among ten individuals with CODAS syndrome. The individuals come from three different ancestral backgrounds (Amish- Swiss from United States, n ¼ 8; Mennonite-German from Canada, n ¼ 1; mixed European from Canada, n ¼ 1). LONP1 encodes Lon protease, a homohexameric enzyme that mediates protein quality control, respiratory-complex assembly, gene expression, and stress responses in mitochondria. All four pathogenic amino acid substitutions cluster within the AAAþ domain at residues near the ATP-bind- ing pocket. In biochemical assays, pathogenic Lon proteins show substrate-specific defects in ATP-dependent proteolysis. When ex- pressed recombinantly in cells, all altered Lon proteins localize to mitochondria. The Old Order Amish Lon variant (LONP1 c.2161C>G[p.Arg721Gly]) homo-oligomerizes poorly in vitro. Lymphoblastoid cell lines generated from affected children have (1) swollen mitochondria with electron-dense inclusions and abnormal inner-membrane morphology; (2) aggregated MT-CO2, the mtDNA-encoded subunit II of cytochrome c oxidase; and (3) reduced spare respiratory capacity, leading to impaired mitochondrial pro- teostasis and function. CODAS syndrome is a distinct, autosomal-recessive, developmental disorder associated with dysfunction of the mitochondrial Lon protease. Introduction Cerebral, ocular, dental, auricular, skeletal (CODAS) syn- drome (MIM 600373) was first described in 1991 as a distinctive constellation of developmental delay, craniofa- cial anomalies, cataracts, ptosis, median nasal groove, delayed tooth eruption, anomalous cusp morphology, malformed helices, hearing loss, short stature, delayed epiphyseal ossification, metaphyseal hip dyplasia, and vertebral coronal clefts.1 CODAS is a rare disease; only three additional cases (from Brazil, Canada, and France) were reported between 1995 and 2010.2–4 Autosomal- recessive inheritance of CODAS syndrome was first sug- gested by Shebib et al.,1 who identified the phenotype within an endogamous Mennonite community, and it was corroborated by the observations of Innes et al.3 We used whole-exome and Sanger sequencing to study ten individuals with CODAS syndrome from three distinct ancestral backgrounds: Amish-Swiss from the United States (n ¼ 8), Mennonite-German from Canada (n ¼ 1), and 1Clinic for Special Children, Strasburg, PA 17579, USA; 2Lancaster General Hosp dations of Behavior Program, Franklin and Marshall College, Lancaster, PA Genetics, New Jersey Medical School, Rutgers, The State University of New Jerse University, Cleveland, OH 44106, USA; 6Department of Molecular Microbiolo Missouri, Columbia, Columbia, MO 65201, USA; 7Auditory Physiology and Psy ton, DE 19803, USA; 8Division of Human Genetics, The Children’s Hospital of P man School of Medicine, University of Pennsylvania, Philadelphia, PA 19104 Research Institute, Cumming School of Medicine, University of Calgary, Ca Pediatrics, Children’s Health Research Institute and Western University, Lond University of Manitoba, Winnipeg, MB R3A 1S1, Canada; 13Department of MB R3A 1S1, Canada 14These authors contributed equally to this work *Correspondence: kstrauss@clinicforspecialchildren.org http://dx.doi.org/10.1016/j.ajhg.2014.12.003. �2015 by The American Societ The Amer mixed European (German, Scottish, Irish, English) from Canada (n ¼ 1). Our analysis revealed four different mutations within LONP1 (RefSeq accession number NM_004793.3 [isoform 1]; [MIM 605490]). These muta- tions were either homozygous (LONP1 c.2161C>G [p.Arg721Gly] and LONP1 c.2026C>T [p.Pro676Ser]) or compound heterozygous (LONP1 c.1892C>A/c.2171C>T [p.Ser631Tyr/p.Ala724Val]) (Table 1). In humans, LONP1 encodes the ATP-dependent mitochondrial Lon protease, which belongs to the AAAþ superfamily of ATPases associ- ated with various cellular activities. AAAþ proteins share a conserved ATPase module and mediate diverse cellular pro- cesses such as DNA replication, membrane fusion, signal transduction, and transcriptional regulation.5–8 Lon is an ATP-driven proteolytic machine and is highly conserved from archaea to mammals. In human mito- chondria, Lon is the simplest of four energy-dependent proteolytic systems that include ClpXP, m-AAA, and i-AAA.9 Mutations of CLPP (MIM 601119), which encodes the proteolytic component of ClpXP, are associated with ital, Lancaster, PA 17602, USA; 3Department of Biology and Biological Foun- 17603, USA; 4Department of Microbiology, Biochemistry, and Molecular y, Newark, NJ 07103, USA; 5Department of Chemistry, Case Western Reserve gy and Immunology, Christopher Bond Life Sciences Center, University of choacoustics Research Laboratory, duPont Hospital for Children, Wilming- hiladelphia, Philadelphia, PA 19104, USA; 9Department of Pediatrics, Perel- , USA; 10Department of Medical Genetics and Alberta Children’s Hospital lgary, AB T2N 1N4, Canada; 11Medical Genetics Program, Department of on, ON N6C 2V5, Canada; 12Department of Pediatrics and Child Health, Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, y of Human Genetics. All rights reserved. ican Journal of Human Genetics 96, 121–135, January 8, 2015 121 mailto:kstrauss@clinicforspecialchildren.org http://dx.doi.org/10.1016/j.ajhg.2014.12.003 http://crossmark.crossref.org/dialog/?doi=10.1016/j.ajhg.2014.12.003&domain=pdf Table 1. Four Pathogenic LONP1 Alleles Among Ten Individuals with CODAS Syndrome Origin of Sample n Ancestry LONP1 Variants Previously ReportedAllele 1a Allele 2a Pennsylvania 8 Amish-Swiss Arg721Gly Arg721Gly Canada 1 Mennonite-German Pro676Ser Pro676Ser Shebib et al., 19911 Canada 1 Mixed European Ser631Tyr Ala724Val aIn every case, genotyping showed allele segregation consistent with recessive inheritance. Perrault syndrome type 3 (MIM 614129),10 and mutations of SPG7, which encodes a subunit of m-AAA, are linked to hereditary spastic paraplegia (MIM 607259).11 In humans, Lon functions as a protease, chaperone, and DNA-binding protein. It contributes to mitochondrial proteostasis by eliminating misfolded and oxidatively damaged pro- teins12–14 and supports cell viability during proteotoxic, hypoxic, and endoplasmic-reticulum stress.15–17 Lon also degrades key regulatory proteins of mitochondrial metabolism, for example ALAS-1 and StAR, rate-limiting enzymes of heme and steroid hormone biosynthesis, respectively.18,19 Independent of its proteolytic activity, Lon influences function of the mitochondrial genome and respiratory chain. Its chaperone-like function promotes assembly of respiratory-chain complexes within the mitochondrial inner membrane,16,20,21 and knockout of Lon-encoding genes in yeast and mammals leads to mtDNA deletions and depletion, respectively.15,22,23 Lon binds single- stranded mtDNA in a sequence- and strand-specific manner24–27 and associates preferentially with the mtDNA control region where replication and transcription are initiated.24 Here, we establish a link between LONP1 and CODAS syndrome in humans and discuss potential mechanisms linking the dysfunction of mitochondrial Lon to diverse clinical manifestations. Subjects and Methods Subjects We studied a total of ten individuals from two countries (United States, n ¼ 8 and Canada, n ¼ 2). All had the core develop- mental, ocular, dental, auricular, and skeletal features of CODAS syndrome.1–4 Detailed phenotype studies were conducted on four children homozygous for LONP1 c.2161C>G, a founder mutation among the Pennsylvania Amish (Figures 1 and 2). The study was approved by the institutional review board of Lancaster General Hospital. Parents consented to research on their affected children on their behalf and signed a separate con- sent to reproduce photographs. Affected children were evaluated by routine clinical methods, including bone radiographs, 1.5 T magnetic resonance imaging, slit-lamp and dilated indirect ophthalmoscopy, middle-ear tympanometry, distortion-product otoacoustic emissions, evoked auditory brainstem responses, and direct laryngoscopy. 122 The American Journal of Human Genetics 96, 121–135, January 8 Molecular Genetics Genetic mapping of CODAS syndrome was performed as previ- ously described.28–33 In brief, SNP genotyping was performed with the GeneChip Mapping 10K Assay Kits (Affymetrix). Data were analyzed in Microsoft Excel spreadsheets that were custom formatted at the Clinic for Special Children. SNP positions came from Affymetrix genome-annotation files, and genotype data came from the Affymetrix GeneChip Human Mapping 10K Xba 142 arrays. Data analyses were designed for rapid identification of genomic regions that were identically homozygous between affected individuals. Exome sequencing was performed as previously described.34–36 A solution-hybrid-selection methodology was used for isolating exomic DNA, which was subjected to sequencing on the Illumina GA-II platform. In brief, DNA oligonucleotides, corresponding to 120 bp of target sequence flanked by 15 bp of universal primer sequence, were synthesized in parallel on an Agilent microarray, then were cleaved from the array. The oligonucleotides were amplified by PCR, then transcribed in vitro in the presence of bio- tinylated uridine-50-triphosphate so that single-stranded RNA ‘‘bait’’ would be generated. Genomic DNA was sheared, ligated to Illumina sequencing adapters, and selected for lengths between 200 and 350 bp. This ‘‘pond’’ of DNA was hybridized with an excess of bait in the solution. The ‘‘catch’’ was pulled down by magnetic beads coated with streptavidin, eluted, and sequenced on the Illumina GA-II. Massively parallel sequencing data were processed by the Sequencing Platform at the Broad Institute via two pipelines. The first, called ‘‘Picard,’’ utilized the reads and qualities produced by the Illumina software and produced a single BAM file represent- ing each sample. The second pipeline, called Genome Sequencing Analysis, then performed post-processing and analysis of the data, including SNP identification, small-insertion and -deletion identi- fication, local realignment of insertion- or deletion-containing reads, gene annotation, and filtering with common polymor- phisms (1000 Genomes, dbSNP build 129). The details of the sequencing-data processing have been described elsewhere.37 For inter-exome analyses, variant call data were imported into a File- Maker Pro database, and rare and novel variants shared among affected individuals were queried. Molecular Modeling A homology model of the human mitochondrial Lon hexamer was generated from the crystal structure of B. subtilis Lon (PDB 3M6A).9,38 In brief, we generated a homology model of a single Lon subunit that lacks its amino-terminal mitochondrial-targeting sequence (MTS) (see Figures 3B–3D) by using three independent molecular modeling software tools: Prime (Schrodinger), Modeller (see Web Resources), and the I-TASSER online server at the , 2015 Figure 1. CODAS Phenotype (A and B) A 5 year-old child who is homozygous for LONP1 c.2161C>G (p.Arg721Gly) and who has characteristic physical features of CODAS syndrome. Features include (A) broad skull and flattened midface; ptosis; grooved nasal tip; anteverted nares; and (B) helix hypoplasia (‘‘crumpled’’ ears). (C) Laryngoscopy of an affected child revealed airway narrowing secondary to paretic, atrophic vocal cords (inset: normal vocal cords for comparison). (D) Dense bilateral nuclear cataracts develop rapidly between 2 and 6 months of age in all CODAS-syndrome-affected Amish individuals. (E and F) Scoliosis and genu valgum progress with advancing age. (G) Despite previous reports of developmental delay in children with CODAS syndrome, timely ophthalmologic and audiologic inter- vention appear to be critical developmental determinants; at age 5 years our oldest subject has a vocabulary of more than 100 words, reads at the first grade level, communicates effectively by signing, and has advanced drawing and writing skills. University of Michigan. We obtained the final model of the Lon subunit by averaging three models. We obtained the Lon hexamer by superimposing the structure of one subunit onto the crystal structure of B. subtilis Lon (PDB 3M6A). cDNA Constructs Human LONP1 cDNA (RefSeq NM_004793.3) that was amplified from ARPE-19 (human adherent retinal pigment epithelium-19) total RNA by RT-PCR (SuperScript II reverse transcriptase, Life Tech- nologies; Phusion DNA polymerase, New England BioLabs)35 was cloned into pENTR/D-TOPO (Life Technologies), and CODAS mutations were introduced by site-directed mutagenesis (Quik- Change II, Agilent) (Table S2). For mammalian expression of Lon proteins with C-terminal V5 or FLAG tags, the LONP1 cDNA was cloned in frame into pcDNA3.2/V5-DEST or pcDNA3.2/FLAG- DEST by recombination from pENTR/D-TOPO/LONP1 constructs with LR Clonase II (Life Technologies). For bacterial expression, cDNAs encoding Lon proteins that lacked the first 114 amino acids (containing the predicted MTS) were sub-cloned into pET24c(þ) as well as into pProEX-1; the latter constructs encoded a His tag fused to the N terminus.25 All amino acid numberings were assigned on the basis of the full-length Lon precursor. Purification of Recombinant Lon Wild-type and pathogenic Lon variants were overexpressed in Ro- setta (DE3) (Novagen) and purified as previously described.39–41 The Amer Lon constructs were cloned into either pET24c(þ) (Figure 4) or ProEX (hexa-histidine tag26) (Figure S1). For Figure 4, Lon-positive fractions were pooled and concentrated, and protein concentra- tions were determined via the Bradford assay, for which BSA was used as a standard and UV absorbance was set at 280 nm, and us- ing the equation ε280 ¼ W(5500) þ Y(1490) þ C(125), where ε280 is the molar absorptivity at 280 nm and W, Y, and C are the numbers of tryptophan, tyrosine, and cysteine residues, respectively.42,43 Protein concentrations were further corrected on the basis of pro- tein-band optical density normalized to LonWT in Coomassie- stained SDS-PAGE gels (Image J). Recombinant steroidogenic acute regulatory protein (StAR) and mitochondrial transcription factor A (TFAM) were purified as described previously.41,44 Enzymatic Assays Peptidase and Protease Assays Fmoc-protected amino acids, Boc-Abz, Fmoc-protected Lys-Wang resin, and HBTU (N,N,N0,N0-tetramethyl-O-(1H-benzotriazol-1-yl) uronium hexafluorophosphate, O-(benzotriazol-1-yl)-N,N,N0,N0- tetramethyluronium hexafluorophosphate) were purchased from Advanced ChemTech and NovaBiochem. Peptide substrate S1 was purchased from Genscript. Peptide synthesis of S2 was per- formed with Fmoc solid-phase synthesis methodologies as described previously.40 Peptidase activity was measured on a Fluo- romax-4 spectrofluorimeter (Horiba Group). All reactions mixtures contained 50 mM HEPES (pH 8.0), 5 mM Mg(OAc)2, 1 mM DTT, ican Journal of Human Genetics 96, 121–135, January 8, 2015 123 Figure 2. Radiographic Features of CODAS Syndrome (A) Sagittal T1 (upper) and axial T2 (lower) images of an affected 4-year-old Amish child show mild diffuse cortical atrophy, an immature pattern of myelination, and hypoplasia of the corpus callosum. (B) Skeletal radiographs show metaphyseal dysplasia most evident at hip joints (upper) and valgus knees with both hypoplasia and delayed ossification of epiphyses (lower). (C) Severe scoliosis is evident early in childhood, and coronal clefts are observed at various levels of the vertebral column. (D–G) Prenatal imaging of a 32-week-old fetus with CODAS syndrome revealed polyhydramnios, a two-vessel umbilical cord, midfacial hypoplasia, (D) crumpled helices, (E) a balanced atrioventricular canal, (F) common atrium, omphalocele, (G, arrow) and absence of left lower long bones, which were replaced by a relatively well-formed foot attached directly to the hip. (G, arrowhead) The asterisk in (G) marks the left palpebral fissure for orientation. 250 nM wild-type and mutant Lon proteins, and 0.5 mM S3. To avoid problems from the inner filter effect, we used S3, a mixed substrate consisting of 10% fluorescently labeled peptide S2 [Y(NO2)RGITCSGRQK(Abz)] and 90% non-fluorescent analog of the peptide S1 (YRGITCSGRQK(Ac)). After the mixture equili- brated at 37�C for 1 min, the reaction was initiated with 1 mM ATP. We measured the amount of hydrolyzed peptide by deter- mining the maximum fluorescence generated per micromole of peptide after complete digestion by trypsin. The steady-state rate of the reaction was determined from the tangent of the linear portion of the time course. This rate was converted to an observed rate constant (kobs) by division of the rate by the enzyme concen- tration. At least three identical experiments were performed. Peptidase assays using the fluorescent dipeptide substrate rhodamine 110, bis-(CBZ-L-alanyl-L-alanine amide) (AA2-Rh110, Anaspec), were performed in quadruplicate (20 ml) in 384-well plates as previously described.39 Lon (800 nM monomer) or no- enzyme controls were incubated in a reaction buffer (150 mM NaCl, 50 mM HEPES [pH 8.0], and 10 mM MgCl2) containing AA2-Rh110 (6 mM) and ATP (2 mM) for 3 hr at 37 �C, after which fluorescence was measured at an excitation and emission of 485 and 535 nm, respectively, with a Perkin Elmer Victor3 V. 124 The American Journal of Human Genetics 96, 121–135, January 8 The relative fluorescence units (RFU) of the background (no- enzyme control) measurements were subtracted, and the resultant values were normalized to percent activity of the no-drug reac- tions. Data were fit to 4-parameter dose-response curves using GraphPad Prism 6, and the error bars represent the SD of four repli- cate reactions. At least three identical experiments were per- formed. StAR and TFAM purification and degradation assays were performed as previously described.19,41,44,45 ATPase Assay The ATPase activity of wild-type and mutant Lon variants was measured in quadruplicate reactions (5 ml) in a 384-well plate via the ADP-Glo Assay (Promega) according to the manufacturer’s recommendations. All steps were performed at room temperature. Purified human Lon (400 nM monomer) was incubated in a reaction buffer (40 mM Tris-HCl [pH 7.5], 20 mM MgCl2, and 0.1 mg/ml BSA) with Ultra Pure ATP (1 mM; Promega) for 60 min. ADP-Glo� Reagent (5 ml) was added for 40 min followed by Kinase Detection Substrate (10 ml) for 60 min, which was linear with [ADP]. Relative luminescence units (RLU) were measured on a Perkin Elmer Victor3 V, and background (no-enzyme control) RLU were subtracted. At least three identical experiments were performed. , 2015 Figure 3. Location of Amino Acids That Are Mutated within Mitochondrial Lon in CODAS Syndrome (A) Domain structure of the human Lon subunit. Mitochondrial targeting sequence (MTS), substrate recognition/binding (N) domain, ATPase domain (AAAþ), and a pro- tease domain (P). Red arrows indicate the position of pathogenic CODAS mutations. (B) Homology model of human mitochon- drial Lon. The model shows the position, within a single Lon subunit (shown in amber), of amino acids that are altered in CODAS syndrome. (C and D) The positions of amino acids Pro676 and Arg721 (yellow), which are altered in proteins encoded by CODAS homozygous alleles. ADP is shown occu- pying the ATP- and ADP-binding pocket (a green dashed line represents a proline-induced kink in the helix; a yellow dotted line represents a salt bridge; starbursts represent hydrophobic interac- tions). (C) The Pro676-induced kink pro- motes hydrophobic interactions between Ala670, Leu667, Leu696, and Val716. (D) The positions of Arg721 as well as Pro676, Glu717, and the ATP- and ADP-binding site are located on the same Lon subunit (amber). Arg721 and Glu717 form salt bridges with Glu654 and Lys517, respectively, which are located on the adjacent Lon subunit (green). Cell Culture ARPE-19 cells, human embryonic kidney (HEK293T) cells, HeLa cells (ATCC, Manassas, VA), and mouse auditory sensory epithe- lium (UB/OC-2) cells (gift from Matthew Holley, University of Sheffield, UK) were cultured in DMEM:F12, DMEM, Eagle’s MEM, or MEM with GlutaMax and 50U/ml interferon-gamma (Life Technologies), respectively, each with 10% fetal bovine serum (FBS; Sigma or Life Technologies) at 37�C/5% CO2. Ep- stein-Barr virus (EBV)-transformed B-lymphoblastoid cell lines (LCLs) were generated from two Amish CODAS-syndrome- affected probands homozygous for LONP1 c.2161C>G (p.Arg721- Gly) as well as from their respective heterozygous parents (Lineberger Comprehensive Cancer Center, University of North Carolina). LCLs were cultured in RPMI (ATCC formulation) with 15% FBS at 37�C and 5% CO2. Genotyping of LCLs was performed by Sanger sequencing (The Pennsylvania State University Nucleic Acid Facility). In brief, SuperScript III reverse transcriptase (Life Technologies) and oligo dT primers were used for reverse transcrip- tion of cDNA from LCL total RNA (Trizol, Life Technologies). cDNA was then amplified by the polymerase chain reaction (PCR) with LONP1-specific primers (Table S1) and purified for DNA sequencing with QIAGEN PCR purification spin columns. Lon knockdown in HeLa, PC3, and HEK293T cells was performed as previously described46 (and as also described in Supplemental Data). Mitochondrial Localization Wild-type and mutant Lon variants were overexpressed by tran- sient transfection (FuGENE 6, Promega) in ARPE-19, HeLa, and UB/OC-2 cells cultured on glass coverslips. After 40–48 hr, cells were incubated with MitoTracker Red CMXRos (300 nM) at 37�C/5% CO2 for 30 min, washed with PBS (pH 7.2), and fixed in 4% paraformaldehyde in PBS for 15 min at 37�C. Cells were per- meabilized in acetone at �20�C for 5 min, washed with PBS, and The Amer blocked in 1% bovine serum albumin (BSA) with 0.2% Triton X- 100 in PBS for 10 min. Immunofluorescence detection of Lon- V5 fusion proteins was conducted as described elsewhere.35 Immunoblotting and Immunoprecipitation Cells were lysed in the following buffers as indicated: in RIPA buffer (50 mM Tris-HCl [pH 8.0], 150 mM NaCl, 1.0% octylphe- noxypolyethoxyethanol CA-630, 0.5% sodium deoxycholate, and 0.1% sodium dodecyl sulfate) supplemented with 1 mM EDTA, 2 mM NaF, 1 mM Na3VO4, and protease-inhibitor cocktail (Roche); Triton X-100 buffer (50 mM Tris, [pH 7.5], 300 mM NaCl, and 0.5% Triton X-100) supplemented with Halt protease- and phosphatase-inhibitor cocktail (ThermoScientific); or urea buffer (50 mM triethylammonium bicarbonate buffer [TEAB, Sigma] [pH 8.5], 8 M urea) supplemented with protease-inhibitor cocktail and phosphatase-inhibitor cocktail (Sigma). When the RIPA or Triton X-100 buffer was used, cells were incubated for 15 min on ice, then centrifuged at 14,000 rpm for 15 min at 4�C so that lysates would be cleared. For the urea buffer, cells were son- icated with microprobe for 15 s followed by a 30 s pulse, and this was repeated 3–5 times before centrifugation at 4�C cleared lysates. The protein concentration of cell extracts was measured with the Bradford assay and then normalized.42 Immunoblotting and co-immunoprecipitation were performed as described else- where;35 washes following immunoprecipitation contained 5% Tween-20. The following antibodies were employed: rabbit anti- Lon (1:400, custom made and affinity purified as previously described)24 or rabbit anti-Lon (1:100, Novus, H00009361- D01P); rabbit anti-ClpX (1:3,000, custom made by Dr. Irene Lee); rabbit anti-aconitase (1:200, provided by Dr. Luke Szweda, Okla- homa Medical Research Foundation); mouse anti-ClpP (1:4,000, Abcam 56455); rabbit anti-Pink1 (1:500, Abcam 23707); rabbit anti-mt-cytochrome oxidase subunit II (1:5,000, Abcam 91317); ican Journal of Human Genetics 96, 121–135, January 8, 2015 125 Figure 4. Enzymology, Cellular Expression, and Homo-oligomerization of Pathogenic Lon Proteins (A) ATP-dependent degradation of the fluorogenic peptide substrate S3. The rates of S3 cleavage were determined with recombinant Lon proteins (250 nM monomer) and 0.5 mM of S3 (0.5-fold Km) in buffer. The degradation reactions were initiated by the addition of Mg-ATP at 37�C for 600–900 s. S3 cleavage was monitored as an increase in fluorescence at an excitation and emission of 320 nm and 420 nm, respectively (three or more replicates). Error bars represent SEM. (B) StAR (5.6 mM) or TFAM (5 mM) were combined in buffer with Mg-ATP, and reactions were initiated by the addition of Lon (500 nM) at 37�C. Aliquots were removed at time points, reactions were terminated with 53 reducing sample buffer, and aliquots were analyzed by SDS-PAGE and Coomassie Brilliant Blue staining. ImageJ was used for determining the percent StAR and TFAM degraded after a 60 min incubation period with the same Lon variant (no replicates). (C and D) Overexpression of (C) wild-type Lon-V5 and (D) p.Arg721Gly-V5 (green ¼ V5 immunofluorescence) in ARPE-19 cells counter- stained with Mitotracker Red (red) and DAPI (blue). Nontransfected cells display only red Mitotracker fluorescence. (E) Coimmunoprecipitation (IP) of Lon-V5 by Lon-FLAG (co-overexpressed in HEK293T cell lysates) with anti-FLAG M2 antibody. Immunoblotting (IB) with anti-FLAG and V5 antibodies followed. OXPHOS Antibody Cocktail (1:500, Mitoscience 604); goat anti- actin (1:2,000, 1615; Santa Cruz); mouse anti-b-actin (Sigma, 1:1,000,000); rabbit anti-mt-cytochrome B (1:1,000, Aviva Systems Biology ARP50256_P050); mouse monoclonal anti-V5 (1:5,000, Life Technologies) or mouse anti-FLAG M2 (1:1,000 Sigma); and HRP-conjugated secondary antibodies (Cell Signaling Technologies or Santa Cruz Biotechnology). Immunoblots were developed by chemiluminescence. Transmission Electron Microscopy LCLs were prepared for transmission electron microscopy (TEM) as described previously.39,47 In brief, cells were pelleted by centrifuga- tion at 1,000 rpm for 5 min, washed with PBS, re-pelleted, and re- suspended in 2.5% glutaraldehyde and 4% paraformaldehyde in 0.1 M PBS overnight at 4�C. LCLs were then post-fixed in 1% osmium tetroxide for 1 hr at 20�C and dehydrated with ethanol; infiltration with Spurr’s resin (Electron Microscopy Sciences) fol- lowed. Infiltrated LCLs were collected by centrifugation and embedded directly into BEEM capsules, which were gently centri- fuged so that the cells would be driven to the bottom of the capsule. Ultrathin sections (70-90 nm) were cut with glass knives, 126 The American Journal of Human Genetics 96, 121–135, January 8 mounted on naked copper grids (200 mesh), and stained with a saturated solution of uranyl acetate in 50% ethanol for 20 min, fol- lowed by Reynold’s lead citrate for 15 min. Differences in the fre- quency of abnormal mitochondrial morphology between proband and parent cells were analyzed by two-way ANOVA with Bonfer- roni post-hoc correction (SPSS, IBM). Quantitative PCR (qPCR) Determination of mtDNA Copy Number Genomic DNA (100 ng) from peripheral blood lymphocytes and LCL cells was used for amplification of both the mtDNA- encoded MT-CYB gene and the nuclear-DNA-encoded APP gene as a reference control. Genomic DNA was amplified in reactions (20 ml; Applied Biosystems; universal PCR master mix) with the following Taqman primers. Forward: APP, 50-TTTTTGTGTGCT CTCCCAGGTCT-30; and MT-CYB, 50-GCCTGCCTGATCCTCCAA AT-30). Reverse: APP 50-TGGTCACTGGTTGGTTGGC-30; and MT CYB, 50-AAGGTAGCGGATGATTCAGCC-30. TaqMan probes: APP, 50-[6FAM]CCCTGAACTGCAGATCACCAATGTGGTAG[TAM]-30; and MT-CYB, 50-[6FAM]CACCAGACGCCTCAACCGCCTT [TAM]-30 An AB 7500 RT-PCR system (Applied Biosystems) and , 2015 standard qPCR conditions were employed: 95�C for 10 min, 40 cycles at 95�C for 15 s, and 60�C for 1 min. All reactions were performed in triplicate. Relative quantitation of the mtDNA copy number was calculated by the DDCt method. Mitochondrial Oxygen-Consumption Measurements Proband and parent cells (1 3 105/well) in a 100 ml unbuffered assay medium (DMEM with 11 mM glucose, 1 mM pyruvate, and 2 mM glutamine) were plated in poly-D-lysine-coated Seahorse XF 24-well plates (50 mg/ml) and centrifuged so that cells were adhered to the plates. After the uniformity of cells plated in each well was ensured, plates were transferred to a 37�C incubator without CO2 for 30 min, 575 ml warm assay medium was added, and then plates were left in the incubator for another 20 min. Each plate was transferred to the XF24 Analyzer, and the oxygen consumption rate (OCR) was measured at baseline. Sequential injections were then performed with the following: oligomycin (1 mM), an inhibitor of complex V (the ATP synthase); FCCP (750 nM), an uncoupler that dissipates the membrane potential across the mitochondrial inner membrane; and rotenone (1 mM), an inhibitor of complex I (NADH dehydrogenase), which blocks the transfer of electrons from complex I to ubiquinone. The spare respiratory capacity (SRC), which is the maximum OCR after FCCP injection, was calculated as a percentage of baseline OCR. The mean, standard deviation, and statistical significant difference between proband and mother cells were calculated by Student’s t test in Microsoft Excel. Results Molecular Genetic Studies We initiated a genetic mapping study involving two Amish children who were affected with CODAS syndrome. Affy- metrix GeneChip Mapping 10K DNA microarrays were initially used for genotyping, but no significant regions of homozygosity were detected. Exome sequencing was performed, and variants were filtered so that those with minor allele frequencies (MAFs) > 1% in either 1000 Genomes or the Exome Variant Server or with MAF > 10% in our own internal Amish variant database were excluded. Analysis of variants detected in CODAS-affected individuals revealed a single shared homozygous allele, LONP1 c.2161C>G (p.Arg721Gly), located at physical position chr19:5,694,557 (NCBI GRCh37, hg19) The gene mapped to a region of chromosome 19 with poor coverage on the Affymetrix 10K array, hence our inability to conclusively map the phenotype. The presence of LONP1 c.2161C>G was confirmed by Sanger sequencing, and we verified proper segregation by genotyping parents and siblings. Using an unlabeled probe to genotype 576 Amish controls by high-resolution melt analysis, we found a surprisingly high population allele fre- quency of 5.9%. We next acquired additional CODAS-syn- drome DNA samples accompanied by clinical data from collaborators in Canada and designed primers to amplify and sequence all 18 exons of LONP1. The first sample was from a Mennonite child from Manitoba, Canada, who was the subject of a previous publication.1 She was The Amer found to be homozygous for LONP1 c.2026C>T (p.Pro676Ser). The second sample was from a child of mixed European ancestry (German, Scottish, Irish, and En- glish) from Ontario, and she was found to be compound heterozygous for LONP1 c.1892C>A (p.Ser631Tyr) and c.2171C>T (p.Ala724Val)(Table 1). Phenotype Associated with CODAS Syndrome Eight children from the Old Order Amish community of Pennsylvania were homozygous for LONP1 c.2161C>G (Table 2). Three died of laryngeal obstruction in the first days of life, and one with a tracheostomy died from pneu- monia during early infancy. Four remaining children (ages 0.6–14.2 years, two female) required tracheostomies to sur- vive infancy; only one child was subsequently decannu- lated. All three surviving children who underwent sedated laryngoscopies had paretic, atrophic vocal cords, glottic narrowing, chronic sialorrhea, and swallowing dysfunc- tion. Three of four surviving affected children were nour- ished exclusively by gastrostomy tube. LONP1 c.2161C>G homozygotes had physical features characteristic of CODAS syndrome.1–4 These included broad skull and flattened midface, helix hypoplasia (‘‘crumpled’’ ears), ptosis, grooved nasal tip, anteverted nares, and with advancing age, short stature, scoliosis, genu valgus, and pes valgus (Figures 1A–1F). Teeth erupted late with cusp-tip extensions. Dense bilateral nuclear cata- racts developed rapidly between 2 and 6 months of age. Audiological testing showed impaired tympanic mem- brane mobility (type B pattern), normal otoacoustic emissions and neural synchrony, and low-frequency conductive hearing loss. The two oldest individuals had a mixed pattern involving mild-to-moderate sensory hear- ing loss in the 2–4 kHz range. All affected children had hypotonia, developmental delay, and variable intellectual disability, some of which was remediable. One child with delayed ambulation was able to walk independently when knee-ankle-foot orthoses were applied at age 4 years. This same child had early audi- ological intervention (myringotomy tubes and amplifica- tion). By age 5, he had a vocabulary of more than 100 words, could read at the first grade level, could communi- cate effectively by signing, and had advanced drawing and writing skills (Figure 1G). Neuroimaging of one affected child (age 4 years) was notable for prominent cortical sulci, symmetric ventriculo- megaly, subcortical hypomyelination, and a thin corpus callosum (Figure 2A). Skeletal radiographs showed meta- physeal dysplasia (most evident at hip joints) and both hy- poplasia and delayed ossification of epiphyses (Figure 2B). Two children had cervical radiographs that showed dens hypoplasia and, in one case, synostosis between the odon- toid and C2. Scoliosis was evident within the first few years of life, and coronal clefts could be observed at various levels of the vertebral column (Figure 2C). Two surviving LONP1 c.2161C>G homozygotes had an imperforate anus, in one case associated with omphalocele ican Journal of Human Genetics 96, 121–135, January 8, 2015 127 Table 2. Phenotype of CODAS-Syndrome-Affected Children Homozygous for LONP1 c.2161C>G (p.Arg721Gly) Agea 14.2b 4.6 1.2 0.6 Gender F M F M Cerebral Hypotonia and motor delay þ þ þ þ Intellectual disability þ þ n/a n/a Epilepsy þ � � � Ocular Ptosis þ þ þ þ Cataractsc þ þ þ þ Dental Delayed tooth eruption þ þ n/a n/a Enamel dysplasia þ þ n/a n/a Auricular Scapha and helix dysplasia (‘‘Crumpled’’ Ears) þ þ þ þ Conductive hearing loss n/a þ þ þ Sensorineural hearing loss þ þ � � Skeletal Short stature þ þ þ þ Metaphyseal dysplasia þ þ þ þ Epiphyseal hypoplasia and delayed ossification þ þ þ þ Scoliosis þ þ þ þ Genu valgus þ þ þ þ Pes valgus þ þ n/a n/a Vertebral coronal clefts þ þ n/a n/a Dens hypoplasia þ þ n/a n/a Other Grooved nasal tip þ þ þ þ Anteverted nares þ þ þ þ Tracheostomy þ þ þ þ Vocal-cord paresis � þ þ þ Endocardial cushion (atrial septum defect)d � þ þ � Gastresophageal reflux � þ � þ Gastrostomy tube feeding � þ þ þ Omphalocele � � � þ Imperforate anus þ � � þ Rectovaginal fistula þ � � � Cryptorchidism n/a þ n/a � Tongue hemiatrophy � þ � þ Abbreviations are as follows: F, female; M, male; n/a, not applicable because age or clinical data were not available. aLimited clinical data are available for four of eight children who had CODAS syndrome and died in the perinatal period as a result of respiratory 128 The American Journal of Human Genetics 96, 121–135, January 8 and in the other with a rectovaginal fistula. Omphalocele was also observed in two CODAS-affected siblings who died during infancy. One of two affected males had unilat- eral cryptorchidism. Two affected children had striking hemiatrophy of the tongue, presumably because of hypo- plasia or aplasia of the ipsilateral hypoglossal nerve; both also had vocal cord paresis and chronic sialorrhea. Given the diverse roles of Lon in mitochondrial homeo- stasis, we screened affected individuals for biochemical indices of mitochondrial dysfunction. Plasma lactate (1.4 5 0.9 mM; reference range 1.9 5 0.7 mM) and alanine (332 5 150 mM; reference range 435 5 121 mM)48 were normal in children with CODAS syndrome, as were all urine-tricarboxylic-acid-cycle intermediates, including citrate, alpha-ketoglutarate, succinate, fumarate, and ma- late49 (data not shown). Prenatal Presentation of CODAS Syndrome One unborn female Amish infant with CODAS syndrome had a prenatal ultrasound at 32 weeks’ gestational age. Esti- mated weight (1,516 g; 1st centile) and femur length (6.0 cm, 10th centile) were low; estimated head circumfer- ence was normal (32.1 cm; 93rd centile). The study was notable for the presence of polyhydramnios, a two-vessel umbilical cord, midfacial hypoplasia, crumpled helices, retrognathia, a balanced atrioventricular canal, omphalo- cele (containing intestine, stomach, and part of the liver), and mild thoracic scoliosis and the absence of left lower long bones, which were replaced by a relatively well- formed foot attached directly to the hip (Figures 2D–2G). She was stillborn at 39 weeks’ gestation. Location of Amino Acids Altered in CODAS Syndrome within the Homology Model of Lon The Lon holoenzyme is composed of six identical subunits that each consist of four domains: an MTS, a substrate recognition and binding (N) domain, an ATPase (AAAþ) domain, and a proteolytic (P) domain (Figure 3A). The amino terminal MTS targets the newly synthesized cytosolic polypeptide to the mitochondrial matrix, where removal of the MTS produces mature Lon. The AAAþ domain contains Walker boxes A and B that mediate ATP binding and hydrolysis, and the P domain active site contains a Ser855-Lys896 dyad required for peptide-bond hydrolysis.50 Using the crystal structure of B. subtilis Lon as a tem- plate, we established a homology model of human mito- chondrial Lon, which has a homohexameric ring-shaped arrangement.9 All four Lon variants identified in CODAS complications (these children are not represented in Table 1). Two had ompha- locele. bDied at 14.2 years of age as a result of accidental trauma. cDense nuclear cataracts develop rapidly between 2 and 6 months of age. dTwo surviving children with CODAS syndrome have atrial septal defects; two who died perinatally had atrioventricular canal defects. The total incidence of congenital heart disease was 50%. , 2015 syndrome resulted in amino acid substitutions within the AAAþ domain. Ser631, Arg721, and Ala724 are located at the interface between Lon subunits and mediate subunit interactions. Pro676 is at the ATP-binding site (Figures 3A and 3B). Replacements of Pro676 with serine (p.Pro676Ser) and Arg721with glycine (p.Arg721Gly) are pathogenic in the homozygous state (Table 1), allowing homo-oligomeric composition of the holoenzyme to be inferred. Pro676 is on an alpha-helix near the ATP-binding pocket and is involved in subdomain interactions (Figure 3C); it gener- ates an alpha-helix kink that facilitates hydrophobic inter- actions between Leu667, Ala670, Leu696, and Val716. These interactions most likely maintain the geometry of the ATP-binding pocket. Arg721 contributes to intersubu- nit interactions by forming a salt bridge with Glu654 on the neighboring subunit (Figure 3D). Arg721 is near Glu717, which also forms a salt bridge with Lys517 in the AAAþ domain of the neighboring subunit (Figure 3D). The p.Arg721Gly substitution most likely weakens intersu- bunit interactions and reduces stability of the holoenzyme. Recombinant CODAS Proteins Show Substrate- Specific Defects in Enzymatic Activity Individual CODAS variants were purified as recombinant proteins from bacteria, and their respective peptidase, pro- tease, and ATPase activities were compared with those of wild-type Lon. We measured ATP-dependent peptidase activity by using the fluorogenic decapeptide reporter substrate S3.51 Relative to wild-type protein, all four CODAS variants (p.Ser631Tyr, p.Pro676Ser, p.Arg721Gly, and p.Ala724Val) showed decreased energy-dependent peptidase activity; cleavage of S3 ranged from 19%–39% of wild-type activity (Figure 4A). For the Amish p.Arg721Gly variant, we also determined ATP-independent peptidase activity by using the fluores- cent dipeptide reporter AA2-Rh110, which is cleaved by wild-type Lon.39 In the absence of ATP, p.Arg721Gly and wild-type Lon showed significantly increased cleavage of AA2-Rh110. However, in the presence of ATP, peptidase ac- tivity of wild-type Lon was augmented 4-fold in compari- son to minimal augmentation of p.Arg721Gly, and direct determination of ATPase function showed that p.Arg721- Gly was 80% less active than wild-type Lon (data not shown). Finally, we examined ATP-dependent proteolysis of p.Pro676Ser and p.Arg721Gly by using two purified recom- binant natural Lon substrates: StAR19,44,45 and TFAM41,52 (Figures 4B and S1). Compared to wild-type Lon, which degraded 50%–60% StAR, p.Pro676Ser and p.Arg721Gly degraded only 10%–20% StAR during a 60 min incubation. In contrast, TFAM degradation by both p.Pro676Ser and p.Arg721Gly was comparable to that of the wild-type (50%–55% of TFAM was degraded during 60 min incubation)(Figure 4B). Control experiments confirmed that TFAM degradation was dependent on both ATP and Lon (Figures S1C and S1D). The efficiency with which The Amer p.Pro676Ser and p.Arg721Gly degrade different proteins is probably influenced by the substrate’s secondary or ter- tiary structure (i.e., whether the substrate is tightly or loosely folded). For example, folded and unfolded TFAM have been shown to be in equilibrium with one another.53 These results suggest that in CODAS syndrome, Lon-mediated degradation of some but not all protein sub- strates is impaired. Impaired Homo-oligomerization of Transiently Expressed R721G To characterize the behavior of CODAS variants in cultured cells, we engineered fusion proteins with C-terminal epitope tags to distinguish them from endogenous Lon. When transiently expressed in ARPE-19, UB/OC-2, or HeLa cells, CODAS and wild-type Lon were expressed at similar levels (Figures 4, S2, and S3). All four CODAS vari- ants localized to mitochondria, as determined by double labeling with immunofluorescence and Mitotracker Red (Figures 4, S2, and S3). To determine whether p.Arg721Gly was able to oligo- merize with wild-type or other p.Arg721Gly Lon, we per- formed co-immunoprecipitation experiments by using recombinant proteins with different epitope tags (e.g., FLAG and V5 epitopes) (Figure 4E), and we transiently co-expressed p.Arg721GlyFLAG or wild-typeFLAG with either p.Arg721GlyV5 or wild-typeV5 in HEK293T cells. Both wild- typeFLAG and p.Arg721GlyFLAG co-precipitated wild-typeV5; p.Arg721GlyFLAG co-precipitated significantly less wild- typeV5 (Figure 4E). Co-immunoprecipitation of p.Arg721- GlyV5 by p.Arg721GlyFLAG was barely detected, suggesting impaired homo-oligomeric assembly of p.Arg721Gly or disassembly of an unstable complex during the experi- mental procedure. Normal mtDNA Copy Number and Lon Substrate Levels in Immortalized CODAS Lymphoblast Cell Lines We generated EBV-immortalized B- LCLs from the periph- eral blood of two Amish family trios, each consisting of a homozygous p.Arg721Gly/p.Arg721Gly proband and his or her unaffected heterozygous p.Arg721Gly/wt parents. No difference in LONP1 transcript abundance was observed by RT-PCR in proband versus parent cells (Figure S4). Immunoblot analysis of whole-cell extracts demonstrated that Lon protein abundance was virtually identical in cell lines from the proband, mother, and father (Figure 5A). There was no compensatory upregula- tion of ClpX or ClpP, subunits of another mitochondrial AAAþ protease, and we found no change in the abundance of mitochondrial aconitase or TFAM (Figure 5A), which are natural Lon substrates under conditions of oxidative stress or mtDNA depletion, respectively.13,41 Recent work sug- gests that Lon is involved in degrading the PTEN-induced putative kinase 1 (PINK1).54,55 However, an increase in PINK1 abundance was not detected in proband or parent cells, not even in the presence of proteasome inhibitors ican Journal of Human Genetics 96, 121–135, January 8, 2015 129 Figure 5. Lon Protein Levels and mtDNA Copy Numbers in CO- DAS-Syndrome-Affected Probands Are Similar to Those of Their Unaffected Heterozygous Parents (A) Proteins were extracted from LCLs generated from CODAS- syndrome-affected probands and their parents. Triton X-100 lysis buffer was used, and 40 mg of the extract was immunoblotted for the proteins as shown. (B) Quantitative PCR (qPCR) was performed so that the relative mtDNA copy number in CODAS-syndrome-affected probands and their respective mothers and fathers could be determined with genomic DNA isolated from the peripheral blood lympho- cytes from the family members as shown. qPCR was performed with Taqman probes and primers for the mtDNA-encoded MT- CYB gene and the nuclear DNA-encoded APP gene as a reference. The relative quantity of mtDNA copy number was calculated by the DDCt method. Error bars represent the maximum and mini- mum relative quantitation. MG262 or epoxomicin (data not shown), which stabilize PINK1 in vitro.56,57 In lower and higher eukaryotes, knocking out Lon leads to reduced mtDNA copy-number or mtDNA dele- tions.22,23,58,59 We performed mtDNA copy-number anal- ysis and mtDNA sequencing by using total genomic DNA isolated from peripheral blood lymphocytes obtained from three homozygous p.Arg721Gly probands and their unaffected heterozygous parents. Quantitative PCR showed no consistent difference of mtDNA copy number between probands and parents (Figure 5B). Similar results were obtained with genomic DNA isolated from proband and parental LCLs (Figure S5). Using genomic DNA iso- lated from peripheral blood lymphocytes for deep sequencing of mtDNA exposed no deletions or significant differences between mtDNA sequences of each proband and that individual’s mother (Table S1). Abnormal Mitochondrial Morphology in CODAS Lymphoblast Cell Lines Transmission electron microscopy of proband LCLs re- vealed enlarged mitochondria with swollen intra- or inter- cristal compartments, uniform vesicular structures, and electron-dense intramitochondrial inclusions (Figures 6A and 6C), suggestive of abnormal inner-membrane topol- ogy.60,61 By contrast, parental LCLs showed typical mito- chondrial morphology with well-organized cristae and only occasional swollen cristal spaces (Figures 6B and 6D, and S6). Approximately 43 5 3.4% (SEM) of mitochondria (n ¼ 59) examined from two CODAS probands displayed 130 The American Journal of Human Genetics 96, 121–135, January 8 abnormal mitochondrial morphology (p < 0.001), whereas only 5.2 5 3.7% of paternal (n ¼ 48), 14.1 5 3.3% of maternal (n ¼ 63), and 10.6 5 5.1% (n ¼ 326) of homozy- gous-normal LCL mitochondria were abnormal. There was no significant difference between maternal and paternal mitochondria (p ¼ 0.232). Reduced MT-CO2 and SRC in CODAS Mitochondria Lon has been shown to mediate protein quality control and remodeling of cytochrome c oxidase, which is com- plex IV of the mitochondrial respiratory chain.16,17,21 We examined the abundance of various protein subunits belonging to respiratory chain complexes I-V and detected a specific and substantial reduction of mtDNA-encoded subunit II of cytochrome c oxidase (MT-CO2) in proband cells (Figure 7A). In contrast, proband and parental LCLs had similar abundances of other mtDNA-encoded (MT- CYB, MT-CO1) and nuclear-encoded (NDUFB8, SDHB, UQCRC2, ATP5A1) respiratory chain subunits. One possible explanation for reduced MT-CO2 abun- dance in CODAS cell lines was its poor solubility in the protein extraction buffer employed, which contained 0.5% Triton X-100. To address this, we used cell-lysis buffer with 8 M urea, which effectively solubilizes most aggre- gated proteins and inclusion bodies. When cells were ex- tracted with this strong denaturant, similar abundances of MT-CO2 were observed in proband and parental LCLs. Urea did not alter solubility of other mtDNA- or nuclear- encoded respiratory-chain subunits (Figure 7A). Insolubility of MT-CO2 might be caused by the failure of p.Arg721Gly to properly degrade and/or assemble this sub- unit into complex IV. To determine whether MT-CO2 is a natural substrate of Lon, we knocked down the protease in different cell lines. MT-CO2 is highly hydrophobic and difficult to purify; we therefore used a cellular system rather than directly testing purified Lon activity against that of MT-CO2. Inducible or transient knockdown of Lon was associated with increased MT-CO2 abundance in HeLa, HEK293T, and PC3 cells (Figure 7B). Taken together, these results suggest that MT-CO2 is a natural Lon sub- strate, and that p.Arg721Gly’s failure either to degrade or to properly assemble this protein in CODAS cells results in its aggregation within mitochondria. After observing the reduced solubility of MT-CO2, we investigated the efficiency of mitochondrial respiration. Normalized basal mitochondrial oxygen consumption and ATP-linked respiration were similar in proband and parental cells, but CODAS cells had significantly lower SRC when mitochondrial membrane potential was dissi- pated by the uncoupler FCCP (Figure 7C). Discussion Pathophysiology of CODAS Syndrome Mitochondrial Lon is a multi-functional enzyme with diverse roles that include (1) elimination of misfolded and oxidatively damaged proteins,12–14,20,62 (2) selective , 2015 Figure 6. Abnormal Mitochondrial Ultra- structure in CODAS Cells (A–D) Transmission electron micrographs were generated from cell lines of (A and C) two CODAS-syndrome-affected chil- dren homozygous for p.Arg721Gly (Gly721/Gly721) and (B and D, respectively) their heterozygous unaffected mothers (Arg721/Gly721). (A and C) In CODAS cells, 43% of mitochondria (n ¼ 59 cells, p < 0.001) display abnormal morphology, including swollen intracristal space and/or vesiculated cristae and electron-dense ag- gregations (white arrows). The mitochon- drion below the asterisk in (A0) appears to be in transition toward the vesiculated state when it is compared to the normal mitochondrion to the left of the asterisk. (B and D) The mothers’ cells display normal mitochondrial morphology (e.g., arrowheads); 14% show swelling of the intracristal space (n ¼ 63 cells, p < 0.001) (N ¼ nucleus). Scale bars represent 1 mm in (A–D) and 0.3 mm in magnified regions (A0–D0). The fathers’ cells show normal mitochondrial morphology, and only 5.2% show cristal-space swelling (n ¼ 48) (Figure S6). degradation of key rate-determining mitochondrial pro- teins in response to metabolic flux and cell stress,16– 18,58,63 (3) chaperone-like assembly of protein complexes within the respiratory chain,21,64 and (4) regulation of mitochondrial gene expression.46 Homozygous loss of Lon in mice results in embryonic death associated with poor growth and reduced mtDNA copy number, indicating the vital importance of this protein in mammals.15 In contrast, heterozygous mice mature similarly to wild- type.15 These findings, combined with our enzymatic studies, support the idea that allele combinations giving rise to CODAS syndrome must result in alteration and not complete abrogation of Lon function. We speculate that the complex pathophysiology of CODAS syndrome most likely reflects dysregulation of multiple Lon-dependent proteins and pathways. There are clinical similarities between CODAS syndrome and other mitochondrial disorders, including hypotonia, motor delay, ptosis, and sensorineural hearing loss.48,65–67 However, many classic clinical and biochemical hall- marks48,67,68 of mitochondrial dysfunction were not observed in CODAS children, and several CODAS manifes- tations, including postnatal cataracts, skeletal and dental abnormalities, and conductive hearing loss, are atypical of mitochondrial disease. Such complications might reflect Lon-dependent processes that are currently unknown. Moreover, tissue-specific defects could arise because of dif- ferential histological dependence on Lon as a protease, a The American Journal of Human G chaperone, and/or a DNA-binding protein. For example, some tissues might be more dependent on Lon to degrade oxidatively damaged pro- teins, whereas others might depend critically on as-yet- unknown Lon-mediated pathways. The reduced SRC of CODAS LCLs suggests limited capacity to boost ATP pro- duction in response to high energy demands, such as oc- curs during increased cellular stress and work load. Further investigations using tissues with metabolic demands higher than those of lymphoblastoid cells are needed if we are to delineate the physiological relevance of mecha- nistic defects in mitochondrial function and morphology associated with CODAS syndrome. Pathogenic Amino Acid Changes in CODAS Syndrome Cluster in the AAAþ Domain Our structural analysis of homozygous CODAS variants is consistent with observed alterations of their enzyme activ- ity in vitro. Proline676 is located within an a helix near the ATP-binding pocket and generates a kink in this a helix69 (Figure 3C). Removal of the kink at Pro676 might alter key hydrophobic interactions with four other residues— Ala670 (on the same helix as Pro676), Leu667, Leu696, and Val716—that likely help to shape the active site (Figure 3C). The p.Pro676Ser substitution could disrupt this kink, changing active site geometry or flexibility of the AAAþ module (Figure 3C) and decreasing the protein’s capacity for ATP binding and hydrolysis. Accordingly, we found reduced ATP-dependent peptidase (Figure 4A) and substrate-specific protease (Figure 4B) activities of the p.Pro676Ser variant. enetics 96, 121–135, January 8, 2015 131 Figure 7. Alterations in the Structure and Function of the Electron Transport Chain in Proband LCLs (A) For extraction of cellular proteins, lysis buffer containing 0.1% Triton X-100 (TX100) or 8 M urea (UREA) was used. Immunoblotting was performed for proteins as shown. MT-CO2 abundances were selectively reduced in CODAS cells but were recovered in the presence of 8 M urea. (B) MT-CO2 is a protein substrate of mitochondrial Lon. A HeLa cell line with a stably integrated anhydrotetracycline (ATc)-inducible shRNA targeting Lon was grown in the presence or absence of ATc (500 ng/ml) for 4 days. HEK293T cells and PC3 cells were transiently transfected with a siRNA targeting the 30 untranslated regionof Lon or a control siRNA and cultured for 3 days. Cell extracts (25 mg) were immunoblotted for Lon, MT-CO2, or actin as the loading control. (C) Basal normalized oxygen-consumption rate was established at 18 min (representing 100%), followed by addition of oligomycin (Oligo, 1 mM); FCCP (750 nM); and finally rotenone (1 mM). A Student’s t test showed that cells from CODAS-syndrome-affected probands had significantly (p < 0.05) lower spare mitochondrial respiratory capacity in the presence of the uncoupler FCCP. Error bars represent 5SD. Arg721 also makes an important intersubunit contact by forming a salt bridge with Glu654 from the AAAþ module of the neighboring subunit (Figure 3D). This salt bridge, together with one formed between Glu717 and Lys517, appears to be critical for stabilizing the Lon complex. p.Arg721Gly might thus destabilize the Lon complex and explain reduced homo-oligomerization, ATPase activity, and peptidase activity of this variant (Figure 4). Functional Consequences of Homozygous Lon p.Arg721Gly in Human Cells Our studies using recombinant protein and CODAS LCLs indicate that the p.Arg721Gly pathogenic variant has altered activity. Purified p.Arg721Gly has reduced ATPase and peptidase activities and degrades some natural sub- strates (e.g., TFAM) but not others (e.g., StAR) (Figure 4B). Substrate-restricted protease activity might reflect poor conformational resilience or stability of the p.Arg721Gly holoenzyme, as suggested by co-immunoprecipitation experiments (Figure 4E). In this study, we show that MT-CO2 is upregulated when Lon is knocked down in three different cell types, suggest- ing that MT-CO2 is a Lon substrate. In CODAS cells, we observed insolubility and aggregation of MT-CO2 (Figure 7A) but not MT-CO1, which is similarly hydropho- bic. Previous work shows that wild-type Lon plays a role in the assembly of the cytochrome c oxidase complex and in its remodeling the during adaptive switching between normoxic and hypoxic conditions.16,17,21 Work in yeast and mammalian cells suggests that Lon promotes cyto- chrome c oxidase assembly through a chaperone function independent of its protease activity.16,21 132 The American Journal of Human Genetics 96, 121–135, January 8 The failure of p.Arg721Gly to degrade or chaperone MT-CO2 (and possibly other inner-membrane proteins yet to be identified) might promote protein aggregation and explain the morphological and bioenergetic changes observed in CODAS cell mitochondria (Figures 6 and 7). In this study, we examined immortalized B-LCLs generated from CODAS-syndrome-affected probands and their par- ents. Because this cell type is not principally affected in CODAS syndrome, future studies are aimed at employing primary cells and tissue samples as well as cell lines engi- neered with specific CODAS LONP1 mutations. Natural History and Treatment of CODAS Syndrome The high allele frequency (MAF ¼ 5.9%, carrier fre- quency ¼ 11.8%) of LONP1 c.2161C>G among the Amish was unexpected. Of all the known pathogenic alleles in the Lancaster County Amish population, the LONP1 variant has the second-highest carrier frequency, eclipsed only by that of a founder mutation for Ellis-van Creveld syndrome (EVC; MIM 225500; RefSeq NM_153717.2); c.1886þ5G>T; carrier frequency ¼ 16.2%). The eight known Amish CODAS-syndrome-affected individuals comprise six sibships. In two sibships, three affected children died shortly after birth. One affected fetus was stillborn. Given the high allele frequency and perinatal mortality, the Amish CODAS variant might result in a high rate of fetal demise. This could explain an apparent discrepancy between the relatively low birth incidence of CODAS syndrome and the high allele frequency of LONP1 c.2161C>G in the Amish population. Accurate prenatal diagnosis of CODAS syndrome has important clinical implications (Figure 2). Newborns , 2015 with CODAS syndrome are likely to struggle through the neonatal transition. For syndrome variants (e.g., p.Arg721- Gly) associated with neonatal vocal-cord paresis, the ability to anticipate and prevent perinatal asphyxiation can be life saving. Three of eight (38%) Amish CODAS- syndrome-affected individuals died of airway complica- tions shortly after birth, and one was stillborn. Surviving children required intubation, mechanical respiratory support, and tracheostomy. Incomplete cardiac septation occurred in 50% of p.Arg721Gly homozygotes and mani- fested in two cases (25%) as a complete atrioventricular defect. Thus, supportive cardiopulmonary care for babies with a suspected or confirmed diagnosis of CODAS syndrome requires tertiary-level subspecialty services and must begin immediately after birth. A careful, staged approach to subsequent management can allow for rela- tively good cognitive outcome and for independence in affected individuals. Supplemental Data Supplemental Data include six figures and two tables and can be found with this article online at http://dx.doi.org/10.1016/ j.ajhg.2014.12.003. Acknowledgments We thank Kimberly Jacob, Prashant Patel, Ana Quiroz, and Jennie Silver for technical assistance and D. Streck from the Molecular Genetics Laboratory at the Rutgers New Jersey Medical School for next-generation mtDNA sequencing. R.N.J., E.G.P., and K.A.S. were supported by Howard Hughes Medical Institute undergradu- ate science education awards 52006294 and 52007538. R.N.J. was also supported by the Center for Research on Women and Newborn Health and by ConnectCare3. I.L. was supported by National Science Foundation grant CHE-1213175. C.K.S. was supported in part by National Institutes of Health grants R01GM084039 and 1R21NS067668 and also by the NJ Health Science Foundation. The Clinic for Special Children is supported by charitable donations from the communities it serves. Received: November 16, 2014 Accepted: December 5, 2014 Published: January 8, 2015 Web Resources The URLs for data presented herein are as follows: I-Tasser Online Protein Structure and Function Predictions, Zhang Laboratory, University of Michigan, http://zhanglab.ccmb.med. umich.edu/I-TASSER/ Modeller, https://salilab.org/modeller/ Online Mendelian Inheritance in Man (OMIM), http://www. omim.org RefSeq: NCBI Reference Sequence Database, http://www.ncbi.nlm. nih.gov/refseq/ Sali Laboratory, University of California San Francisco, http:// salilab.org/index.html SAMtools, http://samtools.sourceforge.net/ The Amer References 1. Shebib, S.M., Reed, M.H., Shuckett, E.P., Cross, H.G., Perry, J.B., and Chudley, A.E. (1991). Newly recognized syndrome of cerebral, ocular, dental, auricular, skeletal anomalies: CODAS syndrome—a case report. Am. J. Med. Genet. 40, 88–93. 2. de Almeida, J.C., Vargas, F.R., Barbosa-Neto, J.G., and Llerena, J.C., Jr. (1995). CODAS syndrome: a new distinct MCA/MR syndrome with radiological changes of spondyloe- piphyseal dysplasia. Another case report. Am. J. Med. Genet. 55, 19–20. 3. Innes, A.M., Chudley, A.E., Reed, M.H., Shuckett, E.P., Hildes- Ripstein, G.E., and Greenberg, C.R. (2001). Third case of cerebral, ocular, dental, auricular, skeletal anomalies (CODAS) syndrome, further delineating a new malformation syndrome: first report of an affected male and review of literature. Am. J. Med. Genet. 102, 44–47. 4. Marlin, S., Ducou Le Pointe, H., Le Merrer, M., Portnoi, M.F., Chantot, S., Jonard, L., Mantel-Guiochon, A., Siffroi, J.P., Garabedian, E.N., and Denoyelle, F. (2010). Fourth case of cerebral, ocular, dental, auricular, skeletal syndrome (CODAS), description of new features and molecular analysis. Am. J. Med. Genet. A. 152A, 1510–1514. 5. Lee, I., and Suzuki, C.K. (2008). Functional mechanics of the ATP-dependent Lon protease- lessons from endogenous pro- tein and synthetic peptide substrates. Biochim. Biophys. Acta 1784, 727–735. 6. Hanson, P.I., and Whiteheart, S.W. (2005). AAAþ proteins: have engine, will work. Nat. Rev. Mol. Cell Biol. 6, 519–529. 7. Neuwald, A.F., Aravind, L., Spouge, J.L., and Koonin, E.V. (1999). AAAþ: A class of chaperone-like ATPases associated with the assembly, operation, and disassembly of protein com- plexes. Genome Res. 9, 27–43. 8. Truscott, K.N., Lowth, B.R., Strack, P.R., and Dougan, D.A. (2010). Diverse functions of mitochondrial AAAþ proteins: protein activation, disaggregation, and degradation. Biochem. Cell Biol. 88, 97–108. 9. Venkatesh, S., Lee, J., Singh, K., Lee, I., and Suzuki, C.K. (2012). Multitasking in the mitochondrion by the ATP-depen- dent Lon protease. Biochim. Biophys. Acta 1823, 56–66. 10. Jenkinson, E.M., Rehman, A.U., Walsh, T., Clayton-Smith, J., Lee, K., Morell, R.J., Drummond, M.C., Khan, S.N., Naeem, M.A., Rauf, B., et al.; University of Washington Center for Mendelian Genomics (2013). Perrault syndrome is caused by recessive mutations in CLPP, encoding a mitochondrial ATP- dependent chambered protease. Am. J. Hum. Genet. 92, 605–613. 11. Casari, G., De Fusco, M., Ciarmatori, S., Zeviani, M., Mora, M., Fernandez, P., De Michele, G., Filla, A., Cocozza, S., Marconi, R., et al. (1998). Spastic paraplegia and OXPHOS impairment caused by mutations in paraplegin, a nuclear-encoded mito- chondrial metalloprotease. Cell 93, 973–983. 12. Bota, D.A., Van Remmen, H., and Davies, K.J. (2002). Modula- tion of Lon protease activity and aconitase turnover during aging and oxidative stress. FEBS Lett. 532, 103–106. 13. Bota, D.A., and Davies, K.J. (2002). Lon protease preferentially degrades oxidized mitochondrial aconitase by an ATP-stimu- lated mechanism. Nat. Cell Biol. 4, 674–680. 14. Ngo, J.K., Pomatto, L.C., and Davies, K.J. (2013). Upregulation of the mitochondrial Lon Protease allows adaptation to acute oxidative stress but dysregulation is associated with chronic stress, disease, and aging. Redox Biol. 1, 258–264. ican Journal of Human Genetics 96, 121–135, January 8, 2015 133 http://dx.doi.org/10.1016/j.ajhg.2014.12.003 http://dx.doi.org/10.1016/j.ajhg.2014.12.003 http://zhanglab.ccmb.med.umich.edu/I-TASSER/ http://zhanglab.ccmb.med.umich.edu/I-TASSER/ https://salilab.org/modeller/ http://www.omim.org http://www.omim.org http://www.ncbi.nlm.nih.gov/refseq/ http://www.ncbi.nlm.nih.gov/refseq/ http://salilab.org/index.html http://salilab.org/index.html http://samtools.sourceforge.net/ 15. Quiros, P.M., Espanol, Y., Acin-Perez, R., Rodriguez, F., Barcena, C., Watanabe, K., Calvo, E., Loureiro, M., Fernan- dez-Garcia, M.S., Fueyo, A., et al. (2014). ATP-dependent Lon protease controls tumor bioenergetics by reprogramming mitochondrial activity. Cell Rep. 8, 542–556. 16. Hori, O., Ichinoda, F., Tamatani, T., Yamaguchi, A., Sato, N., Ozawa, K., Kitao, Y., Miyazaki, M., Harding, H.P., Ron, D., et al. (2002). Transmission of cell stress from endoplasmic reticulum to mitochondria: enhanced expression of Lon pro- tease. J. Cell Biol. 157, 1151–1160. 17. Fukuda, R., Zhang, H., Kim, J.W., Shimoda, L., Dang, C.V., and Semenza, G.L. (2007). HIF-1 regulates cytochrome oxidase subunits to optimize efficiency of respiration in hypoxic cells. Cell 129, 111–122. 18. Tian, Q., Li, T., Hou, W., Zheng, J., Schrum, L.W., and Bonkovsky, H.L. (2011). Lon peptidase 1 (LONP1)-dependent breakdown of mitochondrial 5-aminolevulinic acid synthase protein by heme in human liver cells. J. Biol. Chem. 286, 26424–26430. 19. Granot, Z., Kobiler, O., Melamed-Book, N., Eimerl, S., Bahat, A., Lu, B., Braun, S., Maurizi, M.R., Suzuki, C.K., Oppenheim, A.B., and Orly, J. (2007). Turnover of mitochondrial steroido- genic acute regulatory (StAR) protein by Lon protease: the un- expected effect of proteasome inhibitors. Mol. Endocrinol. 21, 2164–2177. 20. Wohlever, M.L., Baker, T.A., and Sauer, R.T. (2014). Roles of the N domain of the AAAþ Lon protease in substrate recognition, allosteric regulation and chaperone activity. Mol. Microbiol. 91, 66–78. 21. Rep, M., van Dijl, J.M., Suda, K., Schatz, G., Grivell, L.A., and Suzuki, C.K. (1996). Promotion of mitochondrial membrane complex assembly by a proteolytically inactive yeast Lon. Science 274, 103–106. 22. Suzuki, C.K., Suda, K., Wang, N., and Schatz, G. (1994). Requirement for the yeast gene LON in intramitochondrial proteolysis and maintenance of respiration. Science 264, 273–276. 23. Van Dyck, L., Pearce, D.A., and Sherman, F. (1994). PIM1 encodes a mitochondrial ATP-dependent protease that is required for mitochondrial function in the yeast Saccharo- myces cerevisiae. J. Biol. Chem. 269, 238–242. 24. Lu, B., Liu, T., Crosby, J.A., Thomas-Wohlever, J., Lee, I., and Suzuki, C.K. (2003). The ATP-dependent Lon protease of Mus musculus is a DNA-binding protein that is functionally conserved between yeast and mammals. Gene 306, 45–55. 25. Fu, G.K., and Markovitz, D.M. (1998). The human LON prote- ase binds to mitochondrial promoters in a single-stranded, site- specific, strand-specific manner. Biochemistry 37, 1905–1909. 26. Liu, T., Lu, B., Lee, I., Ondrovicová, G., Kutejová, E., and Suzuki, C.K. (2004). DNA and RNA binding by the mitochon- drial lon protease is regulated by nucleotide and protein sub- strate. J. Biol. Chem. 279, 13902–13910. 27. Chen, S.H., Suzuki, C.K., and Wu, S.H. (2008). Thermody- namic characterization of specific interactions between the human Lon protease and G-quartet DNA. Nucleic Acids Res. 36, 1273–1287. 28. Lohr, N.J., Molleston, J.P., Strauss, K.A., Torres-Martinez, W., Sherman, E.A., Squires, R.H., Rider, N.L., Chikwava, K.R., Cummings, O.W., Morton, D.H., and Puffenberger, E.G. (2010). Human ITCH E3 ubiquitin ligase deficiency causes syndromic multisystem autoimmune disease. Am. J. Hum. Genet. 86, 447–453. 134 The American Journal of Human Genetics 96, 121–135, January 8 29. Puffenberger, E.G., Hu-Lince, D., Parod, J.M., Craig, D.W., Dobrin, S.E., Conway, A.R., Donarum, E.A., Strauss, K.A., Dunckley, T., Cardenas, J.F., et al. (2004). Mapping of sudden infant death with dysgenesis of the testes syndrome (SIDDT) by a SNP genome scan and identification of TSPYL loss of function. Proc. Natl. Acad. Sci. USA 101, 11689–11694. 30. Puffenberger, E.G., Strauss, K.A., Ramsey, K.E., Craig, D.W., Stephan, D.A., Robinson, D.L., Hendrickson, C.L., Gottlieb, S., Ramsay, D.A., Siu, V.M., et al. (2007). Polyhydramnios, megalencephaly and symptomatic epilepsy caused by a homozygous 7-kilobase deletion in LYK5. Brain 130, 1929–1941. 31. Sherman, E.A., Strauss, K.A., Tortorelli, S., Bennett, M.J., Knerr, I., Morton, D.H., and Puffenberger, E.G. (2008). Genetic mapping of glutaric aciduria, type 3, to chromosome 7 and identification of mutations in c7orf10. Am. J. Hum. Genet. 83, 604–609. 32. Strauss, K.A., Puffenberger, E.G., Bunin, N., Rider, N.L., Morton, M.C., Eastman, J.T., 3rd, and Morton, D.H. (2008). Clinical application of DNA microarrays: molecular diagnosis and HLA matching of an Amish child with severe combined immune deficiency. Clin. Immunol. 128, 31–38. 33. Strauss, K.A., Puffenberger, E.G., Craig, D.W., Panganiban, C.B., Lee, A.M., Hu-Lince, D., Stephan, D.A., and Morton, D.H. (2005). Genome-wide SNP arrays as a diagnostic tool: clinical description, genetic mapping, and molecular charac- terization of Salla disease in an Old Order Mennonite popula- tion. Am. J. Med. Genet. A. 138A, 262–267. 34. Gnirke, A., Melnikov, A., Maguire, J., Rogov, P., LeProust, E.M., Brockman, W., Fennell, T., Giannoukos, G., Fisher, S., Russ, C., et al. (2009). Solution hybrid selection with ultra-long oligo- nucleotides for massively parallel targeted sequencing. Nat. Biotechnol. 27, 182–189. 35. Puffenberger, E.G., Jinks, R.N., Sougnez, C., Cibulskis, K., Willert, R.A., Achilly, N.P., Cassidy, R.P., Fiorentini, C.J., Heiken, K.F., Lawrence, J.J., et al. (2012). Genetic mapping and exome sequencing identify variants associated with five novel diseases. PLoS ONE 7, e28936. 36. Puffenberger, E.G., Jinks, R.N., Wang, H., Xin, B., Fiorentini, C., Sherman, E.A., Degrazio, D., Shaw, C., Sougnez, C., Cibul- skis, K., et al. (2012). A homozygous missense mutation in HERC2 associated with global developmental delay and autism spectrum disorder. Hum. Mutat. 33, 1639–1646. 37. DePristo, M.A., Banks, E., Poplin, R., Garimella, K.V., Maguire, J.R., Hartl, C., Philippakis, A.A., del Angel, G., Rivas, M.A., Hanna, M., et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498. 38. Duman, R.E., and Löwe, J. (2010). Crystal structures of Bacillus subtilis Lon protease. J. Mol. Biol. 401, 653–670. 39. Bernstein, S.H., Venkatesh, S., Li, M., Lee, J., Lu, B., Hilchey, S.P., Morse, K.M., Metcalfe, H.M., Skalska, J., Andreeff, M., et al. (2012). The mitochondrial ATP-dependent Lon protease: a novel target in lymphoma death mediated by the synthetic triterpenoid CDDO and its derivatives. Blood 119, 3321–3329. 40. Lee, I., and Berdis, A.J. (2001). Adenosine triphosphate-depen- dent degradation of a fluorescent lambda N substrate mimic by Lon protease. Anal. Biochem. 291, 74–83. 41. Lu, B., Lee, J., Nie, X., Li, M., Morozov, Y.I., Venkatesh, S., Bogenhagen, D.F., Temiakov, D., and Suzuki, C.K. (2013). Phosphorylation of human TFAM in mitochondria impairs , 2015 DNA binding and promotes degradation by the AAAþ Lon protease. Mol. Cell 49, 121–132. 42. Bradford, M.M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. 43. Pace, C.N., Vajdos, F., Fee, L., Grimsley, G., and Gray, T. (1995). How to measure and predict the molar absorption coefficient of a protein. Protein Sci. 4, 2411–2423. 44. Ondrovicová, G., Liu, T., Singh, K., Tian, B., Li, H., Gakh, O., Perecko, D., Janata, J., Granot, Z., Orly, J., et al. (2005). Cleav- age site selection within a folded substrate by the ATP- dependent lon protease. J. Biol. Chem. 280, 25103–25110. 45. Fishovitz, J., Li, M., Frase, H., Hudak, J., Craig, S., Ko, K., Berdis, A.J., Suzuki, C.K., and Lee, I. (2011). Active-site- directed chemical tools for profiling mitochondrial Lon prote- ase. ACS Chem. Biol. 6, 781–788. 46. Lu, B., Yadav, S., Shah, P.G., Liu, T., Tian, B., Pukszta, S., Villaluna, N., Kutejová, E., Newlon, C.S., Santos, J.H., and Suzuki, C.K. (2007). Roles for the human ATP-dependent Lon protease in mitochondrial DNA maintenance. J. Biol. Chem. 282, 17363–17374. 47. di Sant’Agnese, P.A., and De Mesy-Jensen, K.L. (1984). Diag- nostic electron microscopy on reembedded (‘‘popped off’’) areas of large Spurr epoxy sections. Ultrastruct. Pathol. 6, 247–253. 48. DiMauro, S., Schon, E.A., Carelli, V., and Hirano, M. (2013). The clinical maze of mitochondrial neurology. Nat Rev Neurol 9, 429–444. 49. Barshop, B.A. (2004). Metabolomic approaches to mitochon- drial disease: correlation of urine organic acids. Mitochon- drion 4, 521–527. 50. Botos, I., Melnikov, E.E., Cherry, S., Tropea, J.E., Khalatova, A.G., Rasulova, F., Dauter, Z., Maurizi, M.R., Rotanova, T.V., Wlodawer, A., and Gustchina, A. (2004). The catalytic domain of Escherichia coli Lon protease has a unique fold and a Ser-Lys dyad in the active site. J. Biol. Chem. 279, 8140–8148. 51. Thomas-Wohlever, J., and Lee, I. (2002). Kinetic characteriza- tion of the peptidase activity of Escherichia coli Lon reveals the mechanistic similarities in ATP-dependent hydrolysis of peptide and protein substrates. Biochemistry 41, 9418–9425. 52. Matsushima, Y., Goto, Y., and Kaguni, L.S. (2010). Mitochon- drial Lon protease regulates mitochondrial DNA copy number and transcription by selective degradation of mitochondrial transcription factor A (TFAM). Proc. Natl. Acad. Sci. USA 107, 18410–18415. 53. Wong, T.S., Rajagopalan, S., Freund, S.M., Rutherford, T.J., Andreeva, A., Townsley, F.M., Petrovich, M., and Fersht, A.R. (2009). Biophysical characterizations of human mitochon- drial transcription factor A and its binding to tumor suppres- sor p53. Nucleic Acids Res. 37, 6765–6783. 54. Jin, S.M., and Youle, R.J. (2013). The accumulation of mis- folded proteins in the mitochondrial matrix is sensed by PINK1 to induce PARK2/Parkin-mediated mitophagy of polar- ized mitochondria. Autophagy 9, 1750–1757. 55. Thomas, R.E., Andrews, L.A., Burman, J.L., Lin, W.Y., and Pallanck, L.J. (2014). PINK1-Parkin pathway activity is regu- The Amer lated by degradation of PINK1 in the mitochondrial matrix. PLoS Genet. 10, e1004279. 56. Um, J.W., Stichel-Gunkel, C., Lübbert, H., Lee, G., and Chung, K.C. (2009). Molecular interaction between parkin and PINK1 in mammalian neuronal cells. Mol. Cell. Neurosci. 40, 421–432. 57. Lin, W., and Kang, U.J. (2008). Characterization of PINK1 pro- cessing, stability, and subcellular localization. J. Neurochem. 106, 464–474. 58. Quirós, P.M., Español, Y., Acı́n-Pérez, R., Rodrı́guez, F., Bárcena, C., Watanabe, K., Calvo, E., Loureiro, M., Fernán- dez-Garcı́a, M.S., Fueyo, A., et al. (2014). ATP-dependent Lon protease controls tumor bioenergetics by reprogramming mitochondrial activity. Cell Rep. 8, 542–556. 59. Guha, S., López-Maury, L., Shaw, M., Bähler, J., Norbury, C.J., and Agashe, V.R. (2011). Transcriptional and cellular re- sponses to defective mitochondrial proteolysis in fission yeast. J. Mol. Biol. 408, 222–237. 60. Mannella, C.A. (2006). Structure and dynamics of the mito- chondrial inner membrane cristae. Biochim. Biophys. Acta 1763, 542–548. 61. Scorrano, L., Ashiya, M., Buttle, K., Weiler, S., Oakes, S.A., Mannella, C.A., and Korsmeyer, S.J. (2002). A distinct pathway remodels mitochondrial cristae and mobilizes cytochrome c during apoptosis. Dev. Cell 2, 55–67. 62. Smakowska, E., Czarna, M., and Janska, H. (2014). Mitochon- drial ATP-dependent proteases in protection against accumu- lation of carbonylated proteins. Mitochondrion 19 (Pt B), 245–251. 63. Teng, H., Wu, B., Zhao, K., Yang, G., Wu, L., and Wang, R. (2013). Oxygen-sensitive mitochondrial accumulation of cys- tathionine b-synthase mediated by Lon protease. Proc. Natl. Acad. Sci. USA 110, 12679–12684. 64. Solheim, C., Li, L., Hatzopoulos, P., and Millar, A.H. (2012). Loss of Lon1 in Arabidopsis changes the mitochondrial prote- ome leading to altered metabolite profiles and growth retarda- tion without an accumulation of oxidative damage. Plant Physiol. 160, 1187–1203. 65. Di Fonzo, A., Ronchi, D., Lodi, T., Fassone, E., Tigano, M., Lamperti, C., Corti, S., Bordoni, A., Fortunato, F., Nizzardo, M., et al. (2009). The mitochondrial disulfide relay system pro- tein GFER is mutated in autosomal-recessive myopathy with cataract and combined respiratory-chain deficiency. Am. J. Hum. Genet. 84, 594–604. 66. Hirano, M., Nishigaki, Y., and Martı́, R. (2004). Mitochondrial neurogastrointestinal encephalomyopathy (MNGIE): a dis- ease of two genomes. Neurologist 10, 8–17. 67. Koopman, W.J., Distelmaier, F., Smeitink, J.A., and Willems, P.H. (2013). OXPHOS mutations and neurodegeneration. EMBO J. 32, 9–29. 68. Breuer, M.E., Koopman, W.J., Koene, S., Nooteboom, M., Rodenburg, R.J., Willems, P.H., and Smeitink, J.A. (2013). The role of mitochondrial OXPHOS dysfunction in the devel- opment of neurologic diseases. Neurobiol. Dis. 51, 27–34. 69. Wilman, H.R., Shi, J., and Deane, C.M. (2014). Helix kinks are equally prevalent in soluble and membrane proteins. Proteins 82, 1960–1970. ican Journal of Human Genetics 96, 121–135, January 8, 2015 135 CODAS Syndrome Is Associated with Mutations of LONP1, Encoding Mitochondrial AAA+ Lon Protease Introduction Subjects and Methods Subjects Molecular Genetics Molecular Modeling cDNA Constructs Purification of Recombinant Lon Enzymatic Assays Peptidase and Protease Assays ATPase Assay Cell Culture Mitochondrial Localization Immunoblotting and Immunoprecipitation Transmission Electron Microscopy Quantitative PCR (qPCR) Determination of mtDNA Copy Number Mitochondrial Oxygen-Consumption Measurements Results Molecular Genetic Studies Phenotype Associated with CODAS Syndrome Prenatal Presentation of CODAS Syndrome Location of Amino Acids Altered in CODAS Syndrome within the Homology Model of Lon Recombinant CODAS Proteins Show Substrate-Specific Defects in Enzymatic Activity Impaired Homo-oligomerization of Transiently Expressed R721G Normal mtDNA Copy Number and Lon Substrate Levels in Immortalized CODAS Lymphoblast Cell Lines Abnormal Mitochondrial Morphology in CODAS Lymphoblast Cell Lines Reduced MT-CO2 and SRC in CODAS Mitochondria Discussion Pathophysiology of CODAS Syndrome Pathogenic Amino Acid Changes in CODAS Syndrome Cluster in the AAA+ Domain Functional Consequences of Homozygous Lon p.Arg721Gly in Human Cells Natural History and Treatment of CODAS Syndrome Supplemental Data Acknowledgments Web Resources References