The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing.

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

A Genocentric Approach to Discovery of Mendelian Disorders.

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Bi-allelic Variants in TONSL Cause SPONASTRIME Dysplasia and a Spectrum of Skeletal Dysplasia Phenotypes.

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.

Retrospective analysis of a clinical exome sequencing cohort reveals the mutational spectrum and identifies candidate disease-associated loci for BAFopathies.

PURPOSE: BRG1/BRM-associated factor (BAF) complex is a chromatin remodeling complex that plays a critical role in gene regulation. Defects in the genes encoding BAF subunits lead to BAFopathies, a group of neurodevelopmental disorders with extensive locus and phenotypic heterogeneity. METHODS: We retrospectively analyzed data from 16,243 patients referred for clinical exome sequencing (ES) with a focus on the BAF complex. We applied a genotype-first approach, combining predicted genic constraints to propose candidate BAFopathy genes. RESULTS: We identified 127 patients carrying pathogenic variants, likely pathogenic variants, or de novo variants of unknown clinical significance in 11 known BAFopathy genes. Those include 34 patients molecularly diagnosed using ES reanalysis with new gene-disease evidence (n = 21) or variant reclassifications in known BAFopathy genes (n = 13). We also identified de novo or predicted loss-of-function variants in 4 candidate BAFopathy genes, including ACTL6A, BICRA (implicated in Coffin-Siris syndrome during this study), PBRM1, and SMARCC1. CONCLUSION: We report the mutational spectrum of BAFopathies in an ES cohort. A genotype-driven and pathway-based reanalysis of ES data identified new evidence for candidate genes involved in BAFopathies. Further mechanistic and phenotypic characterization of additional patients are warranted to confirm their roles in human disease and to delineate their associated phenotypic spectrums.

Functional and structural analysis of cytokine-selective IL6ST defects that cause recessive hyper-IgE syndrome.

BACKGROUND: Biallelic variants in IL6ST, encoding GP130, cause a recessive form of hyper-IgE syndrome (HIES) characterized by high IgE level, eosinophilia, defective acute phase response, susceptibility to bacterial infections, and skeletal abnormalities due to cytokine-selective loss of function in GP130, with defective IL-6 and IL-11 and variable oncostatin M (OSM) and IL-27 levels but sparing leukemia inhibitory factor (LIF) signaling. OBJECTIVE: Our aim was to understand the functional and structural impact of recessive HIES-associated IL6ST variants. METHODS: We investigated a patient with HIES by using exome, genome, and RNA sequencing. Functional assays assessed IL-6, IL-11, IL-27, OSM, LIF, CT-1, CLC, and CNTF signaling. Molecular dynamics simulations and structural modeling of GP130 cytokine receptor complexes were performed. RESULTS: We identified a patient with compound heterozygous novel missense variants in IL6ST (p.Ala517Pro and the exon-skipping null variant p.Gly484_Pro518delinsArg). The p.Ala517Pro variant resulted in a more profound IL-6- and IL-11-dominated signaling defect than did the previously identified recessive HIES IL6ST variants p.Asn404Tyr and p.Pro498Leu. Molecular dynamics simulations suggested that the p.Ala517Pro and p.Asn404Tyr variants result in increased flexibility of the extracellular membrane-proximal domains of GP130. We propose a structural model that explains the cytokine selectivity of pathogenic IL6ST variants that result in recessive HIES. The variants destabilized the conformation of the hexameric cytokine receptor complexes, whereas the trimeric LIF-GP130-LIFR complex remained stable through an additional membrane-proximal interaction. Deletion of this membrane-proximal interaction site in GP130 consequently caused additional defective LIF signaling and Stüve-Wiedemann syndrome. CONCLUSION: Our data provide a structural basis to understand clinical phenotypes in patients with IL6ST variants.

De Novo Missense Variants in TRAF7 Cause Developmental Delay, Congenital Anomalies, and Dysmorphic Features.

