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.

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.

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.

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.

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.

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.

Why are there hotspot mutations in the TP53 gene in human cancers?

The p53 gene contains homozygous mutations in ~50-60% of human cancers. About 90% of these mutations encode missense mutant proteins that span ~190 different codons localized in the DNA-binding domain of the gene and protein. These mutations produce a protein with a reduced capacity to bind to a specific DNA sequence that regulates the p53 transcriptional pathway. Eight of these mutations are localized in codons that account for ~28% of the total p53 mutations and these alleles appear to be selected for preferentially in human cancers of many tissue types. This article explores the question 'Why are there hotspot mutations in the p53 gene in human cancers?' Four possible reasons for this are considered; (1) the hotspot mutant alleles produce a protein that has a highly altered structure, (2) environmental mutagens produce allele-specific changes in the p53 gene, (3) these mutations arise at selected sites in the gene due to a specific DNA sequence, such as a methylated cytosine residue in a CpG dinucleotide, which has a higher mutation rate changing C to T nucleotides, (4) along with the observed change in mutant p53 proteins, which produce a loss of function (DNA binding and transcription), some mutant proteins have an allele-specific gain of function that promotes cancer. Evidence is presented that demonstrates the first three possibilities all contribute some property to this list of hotspot mutations. The fourth possibility remains to be tested.

Role of WNT10A in failure of tooth development in humans and zebrafish.

BACKGROUND: Oligodontia is a severe form of tooth agenesis characterized by the absence of six or more permanent teeth. Oligodontia has complex etiology and variations in numerous genes have been suggested as causal for the condition. METHODS: We applied whole-exome sequencing (WES) to identify the cause of oligodontia in a 9-year-old girl missing 11 permanent teeth. Protein modeling and functional analysis in zebrafish were also performed to understand the impact of identified variants on the phenotype. RESULTS: We identified a novel compound heterozygous missense mutation in WNT10A (c.637G>A:p.Gly213Ser and c.1070C>T:p.Thr357Ile) as the likely cause of autosomal recessive oligodontia in the child. Affected residues are located in conserved regions and variants are predicted to be highly deleterious for potentially destabilizing the protein fold and inhibiting normal protein function. Functional studies in zebrafish embryos showed that wnt10a is expressed in the craniofacies at critical time points for tooth development, and that perturbations of wnt10a expression impaired normal tooth development and arrested tooth development at 5 days postfertilization (dpf). Furthermore, mRNA expression levels of additional tooth development genes were directly correlated with wnt10a expression; expression of msx1, dlx2b, eda, and axin2 was decreased upon wnt10a knockdown, and increased upon wnt10a overexpression. CONCLUSIONS: Our results reveal a novel compound heterozygous variant in WNT10A as pathogenic for oligodontia, and demonstrate that perturbations of wnt10a expression in zebrafish may directly and/or indirectly affect tooth development recapitulating the agenesis phenotype observed in humans.

Real-time PyMOL visualization for Rosetta and PyRosetta.

Computational structure prediction and design of proteins and protein-protein complexes have long been inaccessible to those not directly involved in the field. A key missing component has been the ability to visualize the progress of calculations to better understand them. Rosetta is one simulation suite that would benefit from a robust real-time visualization solution. Several tools exist for the sole purpose of visualizing biomolecules; one of the most popular tools, PyMOL (Schrödinger), is a powerful, highly extensible, user friendly, and attractive package. Integrating Rosetta and PyMOL directly has many technical and logistical obstacles inhibiting usage. To circumvent these issues, we developed a novel solution based on transmitting biomolecular structure and energy information via UDP sockets. Rosetta and PyMOL run as separate processes, thereby avoiding many technical obstacles while visualizing information on-demand in real-time. When Rosetta detects changes in the structure of a protein, new coordinates are sent over a UDP network socket to a PyMOL instance running a UDP socket listener. PyMOL then interprets and displays the molecule. This implementation also allows remote execution of Rosetta. When combined with PyRosetta, this visualization solution provides an interactive environment for protein structure prediction and design.