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.

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.

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.

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.

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.

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.

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.

Novel variants in TUBA1A cause congenital fibrosis of the extraocular muscles with or without malformations of cortical brain development

Variants in multiple tubulin genes have been implicated in neurodevelopmental disorders, including malformations of cortical development (MCD) and congenital fibrosis of the extraocular muscles (CFEOM). Distinct missense variants in the beta-tubulin encoding genes TUBB3 and TUBB2B cause MCD, CFEOM, or both, suggesting substitution-specific mechanisms. Variants in the alpha tubulin-encoding gene TUBA1A have been associated with MCD, but not with CFEOM. Using exome sequencing (ES) and genome sequencing (GS), we identified 3 unrelated probands with CFEOM who harbored novel heterozygous TUBA1A missense variants c.1216C>G, p.(His406Asp); c.467G>A, p.(Arg156His); and c.1193T>G, p.(Met398Arg). MRI revealed small oculomotor-innervated muscles and asymmetrical caudate heads and lateral ventricles with or without corpus callosal thinning. Two of the three probands had MCD. Mutated amino acid residues localize either to the longitudinal interface at which α and β tubulins heterodimerize (Met398, His406) or to the lateral interface at which tubulin protofilaments interact (Arg156), and His406 interacts with the motor domain of kinesin-1. This series of individuals supports TUBA1A variants as a cause of CFEOM and expands our knowledge of tubulinopathies.

Phosphorylation state of pyruvate dehydrogenase and metabolite levels in turkey skeletal muscle in normal and pale, soft, exudative meats

1. Turkey production has increased dramatically as genetic selection has succeeded in increasing body weight and muscle yield to fulfil increasing consumer demand. However, producing fast-growing, heavily muscled birds is linked to increased heat stress susceptibility and can result in pale, soft, exudative (PSE) meat. Previous studies indicated that pyruvate dehydrogenase kinase 4 (PDK4) is significantly reduced in PSE samples, suggesting this as a candidate gene associated with the development of this problem.2. The objective of this study was to determine whether pre-market thermal challenge results in PSE meat as a result of differential expression of PDK4. Two genetic lines of turkeys were used in this study; the Randombred Control Line 2 (RBC2) and a commercial line. Turkeys were exposed to a pre-market thermal challenge of 12 h at 35°C followed by 12 h at 27°C for 5 d. Birds were slaughtered and processed according to industry standards. Pectoralis major samples were categorised as PSE or normal based on marinade uptake and cook loss indicators. In the first experiment, the relative expression of pyruvate dehydrogenase (PDH) and the phosphorylation state of PDH in normal and PSE turkey meat were analysed by western blotting. In the second experiment, the same samples were used to measure metabolite levels at 5 min post-mortem, comparing the normal to the PSE samples.3. The results of the first experiment showed that PSE samples had significantly lower total PDH (P = 0.029) compared to normal meat. However, there was no significant difference in the degree of phosphorylation of sites 1, 2 or 3. In the second experiment, there were no significant differences in glycogen, lactate, glycolytic potential or ATP when comparing PSE to control samples.4. These results suggested that a reduction in PDK4 expression alone does not explain the development of PSE meat.

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.

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.

Muscle from grass- and grain-fed cattle differs energetically

Insufficient acidification results in dark, firm, and dry beef. While this defect is often indicative of a stress event antemortem, muscle tissue may change in response to feeding regime. Longissimus dorsi muscle samples from 10 grain-fed and 10 grass-fed market weight, angus-crossbred beef cattle were collected postmortem. Lower (P < .05) L* and a* values were recorded for steaks from grass-fed cattle. Higher (P < .05) ultimate pH values were noted in lean of grass-fed cattle compared to grain-fed cattle, yet differences in lactate, glycogen and glucose were not detected. Further, increased (P < .05) ultimate pH values and lower (P < .05) lactate accumulations were noted when samples from grass-fed cattle were subjected to an in vitro glycolysis system. Muscle from grass-fed beef possessed nearly two-fold more (P < .05) succinate dehydrogenase and (P < .001) myoglobin than that of grain-fed cattle. These data show lean from grass-fed beef has greater enzymes reflective of oxidative metabolism and suggest dark lean from grass-fed cattle may be a function of more oxidative metabolism rather than a stress-related event antemortem.

Muscle from grass- and grain-fed cattle differs energetically

Insufficient acidification results in dark, firm, and dry beef. While this defect is often indicative of a stress event antemortem, muscle tissue may change in response to feeding regime. Longissimus dorsi muscle samples from 10 grain-fed and 10 grass-fed market weight, angus-crossbred beef cattle were collected postmortem. Lower (P < .05) L* and a* values were recorded for steaks from grass-fed cattle. Higher (P < .05) ultimate pH values were noted in lean of grass-fed cattle compared to grain-fed cattle, yet differences in lactate, glycogen and glucose were not detected. Further, increased (P < .05) ultimate pH values and lower (P < .05) lactate accumulations were noted when samples from grass-fed cattle were subjected to an in vitro glycolysis system. Muscle from grass-fed beef possessed nearly two-fold more (P < .05) succinate dehydrogenase and (P < .001) myoglobin than that of grain-fed cattle. These data show lean from grass-fed beef has greater enzymes reflective of oxidative metabolism and suggest dark lean from grass-fed cattle may be a function of more oxidative metabolism rather than a stress-related event antemortem.

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.

Chronic activation of AMP-activated protein kinase increases monocarboxylate transporter 2 and 4 expression in skeletal muscle

Acute activation of AMP-activated protein kinase (AMPK) increases monocarboxylate transporter (MCT) expression in skeletal muscle. However, the impact of chronic activation of AMPK on MCT expression in skeletal muscle is unknown. To investigate, MCT1, MCT2, and MCT4 mRNA expression and protein abundance were measured in the longissimus lumborum (glycolytic), masseter (oxidative), and heart from wild-type (control) and AMPK γ3 pigs. The AMPK γ3 gain in function mutation results in AMPK being constitutively active in glycolytic skeletal muscle and increases energy producing pathways. The MCT1 and MCT2 mRNA expression in muscle was lower ( < 0.05) from both wild-type and AMPK γ3 animals compared to other tissues. However, in both genotypes, MCT1 and MCT2 mRNA expression was greater ( < 0.05) in the masseter than the longissimus lumborum. The MCT1 protein was not detected in skeletal muscle, but MCT2 was greater ( < 0.05) in muscles with an oxidative muscle phenotype. Monocarboxylate transporter 2 was also detected in muscle mitochondria and may explain the differences between muscles. The MCT4 mRNA expression was intermediate among all tissues tested and greater ( < 0.05) in the longissimus lumborum than the masseter. Furthermore, MCT4 protein expression in the longissimus lumborum from AMPK γ3 animals was greater ( < 0.05) than in the longissimus lumborum from wild-type animals. In totality, these data indicate that chronic AMPK activation simultaneously increases MCT2 and MCT4 expression in skeletal muscle.

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.