A series of simplex cases have been reported under various diagnoses sharing early aging, especially evident in congenitally decreased subcutaneous fat tissue and sparse hair, bone dysplasia of the skull and fingers, a distinctive facial gestalt, and prenatal and postnatal growth retardation. For historical reasons, we suggest naming the entity Fontaine syndrome. Exome sequencing of four unrelated affected individuals showed that all carried the de novo missense variant c.649C>T (p.Arg217Cys) or c.650G>A (p.Arg217His) in SLC25A24, a solute carrier 25 family member coding for calcium-binding mitochondrial carrier protein (SCaMC-1, also known as SLC25A24). SLC25A24 allows an electro-neutral and reversible exchange of ATP-Mg and phosphate between the cytosol and mitochondria, which is required for maintaining optimal adenine nucleotide levels in the mitochondrial matrix. Molecular dynamic simulation studies predict that p.Arg217Cys and p.Arg217His narrow the substrate cavity of the protein and disrupt transporter dynamics. SLC25A24-mutant fibroblasts and cells expressing p.Arg217Cys or p.Arg217His variants showed altered mitochondrial morphology, a decreased proliferation rate, increased mitochondrial membrane potential, and decreased ATP-linked mitochondrial oxygen consumption. The results suggest that the SLC25A24 mutations lead to impaired mitochondrial ATP synthesis and cause hyperpolarization and increased proton leak in association with an impaired energy metabolism. Our findings identify SLC25A24 mutations affecting codon 217 as the underlying genetic cause of human progeroid Fontaine syndrome.

By using whole-exome sequencing, we identified a homozygous guanine-to-adenine transition at the invariant -1 position of the acceptor site of intron 1 (c.97-1G>A) in solute carrier organic anion transporter family member 2A1 (SLCO2A1), which encodes a prostaglandin transporter protein, as the causative mutation in a single individual with primary hypertrophic osteoarthropathy (PHO) from a consanguineous family. In two other affected individuals with PHO from two unrelated nonconsanguineous families, we identified two different compound heterozygous mutations by using Sanger sequencing. These findings confirm that SLCO2A1 mutations inactivate prostaglandin E(2) (PGE(2)) transport, and they indicate that mutations in SLCO2A1 are the pathogenic cause of PHO. Moreover, this study might also help to explain the cause of secondary hypertrophic osteoarthropathy.

Secreted protein, acidic, cysteine-rich (SPARC) is a glycoprotein that binds to collagen type I and other proteins in the extracellular matrix. Using whole-exome sequencing to identify the molecular defect in two unrelated girls with severe bone fragility and a clinical diagnosis of osteogenesis imperfecta type IV, we identified two homozygous variants in SPARC (GenBank: NM_003118.3; c.497G>A [p.Arg166His] in individual 1; c.787G>A [p.Glu263Lys] in individual 2). Published modeling and site-directed mutagenesis studies had previously shown that the residues substituted by these mutations form an intramolecular salt bridge in SPARC and are essential for the binding of SPARC to collagen type I. The amount of SPARC secreted by skin fibroblasts was reduced in individual 1 but appeared normal in individual 2. The migration of collagen type I alpha chains produced by these fibroblasts was mildly delayed on SDS-PAGE gel, suggesting some overmodification of collagen during triple helical formation. Pulse-chase experiments showed that collagen type I secretion was mildly delayed in skin fibroblasts from both individuals. Analysis of an iliac bone sample from individual 2 showed that trabecular bone was hypermineralized on the material level. In conclusion, these observations show that homozygous mutations in SPARC can give rise to severe bone fragility in humans.