TRAF7 is a multi-functional protein involved in diverse signaling pathways and cellular processes. The phenotypic consequence of germline TRAF7 variants remains unclear. Here we report missense variants in TRAF7 in seven unrelated individuals referred for clinical exome sequencing. The seven individuals share substantial phenotypic overlap, with developmental delay, congenital heart defects, limb and digital anomalies, and dysmorphic features emerging as key unifying features. The identified variants are de novo in six individuals and comprise four distinct missense changes, including a c.1964G>A (p.Arg655Gln) variant that is recurrent in four individuals. These variants affect evolutionarily conserved amino acids and are located in key functional domains. Gene-specific mutation rate analysis showed that the occurrence of the de novo variants in TRAF7 (p = 2.6 × 10-3) and the recurrent de novo c.1964G>A (p.Arg655Gln) variant (p = 1.9 × 10-8) in our exome cohort was unlikely to have occurred by chance. In vitro analyses of the observed TRAF7 mutations showed reduced ERK1/2 phosphorylation. Our findings suggest that missense mutations in TRAF7 are associated with a multisystem disorder and provide evidence of a role for TRAF7 in human development.

Effect of Genetic Diagnosis on Patients with Previously Undiagnosed Disease.

BACKGROUND: Many patients remain without a diagnosis despite extensive medical evaluation. The Undiagnosed Diseases Network (UDN) was established to apply a multidisciplinary model in the evaluation of the most challenging cases and to identify the biologic characteristics of newly discovered diseases. The UDN, which is funded by the National Institutes of Health, was formed in 2014 as a network of seven clinical sites, two sequencing cores, and a coordinating center. Later, a central biorepository, a metabolomics core, and a model organisms screening center were added. METHODS: We evaluated patients who were referred to the UDN over a period of 20 months. The patients were required to have an undiagnosed condition despite thorough evaluation by a health care provider. We determined the rate of diagnosis among patients who subsequently had a complete evaluation, and we observed the effect of diagnosis on medical care. RESULTS: A total of 1519 patients (53% female) were referred to the UDN, of whom 601 (40%) were accepted for evaluation. Of the accepted patients, 192 (32%) had previously undergone exome sequencing. Symptoms were neurologic in 40% of the applicants, musculoskeletal in 10%, immunologic in 7%, gastrointestinal in 7%, and rheumatologic in 6%. Of the 382 patients who had a complete evaluation, 132 received a diagnosis, yielding a rate of diagnosis of 35%. A total of 15 diagnoses (11%) were made by clinical review alone, and 98 (74%) were made by exome or genome sequencing. Of the diagnoses, 21% led to recommendations regarding changes in therapy, 37% led to changes in diagnostic testing, and 36% led to variant-specific genetic counseling. We defined 31 new syndromes. CONCLUSIONS: The UDN established a diagnosis in 132 of the 382 patients who had a complete evaluation, yielding a rate of diagnosis of 35%. (Funded by the National Institutes of Health Common Fund.).

PPP3CA truncating variants clustered in the regulatory domain cause early-onset refractory epilepsy.

PPP3CA encodes the catalytic subunit of calcineurin, a calcium-calmodulin-regulated serine-threonine phosphatase. Loss-of-function (LoF) variants in the catalytic domain have been associated with epilepsy, while gain-of-function (GoF) variants in the auto-inhibitory domain cause multiple congenital abnormalities. We herein report five new patients with de novo PPP3CA variants. Interestingly, the two frameshift variants in this study and the six truncating variants reported previously are all located within a 26-amino acid region in the regulatory domain (RD). Patients with a truncating variant had more severe earlier onset seizures compared to patients with a LoF missense variant, while autism spectrum disorder was a more frequent feature in the latter. Expression studies of a truncating variant showed apparent RNA expression from the mutant allele, but no detectable mutant protein. Our data suggest that PPP3CA truncating variants clustered in the RD, causing more severe early-onset refractory epilepsy and representing a type of variants distinct from LoF or GoF missense variants.

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.

Heterozygous variants in SPTBN1 cause intellectual disability and autism.

