Next-generation phenotyping integrated in a national framework for patients with ultrarare disorders improves genetic diagnostics and yields new molecular findings.

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams.

GestaltMatcher Database - A global reference for facial phenotypic variability in rare human diseases.

The most important factor that complicates the work of dysmorphologists is the significant phenotypic variability of the human face. Next-Generation Phenotyping (NGP) tools that assist clinicians with recognizing characteristic syndromic patterns are particularly challenged when confronted with patients from populations different from their training data. To that end, we systematically analyzed the impact of genetic ancestry on facial dysmorphism. For that purpose, we established the GestaltMatcher Database (GMDB) as a reference dataset for medical images of patients with rare genetic disorders from around the world. We collected 10,980 frontal facial images - more than a quarter previously unpublished - from 8,346 patients, representing 581 rare disorders. Although the predominant ancestry is still European (67%), data from underrepresented populations have been increased considerably via global collaborations (19% Asian and 7% African). This includes previously unpublished reports for more than 40% of the African patients. The NGP analysis on this diverse dataset revealed characteristic performance differences depending on the composition of training and test sets corresponding to genetic relatedness. For clinical use of NGP, incorporating non-European patients resulted in a profound enhancement of GestaltMatcher performance. The top-5 accuracy rate increased by +11.29%. Importantly, this improvement in delineating the correct disorder from a facial portrait was achieved without decreasing the performance on European patients. By design, GMDB complies with the FAIR principles by rendering the curated medical data findable, accessible, interoperable, and reusable. This means GMDB can also serve as data for training and benchmarking. In summary, our study on facial dysmorphism on a global sample revealed a considerable cross ancestral phenotypic variability confounding NGP that should be counteracted by international efforts for increasing data diversity. GMDB will serve as a vital reference database for clinicians and a transparent training set for advancing NGP technology.

REEV: review, evaluate and explain variants.