Skeletal dysplasias are highly variable Mendelian phenotypes. Molecular diagnosis of skeletal dysplasias is complicated by their extreme clinical and genetic heterogeneity. We describe a clinically recognizable autosomal-recessive disorder in four affected siblings from a consanguineous Saudi family, comprising progressive spondyloepimetaphyseal dysplasia, short stature, facial dysmorphism, short fourth metatarsals, and intellectual disability. Combined autozygome/exome analysis identified a homozygous frameshift mutation in RSPRY1 with resulting nonsense-mediated decay. Using a gene-centric "matchmaking" system, we were able to identify a Peruvian simplex case subject whose phenotype is strikingly similar to the original Saudi family and whose exome sequencing had revealed a likely pathogenic homozygous missense variant in the same gene. RSPRY1 encodes a hypothetical RING and SPRY domain-containing protein of unknown physiological function. However, we detect strong RSPRY1 protein localization in murine embryonic osteoblasts and periosteal cells during primary endochondral ossification, consistent with a role in bone development. This study highlights the role of gene-centric matchmaking tools to establish causal links to genes, especially for rare or previously undescribed clinical entities.

The natriuretic peptide signaling pathway has been implicated in many cellular processes, including endochondral ossification and bone growth. More precisely, different mutations in the NPR-B receptor and the CNP ligand have been identified in individuals with either short or tall stature. In this study we show that the NPR-C receptor (encoded by NPR3) is also important for the regulation of linear bone growth. We report four individuals, originating from three different families, with a phenotype characterized by tall stature, long digits, and extra epiphyses in the hands and feet. In addition, aortic dilatation was observed in two of these families. In each affected individual, we identified a bi-allelic loss-of-function mutation in NPR3. The missense mutations (c.442T>C [p.Ser148Pro] and c.1088A>T [p.Asp363Val]) resulted in intracellular retention of the NPR-C receptor and absent localization on the plasma membrane, whereas the nonsense mutation (c.1524delC [p.Tyr508∗]) resulted in nonsense-mediated mRNA decay. Biochemical analysis of plasma from two affected and unrelated individuals revealed a reduced NTproNP/NP ratio for all ligands and also high cGMP levels. These data strongly suggest a reduced clearance of natriuretic peptides by the defective NPR-C receptor and consequently increased activity of the NPR-A/B receptors. In conclusion, this study demonstrates that loss-of-function mutations in NPR3 result in increased NPR-A/B signaling activity and cause a phenotype marked by enhanced bone growth and cardiovascular abnormalities.

Ribosomal protein (RP) gene mutations, mostly associated with inherited or acquired bone marrow failure, are believed to drive disease by slowing the rate of protein synthesis. Here de novo missense mutations in the RPS23 gene, which codes for uS12, are reported in two unrelated individuals with microcephaly, hearing loss, and overlapping dysmorphic features. One individual additionally presents with intellectual disability and autism spectrum disorder. The amino acid substitutions lie in two highly conserved loop regions of uS12 with known roles in maintaining the accuracy of mRNA codon translation. Primary cells revealed one substitution severely impaired OGFOD1-dependent hydroxylation of a neighboring proline residue resulting in 40S ribosomal subunits that were blocked from polysome formation. The other disrupted a predicted pi-pi stacking interaction between two phenylalanine residues leading to a destabilized uS12 that was poorly tolerated in 40S subunit biogenesis. Despite no evidence of a reduction in the rate of mRNA translation, these uS12 variants impaired the accuracy of mRNA translation and rendered cells highly sensitive to oxidative stress. These discoveries describe a ribosomopathy linked to uS12 and reveal mechanistic distinctions between RP gene mutations driving hematopoietic disease and those resulting in developmental disorders.

SQSTM1 (sequestosome 1; also known as p62) encodes a multidomain scaffolding protein involved in various key cellular processes, including the removal of damaged mitochondria by its function as a selective autophagy receptor. Heterozygous variants in SQSTM1 have been associated with Paget disease of the bone and might contribute to neurodegeneration in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Using exome sequencing, we identified three different biallelic loss-of-function variants in SQSTM1 in nine affected individuals from four families with a childhood- or adolescence-onset neurodegenerative disorder characterized by gait abnormalities, ataxia, dysarthria, dystonia, vertical gaze palsy, and cognitive decline. We confirmed absence of the SQSTM1/p62 protein in affected individuals' fibroblasts and found evidence of a defect in the early response to mitochondrial depolarization and autophagosome formation. Our findings expand the SQSTM1-associated phenotypic spectrum and lend further support to the concept of disturbed selective autophagy pathways in neurodegenerative diseases.