Spectrins are common components of cytoskeletons, binding to cytoskeletal elements and the plasma membrane, allowing proper localization of essential membrane proteins, signal transduction, and cellular scaffolding. Spectrins are assembled from α and ß subunits, encoded by SPTA1 and SPTAN1 (α) and SPTB, SPTBN1, SPTBN2, SPTBN4, and SPTBN5 (ß). Pathogenic variants in various spectrin genes are associated with erythroid cell disorders (SPTA1, SPTB) and neurologic disorders (SPTAN1, SPTBN2, and SPTBN4), but no phenotypes have been definitively associated with variants in SPTBN1 or SPTBN5. Through exome sequencing and case matching, we identified seven unrelated individuals with heterozygous SPTBN1 variants: two with de novo missense variants and five with predicted loss-of-function variants (found to be de novo in two, while one was inherited from a mother with a history of learning disabilities). Common features include global developmental delays, intellectual disability, and behavioral disturbances. Autistic features (4/6) and epilepsy (2/7) or abnormal electroencephalogram without overt seizures (1/7) were present in a subset. Identification of loss-of-function variants suggests a haploinsufficiency mechanism, but additional functional studies are required to fully elucidate disease pathogenesis. Our findings support the essential roles of SPTBN1 in human neurodevelopment and expand the knowledge of human spectrinopathy disorders.

Biallelic Variants in OTUD6B Cause an Intellectual Disability Syndrome Associated with Seizures and Dysmorphic Features.

Ubiquitination is a posttranslational modification that regulates many cellular processes including protein degradation, intracellular trafficking, cell signaling, and protein-protein interactions. Deubiquitinating enzymes (DUBs), which reverse the process of ubiquitination, are important regulators of the ubiquitin system. OTUD6B encodes a member of the ovarian tumor domain (OTU)-containing subfamily of deubiquitinating enzymes. Herein, we report biallelic pathogenic variants in OTUD6B in 12 individuals from 6 independent families with an intellectual disability syndrome associated with seizures and dysmorphic features. In subjects with predicted loss-of-function alleles, additional features include global developmental delay, microcephaly, absent speech, hypotonia, growth retardation with prenatal onset, feeding difficulties, structural brain abnormalities, congenital malformations including congenital heart disease, and musculoskeletal features. Homozygous Otud6b knockout mice were subviable, smaller in size, and had congenital heart defects, consistent with the severity of loss-of-function variants in humans. Analysis of peripheral blood mononuclear cells from an affected subject showed reduced incorporation of 19S subunits into 26S proteasomes, decreased chymotrypsin-like activity, and accumulation of ubiquitin-protein conjugates. Our findings suggest a role for OTUD6B in proteasome function, establish that defective OTUD6B function underlies a multisystemic human disorder, and provide additional evidence for the emerging relationship between the ubiquitin system and human disease.

Resolution of Disease Phenotypes Resulting from Multilocus Genomic Variation.

BACKGROUND: Whole-exome sequencing can provide insight into the relationship between observed clinical phenotypes and underlying genotypes. METHODS: We conducted a retrospective analysis of data from a series of 7374 consecutive unrelated patients who had been referred to a clinical diagnostic laboratory for whole-exome sequencing; our goal was to determine the frequency and clinical characteristics of patients for whom more than one molecular diagnosis was reported. The phenotypic similarity between molecularly diagnosed pairs of diseases was calculated with the use of terms from the Human Phenotype Ontology. RESULTS: A molecular diagnosis was rendered for 2076 of 7374 patients (28.2%); among these patients, 101 (4.9%) had diagnoses that involved two or more disease loci. We also analyzed parental samples, when available, and found that de novo variants accounted for 67.8% (61 of 90) of pathogenic variants in autosomal dominant disease genes and 51.7% (15 of 29) of pathogenic variants in X-linked disease genes; both variants were de novo in 44.7% (17 of 38) of patients with two monoallelic variants. Causal copy-number variants were found in 12 patients (11.9%) with multiple diagnoses. Phenotypic similarity scores were significantly lower among patients in whom the phenotype resulted from two distinct mendelian disorders that affected different organ systems (50 patients) than among patients with disorders that had overlapping phenotypic features (30 patients) (median score, 0.21 vs. 0.36; P=1.77×10-7). CONCLUSIONS: In our study, we found multiple molecular diagnoses in 4.9% of cases in which whole-exome sequencing was informative. Our results show that structured clinical ontologies can be used to determine the degree of overlap between two mendelian diseases in the same patient; the diseases can be distinct or overlapping. Distinct disease phenotypes affect different organ systems, whereas overlapping disease phenotypes are more likely to be caused by two genes encoding proteins that interact within the same pathway. (Funded by the National Institutes of Health and the Ting Tsung and Wei Fong Chao Foundation.).

