Federated analysis of autosomal recessive coding variants in 29,745 developmental disorder patients from diverse populations.

Autosomal recessive coding variants are well-known causes of rare disorders. We quantified the contribution of these variants to developmental disorders in a large, ancestrally diverse cohort comprising 29,745 trios, of whom 20.4% had genetically inferred non-European ancestries. The estimated fraction of patients attributable to exome-wide autosomal recessive coding variants ranged from ~2-19% across genetically inferred ancestry groups and was significantly correlated with average autozygosity. Established autosomal recessive developmental disorder-associated (ARDD) genes explained 84.0% of the total autosomal recessive coding burden, and 34.4% of the burden in these established genes was explained by variants not already reported as pathogenic in ClinVar. Statistical analyses identified two novel ARDD genes: KBTBD2 and ZDHHC16. This study expands our understanding of the genetic architecture of developmental disorders across diverse genetically inferred ancestry groups and suggests that improving strategies for interpreting missense variants in known ARDD genes may help diagnose more patients than discovering the remaining genes.

Characterization of a missense variant in COG5 in a Tunisian patient with COG5-CDG syndrome and insights into the effect of non-synonymous variants on COG5 protein.

The clinical diagnosis of patients with multisystem involvement including a pronounced neurologic damage is challenging. High-throughput sequencing methods remains crucial to provide an accurate diagnosis. In this study, we reported a Tunisian patient manifesting hypotonia and global developmental delay with visual and skin abnormalities. Exome sequencing was conducted followed by segregation analysis and, subsequently additional investigations. In silico analysis of non-synonymous variants (nsSNPs) described in COG5 in conserved positions was made. Results revealed a homozygous missense variant c.298 C > T (p.Leu100Phe) in the COG5 inherited from both parents. This variant altered both protein solubility and stability, in addition to a putative disruption of the COG5-COG7 interaction. This disruption has been confirmed using patient-derived cells in vitro in a COG5 co-immuno-precipitation, where interaction with binding partner COG7 was abrogated. Hence, we established the COG5-CDG diagnosis. Clinically, the patient shared common features with the already described cases with the report of the ichtyosis as a new manifestation. Conversely, the CADD scoring revealed 19 putatively pathogenic nsSNPs (Minor Allele Frequency MAF < 0.001, CADD > 30), 11 of which had a significant impact on the solubility and/or stability of COG5. These properties seem to be disrupted by six of the seven missense COG5-CDG variants. In conclusion, our study expands the genetic and phenotypic spectrum of COG5-CDG disease and highlight the utility of the next generation sequencing as a powerful tool in accurate diagnosis. Our results shed light on a likely molecular mechanism underlying the pathogenic effect of missense COG5 variants, which is the alteration of COG5 stability and solubility.

DISP1 deficiency: Monoallelic and biallelic variants cause a spectrum of midline craniofacial malformations.

PURPOSE: DISP1 encodes a transmembrane protein that regulates the secretion of the morphogen, Sonic hedgehog, a deficiency of which is a major cause of holoprosencephaly (HPE). This disorder covers a spectrum of brain and midline craniofacial malformations. The objective of the present study was to better delineate the clinical phenotypes associated with division transporter dispatched-1 (DISP1) variants. METHODS: This study was based on the identification of at least 1 pathogenic variant of the DISP1 gene in individuals for whom detailed clinical data were available. RESULTS: A total of 23 DISP1 variants were identified in heterozygous, compound heterozygous or homozygous states in 25 individuals with midline craniofacial defects. Most cases were minor forms of HPE, with craniofacial features such as orofacial cleft, solitary median maxillary central incisor, and congenital nasal pyriform aperture stenosis. These individuals had either monoallelic loss-of-function variants or biallelic missense variants in DISP1. In individuals with severe HPE, the DISP1 variants were commonly found associated with a variant in another HPE-linked gene (ie, oligogenic inheritance). CONCLUSION: The genetic findings we have acquired demonstrate a significant involvement of DISP1 variants in the phenotypic spectrum of midline defects. This underlines its importance as a crucial element in the efficient secretion of Sonic hedgehog. We also demonstrated that the very rare solitary median maxillary central incisor and congenital nasal pyriform aperture stenosis combination is part of the DISP1-related phenotype. The present study highlights the clinical risks to be flagged up during genetic counseling after the discovery of a pathogenic DISP1 variant.

Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish.

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

Deep transcriptome profiling reveals limited conservation of A-to-I RNA editing in Xenopus.

