The Mediator complex regulates protein-coding gene transcription by coordinating the interaction of upstream enhancers with the basal transcription machinery at the promoter. Pathogenic variants in Mediator subunits typically lead to neurodevelopmental or neurodegenerative disorders with variable clinical presentations, designated as MEDopathies. Here, we report the identification of 25 individuals from 18 families with bi-allelic MED16 variants who have a multiple congenital anomalies (MCAs)-intellectual disability syndrome. Intellectual disability, speech delay, and/or motor delay of variable severity were constant and associated with variable combinations of craniofacial defects (micro/retrognathia, cleft palate, and preauricular tags), anomalies of the extremities, and heart defects (predominantly tetralogy of Fallot). Visual impairment, deafness, and magnetic resonance imaging (MRI) abnormalities were also frequent. The 26 variants identified were comprised of eight predicted protein-truncating (three intragenic deletions, two frameshifts, and one nonsense and two essential splice site variants) and 18 missense or in-frame duplication variants affecting conserved residues, without clear correlation between phenotypic severity and variant type combination. Three-dimensional modeling indicated that the missense and duplication variants likely have a destabilizing effect on the structural elements of the protein. Immunofluorescence assays demonstrated protein mislocalization from the nucleus to the cytoplasm for 16 of the 17 variants studied. Homozygous mutant med16 zebrafish presented growth delay and increased mortality compared with wild-type fish, and Med16 knockout mice are preweaning lethal, highlighting the conserved requirement of MED16 for development. Overall, we describe an autosomal recessive MCAs-intellectual disability MEDopathy, emphasizing the importance of Mediator during neurodevelopment and suggesting that some tissues are particularly sensitive to the loss of certain subunits.

PURPOSE. To functionally evaluate novel human sequence-derived candidate genes and variants for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). METHODS. Through exome and genome sequencing of a genetically unsolved human oCCDD cohort, we previously reported the identification of variants in many candidate genes. Here, we describe a parallel study that prioritized a subset of these genes (43 human genes, 57 zebrafish genes) using a G0 CRISPR/Cas9-based knockout assay in zebrafish and generated F2 germline mutants for 17. We tested the functionality of variants of uncertain significance in known and novel candidate transcription factor-encoding genes through protein binding microarrays. RESULTS. We first demonstrated the feasibility of the G0 screen by targeting known oCCDD genes phox2a and mafba. Approximately 70% to 90% of gene-targeted G0 zebrafish embryos recapitulated germline homozygous null-equivalent phenotypes. Using this approach, we then identified three novel candidate oCCDD genes (SEMA3F, OLIG2, and FRMD4B) with putative contributions to human and zebrafish cranial motor development. In addition, protein binding microarrays demonstrated reduced or abolished DNA binding of human variants of uncertain significance in known and novel sequence-derived transcription factors PHOX2A (p.(Trp137Cys)), MAFB (p.(Glu223Lys)), and OLIG2 (p.(Arg156Leu)). CONCLUSIONS. This study nominates three strong novel candidate oCCDD genes (SEMA3F, OLIG2, and FRMD4B) and supports the functionality and putative pathogenicity of transcription factor candidate variants PHOX2A p.(Trp137Cys), MAFB p.(Glu223Lys), and OLIG2 p.(Arg156Leu). Our findings support that G0 loss-of-function screening in zebrafish can be coupled with human sequence analysis and protein binding microarrays to aid in prioritizing oCCDD candidate genes/variants.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3 , encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3 / TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish ( srpk3 −/− ; ttn.1 +/− ) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

CHASERR encodes a human long noncoding RNA (lncRNA) adjacent to CHD2, a coding gene in which de novo loss-of-function variants cause developmental and epileptic encephalopathy. Here, we report our findings in three unrelated children with a syndromic, early-onset neurodevelopmental disorder, each of whom had a de novo deletion in the CHASERR locus. The children had severe encephalopathy, shared facial dysmorphisms, cortical atrophy, and cerebral hypomyelination - a phenotype that is distinct from the phenotypes of patients with CHD2 haploinsufficiency. We found that the CHASERR deletion results in increased CHD2 protein abundance in patient-derived cell lines and increased expression of the CHD2 transcript in cis. These findings indicate that CHD2 has bidirectional dosage sensitivity in human disease, and we recommend that other lncRNA-encoding genes be evaluated, particularly those upstream of genes associated with mendelian disorders. (Funded by the National Human Genome Research Institute and others.) Deletions in a gene (CHASERR) that encodes a long noncoding RNA cause severe developmental delay and an increase in expression of a neighboring gene, CHD2.

Background Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined.Methods We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center.Results We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments.Conclusions The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.) Genetic diagnosis of rare diseases is made through a variety of methods. This study gauged the diagnostic yield of genome sequencing after negative results from exome sequencing and other methods.

Genes encoding long non-coding RNAs (lncRNAs) comprise a large fraction of the human genome, yet haploinsufficiency of a lncRNA has not been shown to cause a Mendelian disease. CHASERR is a highly conserved human lncRNA adjacent to CHD2-a coding gene in which de novo loss-of-function variants cause developmental and epileptic encephalopathy. Here we report three unrelated individuals each harboring an ultra-rare heterozygous de novo deletion in the CHASERR locus. We report similarities in severe developmental delay, facial dysmorphisms, and cerebral dysmyelination in these individuals, distinguishing them from the phenotypic spectrum of CHD2 haploinsufficiency. We demonstrate reduced CHASERR mRNA expression and corresponding increased CHD2 mRNA and protein in whole blood and patient-derived cell lines-specifically increased expression of the CHD2 allele in cis with the CHASERR deletion, as predicted from a prior mouse model of Chaserr haploinsufficiency. We show for the first time that de novo structural variants facilitated by Alu-mediated non-allelic homologous recombination led to deletion of a non-coding element (the lncRNA CHASERR) to cause a rare syndromic neurodevelopmental disorder. We also demonstrate that CHD2 has bidirectional dosage sensitivity in human disease. This work highlights the need to carefully evaluate other lncRNAs, particularly those upstream of genes associated with Mendelian disorders.

The ADAT2/ADAT3 complex catalyzes the adenosine to inosine modification at the wobble position of eukaryotic tRNAs. Mutations in ADAT3 , the catalytically inactive subunit of the ADAT2/ADAT3 complex, have been identified in patients presenting with severe neurodevelopmental disorders (NDDs). Yet, the physiological function of ADAT2/ADAT3 complex during brain development remains totally unknown. Here we showed that maintaining a proper level of ADAT2/ADAT3 catalytic activity is required for correct radial migration of projection neurons in the developing mouse cortex. In addition, we not only reported 20 new NDD patients carrying biallelic variants in ADAT3 but also deeply characterized the impact of those variants on ADAT2/ADAT3 structure, biochemical properties, enzymatic activity and tRNAs editing and abundance. We demonstrated that all the identified variants alter both the abundance and the activity of the complex leading to a significant decrease of I 34 with direct consequence on their steady-state. Using in vivo complementation assays, we correlated the severity of the migration phenotype with the degree of the loss of function caused by the variants. Altogether, our results indicate a critical role of ADAT2/ADAT3 during cortical development and provide cellular and molecular insights into the pathogenicity of ADAT3-related neurodevelopmental disorder.

Purpose: Fem1 homolog B (FEM1B) acts as a substrate recognition subunit for ubiquitin ligase complexes belonging to the CULLIN 2-based E3 family. Several biological functions have been proposed for FEM1B, including a structurally resolved function as a sensor for redox cell status by controlling mitochondrial activity, but its implication in human disease remains elusive. Methods: To understand the involvement of FEM1B in human disease, we made use of Matchmaker exchange platforms to identify individuals with de novo variants in FEM1B and performed their clinical evaluation. We performed functional validation using primary neuronal cultures and in utero electroporation assays, as well as experiments on patient's cells. Results: Five individuals with a recurrent de novo missense variant in FEM1B were identified: NM_015322.5:c.377G>A NP_056137.1:p.(Arg126Gln) (FEM1BR126Q). Affected individuals shared a severe neurodevelopmental disorder with behavioral phenotypes and a variable set of malformations, including brain anomalies, clubfeet, skeletal abnormalities, and facial dysmorphism. Overexpression of the FEM1BR126Q variant but not FEM1B wild-type protein, during mouse brain development, resulted in delayed neuronal migration of the target cells. In addition, the individuals' cells exhibited signs of oxidative stress and induction of type I interferon signaling. Conclusion: Overall, our data indicate that p.(Arg126Gln) induces aberrant FEM1B activation, resulting in a gain-of-function mechanism associated with a severe syndromic developmental disorder in humans.

