A complete pedigree-based graph workflow for rare candidate variant analysis.

Methods that use a linear genome reference for genome sequencing data analysis are reference-biased. In the field of clinical genetics for rare diseases, a resulting reduction in genotyping accuracy in some regions has likely prevented the resolution of some cases. Pangenome graphs embed population variation into a reference structure. Although pangenome graphs have helped to reduce reference mapping bias, further performance improvements are possible. We introduce VG-Pedigree, a pedigree-aware workflow based on the pangenome-mapping tool of Giraffe and the variant calling tool DeepTrio using a specially trained model for Giraffe-based alignments. We demonstrate mapping and variant calling improvements in both single-nucleotide variants (SNVs) and insertion and deletion (indel) variants over those produced by alignments created using BWA-MEM to a linear-reference and Giraffe mapping to a pangenome graph containing data from the 1000 Genomes Project. We have also adapted and upgraded deleterious-variant (DV) detecting methods and programs into a streamlined workflow. We used these workflows in combination to detect small lists of candidate DVs among 15 family quartets and quintets of the Undiagnosed Diseases Program (UDP). All candidate DVs that were previously diagnosed using the Mendelian models covered by the previously published methods were recapitulated by these workflows. The results of these experiments indicate that a slightly greater absolute count of DVs are detected in the proband population than in their matched unaffected siblings.

Novel CUL3 Variant Causing Familial Hyperkalemic Hypertension Impairs Regulation and Function of Ubiquitin Ligase Activity.

Familial hyperkalemic hypertension is caused by pathogenic variants in genes of the CUL3 (cullin-3)-KLHL3 (kelch-like-family-member-3)-WNK (with no-lysine [K] kinase) pathway, manifesting clinically as hyperkalemia, metabolic acidosis, and high systolic blood pressure. The ubiquitin E3 ligase CUL3-KLHL3 targets WNK kinases for degradation to limit activation of the thiazide-sensitive NCC (Na-Cl cotransporter). All known variants in CUL3 lead to exon 9 skipping (CUL3Δ9) and typically result in severe familial hyperkalemic hypertension and growth disturbances in patients. Whether other variants in CUL3 cause familial hyperkalemic hypertension is unknown. Here, we identify a novel de novo heterozygous CUL3 variant (CUL3Δ474-477) in a pediatric familial hyperkalemic hypertension patient with multiple congenital anomalies and reveal molecular mechanisms by which CUL3Δ474-477 leads to dysregulation of the CUL3-KLHL3-WNK signaling axis. Using patient-derived urinary extracellular vesicles and dermal fibroblasts, in vitro assays, and cultured kidney cells, we demonstrate that CUL3Δ474-477 causes reduced total CUL3 levels due to increased autoubiquitination. The CUL3Δ474-477 that escapes autodegradation shows enhanced modification with NEDD8 (neural precursor cell expressed developmentally down-regulated protein 8) and increased formation of CUL3-KLHL3 complexes that are impaired in ubiquitinating WNK4. Proteomic analysis of CUL3 complexes revealed that, in addition to increased KLHL3 binding, the CUL3Δ474-477 variant also exhibits increased interactions with other BTB (Bric-a-brac, Tramtrack, and Broad complex) substrate adaptors, providing a rationale for the patient's diverse phenotypes. We conclude that the pathophysiological effects of CUL3Δ474-477 are caused by reduced CUL3 levels and formation of catalytically impaired CUL3 ligase complexes.

FOXR1 regulates stress response pathways and is necessary for proper brain development.