In the era of high throughput sequencing, special software is required for the clinical evaluation of genetic variants. We developed REEV (Review, Evaluate and Explain Variants), a user-friendly platform for clinicians and researchers in the field of rare disease genetics. Supporting data was aggregated from public data sources. We compared REEV with seven other tools for clinical variant evaluation. REEV (semi-)automatically fills individual ACMG criteria facilitating variant interpretation. REEV can store disease and phenotype data related to a case to use these for phenotype similarity measures. Users can create public permanent links for individual variants that can be saved as browser bookmarks and shared. REEV may help in the fast diagnostic assessment of genetic variants in a clinical as well as in a research context. REEV (https://reev.bihealth.org/) is free and open to all users and there is no login requirement.

Rare Genetic Developmental Disabilities: Mabry Syndrome (MIM 239300) Index Cases and Glycophosphatidylinositol (GPI) Disorders.

The case report by Mabry et al. (1970) of a family with four children with elevated tissue non-specific alkaline phosphatase, seizures and profound developmental disability, became the basis for phenotyping children with the features that became known as Mabry syndrome. Aside from improvements in the services available to patients and families, however, the diagnosis and treatment of this, and many other developmental disabilities, did not change significantly until the advent of massively parallel sequencing. As more patients with features of the Mabry syndrome were identified, exome and genome sequencing were used to identify the glycophosphatidylinositol (GPI) biosynthesis disorders (GPIBDs) as a group of congenital disorders of glycosylation (CDG). Biallelic variants of the phosphatidylinositol glycan (PIG) biosynthesis, type V (PIGV) gene identified in Mabry syndrome became evidence of the first in a phenotypic series that is numbered HPMRS1-6 in the order of discovery. HPMRS1 [MIM: 239300] is the phenotype resulting from inheritance of biallelic PIGV variants. Similarly, HPMRS2 (MIM 614749), HPMRS5 (MIM 616025) and HPMRS6 (MIM 616809) result from disruption of the PIGO, PIGW and PIGY genes expressed in the endoplasmic reticulum. By contrast, HPMRS3 (MIM 614207) and HPMRS4 (MIM 615716) result from disruption of post attachment to proteins PGAP2 (HPMRS3) and PGAP3 (HPMRS4). The GPI biosynthesis disorders (GPIBDs) are currently numbered GPIBD1-21. Working with Dr. Mabry, in 2020, we were able to use improved laboratory diagnostics to complete the molecular diagnosis of patients he had originally described in 1970. We identified biallelic variants of the PGAP2 gene in the first reported HPMRS patients. We discuss the longevity of the Mabry syndrome index patients in the context of the utility of pyridoxine treatment of seizures and evidence for putative glycolipid storage in patients with HPMRS3. From the perspective of the laboratory innovations made that enabled the identification of the HPMRS phenotype in Dr. Mabry's patients, the need for treatment innovations that will benefit patients and families affected by developmental disabilities is clear.

Next-generation phenotyping in Nigerian children with Cornelia de Lange syndrome.

Next-generation phenotyping (NGP) can be used to compute the similarity of dysmorphic patients to known syndromic diseases. So far, the technology has been evaluated in variant prioritization and classification, providing evidence for pathogenicity if the phenotype matched with other patients with a confirmed molecular diagnosis. In a Nigerian cohort of individuals with facial dysmorphism, we used the NGP tool GestaltMatcher to screen portraits prior to genetic testing and subjected individuals with high similarity scores to exome sequencing (ES). Here, we report on two individuals with global developmental delay, pulmonary artery stenosis, and genital and limb malformations for whom GestaltMatcher yielded Cornelia de Lange syndrome (CdLS) as the top hit. ES revealed a known pathogenic nonsense variant, NM_133433.4: c.598C>T; p.(Gln200*), as well as a novel frameshift variant c.7948dup; p.(Ile2650Asnfs*11) in NIPBL. Our results suggest that NGP can be used as a screening tool and thresholds could be defined for achieving high diagnostic yields in ES. Training the artificial intelligence (AI) with additional cases of the same ethnicity might further increase the positive predictive value of GestaltMatcher.

Enhancing Variant Prioritization in VarFish through On-Premise Computational Facial Analysis.

Genomic variant prioritization is crucial for identifying disease-associated genetic variations. Integrating facial and clinical feature analyses into this process enhances performance. This study demonstrates the integration of facial analysis (GestaltMatcher) and Human Phenotype Ontology analysis (CADA) within VarFish, an open-source variant analysis framework. Challenges related to non-open-source components were addressed by providing an open-source version of GestaltMatcher, facilitating on-premise facial analysis to address data privacy concerns. Performance evaluation on 163 patients recruited from a German multi-center study of rare diseases showed PEDIA's superior accuracy in variant prioritization compared to individual scores. This study highlights the importance of further benchmarking and future integration of advanced facial analysis approaches aligned with ACMG guidelines to enhance variant classification.

GestaltMatcher Database - A global reference for the facial phenotypic variability of rare human diseases.

Dysmorphologists sometimes encounter challenges in recognizing disorders due to phenotypic variability influenced by factors such as age and ethnicity. Moreover, the performance of Next Generation Phenotyping Tools such as GestaltMatcher is dependent on the diversity of the training set. Therefore, we developed GestaltMatcher Database (GMDB) - a global reference for the phenotypic variability of rare diseases that complies with the FAIR-principles. We curated dysmorphic patient images and metadata from 2,224 publications, transforming GMDB into an online dynamic case report journal. To encourage clinicians worldwide to contribute, each case can receive a Digital Object Identifier (DOI), making it a citable micro-publication. This resulted in a collection of 2,312 unpublished images, partly with longitudinal data. We have compiled a collection of 10,189 frontal images from 7,695 patients representing 683 disorders. The web interface enables gene- and phenotype-centered queries for registered users (https://db.gestaltmatcher.org/). Despite the predominant European ancestry of most patients (59%), our global collaborations have facilitated the inclusion of data from frequently underrepresented ethnicities, with 17% Asian, 4% African, and 6% with other ethnic backgrounds. The analysis has revealed a significant enhancement in GestaltMatcher performance across all ethnic groups, incorporating non-European ethnicities, showcasing a remarkable increase in Top-1-Accuracy by 31.56% and Top-5-Accuracy by 12.64%. Importantly, this improvement was achieved without altering the performance metrics for European patients. GMDB addresses dysmorphology challenges by representing phenotypic variability and including underrepresented groups, enhancing global diagnostic rates and serving as a vital clinician reference database.

Biallelic variants in PIGN cause Fryns syndrome, multiple congenital anomalies-hypotonia-seizures syndrome, and neurologic phenotypes: A genotype-phenotype correlation study.

PURPOSE: Biallelic PIGN variants have been described in Fryns syndrome, multiple congenital anomalies-hypotonia-seizure syndrome (MCAHS), and neurologic phenotypes. The full spectrum of clinical manifestations in relation to the genotypes is yet to be reported. METHODS: Genotype and phenotype data were collated and analyzed for 61 biallelic PIGN cases: 21 new and 40 previously published cases. Functional analysis was performed for 2 recurrent variants (c.2679C>G p.Ser893Arg and c.932T>G p.Leu311Trp). RESULTS: Biallelic-truncating variants were detected in 16 patients-10 with Fryns syndrome, 1 with MCAHS1, 2 with Fryns syndrome/MCAHS1, and 3 with neurologic phenotype. There was an increased risk of prenatal or neonatal death within this group (6 deaths were in utero or within 2 months of life; 6 pregnancies were terminated). Incidence of polyhydramnios, congenital anomalies (eg, diaphragmatic hernia), and dysmorphism was significantly increased. Biallelic missense or mixed genotype were reported in the remaining 45 cases-32 showed a neurologic phenotype and 12 had MCAHS1. No cases of diaphragmatic hernia or abdominal wall defects were seen in this group except patient 1 in which we found the missense variant p.Ser893Arg to result in functionally null alleles, suggesting the possibility of an undescribed functionally important region in the final exon. For all genotypes, there was complete penetrance for developmental delay and near-complete penetrance for seizures and hypotonia in patients surviving the neonatal period. CONCLUSION: We have expanded the described spectrum of phenotypes and natural history associated with biallelic PIGN variants. Our study shows that biallelic-truncating variants usually result in the more severe Fryns syndrome phenotype, but neurologic problems, such as developmental delay, seizures, and hypotonia, present across all genotypes. Functional analysis should be considered when the genotypes do not correlate with the predicted phenotype because there may be other functionally important regions in PIGN that are yet to be discovered.

Next-generation phenotyping contributing to the identification of a 4.7 kb deletion in KANSL1 causing Koolen-de Vries syndrome.

Next-generation phenotyping (NGP) is an application of advanced methods of computer vision on medical imaging data such as portrait photos of individuals with rare disorders. NGP on portraits results in gestalt scores that can be used for the selection of appropriate genetic tests, and for the interpretation of the molecular data. Here, we report on an exceptional case of a young girl that was presented at the age of 8 and 15 and enrolled in NGP diagnostics on the latter occasion. The girl had clinical features associated with Koolen-de Vries syndrome (KdVS) and a suggestive facial gestalt. However, chromosomal microarray (CMA), Sanger sequencing, multiplex ligation-dependent probe analysis (MLPA), and trio exome sequencing remained inconclusive. Based on the highly indicative gestalt score for KdVS, the decision was made to perform genome sequencing to also evaluate noncoding variants. This analysis revealed a 4.7 kb de novo deletion partially affecting intron 6 and exon 7 of the KANSL1 gene. This is the smallest reported structural variant to date for this phenotype. The case illustrates how NGP can be integrated into the iterative diagnostic process of test selection and interpretation of sequencing results.

GestaltMatcher facilitates rare disease matching using facial phenotype descriptors.

Many monogenic disorders cause a characteristic facial morphology. Artificial intelligence can support physicians in recognizing these patterns by associating facial phenotypes with the underlying syndrome through training on thousands of patient photographs. However, this 'supervised' approach means that diagnoses are only possible if the disorder was part of the training set. To improve recognition of ultra-rare disorders, we developed GestaltMatcher, an encoder for portraits that is based on a deep convolutional neural network. Photographs of 17,560 patients with 1,115 rare disorders were used to define a Clinical Face Phenotype Space, in which distances between cases define syndromic similarity. Here we show that patients can be matched to others with the same molecular diagnosis even when the disorder was not included in the training set. Together with mutation data, GestaltMatcher could not only accelerate the clinical diagnosis of patients with ultra-rare disorders and facial dysmorphism but also enable the delineation of new phenotypes.

PIGN encephalopathy: Characterizing the epileptology.

OBJECTIVE: Epilepsy is common in patients with PIGN diseases due to biallelic variants; however, limited epilepsy phenotyping data have been reported. We describe the epileptology of PIGN encephalopathy. METHODS: We recruited patients with epilepsy due to biallelic PIGN variants and obtained clinical data regarding age at seizure onset/offset and semiology, development, medical history, examination, electroencephalogram, neuroimaging, and treatment. Seizure and epilepsy types were classified. RESULTS: Twenty six patients (13 female) from 26 families were identified, with mean age 7 years (range = 1 month to 21 years; three deceased). Abnormal development at seizure onset was present in 25 of 26. Developmental outcome was most frequently profound (14/26) or severe (11/26). Patients presented with focal motor (12/26), unknown onset motor (5/26), focal impaired awareness (1/26), absence (2/26), myoclonic (2/26), myoclonic-atonic (1/26), and generalized tonic-clonic (2/26) seizures. Twenty of 26 were classified as developmental and epileptic encephalopathy (DEE): 55% (11/20) focal DEE, 30% (6/20) generalized DEE, and 15% (3/20) combined DEE. Six had intellectual disability and epilepsy (ID+E): two generalized and four focal epilepsy. Mean age at seizure onset was 13 months (birth to 10 years), with a lower mean onset in DEE (7 months) compared with ID+E (33 months). Patients with DEE had drug-resistant epilepsy, compared to 4/6 ID+E patients, who were seizure-free. Hyperkinetic movement disorder occurred in 13 of 26 patients. Twenty-seven of 34 variants were novel. Variants were truncating (n = 7), intronic and predicted to affect splicing (n = 7), and missense or inframe indels (n = 20, of which 11 were predicted to affect splicing). Seven variants were recurrent, including p.Leu311Trp in 10 unrelated patients, nine with generalized seizures, accounting for nine of the 11 patients in this cohort with generalized seizures. SIGNIFICANCE: PIGN encephalopathy is a complex autosomal recessive disorder associated with a wide spectrum of epilepsy phenotypes, typically with substantial profound to severe developmental impairment.

Induction of Rosette-to-Lumen stage embryoids using reprogramming paradigms in ESCs.

Blastocyst-derived stem cell lines were shown to self-organize into embryo-like structures in 3D cell culture environments. Here, we provide evidence that embryo-like structures can be generated solely based on transcription factor-mediated reprogramming of embryonic stem cells in a simple 3D co-culture system. Embryonic stem cells in these cultures self-organize into elongated, compartmentalized embryo-like structures reflecting aspects of the inner regions of the early post-implantation embryo. Single-cell RNA-sequencing reveals transcriptional profiles resembling epiblast, primitive-/visceral endoderm, and extraembryonic ectoderm of early murine embryos around E4.5-E5.5. In this stem cell-based embryo model, progression from rosette formation to lumenogenesis accompanied by progression from naïve- to primed pluripotency was observed within Epi-like cells. Additionally, lineage specification of primordial germ cells and distal/anterior visceral endoderm-like cells was observed in epiblast- or visceral endoderm-like compartments, respectively. The system presented in this study allows for fast and reproducible generation of embryo-like structures, providing an additional tool to study aspects of early embryogenesis.

CADA: phenotype-driven gene prioritization based on a case-enriched knowledge graph.

Many rare syndromes can be well described and delineated from other disorders by a combination of characteristic symptoms. These phenotypic features are best documented with terms of the Human Phenotype Ontology (HPO), which are increasingly used in electronic health records (EHRs), too. Many algorithms that perform HPO-based gene prioritization have also been developed; however, the performance of many such tools suffers from an over-representation of atypical cases in the medical literature. This is certainly the case if the algorithm cannot handle features that occur with reduced frequency in a disorder. With Cada, we built a knowledge graph based on both case annotations and disorder annotations. Using network representation learning, we achieve gene prioritization by link prediction. Our results suggest that Cada exhibits superior performance particularly for patients that present with the pathognomonic findings of a disease. Additionally, information about the frequency of occurrence of a feature can readily be incorporated, when available. Crucial in the design of our approach is the use of the growing amount of phenotype-genotype information that diagnostic labs deposit in databases such as ClinVar. By this means, Cada is an ideal reference tool for differential diagnostics in rare disorders that can also be updated regularly.

TBK1 and TNFRSF13B mutations and an autoinflammatory disease in a child with lethal COVID-19.

Among children, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections are typically mild. Here, we describe the case of a 3.5-year-old girl with an unusually severe presentation of coronavirus disease (COVID-19). The child had an autoinflammatory disorder of unknown etiology, which had been treated using prednisolone and methotrexate, and her parents were half cousins of Turkish descent. After 5 days of nonspecific viral infection symptoms, tonic-clonic seizures occurred followed by acute cardiac insufficiency, multi-organ insufficiency, and ultimate death. Trio exome sequencing identified a homozygous splice-variant in the gene TBK1, and a homozygous missense variant in the gene TNFRSF13B. Heterozygous deleterious variants in the TBK1 gene have been associated with severe COVID-19, and the variant in the TNFRSF13B gene has been associated with common variable immunodeficiency (CVID). We suggest that the identified variants, the autoinflammatory disorder and its treatment, or a combination of these factors probably predisposed to lethal COVID-19 in the present case.

A CRISPR-Cas9-engineered mouse model for GPI-anchor deficiency mirrors human phenotypes and exhibits hippocampal synaptic dysfunctions.

Pathogenic germline mutations in PIGV lead to glycosylphosphatidylinositol biosynthesis deficiency (GPIBD). Individuals with pathogenic biallelic mutations in genes of the glycosylphosphatidylinositol (GPI)-anchor pathway exhibit cognitive impairments, motor delay, and often epilepsy. Thus far, the pathophysiology underlying the disease remains unclear, and suitable rodent models that mirror all symptoms observed in human patients have not been available. Therefore, we used CRISPR-Cas9 to introduce the most prevalent hypomorphic missense mutation in European patients, Pigv:c.1022C > A (p.A341E), at a site that is conserved in mice. Mirroring the human pathology, mutant Pigv341E mice exhibited deficits in motor coordination, cognitive impairments, and alterations in sociability and sleep patterns, as well as increased seizure susceptibility. Furthermore, immunohistochemistry revealed reduced synaptophysin immunoreactivity in Pigv341E mice, and electrophysiology recordings showed decreased hippocampal synaptic transmission that could underlie impaired memory formation. In single-cell RNA sequencing, Pigv341E-hippocampal cells exhibited changes in gene expression, most prominently in a subtype of microglia and subicular neurons. A significant reduction in Abl1 transcript levels in several cell clusters suggested a link to the signaling pathway of GPI-anchored ephrins. We also observed elevated levels of Hdc transcripts, which might affect histamine metabolism with consequences for circadian rhythm. This mouse model will not only open the doors to further investigation into the pathophysiology of GPIBD, but will also deepen our understanding of the role of GPI-anchor-related pathways in brain development.

Evidence of the milder phenotypic spectrum of c.1582G>A PIGT variant: Delineation based on seven novel Polish patients.

PIGT is one of over 29 glycosylphosphatidylinositol biosynthesis defect genes. Mutations cause genetically determined disorders characterized mainly by epilepsy with fever-sensitivity, central hypotonia, psychomotor delay and congenital malformations. The disease is known as multiple congenital anomalies-hypotonia-seizures syndrome 3 (MCAHS3) or glycosylphosphatidylinositol biosynthesis defect-7. Twenty-eight cases have been reported until today. We present seven novel Polish patients, all harboring 1582G>A variant in a homozygous or compound heterozygous state which seems to cause a milder phenotype of the disease.

Lessons learned from 40 novel PIGA patients and a review of the literature.

OBJECTIVE: To define the phenotypic spectrum of phosphatidylinositol glycan class A protein (PIGA)-related congenital disorder of glycosylation (PIGA-CDG) and evaluate genotype-phenotype correlations. METHODS: Our cohort encompasses 40 affected males with a pathogenic PIGA variant. We performed a detailed phenotypic assessment, and in addition, we reviewed the available clinical data of 36 previously published cases and assessed the variant pathogenicity using bioinformatical approaches. RESULTS: Most individuals had hypotonia, moderate to profound global developmental delay, and intractable seizures. We found that PIGA-CDG spans from a pure neurological phenotype at the mild end to a Fryns syndrome-like phenotype. We found a high frequency of cardiac anomalies including structural anomalies and cardiomyopathy, and a high frequency of spontaneous death, especially in childhood. Comparative bioinformatical analysis of common variants, found in the healthy population, and pathogenic variants, identified in affected individuals, revealed a profound physiochemical dissimilarity of the substituted amino acids in variant constrained regions of the protein. SIGNIFICANCE: Our comprehensive analysis of the largest cohort of published and novel PIGA patients broadens the spectrum of PIGA-CDG. Our genotype-phenotype correlation facilitates the estimation on pathogenicity of variants with unknown clinical significance and prognosis for individuals with pathogenic variants in PIGA.

VarFish: comprehensive DNA variant analysis for diagnostics and research.

VarFish is a user-friendly web application for the quality control, filtering, prioritization, analysis, and user-based annotation of DNA variant data with a focus on rare disease genetics. It is capable of processing variant call files with single or multiple samples. The variants are automatically annotated with population frequencies, molecular impact, and presence in databases such as ClinVar. Further, it provides support for pathogenicity scores including CADD, MutationTaster, and phenotypic similarity scores. Users can filter variants based on these annotations and presumed inheritance pattern and sort the results by these scores. Variants passing the filter are listed with their annotations and many useful link-outs to genome browsers, other gene/variant data portals, and external tools for variant assessment. VarFish allows users to create their own annotations including support for variant assessment following ACMG-AMP guidelines. In close collaboration with medical practitioners, VarFish was designed for variant analysis and prioritization in diagnostic and research settings as described in the software's extensive manual. The user interface has been optimized for supporting these protocols. Users can install VarFish on their own in-house servers where it provides additional lab notebook features for collaborative analysis and allows re-analysis of cases, e.g. after update of genotype or phenotype databases.

A Novel Mutation in PIGA Associated with Multiple Congenital Anomalies-Hypotonia-Seizure Syndrome 2 (MCAHS2) in a Boy with a Combination of Severe Epilepsy and Gingival Hyperplasia.

Multiple congenital anomalies-hypotonia-seizures syndrome 2 (MCAHS2) is a rare disease caused by mutations in the X chromosomal PIGA gene. Clinically it is characterized by early-onset epilepsy, hypotonia, dysmorphic features, and variable congenital anomalies. PIGA codes for the phosphatidylinositol glycan-class A protein, which forms a subunit of an enzymatic complex involved in glycophosphatidylinositol (GPI) biosynthesis. We present a new case of MCAHS2 and perform a comprehensive review of the available literature to delineate the phenotypical traits associated with germline PIGA mutations. Furthermore, we provide functional evidence of pathogenicity of the novel missense mutation, c.154C>T; (p.His52Tyr), in the PIGA gene causative of MCAHS2 in our patient. By flow cytometry, we observed reduced expression of GPI-anchored surface proteins in patient granulocytes compared to control samples, proving GPI-biogenesis impairment. The patient's severe epilepsy with several daily attacks was refractory to treatment, but the frequency of seizures reduced temporarily under triple therapy with perampanel, rufinamide and vigabatrin. Our study delineates the known MCAHS2 phenotype and discusses challenges of diagnosis and clinical management in this complex, rare disease. Furthermore, we present a novel mutation with functional evidence of pathogenicity.