Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease.

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

Rare germline heterozygous missense variants in BRCA1-associated protein 1, BAP1, cause a syndromic neurodevelopmental disorder.

Nuclear deubiquitinase BAP1 (BRCA1-associated protein 1) is a core component of multiprotein complexes that promote transcription by reversing the ubiquitination of histone 2A (H2A). BAP1 is a tumor suppressor whose germline loss-of-function variants predispose to cancer. To our knowledge, there are very rare examples of different germline variants in the same gene causing either a neurodevelopmental disorder (NDD) or a tumor predisposition syndrome. Here, we report a series of 11 de novo germline heterozygous missense BAP1 variants associated with a rare syndromic NDD. Functional analysis showed that most of the variants cannot rescue the consequences of BAP1 inactivation, suggesting a loss-of-function mechanism. In T cells isolated from two affected children, H2A deubiquitination was impaired. In matching peripheral blood mononuclear cells, histone H3 K27 acetylation ChIP-seq indicated that these BAP1 variants induced genome-wide chromatin state alterations, with enrichment for regulatory regions surrounding genes of the ubiquitin-proteasome system (UPS). Altogether, these results define a clinical syndrome caused by rare germline missense BAP1 variants that alter chromatin remodeling through abnormal histone ubiquitination and lead to transcriptional dysregulation of developmental genes.

A concurrent dual analysis of genomic data augments diagnoses: Experiences of 2 clinical sites in the Undiagnosed Diseases Network.

PURPOSE: Next-generation sequencing (NGS) has revolutionized the diagnostic process for rare/ultrarare conditions. However, diagnosis rates differ between analytical pipelines. In the National Institutes of Health-Undiagnosed Diseases Network (UDN) study, each individual's NGS data are concurrently analyzed by the UDN sequencing core laboratory and the clinical sites. We examined the outcomes of this practice. METHODS: A retrospective review was performed at 2 UDN clinical sites to compare the variants and diagnoses/candidate genes identified with the dual analyses of the NGS data. RESULTS: In total, 95 individuals had 100 diagnoses/candidate genes. There was 59% concordance between the UDN sequencing core laboratories and the clinical sites in identifying diagnoses/candidate genes. The core laboratory provided more diagnoses, whereas the clinical sites prioritized more research variants/candidate genes (P < .001). The clinical sites solely identified 15% of the diagnoses/candidate genes. The differences between the 2 pipelines were more often because of variant prioritization disparities than variant detection. CONCLUSION: The unique dual analysis of NGS data in the UDN synergistically enhances outcomes. The core laboratory provided a clinical analysis with more diagnoses and the clinical sites prioritized more research variants/candidate genes. Implementing such concurrent dual analyses in other genomic research studies and clinical settings can improve both variant detection and prioritization.

The benefit of diagnostic whole genome sequencing in schizophrenia and other psychotic disorders.

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families.

Pre-detection history of extensively drug-resistant tuberculosis in KwaZulu-Natal, South Africa.

Antimicrobial-resistant (AMR) infections pose a major threat to global public health. Similar to other AMR pathogens, both historical and ongoing drug-resistant tuberculosis (TB) epidemics are characterized by transmission of a limited number of predominant Mycobacterium tuberculosis (Mtb) strains. Understanding how these predominant strains achieve sustained transmission, particularly during the critical period before they are detected via clinical or public health surveillance, can inform strategies for prevention and containment. In this study, we employ whole-genome sequence (WGS) data from TB clinical isolates collected in KwaZulu-Natal, South Africa to examine the pre-detection history of a successful strain of extensively drug-resistant (XDR) TB known as LAM4/KZN, first identified in a widely reported cluster of cases in 2005. We identify marked expansion of this strain concurrent with the onset of the generalized HIV epidemic 12 y prior to 2005, localize its geographic origin to a location in northeastern KwaZulu-Natal ∼400 km away from the site of the 2005 outbreak, and use protein structural modeling to propose a mechanism for how strain-specific rpoB mutations offset fitness costs associated with rifampin resistance in LAM4/KZN. Our findings highlight the importance of HIV coinfection, high preexisting rates of drug-resistant TB, human migration, and pathoadaptive evolution in the emergence and dispersal of this critical public health threat. We propose that integrating whole-genome sequencing into routine public health surveillance can enable the early detection and local containment of AMR pathogens before they achieve widespread dispersal.

