Generating Clinical-Grade Gene-Disease Validity Classifications Through the ClinGen Data Platforms.

Clinical genetic laboratories must have access to clinically validated biomedical data for precision medicine. A lack of accessibility, normalized structure, and consistency in evaluation complicates interpretation of disease causality, resulting in confusion in assessing the clinical validity of genes and genetic variants for diagnosis. A key goal of the Clinical Genome Resource (ClinGen) is to fill the knowledge gap concerning the strength of evidence supporting the role of a gene in a monogenic disease, which is achieved through a process known as Gene-Disease Validity curation. Here we review the work of ClinGen in developing a curation infrastructure that supports the standardization, harmonization, and dissemination of Gene-Disease Validity data through the creation of frameworks and the utilization of common data standards. This infrastructure is based on several applications, including the ClinGen GeneTracker, Gene Curation Interface, Data Exchange, GeneGraph, and website.

The landscape of regional missense mutational intolerance quantified from 125,748 exomes.

Missense variants can have a range of functional impacts depending on factors such as the specific amino acid substitution and location within the gene. To interpret their deleteriousness, studies have sought to identify regions within genes that are specifically intolerant of missense variation 1-12 . Here, we leverage the patterns of rare missense variation in 125,748 individuals in the Genome Aggregation Database (gnomAD) 13 against a null mutational model to identify transcripts that display regional differences in missense constraint. Missense-depleted regions are enriched for ClinVar 14 pathogenic variants, de novo missense variants from individuals with neurodevelopmental disorders (NDDs) 15,16 , and complex trait heritability. Following ClinGen calibration recommendations for the ACMG/AMP guidelines, we establish that regions with less than 20% of their expected missense variation achieve moderate support for pathogenicity. We create a missense deleteriousness metric (MPC) that incorporates regional constraint and outperforms other deleteriousness scores at stratifying case and control de novo missense variation, with a strong enrichment in NDDs. These results provide additional tools to aid in missense variant interpretation.

Genetic sex validation for sample tracking in next-generation sequencing clinical testing.

OBJECTIVE: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. RESULTS: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors (49.09%), samples from transgender participants (3.64%) and stem cell or bone marrow transplant patients (7.27%) along with undetermined sample mix-ups (40%) for which sample swaps occurred prior to arrival at genome centers, however the exact cause of the events at the sampling sites resulting in the mix-ups were not able to be determined.

Heterozygous MAP3K20 variants cause ectodermal dysplasia, craniosynostosis, sensorineural hearing loss, and limb anomalies.

Biallelic pathogenic variants in MAP3K20, which encodes a mitogen-activated protein kinase, are a rare cause of split-hand foot malformation (SHFM), hearing loss, and nail abnormalities or congenital myopathy. However, heterozygous variants in this gene have not been definitively associated with a phenotype. Here, we describe the phenotypic spectrum associated with heterozygous de novo variants in the linker region between the kinase domain and leucine zipper domain of MAP3K20. We report five individuals with diverse clinical features, including craniosynostosis, limb anomalies, sensorineural hearing loss, and ectodermal dysplasia-like phenotypes who have heterozygous de novo variants in this specific region of the gene. These individuals exhibit both shared and unique clinical manifestations, highlighting the complexity and variability of the disorder. We propose that the involvement of MAP3K20 in endothelial-mesenchymal transition provides a plausible etiology of these features. Together, these findings characterize a disorder that both expands the phenotypic spectrum associated with MAP3K20 and highlights the need for further studies on its role in early human development.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing.

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing (ES), the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare disease. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery which should in turn increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints, and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks like Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Detecting missed diagnoses of spinal muscular atrophy in genome, exome, and panel sequencing datasets.