The forkhead box (Fox) family of transcription factors are highly conserved and play essential roles in a wide range of cellular and developmental processes. We report an individual with severe neurological symptoms including postnatal microcephaly, progressive brain atrophy and global developmental delay associated with a de novo missense variant (M280L) in the FOXR1 gene. At the protein level, M280L impaired FOXR1 expression and induced a nuclear aggregate phenotype due to protein misfolding and proteolysis. RNAseq and pathway analysis showed that FOXR1 acts as a transcriptional activator and repressor with central roles in heat shock response, chaperone cofactor-dependent protein refolding and cellular response to stress pathways. Indeed, FOXR1 expression is increased in response to cellular stress, a process in which it directly controls HSPA6, HSPA1A and DHRS2 transcripts. The M280L mutant compromises FOXR1's ability to respond to stress, in part due to impaired regulation of downstream target genes that are involved in the stress response pathway. Quantitative PCR of mouse embryo tissues show Foxr1 expression in the embryonic brain. Using CRISPR/Cas9 gene editing, we found that deletion of mouse Foxr1 leads to a severe survival deficit while surviving newborn Foxr1 knockout mice have reduced body weight. Further examination of newborn Foxr1 knockout brains revealed a decrease in cortical thickness and enlarged ventricles compared to littermate wild-type mice, suggesting that loss of Foxr1 leads to atypical brain development. Combined, these results suggest FOXR1 plays a role in cellular stress response pathways and is necessary for normal brain development.

Compound heterozygous KCTD7 variants in progressive myoclonus epilepsy.

KCTD7 is a member of the potassium channel tetramerization domain-containing protein family and has been associated with progressive myoclonic epilepsy (PME), characterized by myoclonus, epilepsy, and neurological deterioration. Here we report four affected individuals from two unrelated families in which we identified KCTD7 compound heterozygous single nucleotide variants through exome sequencing. RNAseq was used to detect a non-annotated splicing junction created by a synonymous variant in the second family. Whole-cell patch-clamp analysis of neuroblastoma cells overexpressing the patients' variant alleles demonstrated aberrant potassium regulation. While all four patients experienced many of the common clinical features of PME, they also showed variable phenotypes not previously reported, including dysautonomia, brain pathology findings including a significantly reduced thalamus, and the lack of myoclonic seizures. To gain further insight into the pathogenesis of the disorder, zinc finger nucleases were used to generate kctd7 knockout zebrafish. Kctd7 homozygous mutants showed global dysregulation of gene expression and increased transcription of c-fos, which has previously been correlated with seizure activity in animal models. Together these findings expand the known phenotypic spectrum of KCTD7-associated PME, report a new animal model for future studies, and contribute valuable insights into the disease.

GGPS1 Mutations Cause Muscular Dystrophy/Hearing Loss/Ovarian Insufficiency Syndrome.

OBJECTIVE: A hitherto undescribed phenotype of early onset muscular dystrophy associated with sensorineural hearing loss and primary ovarian insufficiency was initially identified in 2 siblings and in subsequent patients with a similar constellation of findings. The goal of this study was to understand the genetic and molecular etiology of this condition. METHODS: We applied whole exome sequencing (WES) superimposed on shared haplotype regions to identify the initial biallelic variants in GGPS1 followed by GGPS1 Sanger sequencing or WES in 5 additional families with the same phenotype. Molecular modeling, biochemical analysis, laser membrane injury assay, and the generation of a Y259C knock-in mouse were done. RESULTS: A total of 11 patients in 6 families carrying 5 different biallelic pathogenic variants in specific domains of GGPS1 were identified. GGPS1 encodes geranylgeranyl diphosphate synthase in the mevalonate/isoprenoid pathway, which catalyzes the synthesis of geranylgeranyl pyrophosphate, the lipid precursor of geranylgeranylated proteins including small guanosine triphosphatases. In addition to proximal weakness, all but one patient presented with congenital sensorineural hearing loss, and all postpubertal females had primary ovarian insufficiency. Muscle histology was dystrophic, with ultrastructural evidence of autophagic material and large mitochondria in the most severe cases. There was delayed membrane healing after laser injury in patient-derived myogenic cells, and a knock-in mouse of one of the mutations (Y259C) resulted in prenatal lethality. INTERPRETATION: The identification of specific GGPS1 mutations defines the cause of a unique form of muscular dystrophy with hearing loss and ovarian insufficiency and points to a novel pathway for this clinical constellation. ANN NEUROL 2020;88:332-347.