BACKGROUND: Xenopus has served as a valuable model system for biomedical research over the past decades. Notably, ADAR was first detected in frog oocytes and embryos as an activity that unwinds RNA duplexes. However, the scope of A-to-I RNA editing by the ADAR enzymes in Xenopus remains underexplored. RESULTS: Here, we identify millions of editing events in Xenopus with high accuracy and systematically map the editome across developmental stages, adult organs, and species. We report diverse spatiotemporal patterns of editing with deamination activity highest in early embryogenesis before zygotic genome activation and in the ovary. Strikingly, editing events are poorly conserved across different Xenopus species. Even sites that are detected in both X. laevis and X. tropicalis show largely divergent editing levels or developmental profiles. In protein-coding regions, only a small subset of sites that are found mostly in the brain are well conserved between frogs and mammals. CONCLUSIONS: Collectively, our work provides fresh insights into ADAR activity in vertebrates and suggest that species-specific editing may play a role in each animal's unique physiology or environmental adaptation.

SMCHD1 has separable roles in chromatin architecture and gene silencing that could be targeted in disease.

The interplay between 3D chromatin architecture and gene silencing is incompletely understood. Here, we report a novel point mutation in the non-canonical SMC protein SMCHD1 that enhances its silencing capacity at endogenous developmental targets. Moreover, it also results in enhanced silencing at the facioscapulohumeral muscular dystrophy associated macrosatellite-array, D4Z4, resulting in enhanced repression of DUX4 encoded by this repeat. Heightened SMCHD1 silencing perturbs developmental Hox gene activation, causing a homeotic transformation in mice. Paradoxically, the mutant SMCHD1 appears to enhance insulation against other epigenetic regulators, including PRC2 and CTCF, while depleting long range chromatin interactions akin to what is observed in the absence of SMCHD1. These data suggest that SMCHD1's role in long range chromatin interactions is not directly linked to gene silencing or insulating the chromatin, refining the model for how the different levels of SMCHD1-mediated chromatin regulation interact to bring about gene silencing in normal development and disease.

The Opto-inflammasome in zebrafish as a tool to study cell and tissue responses to speck formation and cell death.

The inflammasome is a conserved structure for the intracellular detection of danger or pathogen signals. As a large intracellular multiprotein signaling platform, it activates downstream effectors that initiate a rapid necrotic programmed cell death (PCD) termed pyroptosis and activation and secretion of pro-inflammatory cytokines to warn and activate surrounding cells. However, inflammasome activation is difficult to control experimentally on a single-cell level using canonical triggers. We constructed Opto-ASC, a light-responsive form of the inflammasome adaptor protein ASC (Apoptosis-Associated Speck-Like Protein Containing a CARD) which allows tight control of inflammasome formation in vivo. We introduced a cassette of this construct under the control of a heat shock element into zebrafish in which we can now induce ASC inflammasome (speck) formation in individual cells of the skin. We find that cell death resulting from ASC speck formation is morphologically distinct from apoptosis in periderm cells but not in basal cells. ASC-induced PCD can lead to apical or basal extrusion from the periderm. The apical extrusion in periderm cells depends on Caspb and triggers a strong Ca2+ signaling response in nearby cells.

In skeletal muscle and neural crest cells, SMCHD1 regulates biological pathways relevant for Bosma syndrome and facioscapulohumeral dystrophy phenotype.

Many genetic syndromes are linked to mutations in genes encoding factors that guide chromatin organization. Among them, several distinct rare genetic diseases are linked to mutations in SMCHD1 that encodes the structural maintenance of chromosomes flexible hinge domain containing 1 chromatin-associated factor. In humans, its function as well as the impact of its mutations remains poorly defined. To fill this gap, we determined the episignature associated with heterozygous SMCHD1 variants in primary cells and cell lineages derived from induced pluripotent stem cells for Bosma arhinia and microphthalmia syndrome (BAMS) and type 2 facioscapulohumeral dystrophy (FSHD2). In human tissues, SMCHD1 regulates the distribution of methylated CpGs, H3K27 trimethylation and CTCF at repressed chromatin but also at euchromatin. Based on the exploration of tissues affected either in FSHD or in BAMS, i.e. skeletal muscle fibers and neural crest stem cells, respectively, our results emphasize multiple functions for SMCHD1, in chromatin compaction, chromatin insulation and gene regulation with variable targets or phenotypical outcomes. We concluded that in rare genetic diseases, SMCHD1 variants impact gene expression in two ways: (i) by changing the chromatin context at a number of euchromatin loci or (ii) by directly regulating some loci encoding master transcription factors required for cell fate determination and tissue differentiation.

Inactivation of DRG1, encoding a translation factor GTPase, causes a recessive neurodevelopmental disorder.