Spinal muscular atrophy (SMA) is a genetic disorder that causes progressive degeneration of lower motor neurons and the subsequent loss of muscle function throughout the body. It is the second most common recessive disorder in individuals of European descent and is present in all populations. Accurate tools exist for diagnosing SMA from short read and long read genome sequencing data. However, there are no publicly available tools for GRCh38-aligned data from panel or exome sequencing assays which continue to be used as first line tests for neuromuscular disorders. We therefore developed and extensively validated a new tool - SMA Finder - that can diagnose SMA not only in genome, but also exome and targeted sequencing samples aligned to GRCh37, GRCh38, or T2T-CHM13. It works by evaluating aligned reads that overlap the c.840 position of SMN1 and SMN2 in order to detect the most common molecular causes of SMA. We applied SMA Finder to 16,626 exomes and 3,911 genomes from heterogeneous rare disease cohorts sequenced at the Broad Institute Center for Mendelian Genomics as well as 1,157 exomes and 8,762 targeted sequencing samples from Tartu University Hospital. SMA Finder correctly identified all 16 known SMA cases and reported nine novel diagnoses which have since been confirmed by clinical testing, with another four novel diagnoses undergoing validation. Notably, out of the 29 total SMA positive cases, 21 had an initial clinical diagnosis of muscular dystrophy, congenital myasthenic syndrome, or congenital myopathy. This underscored the frequency with which SMA can be misdiagnosed as other neuromuscular disorders and confirmed the utility of using SMA Finder to reanalyze phenotypically diverse neuromuscular disease cohorts. Finally, we evaluated SMA Finder on 198,868 individuals that had both exome and genome sequencing data within the UK Biobank (UKBB) and found that SMA Finder's overall false positive rate was less than 1 / 200,000 exome samples, and its positive predictive value (PPV) was 96%. We also observed 100% concordance between UKBB exome and genome calls. This analysis showed that, even though it is located within a segmental duplication, the most common causal variant for SMA can be detected with comparable accuracy to monogenic disease variants in non-repetitive regions. Additionally, the high PPV demonstrated by SMA Finder, the existence of treatment options for SMA in which early diagnosis is imperative for therapeutic benefit, as well as widespread availability of clinical confirmatory testing for SMA, may warrant the addition of SMN1 to the ACMG list of genes with reportable secondary findings after genome and exome sequencing.

A gene pathogenicity tool "GenePy" identifies missed biallelic diagnoses in the 100,000 Genomes Project.

PURPOSE: The 100,000 Genomes Project diagnosed a quarter of affected participants, but 26% of diagnoses were not on the applied gene panel(s); with many being de novo variants. Assessing biallelic variants without a gene panel is more challenging. METHODS: We sought to identify missed biallelic diagnoses using GenePy, which incorporates allele frequency, zygosity, and a user-defined deleterious metric, generating an aggregate GenePy score per gene, per participant. We calculated GenePy scores for 2862 recessive disease genes in 78,216 100,000 Genomes Project participants. For each gene, we ranked participant GenePy scores and scrutinized affected participants without a diagnosis, whose scores ranked among the top 5 for each gene. In cases which participant phenotypes overlapped with the disease gene of interest, we extracted rare variants and applied phase, ClinVar, and ACMG classification. RESULTS: 3184 affected individuals without a molecular diagnosis had a top-5-ranked GenePy score and 682 of 3184 (21%) had phenotypes overlapping with a top-ranking gene. In 122 of 669 (18%) phenotype-matched cases (excluding 13 withdrawn participants), we identified a putative missed diagnosis (2.2% of all undiagnosed participants). A further 334 of 669 (50%) cases have a possible missed diagnosis but require functional validation. CONCLUSION: Applying GenePy at scale has identified 456 potential diagnoses, demonstrating the value of novel diagnostic strategies.

Rare ACTN2 Frameshift Variants Resulting in Protein Extension Cause Distal Myopathy and Hypertrophic Cardiomyopathy through Protein Aggregation.

Distal myopathies are a group of rare, inherited muscular disorders characterized by progressive loss of muscle fibers that begins in the distal parts of arms and legs. Recently, variants in a new disease gene, ACTN2 , have been shown to cause distal myopathy. ACTN2 , a gene previously only associated with cardiomyopathies, encodes alpha-actinin-2, a protein expressed in both cardiac and skeletal sarcomeres. The primary function of alpha-actinin-2 is to link actin and titin to the sarcomere Z-disk. New ACTN2 variants are continuously discovered, however, the clinical significance of many variants remains unknown. Thus, lack of clear genotype-phenotype correlations in ACTN2 -related diseases, actininopathies, persists. Objective: The objective of the study is to characterize the pathomechanisms underlying actininopathies. Methods: Functional characterization in C2C12 cell models of several ACTN2 variants is conducted, including frameshift and missense variants associated with dominant actininopathies. We assess the genotype-phenotype correlations of actininopathies using clinical data from several patients carrying these variants. Results: The results show that the missense variants associated with a recessive form of actininopathy do not cause detectable alpha-actinin-2 aggregates in the cell model. Conversely, dominant frameshift variants causing a protein extension do produce alpha-actinin-2 aggregates. Interpretation: The results suggest that alpha-actinin-2 aggregation is the disease mechanism underlying some dominant actininopathies, and thus we recommend that protein-extending frameshift variants in ACTN2 should be classified as pathogenic. However, this mechanism is likely elicited by only a limited number of variants. Alternative functional characterization methods should be explored to further investigate other molecular mechanisms underlying actininopathies.

A harmonized public resource of deeply sequenced diverse human genomes.

Underrepresented populations are often excluded from genomic studies due in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high quality set of 4,094 whole genomes from HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also demonstrate substantial added value from this dataset compared to the prior versions of the component resources, typically combined via liftover and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared to previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Will variants of uncertain significance still exist in 2030?

In 2020, the National Human Genome Research Institute (NHGRI) made ten "bold predictions," including that "the clinical relevance of all encountered genomic variants will be readily predictable, rendering the diagnostic designation 'variant of uncertain significance (VUS)' obsolete." We discuss the prospects for this prediction, arguing that many, if not most, VUS in coding regions will be resolved by 2030. We outline a confluence of recent changes making this possible, especially advances in the standards for variant classification that better leverage diverse types of evidence, improvements in computational variant effect predictor performance, scalable multiplexed assays of variant effect capable of saturating the genome, and data-sharing efforts that will maximize the information gained from each new individual sequenced and variant interpreted. We suggest that clinicians and researchers can realize a future where VUSs have largely been eliminated, in line with the NHGRI's bold prediction. The length of time taken to reach this future, and thus whether we are able to achieve the goal of largely eliminating VUSs by 2030, is largely a consequence of the choices made now and in the next few years. We believe that investing in eliminating VUSs is worthwhile, since their predominance remains one of the biggest challenges to precision genomic medicine.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both copies of a gene are impacted by a damaging genetic variant. When a patient carries two potentially causal variants in a gene, accurate diagnosis requires determining that these variants occur on different copies of the chromosome (that is, are in trans) rather than on the same copy (that is, in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. Here we developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in the Genome Aggregation Database (v2, n = 125,748 exomes). Our approach estimates phase with 96% accuracy, both in trio data and in patients with Mendelian conditions and presumed causal compound heterozygous variants. We provide a public resource of phasing estimates for coding variants and counts per gene of rare variants in trans that can aid interpretation of rare co-occurring variants in the context of recessive disease.

A genomic mutational constraint map using variation in 76,156 human genomes.

The depletion of disruptive variation caused by purifying natural selection (constraint) has been widely used to investigate protein-coding genes underlying human disorders1-4, but attempts to assess constraint for non-protein-coding regions have proved more difficult. Here we aggregate, process and release a dataset of 76,156 human genomes from the Genome Aggregation Database (gnomAD)-the largest public open-access human genome allele frequency reference dataset-and use it to build a genomic constraint map for the whole genome (genomic non-coding constraint of haploinsufficient variation (Gnocchi)). We present a refined mutational model that incorporates local sequence context and regional genomic features to detect depletions of variation. As expected, the average constraint for protein-coding sequences is stronger than that for non-coding regions. Within the non-coding genome, constrained regions are enriched for known regulatory elements and variants that are implicated in complex human diseases and traits, facilitating the triangulation of biological annotation, disease association and natural selection to non-coding DNA analysis. More constrained regulatory elements tend to regulate more constrained protein-coding genes, which in turn suggests that non-coding constraint can aid the identification of constrained genes that are as yet unrecognized by current gene constraint metrics. We demonstrate that this genome-wide constraint map improves the identification and interpretation of functional human genetic variation.

Recommendations for risk allele evidence curation, classification, and reporting from the ClinGen Low Penetrance/Risk Allele Working Group.

PURPOSE: Genetic variants at the low end of the penetrance spectrum have historically been challenging to interpret because their high population frequencies exceed the disease prevalence of the associated condition, leading to a lack of clear segregation between the variant and disease. There is currently substantial variation in the classification of these variants, and no formal classification framework has been widely adopted. The Clinical Genome Resource Low Penetrance/Risk Allele Working Group was formed to address these challenges and promote harmonization within the clinical community. METHODS: The work presented here is the product of internal and community Likert-scaled surveys in combination with expert consensus within the Working Group. RESULTS: We formally recognize risk alleles and low-penetrance variants as distinct variant classes from those causing highly penetrant disease that require special considerations regarding their clinical classification and reporting. First, we provide a preferred terminology for these variants. Second, we focus on risk alleles and detail considerations for reviewing relevant studies and present a framework for the classification these variants. Finally, we discuss considerations for clinical reporting of risk alleles. CONCLUSION: These recommendations support harmonized interpretation, classification, and reporting of variants at the low end of the penetrance spectrum.

Toward robust clinical genome interpretation: Developing a consistent terminology to characterize Mendelian disease-gene relationships-allelic requirement, inheritance modes, and disease mechanisms.

PURPOSE: The terminology used for gene-disease curation and variant annotation to describe inheritance, allelic requirement, and both sequence and functional consequences of a variant is currently not standardized. There is considerable discrepancy in the literature and across clinical variant reporting in the derivation and application of terms. Here, we standardize the terminology for the characterization of disease-gene relationships to facilitate harmonized global curation and to support variant classification within the ACMG/AMP framework. METHODS: Terminology for inheritance, allelic requirement, and both structural and functional consequences of a variant used by Gene Curation Coalition members and partner organizations was collated and reviewed. Harmonized terminology with definitions and use examples was created, reviewed, and validated. RESULTS: We present a standardized terminology to describe gene-disease relationships, and to support variant annotation. We demonstrate application of the terminology for classification of variation in the ACMG SF 2.0 genes recommended for reporting of secondary findings. Consensus terms were agreed and formalized in both Sequence Ontology (SO) and Human Phenotype Ontology (HPO) ontologies. Gene Curation Coalition member groups intend to use or map to these terms in their respective resources. CONCLUSION: The terminology standardization presented here will improve harmonization, facilitate the pooling of curation datasets across international curation efforts and, in turn, improve consistency in variant classification and genetic test interpretation.

ClinGen guidance for use of the PP1/BS4 co-segregation and PP4 phenotype specificity criteria for sequence variant pathogenicity classification.

The 2015 American College of Medical Genetics and Genomics and the Association for Molecular Pathology variant classification publication established a standard employed internationally to guide laboratories in variant assessment. Those recommendations included both pathogenic (PP1) and benign (BS4) criteria for evaluating the inheritance patterns of variants, but details of how to apply those criteria at appropriate evidence levels were sparse. Several publications have since attempted to provide additional guidance, but anecdotally, this issue is still challenging. Additionally, it is not clear that those prior efforts fully distinguished disease-gene identification considerations from variant pathogenicity considerations nor did they address autosomal-recessive and X-linked inheritance. Here, we have taken a mixed inductive and deductive approach to this problem using real diseases as examples. We have developed a practical heuristic for genetic co-segregation evidence and have also determined that the specific phenotype criterion (PP4) is inseparably coupled to the co-segregation criterion. We have also determined that negative evidence at one locus constitutes positive evidence for other loci for disorders with locus heterogeneity. Finally, we provide a points-based system for evaluating phenotype and co-segregation as evidence types to support or refute a locus and show how that can be integrated into the Bayesian framework now used for variant classification and consistent with the 2015 guidelines.

A Panel-Agnostic Strategy 'HiPPo' Improves Diagnostic Efficiency in the UK Genomic Medicine Service.

Genome sequencing is available as a clinical test in the UK through the Genomic Medicine Service (GMS). The GMS analytical strategy predominantly filters genome data on preselected gene panels. Whilst this reduces variants requiring assessment by reporting laboratories, pathogenic variants outside applied panels may be missed, and variants in genes without established disease-gene relationships are largely ignored. This study compares the analysis of a research exome to a GMS clinical genome for the same patients. For the research exome, we applied a panel-agnostic approach filtering for variants with High Pathogenic Potential (HiPPo) using ClinVar, allele frequency, and in silico prediction tools. We then restricted HiPPo variants to Gene Curation Coalition (GenCC) disease genes. These results were compared with the GMS genome panel-based approach. Twenty-four participants from eight families underwent parallel research exome and GMS genome sequencing. Exome HiPPo analysis identified a similar number of variants as the GMS panel-based approach. GMS genome analysis returned two pathogenic variants and one de novo variant. Exome HiPPo analysis returned the same variants plus an additional pathogenic variant and three further de novo variants in novel genes, where case series are underway. When HiPPo was restricted to GenCC disease genes, statistically fewer variants required assessment to identify more pathogenic variants than reported by the GMS, giving a diagnostic rate per variant assessed of 20% for HiPPo versus 3% for the GMS. With UK plans to sequence 5 million genomes, strategies are needed to optimise genome analysis beyond gene panels whilst minimising the burden of variants requiring clinical assessment.

Genetic Sex Validation for Sample Tracking in Clinical Testing.

Objective: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. Results: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors, samples from transgender participants and stem cell or bone marrow transplant patients along with undetermined sample mix-ups.

Aligning NIH's existing data use restrictions to the GA4GH DUO standard.

It is widely accepted that large-scale genomic data (e.g., whole-genome sequencing, whole-exome sequencing, and genome-wide association study data) be shared through a controlled-access mechanism. This protects the privacy of research participants and ensures downstream uses of data align with participants' informed consent regarding future sharing of their data. In 2019, GA4GH approved the Data Use Ontology (DUO) standard to define data use terms with machine-readable representations to represent how a dataset can be used. We endeavored to determine the parity of existing data use restrictions ("Data Use Limitations" [DULs]) for datasets registered in the National Institutes of Health database for Genotypes and Phenotypes (dbGaP) with the DUO standard. We found substantial (93%) parity between the dbGaP DULs (n = 3,575) and DUO. This study demonstrates the comprehensiveness of the DUO standard and encourages data stewards to standardize data use restrictions in machine-readable formats to facilitate data sharing.