Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Missense variants in DPYSL5 cause a neurodevelopmental disorder with corpus callosum agenesis and cerebellar abnormalities.

The collapsin response mediator protein (CRMP) family proteins are intracellular mediators of neurotrophic factors regulating neurite structure/spine formation and are essential for dendrite patterning and directional axonal pathfinding during brain developmental processes. Among this family, CRMP5/DPYSL5 plays a significant role in neuronal migration, axonal guidance, dendrite outgrowth, and synapse formation by interacting with microtubules. Here, we report the identification of missense mutations in DPYSL5 in nine individuals with brain malformations, including corpus callosum agenesis and/or posterior fossa abnormalities, associated with variable degrees of intellectual disability. A recurrent de novo p.Glu41Lys variant was found in eight unrelated patients, and a p.Gly47Arg variant was identified in one individual from the first family reported with Ritscher-Schinzel syndrome. Functional analyses of the two missense mutations revealed impaired dendritic outgrowth processes in young developing hippocampal primary neuronal cultures. We further demonstrated that these mutations, both located in the same loop on the surface of DPYSL5 monomers and oligomers, reduced the interaction of DPYSL5 with neuronal cytoskeleton-associated proteins MAP2 and ßIII-tubulin. Our findings collectively indicate that the p.Glu41Lys and p.Gly47Arg variants impair DPYSL5 function on dendritic outgrowth regulation by preventing the formation of the ternary complex with MAP2 and ßIII-tubulin, ultimately leading to abnormal brain development. This study adds DPYSL5 to the list of genes implicated in brain malformation and in neurodevelopmental disorders.

A form of muscular dystrophy associated with pathogenic variants in JAG2.

JAG2 encodes the Notch ligand Jagged2. The conserved Notch signaling pathway contributes to the development and homeostasis of multiple tissues, including skeletal muscle. We studied an international cohort of 23 individuals with genetically unsolved muscular dystrophy from 13 unrelated families. Whole-exome sequencing identified rare homozygous or compound heterozygous JAG2 variants in all 13 families. The identified bi-allelic variants include 10 missense variants that disrupt highly conserved amino acids, a nonsense variant, two frameshift variants, an in-frame deletion, and a microdeletion encompassing JAG2. Onset of muscle weakness occurred from infancy to young adulthood. Serum creatine kinase (CK) levels were normal or mildly elevated. Muscle histology was primarily dystrophic. MRI of the lower extremities revealed a distinct, slightly asymmetric pattern of muscle involvement with cores of preserved and affected muscles in quadriceps and tibialis anterior, in some cases resembling patterns seen in POGLUT1-associated muscular dystrophy. Transcriptome analysis of muscle tissue from two participants suggested misregulation of genes involved in myogenesis, including PAX7. In complementary studies, Jag2 downregulation in murine myoblasts led to downregulation of multiple components of the Notch pathway, including Megf10. Investigations in Drosophila suggested an interaction between Serrate and Drpr, the fly orthologs of JAG1/JAG2 and MEGF10, respectively. In silico analysis predicted that many Jagged2 missense variants are associated with structural changes and protein misfolding. In summary, we describe a muscular dystrophy associated with pathogenic variants in JAG2 and evidence suggests a disease mechanism related to Notch pathway dysfunction.

A recurrent missense variant in the E3 ubiquitin ligase substrate recognition subunit FEM1B causes a rare syndromic 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.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation 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.

Recessive variants in COL25A1 gene as novel cause of arthrogryposis multiplex congenita with ocular congenital cranial dysinnervation disorder.

A proper interaction between muscle-derived collagen XXV and its motor neuron-derived receptors protein tyrosine phosphatases σ and δ (PTP σ/δ) is indispensable for intramuscular motor innervation. Despite this, thus far, pathogenic recessive variants in the COL25A1 gene had only been detected in a few patients with isolated ocular congenital cranial dysinnervation disorders. Here we describe five patients from three unrelated families with recessive missense and splice site COL25A1 variants presenting with a recognizable phenotype characterized by arthrogryposis multiplex congenita with or without an ocular congenital cranial dysinnervation disorder phenotype. The clinical features of the older patients remained stable over time, without central nervous system involvement. This study extends the phenotypic and genotypic spectrum of COL25A1 related conditions, and further adds to our knowledge of the complex process of intramuscular motor innervation. Our observations indicate a role for collagen XXV in regulating the appropriate innervation not only of extraocular muscles, but also of bulbar, axial, and limb muscles in the human.