Truncating variants in the SHANK1 gene are associated with a spectrum of neurodevelopmental disorders.

PURPOSE: In this study, we aimed to characterize the clinical phenotype of a SHANK1-related disorder and define the functional consequences of SHANK1 truncating variants. METHODS: Exome sequencing (ES) was performed for six individuals who presented with neurodevelopmental disorders. Individuals were ascertained with the use of GeneMatcher and Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources (DECIPHER). We evaluated potential nonsense-mediated decay (NMD) of two variants by making knock-in cell lines of endogenous truncated SHANK1, and expressed the truncated SHANK1 complementary DNA (cDNA) in HEK293 cells and cultured hippocampal neurons to examine the proteins. RESULTS: ES detected de novo truncating variants in SHANK1 in six individuals. Evaluation of NMD resulted in stable transcripts, and the truncated SHANK1 completely lost binding with Homer1, a linker protein that binds to the C-terminus of SHANK1. These variants may disrupt protein-protein networks in dendritic spines. Dispersed localization of the truncated SHANK1 variants within the spine and dendritic shaft was also observed when expressed in neurons, indicating impaired synaptic localization of truncated SHANK1. CONCLUSION: This report expands the clinical spectrum of individuals with truncating SHANK1 variants and describes the impact these variants may have on the pathophysiology of neurodevelopmental disorders.

Defining the genotypic and phenotypic spectrum of X-linked MSL3-related disorder.

PURPOSE: We sought to delineate the genotypic and phenotypic spectrum of female and male individuals with X-linked, MSL3-related disorder (Basilicata-Akhtar syndrome). METHODS: Twenty-five individuals (15 males, 10 females) with causative variants in MSL3 were ascertained through exome or genome sequencing at ten different sequencing centers. RESULTS: We identified multiple variant types in MSL3 (ten nonsense, six frameshift, four splice site, three missense, one in-frame-deletion, one multi-exon deletion), most proven to be de novo, and clustering in the terminal eight exons suggesting that truncating variants in the first five exons might be compensated by an alternative MSL3 transcript. Three-dimensional modeling of missense and splice variants indicated that these have a deleterious effect. The main clinical findings comprised developmental delay and intellectual disability ranging from mild to severe. Autism spectrum disorder, muscle tone abnormalities, and macrocephaly were common as well as hearing impairment and gastrointestinal problems. Hypoplasia of the cerebellar vermis emerged as a consistent magnetic resonance image (MRI) finding. Females and males were equally affected. Using facial analysis technology, a recognizable facial gestalt was determined. CONCLUSION: Our aggregated data illustrate the genotypic and phenotypic spectrum of X-linked, MSL3-related disorder (Basilicata-Akhtar syndrome). Our cohort improves the understanding of disease related morbidity and allows us to propose detailed surveillance guidelines for affected individuals.

CSNK2B: A broad spectrum of neurodevelopmental disability and epilepsy severity.

CSNK2B has recently been implicated as a disease gene for neurodevelopmental disability (NDD) and epilepsy. Information about developmental outcomes has been limited by the young age and short follow-up for many of the previously reported cases, and further delineation of the spectrum of associated phenotypes is needed. We present 25 new patients with variants in CSNK2B and refine the associated NDD and epilepsy phenotypes. CSNK2B variants were identified by research or clinical exome sequencing, and investigators from different centers were connected via GeneMatcher. Most individuals had developmental delay and generalized epilepsy with onset in the first 2 years. However, we found a broad spectrum of phenotypic severity, ranging from early normal development with pharmacoresponsive seizures to profound intellectual disability with intractable epilepsy and recurrent refractory status epilepticus. These findings suggest that CSNK2B should be considered in the diagnostic evaluation of patients with a broad range of NDD with treatable or intractable seizures.

Genetic testing in individuals with cerebral palsy.