PURPOSE: Developmentally regulated Guanosine-5'-triphosphate-binding protein 1 (DRG1) is a highly conserved member of a class of GTPases implicated in translation. Although the expression of mammalian DRG1 is elevated in the central nervous system during development, and its function has been implicated in fundamental cellular processes, no pathogenic germline variants have yet been identified. Here, we characterize the clinical and biochemical consequences of DRG1 variants. METHODS: We collate clinical information of 4 individuals with germline DRG1 variants and use in silico, in vitro, and cell-based studies to study the pathogenicity of these alleles. RESULTS: We identified private germline DRG1 variants, including 3 stop-gained p.Gly54∗, p.Arg140∗, p.Lys263∗, and a p.Asn248Phe missense variant. These alleles are recessively inherited in 4 affected individuals from 3 distinct families and cause a neurodevelopmental disorder with global developmental delay, primary microcephaly, short stature, and craniofacial anomalies. We show that these loss-of-function variants (1) severely disrupt DRG1 messenger RNA/protein stability in patient-derived fibroblasts, (2) impair its GTPase activity, and (3) compromise its binding to partner protein ZC3H15. Consistent with the importance of DRG1 in humans, targeted inactivation of mouse Drg1 resulted in preweaning lethality. CONCLUSION: Our work defines a new Mendelian disorder of DRG1 deficiency. This study highlights DRG1's importance for normal mammalian development and underscores the significance of translation factor GTPases in human physiology and homeostasis.

RAF1 deficiency causes a lethal syndrome that underscores RTK signaling during embryogenesis.

Somatic and germline gain-of-function point mutations in RAF, one of the first oncogenes to be discovered in humans, delineate a group of tumor-prone syndromes known as the RASopathies. In this study, we document the first human phenotype resulting from the germline loss-of-function of the proto-oncogene RAF1 (a.k.a. CRAF). In a consanguineous family, we uncovered a homozygous p.Thr543Met variant segregating with a neonatal lethal syndrome with cutaneous, craniofacial, cardiac, and limb anomalies. Structure-based prediction and functional tests using human knock-in cells showed that threonine 543 is essential to: (i) ensure RAF1's stability and phosphorylation, (ii) maintain its kinase activity toward substrates of the MAPK pathway, and (iii) protect from stress-induced apoptosis mediated by ASK1. In Xenopus embryos, mutant RAF1T543M failed to phenocopy the effects of normal and overactive FGF/MAPK signaling, confirming its hypomorphic activity. Collectively, our data disclose the genetic and molecular etiology of a novel lethal syndrome with progeroid features, highlighting the importance of RTK signaling for human development and homeostasis.

A homozygous p.Leu813Pro gain-of-function NLRP1 variant causes phenotypes of different severity in two siblings.

BACKGROUND: A trio exome sequencing study identified a previously unreported NLRP1 gene variant resulting in a p.Leu813Pro substitution of the LRR (leucine-rich repeats) domain of the NLRP1 protein (NACHT, LRR and PYD domains-containing protein 1). This homozygous mutation was shared by two sisters with different clinical presentation: the younger sister had generalized inflammatory nodules with keratotic plugs, clinically resembling multiple keratoacanthomas, while the older had manifestations of familial keratosis lichenoides chronica. OBJECTIVES: To analyse the consequences of this NLRP1 variant in two siblings with a different clinical spectrum of severity. METHODS: To demonstrate the pathogenicity, p.Leu813Pro was recombinantly expressed, and its effect on inflammasome assembly was assessed. Exome sequencing and RNA-Seq were performed to identify factors with potentially modifying effects on the severity of the skin manifestation between each sibling. RESULTS: The variant p.Leu813Pro triggered activation of the NLRP1 inflammasome leading to ASC (apoptosis-associated speck-like protein containing a CARD) speck formation and interleukin (IL)-1ß release. The more severely affected sister had several additional genomic variants associated with atopy and psoriasis that were not present in her sibling. IL-5 and IL-17 emerged as dominant cytokines driving prominent inflammation in the skin of the severely affected sibling. CONCLUSIONS: To the best of our knowledge, this is the first report of a NLRP1 variant that leads to a different clinical spectrum of severity within the same sibship. IL-5 and IL-17 were the main cytokines expressed in the inflammatory lesions of the severely affected patient and might be regarded as disease modifying factors, and therefore may be considered as therapeutic targets.

A progeroid syndrome caused by a deep intronic variant in TAPT1 is revealed by RNA/SI-NET sequencing.