Correction: Clinical experience with carrier screening in a general population: support for a comprehensive pan-ethnic approach.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Clinical experience with carrier screening in a general population: support for a comprehensive pan-ethnic approach.

PURPOSE: To present results from a large cohort of individuals receiving expanded carrier screening (CS) in the United States. METHODS: Single-gene disorder carrier status for 381,014 individuals was determined using next-generation sequencing (NGS) based CS for up to 274 genes. Detection rates were compared with literature-reported values derived from disease prevalence and carrier frequencies. Combined theoretical affected pregnancy rates for the 274 screened disorders were calculated. RESULTS: For Ashkenazi Jewish (AJ) diseases, 81.6% (4434/5435) of carriers identified did not report AJ ancestry. For cystic fibrosis, 44.0% (6260/14,229) of carriers identified had a variant not on the standard genotyping panel. Individuals at risk of being a silent spinal muscular atrophy carrier, not detectable by standard screening, comprised 1/39 (8763/344,407) individuals. For fragile X syndrome, compared with standard premutation screening, AGG interruption analysis modified risk in 83.2% (1128/1356) premutation carriers. Assuming random pairing across the study population, approximately 1/175 pregnancies would be affected by a disorder in the 274-gene screening panel. CONCLUSION: Compared with standard screening, NGS-based CS provides additional information that may impact reproductive choices. Pan-ethnic CS leads to substantially increased identification of at-risk couples. These data support offering NGS-based CS to all reproductive-aged women.

Low-level parental somatic mosaic SNVs in exomes from a large cohort of trios with diverse suspected Mendelian conditions.

PURPOSE: The goal of this study was to assess the scale of low-level parental mosaicism in exome sequencing (ES) databases. METHODS: We analyzed approximately 2000 family trio ES data sets from the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) and Baylor Genetics (BG). Among apparent de novo single-nucleotide variants identified in the affected probands, we selected rare unique variants with variant allele fraction (VAF) between 30% and 70% in the probands and lower than 10% in one of the parents. RESULTS: Of 102 candidate mosaic variants validated using amplicon-based next-generation sequencing, droplet digital polymerase chain reaction, or blocker displacement amplification, 27 (26.4%) were confirmed to be low- (VAF between 1% and 10%) or very low (VAF <1%) level mosaic. Detection precision in parental samples with two or more alternate reads was 63.6% (BHCMG) and 43.6% (BG). In nine investigated individuals, we observed variability of mosaic ratios among blood, saliva, fibroblast, buccal, hair, and urine samples. CONCLUSION: Our computational pipeline enables robust discrimination between true and false positive candidate mosaic variants and efficient detection of low-level mosaicism in ES samples. We confirm that the presence of two or more alternate reads in the parental sample is a reliable predictor of low-level parental somatic mosaicism.

CNVs cause autosomal recessive genetic diseases with or without involvement of SNV/indels.

PURPOSE: Improved resolution of molecular diagnostic technologies enabled detection of smaller sized exonic level copy-number variants (CNVs). The contribution of CNVs to autosomal recessive (AR) conditions may be better recognized using a large clinical cohort. METHODS: We retrospectively investigated the CNVs' contribution to AR conditions in cases subjected to chromosomal microarray analysis (CMA, N = ~70,000) and/or clinical exome sequencing (ES, N = ~12,000) at Baylor Genetics; most had pediatric onset neurodevelopmental disorders. RESULTS: CNVs contributed to biallelic variations in 87 cases, including 81 singletons and three affected sibling pairs. Seventy cases had CNVs affecting both alleles, and 17 had a CNV and a single-nucleotide variant (SNV)/indel in trans. In total, 94.3% of AR-CNVs affected one gene; among these 41.4% were single-exon and 35.0% were multiexon partial-gene events. Sixty-nine percent of homozygous AR-CNVs were embedded in homozygous genomic intervals. Five cases had large deletions unmasking an SNV/indel on the intact allele for a recessive condition, resulting in multiple molecular diagnoses. CONCLUSIONS: AR-CNVs are often smaller in size, transmitted through generations, and underrecognized due to limitations in clinical CNV detection methods. Our findings from a large clinical cohort emphasized integrated CNV and SNV/indel analyses for precise clinical and molecular diagnosis especially in the context of genomic disorders.

