ETV4-Dependent Transcriptional Plasticity Maintains MYC Expression and Results in IMiD Resistance in Multiple Myeloma

Immunomodulatory drugs (IMiD) are a backbone therapy for multiple myeloma (MM). Despite their efficacy, most patients develop resistance, and the mechanisms are not fully defined. Here, we show that IMiD responses are directed by IMiD-dependent degradation of IKZF1 and IKZF3 that bind to enhancers necessary to sustain the expression of MYC and other myeloma oncogenes. IMiD treatment universally depleted chromatin-bound IKZF1, but eviction of P300 and BRD4 coactivators only occurred in IMiD-sensitive cells. IKZF1-bound enhancers overlapped other transcription factor binding motifs, including ETV4. Chromatin immunoprecipitation sequencing showed that ETV4 bound to the same enhancers as IKZF1, and ETV4 CRISPR/Cas9-mediated ablation resulted in sensitization of IMiD-resistant MM. ETV4 expression is associated with IMiD resistance in cell lines, poor prognosis in patients, and is upregulated at relapse. These data indicate that ETV4 alleviates IKZF1 and IKZF3 dependency in MM by maintaining oncogenic enhancer activity and identify transcriptional plasticity as a previously unrecognized mechanism of IMiD resistance.

SIGNIFICANCE

We show that IKZF1-bound enhancers are critical for IMiD efficacy and that the factor ETV4 can bind the same enhancers and substitute for IKZF1 and mediate IMiD resistance by maintaining MYC and other oncogenes. These data implicate transcription factor redundancy as a previously unrecognized mode of IMiD resistance in MM. See related article by Welsh, Barwick, et al., p. 34. See related commentary by Yun and Cleveland, p. 5. This article is featured in Selected Articles from This Issue, p. 4.

The best of both worlds: Blending cutting-edge research with clinical processes for a productive exome clinic

Genomic medicine has been transformed by next-generation sequencing (NGS), inclusive of exome sequencing (ES) and genome sequencing (GS). Currently, ES is offered widely in clinical settings, with a less prevalent alternative model consisting of hybrid programs that incorporate research ES along with clinical patient workflows. We were among the earliest to implement a hybrid ES clinic, have provided diagnoses to 45% of probands, and have identified several novel candidate genes. Our program is enabled by a cost-effective investment by the health system and is unique in encompassing all the processes that have been variably included in other hybrid/clinical programs. These include careful patient selection, utilization of a phenotype-agnostic bioinformatics pipeline followed by manual curation of variants and phenotype integration by clinicians, close collaborations between the clinicians and the bioinformatician, pursuit of interesting variants, communication of results to patients in categories that are predicated upon the certainty of a diagnosis, and tracking changes in results over time and the underlying mechanisms for such changes. Due to its effectiveness, scalability to GS and its resource efficiency, specific elements of our paradigm can be incorporated into existing clinical settings, or the entire hybrid model can be implemented within health systems that have genomic medicine programs, to provide NGS in a scientifically rigorous, yet pragmatic setting.

Multiomic Analysis Identifies a High-Risk Metabolic and TME Depleted Signature that Predicts Early Clinical Failure in DLBCL

PURPOSE

60-70% of newly diagnosed diffuse large B-cell lymphoma (DLBCL) patients avoid events within 24 months of diagnosis (EFS24) and the remainder have poor outcomes. Recent genetic and molecular classification of DLBCL has advanced our knowledge of disease biology, yet were not designed to predict early events and guide anticipatory selection of novel therapies. To address this unmet need, we used an integrative multiomic approach to identify a signature at diagnosis that will identify DLBCL at high risk of early clinical failure.

PATIENTS AND METHODS

Tumor biopsies from 444 newly diagnosed DLBCL were analyzed by WES and RNAseq. A combination of weighted gene correlation network analysis and differential gene expression analysis followed by integration with clinical and genomic data was used to identify a multiomic signature associated with high risk of early clinical failure.

RESULTS

Current DLBCL classifiers are unable to discriminate cases who fail EFS24. We identified a high risk RNA signature that had a hazard ratio (HR, 18.46 [95% CI 6.51-52.31] P < .001) in a univariate model, which did not attenuate after adjustment for age, IPI and COO (HR, 20.8 [95% CI, 7.14-61.09] P < .001). Further analysis revealed the signature was associated with metabolic reprogramming and a depleted immune microenvironment. Finally, WES data was integrated into the signature and we found that inclusion of ARID1A mutations resulted in identification of 45% of cases with an early clinical failure which was validated in external DLBCL cohorts.

CONCLUSION

This novel and integrative approach is the first to identify a signature at diagnosis that will identify DLBCL at high risk for early clinical failure and may have significant implications for design of therapeutic options.

Whole-genome analysis identifies novel drivers and high-risk double-hit events in relapsed/refractory myeloma

Large-scale analyses of genomic data from patients with newly diagnosed multiple myeloma (ndMM) have been undertaken, however, large-scale analysis of relapsed/refractory MM (rrMM) has not been performed. We hypothesize that somatic variants chronicle the therapeutic exposures and clonal structure of myeloma from ndMM to rrMM stages. We generated whole-genome sequencing (WGS) data from 418 tumors (386 patients) derived from 6 rrMM clinical trials and compared them with WGS from 198 unrelated patients with ndMM in a population-based case-control fashion. We identified significantly enriched events at the rrMM stage, including drivers (DUOX2, EZH2, TP53), biallelic inactivation (TP53), noncoding mutations in bona fide drivers (TP53BP1, BLM), copy number aberrations (CNAs; 1qGain, 17pLOH), and double-hit events (Amp1q-ISS3, 1qGain-17p loss-of-heterozygosity). Mutational signature analysis identified a subclonal defective mismatch repair signature enriched in rrMM and highly active in high mutation burden tumors, a likely feature of therapy-associated expanding subclones. Further analysis focused on the association of genomic aberrations enriched at different stages of resistance to immunomodulatory agent (IMiD)-based therapy. This analysis revealed that TP53, DUOX2, 1qGain, and 17p loss-of-heterozygosity increased in prevalence from ndMM to lenalidomide resistant (LENR) to pomalidomide resistant (POMR) stages, whereas enrichment of MAML3 along with immunoglobulin lambda (IGL) and MYC translocations distinguished POM from the LEN subgroup. Genomic drivers associated with rrMM are those that confer clonal selective advantage under therapeutic pressure. Their role in therapy evasion should be further evaluated in longitudinal patient samples, to confirm these associations with the evolution of clinical resistance and to identify molecular subsets of rrMM for the development of targeted therapies.

The microRNA processor DROSHA is a candidate gene for a severe progressive neurological disorder