Exome sequencing has introduced a paradigm shift for the identification of germline variations responsible for Mendelian diseases. However, non-coding regions, which make up 98% of the genome, cannot be captured. The lack of functional annotation for intronic and intergenic variants makes RNA-seq a powerful companion diagnostic. Here, we illustrate this point by identifying six patients with a recessive Osteogenesis Imperfecta (OI) and neonatal progeria syndrome. By integrating homozygosity mapping and RNA-seq, we delineated a deep intronic TAPT1 mutation (c.1237-52 G>A) that segregated with the disease. Using SI-NET-seq, we document that TAPT1's nascent transcription was not affected in patients' fibroblasts, indicating instead that this variant leads to an alteration of pre-mRNA processing. Predicted to serve as an alternative splicing branchpoint, this mutation enhances TAPT1 exon 12 skipping, creating a protein-null allele. Additionally, our study reveals dysregulation of pathways involved in collagen and extracellular matrix biology in disease-relevant cells. Overall, our work highlights the power of transcriptomic approaches in deciphering the repercussions of non-coding variants, as well as in illuminating the molecular mechanisms of human diseases.

Omics profiling identifies the regulatory functions of the MAPK/ERK pathway in nephron progenitor metabolism.

Nephron endowment is defined by fetal kidney growth and crucially dictates renal health in adults. Defects in the molecular regulation of nephron progenitors contribute to only a fraction of reduced nephron mass cases, suggesting alternative causative mechanisms. The importance of MAPK/ERK activation in nephron progenitor maintenance has been previously demonstrated, and here, we characterized the metabolic consequences of MAPK/ERK deficiency. Liquid chromatography/mass spectrometry-based metabolomics profiling identified 42 reduced metabolites, of which 26 were supported by in vivo transcriptional changes in MAPK/ERK-deficient nephron progenitors. Among these, mitochondria, ribosome and amino acid metabolism, together with diminished pyruvate and proline metabolism, were the most affected pathways. In vitro cultures of mouse kidneys demonstrated a dosage-specific function for pyruvate in controlling the shape of the ureteric bud tip, a regulatory niche for nephron progenitors. In vivo disruption of proline metabolism caused premature nephron progenitor exhaustion through their accelerated differentiation in pyrroline-5-carboxylate reductases 1 (Pycr1) and 2 (Pycr2) double-knockout kidneys. Pycr1/Pycr2-deficient progenitors showed normal cell survival, indicating no changes in cellular stress. Our results suggest that MAPK/ERK-dependent metabolism functionally participates in nephron progenitor maintenance by monitoring pyruvate and proline biogenesis in developing kidneys.

INTS13 variants causing a recessive developmental ciliopathy disrupt assembly of the Integrator complex.

Oral-facial-digital (OFD) syndromes are a heterogeneous group of congenital disorders characterized by malformations of the face and oral cavity, and digit anomalies. Mutations within 12 cilia-related genes have been identified that cause several types of OFD, suggesting that OFDs constitute a subgroup of developmental ciliopathies. Through homozygosity mapping and exome sequencing of two families with variable OFD type 2, we identified distinct germline variants in INTS13, a subunit of the Integrator complex. This multiprotein complex associates with RNA Polymerase II and cleaves nascent RNA to modulate gene expression. We determined that INTS13 utilizes its C-terminus to bind the Integrator cleavage module, which is disrupted by the identified germline variants p.S652L and p.K668Nfs*9. Depletion of INTS13 disrupts ciliogenesis in human cultured cells and causes dysregulation of a broad collection of ciliary genes. Accordingly, its knockdown in Xenopus embryos leads to motile cilia anomalies. Altogether, we show that mutations in INTS13 cause an autosomal recessive ciliopathy, which reveals key interactions between components of the Integrator complex.

DPP9 deficiency: An inflammasomopathy that can be rescued by lowering NLRP1/IL-1 signaling.

Dipeptidyl peptidase 9 (DPP9) is a direct inhibitor of NLRP1, but how it affects inflammasome regulation in vivo is not yet established. Here, we report three families with immune-associated defects, poor growth, pancytopenia, and skin pigmentation abnormalities that segregate with biallelic DPP9 rare variants. Using patient-derived primary cells and biochemical assays, these variants were shown to behave as hypomorphic or knockout alleles that failed to repress NLRP1. The removal of a single copy of Nlrp1a/b/c, Asc, Gsdmd, or Il-1r, but not Il-18, was sufficient to rescue the lethality of Dpp9 mutant neonates in mice. Similarly, dpp9 deficiency was partially rescued by the inactivation of asc, an obligate downstream adapter of the NLRP1 inflammasome, in zebrafish. These experiments suggest that the deleterious consequences of DPP9 deficiency were mostly driven by the aberrant activation of the canonical NLRP1 inflammasome and IL-1ß signaling. Collectively, our results delineate a Mendelian disorder of DPP9 deficiency driven by increased NLRP1 activity as demonstrated in patient cells and in two animal models of the disease.