Biallelic loss of function variants in PPP1R21 cause a neurodevelopmental syndrome with impaired endocytic function.

Next-generation sequencing (NGS) has been instrumental in solving the genetic basis of rare inherited diseases, especially neurodevelopmental syndromes. However, functional workup is essential for precise phenotype definition and to understand the underlying disease mechanisms. Using whole exome (WES) and whole genome sequencing (WGS) in four independent families with hypotonia, neurodevelopmental delay, facial dysmorphism, loss of white matter, and thinning of the corpus callosum, we identified four previously unreported homozygous truncating PPP1R21 alleles: c.347delT p.(Ile116Lysfs*25), c.2170_2171insGGTA p.(Ile724Argfs*8), c.1607dupT p.(Leu536Phefs*7), c.2063delA p.(Lys688Serfs*26) and found that PPP1R21 was absent in fibroblasts of an affected individual, supporting the allele's loss of function effect. PPP1R21 function had not been studied except that a large scale affinity proteomics approach suggested an interaction with PIBF1 defective in Joubert syndrome. Our co-immunoprecipitation studies did not confirm this but in contrast defined the localization of PPP1R21 to the early endosome. Consistent with the subcellular expression pattern and the clinical phenotype exhibiting features of storage diseases, we found patient fibroblasts exhibited a delay in clearance of transferrin-488 while uptake was normal. In summary, we delineate a novel neurodevelopmental syndrome caused by biallelic PPP1R21 loss of function variants, and suggest a role of PPP1R21 within the endosomal sorting process or endosome maturation pathway.

Bi-allelic Mutations in PKD1L1 Are Associated with Laterality Defects in Humans.

Disruption of the establishment of left-right (L-R) asymmetry leads to situs anomalies ranging from situs inversus totalis (SIT) to situs ambiguus (heterotaxy). The genetic causes of laterality defects in humans are highly heterogeneous. Via whole-exome sequencing (WES), we identified homozygous mutations in PKD1L1 from three affected individuals in two unrelated families. PKD1L1 encodes a polycystin-1-like protein and its loss of function is known to cause laterality defects in mouse and medaka fish models. Family 1 had one fetus and one deceased child with heterotaxy and complex congenital heart malformations. WES identified a homozygous splicing mutation, c.6473+2_6473+3delTG, which disrupts the invariant splice donor site in intron 42, in both affected individuals. In the second family, a homozygous c.5072G>C (p.Cys1691Ser) missense mutation was detected in an individual with SIT and congenital heart disease. The p.Cys1691Ser substitution affects a highly conserved cysteine residue and is predicted by molecular modeling to disrupt a disulfide bridge essential for the proper folding of the G protein-coupled receptor proteolytic site (GPS) motif. Damaging effects associated with substitutions of this conserved cysteine residue in the GPS motif have also been reported in other genes, namely GPR56, BAI3, and PKD1 in human and lat-1 in C. elegans, further supporting the likely pathogenicity of p.Cys1691Ser in PKD1L1. The identification of bi-allelic PKD1L1 mutations recapitulates previous findings regarding phenotypic consequences of loss of function of the orthologous genes in mice and medaka fish and further expands our understanding of genetic contributions to laterality defects in humans.

De Novo Truncating Variants in SON Cause Intellectual Disability, Congenital Malformations, and Failure to Thrive.

SON is a key component of the spliceosomal complex and a critical mediator of constitutive and alternative splicing. Additionally, SON has been shown to influence cell-cycle progression, genomic integrity, and maintenance of pluripotency in stem cell populations. The clear functional relevance of SON in coordinating essential cellular processes and its presence in diverse human tissues suggests that intact SON might be crucial for normal growth and development. However, the phenotypic effects of deleterious germline variants in SON have not been clearly defined. Herein, we describe seven unrelated individuals with de novo variants in SON and propose that deleterious variants in SON are associated with a severe multisystem disorder characterized by developmental delay, persistent feeding difficulties, and congenital malformations, including brain anomalies.

Recurrent De Novo and Biallelic Variation of ATAD3A, Encoding a Mitochondrial Membrane Protein, Results in Distinct Neurological Syndromes.

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of borR534W, the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of borWT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations.