Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability

The origin and evolution of Wnt signalling

The Wnt signal transduction pathway has essential roles in the formation of the primary body axis during development, cellular differentiation and tissue homeostasis. This animal-specific pathway has been studied extensively in contexts ranging from developmental biology to medicine for more than 40 years. Despite its physiological importance, an understanding of the evolutionary origin and primary function of Wnt signalling has begun to emerge only recently. Recent studies on very basal metazoan species have shown high levels of conservation of components of both canonical and non-canonical Wnt signalling pathways. Furthermore, some pathway proteins have been described also in non-animal species, suggesting that recruitment and functional adaptation of these factors has occurred in metazoans. In this Review, we summarize the current state of research regarding the evolutionary origin of Wnt signalling, its ancestral function and the characteristics of the primal Wnt ligand, with emphasis on the importance of genomic studies in various pre-metazoan and basal metazoan species.

Craniofacial studies in chicken embryos confirm the pathogenicity of human FZD2 variants associated with Robinow syndrome

Robinow syndrome is a rare disease caused by variants of seven WNT pathway genes. Craniofacial features include widening of the nasal bridge and jaw hypoplasia. We used the chicken embryo to test whether two missense human FZD2 variants (1301G>T, p.Gly434Val; 425C>T, p.Pro142Lys) were sufficient to change frontonasal mass development. In vivo, the overexpression of retroviruses with wild-type or variant human FZD2 inhibited upper beak ossification. In primary cultures, wild-type and variant human FZD2 significantly inhibited chondrogenesis, with the 425C>T variant significantly decreasing activity of a SOX9 luciferase reporter compared to that for the wild type or 1301G>T. Both variants also increased nuclear shuttling of β-catenin (CTNNB1) and increased the expression of TWIST1, which are inhibitory to chondrogenesis. In canonical WNT luciferase assays using frontonasal mass cells, the variants had dominant-negative effects on wild-type FZD2. In non-canonical assays, the 425C>T variant failed to activate the reporter above control levels and was unresponsive to exogenous WNT5A. This is the first single amino acid change to selectively alter ligand binding in a FZD receptor. Therefore, FZD2 missense variants are pathogenic and could lead to the altered craniofacial morphogenesis seen in Robinow syndrome.

NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy

BACKGROUND

NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects.

METHODS

We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences.

RESULTS

Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model.

CONCLUSION

Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

ROR2-Related Skeletal Dysplasia Reveals Disrupted Chondrocyte Polarity through Modulation of BMP/TGF-β Signaling

Genetic studies have shown that Robinow syndrome (RS), a rare skeletal dysplasia, is caused by ROR2 mutation. However, the cell origin and molecular mechanisms underlying this disease remain elusive. We established a conditional knockout system by crossing Prx1cre and Osxcre with Ror2 flox/flox mice. and conducted histological and immunofluorescence analyses to investigate the phenotypes during skeletal development. In the Prx1cre line, we observed RS-like skeletal abnormities, including short stature and an arched skull. Additionally, we found inhibition of chondrocyte differentiation and proliferation. In the Osxcre line, loss of ROR2 in osteoblast lineage cells led to reduced osteoblast differentiation during both embryonic and postnatal stages. Furthermore, ROR2 mutant mice exhibited increased adipogenesis in the bone marrow compared to their littermate controls. To further explore the underlying mechanisms, bulk RNA-seq analysis of Prx1cre; Ror2 flox/flox embryos was performed, results revealed decreased BMP/TGF-β signaling. Immunofluorescence analysis further confirmed the decreased expression of p-smad1/5/8, accompanied by disrupted cell polarity in the developing growth plate. Pharmacological treatment using FK506 partially rescued the skeletal dysplasia and resulted in increased mineralization and osteoblast differentiation. By modeling the phenotype of RS in mice, our findings provide evidence for the involvement of mesenchymal progenitors as the cell origin and highlight the molecular mechanism of BMP/TGF-β signaling in skeletal dysplasia.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

FZD2 regulates limb development by mediating β-catenin-dependent and -independent Wnt signaling pathways

Human Robinow syndrome (RS) and dominant omodysplasia type 2 (OMOD2), characterized by skeletal limb and craniofacial defects, are associated with heterozygous mutations in the Wnt receptor FZD2. However, as FZD2 can activate both canonical and non-canonical Wnt pathways, its precise functions and mechanisms of action in limb development are unclear. To address these questions, we generated mice harboring a single-nucleotide insertion in Fzd2 (Fzd2em1Smill), causing a frameshift mutation in the final Dishevelled-interacting domain. Fzd2em1Smill mutant mice had shortened limbs, resembling those of RS and OMOD2 patients, indicating that FZD2 mutations are causative. Fzd2em1Smill mutant embryos displayed decreased canonical Wnt signaling in developing limb mesenchyme and disruption of digit chondrocyte elongation and orientation, which is controlled by the β-catenin-independent WNT5A/planar cell polarity (PCP) pathway. In line with these observations, we found that disruption of FZD function in limb mesenchyme caused formation of shortened bone elements and defects in Wnt/β-catenin and WNT5A/PCP signaling. These findings indicate that FZD2 controls limb development by mediating both canonical and non-canonical Wnt pathways and reveal causality of pathogenic FZD2 mutations in RS and OMOD2 patients.

Successful therapeutic intervention in new mouse models of frizzled 2-associated congenital malformations

Frizzled 2 (FZD2) is a transmembrane Wnt receptor. We previously identified a pathogenic human FZD2 variant in individuals with FZD2-associated autosomal dominant Robinow syndrome. The variant encoded a protein with a premature stop and loss of 17 amino acids, including a region of the consensus dishevelled-binding sequence. To model this variant, we used zygote microinjection and i-GONAD-based CRISPR/Cas9-mediated genome editing to generate a mouse allelic series. Embryos mosaic for humanized Fzd2W553* knock-in exhibited cleft palate and shortened limbs, consistent with patient phenotypes. We also generated two germline mouse alleles with small deletions: Fzd2D3 and Fzd2D4. Homozygotes for each allele exhibit a highly penetrant cleft palate phenotype, shortened limbs compared with wild type and perinatal lethality. Fzd2D4 craniofacial tissues indicated decreased canonical Wnt signaling. In utero treatment with IIIC3a (a DKK inhibitor) normalized the limb lengths in Fzd2D4 homozygotes. The in vivo replication represents an approach for further investigating the mechanism of FZD2 phenotypes and demonstrates the utility of CRISPR knock-in mice as a tool for investigating the pathogenicity of human genetic variants. We also present evidence for a potential therapeutic intervention.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Obstetrical Challenges in Robinow Syndrome

Robinow syndrome is a genetically heterogenous syndrome that exhibits great pleiotropy, involving skeletal genital, cardiac, and craniofacial developmental anomalies. Fertility is not always compromised, and many individuals may be able to have a healthy pregnancy. Similar to other more common skeletal dysplasias and growth disorders such as achondroplasia, there are several challenges to be addressed in managing physiologic differences that occur in the context of pregnancy, and published literature centers on pregnant people with achondroplasia. We present a patient with Robinow syndrome (ROR2 variant), follow her clinical course through three of her pregnancies (one 20-week loss followed by two preterm cesarean deliveries at 36-week gestation), and highlight the major obstetrical considerations in her individualized care.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the non-canonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable AR trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome. This article is protected by copyright. All rights reserved.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.