Loss of FOCAD, operating via the SKI messenger RNA surveillance pathway, causes a pediatric syndrome with liver cirrhosis.

Cirrhosis is usually a late-onset and life-threatening disease characterized by fibrotic scarring and inflammation that disrupts liver architecture and function. While it is typically the result of alcoholism or hepatitis viral infection in adults, its etiology in infants is much less understood. In this study, we report 14 children from ten unrelated families presenting with a syndromic form of pediatric liver cirrhosis. By genome/exome sequencing, we found recessive variants in FOCAD segregating with the disease. Zebrafish lacking focad phenocopied the human disease, revealing a signature of altered messenger RNA (mRNA) degradation processes in the liver. Using patient's primary cells and CRISPR-Cas9-mediated inactivation in human hepatic cell lines, we found that FOCAD deficiency compromises the SKI mRNA surveillance pathway by reducing the levels of the RNA helicase SKIC2 and its cofactor SKIC3. FOCAD knockout hepatocytes exhibited lowered albumin expression and signs of persistent injury accompanied by CCL2 overproduction. Our results reveal the importance of FOCAD in maintaining liver homeostasis and disclose a possible therapeutic intervention point via inhibition of the CCL2/CCR2 signaling axis.

R-SPONDIN2+ mesenchymal cells form the bud tip progenitor niche during human lung development.

The human respiratory epithelium is derived from a progenitor cell in the distal buds of the developing lung. These "bud tip progenitors" are regulated by reciprocal signaling with surrounding mesenchyme; however, mesenchymal heterogeneity and function in the developing human lung are poorly understood. We interrogated single-cell RNA sequencing data from multiple human lung specimens and identified a mesenchymal cell population present during development that is highly enriched for expression of the WNT agonist RSPO2, and we found that the adjacent bud tip progenitors are enriched for the RSPO2 receptor LGR5. Functional experiments using organoid models, explant cultures, and FACS-isolated RSPO2+ mesenchyme show that RSPO2 is a critical niche cue that potentiates WNT signaling in bud tip progenitors to support their maintenance and multipotency.

Direct identification of A-to-I editing sites with nanopore native RNA sequencing.

Inosine is a prevalent RNA modification in animals and is formed when an adenosine is deaminated by the ADAR family of enzymes. Traditionally, inosines are identified indirectly as variants from Illumina RNA-sequencing data because they are interpreted as guanosines by cellular machineries. However, this indirect method performs poorly in protein-coding regions where exons are typically short, in non-model organisms with sparsely annotated single-nucleotide polymorphisms, or in disease contexts where unknown DNA mutations are pervasive. Here, we show that Oxford Nanopore direct RNA sequencing can be used to identify inosine-containing sites in native transcriptomes with high accuracy. We trained convolutional neural network models to distinguish inosine from adenosine and guanosine, and to estimate the modification rate at each editing site. Furthermore, we demonstrated their utility on the transcriptomes of human, mouse and Xenopus. Our approach expands the toolkit for studying adenosine-to-inosine editing and can be further extended to investigate other RNA modifications.

RABENOSYN separation-of-function mutations uncouple endosomal recycling from lysosomal degradation, causing a distinct Mendelian disorder.

Rabenosyn (RBSN) is a conserved endosomal protein necessary for regulating internalized cargo. Here, we present clinical, genetic, cellular and biochemical evidence that two distinct RBSN missense variants are responsible for a novel Mendelian disorder consisting of progressive muscle weakness, facial dysmorphisms, ophthalmoplegia and intellectual disability. Using exome sequencing, we identified recessively acting germline alleles p.Arg180Gly and p.Gly183Arg, which are both situated in the FYVE domain of RBSN. We find that these variants abrogate binding to its cognate substrate phosphatidylinositol 3-phosphate (PI3P) and thus prevent its translocation to early endosomes. Although the endosomal recycling pathway was unaltered, mutant p.Gly183Arg patient fibroblasts show accumulation of cargo tagged for lysosomal degradation. Our results suggest that these variants are separation-of-function alleles, which cause a delay in endosomal maturation without affecting cargo recycling. We conclude that distinct germline mutations in RBSN cause non-overlapping phenotypes with specific and discrete endolysosomal cellular defects.