DROSHA encodes a ribonuclease that is a subunit of the Microprocessor complex and is involved in the first step of microRNA (miRNA) biogenesis. To date, DROSHA has not yet been associated with a Mendelian disease. Here, we describe two individuals with profound intellectual disability, epilepsy, white matter atrophy, microcephaly and dysmorphic features, who carry damaging de novo heterozygous variants in DROSHA. DROSHA is constrained for missense variants and moderately intolerant to loss-of-function (o/e = 0.24). The loss of the fruit fly ortholog drosha causes developmental arrest and death in third instar larvae, a severe reduction in brain size and loss of imaginal discs in the larva. Loss of drosha in eye clones causes small and rough eyes in adult flies. One of the identified DROSHA variants (p.Asp1219Gly) behaves as a strong loss-of-function allele in flies, while another variant (p.Arg1342Trp) is less damaging in our assays. In worms, a knock-in that mimics the p.Asp1219Gly variant at a worm equivalent residue causes loss of miRNA expression and heterochronicity, a phenotype characteristic of the loss of miRNA. Together, our data show that the DROSHA variants found in the individuals presented here are damaging based on functional studies in model organisms and likely underlie the severe phenotype involving the nervous system.

Expanding the phenotypic spectrum of ARCN1-related syndrome

PURPOSE

This study aimed to describe the phenotypic and molecular characteristics of ARCN1-related syndrome.

METHODS

Patients with ARCN1 variants were identified, and clinician researchers were connected using GeneMatcher and physician referrals. Clinical histories were collected from each patient.

RESULTS

In total, we identified 14 cases of ARCN1-related syndrome, (9 pediatrics, and 5 fetal cases from 3 families). The clinical features these newly identified cases were compared to 6 previously reported cases for a total of 20 cases. Intrauterine growth restriction, micrognathia, and short stature were present in all patients. Other common features included prematurity (11/15, 73.3%), developmental delay (10/14, 71.4%), genitourinary malformations in males (6/8, 75%), and microcephaly (12/15, 80%). Novel features of ARCN1-related syndrome included transient liver dysfunction and specific glycosylation abnormalities during illness, giant cell hepatitis, hepatoblastoma, cataracts, and lethal skeletal manifestations. Developmental delay was seen in 73% of patients, but only 3 patients had intellectual disability, which is less common than previously reported.

CONCLUSION

ARCN1-related syndrome presents with a wide clinical spectrum ranging from a severe embryonic lethal syndrome to a mild syndrome with intrauterine growth restriction, micrognathia, and short stature without intellectual disability. Patients with ARCN1-related syndrome should be monitored for liver dysfunction during illness, cataracts, and hepatoblastoma. Additional research to further define the phenotypic spectrum and possible genotype-phenotype correlations are required.

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.

Integrative multi-omics identifies high risk multiple myeloma subgroup associated with significant DNA loss and dysregulated DNA repair and cell cycle pathways

Background

Despite significant therapeutic advances in improving lives of multiple myeloma (MM) patients, it remains mostly incurable, with patients ultimately becoming refractory to therapies. MM is a genetically heterogeneous disease and therapeutic resistance is driven by a complex interplay of disease pathobiology and mechanisms of drug resistance. We applied a multi-omics strategy using tumor-derived gene expression, single nucleotide variant, copy number variant, and structural variant profiles to investigate molecular subgroups in 514 newly diagnosed MM (NDMM) samples and identified 12 molecularly defined MM subgroups (MDMS1-12) with distinct genomic and transcriptomic features.

Results

Our integrative approach let us identify NDMM subgroups with transversal profiles to previously described ones, based on single data types, which shows the impact of this approach for disease stratification. One key novel subgroup is our MDMS8, associated with poor clinical outcome [median overall survival, 38 months (global log-rank p-value < 1 × 10⁻⁶)], which uniquely presents a broad genomic loss (> 9% of entire genome, t-test p value < 1e−5) driving dysregulation of various transcriptional programs affecting DNA repair and cell cycle/mitotic processes. This subgroup was validated on multiple independent datasets, and a master regulator analyses identified transcription factors controlling MDMS8 transcriptomic profile, including CKS1B and PRKDC among others, which are regulators of the DNA repair and cell cycle pathways.

Conclusion

Using multi-omics unsupervised clustering we were able to discover a new high-risk multiple myeloma patient segment. This high-risk group presents diverse previously known genetic markers, but also a new characteristic defined by accumulation of genomic loss which seems to drive transcriptional dysregulation of cell cycle, DNA repair and DNA damage. Finally, our work identified various master regulators, including E2F2 and CKS1B as the genes controlling these key biological pathways.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01140-5.

Rare deleterious de novo missense variants in Rnf2/Ring2 are associated with a neurodevelopmental disorder with unique clinical features

The Polycomb group (PcG) gene RNF2 (RING2) encodes a catalytic subunit of the Polycomb repressive complex 1 (PRC1), an evolutionarily conserved machinery that post-translationally modifies chromatin to maintain epigenetic transcriptional repressive states of target genes including Hox genes. Here, we describe two individuals, each with rare de novo missense variants in RNF2. Their phenotypes include intrauterine growth retardation, severe intellectual disabilities, behavioral problems, seizures, feeding difficulties and dysmorphic features. Population genomics data suggest that RNF2 is highly constrained for loss-of-function (LoF) and missense variants, and both p.R70H and p.S82R variants have not been reported to date. Structural analyses of the two alleles indicate that these changes likely impact the interaction between RNF2 and BMI1, another PRC1 subunit or its substrate Histone H2A, respectively. Finally, we provide functional data in Drosophila that these two missense variants behave as LoF alleles in vivo. The evidence provide support for deleterious alleles in RNF2 being associated with a new and recognizable genetic disorder. This tentative gene-disease association in addition to the 12 previously identified disorders caused by PcG genes attests to the importance of these chromatin regulators in Mendelian disorders.

TSPEAR variants are primarily associated with ectodermal dysplasia and tooth agenesis but not hearing loss: A novel cohort study

Biallelic loss‐of‐function variants in the thrombospondin‐type laminin G domain and epilepsy‐associated repeats (TSPEAR) gene have recently been associated with ectodermal dysplasia and hearing loss. The first reports describing a TSPEAR disease association identified this gene is a cause of nonsyndromic hearing loss, but subsequent reports involving additional affected families have questioned this evidence and suggested a stronger association with ectodermal dysplasia. To clarify genotype–phenotype associations for TSPEAR variants, we characterized 13 individuals with biallelic TSPEAR variants. Individuals underwent either exome sequencing or panel‐based genetic testing. Nearly all of these newly reported individuals (11/13) have phenotypes that include tooth agenesis or ectodermal dysplasia, while three newly reported individuals have hearing loss. Of the individuals displaying hearing loss, all have additional variants in other hearing‐loss‐associated genes, specifically TMPRSS3, GJB2, and GJB6, that present competing candidates for their hearing loss phenotype. When presented alongside previous reports, the overall evidence supports the association of TSPEAR variants with ectodermal dysplasia and tooth agenesis features but creates significant doubt as to whether TSPEAR variants are a monogenic cause of hearing loss. Further functional evidence is needed to evaluate this phenotypic association.

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.

De novo variants in TCF7L2 are associated with a syndromic neurodevelopmental disorder