The periosteum contributes to bone repair and maintenance of cortical bone mass. In contrast to the understanding of bone development within the epiphyseal growth plate, factors that regulate periosteal osteogenesis have not been studied as intensively. Osteofibrous dysplasia (OFD) is a congenital disorder of osteogenesis and is typically sporadic and characterized by radiolucent lesions affecting the cortical bone immediately under the periosteum of the tibia and fibula. We identified germline mutations in MET, encoding a receptor tyrosine kinase, that segregate with an autosomal-dominant form of OFD in three families and a mutation in a fourth affected subject from a simplex family and with bilateral disease. Mutations identified in all families with dominant inheritance and in the one simplex subject with bilateral disease abolished the splice inclusion of exon 14 in MET transcripts, which resulted in a MET receptor (MET(Δ14)) lacking a cytoplasmic juxtamembrane domain. Splice exclusion of this domain occurs during normal embryonic development, and forced induction of this exon-exclusion event retarded osteoblastic differentiation in vitro and inhibited bone-matrix mineralization. In an additional subject with unilateral OFD, we identified a somatic MET mutation, also affecting exon 14, that substituted a tyrosine residue critical for MET receptor turnover and, as in the case of the MET(Δ14) mutations, had a stabilizing effect on the mature protein. Taken together, these data show that aberrant MET regulation via the juxtamembrane domain subverts core MET receptor functions that regulate osteogenesis within cortical diaphyseal bone.

Disorders of Golgi homeostasis form an emerging group of genetic defects. The highly heterogeneous clinical spectrum is not explained by our current understanding of the underlying cell-biological processes in the Golgi. Therefore, uncovering genetic defects and annotating gene function are challenging. Exome sequencing in a family with three siblings affected by abnormal Golgi glycosylation revealed a homozygous missense mutation, c.92T>C (p.Leu31Ser), in coiled-coil domain containing 115 (CCDC115), the function of which is unknown. The same mutation was identified in three unrelated families, and in one family it was compound heterozygous in combination with a heterozygous deletion of CCDC115. An additional homozygous missense mutation, c.31G>T (p.Asp11Tyr), was found in a family with two affected siblings. All individuals displayed a storage-disease-like phenotype involving hepatosplenomegaly, which regressed with age, highly elevated bone-derived alkaline phosphatase, elevated aminotransferases, and elevated cholesterol, in combination with abnormal copper metabolism and neurological symptoms. Two individuals died of liver failure, and one individual was successfully treated by liver transplantation. Abnormal N- and mucin type O-glycosylation was found on serum proteins, and reduced metabolic labeling of sialic acids was found in fibroblasts, which was restored after complementation with wild-type CCDC115. PSI-BLAST homology detection revealed reciprocal homology with Vma22p, the yeast V-ATPase assembly factor located in the endoplasmic reticulum (ER). Human CCDC115 mainly localized to the ERGIC and to COPI vesicles, but not to the ER. These data, in combination with the phenotypic spectrum, which is distinct from that associated with defects in V-ATPase core subunits, suggest a more general role for CCDC115 in Golgi trafficking. Our study reveals CCDC115 deficiency as a disorder of Golgi homeostasis that can be readily identified via screening for abnormal glycosylation in plasma.