Deficiency in the endocytic adaptor proteins PHETA1/2 impairs renal and craniofacial development.

A critical barrier in the treatment of endosomal and lysosomal diseases is the lack of understanding of the in vivo functions of the putative causative genes. We addressed this by investigating a key pair of endocytic adaptor proteins, PH domain-containing endocytic trafficking adaptor 1 and 2 (PHETA1/2; also known as FAM109A/B, Ses1/2, IPIP27A/B), which interact with the protein product of OCRL, the causative gene for Lowe syndrome. Here, we conducted the first study of PHETA1/2 in vivo, utilizing the zebrafish system. We found that impairment of both zebrafish orthologs, pheta1 and pheta2, disrupted endocytosis and ciliogenesis in renal tissues. In addition, pheta1/2 mutant animals exhibited reduced jaw size and delayed chondrocyte differentiation, indicating a role in craniofacial development. Deficiency of pheta1/2 resulted in dysregulation of cathepsin K, which led to an increased abundance of type II collagen in craniofacial cartilages, a marker of immature cartilage extracellular matrix. Cathepsin K inhibition rescued the craniofacial phenotypes in the pheta1/2 double mutants. The abnormal renal and craniofacial phenotypes in the pheta1/2 mutant animals were consistent with the clinical presentation of a patient with a de novo arginine (R) to cysteine (C) variant (R6C) of PHETA1. Expressing the patient-specific variant in zebrafish exacerbated craniofacial deficits, suggesting that the R6C allele acts in a dominant-negative manner. Together, these results provide insights into the in vivo roles of PHETA1/2 and suggest that the R6C variant is contributory to the pathogenesis of disease in the patient.This article has an associated First Person interview with the first author of the paper.

De Novo Variants in WDR37 Are Associated with Epilepsy, Colobomas, Dysmorphism, Developmental Delay, Intellectual Disability, and Cerebellar Hypoplasia.

WD40 repeat-containing proteins form a large family of proteins present in all eukaryotes. Here, we identified five pediatric probands with de novo variants in WDR37, which encodes a member of the WD40 repeat protein family. Two probands shared one variant and the others have variants in nearby amino acids outside the WD40 repeats. The probands exhibited shared phenotypes of epilepsy, colobomas, facial dysmorphology reminiscent of CHARGE syndrome, developmental delay and intellectual disability, and cerebellar hypoplasia. The WDR37 protein is highly conserved in vertebrate and invertebrate model organisms and is currently not associated with a human disease. We generated a null allele of the single Drosophila ortholog to gain functional insights and replaced the coding region of the fly gene CG12333/wdr37 with GAL4. These flies are homozygous viable but display severe bang sensitivity, a phenotype associated with seizures in flies. Additionally, the mutant flies fall when climbing the walls of the vials, suggesting a defect in grip strength, and repeat the cycle of climbing and falling. Similar to wall clinging defect, mutant males often lose grip of the female abdomen during copulation. These phenotypes are rescued by using the GAL4 in the CG12333/wdr37 locus to drive the UAS-human reference WDR37 cDNA. The two variants found in three human subjects failed to rescue these phenotypes, suggesting that these alleles severely affect the function of this protein. Taken together, our data suggest that variants in WDR37 underlie a novel syndromic neurological disorder.

Lysosomal Storage and Albinism Due to Effects of a De Novo CLCN7 Variant on Lysosomal Acidification.