TCF7L2 encodes transcription factor 7-like 2 (OMIM 602228), a key mediator of the evolutionary conserved canonical Wnt signaling pathway. Although several large-scale sequencing studies have implicated TCF7L2 in intellectual disability and autism, both the genetic mechanism and clinical phenotype have remained incompletely characterized. We present here a comprehensive genetic and phenotypic description of 11 individuals who have been identified to carry de novo variants in TCF7L2, both truncating and missense. Missense variation is clustered in or near a high mobility group box domain, involving this region in these variants' pathogenicity. All affected individuals present with developmental delays in childhood, but most ultimately achieved normal intelligence or had only mild intellectual disability. Myopia was present in approximately half of the individuals, and some individuals also possessed dysmorphic craniofacial features, orthopedic abnormalities, or neuropsychiatric comorbidities including autism and attention-deficit/hyperactivity disorder (ADHD). We thus present an initial clinical and genotypic spectrum associated with variation in TCF7L2, which will be important in informing both medical management and future research.

Clinical sites of the Undiagnosed Diseases Network: unique contributions to genomic medicine and science

PURPOSE

The NIH Undiagnosed Diseases Network (UDN) evaluates participants with disorders that have defied diagnosis, applying personalized clinical and genomic evaluations and innovative research. The clinical sites of the UDN are essential to advancing the UDN mission; this study assesses their contributions relative to standard clinical practices.

METHODS

We analyzed retrospective data from four UDN clinical sites, from July 2015 to September 2019, for diagnoses, new disease gene discoveries and the underlying investigative methods.

RESULTS

Of 791 evaluated individuals, 231 received 240 diagnoses and 17 new disease-gene associations were recognized. Straightforward diagnoses on UDN exome and genome sequencing occurred in 35% (84/240). We considered these tractable in standard clinical practice, although genome sequencing is not yet widely available clinically. The majority (156/240, 65%) required additional UDN-driven investigations, including 90 diagnoses that occurred after prior nondiagnostic exome sequencing and 45 diagnoses (19%) that were nongenetic. The UDN-driven investigations included complementary/supplementary phenotyping, innovative analyses of genomic variants, and collaborative science for functional assays and animal modeling.

CONCLUSION

Investigations driven by the clinical sites identified diagnostic and research paradigms that surpass standard diagnostic processes. The new diagnoses, disease gene discoveries, and delineation of novel disorders represent a model for genomic medicine and science.

De novo variants in SNAP25 cause an early-onset developmental and epileptic encephalopathy

PURPOSE

This study aims to provide a comprehensive description of the phenotypic and genotypic spectrum of SNAP25 developmental and epileptic encephalopathy (SNAP25-DEE) by reviewing newly identified and previously reported individuals.

METHODS

Individuals harboring heterozygous missense or loss-of-function variants in SNAP25 were assembled through collaboration with international colleagues, matchmaking platforms, and literature review. For each individual, detailed phenotyping, classification, and structural modeling of the identified variant were performed.

RESULTS

The cohort comprises 23 individuals with pathogenic or likely pathogenic de novo variants in SNAP25. Intellectual disability and early-onset epilepsy were identified as the core symptoms of SNAP25-DEE, with recurrent findings of movement disorders, cerebral visual impairment, and brain atrophy. Structural modeling for all variants predicted possible functional defects concerning SNAP25 or impaired interaction with other components of the SNARE complex.

CONCLUSION

We provide a comprehensive description of SNAP25-DEE with intellectual disability and early-onset epilepsy mostly occurring before the age of two years. These core symptoms and additional recurrent phenotypes show an overlap to genes encoding other components or associated proteins of the SNARE complex such as STX1B, STXBP1, or VAMP2. Thus, these findings advance the concept of a group of neurodevelopmental disorders that may be termed "SNAREopathies."

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.

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.

Alternative transcripts in variant interpretation: the potential for missed diagnoses and misdiagnoses

PURPOSE

Guidelines by professional organizations for assessing variant pathogenicity include the recommendation to utilize biologically relevant transcripts; however, there is variability in transcript selection by laboratories.

METHODS

We describe three patients whose genomic results were incorrect, because alternative transcripts and tissue expression patterns were not considered by the commercial laboratories.

RESULTS

In individual 1, a pathogenic coding variant in a brain-expressed isoform of CKDL5 was missed twice on sequencing, because the variant was intronic in the transcripts considered in analysis. In individual 2, a microdeletion affecting KMT2C was not reported on microarray, since deletions of proximal exons in this gene are seen in healthy individuals; however, this individual had a more distal deletion involving the brain-expressed KMT2C isoform, giving her a diagnosis of Kleefstra syndrome. Individual 3 was reported to have a pathogenic variant in exon 10 of OFD1 on exome, but had no typical features of the OFD1-related disorders. Since exon 10 is spliced from the more biologically relevant transcripts of OFD1, it was determined that he did not have an OFD1 disorder.

CONCLUSION

These examples illustrate the importance of considering alternative transcripts as a potential confounder when genetic results are negative or discordant with the phenotype.

De novo EIF2AK1 and EIF2AK2 Variants Are Associated with Developmental Delay, Leukoencephalopathy, and Neurologic Decompensation

EIF2AK1 and EIF2AK2 encode members of the eukaryotic translation initiation factor 2 alpha kinase (EIF2AK) family that inhibits protein synthesis in response to physiologic stress conditions. EIF2AK2 is also involved in innate immune response and the regulation of signal transduction, apoptosis, cell proliferation, and differentiation. Despite these findings, human disorders associated with deleterious variants in EIF2AK1 and EIF2AK2 have not been reported. Here, we describe the identification of nine unrelated individuals with heterozygous de novo missense variants in EIF2AK1 (1/9) or EIF2AK2 (8/9). Features seen in these nine individuals include white matter alterations (9/9), developmental delay (9/9), impaired language (9/9), cognitive impairment (8/9), ataxia (6/9), dysarthria in probands with verbal ability (6/9), hypotonia (7/9), hypertonia (6/9), and involuntary movements (3/9). Individuals with EIF2AK2 variants also exhibit neurological regression in the setting of febrile illness or infection. We use mammalian cell lines and proband-derived fibroblasts to further confirm the pathogenicity of variants in these genes and found reduced kinase activity. EIF2AKs phosphorylate eukaryotic translation initiation factor 2 subunit 1 (EIF2S1, also known as EIF2α), which then inhibits EIF2B activity. Deleterious variants in genes encoding EIF2B proteins cause childhood ataxia with central nervous system hypomyelination/vanishing white matter (CACH/VWM), a leukodystrophy characterized by neurologic regression in the setting of febrile illness and other stressors. Our findings indicate that EIF2AK2 missense variants cause a neurodevelopmental syndrome that may share phenotypic and pathogenic mechanisms with CACH/VWM.

Modification of the PROM1 disease phenotype by a mutation in ABCA4

Background: The extensive phenotypic heterogeneity of monogenic diseases can be largely traced to intragenic variation; however, recent advances in clinical detection and gene sequencing have uncovered the emerging role of non-allelic variation (i.e. genetic trans-modifiers) in shaping disease phenotypes. Identifying these associations are not only of significant diagnostic value, but also provides scientific insight into the expanded molecular etiology of rare diseases. This reports describes the discordant clinical manifestation of a family segregating mutations in ABCA4 and PROM1. Methods: Three patients across a two generation family underwent multimodal imaging and functional testing of the retina including color photography, fundus autofluorescence (AF), spectral domain-optical coherence tomography (SD-OCT) and full-field electroretinography (ffERG). Genetic characterization was carried out by direct Sanger and whole exome sequencing. Results: Clinical examination revealed similar retinal degenerative phenotypes in the proband and her mother. Despite being younger, the proband's phenotype was more advanced and exhibited additional features related to Stargardt disease not found in the mother. Whole exome sequencing identified a pathogenic missense variant in PROM1, c.400C > T, p.(Arg134Cys), as the underlying cause of retinal disease in both the proband and mother. Sequencing of the ABCA4 locus uncovered a single disease-causing variant, c.5714 + 5G > A in the daughter segregating from the father who, surprisingly, also exhibited very subtle disease changes associated with STGD1 despite being a heterozygous carrier. Conclusions: Harboring an additional heterozygous ABCA4 mutation increases severity and confers STGD1-like features in patients with PROM1 disease which provides supporting evidence for their shared pathophysiology and potential treatment prospects.

Heterozygous variants in MYBPC1 are associated with an expanded neuromuscular phenotype beyond arthrogryposis

Encoding the slow skeletal muscle isoform of myosin binding protein-C, MYBPC1 is associated with autosomal dominant and recessive forms of arthrogryposis. The authors describe a novel association for MYBPC1 in four patients from three independent families with skeletal muscle weakness, myogenic tremors, and hypotonia with gradual clinical improvement. The patients carried one of two de novo heterozygous variants in MYBPC1, with the p.Leu263Arg variant seen in three individuals and the p.Leu259Pro variant in one individual. Both variants are absent from controls, well conserved across vertebrate species, predicted to be damaging, and located in the M-motif. Protein modeling studies suggested that the p.Leu263Arg variant affects the stability of the M-motif, whereas the p.Leu259Pro variant alters its structure. In vitro biochemical and kinetic studies demonstrated that the p.Leu263Arg variant results in decreased binding of the M-motif to myosin, which likely impairs the formation of actomyosin cross-bridges during muscle contraction. Collectively, our data substantiate that damaging variants in MYBPC1 are associated with a new form of an early-onset myopathy with tremor, which is a defining and consistent characteristic in all affected individuals, with no contractures. Recognition of this expanded myopathic phenotype can enable identification of individuals with MYBPC1 variants without arthrogryposis.

Autism and developmental disability caused by KCNQ3 gain-of-function variants

OBJECTIVE

Recent reports have described single individuals with neurodevelopmental disability (NDD) harboring heterozygous KCNQ3 de novo variants (DNVs). We sought to assess whether pathogenic variants in KCNQ3 cause NDD and to elucidate the associated phenotype and molecular mechanisms.

METHODS

Patients with NDD and KCNQ3 DNVs were identified through an international collaboration. Phenotypes were characterized by clinical assessment, review of charts, electroencephalographic (EEG) recordings, and parental interview. Functional consequences of variants were analyzed in vitro by patch-clamp recording.

RESULTS

Eleven patients were assessed. They had recurrent heterozygous DNVs in KCNQ3 affecting residues R230 (R230C, R230H, R230S) and R227 (R227Q). All patients exhibited global developmental delay within the first 2 years of life. Most (8/11, 73%) were nonverbal or had a few words only. All patients had autistic features, and autism spectrum disorder (ASD) was diagnosed in 5 of 11 (45%). EEGs performed before 10 years of age revealed frequent sleep-activated multifocal epileptiform discharges in 8 of 11 (73%). For 6 of 9 (67%) recorded between 1.5 and 6 years of age, spikes became near-continuous during sleep. Interestingly, most patients (9/11, 82%) did not have seizures, and no patient had seizures in the neonatal period. Voltage-clamp recordings of the mutant KCNQ3 channels revealed gain-of-function (GoF) effects.

INTERPRETATION

Specific GoF variants in KCNQ3 cause NDD, ASD, and abundant sleep-activated spikes. This new phenotype contrasts both with self-limited neonatal epilepsy due to KCNQ3 partial loss of function, and with the neonatal or infantile onset epileptic encephalopathies due to KCNQ2 GoF. ANN NEUROL 2019;86:181-192.

Late-onset pattern macular dystrophy mimicking ABCA4 and PRPH2 disease is caused by a homozygous frameshift mutation in ROM1

ROM1 (retinal outer segment membrane protein 1) is a 351-amino acid integral membrane protein on Chromosome 11q, with high structural similarity to PRPH2/RDS. Localized at the rims of photoreceptor outer segments (OSs), it is required for the maintenance of OS structure. Here, we describe a case with a phenotypic manifestation of a homozygous single-base pair deletion, c.712delC (p.Leu238Cysfs*78) in the ROM1 gene, resulting in early termination at exon 2. The variant was detected by whole-exome sequencing (WES) in a 63-yr-old Caucasian woman with late-onset pattern macular dystrophy. Notably, although the phenotype resembles those caused by pathogenic variants in ABCA4 or RDS/PRPH2, no pathogenic variants in these, or any other plausible candidate genes, were identified by WES. Clinical features include the presence of hyperautofluorescent flecks, relative sparing of the central macula, and preserved visual acuity. Reduced visual sensitivity was detected among flecked regions in the retina; however, full-field electroretinogram testing revealed no generalized cone dysfunction. The described first case of the complete loss of ROM1 protein function in the retina suggests its sufficiency for late-onset macular dystrophy. ROM1 and PRPH2 pattern macular dystrophies exhibit phenotype overlap, which may be attributable to their shared role in maintenance of the photoreceptor outer segment structure.

Missense Variants in the Histone Acetyltransferase Complex Component Gene TRRAP Cause Autism and Syndromic Intellectual Disability

Acetylation of the lysine residues in histones and other DNA-binding proteins plays a major role in regulation of eukaryotic gene expression. This process is controlled by histone acetyltransferases (HATs/KATs) found in multiprotein complexes that are recruited to chromatin by the scaffolding subunit transformation/transcription domain-associated protein (TRRAP). TRRAP is evolutionarily conserved and is among the top five genes intolerant to missense variation. Through an international collaboration, 17 distinct de novo or apparently de novo variants were identified in TRRAP in 24 individuals. A strong genotype-phenotype correlation was observed with two distinct clinical spectra. The first is a complex, multi-systemic syndrome associated with various malformations of the brain, heart, kidneys, and genitourinary system and characterized by a wide range of intellectual functioning; a number of affected individuals have intellectual disability (ID) and markedly impaired basic life functions. Individuals with this phenotype had missense variants clustering around the c.3127G>A p.(Ala1043Thr) variant identified in five individuals. The second spectrum manifested with autism spectrum disorder (ASD) and/or ID and epilepsy. Facial dysmorphism was seen in both groups and included upslanted palpebral fissures, epicanthus, telecanthus, a wide nasal bridge and ridge, a broad and smooth philtrum, and a thin upper lip. RNA sequencing analysis of skin fibroblasts derived from affected individuals skin fibroblasts showed significant changes in the expression of several genes implicated in neuronal function and ion transport. Thus, we describe here the clinical spectrum associated with TRRAP pathogenic missense variants, and we suggest a genotype-phenotype correlation useful for clinical evaluation of the pathogenicity of the variants.

Whole-exome sequencing in the evaluation of fetal structural anomalies: a prospective cohort study

BACKGROUND

Identification of chromosomal aneuploidies and copy number variants that are associated with fetal structural anomalies has substantial value. Although whole-exome sequencing (WES) has been applied to case series of a few selected prenatal cases, its value in routine clinical settings has not been prospectively assessed in a large unselected cohort of fetuses with structural anomalies. We therefore aimed to determine the incremental diagnostic yield (ie, the added value) of WES following uninformative results of standard investigations with karyotype testing and chromosomal microarray in an unselected cohort of sequential pregnancies showing fetal structural anomalies.

METHODS

In this prospective cohort study, the parents of fetuses who were found to have a structural anomaly in a prenatal ultrasound were screened for possible participation in the study. These participants were predominantly identified in or were referred to the Columbia University Carmen and John Thain Center for Prenatal Pediatrics (New York, NY, USA). Fetuses with confirmed aneuploidy or a causal pathogenic copy number variant were excluded from WES analyses. By use of WES of the fetuses and parents (parent-fetus trios), we identified genetic variants that indicated an underlying cause (diagnostic genetic variants) and genetic variants that met the criteria of bioinformatic signatures that had previously been described to be significantly enriched among diagnostic genetic variants.

FINDINGS

Between April 24, 2015, and April 19, 2017, 517 sequentially identified pregnant women found to have fetuses with a structural anomaly were screened for their eligibility for inclusion in our study. 71 (14%) couples declined testing, 87 (17%) trios were missing at least one DNA sample (from either parent or the fetus), 69 (13%) trios had a clinically relevant abnormal karyotype or chromosomal microarray finding, 51 (10%) couples did not consent to WES or withdrew consent, and five (1%) samples were not of good enough quality for analysis. DNA samples from 234 (45%) eligible trios were therefore used for analysis of the primary outcome. By use of trio sequence data, we identified diagnostic genetic variants in 24 (10%) families. Mutations with bioinformatic signatures that were indicative of pathogenicity but with insufficient evidence to be considered diagnostic were also evaluated; 46 (20%) of the 234 fetuses assessed were found to have such signatures.

INTERPRETATION

Our analysis of WES data in a prospective cohort of unselected fetuses with structural anomalies shows the value added by WES following the use of routine genetic tests. Our findings suggest that, in cases of fetal anomalies in which assessment with karyotype testing and chromosomal microarray fail to determine the underlying cause of a structural anomaly, WES can add clinically relevant information that could assist current management of a pregnancy. The unique challenges of WES-based prenatal diagnostics require analysis by a multidisciplinary team of perinatal practitioners and laboratory specialists.

FUNDING

Institute for Genomic Medicine (Columbia University Irving Medical Center).

Improved Pathogenic Variant Localization via a Hierarchical Model of Sub-regional Intolerance

Different parts of a gene can be of differential importance to development and health. This regional heterogeneity is also apparent in the distribution of disease-associated mutations, which often cluster in particular regions of disease-associated genes. The ability to precisely estimate functionally important sub-regions of genes will be key in correctly deciphering relationships between genetic variation and disease. Previous methods have had some success using standing human variation to characterize this variability in importance by measuring sub-regional intolerance, i.e., the depletion in functional variation from expectation within a given region of a gene. However, the ability to precisely estimate local intolerance was restricted by the fact that only information within a given sub-region is used, leading to instability in local estimates, especially for small regions. We show that borrowing information across regions using a Bayesian hierarchical model stabilizes estimates, leading to lower variability and improved predictive utility. Specifically, our approach more effectively identifies regions enriched for ClinVar pathogenic variants. We also identify significant correlations between sub-region intolerance and the distribution of pathogenic variation in disease-associated genes, with AUCs for classifying de novo missense variants in Online Mendelian Inheritance in Man (OMIM) genes of up to 0.86 using exonic sub-regions and 0.91 using sub-regions defined by protein domains. This result immediately suggests that considering the intolerance of regions in which variants are found may improve diagnostic interpretation. We also illustrate the utility of integrating regional intolerance into gene-level disease association tests with a study of known disease-associated genes for epileptic encephalopathy.

Expanding the Spectrum of BAF-Related Disorders: De Novo Variants in SMARCC2 Cause a Syndrome with Intellectual Disability and Developmental Delay

SMARCC2 (BAF170) is one of the invariable core subunits of the ATP-dependent chromatin remodeling BAF (BRG1-associated factor) complex and plays a crucial role in embryogenesis and corticogenesis. Pathogenic variants in genes encoding other components of the BAF complex have been associated with intellectual disability syndromes. Despite its significant biological role, variants in SMARCC2 have not been directly associated with human disease previously. Using whole-exome sequencing and a web-based gene-matching program, we identified 15 individuals with variable degrees of neurodevelopmental delay and growth retardation harboring one of 13 heterozygous variants in SMARCC2, most of them novel and proven de novo. The clinical presentation overlaps with intellectual disability syndromes associated with other BAF subunits, such as Coffin-Siris and Nicolaides-Baraitser syndromes and includes prominent speech impairment, hypotonia, feeding difficulties, behavioral abnormalities, and dysmorphic features such as hypertrichosis, thick eyebrows, thin upper lip vermilion, and upturned nose. Nine out of the fifteen individuals harbor variants in the highly conserved SMARCC2 DNA-interacting domains (SANT and SWIRM) and present with a more severe phenotype. Two of these individuals present cardiac abnormalities. Transcriptomic analysis of fibroblasts from affected individuals highlights a group of differentially expressed genes with possible roles in regulation of neuronal development and function, namely H19, SCRG1, RELN, and CACNB4. Our findings suggest a novel SMARCC2-related syndrome that overlaps with neurodevelopmental disorders associated with variants in BAF-complex subunits.

A comprehensive iterative approach is highly effective in diagnosing individuals who are exome negative

PURPOSE

Sixty to seventy-five percent of individuals with rare and undiagnosed phenotypes remain undiagnosed after exome sequencing (ES). With standard ES reanalysis resolving 10-15% of the ES negatives, further approaches are necessary to maximize diagnoses in these individuals.

METHODS

In 38 ES negative patients an individualized genomic-phenotypic approach was employed utilizing (1) phenotyping; (2) reanalyses of FASTQ files, with innovative bioinformatics; (3) targeted molecular testing; (4) genome sequencing (GS); and (5) conferring of clinical diagnoses when pathognomonic clinical findings occurred.

RESULTS

Certain and highly likely diagnoses were made in 18/38 (47%) individuals, including identifying two new developmental disorders. The majority of diagnoses (>70%) were due to our bioinformatics, phenotyping, and targeted testing identifying variants that were undetected or not prioritized on prior ES. GS diagnosed 3/18 individuals with structural variants not amenable to ES. Additionally, tentative diagnoses were made in 3 (8%), and in 5 individuals (13%) candidate genes were identified. Overall, diagnoses/potential leads were identified in 26/38 (68%).

CONCLUSIONS

Our comprehensive approach to ES negatives maximizes the ES and clinical data for both diagnoses and candidate gene identification, without GS in the majority. This iterative approach is cost-effective and is pertinent to the current conundrum of ES negatives.

Loss of tubulin deglutamylase CCP1 causes infantile-onset neurodegeneration

A set of glutamylases and deglutamylases controls levels of tubulin polyglutamylation, a prominent post-translational modification of neuronal microtubules. Defective tubulin polyglutamylation was first linked to neurodegeneration in the Purkinje cell degeneration (pcd) mouse, which lacks deglutamylase CCP1, displays massive cerebellar atrophy, and accumulates abnormally glutamylated tubulin in degenerating neurons. We found biallelic rare and damaging variants in the gene encoding CCP1 in 13 individuals with infantile-onset neurodegeneration and confirmed the absence of functional CCP1 along with dysregulated tubulin polyglutamylation. The human disease mainly affected the cerebellum, spinal motor neurons, and peripheral nerves. We also demonstrate previously unrecognized peripheral nerve and spinal motor neuron degeneration in pcd mice, which thus recapitulated key features of the human disease. Our findings link human neurodegeneration to tubulin polyglutamylation, entailing this post-translational modification as a potential target for drug development for neurodegenerative disorders.

De novo mutations in the GTP/GDP-binding region of RALA, a RAS-like small GTPase, cause intellectual disability and developmental delay

Mutations that alter signaling of RAS/MAPK-family proteins give rise to a group of Mendelian diseases known as RASopathies. However, among RASopathies, the matrix of genotype-phenotype relationships is still incomplete, in part because there are many RAS-related proteins and in part because the phenotypic consequences may be variable and/or pleiotropic. Here, we describe a cohort of ten cases, drawn from six clinical sites and over 16,000 sequenced probands, with de novo protein-altering variation in RALA, a RAS-like small GTPase. All probands present with speech and motor delays, and most have intellectual disability, low weight, short stature, and facial dysmorphism. The observed rate of de novo RALA variants in affected probands is significantly higher (p = 4.93 x 10⁻¹¹) than expected from the estimated random mutation rate. Further, all de novo variants described here affect residues within the GTP/GDP-binding region of RALA; in fact, six alleles arose at only two codons, Val25 and Lys128. The affected residues are highly conserved across both RAL- and RAS-family genes, are devoid of variation in large human population datasets, and several are homologous to positions at which disease-associated variants have been observed in other GTPase genes. We directly assayed GTP hydrolysis and RALA effector-protein binding of the observed variants, and found that all but one tested variant significantly reduced both activities compared to wild-type. The one exception, S157A, reduced GTP hydrolysis but significantly increased RALA-effector binding, an observation similar to that seen for oncogenic RAS variants. These results show the power of data sharing for the interpretation and analysis of rare variation, expand the spectrum of molecular causes of developmental disability to include RALA, and provide additional insight into the pathogenesis of human disease caused by mutations in small GTPases.

While many causes of developmental disabilities have been identified, a large number of affected children cannot be diagnosed despite extensive medical testing. Previously unknown genetic factors are likely to be the culprits in many of these cases. Using DNA sequencing, and by sharing information among many doctors and researchers, we have identified a set of individuals with developmental problems who all have changes to the same gene, RALA. The affected individuals all have similar symptoms, including intellectual disability, speech delay (or no speech), and problems with motor skills like walking. In nearly all of these cases (10 of 11), the genetic change found in the child was not inherited from either parent. The locations and biological properties of these changes suggest that they are likely to disrupt the normal functions of RALA. Functional experiments also show that the genetic changes found in these individuals alter two key functions of RALA. Together, we have provided evidence that genetic changes in RALA can cause developmental disabilities. These results will allow doctors and researchers to identify additional children with the same condition, providing a clinical diagnosis to these families and leading to new research opportunities.

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.

Further evidence for the involvement of EFL1 in a Shwachman-Diamond-like syndrome and expansion of the phenotypic features

Recent evidence has implicated EFL1 in a phenotype overlapping Shwachman-Diamond syndrome (SDS), with the functional interplay between EFL1 and the previously known causative gene SBDS accounting for the similarity in clinical features. Relatively little is known about the phenotypes associated with pathogenic variants in the EFL1 gene, but the initial indication was that phenotypes may be more severe, when compared with SDS. We report a pediatric patient who presented with a metaphyseal dysplasia and was found to have biallelic variants in EFL1 on reanalysis of trio whole-exome sequencing data. The variant had not been initially reported because of the research laboratory's focus on de novo variants. Subsequent phenotyping revealed variability in her manifestations. Although her metaphyseal abnormalities were more severe than in the original reported cohort with EFL1 variants, the bone marrow abnormalities were generally mild, and there was equivocal evidence for pancreatic insufficiency. Despite the limited number of reported patients, variants in EFL1 appear to cause a broader spectrum of symptoms that overlap with those seen in SDS. Our report adds to the evidence of EFL1 being associated with an SDS-like phenotype and provides information adding to our understanding of the phenotypic variability of this disorder. Our report also highlights the value of exome data reanalysis when a diagnosis is not initially apparent.

IRF2BPL Is Associated with Neurological Phenotypes

Interferon regulatory factor 2 binding protein-like (IRF2BPL) encodes a member of the IRF2BP family of transcriptional regulators. Currently the biological function of this gene is obscure, and the gene has not been associated with a Mendelian disease. Here we describe seven individuals who carry damaging heterozygous variants in IRF2BPL and are affected with neurological symptoms. Five individuals who carry IRF2BPL nonsense variants resulting in a premature stop codon display severe neurodevelopmental regression, hypotonia, progressive ataxia, seizures, and a lack of coordination. Two additional individuals, both with missense variants, display global developmental delay and seizures and a relatively milder phenotype than those with nonsense alleles. The IRF2BPL bioinformatics signature based on population genomics is consistent with a gene that is intolerant to variation. We show that the fruit-fly IRF2BPL ortholog, called pits (protein interacting with Ttk69 and Sin3A), is broadly detected, including in the nervous system. Complete loss of pits is lethal early in development, whereas partial knockdown with RNA interference in neurons leads to neurodegeneration, revealing a requirement for this gene in proper neuronal function and maintenance. The identified IRF2BPL nonsense variants behave as severe loss-of-function alleles in this model organism, and ectopic expression of the missense variants leads to a range of phenotypes. Taken together, our results show that IRF2BPL and pits are required in the nervous system in humans and flies, and their loss leads to a range of neurological phenotypes in both species.

Exome sequencing of an adolescent with nonalcoholic fatty liver disease identifies a clinically actionable case of Wilson disease

Diagnostic whole-exome sequencing has proven highly successful in a range of rare diseases, particularly early-onset genetic conditions. In more common conditions, however, exome sequencing for diagnostic purposes remains the exception. Here we describe a patient initially diagnosed with a common, complex liver disease, nonalcoholic fatty liver disease (NAFLD), who was determined to have Wilson disease (WD) upon research-related exome sequencing. The patient presented as a 14.5-yr-old adolescent with chronically elevated aminotransferases, normal ceruloplasmin, and histologic examination consistent with NAFLD with advanced fibrosis. He was enrolled in a large longitudinal study of patients with NAFLD and was found to have WD by exome sequencing performed 4 yr later. This new diagnosis, confirmed clinically by 24 h urine copper quantification, led to a change in the therapy from lifestyle counseling to directed treatment with d-penicillamine, a copper chelating agent. In this case, the likelihood of making the correct diagnosis and thereby choosing the appropriate treatment was increased by exome sequencing and careful interpretation. This example illustrates the utility of exome sequencing diagnostically in more common conditions not currently considered as targets for genome-wide evaluation and adds to a growing body of evidence that patients diagnosed with more common conditions often in fact have rarer genetically determined syndromes that have escaped clinical detection.

Functional variants in TBX2 are associated with a syndromic cardiovascular and skeletal developmental disorder

The 17 genes of the T-box family are transcriptional regulators that are involved in all stages of embryonic development, including craniofacial, brain, heart, skeleton and immune system. Malformation syndromes have been linked to many of the T-box genes. For example, haploinsufficiency of TBX1 is responsible for many structural malformations in DiGeorge syndrome caused by a chromosome 22q11.2 deletion. We report four individuals with an overlapping spectrum of craniofacial dysmorphisms, cardiac anomalies, skeletal malformations, immune deficiency, endocrine abnormalities and developmental impairments, reminiscent of DiGeorge syndrome, who are heterozygotes for TBX2 variants. The p.R20Q variant is shared by three affected family members in an autosomal dominant manner; the fourth unrelated individual has a de novo p.R305H mutation. Bioinformatics analyses indicate that these variants are rare and predict them to be damaging. In vitro transcriptional assays in cultured cells show that both variants result in reduced transcriptional repressor activity of TBX2. We also show that the variants result in reduced protein levels of TBX2. Heterologous over-expression studies in Drosophila demonstrate that both p.R20Q and p.R305H function as partial loss-of-function alleles. Hence, these and other data suggest that TBX2 is a novel candidate gene for a new multisystem malformation disorder.

Looking beyond the exome: a phenotype-first approach to molecular diagnostic resolution in rare and undiagnosed diseases

PurposeTo describe examples of missed pathogenic variants on whole-exome sequencing (WES) and the importance of deep phenotyping for further diagnostic testing.MethodsGuided by phenotypic information, three children with negative WES underwent targeted single-gene testing.ResultsIndividual 1 had a clinical diagnosis consistent with infantile systemic hyalinosis, although WES and a next-generation sequencing (NGS)-based ANTXR2 test were negative. Sanger sequencing of ANTXR2 revealed a homozygous single base pair insertion, previously missed by the WES variant caller software. Individual 2 had neurodevelopmental regression and cerebellar atrophy, with no diagnosis on WES. New clinical findings prompted Sanger sequencing and copy number testing of PLA2G6. A novel homozygous deletion of the noncoding exon 1 (not included in the WES capture kit) was detected, with extension into the promoter, confirming the clinical suspicion of infantile neuroaxonal dystrophy. Individual 3 had progressive ataxia, spasticity, and magnetic resonance image changes of vanishing white matter leukoencephalopathy. An NGS leukodystrophy gene panel and WES showed a heterozygous pathogenic variant in EIF2B5; no deletions/duplications were detected. Sanger sequencing of EIF2B5 showed a frameshift indel, probably missed owing to failure of alignment.ConclusionThese cases illustrate potential pitfalls of WES/NGS testing and the importance of phenotype-guided molecular testing in yielding diagnoses.

A Novel Mutation in Junctional Plakoglobin Causing Lethal Congenital Epidermolysis Bullosa

We report a case of neonatal generalized erythema and epidermolysis resulting from a novel mutation in the junctional plakoglobin gene causing truncation of the plakoglobin protein. Expedited genetic testing enabled diagnosis while the patient was in the neonatal intensive care unit, providing valuable information for the clinicians and family.

De Novo Mutations in PPP3CA Cause Severe Neurodevelopmental Disease with Seizures

Exome sequencing has readily enabled the discovery of the genetic mutations responsible for a wide range of diseases. This success has been particularly remarkable in the severe epilepsies and other neurodevelopmental diseases for which rare, often de novo, mutations play a significant role in disease risk. Despite significant progress, the high genetic heterogeneity of these disorders often requires large sample sizes to identify a critical mass of individuals with disease-causing mutations in a single gene. By pooling genetic findings across multiple studies, we have identified six individuals with severe developmental delay (6/6), refractory seizures (5/6), and similar dysmorphic features (3/6), each harboring a de novo mutation in PPP3CA. PPP3CA encodes the alpha isoform of a subunit of calcineurin. Calcineurin encodes a calcium- and calmodulin-dependent serine/threonine protein phosphatase that plays a role in a wide range of biological processes, including being a key regulator of synaptic vesicle recycling at nerve terminals. Five individuals with de novo PPP3CA mutations were identified among 4,760 trio probands with neurodevelopmental diseases; this is highly unlikely to occur by chance (p = 1.2 × 10^-8) given the size and mutability of the gene. Additionally, a sixth individual with a de novo mutation in PPP3CA was connected to this study through GeneMatcher. Based on these findings, we securely implicate PPP3CA in early-onset refractory epilepsy and further support the emerging role for synaptic dysregulation in epilepsy.

A Recurrent De Novo Variant in NACC1 Causes a Syndrome Characterized by Infantile Epilepsy, Cataracts, and Profound Developmental Delay

Whole-exome sequencing (WES) has increasingly enabled new pathogenic gene variant identification for undiagnosed neurodevelopmental disorders and provided insights into both gene function and disease biology. Here, we describe seven children with a neurodevelopmental disorder characterized by microcephaly, profound developmental delays and/or intellectual disability, cataracts, severe epilepsy including infantile spasms, irritability, failure to thrive, and stereotypic hand movements. Brain imaging in these individuals reveals delay in myelination and cerebral atrophy. We observe an identical recurrent de novo heterozygous c.892C>T (p.Arg298Trp) variant in the nucleus accumbens associated 1 (NACC1) gene in seven affected individuals. One of the seven individuals is mosaic for this variant. NACC1 encodes a transcriptional repressor implicated in gene expression and has not previously been associated with germline disorders. The probability of finding the same missense NACC1 variant by chance in 7 out of 17,228 individuals who underwent WES for diagnoses of neurodevelopmental phenotypes is extremely small and achieves genome-wide significance (p = 1.25 × 10-14). Selective constraint against missense variants in NACC1 makes this excess of an identical missense variant in all seven individuals more remarkable. Our findings are consistent with a germline recurrent mutational hotspot associated with an allele-specific neurodevelopmental phenotype in NACC1.

A Syndromic Neurodevelopmental Disorder Caused by De Novo Variants in EBF3

Early B cell factor 3 (EBF3) is a member of the highly evolutionarily conserved Collier/Olf/EBF (COE) family of transcription factors. Prior studies on invertebrate and vertebrate animals have shown that EBF3 homologs are essential for survival and that loss-of-function mutations are associated with a range of nervous system developmental defects, including perturbation of neuronal development and migration. Interestingly, aristaless-related homeobox (ARX), a homeobox-containing transcription factor critical for the regulation of nervous system development, transcriptionally represses EBF3 expression. However, human neurodevelopmental disorders related to EBF3 have not been reported. Here, we describe three individuals who are affected by global developmental delay, intellectual disability, and expressive speech disorder and carry de novo variants in EBF3. Associated features seen in these individuals include congenital hypotonia, structural CNS malformations, ataxia, and genitourinary abnormalities. The de novo variants affect a single conserved residue in a zinc finger motif crucial for DNA binding and are deleterious in a fly model. Our findings indicate that mutations in EBF3 cause a genetic neurodevelopmental syndrome and suggest that loss of EBF3 function might mediate a subset of neurologic phenotypes shared by ARX-related disorders, including intellectual disability, abnormal genitalia, and structural CNS malformations.

De Novo Truncating Variants in ASXL2 Are Associated with a Unique and Recognizable Clinical Phenotype

The ASXL genes (ASXL1, ASXL2, and ASXL3) participate in body patterning during embryogenesis and encode proteins involved in epigenetic regulation and assembly of transcription factors to specific genomic loci. Germline de novo truncating variants in ASXL1 and ASXL3 have been respectively implicated in causing Bohring-Opitz and Bainbridge-Ropers syndromes, which result in overlapping features of severe intellectual disability and dysmorphic features. ASXL2 has not yet been associated with a human Mendelian disorder. In this study, we performed whole-exome sequencing in six unrelated probands with developmental delay, macrocephaly, and dysmorphic features. All six had de novo truncating variants in ASXL2. A careful review enabled the recognition of a specific phenotype consisting of macrocephaly, prominent eyes, arched eyebrows, hypertelorism, a glabellar nevus flammeus, neonatal feeding difficulties, hypotonia, and developmental disabilities. Although overlapping features with Bohring-Opitz and Bainbridge-Ropers syndromes exist, features that distinguish the ASXL2-associated condition from ASXL1- and ASXL3-related disorders are macrocephaly, absence of growth retardation, and more variability in the degree of intellectual disabilities. We were also able to demonstrate with mRNA studies that these variants are likely to exert a dominant-negative effect, given that both alleles are expressed in blood and the mutated ASXL2 transcripts escape nonsense-mediated decay. In conclusion, de novo truncating variants in ASXL2 underlie a neurodevelopmental syndrome with a clinically recognizable phenotype. This report expands the germline disorders that are linked to the ASXL genes.

Mapping H4K20me3 onto the chromatin landscape of senescent cells indicates a function in control of cell senescence and tumor suppression through preservation of genetic and epigenetic stability

BACKGROUND

Histone modification H4K20me3 and its methyltransferase SUV420H2 have been implicated in suppression of tumorigenesis. The underlying mechanism is unclear, although H4K20me3 abundance increases during cellular senescence, a stable proliferation arrest and tumor suppressor process, triggered by diverse molecular cues, including activated oncogenes. Here, we investigate the function of H4K20me3 in senescence and tumor suppression.

RESULTS

Using immunofluorescence and ChIP-seq we determine the distribution of H4K20me3 in proliferating and senescent human cells. Altered H4K20me3 in senescence is coupled to H4K16ac and DNA methylation changes in senescence. In senescent cells, H4K20me3 is especially enriched at DNA sequences contained within specialized domains of senescence-associated heterochromatin foci (SAHF), as well as specific families of non-genic and genic repeats. Altered H4K20me3 does not correlate strongly with changes in gene expression between proliferating and senescent cells; however, in senescent cells, but not proliferating cells, H4K20me3 enrichment at gene bodies correlates inversely with gene expression, reflecting de novo accumulation of H4K20me3 at repressed genes in senescent cells, including at genes also repressed in proliferating cells. Although elevated SUV420H2 upregulates H4K20me3, this does not accelerate senescence of primary human cells. However, elevated SUV420H2/H4K20me3 reinforces oncogene-induced senescence-associated proliferation arrest and slows tumorigenesis in vivo.

CONCLUSIONS

These results corroborate a role for chromatin in underpinning the senescence phenotype but do not support a major role for H4K20me3 in initiation of senescence. Rather, we speculate that H4K20me3 plays a role in heterochromatinization and stabilization of the epigenome and genome of pre-malignant, oncogene-expressing senescent cells, thereby suppressing epigenetic and genetic instability and contributing to long-term senescence-mediated tumor suppression.

Subtelomeric p53 binding prevents accumulation of DNA damage at human telomeres

Telomeres and tumor suppressor protein TP53 (p53) function in genome protection, but a direct role of p53 at telomeres has not yet been described. Here, we have identified non-canonical p53-binding sites within the human subtelomeres that suppress the accumulation of DNA damage at telomeric repeat DNA. These non-canonical subtelomeric p53-binding sites conferred transcription enhancer-like functions that include an increase in local histone H3K9 and H3K27 acetylation and stimulation of subtelomeric transcripts, including telomere repeat-containing RNA (TERRA). p53 suppressed formation of telomere-associated γH2AX and prevented telomere DNA degradation in response to DNA damage stress. Our findings indicate that p53 provides a direct chromatin-associated protection to human telomeres, as well as other fragile genomic sites. We propose that p53-associated chromatin modifications enhance local DNA repair or protection to provide a previously unrecognized tumor suppressor function of p53.

Subtelomeric CTCF and cohesin binding site organization using improved subtelomere assemblies and a novel annotation pipeline

Mapping genome-wide data to human subtelomeres has been problematic due to the incomplete assembly and challenges of low-copy repetitive DNA elements. Here, we provide updated human subtelomere sequence assemblies that were extended by filling telomere-adjacent gaps using clone-based resources. A bioinformatic pipeline incorporating multiread mapping for annotation of the updated assemblies using short-read data sets was developed and implemented. Annotation of subtelomeric sequence features as well as mapping of CTCF and cohesin binding sites using ChIP-seq data sets from multiple human cell types confirmed that CTCF and cohesin bind within 3 kb of the start of terminal repeat tracts at many, but not all, subtelomeres. CTCF and cohesin co-occupancy were also enriched near internal telomere-like sequence (ITS) islands and the nonterminal boundaries of subtelomere repeat elements (SREs) in transformed lymphoblastoid cell lines (LCLs) and human embryonic stem cell (ES) lines, but were not significantly enriched in the primary fibroblast IMR90 cell line. Subtelomeric CTCF and cohesin sites predicted by ChIP-seq using our bioinformatics pipeline (but not predicted when only uniquely mapping reads were considered) were consistently validated by ChIP-qPCR. The colocalized CTCF and cohesin sites in SRE regions are candidates for mediating long-range chromatin interactions in the transcript-rich SRE region. A public browser for the integrated display of short-read sequence-based annotations relative to key subtelomere features such as the start of each terminal repeat tract, SRE identity and organization, and subtelomeric gene models was established.