Spermatogenesis-associated 5 like 1 (SPATA5L1) represents an orphan gene encoding a protein of unknown function. We report 28 bi-allelic variants in SPATA5L1 associated with sensorineural hearing loss in 47 individuals from 28 (26 unrelated) families. In addition, 25/47 affected individuals (53%) presented with microcephaly, developmental delay/intellectual disability, cerebral palsy, and/or epilepsy. Modeling indicated damaging effect of variants on the protein, largely via destabilizing effects on protein domains. Brain imaging revealed diminished cerebral volume, thin corpus callosum, and periventricular leukomalacia, and quantitative volumetry demonstrated significantly diminished white matter volumes in several individuals. Immunofluorescent imaging in rat hippocampal neurons revealed localization of Spata5l1 in neuronal and glial cell nuclei and more prominent expression in neurons. In the rodent inner ear, Spata5l1 is expressed in the neurosensory hair cells and inner ear supporting cells. Transcriptomic analysis performed with fibroblasts from affected individuals was able to distinguish affected from controls by principal components. Analysis of differentially expressed genes and networks suggested a role for SPATA5L1 in cell surface adhesion receptor function, intracellular focal adhesions, and DNA replication and mitosis. Collectively, our results indicate that bi-allelic SPATA5L1 variants lead to a human disease characterized by sensorineural hearing loss (SNHL) with or without a nonprogressive mixed neurodevelopmental phenotype.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Colony stimulating factor 1 receptor (CSF1R) plays key roles in regulating development and function of the monocyte/macrophage lineage, including microglia and osteoclasts. Mono-allelic mutations of CSF1R are known to cause hereditary diffuse leukoencephalopathy with spheroids (HDLS), an adult-onset progressive neurodegenerative disorder. Here, we report seven affected individuals from three unrelated families who had bi-allelic CSF1R mutations. In addition to early-onset HDLS-like neurological disorders, they had brain malformations and skeletal dysplasia compatible to dysosteosclerosis (DOS) or Pyle disease. We identified five CSF1R mutations that were homozygous or compound heterozygous in these affected individuals. Two of them were deep intronic mutations resulting in abnormal inclusion of intron sequences in the mRNA. Compared with Csf1r-null mice, the skeletal and neural phenotypes of the affected individuals appeared milder and variable, suggesting that at least one of the mutations in each affected individual is hypomorphic. Our results characterized a unique human skeletal phenotype caused by CSF1R deficiency and implied that bi-allelic CSF1R mutations cause a spectrum of neurological and skeletal disorders, probably depending on the residual CSF1R function.

Replicating the human genome efficiently and accurately is a daunting challenge involving the duplication of upward of three billion base pairs. At the core of the complex machinery that achieves this task are three members of the B family of DNA polymerases: DNA polymerases α, δ, and ε. Collectively these multimeric polymerases ensure DNA replication proceeds at optimal rates approaching 2 × 103 nucleotides/min with an error rate of less than one per million nucleotides polymerized. The majority of DNA replication of undamaged DNA is conducted by DNA polymerases δ and ε. The DNA polymerase α-primase complex performs limited synthesis to initiate the replication process, along with Okazaki-fragment synthesis on the discontinuous lagging strand. An increasing number of human disorders caused by defects in different components of the DNA-replication apparatus have been described to date. These are clinically diverse and involve a wide range of features, including variable combinations of growth delay, immunodeficiency, endocrine insufficiencies, lipodystrophy, and cancer predisposition. Here, by using various complementary approaches, including classical linkage analysis, targeted next-generation sequencing, and whole-exome sequencing, we describe distinct missense and splice-impacting mutations in POLA1 in five unrelated families presenting with an X-linked syndrome involving intellectual disability, proportionate short stature, microcephaly, and hypogonadism. POLA1 encodes the p180 catalytic subunit of DNA polymerase α-primase. A range of replicative impairments could be demonstrated in lymphoblastoid cell lines derived from affected individuals. Our findings describe the presentation of pathogenic mutations in a catalytic component of a B family DNA polymerase member, DNA polymerase α.

Human C2orf69 is an evolutionarily conserved gene whose function is unknown. Here, we report eight unrelated families from which 20 children presented with a fatal syndrome consisting of severe autoinflammation and progredient leukoencephalopathy with recurrent seizures; 12 of these subjects, whose DNA was available, segregated homozygous loss-of-function C2orf69 variants. C2ORF69 bears homology to esterase enzymes, and orthologs can be found in most eukaryotic genomes, including that of unicellular phytoplankton. We found that endogenous C2ORF69 (1) is loosely bound to mitochondria, (2) affects mitochondrial membrane potential and oxidative respiration in cultured neurons, and (3) controls the levels of the glycogen branching enzyme 1 (GBE1) consistent with a glycogen-storage-associated mitochondriopathy. We show that CRISPR-Cas9-mediated inactivation of zebrafish C2orf69 results in lethality by 8 months of age due to spontaneous epileptic seizures, which is preceded by persistent brain inflammation. Collectively, our results delineate an autoinflammatory Mendelian disorder of C2orf69 deficiency that disrupts the development/homeostasis of the immune and central nervous systems.

SPONASTRIME dysplasia is an autosomal-recessive spondyloepimetaphyseal dysplasia characterized by spine (spondylar) abnormalities, midface hypoplasia with a depressed nasal bridge, metaphyseal striations, and disproportionate short stature. Scoliosis, coxa vara, childhood cataracts, short dental roots, and hypogammaglobulinemia have also been reported in this disorder. Although an autosomal-recessive inheritance pattern has been hypothesized, pathogenic variants in a specific gene have not been discovered in individuals with SPONASTRIME dysplasia. Here, we identified bi-allelic variants in TONSL, which encodes the Tonsoku-like DNA repair protein, in nine subjects (from eight families) with SPONASTRIME dysplasia, and four subjects (from three families) with short stature of varied severity and spondylometaphyseal dysplasia with or without immunologic and hematologic abnormalities, but no definitive metaphyseal striations at diagnosis. The finding of early embryonic lethality in a Tonsl-/- murine model and the discovery of reduced length, spinal abnormalities, reduced numbers of neutrophils, and early lethality in a tonsl-/- zebrafish model both support the hypomorphic nature of the identified TONSL variants. Moreover, functional studies revealed increased amounts of spontaneous replication fork stalling and chromosomal aberrations, as well as fewer camptothecin (CPT)-induced RAD51 foci in subject-derived cell lines. Importantly, these cellular defects were rescued upon re-expression of wild-type (WT) TONSL; this rescue is consistent with the hypothesis that hypomorphic TONSL variants are pathogenic. Overall, our studies in humans, mice, zebrafish, and subject-derived cell lines confirm that pathogenic variants in TONSL impair DNA replication and homologous recombination-dependent repair processes, and they lead to a spectrum of skeletal dysplasia phenotypes with numerous extra-skeletal manifestations.

Chromatin is essentially an array of nucleosomes, each of which consists of the DNA double-stranded fiber wrapped around a histone octamer. This organization supports cellular processes such as DNA replication, DNA transcription, and DNA repair in all eukaryotes. Human histone H4 is encoded by fourteen canonical histone H4 genes, all differing at the nucleotide level but encoding an invariant protein. Here, we present a cohort of 29 subjects with de novo missense variants in six H4 genes (H4C3, H4C4, H4C5, H4C6, H4C9, and H4C11) identified by whole-exome sequencing and matchmaking. All individuals present with neurodevelopmental features of intellectual disability and motor and/or gross developmental delay, while non-neurological features are more variable. Ten amino acids are affected, six recurrently, and are all located within the H4 core or C-terminal tail. These variants cluster to specific regions of the core H4 globular domain, where protein-protein interactions occur with either other histone subunits or histone chaperones. Functional consequences of the identified variants were evaluated in zebrafish embryos, which displayed abnormal general development, defective head organs, and reduced body axis length, providing compelling evidence for the causality of the reported disorder(s). While multiple developmental syndromes have been linked to chromatin-associated factors, missense-bearing histone variants (e.g., H3 oncohistones) are only recently emerging as a major cause of pathogenicity. Our findings establish a broader involvement of H4 variants in developmental syndromes.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.