Major brain malformations: corpus callosum dysgenesis, agenesis of septum pellucidum and polymicrogyria in patients with BCORL1-related disorders.

OBJECTIVE: BCORL1, a transcriptional co-repressor, has a role in cortical migration, neuronal differentiation, maturation, and cerebellar development. We describe BCORL1 as a new genetic cause for major brain malformations. METHODS AND RESULTS: We report three patients from two unrelated families with neonatal onset intractable epilepsy and profound global developmental delay. Brain MRI of two siblings from the first family depicted hypoplastic corpus callosum and septal agenesis (ASP) in the older brother and unilateral perisylvian polymicrogyria (PMG) in the younger one. MRI of the patient from the second family demonstrated complete agenesis of corpus callosum (CC). Whole Exome Sequencing revealed a novel hemizygous variant in NM_021946.5 (BCORL1):c.796C>T (p.Pro266Ser) in the two siblings from the first family and the NM_021946.5 (BCORL1): c.3376G>A; p.Asp1126Asn variant in the patient from the second family, both variants inherited from healthy mothers. We reviewed the patients' charts and MRIs and compared the phenotype to the other published BCORL1-related cases. Brain malformations have not been previously described in association with the BCORL1 phenotype. We discuss the potential influence of BCORL1 on brain development. CONCLUSIONS: We suggest that BCORL1 variants present with a spectrum of neurodevelopmental disorders and can lead to major brain malformations originating at different stages of fetal development. We suggest adding BCORL1 to the genetic causes of PMG, ASP, and CC dysgenesis.

Cochlear nerve deficiency in SOX11-related Coffin-Siris syndrome.

The phenotypic spectrum of SOX11-related Coffin-Siris syndrome (CSS) is expanding with reports of new associations. SOX11 is implicated in neurogenesis and inner ear development. Cochlear nerve deficiency, absence or hypoplasia, is commonly associated with cochlear canal stenosis or with CHARGE syndrome, a monogenic condition that affects inner ear development. SOX11 is a transcription factor essential for neuronal identity, highly correlated with the expression of CHD7, which regulates SOX11. We present two unrelated probands, each with novel de novo SOX11 likely pathogenic variants and phenotypic manifestations of CSS including global developmental delay, growth deficiency, and hypoplastic nails. They have unilateral sensorineural hearing loss due to cochlear nerve deficiency confirmed on MRI. SOX11 is implicated in sensory neuron survival and maturation. It is highly expressed in the developing inner ear. Homozygous ablation of SOX11 in a mouse model resulted in a reduction in sensory neuron survival and decreased axonal growth. A heterozygous knockout mice model had hearing impairment with grossly normal inner ear structures like the two probands reported. We propose cochlear nerve deficiency as a new phenotypic feature of SOX11-related CSS. Magnetic resonance imaging is useful in delineating the cochlear nerve deficiency and other CSS-related brain malformations.

Alternative genomic diagnoses for individuals with a clinical diagnosis of Dubowitz syndrome.