Abstract Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 – as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Purpose To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. Results Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. Conclusion This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Background Causal variants underlying rare disorders may remain elusive even after expansive gene panels or exome sequencing (ES). Clinicians and researchers may then turn to genome sequencing (GS), though the added value of this technique and its optimal use remain poorly defined. We therefore investigated the advantages of GS within a phenotypically diverse cohort. Methods GS was performed for 744 individuals with rare disease who were genetically undiagnosed. Analysis included review of single nucleotide, indel, structural, and mitochondrial variants. Results We successfully solved 218/744 (29.3%) cases using GS, with most solves involving established disease genes (157/218, 72.0%). Of all solved cases, 148 (67.9%) had previously had non-diagnostic ES. We systematically evaluated the 218 causal variants for features requiring GS to identify and 61/218 (28.0%) met these criteria, representing 8.2% of the entire cohort. These included small structural variants (13), copy neutral inversions and complex rearrangements (8), tandem repeat expansions (6), deep intronic variants (15), and coding variants that may be more easily found using GS related to uniformity of coverage (19). Conclusion We describe the diagnostic yield of GS in a large and diverse cohort, illustrating several types of pathogenic variation eluding ES or other techniques. Our results reveal a higher diagnostic yield of GS, supporting the utility of a genome-first approach, with consideration of GS as a secondary or tertiary test when higher-resolution structural variant analysis is needed or there is a strong clinical suspicion for a condition and prior targeted genetic testing has been negative.

PurposeFINCA disease (Fibrosis, Neurodegeneration and Cerebral Angiomatosis, OMIM 618278) is an infantile-onset neurodevelopmental and multiorgan disease. Since our initial report in 2018, additional patients have been described. FINCA is the first human disease caused by recessive variants in the highly conserved NHLRC2 gene. Our previous studies have shown that Nhlrc2-null mouse embryos die during gastrulation, indicating the essential role of the protein in embryonic development. Defect in NHLRC2 leads to cerebral neurodegeneration and severe pulmonary, hepatic and cardiac fibrosis. Despite having a structure suggestive of an enzymatic role and the clinical importance of NHLRC2 in multiple organs, the specific physiological role of the protein is unknown. MethodsThe clinical histories of five novel FINCA patients diagnosed with whole exome sequencing were reviewed. Segregation analysis of the biallelic, potentially pathogenic NHLRC2 variants was performed using Sanger sequencing. Studies on neuropathology and NHLRC2 expression in different brain regions were performed on autopsy samples of three previously described deceased FINCA patients. ResultsOne patient was homozygous for the pathogenic variant c.442G > T, while the other four were compound heterozygous for this variant and two other pathogenic NHLRC2 gene variants. All five patients presented with multiorgan dysfunction with neurodevelopmental delay, recurrent infections and macrocytic anemia as key features. Interstitial lung disease was pronounced in infancy but often stabilized. Autopsy samples revealed widespread, albeit at a lower intensity than the control, NHLRC2 expression in the brain. ConclusionThis report expands on the characteristic clinical features of FINCA disease. Presentation is typically in infancy, and although patients can live to late adulthood, the key clinical and histopathological features are fibrosis, infection susceptibility/immunodeficiency/intellectual disability, neurodevelopmental disorder/neurodegeneration and chronic anemia/cerebral angiomatosis (hence the acronym FINCA) that enable an early diagnosis confirmed by genetic investigations.

BackgroundO’Donnell-Luria-Rodan syndrome (ODLURO) is an autosomal-dominant neurodevelopmental disorder caused by pathogenic, mostly truncating variants in KMT2E. It was first described by O’Donnell-Luria et al in 2019 in a cohort of 38 patients. Clinical features encompass macrocephaly, mild intellectual disability (ID), autism spectrum disorder (ASD) susceptibility and seizure susceptibility.MethodsAffected individuals were ascertained at paediatric and genetic centres in various countries by diagnostic chromosome microarray or exome/genome sequencing. Patients were collected into a case cohort and were systematically phenotyped where possible.ResultsWe report 18 additional patients from 17 families with genetically confirmed ODLURO. We identified 15 different heterozygous likely pathogenic or pathogenic sequence variants (14 novel) and two partial microdeletions of KMT2E. We confirm and refine the phenotypic spectrum of the KMT2E-related neurodevelopmental disorder, especially concerning cognitive development, with rather mild ID and macrocephaly with subtle facial features in most patients. We observe a high prevalence of ASD in our cohort (41%), while seizures are present in only two patients. We extend the phenotypic spectrum by sleep disturbances.ConclusionOur study, bringing the total of known patients with ODLURO to more than 60 within 2 years of the first publication, suggests an unexpectedly high relative frequency of this syndrome worldwide. It seems likely that ODLURO, although just recently described, is among the more common single-gene aetiologies of neurodevelopmental delay and ASD. We present the second systematic case series of patients with ODLURO, further refining the mutational and phenotypic spectrum of this not-so-rare syndrome.