Improved exome prioritization of disease genes through cross-species phenotype comparison

Addressing diagnostic gaps and priorities of the global rare diseases community: Recommendations from the IRDiRC diagnostics scientific committee

The International Rare Diseases Research Consortium (IRDiRC) Diagnostic Scientific Committee (DSC) is charged with discussion and contribution to progress on diagnostic aspects of the IRDiRC core mission. Specifically, IRDiRC goals include timely diagnosis, use of globally coordinated diagnostic pipelines, and assessing the impact of rare diseases on affected individuals. As part of this mission, the DSC endeavored to create a list of research priorities to achieve these goals. We present a discussion of those priorities along with aspects of current, global rare disease needs and opportunities that support our prioritization. In support of this discussion, we also provide clinical vignettes illustrating real-world examples of diagnostic challenges.

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.

Lessons and pitfalls of whole genome sequencing

Whole-genome sequencing (WGS) has recently become the first-line genetic investigation for many suspected genetic neurological disorders. While its diagnostic capabilities are innumerable, as with any test, it has its limitations. Clinicians should be aware of where WGS is extremely reliable (detecting single-nucleotide variants), where its reliability is much improved (detecting copy number variants and small repeat expansions) and where it may miss/misinterpret a variant (large repeat expansions, balanced structural variants or low heteroplasmy mitochondrial DNA variants). Bioinformatic technology and virtual gene panels are constantly evolving, and it is important to know what genes and what types of variant are being tested; the current National Health Service Genomic Medicine Service WGS offers more than early iterations of the 100 000 Genomes Project analysis. Close communication between clinician and laboratory, ideally through a multidisciplinary team meeting, is encouraged where there is diagnostic uncertainty.

Multi-gene panel sequencing in highly consanguineous families and patients with congenital forms of skeletal dysplasias

Skeletal dysplasias (SKDs) are a heterogeneous group of more than 750 genetic disorders characterized by abnormal development, growth, and maintenance of bones or cartilage in the human skeleton. SKDs are often caused by variants in early patterning genes and in many cases part of multiple malformation syndromes and occur in combination with non-skeletal phenotypes. The aim of this study was to investigate the underlying genetic cause of congenital SKDs in highly consanguineous Pakistani families, as well as in sporadic and familial SKD cases from India using multigene panel sequencing analysis. Therefore, we performed panel sequencing of 386 bone-related genes in 7 highly consanguineous families from Pakistan and 27 cases from India affected with SKDs. In the highly consanguineous families, we were able to identify the underlying genetic cause in five out of seven families, resulting in a diagnostic yield of 71%. Whereas, in the sporadic and familial SKD cases, we identified 12 causative variants, corresponding to a diagnostic yield of 44%. The genetic heterogeneity in our cohorts was very high and we were able to detect various types of variants, including missense, nonsense, and frameshift variants, across multiple genes known to cause different types of SKDs. In conclusion, panel sequencing proved to be a highly effective way to decipher the genetic basis of SKDs in highly consanguineous families as well as sporadic and or familial cases from South Asia. Furthermore, our findings expand the allelic spectrum of skeletal dysplasias.

Towards a standard benchmark for variant and gene prioritisation algorithms: PhEval - Phenotypic inference Evaluation framework

Background

Computational approaches to support rare disease diagnosis are challenging to build, requiring the integration of complex data types such as ontologies, gene-to-phenotype associations, and cross-species data into variant and gene prioritisation algorithms (VGPAs). However, the performance of VGPAs has been difficult to measure and is impacted by many factors, for example, ontology structure, annotation completeness or changes to the underlying algorithm. Assertions of the capabilities of VGPAs are often not reproducible, in part because there is no standardised, empirical framework and openly available patient data to assess the efficacy of VGPAs - ultimately hindering the development of effective prioritisation tools.

Results

In this paper, we present our benchmarking tool, PhEval, which aims to provide a standardised and empirical framework to evaluate phenotype-driven VGPAs. The inclusion of standardised test corpora and test corpus generation tools in the PhEval suite of tools allows open benchmarking and comparison of methods on standardised data sets.

Conclusions

PhEval and the standardised test corpora solve the issues of patient data availability and experimental tooling configuration when benchmarking and comparing rare disease VGPAs. By providing standardised data on patient cohorts from real-world case-reports and controlling the configuration of evaluated VGPAs, PhEval enables transparent, portable, comparable and reproducible benchmarking of VGPAs. As these tools are often a key component of many rare disease diagnostic pipelines, a thorough and standardised method of assessment is essential for improving patient diagnosis and care.

A corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

Prioritizing genomic variants through neuro-symbolic, knowledge-enhanced learning

MOTIVATION

Whole-exome and genome sequencing have become common tools in diagnosing patients with rare diseases. Despite their success, this approach leaves many patients undiagnosed. A common argument is that more disease variants still await discovery, or the novelty of disease phenotypes results from a combination of variants in multiple disease-related genes. Interpreting the phenotypic consequences of genomic variants relies on information about gene functions, gene expression, physiology, and other genomic features. Phenotype-based methods to identify variants involved in genetic diseases combine molecular features with prior knowledge about the phenotypic consequences of altering gene functions. While phenotype-based methods have been successfully applied to prioritizing variants, such methods are based on known gene-disease or gene-phenotype associations as training data and are applicable to genes that have phenotypes associated, thereby limiting their scope. In addition, phenotypes are not assigned uniformly by different clinicians, and phenotype-based methods need to account for this variability.

RESULTS

We developed an Embedding-based Phenotype Variant Predictor (EmbedPVP), a computational method to prioritize variants involved in genetic diseases by combining genomic information and clinical phenotypes. EmbedPVP leverages a large amount of background knowledge from human and model organisms about molecular mechanisms through which abnormal phenotypes may arise. Specifically, EmbedPVP incorporates phenotypes linked to genes, functions of gene products, and the anatomical site of gene expression, and systematically relates them to their phenotypic effects through neuro-symbolic, knowledge-enhanced machine learning. We demonstrate EmbedPVP's efficacy on a large set of synthetic genomes and genomes matched with clinical information.

AVAILABILITY AND IMPLEMENTATION

EmbedPVP and all evaluation experiments are freely available at https://github.com/bio-ontology-research-group/EmbedPVP.

Significance Associated with Phenotype Score Aids in Variant Prioritization for Exome Sequencing Analysis

Several in silico annotation-based methods have been developed to prioritize variants in exome sequencing analysis. This study introduced a novel metric Significance Associated with Phenotypes (SAP) score, which generates a statistical score by comparing an individual's observed phenotypes against existing gene-phenotype associations. To evaluate the SAP score, a retrospective analysis was performed on 219 exomes. Among them, 82 family-based and 35 singleton exomes had at least one disease-causing variant that explained the patient's clinical features. SAP scores were calculated, and the rank of the disease-causing variant was compared with a known method, Exomiser. Using the SAP score, the known causative variant was ranked in the top 10 retained variants for 94% (77 of 82) of the family-based exomes and in first place for 73% of these cases. For singleton exomes, the SAP score analysis ranked the known pathogenic variants within the top 10 for 80% (28 of 35) of cases. The SAP score, which is independent of detected variants, demonstrates comparable performance with Exomiser, which considers both phenotype and variant-level evidence simultaneously. Among 102 cases with negative results or variants of uncertain significance, SAP score analysis revealed two cases with a potential new diagnosis based on rank. The SAP score, a phenotypic quantitative metric, can be used in conjunction with standard variant filtration and annotation to enhance variant prioritization in exome analysis.

A de novo variant in PAK2 detected in an individual with Knobloch type 2 syndrome

P21-activated kinase 2 (PAK2) is a serine/threonine kinase essential for a variety of cellular processes including signal transduction, cellular survival, proliferation, and migration. A recent report proposed monoallelic PAK2 variants cause Knobloch syndrome type 2 (KNO2)-a developmental disorder primarily characterized by ocular anomalies. Here, we identified a novel de novo heterozygous missense variant in PAK2, NM_002577.4:c.1273G>A, p.(D425N), by whole genome sequencing in an individual with features consistent with KNO2. Notable clinical phenotypes include global developmental delay, congenital retinal detachment, mild cerebral ventriculomegaly, hypotonia, FTT, pyloric stenosis, feeding intolerance, patent ductus arteriosus, and mild facial dysmorphism. The p.(D425N) variant lies within the protein kinase domain and is predicted to be functionally damaging by in silico analysis. Previous clinical genetic testing did not report this variant due to unknown relevance of PAK2 variants at the time of testing, highlighting the importance of reanalysis. Our findings also substantiate the candidacy of PAK2 variants in KNO2 and expand the KNO2 clinical spectrum.

A Rare Inherited Bone Marrow Failure Syndrome Disclosed by Reanalysis of the Exome Data of a Patient Evaluated for Cytopenia and Dysmorphic Features

BACKGROUND

Multisystemic findings of inherited bone marrow failure syndromes may cause difficulty in diagnosis. Exome sequencing (ES) helps to define the etiology of rare diseases and reanalysis offers a valuable new diagnostic approach. Herein, we present the clinical and molecular characteristics of a girl who was referred for cytopenia and frequent infections.

CASE REPORT

A 5-year-old girl with cytopenia, dysmorphism, short stature, developmental delay, and myopia was referred for genetic counseling. Reanalysis of the ES data revealed a homozygous splice-site variant in the DNAJC21 (NM_001012339.3:c.983+1G>A), causing Shwachman-Diamond Syndrome (SDS). It was shown by the RNA sequencing that exon 7 was skipped, causing an 88-nucleotide deletion.

CONCLUSIONS

Precise genetic diagnosis enables genetic counseling and improves patient management by avoiding inappropriate treatment and unnecessary testing. This report would contribute to the clinical and molecular understanding of this rare type of SDS caused by DNAJC21 variants and expand the phenotypic features of this condition.

Universal Exome Sequencing in Critically Ill Adults: A Diagnostic Yield of 25% and Race-Based Disparities in Access to Genetic Testing

Numerous studies have underscored the diagnostic and therapeutic potential of exome or genome sequencing in critically ill pediatric populations. However, an equivalent investigation in critically ill adults remains conspicuously absent. We retrospectively analyzed whole exome sequencing (WES) data available through the PennMedicine Biobank (PMBB) from all 365 young adult patients, aged 18-40 years, with intensive care unit (ICU) admissions at the University of Pennsylvania Health System who met inclusion criteria for our study. For each participant, two Medical Genetics and Internal Medicine-trained clinicians reviewed WES reports and patient charts for variant classification, result interpretation, and identification of genetic diagnoses related to their critical illness. Of the 365 individuals in our study, 90 (24.7%) were found to have clearly diagnostic results on WES; an additional 40 (11.0%) had a suspicious variant of uncertain significance (VUS) identified; and an additional 16 (4.4%) had a medically actionable incidental finding. The diagnostic rate of exome sequencing did not decrease with increasing patient age. Affected genes were primarily involved in cardiac function (18.8%), vascular health (16.7%), cancer (16.7%), and pulmonary disease (11.5%). Only half of all diagnostic findings were known and documented in the patient chart at the time of ICU admission. Significant disparities emerged in subgroup analysis by EHR-reported race, with genetic diagnoses known/documented for 63.5% of White patients at the time of ICU admission but only for 28.6% of Black or Hispanic patients. There was a trend towards patients with undocumented genetic diagnoses having a 66% increased mortality rate, making these race-based disparities in genetic diagnosis even more concerning. Altogether, universal exome sequencing in ICU-admitted adult patients was found to yield a new definitive diagnosis in 11.2% of patients. Of these diagnoses, 76.6% conferred specific care-altering medical management recommendations. Our study suggests that the diagnostic utility of exome sequencing in critically ill young adults is similar to that observed in neonatal and pediatric populations and is age-independent. The high diagnostic rate and striking race-based disparities we find in genetic diagnoses argue for broad and universal approaches to genetic testing for critically ill adults. The widespread implementation of comprehensive genetic sequencing in the adult population promises to enhance medical care for all individuals and holds the potential to rectify disparities in genetic testing referrals, ultimately promoting more equitable healthcare delivery.

Prioritization of oligogenic variant combinations in whole exomes

MOTIVATION

Whole exome sequencing (WES) has emerged as a powerful tool for genetic research, enabling the collection of a tremendous amount of data about human genetic variation. However, properly identifying which variants are causative of a genetic disease remains an important challenge, often due to the number of variants that need to be screened. Expanding the screening to combinations of variants in two or more genes, as would be required under the oligogenic inheritance model, simply blows this problem out of proportion.

RESULTS

We present here the High-throughput oligogenic prioritizer (Hop), a novel prioritization method that uses direct oligogenic information at the variant, gene and gene pair level to detect digenic variant combinations in WES data. This method leverages information from a knowledge graph, together with specialized pathogenicity predictions in order to effectively rank variant combinations based on how likely they are to explain the patient's phenotype. The performance of Hop is evaluated in cross-validation on 36 120 synthetic exomes for training and 14 280 additional synthetic exomes for independent testing. Whereas the known pathogenic variant combinations are found in the top 20 in approximately 60% of the cross-validation exomes, 71% are found in the same ranking range when considering the independent set. These results provide a significant improvement over alternative approaches that depend simply on a monogenic assessment of pathogenicity, including early attempts for digenic ranking using monogenic pathogenicity scores.

AVAILABILITY AND IMPLEMENTATION

Hop is available at https://github.com/oligogenic/HOP.

Myocilin Mutation N480K Leads to Early Onset Juvenile and Adult-onset Primary Open Angle Glaucoma in a Six Generation Family

PRCIS

A pathogenic autosomal dominant MYOC mutation N480K detected in 6 generations of an Indian family is primarily responsible for juvenile open angle glaucoma (JOAG) and adult-onset primary open angle glaucoma (POAG), emphasizing the importance of screening this mutation at a younger age.

PURPOSE

To screen myocilin mutations in a large South Indian family with early-onset JOAG and adult-onset POAG.

METHODS

In a large South Indian family with 20 members, 8 members diagnosed as JOAG, 7 members as POAG, 4 members as JOAG suspect, and 1 member as POAG suspect were screened for myocilin ( MYOC) mutations using Sanger sequencing. Whole exome sequencing was performed on clinically suspected JOAG/POAG individuals.

RESULTS

Myocilin gene mutation N480K (c.1440C>G) was detected in 20 family members, including proband, of whom 8 were JOAG and 7 were POAG patients, 3 were JOAG suspects, and 2 were unaffected. Among the unaffected carriers, 1 was less than 5 years old, and another was 25 years old. The earliest to develop the disease was a 10-year-old child. The penetrance of the mutation was 95% over 10 years of age. This family had JOAG/POAG suspects with no N480K MYOC mutation, and they were further screened for other mutations using whole-exome sequencing. Polymorphisms CYP1B1 L432V and MYOC R76K were detected in 3 JOAG/POAG suspects, and among these 3, one had another CYP1B1 polymorphic variant R368H. The presence of the CYP1B1 polymorphism along with an MYOC polymorphic variant among the JOAG/POAG suspects needs additional studies to explore their combined role in the onset of glaucoma.

CONCLUSIONS

This study reveals that MYOC mutation is primarily responsible for JOAG and adult-onset POAG in a family, emphasizing the importance of screening for this mutation at a younger age for early treatment.

Enhancing phenotype recognition in clinical notes using large language models: PhenoBCBERT and PhenoGPT

To enhance phenotype recognition in clinical notes of genetic diseases, we developed two models-PhenoBCBERT and PhenoGPT-for expanding the vocabularies of Human Phenotype Ontology (HPO) terms. While HPO offers a standardized vocabulary for phenotypes, existing tools often fail to capture the full scope of phenotypes due to limitations from traditional heuristic or rule-based approaches. Our models leverage large language models to automate the detection of phenotype terms, including those not in the current HPO. We compare these models with PhenoTagger, another HPO recognition tool, and found that our models identify a wider range of phenotype concepts, including previously uncharacterized ones. Our models also show strong performance in case studies on biomedical literature. We evaluate the strengths and weaknesses of BERT- and GPT-based models in aspects such as architecture and accuracy. Overall, our models enhance automated phenotype detection from clinical texts, improving downstream analyses on human diseases.

Using multi-scale genomics to associate poorly annotated genes with rare diseases

BACKGROUND

Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes.

METHODS

EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool.

RESULTS

In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes.

CONCLUSIONS

We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

MSeqDR Quick-Mitome (QM): Combining Phenotype-Guided Variant Interpretation and Machine Learning Classifiers to Aid Primary Mitochondrial Disease Genetic Diagnosis

The international Mitochondrial Disease Sequence Data Resource Consortium (MSeqDR) Quick-Mitome (QM) is a web-based platform enabling automated variant interpretation of whole-exome sequencing (WES) datasets for the genetic diagnosis of primary mitochondrial diseases (PMD). Designed specifically to address the unique dual genome nature of PMD etiologies, QM includes features for both nuclear and mitochondrial DNA (mtDNA) genome analysis. QM requires VCF variant lists, HPO ID clinical phenotypes, and pedigree files for multiple-sample VCF inputs. QM maps phenotypes to HPO terms before analysis. QM analysis requires 2 to 20 min for 100,000 variants on an 8-vCPU AWS server using Exomiser's "PASS_ONLY" mode for nuclear variants. QM ranks variants based on allele frequency, phenotype-gene association, functional impact, and inheritance mode. Variants are further annotated with multiple data sources such as OMIM, ClinVar, dbNSFP, gnoMAD, MITOMAP, and MSeqDR. In addition to standard Exomiser results, QM generates an Analysis Report and QM Integrated Report with add-on mtDNA-specific analyses, including haplogroup prediction with Phy-Mer, heteroplasmy calculation, and mvTool annotations. We developed the Mitochondrial Disease Variant (MDV) classifier using XGBoost to predict variant pathogenicity for PMD. The MDV classifier was trained on >120 features and performance benchmarking showed that it correctly classified >98% of nuclear gene variants as being pathogenic or benign, and predicted PMD-causing variants with 94% precision. The MSeqDR QM server is an open-access resource for phenotype-driven dual-genome analyses for PMD diagnosis by the global mitochondrial disease community. It is publicly available for non-commercial, non-clinical research use at https://mseqdr.org/quickmitome.php. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Standardizing clinical phenotypes into human phenotype ontology (HPO) terms as the phenotype input for Quick-Mitome (QM) Basic Protocol 2: Prepare the pedigree input for multiple-sample VCF Basic Protocol 3: Quick-Mitome (QM) analysis Basic Protocol 4: Reviewing and understanding the QM Integrated Report and Analysis Report.

Prenatal Exome Sequencing Analysis in Fetuses with Various Ultrasound Findings

OBJECTIVES

To evaluate the use of Exome Sequencing (ES) for the detection of genome-wide Copy Number Variants (CNVs) and the frequency of SNVs-InDels in selected genes related to developmental disorders in a cohort of consecutive pregnancies undergoing invasive diagnostic procedures for minor or simple ultrasound findings with no indication of ES.

METHODS

Women undergoing invasive diagnostic testing (chorionic villus sampling or amniocentesis) for QF-PCR and chromosomal microarray analysis (CMA) due to prenatal ultrasound findings without an indication for ES were selected over a five-month period (May-September 2021). ES was performed to compare the efficiency of genome-wide CNV detection against CMA analysis and to detect monogenic disorders. Virtual gene panels were selected to target genes related to ultrasound findings and bioinformatic analysis was performed, prioritizing variants based on the corresponding HPO terms. The broad Fetal Gene panel for developmental disorders developed by the PAGE group was also included in the analysis.

RESULTS

A total of 59 out of 61 women consented to participate in this study. There were 36 isolated major fetal anomalies, 11 aneuploidy markers, 6 minor fetal anomalies, 4 multiple anomalies, and 2 other ultrasound signs. Following QF-PCR analysis, two uncultured samples were excluded from this study, and six (10%) common chromosome aneuploidies were detected. In the remaining 51 cases, no pathogenic CNVs were detected at CMA, nor were any pathogenic variants observed in gene panels only targeting the ultrasound indications. Two (3.9%) monogenic diseases, apparently unrelated to the fetal phenotype, were detected: blepharo-cheilo-odontic syndrome (spina bifida) and Duchenne muscular dystrophy (pyelocaliceal dilation).

CONCLUSIONS

In our series of pregnancies with ultrasound findings, common aneuploidies were the only chromosomal abnormalities present, which were detected in 10% of cases. ES CNV analysis was concordant with CMA results in all cases. No additional findings were provided by only targeting selected genes based on ultrasound findings. Broadening the analysis to a larger number of genes involved in fetal developmental disorders revealed monogenic diseases in 3.9% of cases, which, although apparently not directly related to the indications, were clinically relevant.

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND

Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome.

METHODS

We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants.

RESULTS

Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving.

CONCLUSIONS

Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Enhancing Phenotype Recognition in Clinical Notes Using Large Language Models: PhenoBCBERT and PhenoGPT

To enhance phenotype recognition in clinical notes of genetic diseases, we developed two models - PhenoBCBERT and PhenoGPT - for expanding the vocabularies of Human Phenotype Ontology (HPO) terms. While HPO offers a standardized vocabulary for phenotypes, existing tools often fail to capture the full scope of phenotypes, due to limitations from traditional heuristic or rule-based approaches. Our models leverage large language models (LLMs) to automate the detection of phenotype terms, including those not in the current HPO. We compared these models to PhenoTagger, another HPO recognition tool, and found that our models identify a wider range of phenotype concepts, including previously uncharacterized ones. Our models also showed strong performance in case studies on biomedical literature. We evaluated the strengths and weaknesses of BERT-based and GPT-based models in aspects such as architecture and accuracy. Overall, our models enhance automated phenotype detection from clinical texts, improving downstream analyses on human diseases.

Heterozygous nonsense variants in laminin subunit 3α resulting in Ebstein's anomaly

Ebstein's anomaly is a rare congenital heart disease characterized by tricuspid valve downward displacement and is associated with additional cardiac phenotypes such as left ventricle non-compaction. The genetic basis of Ebstein's anomaly has yet to be fully elucidated, although several genes (e.g., NKX2-5, MYH7, TPM1, and FLNA) may contribute to Ebstein's anomaly. Here, in two Ebstein's anomaly families (a three-generation family and a trio), we identified independent heterozygous nonsense variants in laminin subunit 3 α (LAMA3), cosegregated with phenotypes in families with reduced penetrance. Furthermore, knocking out Lama3 in mice revealed that haploinsufficiency of Lama3 led to Ebstein's malformation of the tricuspid valve and an abnormal basement membrane structure. In conclusion, we identified a novel gene-disease association of LAMA3 implicated in Ebstein's anomaly, and the findings extended our understanding of the role of the extracellular matrix in Ebstein's anomaly etiology.

diseaseGPS: auxiliary diagnostic system for genetic disorders based on genotype and phenotype

SUMMARY

The next-generation sequencing brought opportunities for the diagnosis of genetic disorders due to its high-throughput capabilities. However, the majority of existing methods were limited to only sequencing candidate variants, and the process of linking these variants to a diagnosis of genetic disorders still required medical professionals to consult databases. Therefore, we introduce diseaseGPS, an integrated platform for the diagnosis of genetic disorders that combines both phenotype and genotype data for analysis. It offers not only a user-friendly GUI web application for those without a programming background but also scripts that can be executed in batch mode for bioinformatics professionals. The genetic and phenotypic data are integrated using the ACMG-Bayes method and a novel phenotypic similarity method, to prioritize the results of genetic disorders. diseaseGPS was evaluated on 6085 cases from Deciphering Developmental Disorders project and 187 cases from Shanghai Children's hospital. The results demonstrated that diseaseGPS performed better than other commonly used methods.

AVAILABILITY AND IMPLEMENTATION

diseaseGPS is available to freely accessed at https://diseasegps.sjtu.edu.cn with source code at https://github.com/BioHuangDY/diseaseGPS.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

Applications of artificial intelligence in clinical laboratory genomics

The transition from analog to digital technologies in clinical laboratory genomics is ushering in an era of "big data" in ways that will exceed human capacity to rapidly and reproducibly analyze those data using conventional approaches. Accurately evaluating complex molecular data to facilitate timely diagnosis and management of genomic disorders will require supportive artificial intelligence methods. These are already being introduced into clinical laboratory genomics to identify variants in DNA sequencing data, predict the effects of DNA variants on protein structure and function to inform clinical interpretation of pathogenicity, link phenotype ontologies to genetic variants identified through exome or genome sequencing to help clinicians reach diagnostic answers faster, correlate genomic data with tumor staging and treatment approaches, utilize natural language processing to identify critical published medical literature during analysis of genomic data, and use interactive chatbots to identify individuals who qualify for genetic testing or to provide pre-test and post-test education. With careful and ethical development and validation of artificial intelligence for clinical laboratory genomics, these advances are expected to significantly enhance the abilities of geneticists to translate complex data into clearly synthesized information for clinicians to use in managing the care of their patients at scale.

Multi-domain knowledge graph embeddings for gene-disease association prediction

BACKGROUND

Predicting gene-disease associations typically requires exploring diverse sources of information as well as sophisticated computational approaches. Knowledge graph embeddings can help tackle these challenges by creating representations of genes and diseases based on the scientific knowledge described in ontologies, which can then be explored by machine learning algorithms. However, state-of-the-art knowledge graph embeddings are produced over a single ontology or multiple but disconnected ones, ignoring the impact that considering multiple interconnected domains can have on complex tasks such as gene-disease association prediction.

RESULTS

We propose a novel approach to predict gene-disease associations using rich semantic representations based on knowledge graph embeddings over multiple ontologies linked by logical definitions and compound ontology mappings. The experiments showed that considering richer knowledge graphs significantly improves gene-disease prediction and that different knowledge graph embeddings methods benefit more from distinct types of semantic richness.

CONCLUSIONS

This work demonstrated the potential for knowledge graph embeddings across multiple and interconnected biomedical ontologies to support gene-disease prediction. It also paved the way for considering other ontologies or tackling other tasks where multiple perspectives over the data can be beneficial. All software and data are freely available.

Open-Source Artificial Intelligence System Supports Diagnosis of Mendelian Diseases in Acutely Ill Infants

Mendelian disorders are prevalent in neonatal and pediatric intensive care units and are a leading cause of morbidity and mortality in these settings. Current diagnostic pipelines that integrate phenotypic and genotypic data are expert-dependent and time-intensive. Artificial intelligence (AI) tools may help address these challenges. Dx29 is an open-source AI tool designed for use by clinicians. It analyzes the patient's phenotype and genotype to generate a ranked differential diagnosis. We used Dx29 to retrospectively analyze 25 acutely ill infants who had been diagnosed with a Mendelian disorder, using a targeted panel of ~5000 genes. For each case, a trio (proband and both parents) file containing gene variant information was analyzed, alongside patient phenotype, which was provided to Dx29 by three approaches: (1) AI extraction from medical records, (2) AI extraction with manual review/editing, and (3) manual entry. We then identified the rank of the correct diagnosis in Dx29's differential diagnosis. With these three approaches, Dx29 ranked the correct diagnosis in the top 10 in 92-96% of cases. These results suggest that non-expert use of Dx29's automated phenotyping and subsequent data analysis may compare favorably to standard workflows utilized by bioinformatics experts to analyze genomic data and diagnose Mendelian diseases.

Krawitz PDMDPP. The future role of facial image analysis in ACMG classification guidelines

The use of next-generation sequencing (NGS) has dramatically improved the diagnosis of rare diseases. However, the analysis of genomic data has become complex with the increasing detection of variants by exome and genome sequencing. The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) developed a 5-tier classification scheme in 2015 for variant interpretation, that has since been widely adopted. Despite efforts to minimise discrepancies in the application of these criteria, inconsistencies still occur. Further specifications for individual genes were developed by Variant Curation Expert Panels (VCEPs) of the Clinical Genome Resource (ClinGen) consortium, that also take into consideration gene or disease specific features. For instance, in disorders with a highly characerstic facial gestalt a "phenotypic match" (PP4) has higher pathogenic evidence than e.g. in a non-syndromic form of intellectual disability. With computational approaches for quantifying the similarity of dysmorphic features results of such analysis can now be used in a refined Bayesian framework for the ACMG/AMP criteria.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

GA4GH Phenopackets: A Practical Introduction

The Global Alliance for Genomics and Health (GA4GH) is developing a suite of coordinated standards for genomics for healthcare. The Phenopacket is a new GA4GH standard for sharing disease and phenotype information that characterizes an individual person, linking that individual to detailed phenotypic descriptions, genetic information, diagnoses, and treatments. A detailed example is presented that illustrates how to use the schema to represent the clinical course of a patient with retinoblastoma, including demographic information, the clinical diagnosis, phenotypic features and clinical measurements, an examination of the extirpated tumor, therapies, and the results of genomic analysis. The Phenopacket Schema, together with other GA4GH data and technical standards, will enable data exchange and provide a foundation for the computational analysis of disease and phenotype information to improve our ability to diagnose and conduct research on all types of disorders, including cancer and rare diseases.

Targeted Sequencing Approach and Its Clinical Applications for the Molecular Diagnosis of Human Diseases

The outbreak of COVID-19 has positively impacted the NGS market recently. Targeted sequencing (TS) has become an important routine technique in both clinical and research settings, with advantages including high confidence and accuracy, a reasonable turnaround time, relatively low cost, and fewer data burdens with the level of bioinformatics or computational demand. Since there are no clear consensus guidelines on the wide range of next-generation sequencing (NGS) platforms and techniques, there is a vital need for researchers and clinicians to develop efficient approaches, especially for the molecular diagnosis of diseases in the emergency of the disease and the global pandemic outbreak of COVID-19. In this review, we aim to summarize different methods of TS, demonstrate parameters for TS assay designs, illustrate different TS panels, discuss their limitations, and present the challenges of TS concerning their clinical application for the molecular diagnosis of human diseases.

The frequency of somatic mutations in cancer predicts the phenotypic relevance of germline mutations

Genomic sequence mutations can be pathogenic in both germline and somatic cells. Several authors have observed that often the same genes are involved in cancer when mutated in somatic cells and in genetic diseases when mutated in the germline. Recent advances in high-throughput sequencing techniques have provided us with large databases of both types of mutations, allowing us to investigate this issue in a systematic way. Hence, we applied a machine learning based framework to this problem, comparing multiple models. The models achieved significant predictive power as shown by both cross-validation and their application to recently discovered gene/phenotype associations not used for training. We found that genes characterized by high frequency of somatic mutations in the most common cancers and ancient evolutionary age are most likely to be involved in abnormal phenotypes and diseases. These results suggest that the combination of tolerance for mutations at the cell viability level (measured by the frequency of somatic mutations in cancer) and functional relevance (demonstrated by evolutionary conservation) are the main predictors of disease genes. Our results thus confirm the deep relationship between pathogenic mutations in somatic and germline cells, provide new insight into the common origin of cancer and genetic diseases, and can be used to improve the identification of new disease genes.

Generalisable long COVID subtypes: findings from the NIH N3C and RECOVER programmes

BACKGROUND

Stratification of patients with post-acute sequelae of SARS-CoV-2 infection (PASC, or long COVID) would allow precision clinical management strategies. However, long COVID is incompletely understood and characterised by a wide range of manifestations that are difficult to analyse computationally. Additionally, the generalisability of machine learning classification of COVID-19 clinical outcomes has rarely been tested.

METHODS

We present a method for computationally modelling PASC phenotype data based on electronic healthcare records (EHRs) and for assessing pairwise phenotypic similarity between patients using semantic similarity. Our approach defines a nonlinear similarity function that maps from a feature space of phenotypic abnormalities to a matrix of pairwise patient similarity that can be clustered using unsupervised machine learning.

FINDINGS

We found six clusters of PASC patients, each with distinct profiles of phenotypic abnormalities, including clusters with distinct pulmonary, neuropsychiatric, and cardiovascular abnormalities, and a cluster associated with broad, severe manifestations and increased mortality. There was significant association of cluster membership with a range of pre-existing conditions and measures of severity during acute COVID-19. We assigned new patients from other healthcare centres to clusters by maximum semantic similarity to the original patients, and showed that the clusters were generalisable across different hospital systems. The increased mortality rate originally identified in one cluster was consistently observed in patients assigned to that cluster in other hospital systems.

INTERPRETATION

Semantic phenotypic clustering provides a foundation for assigning patients to stratified subgroups for natural history or therapy studies on PASC.

FUNDING

NIH (TR002306/OT2HL161847-01/OD011883/HG010860), U.S.D.O.E. (DE-AC02-05CH11231), Donald A. Roux Family Fund at Jackson Laboratory, Marsico Family at CU Anschutz.

PhenoExam: gene set analyses through integration of different phenotype databases

BACKGROUND

Gene set enrichment analysis (detecting phenotypic terms that emerge as significant in a set of genes) plays an important role in bioinformatics focused on diseases of genetic basis. To facilitate phenotype-oriented gene set analysis, we developed PhenoExam, a freely available R package for tool developers and a web interface for users, which performs: (1) phenotype and disease enrichment analysis on a gene set; (2) measures statistically significant phenotype similarities between gene sets and (3) detects significant differential phenotypes or disease terms across different databases.

RESULTS

PhenoExam generates sensitive and accurate phenotype enrichment analyses. It is also effective in segregating gene sets or Mendelian diseases with very similar phenotypes. We tested the tool with two similar diseases (Parkinson and dystonia), to show phenotype-level similarities but also potentially interesting differences. Moreover, we used PhenoExam to validate computationally predicted new genes potentially associated with epilepsy.

CONCLUSIONS

We developed PhenoExam, a freely available R package and Web application, which performs phenotype enrichment and disease enrichment analysis on gene set G, measures statistically significant phenotype similarities between pairs of gene sets G and G' and detects statistically significant exclusive phenotypes or disease terms, across different databases. We proved with simulations and real cases that it is useful to distinguish between gene sets or diseases with very similar phenotypes. Github R package URL is https://github.com/alexcis95/PhenoExam . Shiny App URL is https://alejandrocisterna.shinyapps.io/phenoexamweb/ .

AutoCaSc: Prioritizing candidate genes for neurodevelopmental disorders

Routine exome sequencing (ES) in individuals with neurodevelopmental disorders (NDD) remains inconclusive in >50% of the cases. Research analysis of unsolved cases can identify novel candidate genes but is time-consuming, subjective, and hard to compare between labs. The field, therefore, requires automated and standardized assessment methods to prioritize candidates for matchmaking. We developed AutoCaSc (https://autocasc.uni-leipzig.de) based on our candidate scoring scheme. We validated our approach using synthetic trios and real in-house trio ES data. AutoCaSc consistently (94.5% of all cases) scored the relevant variants in valid novel NDD genes in the top three ranks. In 93 real trio exomes, AutoCaSc identified most (97.5%) previously manually scored variants while evaluating additional high-scoring variants missed in manual evaluation. It identified candidate variants in previously undescribed NDD candidate genes (CNTN2, DLGAP1, SMURF1, NRXN3, and PRICKLE1). AutoCaSc enables anybody to quickly screen a variant for its plausibility in NDD. After contributing >40 descriptions of NDD-associated genes, we provide usage recommendations based on our extensive experience. Our implementation is capable of pipeline integration and therefore allows the screening of large cohorts for candidate genes. AutoCaSc empowers even small labs to a standardized matchmaking collaboration and to contribute to the ongoing identification of novel NDD entities.

Phenotype-aware prioritisation of rare Mendelian disease variants

A molecular diagnosis from the analysis of sequencing data in rare Mendelian diseases has a huge impact on the management of patients and their families. Numerous patient phenotype-aware variant prioritisation (VP) tools have been developed to help automate this process, and shorten the diagnostic odyssey, but performance statistics on real patient data are limited. Here we identify, assess, and compare the performance of all up-to-date, freely available, and programmatically accessible tools using a whole-exome, retinal disease dataset from 134 individuals with a molecular diagnosis. All tools were able to identify around two-thirds of the genetic diagnoses as the top-ranked candidate, with LIRICAL performing best overall. Finally, we discuss the challenges to overcome most cases remaining undiagnosed after current, state-of-the-art practices.

Whole-exome sequencing reveals PSEN1 and ATP7B combined variants as a possible cause of early-onset Lewy body dementia: a case study of genotype-phenotype correlation

Dementia with Lewy bodies is a neurodegenerative disease, sharing features with Parkinson's and Alzheimer's diseases. We report a case of a patient dementia with Lewy bodies carrying combined PSEN1 and ATP7B mutations. A man developed dementia with Lewy bodies starting at the age of 60 years. CSF biomarkers were of Alzheimer's disease and DaTSCAN was abnormal. Whole-exome sequencing revealed a heterozygous p.Ile408Thr PSEN1 variant and a homozygous p.Arg616Trp ATP7B variant. This case reinstates the need of considering ATP7B mutations when evaluating a patient with parkinsonism and supports p.Ile408Thr as a pathogenic PSEN1 variant.

A catalog of the genetic causes of hereditary angioedema in the Canary Islands (Spain)

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

An additional whole-exome sequencing study in 102 panel-undiagnosed patients: A retrospective study in a Chinese craniosynostosis cohort

Craniosynostosis (CRS) is a disease with prematurely fused cranial sutures. In the last decade, the whole-exome sequencing (WES) was widely used in Caucasian populations. The WES largely contributed in genetic diagnosis and exploration on new genetic mechanisms of CRS. In this study, we enrolled 264 CRS patients in China. After a 17-gene-panel sequencing designed in the previous study, 139 patients were identified with pathogenic/likely pathogenic (P/LP) variants according to the ACMG guideline as positive genetic diagnosis. WES was then performed on 102 patients with negative genetic diagnosis by panel. Ten P/LP variants were additionally identified in ten patients, increasing the genetic diagnostic yield by 3.8% (10/264). The novel variants in ANKH, H1-4, EIF5A, SOX6, and ARID1B expanded the mutation spectra of CRS. Then we designed a compatible research pipeline (RP) for further exploration. The RP could detect all seven P/LP SNVs and InDels identified above, in addition to 15 candidate variants found in 13 patients with worthy of further study. In sum, the 17-gene panel and WES identified positive genetic diagnosis for 56.4% patients (149/264) in 16 genes. At last, in our estimation, the genetic testing strategy of “Panel-first” saves 24.3% of the cost compared with “WES only”, suggesting the “Panel-first” is an economical strategy.

PheNominal: an EHR-integrated web application for structured deep phenotyping at the point of care

BACKGROUND

Clinical phenotype information greatly facilitates genetic diagnostic interpretations pipelines in disease. While post-hoc extraction using natural language processing on unstructured clinical notes continues to improve, there is a need to improve point-of-care collection of patient phenotypes. Therefore, we developed "PheNominal", a point-of-care web application, embedded within Epic electronic health record (EHR) workflows, to permit capture of standardized phenotype data.

METHODS

Using bi-directional web services available within commercial EHRs, we developed a lightweight web application that allows users to rapidly browse and identify relevant terms from the Human Phenotype Ontology (HPO). Selected terms are saved discretely within the patient's EHR, permitting reuse both in clinical notes as well as in downstream diagnostic and research pipelines.

RESULTS

In the 16 months since implementation, PheNominal was used to capture discrete phenotype data for over 1500 individuals and 11,000 HPO terms during clinic and inpatient encounters for a genetic diagnostic consultation service within a quaternary-care pediatric academic medical center. An average of 7 HPO terms were captured per patient. Compared to a manual workflow, the average time to enter terms for a patient was reduced from 15 to 5 min per patient, and there were fewer annotation errors.

CONCLUSIONS

Modern EHRs support integration of external applications using application programming interfaces. We describe a practical application of these interfaces to facilitate deep phenotype capture in a discrete, structured format within a busy clinical workflow. Future versions will include a vendor-agnostic implementation using FHIR. We describe pilot efforts to integrate structured phenotyping through controlled dictionaries into diagnostic and research pipelines, reducing manual effort for phenotype documentation and reducing errors in data entry.

Generalizable Long COVID Subtypes: Findings from the NIH N3C and RECOVER Programs

Accurate stratification of patients with post-acute sequelae of SARS-CoV-2 infection (PASC, or long COVID) would allow precision clinical management strategies. However, the natural history of long COVID is incompletely understood and characterized by an extremely wide range of manifestations that are difficult to analyze computationally. In addition, the generalizability of machine learning classification of COVID-19 clinical outcomes has rarely been tested. We present a method for computationally modeling PASC phenotype data based on electronic healthcare records (EHRs) and for assessing pairwise phenotypic similarity between patients using semantic similarity. Our approach defines a nonlinear similarity function that maps from a feature space of phenotypic abnormalities to a matrix of pairwise patient similarity that can be clustered using unsupervised machine learning procedures. Using k-means clustering of this similarity matrix, we found six distinct clusters of PASC patients, each with distinct profiles of phenotypic abnormalities. There was a significant association of cluster membership with a range of pre-existing conditions and with measures of severity during acute COVID-19. Two of the clusters were associated with severe manifestations and displayed increased mortality. We assigned new patients from other healthcare centers to one of the six clusters on the basis of maximum semantic similarity to the original patients. We show that the identified clusters were generalizable across different hospital systems and that the increased mortality rate was consistently observed in two of the clusters. Semantic phenotypic clustering can provide a foundation for assigning patients to stratified subgroups for natural history or therapy studies on PASC.

Germline mosaicism of a missense variant in KCNC2 in a multiplex family with autism and epilepsy characterized by long-read sequencing

Currently, protein-coding de novo variants and large copy number variants have been identified as important for ~30% of individuals with autism. One approach to identify relevant variation in individuals who lack these types of events is by utilizing newer genomic technologies. In this study, highly accurate PacBio HiFi long-read sequencing was applied to a family with autism, epileptic encephalopathy, cognitive impairment, and mild dysmorphic features (two affected female siblings, unaffected parents, and one unaffected male sibling) with no known clinical variant. From our long-read sequencing data, a de novo missense variant in the KCNC2 gene (encodes Kv3.2) was identified in both affected children. This variant was phased to the paternal chromosome of origin and is likely a germline mosaic. In silico assessment revealed the variant was not in controls, highly conserved, and predicted damaging. This specific missense variant (Val473Ala) has been shown in both an ortholog and paralog of Kv3.2 to accelerate current decay, shift the voltage dependence of activation, and prevent the channel from entering a long-lasting open state. Seven additional missense variants have been identified in other individuals with neurodevelopmental disorders (p = 1.03 × 10^-5 ). KCNC2 is most highly expressed in the brain; in particular, in the thalamus and is enriched in GABAergic neurons. Long-read sequencing was useful in discovering the relevant variant in this family with autism that had remained a mystery for several years and will potentially have great benefits in the clinic once it is widely available.

New Developments and Possibilities in Reanalysis and Reinterpretation of Whole Exome Sequencing Datasets for Unsolved Rare Diseases Using Machine Learning Approaches

Rare diseases impact the lives of 300 million people in the world. Rapid advances in bioinformatics and genomic technologies have enabled the discovery of causes of 20–30% of rare diseases. However, most rare diseases have remained as unsolved enigmas to date. Newer tools and availability of high throughput sequencing data have enabled the reanalysis of previously undiagnosed patients. In this review, we have systematically compiled the latest developments in the discovery of the genetic causes of rare diseases using machine learning methods. Importantly, we have detailed methods available to reanalyze existing whole exome sequencing data of unsolved rare diseases. We have identified different reanalysis methodologies to solve problems associated with sequence alterations/mutations, variation re-annotation, protein stability, splice isoform malfunctions and oligogenic analysis. In addition, we give an overview of new developments in the field of rare disease research using whole genome sequencing data and other omics.

Prenatal phenotyping: A community effort to enhance the Human Phenotype Ontology

Technological advances in both genome sequencing and prenatal imaging are increasing our ability to accurately recognize and diagnose Mendelian conditions prenatally. Phenotype-driven early genetic diagnosis of fetal genetic disease can help to strategize treatment options and clinical preventive measures during the perinatal period, to plan in utero therapies, and to inform parental decision-making. Fetal phenotypes of genetic diseases are often unique and at present are not well understood; more comprehensive knowledge about prenatal phenotypes and computational resources have an enormous potential to improve diagnostics and translational research. The Human Phenotype Ontology (HPO) has been widely used to support diagnostics and translational research in human genetics. To better support prenatal usage, the HPO consortium conducted a series of workshops with a group of domain experts in a variety of medical specialties, diagnostic techniques, as well as diseases and phenotypes related to prenatal medicine, including perinatal pathology, musculoskeletal anomalies, neurology, medical genetics, hydrops fetalis, craniofacial malformations, cardiology, neonatal-perinatal medicine, fetal medicine, placental pathology, prenatal imaging, and bioinformatics. We expanded the representation of prenatal phenotypes in HPO by adding 95 new phenotype terms under the Abnormality of prenatal development or birth (HP:0001197) grouping term, and revised definitions, synonyms, and disease annotations for most of the 152 terms that existed before the beginning of this effort. The expansion of prenatal phenotypes in HPO will support phenotype-driven prenatal exome and genome sequencing for precision genetic diagnostics of rare diseases to support prenatal care.

seqr: A web-based analysis and collaboration tool for rare disease genomics

Exome and genome sequencing have become the tools of choice for rare disease diagnosis, leading to large amounts of data available for analyses. To identify causal variants in these datasets, powerful filtering and decision support tools that can be efficiently used by clinicians and researchers are required. To address this need, we developed seqr - an open-source, web-based tool for family-based monogenic disease analysis that allows researchers to work collaboratively to search and annotate genomic callsets. To date, seqr is being used in several research pipelines and one clinical diagnostic lab. In our own experience through the Broad Institute Center for Mendelian Genomics, seqr has enabled analyses of over 10,000 families, supporting the diagnosis of more than 3,800 individuals with rare disease and discovery of over 300 novel disease genes. Here, we describe a framework for genomic analysis in rare disease that leverages seqr's capabilities for variant filtration, annotation, and causal variant identification, as well as support for research collaboration and data sharing. The seqr platform is available as open source software, allowing low-cost participation in rare disease research, and a community effort to support diagnosis and gene discovery in rare disease.

Innovative in Silico Approaches for Characterization of Genes and Proteins

Bioinformatics is an amalgamation of biology, mathematics and computer science. It is a science which gathers the information from biology in terms of molecules and applies the informatic techniques to the gathered information for understanding and organizing the data in a useful manner. With the help of bioinformatics, the experimental data generated is stored in several databases available online like nucleotide database, protein databases, GENBANK and others. The data stored in these databases is used as reference for experimental evaluation and validation. Till now several online tools have been developed to analyze the genomic, transcriptomic, proteomics, epigenomics and metabolomics data. Some of them include Human Splicing Finder (HSF), Exonic Splicing Enhancer Mutation taster, and others. A number of SNPs are observed in the non-coding, intronic regions and play a role in the regulation of genes, which may or may not directly impose an effect on the protein expression. Many mutations are thought to influence the splicing mechanism by affecting the existing splice sites or creating a new sites. To predict the effect of mutation (SNP) on splicing mechanism/signal, HSF was developed. Thus, the tool is helpful in predicting the effect of mutations on splicing signals and can provide data even for better understanding of the intronic mutations that can be further validated experimentally. Additionally, rapid advancement in proteomics have steered researchers to organize the study of protein structure, function, relationships, and dynamics in space and time. Thus the effective integration of all of these technological interventions will eventually lead to steering up of next-generation systems biology, which will provide valuable biological insights in the field of research, diagnostic, therapeutic and development of personalized medicine.

Deep phenotyping: symptom annotation made simple with SAMS

Precision medicine needs precise phenotypes. The Human Phenotype Ontology (HPO) uses clinical signs instead of diagnoses and has become the standard annotation for patients’ phenotypes when describing single gene disorders. Use of the HPO beyond human genetics is however still limited. With SAMS (Symptom Annotation Made Simple), we want to bring sign-based phenotyping to routine clinical care, to hospital patients as well as to outpatients. Our web-based application provides access to three widely used annotation systems: HPO, OMIM, Orphanet. Whilst data can be stored in our database, phenotypes can also be imported and exported as Global Alliance for Genomics and Health (GA4GH) Phenopackets without using the database. The web interface can easily be integrated into local databases, e.g. clinical information systems. SAMS offers users to share their data with others, empowering patients to record their own signs and symptoms (or those of their children) and thus provide their doctors with additional information. We think that our approach will lead to better characterised patients which is not only helpful for finding disease mutations but also to better understand the pathophysiology of diseases and to recruit patients for studies and clinical trials. SAMS is freely available at https://www.genecascade.org/SAMS/.

Deep phenotyping and whole-exome sequencing improved the diagnostic yield for nuclear pedigrees with neurodevelopmental disorders

Background

Neurodevelopmental disorders, a group of early‐onset neurological disorders with significant clinical and genetic heterogeneity, remain a diagnostic challenge for clinical genetic evaluation. Therefore, we assessed the diagnostic yield by combining standard phenotypes and whole‐exome sequencing in families with these disorders that were “not yet diagnosed” by the traditional testing methods.

Methods

Using a standardized vocabulary of phenotypic abnormalities from human phenotype ontology (HPO), we performed deep phenotyping for 45 “not yet diagnosed” pedigrees to characterize multiple clinical features extracted from Chinese electronic medical records (EMRs). By matching HPO terms with known human diseases and phenotypes from model organisms, together with whole‐exome sequencing data, we prioritized candidate mutations/genes. We made probable genetic diagnoses for the families.

Results

We obtained a diagnostic yield of 29% (13 out of 45) with probably genetic diagnosis, of which compound heterozygosity and de novo mutations accounted for 77% (10/13) of the diagnosis. Of note, these pedigrees are accompanied by a more significant number of non‐neurological features.

Conclusions

Deep phenotyping and whole‐exome sequencing improve the etiological evaluation for neurodevelopmental disorders in the clinical setting.

AutozygosityMapper: Identification of disease-mutations in consanguineous families

With the shift from SNP arrays to high-throughput sequencing, most researchers studying diseases in consanguineous families do not rely on linkage analysis any longer, but simply search for deleterious variants which are homozygous in all patients. AutozygosityMapper allows the fast and convenient identification of disease mutations in patients from consanguineous pedigrees by focussing on homozygous segments shared by all patients. Users can upload multi-sample VCF files, including WGS data, without any pre-processing. Genome-wide runs of homozygosity and the underlying genotypes are presented in graphical interfaces. AutozygosityMapper extends the functions of its predecessor, HomozygosityMapper, to the search for autozygous regions, in which all patients share the same homozygous genotype. We provide export of VCF files containing only the variants found in homozygous regions, this usually reduces the number of variants by two orders of magnitude. These regions can also directly be analysed with our disease mutation identification tool MutationDistiller. The application comes with simple and intuitive graphical interfaces for data upload, analysis, and results. We kept the structure of HomozygosityMapper so that previous users will find it easy to switch. With AutozygosityMapper, we provide a fast web-based way to identify disease mutations in consanguineous families. AutozygosityMapper is freely available at https://www.genecascade.org/AutozygosityMapper/.

SvAnna: efficient and accurate pathogenicity prediction of coding and regulatory structural variants in long-read genome sequencing

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

Implementation of Exome Sequencing in Prenatal Diagnosis and Impact on Genetic Counseling: The Polish Experience

Background: Despite advances in routine prenatal cytogenetic testing, most anomalous fetuses remain without a genetic diagnosis. Exome sequencing (ES) is a molecular technique that identifies sequence variants across protein-coding regions and is now increasingly used in clinical practice. Fetal phenotypes differ from postnatal and, therefore, prenatal ES interpretation requires a large amount of data deriving from prenatal testing. The aim of our study was to present initial results of the implementation of ES to prenatal diagnosis in Polish patients and to discuss its possible clinical impact on genetic counseling. Methods: In this study we performed a retrospective review of all fetal samples referred to our laboratory for ES from cooperating centers between January 2017 and June 2021. Results: During the study period 122 fetuses were subjected to ES at our institution. There were 52 abnormal ES results: 31 in the group of fetuses with a single organ system anomaly and 21 in the group of fetuses with multisystem anomalies. The difference between groups was not statistically significant. There were 57 different pathogenic or likely pathogenic variants reported in 33 different genes. The most common were missense variants. In 17 cases the molecular diagnosis had an actual clinical impact on subsequent pregnancies or other family members. Conclusions: Exome sequencing increases the detection rate in fetuses with structural anomalies and improves genetic counseling for both the affected couple and their relatives.