Dubowitz syndrome (DubS) is considered a recognizable syndrome characterized by a distinctive facial appearance and deficits in growth and development. There have been over 200 individuals reported with Dubowitz or a "Dubowitz-like" condition, although no single gene has been implicated as responsible for its cause. We have performed exome (ES) or genome sequencing (GS) for 31 individuals clinically diagnosed with DubS. After genome-wide sequencing, rare variant filtering and computational and Mendelian genomic analyses, a presumptive molecular diagnosis was made in 13/27 (48%) families. The molecular diagnoses included biallelic variants in SKIV2L, SLC35C1, BRCA1, NSUN2; de novo variants in ARID1B, ARID1A, CREBBP, POGZ, TAF1, HDAC8, and copy-number variation at1p36.11(ARID1A), 8q22.2(VPS13B), Xp22, and Xq13(HDAC8). Variants of unknown significance in known disease genes, and also in genes of uncertain significance, were observed in 7/27 (26%) additional families. Only one gene, HDAC8, could explain the phenotype in more than one family (N = 2). All but two of the genomic diagnoses were for genes discovered, or for conditions recognized, since the introduction of next-generation sequencing. Overall, the DubS-like clinical phenotype is associated with extensive locus heterogeneity and the molecular diagnoses made are for emerging clinical conditions sharing characteristic features that overlap the DubS phenotype.

Sequential targeted exome sequencing of 1001 patients affected by unexplained limb-girdle weakness.

PURPOSE: Several hundred genetic muscle diseases have been described, all of which are rare. Their clinical and genetic heterogeneity means that a genetic diagnosis is challenging. We established an international consortium, MYO-SEQ, to aid the work-ups of muscle disease patients and to better understand disease etiology. METHODS: Exome sequencing was applied to 1001 undiagnosed patients recruited from more than 40 neuromuscular disease referral centers; standardized phenotypic information was collected for each patient. Exomes were examined for variants in 429 genes associated with muscle conditions. RESULTS: We identified suspected pathogenic variants in 52% of patients across 87 genes. We detected 401 novel variants, 116 of which were recurrent. Variants in CAPN3, DYSF, ANO5, DMD, RYR1, TTN, COL6A2, and SGCA collectively accounted for over half of the solved cases; while variants in newer disease genes, such as BVES and POGLUT1, were also found. The remaining well-characterized unsolved patients (48%) need further investigation. CONCLUSION: Using our unique infrastructure, we developed a pathway to expedite muscle disease diagnoses. Our data suggest that exome sequencing should be used for pathogenic variant detection in patients with suspected genetic muscle diseases, focusing first on the most common disease genes described here, and subsequently in rarer and newly characterized disease genes.

Phenotype delineation of ZNF462 related syndrome.

Zinc finger protein 462 (ZNF462) is a relatively newly discovered vertebrate specific protein with known critical roles in embryonic development in animal models. Two case reports and a case series study have described the phenotype of 10 individuals with ZNF462 loss of function variants. Herein, we present 14 new individuals with loss of function variants to the previous studies to delineate the syndrome of loss of function in ZNF462. Collectively, these 24 individuals present with recurring phenotypes that define a multiple congenital anomaly syndrome. Most have some form of developmental delay (79%) and a minority has autism spectrum disorder (33%). Characteristic facial features include ptosis (83%), down slanting palpebral fissures (58%), exaggerated Cupid's bow/wide philtrum (54%), and arched eyebrows (50%). Metopic ridging or craniosynostosis was found in a third of study participants and feeding problems in half. Other phenotype characteristics include dysgenesis of the corpus callosum in 25% of individuals, hypotonia in half, and structural heart defects in 21%. Using facial analysis technology, a computer algorithm applying deep learning was able to accurately differentiate individuals with ZNF462 loss of function variants from individuals with Noonan syndrome and healthy controls. In summary, we describe a multiple congenital anomaly syndrome associated with haploinsufficiency of ZNF462 that has distinct clinical characteristics and facial features.

Neurodevelopmental disorders associated variants in ADAT3 disrupt the activity of the ADAT2/ADAT3 tRNA deaminase complex and impair neuronal migration.

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 7 new NDD patients carrying biallelic variants in ADAT3 but also deeply characterize 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 expression 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.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders.

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.

Gene identification for ocular congenital cranial motor neuron disorders using human sequencing, zebrafish screening, and protein binding microarrays.

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 identified variants in 80 strong candidate genes. Here, we further prioritized a subset of these (43 human genes, 57 zebrafish genes) using a G0 CRISPR/Cas9-based knockout assay in zebrafish and generated F2 germline mutants for seventeen. 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. 70-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.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

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.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation 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 pedigrees 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.