Optimal lysosome function requires maintenance of an acidic pH maintained by proton pumps in combination with a counterion transporter such as the Cl-/H+ exchanger, CLCN7 (ClC-7), encoded by CLCN7. The role of ClC-7 in maintaining lysosomal pH has been controversial. In this paper, we performed clinical and genetic evaluations of two children of different ethnicities. Both children had delayed myelination and development, organomegaly, and hypopigmentation, but neither had osteopetrosis. Whole-exome and -genome sequencing revealed a de novo c.2144A>G variant in CLCN7 in both affected children. This p.Tyr715Cys variant, located in the C-terminal domain of ClC-7, resulted in increased outward currents when it was heterologously expressed in Xenopus oocytes. Fibroblasts from probands displayed a lysosomal pH approximately 0.2 units lower than that of control cells, and treatment with chloroquine normalized the pH. Primary fibroblasts from both probands also exhibited markedly enlarged intracellular vacuoles; this finding was recapitulated by the overexpression of human p.Tyr715Cys CLCN7 in control fibroblasts, reflecting the dominant, gain-of-function nature of the variant. A mouse harboring the knock-in Clcn7 variant exhibited hypopigmentation, hepatomegaly resulting from abnormal storage, and enlarged vacuoles in cultured fibroblasts. Our results show that p.Tyr715Cys is a gain-of-function CLCN7 variant associated with developmental delay, organomegaly, and hypopigmentation resulting from lysosomal hyperacidity, abnormal storage, and enlarged intracellular vacuoles. Our data supports the hypothesis that the ClC-7 antiporter plays a critical role in maintaining lysosomal pH.

Biallelic HEPHL1 variants impair ferroxidase activity and cause an abnormal hair phenotype.

Maintenance of the correct redox status of iron is functionally important for critical biological processes. Multicopper ferroxidases play an important role in oxidizing ferrous iron, released from the cells, into ferric iron, which is subsequently distributed by transferrin. Two well-characterized ferroxidases, ceruloplasmin (CP) and hephaestin (HEPH) facilitate this reaction in different tissues. Recently, a novel ferroxidase, Hephaestin like 1 (HEPHL1), also known as zyklopen, was identified. Here we report a child with compound heterozygous mutations in HEPHL1 (NM_001098672) who presented with abnormal hair (pili torti and trichorrhexis nodosa) and cognitive dysfunction. The maternal missense mutation affected mRNA splicing, leading to skipping of exon 5 and causing an in-frame deletion of 85 amino acids (c.809_1063del; p.Leu271_ala355del). The paternal mutation (c.3176T>C; p.Met1059Thr) changed a highly conserved methionine that is part of a typical type I copper binding site in HEPHL1. We demonstrated that HEPHL1 has ferroxidase activity and that the patient's two mutations exhibited loss of this ferroxidase activity. Consistent with these findings, the patient's fibroblasts accumulated intracellular iron and exhibited reduced activity of the copper-dependent enzyme, lysyl oxidase. These results suggest that the patient's biallelic variants are loss-of-function mutations. Hence, we generated a Hephl1 knockout mouse model that was viable and had curly whiskers, consistent with the hair phenotype in our patient. These results enhance our understanding of the function of HEPHL1 and implicate altered ferroxidase activity in hair growth and hair disorders.

PARP1 inhibition alleviates injury in ARH3-deficient mice and human cells.

Poly(ADP-ribosyl)ation refers to the covalent attachment of ADP-ribose to protein, generating branched, long chains of ADP-ribose moieties, known as poly(ADP-ribose) (PAR). Poly(ADP-ribose) polymerase 1 (PARP1) is the main polymerase and acceptor of PAR in response to DNA damage. Excessive intracellular PAR accumulation due to PARP1 activation leads cell death in a pathway known as parthanatos. PAR degradation is mainly controlled by poly(ADP-ribose) glycohydrolase (PARG) and ADP-ribose-acceptor hydrolase 3 (ARH3). Our previous results demonstrated that ARH3 confers protection against hydrogen peroxide (H2O2) exposure, by lowering cytosolic and nuclear PAR levels and preventing apoptosis-inducing factor (AIF) nuclear translocation. We identified a family with an ARH3 gene mutation that resulted in a truncated, inactive protein. The 8-year-old proband exhibited a progressive neurodegeneration phenotype. In addition, parthanatos was observed in neurons of the patient's deceased sibling, and an older sibling exhibited a mild behavioral phenotype. Consistent with the previous findings, the patient's fibroblasts and ARH3-deficient mice were more sensitive, respectively, to H2O2 stress and cerebral ischemia/reperfusion-induced PAR accumulation and cell death. Further, PARP1 inhibition alleviated cell death and injury resulting from oxidative stress and ischemia/reperfusion. PARP1 inhibitors may attenuate the progression of neurodegeneration in affected patients with ARH3 deficiency.