AIM To determine which patients with cerebral palsy (CP) should undergo genetic testing, we compared the rate of likely causative genetic variants from whole-exome sequencing in individuals with and without environmental risk factors. METHOD Patients were part of a convenience and physician-referred cohort recruited from a single medical center, and research whole-exome sequencing was completed. Participants were evaluated for the following risk factors: extreme preterm birth, brain bleed or stroke, birth asphyxia, brain malformations, and intrauterine infection. RESULTS A total of 151 unrelated individuals with CP (81 females, 70 males; mean age 25y 7mo [SD 17y 5mo], range 3wks-72y) participated. Causative genetic variants were identified in 14 participants (9.3%). There was no significant difference in diagnostic rate between individuals with risk factors (10 out of 123; 8.1%) and those without (4 out of 28; 14.3%) (Fisher's exact p=0.3). INTERPRETATION While the rate of genetic diagnoses among individuals without risk factors was higher than those with risk factors, the difference was not statistically significant at this sample size. The identification of genetic diagnoses in over 8% of cases with risk factors suggests that these might confer susceptibility to environmental factors, and that further research should include individuals with risk factors. What this paper adds There is no significant difference in diagnostic rate between individuals with and without risk factors. Genetic variants may confer susceptibility to environmental risk factors. Six causative variants were identified in genes not previously associated with cerebral palsy. Global developmental delay/intellectual disability is positively associated with a genetic etiology. Extreme preterm birth, stroke/brain hemorrhage, and older age are negatively associated with a genetic etiology.

A pathogenic variant in the SETBP1 hotspot results in a forme-fruste Schinzel-Giedion syndrome.

Schinzel-Giedion syndrome (SGS; OMIM 269150) is an ultra-rare genetic disorder associated with a distinctive facial gestalt, congenital malformations, severe intellectual disability, and a progressive neurological course. The prognosis for SGS is poor, with survival beyond the first decade rare. Germline, de novo heterozygous variants in the SETBP1 gene cause SGS with the pathogenic variants associated with the SGS phenotype missense and confined to exon 4 of the gene, clustered in a four amino acid (12 bp) hotspot in the SKI homologous region of the SETBP1 protein. We report a patient with a de novo I871S variant within the SKI homologous region, which has been associated with the severe phenotype previously; but our patient has fewer features of SGS and a milder course. This is the first report of a forme-fruste phenotype in a patient with a pathogenic variant within the SGS hotspot on the SETBP1 gene and it highlights the importance of considering atypical clinical presentations in the context of severe ultra-rare genetic disorders.

Somatic SLC35A2 variants in the brain are associated with intractable neocortical epilepsy.

OBJECTIVE: Somatic variants are a recognized cause of epilepsy-associated focal malformations of cortical development (MCD). We hypothesized that somatic variants may underlie a wider range of focal epilepsy, including nonlesional focal epilepsy (NLFE). Through genetic analysis of brain tissue, we evaluated the role of somatic variation in focal epilepsy with and without MCD. METHODS: We identified somatic variants through high-depth exome and ultra-high-depth candidate gene sequencing of DNA from epilepsy surgery specimens and leukocytes from 18 individuals with NLFE and 38 with focal MCD. RESULTS: We observed somatic variants in 5 cases in SLC35A2, a gene associated with glycosylation defects and rare X-linked epileptic encephalopathies. Nonsynonymous variants in SLC35A2 were detected in resected brain, and absent from leukocytes, in 3 of 18 individuals (17%) with NLFE, 1 female and 2 males, with variant allele frequencies (VAFs) in brain-derived DNA of 2 to 14%. Pathologic evaluation revealed focal cortical dysplasia type Ia (FCD1a) in 2 of the 3 NLFE cases. In the MCD cohort, nonsynonymous variants in SCL35A2 were detected in the brains of 2 males with intractable epilepsy, developmental delay, and magnetic resonance imaging suggesting FCD, with VAFs of 19 to 53%; Evidence for FCD was not observed in either brain tissue specimen. INTERPRETATION: We report somatic variants in SLC35A2 as an explanation for a substantial fraction of NLFE, a largely unexplained condition, as well as focal MCD, previously shown to result from somatic mutation but until now only in PI3K-AKT-mTOR pathway genes. Collectively, our findings suggest a larger role than previously recognized for glycosylation defects in the intractable epilepsies. Ann Neurol 2018.

NBEA: Developmental disease gene with early generalized epilepsy phenotypes.

NBEA is a candidate gene for autism, and de novo variants have been reported in neurodevelopmental disease (NDD) cohorts. However, NBEA has not been rigorously evaluated as a disease gene, and associated phenotypes have not been delineated. We identified 24 de novo NBEA variants in patients with NDD, establishing NBEA as an NDD gene. Most patients had epilepsy with onset in the first few years of life, often characterized by generalized seizure types, including myoclonic and atonic seizures. Our data show a broader phenotypic spectrum than previously described, including a myoclonic-astatic epilepsy-like phenotype in a subset of patients. Ann Neurol 2018;84:796-803.

Robust classification of protein variation using structural modelling and large-scale data integration.

Existing methods for interpreting protein variation focus on annotating mutation pathogenicity rather than detailed interpretation of variant deleteriousness and frequently use only sequence-based or structure-based information. We present VIPUR, a computational framework that seamlessly integrates sequence analysis and structural modelling (using the Rosetta protein modelling suite) to identify and interpret deleterious protein variants. To train VIPUR, we collected 9477 protein variants with known effects on protein function from multiple organisms and curated structural models for each variant from crystal structures and homology models. VIPUR can be applied to mutations in any organism's proteome with improved generalized accuracy (AUROC .83) and interpretability (AUPR .87) compared to other methods. We demonstrate that VIPUR's predictions of deleteriousness match the biological phenotypes in ClinVar and provide a clear ranking of prediction confidence. We use VIPUR to interpret known mutations associated with inflammation and diabetes, demonstrating the structural diversity of disrupted functional sites and improved interpretation of mutations associated with human diseases. Lastly, we demonstrate VIPUR's ability to highlight candidate variants associated with human diseases by applying VIPUR to de novo variants associated with autism spectrum disorders.

Influence of surface charge, binding site residues and glycosylation on Thielavia terrestris cutinase biochemical characteristics.

Cutinases are esterases of industrial importance for applications in recycling and surface modification of polyesters. The cutinase from Thielavia terrestris (TtC) is distinct in terms of its ability to retain its stability and activity in acidic pH. Stability and activity in acidic pHs are desirable for esterases as the pH of the reaction tends to go down with the generation of acid. The pH stability and activity are governed by the charged state of the residues involved in catalysis or in substrate binding. In this study, we performed the detailed structural and biochemical characterization of TtC coupled with surface charge analysis to understand its acidic tolerance. The stability of TtC in acidic pH was rationalized by evaluating the contribution of charge interactions to the Gibbs free energy of unfolding at varying pHs. The activity of TtC was found to be limited by substrate binding affinity, which is a function of the surface charge. Additionally, the presence of glycosylation affects the biochemical characteristics of TtC owing to steric interactions with residues involved in substrate binding.

Myosin Mutations and Sudden Sensorineural Hearing Loss: Results of Whole Exome Sequencing.

OBJECTIVE: Idiopathic sudden sensorineural hearing loss (ISSNHL) affects 66,000 patients per year in the United States. Genetic mutations have been associated with progressive hearing loss; however, genetic mutations associated with ISSNHL have not been identified. METHODS: A prospective cohort study of adults older than 18 years presenting with ISSNHL at a tertiary academic medical center. Whole exome sequencing (WES) was conducted using Genome Analysis Toolkit best practices. An automated diagnostic screen employing a variety of models for pathogenicity was conducted across all genes with no specific targets. Candidate pathogenic variants were reviewed by a team of geneticists and clinicians. Variants were crossed-referenced with 92 known hearing loss associated genes. RESULTS: Twenty-nine patients with SSNHL were screened using WES. The average age of patients was 53 ± 17.1 years, and most patients were White (62%) and men (55%). The mean pure tone average was 64.8 ± 31.3 dB for the affected ear. Using a 0.1% allele frequency screen, 12 (41%) cases had a mutation in any of the nine selected myosin genes. When we restrict to singletons (allele frequency = 0%), 21% (n = 6) of cases have qualifying variants, whereas only 3.8% (n = 481) of 12,577 healthy controls carry qualifying variants (p < 0.01). Most mutations (80%) were missense mutations. Of the novel mutations, one was a frameshift mutation, and two were a stop-gained function. Three were missense mutations. CONCLUSION: Myosin mutations may be associated with ISSNHL. However, larger population screening is needed to confirm the association of myosin mutation with ISSNHL and steroid responsiveness.

Ancestry adjustment improves genome-wide estimates of regional intolerance.

Genomic regions subject to purifying selection are more likely to carry disease-causing mutations than regions not under selection. Cross species conservation is often used to identify such regions but with limited resolution to detect selection on short evolutionary timescales such as that occurring in only one species. In contrast, genetic intolerance looks for depletion of variation relative to expectation within a species, allowing species-specific features to be identified. When estimating the intolerance of noncoding sequence, methods strongly leverage variant frequency distributions. As the expected distributions depend on ancestry, if not properly controlled for, ancestral population source may obfuscate signals of selection. We demonstrate that properly incorporating ancestry in intolerance estimation greatly improved variant classification. We provide a genome-wide intolerance map that is conditional on ancestry and likely to be particularly valuable for variant prioritization.

Whole exome sequencing reveals potentially pathogenic variants in a small subset of premenopausal women with idiopathic osteoporosis.

Osteoporosis in premenopausal women with intact gonadal function and no known secondary cause of bone loss is termed idiopathic osteoporosis (IOP). Women with IOP diagnosed in adulthood have profound bone structural deficits and often report adult and childhood fractures, and family history of osteoporosis. Some have very low bone formation rates (BFR/BS) suggesting osteoblast dysfunction. These features led us to investigate potential genetic etiologies of bone fragility. In 75 IOP women (aged 20-49) with low trauma fractures and/or very low BMD who had undergone transiliac bone biopsies, we performed Whole Exome Sequencing (WES) using our variant analysis pipeline to select candidate rare and novel variants likely to affect known disease genes. We ran rare-variant burden analyses on all genes individually and on phenotypically-relevant gene sets. For particular genes implicated in osteoporosis, we also assessed the frequency of all (including common) variants in subjects versus 6540 non-comorbid female controls. The variant analysis pipeline identified 4 women with 4 heterozygous variants in LRP5 and PLS3 that were considered to contribute to osteoporosis. All 4 women had adult fractures, and 3 women also had multiple fractures, childhood fractures and a family history of osteoporosis. Two women presented during pregnancy/lactation. In an additional 4 subjects, 4 different relevant Variants of Uncertain Significance (VUS) were detected in the genes FKBP10, SLC34A3, and HGD. Of the subjects with VUS, 2 had multiple adult fractures, childhood fractures, and presented during pregnancy/lactation, and 2 had nephrolithiasis. BFR/BS varied among the 8 subjects with identified variants; BFR/BS was quite low in those with variants that are likely to have adverse effects on bone formation. The analysis pipeline did not discover candidate variants in COL1A1, COL1A2, WNT, or ALPL. Although we found several novel and rare variants in LRP5, cases did not have an increased burden of common LRP5 variants compared to controls. Cohort-wide collapsing analysis did not reveal any novel disease genes with genome-wide significance for qualifying variants between controls and our 75 cases. In summary, WES revealed likely pathogenic variants or relevant VUS in 8 (11%) of 75 women with IOP. Notably, the genetic variants identified were consistent with the affected women's diagnostic evaluations that revealed histological evidence of low BFR/BS or biochemical evidence of increased bone resorption and urinary calcium excretion. These results, and the fact that the majority of the women had no identifiable genetic etiology, also suggest that the pathogenesis of and mechanisms leading to osteoporosis in this cohort are heterogeneous. Future research is necessary to identify both new genetic and non-genetic etiologies of early-onset osteoporosis.

Serine biosynthesis defect due to haploinsufficiency of PHGDH causes retinal disease.

Macular telangiectasia type 2 (MacTel) is a progressive, late-onset retinal degenerative disease linked to decreased serum levels of serine that elevate circulating levels of a toxic ceramide species, deoxysphingolipids (deoxySLs); however, causal genetic variants that reduce serine levels in patients have not been identified. Here we identify rare, functional variants in the gene encoding the rate-limiting serine biosynthetic enzyme, phosphoglycerate dehydrogenase (PHGDH), as the single locus accounting for a significant fraction of MacTel. Under a dominant collapsing analysis model of a genome-wide enrichment analysis of rare variants predicted to impact protein function in 793 MacTel cases and 17,610 matched controls, the PHGDH gene achieves genome-wide significance (P = 1.2 × 10-13) with variants explaining ~3.2% of affected individuals. We further show that the resulting functional defects in PHGDH cause decreased serine biosynthesis and accumulation of deoxySLs in retinal pigmented epithelial cells. PHGDH is a significant locus for MacTel that explains the typical disease phenotype and suggests a number of potential treatment options.

SSBP1 mutations cause mtDNA depletion underlying a complex optic atrophy disorder.

Inherited optic neuropathies include complex phenotypes, mostly driven by mitochondrial dysfunction. We report an optic atrophy spectrum disorder, including retinal macular dystrophy and kidney insufficiency leading to transplantation, associated with mitochondrial DNA (mtDNA) depletion without accumulation of multiple deletions. By whole-exome sequencing, we identified mutations affecting the mitochondrial single-strand binding protein (SSBP1) in 4 families with dominant and 1 with recessive inheritance. We show that SSBP1 mutations in patient-derived fibroblasts variably affect the amount of SSBP1 protein and alter multimer formation, but not the binding to ssDNA. SSBP1 mutations impaired mtDNA, nucleoids, and 7S-DNA amounts as well as mtDNA replication, affecting replisome machinery. The variable mtDNA depletion in cells was reflected in severity of mitochondrial dysfunction, including respiratory efficiency, OXPHOS subunits, and complex amount and assembly. mtDNA depletion and cytochrome c oxidase-negative cells were found ex vivo in biopsies of affected tissues, such as kidney and skeletal muscle. Reduced efficiency of mtDNA replication was also reproduced in vitro, confirming the pathogenic mechanism. Furthermore, ssbp1 suppression in zebrafish induced signs of nephropathy and reduced optic nerve size, the latter phenotype complemented by WT mRNA but not by SSBP1 mutant transcripts. This previously unrecognized disease of mtDNA maintenance implicates SSBP1 mutations as a cause of human pathology.

New insights into tardive dyskinesia genetics: Implementation of whole-exome sequencing approach.

Tardive dyskinesia (TD) is an adverse movement disorder induced by chronic treatment with antipsychotics drugs. The contribution of common genetic variants to TD susceptibility has been investigated in recent years, but with limited success. The aim of the current study was to investigate the potential contribution of rare variants to TD vulnerability. In order to identify TD risk genes, we performed whole-exome sequencing (WES) and gene-based collapsing analysis focusing on rare (allele frequency < 1%) and putatively deleterious variants (qualifying variants). 82 Jewish schizophrenia patients chronically treated with antipsychotics were included and classified as having severe TD or lack of any abnormal movements based on a rigorous definition of the TD phenotype. First, we performed a case-control, exome-wide collapsing analysis comparing 39 schizophrenia patients with severe TD to 3118 unrelated population controls. Then, we checked the potential top candidate genes among 43 patients without any TD manifestations. All the genes that were found to harbor one or more qualifying variants in patients without any TD features were excluded from the final list of candidate genes. Only one gene, regulating synaptic membrane exocytosis 2 (RIMS2), showed significant enrichment of qualifying variants in TD patients compared with unrelated population controls after correcting for multiple testing (Fisher's exact test p = 5.32E-08, logistic regression p = 2.50E-08). Enrichment was caused by a single variant (rs567070433) due to a frameshift in an alternative transcript of RIMS2. None of the TD negative patients had qualifying variants in this gene. In a validation cohort of 140 schizophrenia patients assessed for TD, the variant was also not detected in any individual. Some potentially suggestive TD genes were detected in the TD cohort and warrant follow-up in future studies. No significant enrichment in previously reported TD candidate genes was identified. To the best of our knowledge, this is the first WES study of TD, demonstrating the potential role of rare loss-of-function variant enrichment in this pharmacogenetic phenotype.