A suite of automated sequence analyses reduces the number of candidate deleterious variants and reveals a difference between probands and unaffected siblings.

PURPOSE: Develop an automated exome analysis workflow that can produce a very small number of candidate variants yet still detect different numbers of deleterious variants between probands and unaffected siblings. METHODS: Ninety-seven outbred nuclear families from the Undiagnosed Diseases Program/Network included single probands and the corresponding unaffected sibling(s). Single-nucleotide polymorphism (SNP) chip and exome analyses were performed on all, with proband and unaffected sibling considered independently as the target. The total burden of candidate genetic variants was summed for probands and siblings over all considered disease models. RESULTS: Exome analysis workflow include automated programs for ethnicity-matched genotype calling, salvage pathway for Mendelian inconsistency, compound heterozygous recessive detection, BAM file regional curation, population frequency filtering, pedigree-aware BAM file noise evaluation, and exon deletion filtration. This workflow relied heavily on BAM file analysis. A greater average pathogenic variant number was found compared with unaffected siblings. This was significant (p < 0.05) when using published recommended thresholds, and implies that causal variants are retained in many probands' lists. CONCLUSION: Using Mendelian and non-Mendelian models, this agnostic exome analysis shows a difference between a small group of probands and their unaffected siblings. This workflow produces candidate lists small enough to pursue with laboratory validation.

Expanding the phenotype of COPA syndrome: a kindred with typical and atypical features.

BACKGROUND: Copa syndrome is a rare autosomal dominant disorder with abnormal intracellular vesicle trafficking. The objective of this work is to expand the knowledge about this disorder by delineating phenotypic features of an unreported COPA family. METHODS AND RESULTS: A heterozygous missense variant (c.698 G>A, p.Arg233His) in COPA was identified in four members of a three-generation kindred with lung, autoimmune and malignant disease of unknown aetiology. Ages of onset were 56, 26, 16 and 1 year, with earlier age of onset in successive generations. Presenting symptoms were cough and dyspnoea. Findings included small lung cysts, follicular bronchiolitis, interstitial lung disease, neuroendocrine cell hyperplasia, rheumatoid arthritis, avascular necrosis and select abnormal autoimmune serologies. Neither alveolar haemorrhage nor glomerular disease were present. Features not previously associated with Copa syndrome included neuromyelitis optica, pulmonary carcinoid tumour, clear cell renal carcinoma, renal cysts, hepatic cysts, nephrolithiasis, pyelonephritis and meningitis. Longitudinal evaluations demonstrated slow progression of lung disease and extrapulmonary cysts. CONCLUSIONS: Worsening severity with successive generations may be observed in Copa syndrome. Extrapulmonary cysts, malignancies, autoimmune neurological disorders and infections are clinical features that may be associated with Copa syndrome. Further studies are indicated to fully define the phenotypic spectrum of this disorder.

Early infantile-onset epileptic encephalopathy 28 due to a homozygous microdeletion involving the WWOX gene in a region of uniparental disomy.

The genetic etiologies of many rare disorders, including early infantile epileptic encephalopathies, are largely undiagnosed. A 6-year-old girl was admitted to the National Institutes of Health Undiagnosed Diseases Program with profound intellectual disability, infantile-onset seizures, chronic respiratory failure, facial dysmorphisms, skeletal abnormalities, and atrial septum defect. A large region of homozygosity was discovered on chromosome 16, spanning 16q22.1-16q24.3' caused by uniparental disomy (UPD) that included a maternally inherited homozygous microdeletion covering exon 6 of WWOX (NM_016373.3). mRNA expression analysis revealed that the deletion led to nonsense-mediated decay of the NM_016373.3 transcript; the exon 6 of an alternative transcript (NM_130791.3), lacking the short-chain dehydrogenase, was utilized. The microdeletion in WWOX explains the seizures and intellectual disability, while pathogenic variants in another gene, HSPG2, are likely responsible for the patient's skeletal abnormalities. This report describes a rare autosomal recessive disorder with multiple genetic etiologies, one of which involves UPD.

Novel variants in SPTAN1 without epilepsy: An expansion of the phenotype.

We describe two unrelated children with de novo variants in the non-erythrocytic alpha-II-spectrin (SPTAN1) gene who have hypoplastic brain structures, intellectual disability, and both fine and gross motor impairments. Using agnostic exome sequencing, we identified a nonsense variant creating a premature stop codon in exon 21 of SPTAN1, and in a second patient we identified an intronic substitution in SPTAN1 prior to exon 50 creating a new donor acceptor site. Neither of these variants has been described previously. Although some of these patients' features are consistent with the known SPTAN1 encephalopathy phenotype, these two children do not have epilepsy, in contrast to reports about nearly every other patient with heterozygous SPTAN1 variants and in all patients with a variant near the C-terminal coding region. Moreover, both children have abnormal thyroid function, which has not been previously reported in association with SPTAN1 variant. We present a detailed discussion of the clinical manifestations of these two unique SPTAN1 variants and provide evidence that both variants result in reduced mRNA expression despite different locations within the gene and clinical phenotypes. These findings expand the motor, cognitive, and behavioral spectrum of the SPTAN1-associated phenotype and invite speculation about underlying pathophysiologies.

Bi-allelic CCDC47 Variants Cause a Disorder Characterized by Woolly Hair, Liver Dysfunction, Dysmorphic Features, and Global Developmental Delay.

Ca2+ signaling is vital for various cellular processes including synaptic vesicle exocytosis, muscle contraction, regulation of secretion, gene transcription, and cellular proliferation. The endoplasmic reticulum (ER) is the largest intracellular Ca2+ store, and dysregulation of ER Ca2+ signaling and homeostasis contributes to the pathogenesis of various complex disorders and Mendelian disease traits. We describe four unrelated individuals with a complex multisystem disorder characterized by woolly hair, liver dysfunction, pruritus, dysmorphic features, hypotonia, and global developmental delay. Through whole-exome sequencing and family-based genomics, we identified bi-allelic variants in CCDC47 that encodes the Ca2+-binding ER transmembrane protein CCDC47. CCDC47, also known as calumin, has been shown to bind Ca2+ with low affinity and high capacity. In mice, loss of Ccdc47 leads to embryonic lethality, suggesting that Ccdc47 is essential for early development. Characterization of cells from individuals with predicted likely damaging alleles showed decreased CCDC47 mRNA expression and protein levels. In vitro cellular experiments showed decreased total ER Ca2+ storage, impaired Ca2+ signaling mediated by the IP3R Ca2+ release channel, and reduced ER Ca2+ refilling via store-operated Ca2+ entry. These results, together with the previously described role of CCDC47 in Ca2+ signaling and development, suggest that bi-allelic loss-of-function variants in CCDC47 underlie the pathogenesis of this multisystem disorder.

KCTD7 deficiency defines a distinct neurodegenerative disorder with a conserved autophagy-lysosome defect.

OBJECTIVE: Several small case series identified KCTD7 mutations in patients with a rare autosomal recessive disorder designated progressive myoclonic epilepsy (EPM3) and neuronal ceroid lipofuscinosis (CLN14). Despite the name KCTD (potassium channel tetramerization domain), KCTD protein family members lack predicted channel domains. We sought to translate insight gained from yeast studies to uncover disease mechanisms associated with deficiencies in KCTD7 of unknown function. METHODS: Novel KCTD7 variants in new and published patients were assessed for disease causality using genetic analyses, cell-based functional assays of patient fibroblasts and knockout yeast, and electron microscopy of patient samples. RESULTS: Patients with KCTD7 mutations can exhibit movement disorders or developmental regression before seizure onset, and are distinguished from similar disorders by an earlier age of onset. Although most published KCTD7 patient variants were excluded from a genome sequence database of normal human variations, most newly identified patient variants are present in this database, potentially challenging disease causality. However, genetic analysis and impaired biochemical interactions with cullin 3 support a causal role for patient KCTD7 variants, suggesting deleterious alleles of KCTD7 and other rare disease variants may be underestimated. Both patient-derived fibroblasts and yeast lacking Whi2 with sequence similarity to KCTD7 have impaired autophagy consistent with brain pathology. INTERPRETATION: Biallelic KCTD7 mutations define a neurodegenerative disorder with lipofuscin and lipid droplet accumulation but without defining features of neuronal ceroid lipofuscinosis or lysosomal storage disorders. KCTD7 deficiency appears to cause an underlying autophagy-lysosome defect conserved in yeast, thereby assigning a biological role for KCTD7. Ann Neurol 2018;84:774-788.

Aggrecan Mutations in Nonfamilial Short Stature and Short Stature Without Accelerated Skeletal Maturation.

Aggrecan, a proteoglycan, is an important component of cartilage extracellular matrix, including that of the growth plate. Heterozygous mutations in ACAN, the gene encoding aggrecan, cause autosomal dominant short stature, accelerated skeletal maturation, and joint disease. The inheritance pattern and the presence of bone age equal to or greater than chronological age have been consistent features, serving as diagnostic clues. From family 1, a 6-year-old boy presented with short stature [height standard deviation score (SDS), -1.75] and bone age advanced by 3 years. There was no family history of short stature (height SDS: father, -0.76; mother, 0.7). Exome sequencing followed by Sanger sequencing identified a de novo novel heterozygous frameshift mutation in ACAN (c.6404delC: p.A2135Dfs). From family 2, a 12-year-old boy was evaluated for short stature (height SDS, -3.9). His bone age at the time of genetic evaluation was approximately 1 year less than his chronological age. Family history was consistent with an autosomal dominant inheritance of short stature, with several affected members also showing early-onset osteoarthritis. Exome sequencing, confirmed by Sanger sequencing, identified a novel nonsense mutation in ACAN (c.4852C>T: p.Q1618X), which cosegregated with the phenotype. In conclusion, patients with ACAN mutations may present with nonfamilial short stature and with bone age less than chronological age. These findings expand the known phenotypic spectrum of heterozygous ACAN mutations and indicate that this diagnosis should be considered in children without a family history of short stature and in children without accelerated skeletal maturation.

Compound heterozygosity for loss-of-function GARS variants results in a multisystem developmental syndrome that includes severe growth retardation.

Aminoacyl-tRNA synthetases (ARSs) are ubiquitously expressed enzymes that ligate amino acids onto tRNA molecules. Genes encoding ARSs have been implicated in myriad dominant and recessive disease phenotypes. Glycyl-tRNA synthetase (GARS) is a bifunctional ARS that charges tRNAGly in the cytoplasm and mitochondria. GARS variants have been associated with dominant Charcot-Marie-Tooth disease but have not been convincingly implicated in recessive phenotypes. Here, we describe a patient from the NIH Undiagnosed Diseases Program with a multisystem, developmental phenotype. Whole-exome sequence analysis revealed that the patient is compound heterozygous for one frameshift (p.Glu83Ilefs*6) and one missense (p.Arg310Gln) GARS variant. Using in vitro and in vivo functional studies, we show that both GARS variants cause a loss-of-function effect: the frameshift variant results in depleted protein levels and the missense variant reduces GARS tRNA charging activity. In support of GARS variant pathogenicity, our patient shows striking phenotypic overlap with other patients having ARS-related recessive diseases, including features associated with variants in both cytoplasmic and mitochondrial ARSs; this observation is consistent with the essential function of GARS in both cellular locations. In summary, our clinical, genetic, and functional analyses expand the phenotypic spectrum associated with GARS variants.