Recommendations for application of the functional evidence PS3/BS3 criterion using the ACMG/AMP sequence variant interpretation framework

Assigning credit where it is due: an information content score to capture the clinical value of multiplexed assays of variant effect

BACKGROUND

A variant can be pathogenic or benign with relation to a human disease. Current classification categories from benign to pathogenic reflect a probabilistic summary of the current understanding. A primary metric of clinical utility for multiplexed assays of variant effect (MAVE) is the number of variants that can be reclassified from uncertain significance (VUS). However, a gap in this measure of utility is that it underrepresents the information gained from MAVEs. The aim of this study was to develop an improved quantification metric for MAVE utility. We propose adopting an information content approach that includes data that does not reclassify variants will better reflect true information gain. We adopted an information content approach to evaluate the information gain, in bits, for MAVEs of BRCA1, PTEN, and TP53. Here, one bit represents the amount of information required to completely classify a single variant starting from no information.

RESULTS

BRCA1 MAVEs produced a total of 831.2 bits of information, 6.58% of the total missense information in BRCA1 and a 22-fold increase over the information that only contributed to VUS reclassification. PTEN MAVEs produced 2059.6 bits of information which represents 32.8% of the total missense information in PTEN and an 85-fold increase over the information that contributed to VUS reclassification. TP53 MAVEs produced 277.8 bits of information which represents 6.22% of the total missense information in TP53 and a 3.5-fold increase over the information that contributed to VUS reclassification.

CONCLUSIONS

An information content approach will more accurately portray information gained through MAVE mapping efforts than by counting the number of variants reclassified. This information content approach may also help define the impact of guideline changes that modify the information definitions used to classify groups of variants.

Calibration of variant effect predictors on genome-wide data masks heterogeneous performance across genes

In silico variant effect predictions are available for nearly all missense variants but played a minimal role in clinical variant classification because they were deemed to provide only supporting evidence. Recently, the ClinGen Sequence Variant Interpretation (SVI) Working Group updated recommendations for variant effect prediction use. By analyzing control pathogenic and benign variants across all genes, they were able to compute evidence strength for predictor score intervals with some intervals generating moderate, strong, or even very strong evidence. However, this genome-wide approach could obscure heterogeneous predictor performance in different genes. We quantified the gene-by-gene performance of two top predictors, REVEL and BayesDel, by analyzing control variants in each predictor score interval in 3,668 disease-relevant genes. Approximately 10% of intervals had sufficient control variants for analysis, and ∼70% of these intervals exceeded the maximum number of incorrect predictions implied by the SVI recommendations. These trending discordant intervals arose owing to the divergence of the gene-specific distribution of predictions from the genome-wide distribution, suggesting that gene-specific calibration is needed in many cases. Approximately 22% of ClinVar missense variants of uncertain significance in genes we analyzed (REVEL = 100,629, BayesDel = 71,928) had predictions in trending discordant intervals. Thus, genome-wide calibrations could result in many variants receiving inappropriate evidence strength. To facilitate a review of the SVI's calibrations, we developed a web application enabling visualization of gene-specific predictions and trending concordant and discordant intervals.

The PS4-likelihood ratio calculator: flexible allocation of evidence weighting for case-control data in variant classification

BACKGROUND

The 2015 American College of Medical Genetics/Association of Molecular Pathology (ACMG/AMP) variant classification framework specifies that case-control observations can be scored as 'strong' evidence (PS4) towards pathogenicity.

METHODS

We developed the PS4-likelihood ratio calculator (PS4-LRCalc) for quantitative evidence assignment based on the observed variant frequencies in cases and controls. Binomial likelihoods are computed for two models, each defined by prespecified OR thresholds. Model 1 represents the hypothesis of association between variant and phenotype (eg, OR≥5) and model 2 represents the hypothesis of non-association (eg, OR≤1).

RESULTS

PS4-LRCalc enables continuous quantitation of evidence for variant classification expressed as a likelihood ratio (LR), which can be log-converted into log LR (evidence points). Using PS4-LRCalc, observed data can be used to quantify evidence towards either pathogenicity or benignity. Variants can also be evaluated against models of different penetrance. The approach is applicable to balanced data sets generated for more common phenotypes and smaller data sets more typical in very rare disease variant evaluation.

CONCLUSION

PS4-LRCalc enables flexible evidence quantitation on a continuous scale for observed case-control data. The converted LR is amenable to incorporation into the now widely used 2018 updated Bayesian ACMG/AMP framework.

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

The recommendation of re-classification of variants of uncertain significance (VUS) in adult genetic disorders patients

Since variants of uncertain significance (VUS) reported in genetic testing cannot be acted upon clinically, this classification may delay or prohibit precise diagnosis and genetic counseling in adult genetic disorders patients. Large-scale analyses about qualitatively distinct lines of evidence used for VUS can make them re-classification more accurately. We analyzed 458 Chinese adult patients WES data, within 15 pathogenic evidence PS1, PS2, PM1, PM6 and PP4 were not used for VUS pathogenic classification, meanwhile the PP3, BP4, PP2 were used much more frequently. The PM2_Supporting was used most widely for all reported variants. There were also 31 null variants (nonsense, frameshift, canonical ±1 or 2 splice sites) which were probably the disease-causing variants of the patients were classified as VUS. By analyzed the evidence used for all VUS we recommend that appropriate genetic counseling, reliable releasing of in-house data, allele frequency comparison between case and control, expanded verification in patient family, co-segregation analysis and functional assays were urgent need to gather more evidence to reclassify VUS. We also found adult patients with nervous system disease were reported the most phenotype-associated VUS and the lower the phenotypic specificity, the more reported VUS. This result emphasized the importance of pretest genetic counseling which would make less reporting of VUS. Our result revealed the characteristics of the pathogenic classification evidence used for VUS in adult genetic disorders patients for the first time, recommend a rules-based process to evaluate the pathogenicity of VUS which could provide a strong basis for accurately evaluating the pathogenicity and clinical grade information of VUS. Meanwhile, we further expanded the genetic spectrum and improve the diagnostic rate of adult genetic disorders.

Prevalence of Germline Pathogenic Variants in Renal Cancer Predisposition Genes in a Population-Based Study of Renal Cell Carcinoma

Renal cell carcinoma (RCC) has been associated with germline pathogenic or likely pathogenic (PLP) variants in recognised cancer susceptibility genes. Studies of RCC using gene panel sequencing have been highly variable in terms of study design, genes included, and reported prevalence of PLP variant carriers (4-26%). Studies that restricted their analysis to established RCC predisposition genes identified variants in 1-6% of cases. This work assessed the prevalence of clinically actionable PLP variants in renal cancer predisposition genes in an Australian population-based sample of RCC cases. Germline DNA from 1029 individuals diagnosed with RCC who were recruited through the Victoria and Queensland cancer registries were screened using a custom amplicon-based panel of 21 genes. Mean age at cancer diagnosis was 60 ± 10 years, and two-thirds (690, 67%) of the participants were men. Eighteen participants (1.7%) were found to carry a PLP variant. Genes with PLP variants included BAP1, FH, FLCN, MITF, MSH6, SDHB, TSC1, and VHL. Most carriers of PLP variants did not report a family history of the disease. Further exploration of the clinical utility of gene panel susceptibility testing for all RCCs is warranted.

Splice site variants in the canonical donor site of MED13L exon 7 lead to intron retention in patients with MED13L syndrome

Pathogenic variants in the MED13L gene are associated with the autosomal dominant MED13L syndrome, which is characterised by global developmental delay and cardiac malformations. We investigated two heterozygous MED13L variants located at the canonical donor splice site motif of exon 7: c.1009+1G>C and c.1009+5G>C. We report that in silico predictions suggested two possible outcomes: exon 7 skipping, resulting in loss of the phosphodegron motif essential for MED13L regulation, or activation of a cryptic donor site in intron 7, leading to intron retention. RNA analysis confirmed that both variants affected the exon 7 splice donor site, resulting in the retention of 73 bp of intron 7. This retention caused a frameshift and premature translation termination, consistent with haploinsufficiency. Our results highlight the importance of combining predictive and experimental approaches to understand the functional impact of splice site variants. These insights into the molecular consequences of MED13L variants provide a deeper understanding of the genetic basis of MED13L syndrome.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Functional characterization of 2,832 JAG1 variants supports reclassification for Alagille syndrome and improves guidance for clinical variant interpretation

Pathogenic variants in the JAG1 gene are a primary cause of the multi-system disorder Alagille syndrome. Although variant detection rates are high for this disease, there is uncertainty associated with the classification of missense variants that leads to reduced diagnostic yield. Consequently, up to 85% of reported JAG1 missense variants have uncertain or conflicting classifications. We generated a library of 2,832 JAG1 nucleotide variants within exons 1-7, a region with a high number of reported missense variants, and designed a high-throughput assay to measure JAG1 membrane expression, a requirement for normal function. After calibration using a set of 175 known or predicted pathogenic and benign variants included within the variant library, 486 variants were characterized as functionally abnormal (n = 277 abnormal and n = 209 likely abnormal), of which 439 (90.3%) were missense. We identified divergent membrane expression occurring at specific residues, indicating that loss of the wild-type residue itself does not drive pathogenicity, a finding supported by structural modeling data and with broad implications for clinical variant classification both for Alagille syndrome and globally across other disease genes. Of 144 uncertain variants reported in patients undergoing clinical or research testing, 27 had functionally abnormal membrane expression, and inclusion of our data resulted in the reclassification of 26 to likely pathogenic. Functional evidence augments the classification of genomic variants, reducing uncertainty and improving diagnostics. Inclusion of this repository of functional evidence during JAG1 variant reclassification will significantly affect resolution of variant pathogenicity, making a critical impact on the molecular diagnosis of Alagille syndrome.

Rapid classification of a novel ALS-causing I149S variant in superoxide dismutase-1

Variants of the oxygen free radical scavenging enzyme superoxide dismutase-1 (SOD1) are associated with the neurodegenerative disease amyotrophic lateral sclerosis (ALS). These variants occur in roughly 20% of familial ALS cases, and 1% of sporadic ALS cases. Here, we identified a novel SOD1 variant in a patient in their 50s who presented with movement deficiencies and neuropsychiatric features. The variant was heterozygous and resulted in the isoleucine at position 149 being substituted with a serine (I149S). In silico analysis predicted the variant to be destabilizing to the SOD1 protein structure. Expression of the SOD1I149S variant with a C-terminal EGFP tag in neuronal-like NSC-34 cells resulted in extensive inclusion formation and reduced cell viability. Immunoblotting revealed that the intramolecular disulphide between Cys57 and Cys146 was fully reduced for SOD1I149S. Furthermore, SOD1I149S was highly susceptible to proteolytic digestion, suggesting a large degree of instability to the protein fold. Finally, fluorescence correlation spectroscopy and native-PAGE of cell lysates showed that SOD1I149S was monomeric in solution in comparison to the dimeric SOD1WT. This experimental data was obtained within 3 months and resulted in the rapid re-classification of the variant from a variant of unknown significance (VUS) to a clinically actionable likely pathogenic variant.

Scalable approaches for generating, validating and incorporating data from high-throughput functional assays to improve clinical variant classification

As the adoption and scope of genetic testing continue to expand, interpreting the clinical significance of DNA sequence variants at scale remains a formidable challenge, with a high proportion classified as variants of uncertain significance (VUSs). Genetic testing laboratories have historically relied, in part, on functional data from academic literature to support variant classification. High-throughput functional assays or multiplex assays of variant effect (MAVEs), designed to assess the effects of DNA variants on protein stability and function, represent an important and increasingly available source of evidence for variant classification, but their potential is just beginning to be realized in clinical lab settings. Here, we describe a framework for generating, validating and incorporating data from MAVEs into a semi-quantitative variant classification method applied to clinical genetic testing. Using single-cell gene expression measurements, cellular evidence models were built to assess the effects of DNA variation in 44 genes of clinical interest. This framework was also applied to models for an additional 22 genes with previously published MAVE datasets. In total, modeling data was incorporated from 24 genes into our variant classification method. These data contributed evidence for classifying 4043 observed variants in over 57,000 individuals. Genetic testing laboratories are uniquely positioned to generate, analyze, validate, and incorporate evidence from high-throughput functional data and ultimately enable the use of these data to provide definitive clinical variant classifications for more patients.

Atypical cancer risk profile in carriers of Italian founder BRCA1 variant p.His1673del: Implications for classification and clinical management

BACKGROUND

BRCA1:c.5017_5019del (p.His1673del) is a founder variant relatively frequent in Northern Italy. Despite previous suggestion of pathogenicity, variant classification in public databases is still conflicting, needing additional evidence.

METHODS

Maximum likelihood penetrance of breast/ovarian and other cancer types was estimated using full pedigree data from 53 informative Italian families. The effect of the variant on BRCA1-ABRAXAS1 interaction was assessed using a GFP-fragment reassembly-based PPI assay. Results were combined with additional data from multiple sources to classify the variant according to ACMG/AMP classification rules specified for BRCA1/2.

RESULTS

Variant-carriers displayed increased risk for ovarian cancer (HR = 33.0, 95% CI = 7.0-155.0; cumulative risk at age 70 = 27.6%, 95% CI = 12.6-40.0%) but not for breast cancer (HR = 0.7, 95% CI = 0.2-2.2). An increased risk of uterine cancer (HR = 8.0, 95% CI = 1.03-61.6) emerged, warranting further evaluation. Likelihood-ratio in favor of pathogenicity was 98898642.82 under assumption of standard BRCA1 breast and ovarian penetrance, and 104240832.84 after excluding breast cancer diagnoses (based on penetrance results). Functional analysis demonstrated that the variant abrogates the BRCA1-ABRAXAS1 binding, supporting the PS3 code assignment within the ACMG/AMP rule-based model. Collectively, these findings allowed to classify the variant as pathogenic.

CONCLUSION

Pathogenicity of BRCA1:c.5017_5019del(p.His1673del) has been confirmed; however, breast cancer risk in Italian families is not increased, unlike in families from other countries and in carriers of most BRCA1 pathogenic variants. The knowledge of atypical risk profiles for this and other variants will pave the way for personalized management based on specific genotype.

Multisite Validation of a Functional Assay to Adjudicate SCN5A Brugada Syndrome-Associated Variants

BACKGROUND

Brugada syndrome is an inheritable arrhythmia condition that is associated with rare, loss-of-function variants in SCN5A. Interpreting the pathogenicity of SCN5A missense variants is challenging, and ≈79% of SCN5A missense variants in ClinVar are currently classified as variants of uncertain significance. Automated patch clamp technology enables high-throughput functional studies of ion channel variants and can provide evidence for variant reclassification.

METHODS

An in vitro SCN5A-Brugada syndrome automated patch clamp assay was independently performed at Vanderbilt University Medical Center and Victor Chang Cardiac Research Institute. The assay was calibrated according to ClinGen Sequence Variant Interpretation recommendations using high-confidence variant controls (n=49). Normal and abnormal ranges of function were established based on the distribution of benign variant assay results. Odds of pathogenicity values were derived from the experimental results according to ClinGen Sequence Variant Interpretation recommendations. The calibrated assay was then used to study SCN5A variants of uncertain significance observed in 4 families with Brugada syndrome and other arrhythmia phenotypes associated with SCN5A loss-of-function.

RESULTS

Variant channel parameters generated independently at the 2 research sites showed strong correlations, including peak INa density (R2=0.86). The assay accurately distinguished benign controls (24/25 concordant variants) from pathogenic controls (23/24 concordant variants). Odds of pathogenicity values were 0.042 for normal function and 24.0 for abnormal function, corresponding to strong evidence for both American College of Medical Genetics and Genomics/Association for Molecular Pathology benign and pathogenic functional criteria (BS3 and PS3, respectively). Application of the assay to 4 clinical SCN5A variants of uncertain significance revealed loss-of-function for 3/4 variants, enabling reclassification to likely pathogenic.

CONCLUSIONS

This validated high-throughput assay provides clinical-grade functional evidence to aid the classification of current and future SCN5A-Brugada syndrome variants of uncertain significance.

DNA methylation analysis in patients with neurodevelopmental disorders improves variant interpretation and reveals complexity

Analysis of genomic DNA methylation by generating epigenetic signature profiles (episignatures) is increasingly being implemented in genetic diagnosis. Here we report our experience using episignature analysis to resolve both uncomplicated and complex cases of neurodevelopmental disorders (NDDs). We analyzed 97 NDDs divided into (1) a validation cohort of 59 patients with likely pathogenic/pathogenic variants characterized by a known episignature and (2) a test cohort of 38 patients harboring variants of unknown significance or unidentified variants. The expected episignature was obtained in most cases with likely pathogenic/pathogenic variants (53/59 [90%]), a revealing exception being the overlapping profile of two SMARCB1 pathogenic variants with ARID1A/B:c.6200, confirmed by the overlapping clinical features. In the test cohort, five cases showed the expected episignature, including (1) novel pathogenic variants in ARID1B and BRWD3; (2) a deletion in ATRX causing MRXFH1 X-linked mental retardation; and (3) confirmed the clinical diagnosis of Cornelia de Lange (CdL) syndrome in mutation-negative CdL patients. Episignatures analysis of the in BAF complex components revealed novel functional protein interactions and common episignatures affecting homologous residues in highly conserved paralogous proteins (SMARCA2 M856V and SMARCA4 M866V). Finally, we also found sex-dependent episignatures in X-linked disorders. Implementation of episignature profiling is still in its early days, but with increasing utilization comes increasing awareness of the capacity of this methodology to help resolve the complex challenges of genetic diagnoses.

Integrative analysis of Lunatic Fringe variants associated with spondylocostal dysostosis type-III

Lunatic Fringe (LFNG) is required for spinal development. Biallelic pathogenic variants cause spondylocostal dysostosis type-III (SCD3), a rare disease generally characterized by malformed, asymmetrical, and attenuated development of the vertebral column and ribs. However, a variety of SCD3 cases reported have presented with additional features such as auditory alterations and digit abnormalities. There has yet to be a single, comprehensive, functional evaluation of causative LFNG variants and such analyses could unveil molecular mechanisms for phenotypic variability in SCD3. Therefore, nine LFNG missense variants associated with SCD3, c.564C>A, c.583T>C, c.842C>A, c.467T>G, c.856C>T, c.601G>A, c.446C>T, c.521G>A, and c.766G>A, were assessed in vitro for subcellular localization and protein processing. Glycosyltransferase activity was quantified for the first time in the c.583T>C, c.842C>A, and c.446C>T variants. Primarily, our results are the first to satisfy American College of Medical Genetics and Genomics PS3 criteria (functional evidence via well-established assay) for the pathogenicity of c.583T>C, c.842C>A, and c.446C>T, and replicate this evidence for the remaining six variants. Secondly, this work indicates that all variants that prevent Golgi localization also lead to impaired protein processing. It appears that the FRINGE domain is responsible for this phenomenon. Thirdly, our data suggests that variant proximity to the catalytic residue may influence whether LFNG is improperly trafficked and/or enzymatically dysfunctional. Finally, the phenotype of the axial skeleton, but not elsewhere, may be modulated in a variant-specific fashion. More reports are needed to continue testing this hypothesis. We anticipate our data will be used as a basis for discussion of genotype-phenotype correlations in SCD3.

Developing a genetic testing panel for evaluation of morbidities in kidney transplant recipients

Cardiovascular disease, infection, malignancy, and thromboembolism are major causes of morbidity and mortality in kidney transplant recipients (KTR). Prospectively identifying monogenic conditions associated with post-transplant complications may enable personalized management. Therefore, we developed a transplant morbidity panel (355 genes) associated with major post-transplant complications including cardiometabolic disorders, immunodeficiency, malignancy, and thrombophilia. This gene panel was then evaluated using exome sequencing data from 1590 KTR. Additionally, genes associated with monogenic kidney and genitourinary disorders along with American College of Medical Genetics (ACMG) secondary findings v3.2 were annotated. Altogether, diagnostic variants in 37 genes associated with Mendelian kidney and genitourinary disorders were detected in 9.9% (158/1590) of KTR; 25.9% (41/158) had not been clinically diagnosed. Moreover, the transplant morbidity gene panel detected diagnostic variants for 56 monogenic disorders in 9.1% KTRs (144/1590). Cardiovascular disease, malignancy, immunodeficiency, and thrombophilia variants were detected in 5.1% (81), 2.1% (34), 1.8% (29) and 0.2% (3) among 1590 KTRs, respectively. Concordant phenotypes were present in half of these cases. Reviewing implications for transplant care, these genetic findings would have allowed physicians to set specific risk factor targets in 6.3% (9/144), arrange intensive surveillance in 97.2% (140/144), utilize preventive measures in 13.2% (19/144), guide disease-specific therapy in 63.9% (92/144), initiate specialty referral in 90.3% (130/144) and alter immunosuppression in 56.9% (82/144). Thus, beyond diagnostic testing for kidney disorders, sequence annotation identified monogenic disorders associated with common post-transplant complications in 9.1% of KTR, with important clinical implications. Incorporating genetic diagnostics for transplant morbidities would enable personalized management in pre- and post-transplant care.

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

A new way of looking at transcription factor assays

AF Rubin discusses a new high-throughput functional assay for transcription factors applied for a deep mutational scanning study of the transcription factor PAX6 by Kudla and colleagues (McDonnell et al, 2023 ) in this issue of Molecular Systems Biology .

Saturation genome editing of BAP1 functionally classifies somatic and germline variants

Many variants that we inherit from our parents or acquire de novo or somatically are rare, limiting the precision with which we can associate them with disease. We performed exhaustive saturation genome editing (SGE) of BAP1, the disruption of which is linked to tumorigenesis and altered neurodevelopment. We experimentally characterized 18,108 unique variants, of which 6,196 were found to have abnormal functions, and then used these data to evaluate phenotypic associations in the UK Biobank. We also characterized variants in a large population-ascertained tumor collection, in cancer pedigrees and ClinVar, and explored the behavior of cancer-associated variants compared to that of variants linked to neurodevelopmental phenotypes. Our analyses demonstrated that disruptive germline BAP1 variants were significantly associated with higher circulating levels of the mitogen IGF-1, suggesting a possible pathological mechanism and therapeutic target. Furthermore, we built a variant classifier with >98% sensitivity and specificity and quantify evidence strengths to aid precision variant interpretation.

GM1 gangliosidosis type II: Results of a 10-year prospective study

PURPOSE

GM1 gangliosidosis (GM1) a lysosomal disorder caused by pathogenic variants in GLB1, is characterized by relentless neurodegeneration. There are no approved treatments.

METHODS

Forty-one individuals with type II (late-infantile and juvenile) GM1 participated in a single-site prospective observational study.

RESULTS

Classification of 37 distinct variants using American College of Medical Genetics and Genomics criteria resulted in the upgrade of 6 and the submission of 4 new variants. In contrast to type I infantile disease, children with type II had normal or near normal hearing and did not have cherry-red maculae or hepatosplenomegaly. Some older children with juvenile onset disease developed thickened aortic and/or mitral valves. Serial magnetic resonance images demonstrated progressive brain atrophy, more pronounced in late infantile patients. Magnetic resonance spectroscopy showed worsening elevation of myo-inositol and deficit of N-acetyl aspartate that were strongly correlated with scores on the Vineland Adaptive Behavior Scale, progressing more rapidly in late infantile compared with juvenile onset disease.

CONCLUSION

Serial phenotyping of type II GM1 patients expands the understanding of disease progression and clarifies common misconceptions about type II patients; these are pivotal steps toward more timely diagnosis and better supportive care. The data amassed through this 10-year effort will serve as a robust comparator for ongoing and future therapeutic trials.

Clinical Impact of Telomere Length Testing for Interstitial Lung Disease

BACKGROUND

Shortened telomere length (TL) is a genomic risk factor for fibrotic interstitial lung disease (ILD), but its role in clinical management is unknown.

RESEARCH QUESTION

What is the clinical impact of TL testing on the management of ILD?

STUDY DESIGN AND METHODS

Patients were evaluated in the Columbia University ILD clinic and underwent Clinical Laboratory Improvement Amendments-certified TL testing by flow cytometry and fluorescence in situ hybridization (FlowFISH) as part of clinical treatment. Short TL was defined as below the 10th age-adjusted percentile for either granulocytes or lymphocytes by FlowFISH. Patients were offered genetic counseling and testing if they had short TL or a family history of ILD. FlowFISH TL was compared with research quantitative polymerase chain reaction (qPCR) TL measurement.

RESULTS

A total of 108 patients underwent TL testing, including those with clinical features of short telomere syndrome such as familial pulmonary fibrosis (50%) or extrapulmonary manifestations in the patient (25%) or a relative (41%). The overall prevalence of short TL was 46% and was similar across clinical ILD diagnoses. The number of short telomere clinical features was independently associated with detecting short TL (OR, 2.00; 95% CI, 1.27-3.32). TL testing led to clinical treatment changes for 35 patients (32%), most commonly resulting in reduction or avoidance of immunosuppression. Of the patients who underwent genetic testing (n = 34), a positive or candidate diagnostic finding in telomere-related genes was identified in 10 patients (29%). Inclusion of TL testing below the 1st percentile helped reclassify eight of nine variants of uncertain significance into actionable findings. The qPCR test correlated with FlowFISH, but age-adjusted percentile cutoffs may not be equivalent between the two assays.

INTERPRETATION

Incorporating TL testing in ILD impacted clinical management and led to the discovery of new actionable genetic variants.

Multiplexed Assays of Variant Effect and Automated Patch-clamping Improve KCNH2-LQTS Variant Classification and Cardiac Event Risk Stratification

Background

Long QT syndrome (LQTS) is a lethal arrhythmia syndrome, frequently caused by rare loss-of-function variants in the potassium channel encoded by KCNH2. Variant classification is difficult, often owing to lack of functional data. Moreover, variant-based risk stratification is also complicated by heterogenous clinical data and incomplete penetrance. Here, we sought to test whether variant-specific information, primarily from high-throughput functional assays, could improve both classification and cardiac event risk stratification in a large, harmonized cohort of KCNH2 missense variant heterozygotes.

Methods

We quantified cell-surface trafficking of 18,796 variants in KCNH2 using a Multiplexed Assay of Variant Effect (MAVE). We recorded KCNH2 current density for 533 variants by automated patch clamping (APC). We calibrated the strength of evidence of MAVE data according to ClinGen guidelines. We deeply phenotyped 1,458 patients with KCNH2 missense variants, including QTc, cardiac event history, and mortality. We correlated variant functional data and Bayesian LQTS penetrance estimates with cohort phenotypes and assessed hazard ratios for cardiac events.

Results

Variant MAVE trafficking scores and APC peak tail currents were highly correlated (Spearman Rank-order ρ = 0.69). The MAVE data were found to provide up to pathogenic very strong evidence for severe loss-of-function variants. In the cohort, both functional assays and Bayesian LQTS penetrance estimates were significantly predictive of cardiac events when independently modeled with patient sex and adjusted QT interval (QTc); however, MAVE data became non-significant when peak-tail current and penetrance estimates were also available. The area under the ROC for 20-year event outcomes based on patient-specific sex and QTc (AUC 0.80 [0.76-0.83]) was improved with prospectively available penetrance scores conditioned on MAVE (AUC 0.86 [0.83-0.89]) or attainable APC peak tail current data (AUC 0.84 [0.81-0.88]).

Conclusion

High throughput KCNH2 variant MAVE data meaningfully contribute to variant classification at scale while LQTS penetrance estimates and APC peak tail current measurements meaningfully contribute to risk stratification of cardiac events in patients with heterozygous KCNH2 missense variants.

Multimodal framework to resolve variants of uncertain significance in TSC2

Efforts to resolve the functional impact of variants of uncertain significance (VUS) have lagged behind the identification of new VUS; as such, there is a critical need for scalable VUS resolution technologies. Computational variant effect predictors (VEPs), once trained, can predict pathogenicity for all missense variants in a gene, set of genes, or the exome. Existing tools have employed information on known pathogenic and benign variants throughout the genome to predict pathogenicity of VUS. We hypothesize that taking a gene-specific approach will improve pathogenicity prediction over globally-trained VEPs. We tested this hypothesis using the gene TSC2, whose loss of function results in tuberous sclerosis, a multisystem mTORopathy affecting about 1 in 6,000 individuals born in the United States. TSC2 has been identified as a high-priority target for VUS resolution, with (1) well-characterized molecular and patient phenotypes associated with loss-of-function variants, and (2) more than 2,700 VUS already documented in ClinVar. We developed Tuberous sclerosis classifier to Resolve variants of Uncertain Significance in T SC2 (TRUST), a machine learning model to predict pathogenicity of TSC2 missense VUS. To test whether these predictions are accurate, we further introduce curated loci prime editing (cliPE) as an accessible strategy for performing scalable multiplexed assays of variant effect (MAVEs). Using cliPE, we tested the effects of more than 200 TSC2 variants, including 106 VUS. It is highly likely this functional data alone would be sufficient to reclassify 92 VUS with most being reclassified as likely benign. We found that TRUST's classifications were correlated with the functional data, providing additional validation for the in silico predictions. We provide our pathogenicity predictions and MAVE data to aid with VUS resolution. In the near future, we plan to host these data on a public website and deposit into relevant databases such as MAVEdb as a community resource. Ultimately, this study provides a framework to complete variant effect maps of TSC1 and TSC2 and adapt this approach to other mTORopathy genes.

Neuropathy target esterase activity defines phenotypes among PNPLA6 disorders

Biallelic pathogenic variants in the PNPLA6 gene cause a broad spectrum of disorders leading to gait disturbance, visual impairment, anterior hypopituitarism and hair anomalies. PNPLA6 encodes neuropathy target esterase (NTE), yet the role of NTE dysfunction on affected tissues in the large spectrum of associated disease remains unclear. We present a systematic evidence-based review of a novel cohort of 23 new patients along with 95 reported individuals with PNPLA6 variants that implicate missense variants as a driver of disease pathogenesis. Measuring esterase activity of 46 disease-associated and 20 common variants observed across PNPLA6-associated clinical diagnoses unambiguously reclassified 36 variants as pathogenic and 10 variants as likely pathogenic, establishing a robust functional assay for classifying PNPLA6 variants of unknown significance. Estimating the overall NTE activity of affected individuals revealed a striking inverse relationship between NTE activity and the presence of retinopathy and endocrinopathy. This phenomenon was recaptured in vivo in an allelic mouse series, where a similar NTE threshold for retinopathy exists. Thus, PNPLA6 disorders, previously considered allelic, are a continuous spectrum of pleiotropic phenotypes defined by an NTE genotype:activity:phenotype relationship. This relationship, and the generation of a preclinical animal model, pave the way for therapeutic trials, using NTE as a biomarker.

ATP1A3 Disease Spectrum Includes Paroxysmal Weakness and Encephalopathy Not Triggered by Fever

Background and Objectives

Heterozygous pathogenic variants in ATP1A3, which encodes the catalytic alpha subunit of neuronal Na+/K+-ATPase, cause primarily neurologic disorders with widely variable features that can include episodic movement deficits. One distinctive presentation of ATP1A3-related disease is recurrent fever-triggered encephalopathy. This can occur with generalized weakness and/or ataxia and is described in the literature as relapsing encephalopathy with cerebellar ataxia. This syndrome displays genotype-phenotype correlation with variants at p.R756 causing temperature sensitivity of ATP1A3. We report clinical and in vitro functional evidence for a similar phenotype not triggered by fever but associated with protein loss-of-function.

Methods

We describe the phenotype of an individual with de novo occurrence of a novel heterozygous ATP1A3 variant, NM_152296.5:c.388_390delGTG; p.(V130del). We confirmed the pathogenicity of p.V130del by cell survival complementation assay in HEK293 cells and then characterized its functional impact on enzymatic ion transport and extracellular sodium binding by two-electrode voltage clamp electrophysiology in Xenopus oocytes. To determine whether variant enzymes reach the cell surface, we surface-biotinylated oocytes expressing N-tagged ATP1A3.

Results

The proband is a 7-year-old boy who has had 2 lifetime episodes of paroxysmal weakness, encephalopathy, and ataxia not triggered by fever. He had speech regression and intermittent hand tremors after the second episode but otherwise spontaneously recovered after episodes and is at present developmentally appropriate. The p.V130del variant was identified on clinical trio exome sequencing, which did not reveal any other variants possibly associated with the phenotype. p.V130del eliminated ATP1A3 function in cell survival complementation assay. In Xenopus oocytes, p.V130del variant Na+/K+-ATPases showed complete loss of ion transport activity and marked abnormalities of extracellular Na+ binding at room temperature. Despite this clear loss-of-function effect, surface biotinylation under the same conditions revealed that p.V130del variant enzymes were still present at the oocyte's cell membrane.

Discussion

This individual's phenotype expands the clinical spectrum of ATP1A3-related recurrent encephalopathy to include presentations without fever-triggered events. The total loss of ion transport function with p.V130del, despite enzyme presence at the cell membrane, indicates that haploinsufficiency can cause relatively mild phenotypes in ATP1A3-related disease.

Cosmic Whirl: Navigating the Comet Trail in DNA: H2AX Phosphorylation and the Enigma of Uncertain Significance Variants

Pathogenic variations in the BRCA2 gene have been detected with the development of next-generation sequencing (NGS)-based hereditary cancer panel testing technology. It also reveals an increasing number of variants of uncertain significance (VUSs). Well-established functional tests are crucial to accurately reclassifying VUSs for effective diagnosis and treatment. We retrospectively analyzed the multi-gene cancer panel results of 922 individuals and performed in silico analysis following ClinVar classification. Then, we selected five breast cancer-diagnosed patients' missense BRCA2 VUSs (T1011R, T1104P/M1168K, R2027K, G2044A, and D2819) for reclassification. The effects of VUSs on BRCA2 function were analyzed using comet and H2AX phosphorylation (γH2AX) assays before and after the treatment of peripheral blood mononuclear cells (PBMCs) of subjects with the double-strand break (DSB) agent doxorubicin (Dox). Before and after Dox-induction, the amount of DNA in the comet tails was similar in VUS carriers; however, notable variations in γH2AX were observed, and according to combined computational and functional analyses, we reclassified T1001R as VUS-intermediate, T1104P/M1168K and D2819V as VUS (+), and R2027K and G2044A as likely benign. These findings highlight the importance of the variability of VUSs in response to DNA damage before and after Dox-induction and suggest that further investigation is needed to understand the underlying mechanisms.

Prioritize Variant Reclassification in Pediatric Long QT Syndrome-Time to Revisit

Congenital long QT syndrome (LQTS) is an inherited arrhythmia syndrome associated with sudden cardiac death. Accurate interpretation and classification of genetic variants in LQTS patients are crucial for effective management. All patients with LQTS with a positive genetic test over the past 18 years (2002-2020) in our single tertiary pediatric cardiac center were identified. Reevaluation of the reported variants in LQTS genes was conducted using the American College of Genetics and Genomics (ACMG) guideline after refinement by the US ClinGen SVI working group and guideline by Walsh et al. on genetic variant reclassification, under multidisciplinary input. Among the 59 variants identified. 18 variants (30.5%) were reclassified. A significant larger portion of variants of unknown significance (VUS) were reclassified compared to likely pathogenic (LP)/pathogenic (P) variants (57.7% vs 9.1%, p < 0.001). The rate of reclassification was significantly higher in the limited/disputed evidence group compared to the definite/moderate evidence group (p = 0.0006). All LP/P variants were downgraded in the limited/disputed evidence group (p = 0.0057). VUS upgrades are associated with VUS located in genes within the definite/moderate evidence group (p = 0.0403) and with VUS present in patients exhibiting higher corrected QT intervals (QTc) (p = 0.0445). A significant number of pediatric LQTS variants were reclassified, particularly for VUS. The strength of the gene-disease association of the genes influences the reclassification performance. The study provides important insights and guidance for pediatricians to seek for reclassification of "outdated variants" in order to facilitate contemporary precision medicine.

Systematic analysis of NDUFAF6 in complex I assembly and mitochondrial disease

Isolated complex I (CI) deficiencies are a major cause of primary mitochondrial disease. A substantial proportion of CI deficiencies are believed to arise from defects in CI assembly factors (CIAFs) that are not part of the CI holoenzyme. The biochemistry of these CIAFs is poorly defined, making their role in CI assembly unclear, and confounding interpretation of potential disease-causing genetic variants. To address these challenges, we devised a deep mutational scanning approach to systematically assess the function of thousands of NDUFAF6 genetic variants. Guided by these data, biochemical analyses and cross-linking mass spectrometry, we discovered that the CIAF NDUFAF6 facilitates incorporation of NDUFS8 into CI and reveal that NDUFS8 overexpression rectifies NDUFAF6 deficiency. Our data further provide experimental support of pathogenicity for seven novel NDUFAF6 variants associated with human pathology and introduce functional evidence for over 5,000 additional variants. Overall, our work defines the molecular function of NDUFAF6 and provides a clinical resource for aiding diagnosis of NDUFAF6-related diseases.

In Vitro Functional Analysis Can Aid Precision Diagnostics of HNF1B-MODY

Precision medicine relies on accurate and consistent classification of sequence variants. A correct diagnosis of hepatocyte nuclear factor (HNF) 1B maturity-onset diabetes of the young, caused by pathogenic variants in the HNF1B gene, is important for optimal disease management and prognosis, and it has implications for genetic counseling and follow-up of at-risk family members. We hypothesized that the functional characterization could provide valuable information to assist the interpretation of pathogenicity of HNF1B variants. Using different in vitro functional assays, variants identified among 313 individuals, suspected to have monogenic diabetes with or without kidney disease, were characterized. The data from the functional assays were subsequently conjugated with obtained clinical, biochemical, and in silico data. Two variants (p.A167P, p.H336Pfs∗22) showed severe loss of function due to impaired transactivation, reduced DNA binding (p.A167P), and mRNA instability (p.A167P). Although both these variant carriers were diagnosed with diabetes, the p.H336Pfs∗22 carrier also had congenital absence of a kidney, which is a characteristic trait for HNF1B maturity-onset diabetes of the young. Functional analysis of the p.A167P variant revealed damaging effects on HNF-1B protein function, which may warrant imaging of the kidneys and/or pancreas. In addition, the current study has generated important data, including evidence supporting the benign functional impact of five variants (p.D82N, p.T88A, p.N394D, p.V458G, and p.T544A), and piloting new approaches that will prove critical for the growth of HNF1B-diabetes diagnosis.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

BACKGROUND

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility.

METHODS

In this study, we leveraged the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein-coding KCNE1 variants.

RESULTS

We comprehensively assayed KCNE1 variant cell surface expression (2554/2709 possible single-amino-acid variants) and function (2534 variants). Our study identified 470 loss- or partial loss-of-surface expression and 574 loss- or partial loss-of-function variants. Of the 574 loss- or partial loss-of-function variants, 152 (26.5%) had reduced cell surface expression, indicating that most functionally deleterious variants affect channel gating. Nonsense variants at residues 56-104 generally had WT-like trafficking scores but decreased functional scores, indicating that the latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation (with > 70% loss-of-function variants) were in predicted close contact with binding partners KCNQ1 or calmodulin. Our functional assay data were consistent with gold standard electrophysiological data (ρ =  - 0.64), population and patient cohorts (32/38 presumed benign or pathogenic variants with consistent scores), and computational predictors (ρ =  - 0.62). Our data provide moderate-strength evidence for the American College of Medical Genetics/Association of Molecular Pathology functional criteria for benign and pathogenic variants.

CONCLUSIONS

Comprehensive variant effect maps of KCNE1 can both provide insight into I Ks channel biology and help reclassify variants of uncertain significance.

Specifications of the ACMG/AMP variant curation guidelines for the analysis of germline ATM sequence variants

The ClinGen Hereditary Breast, Ovarian and Pancreatic Cancer (HBOP) Variant Curation Expert Panel (VCEP) is composed of internationally recognized experts in clinical genetics, molecular biology and variant interpretation. This VCEP made specifications for ACMG/AMP guidelines for the ataxia telangiectasia mutated (ATM) gene according to the Food and Drug Administration (FDA)-approved ClinGen protocol. These gene-specific rules for ATM were modified from the American College of Medical Genetics and Association for Molecular Pathology (ACMG/AMP) guidelines and were tested against 33 ATM variants of various types and classifications in a pilot curation phase. The pilot revealed a majority agreement between the HBOP VCEP classifications and the ClinVar-deposited classifications. Six pilot variants had conflicting interpretations in ClinVar and reevaluation with the VCEP's ATM-specific rules resulted in four that were classified as benign, one as likely pathogenic and one as a variant of uncertain significance (VUS) by the VCEP, improving the certainty of interpretations in the public domain. Overall, 28 the 33 pilot variants were not VUS leading to an 85% classification rate. The ClinGen-approved, modified rules demonstrated value for improved interpretation of variants in ATM.

Genetics of inherited peripheral neuropathies and the next frontier: looking backwards to progress forwards

Inherited peripheral neuropathies (IPNs) encompass a clinically and genetically heterogeneous group of disorders causing length-dependent degeneration of peripheral autonomic, motor and/or sensory nerves. Despite gold-standard diagnostic testing for pathogenic variants in over 100 known associated genes, many patients with IPN remain genetically unsolved. Providing patients with a diagnosis is critical for reducing their 'diagnostic odyssey', improving clinical care, and for informed genetic counselling. The last decade of massively parallel sequencing technologies has seen a rapid increase in the number of newly described IPN-associated gene variants contributing to IPN pathogenesis. However, the scarcity of additional families and functional data supporting variants in potential novel genes is prolonging patient diagnostic uncertainty and contributing to the missing heritability of IPNs. We review the last decade of IPN disease gene discovery to highlight novel genes, structural variation and short tandem repeat expansions contributing to IPN pathogenesis. From the lessons learnt, we provide our vision for IPN research as we anticipate the future, providing examples of emerging technologies, resources and tools that we propose that will expedite the genetic diagnosis of unsolved IPN families.

Analyzing the functional effects of DNA variants with gene editing

Continual advancements in genomics have led to an ever-widening disparity between the rate of discovery of genetic variants and our current understanding of their functions and potential roles in disease. Systematic methods for phenotyping DNA variants are required to effectively translate genomics data into improved outcomes for patients with genetic diseases. To make the biggest impact, these approaches must be scalable and accurate, faithfully reflect disease biology, and define complex disease mechanisms. We compare current methods to analyze the function of variants in their endogenous DNA context using genome editing strategies, such as saturation genome editing, base editing and prime editing. We discuss how these technologies can be linked to high-content readouts to gain deep mechanistic insights into variant effects. Finally, we highlight key challenges that need to be addressed to bridge the genotype to phenotype gap, and ultimately improve the diagnosis and treatment of genetic diseases.

Explaining Alport syndrome-lessons from the adult nephrology clinic

Alport syndrome is a genetic kidney disease that causes worsening of kidney function over time, often progressing to kidney failure. Some types of Alport syndrome cause other symptoms and signs, including hearing loss and eye abnormalities. Research now indicates that Alport syndrome (autosomal dominant inheritance) is the most common form. Alport syndrome can have X-linked or a rare form of autosomal recessive inheritance. Traditionally, a kidney biopsy was used to diagnose Alport syndrome, but genetic testing provides a more precise and less invasive means of diagnosis and reveals the underlying pattern of inheritance. At present, there are no specific curative treatments for Alport syndrome however there is a strong international effort in pursuit of future therapies. Currently, angiotensin-converting enzyme inhibitors (ACEi), or an angiotensin receptor blocker (ARB) if a patient cannot tolerate an ACEi, slow down the progression of kidney disease and can delay the onset of kidney failure by years. There are other potential treatments in research that potentially can help delay the onset of kidney issues. Early treatment of patients and identification of their at-risk relatives is a priority. People living with Alport syndrome and their doctors now benefit from an active international research community working on translating further treatments into clinical practice and providing up-to-date clinical guidelines.

Systematic large-scale application of ClinGen InSiGHT APC -specific ACMG/AMP variant classification criteria substantially alleviates the burden of variants of uncertain significance in ClinVar and LOVD databases

Background

Pathogenic constitutional APC variants underlie familial adenomatous polyposis, the most common hereditary gastrointestinal polyposis syndrome. To improve variant classification and resolve the interpretative challenges of variants of uncertain significance (VUS), APC-specific ACMG/AMP variant classification criteria were developed by the ClinGen-InSiGHT Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP).

Methods

A streamlined algorithm using the APC -specific criteria was developed and applied to assess all APC variants in ClinVar and the InSiGHT international reference APC LOVD variant database.

Results

A total of 10,228 unique APC variants were analysed. Among the ClinVar and LOVD variants with an initial classification of (Likely) Benign or (Likely) Pathogenic, 94% and 96% remained in their original categories, respectively. In contrast, 41% ClinVar and 61% LOVD VUS were reclassified into clinically actionable classes, the vast majority as (Likely) Benign. The total number of VUS was reduced by 37%. In 21 out of 36 (58%) promising APC variants that remained VUS despite evidence for pathogenicity, a data mining-driven work-up allowed their reclassification as (Likely) Pathogenic.

Conclusions

The application of APC -specific criteria substantially reduced the number of VUS in ClinVar and LOVD. The study also demonstrated the feasibility of a systematic approach to variant classification in large datasets, which might serve as a generalisable model for other gene-/disease-specific variant interpretation initiatives. It also allowed for the prioritization of VUS that will benefit from in-depth evidence collection. This subset of APC variants was approved by the VCEP and made publicly available through ClinVar and LOVD for widespread clinical use.

Diagnostic utility and reporting recommendations for clinical DNA methylation episignature testing in genetically undiagnosed rare diseases

PURPOSE

This study aims to assess the diagnostic utility and provide reporting recommendations for clinical DNA methylation episignature testing based on the cohort of patients tested through the EpiSign Clinical Testing Network.

METHODS

The EpiSign assay utilized unsupervised clustering techniques and a support vector machine-based classification algorithm to compare each patient's genome-wide DNA methylation profile with the EpiSign Knowledge Database, yielding the result that was reported. An international working group, representing distinct EpiSign Clinical Testing Network health jurisdictions, collaborated to establish recommendations for interpretation and reporting of episignature testing.

RESULTS

Among 2399 cases analyzed, 1667 cases underwent a comprehensive screen of validated episignatures, imprinting, and promoter regions, resulting in 18.7% (312/1667) positive reports. The remaining 732 referrals underwent targeted episignature analysis for assessment of sequence or copy-number variants (CNVs) of uncertain significance or for assessment of clinical diagnoses without confirmed molecular findings, and 32.4% (237/732) were positive. Cases with detailed clinical information were highlighted to describe various utility scenarios for episignature testing.

CONCLUSION

Clinical DNA methylation testing including episignatures, imprinting, and promoter analysis provided by an integrated network of clinical laboratories enables test standardization and demonstrates significant diagnostic yield and clinical utility beyond DNA sequence analysis in rare diseases.

Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification

CRISPR base editing screens enable analysis of disease-associated variants at scale; however, variable efficiency and precision confounds the assessment of variant-induced phenotypes. Here, we provide an integrated experimental and computational pipeline that improves estimation of variant effects in base editing screens. We use a reporter construct to measure guide RNA (gRNA) editing outcomes alongside their phenotypic consequences and introduce base editor screen analysis with activity normalization (BEAN), a Bayesian network that uses per-guide editing outcomes provided by the reporter and target site chromatin accessibility to estimate variant impacts. BEAN outperforms existing tools in variant effect quantification. We use BEAN to pinpoint common regulatory variants that alter low-density lipoprotein (LDL) uptake, implicating previously unreported genes. Additionally, through saturation base editing of LDLR, we accurately quantify missense variant pathogenicity that is consistent with measurements in UK Biobank patients and identify underlying structural mechanisms. This work provides a widely applicable approach to improve the power of base editing screens for disease-associated variant characterization.

Monogenic lupus: insights into disease pathogenesis and therapeutic opportunities

PURPOSE OF REVIEW

This review aims to provide an overview of the genes and molecular pathways involved in monogenic lupus, the implications for genome diagnosis, and the potential therapies targeting these molecular mechanisms.

RECENT FINDINGS

To date, more than 30 genes have been identified as contributors to monogenic lupus. These genes are primarily related to complement deficiency, activation of the type I interferon (IFN) pathway, disruption of B-cell and T-cell tolerance and metabolic pathways, which reveal the multifaceted nature of systemic lupus erythematosus (SLE) pathogenesis.

SUMMARY

In-depth study of the causes of monogenic lupus can provide valuable insights into of pathogenic mechanisms of SLE, facilitate the identification of effective biomarkers, and aid in developing therapeutic strategies.

The importance of variant reinterpretation in inherited cardiovascular diseases: Establishing the optimal timeframe

Inherited cardiovascular diseases are rare diseases that are difficult to diagnose by non-expert professionals. Genetic analyses play a key role in the diagnosis of these diseases, in which the identification of a pathogenic genetic variant is often a diagnostic criterion. Therefore, genetic variant classification and routine reinterpretation as data become available represent one of the main challenges associated with genetic analyses. Using the genetic variants identified in an inherited cardiovascular diseases unit during a 10-year period, the objectives of this study were: 1) to evaluate the impact of genetic variant reinterpretation, 2) to compare the reclassification rates between different cohorts of cardiac channelopathies and cardiomyopathies, and 3) to establish the most appropriate periodicity for genetic variant reinterpretation. All the evaluated cohorts (full cohort of inherited cardiovascular diseases, cardiomyopathies, cardiac channelopathies, hypertrophic cardiomyopathy, dilated cardiomyopathy, arrhythmogenic cardiomyopathy, Brugada syndrome, long QT syndrome and catecholaminergic polymorphic ventricular tachycardia) showed reclassification rates above 25%, showing even higher reclassification rates when there is definitive evidence of the association between the gene and the disease in the cardiac channelopathies. Evaluation of genetic variant reclassification rates based on the year of the initial classification showed that the most appropriate frequency for the reinterpretation would be 2 years, with the possibility of a more frequent reinterpretation if deemed convenient. To keep genetic variant classifications up to date, genetic counsellors play a critical role in the reinterpretation process, providing clinical evidence that genetic diagnostic laboratories often do not have at their disposal and communicating changes in classification and the potential implications of these reclassifications to patients and relatives.

Assessment of BRCA1 and BRCA2 Germline Variant Data From Patients With Breast Cancer in a Real-World Data Registry

PURPOSE

The emergence of large real-world clinical databases and tools to mine electronic medical records has allowed for an unprecedented look at large data sets with clinical and epidemiologic correlates. In clinical cancer genetics, real-world databases allow for the investigation of prevalence and effectiveness of prevention strategies and targeted treatments and for the identification of barriers to better outcomes. However, real-world data sets have inherent biases and problems (eg, selection bias, incomplete data, measurement error) that may hamper adequate analysis and affect statistical power.

METHODS

Here, we leverage a real-world clinical data set from a large health network for patients with breast cancer tested for variants in BRCA1 and BRCA2 (N = 12,423). We conducted data cleaning and harmonization, cross-referenced with publicly available databases, performed variant reassessment and functional assays, and used functional data to inform a variant's clinical significance applying American College of Medical Geneticists and the Association of Molecular Pathology guidelines.

RESULTS

In the cohort, White and Black patients were over-represented, whereas non-White Hispanic and Asian patients were under-represented. Incorrect or missing variant designations were the most significant contributor to data loss. While manual curation corrected many incorrect designations, a sizable fraction of patient carriers remained with incorrect or missing variant designations. Despite the large number of patients with clinical significance not reported, original reported clinical significance assessments were accurate. Reassessment of variants in which clinical significance was not reported led to a marked improvement in data quality.

CONCLUSION

We identify the most common issues with BRCA1 and BRCA2 testing data entry and suggest approaches to minimize data loss and keep interpretation of clinical significance of variants up to date.

Defining and Reducing Variant Classification Disparities

Background

Multiplexed Assays of Variant Effects (MAVEs) can test all possible single variants in a gene of interest. The resulting saturation-style data may help resolve variant classification disparities between populations, especially for variants of uncertain significance (VUS).

Methods

We analyzed clinical significance classifications in 213,663 individuals of European-like genetic ancestry versus 206,975 individuals of non-European-like genetic ancestry from All of Us and the Genome Aggregation Database. Then, we incorporated clinically calibrated MAVE data into the Clinical Genome Resource's Variant Curation Expert Panel rules to automate VUS reclassification for BRCA1, TP53, and PTEN .

Results

Using two orthogonal statistical approaches, we show a higher prevalence ( p ≤5.95e-06) of VUS in individuals of non-European-like genetic ancestry across all medical specialties assessed in all three databases. Further, in the non-European-like genetic ancestry group, higher rates of Benign or Likely Benign and variants with no clinical designation ( p ≤2.5e-05) were found across many medical specialties, whereas Pathogenic or Likely Pathogenic assignments were higher in individuals of European-like genetic ancestry ( p ≤2.5e-05). Using MAVE data, we reclassified VUS in individuals of non-European-like genetic ancestry at a significantly higher rate in comparison to reclassified VUS from European-like genetic ancestry ( p =9.1e-03) effectively compensating for the VUS disparity. Further, essential code analysis showed equitable impact of MAVE evidence codes but inequitable impact of allele frequency ( p =7.47e-06) and computational predictor ( p =6.92e-05) evidence codes for individuals of non-European-like genetic ancestry.

Conclusions

Generation of saturation-style MAVE data should be a priority to reduce VUS disparities and produce equitable training data for future computational predictors.

Clinical interpretation of KCNH2 variants using a robust PS3/BS3 functional patch-clamp assay

Long QT syndrome (LQTS), caused by the dysfunction of cardiac ion channels, increases the risk of sudden death in otherwise healthy young people. For many variants in LQTS genes, there is insufficient evidence to make a definitive genetic diagnosis. We have established a robust functional patch-clamp assay to facilitate classification of missense variants in KCNH2, one of the key LQTS genes. A curated set of 30 benign and 30 pathogenic missense variants were used to establish the range of normal and abnormal function. The extent to which variants reduced protein function was quantified using Z scores, the number of standard deviations from the mean of the normalized current density of the set of benign variant controls. A Z score of -2 defined the threshold for abnormal loss of function, which corresponds to 55% wild-type function. More extreme Z scores were observed for variants with a greater loss-of-function effect. We propose that the Z score for each variant can be used to inform the application and weighting of abnormal and normal functional evidence criteria (PS3 and BS3) within the American College of Medical Genetics and Genomics variant classification framework. The validity of this approach was demonstrated using a series of 18 KCNH2 missense variants detected in a childhood onset LQTS cohort, where the level of function assessed using our assay correlated to the Schwartz score (a scoring system used to quantify the probability of a clinical diagnosis of LQTS) and the length of the corrected QT (QTc) interval.

Comparison of literature mining tools for variant classification: Through the lens of 50 RYR1 variants

PURPOSE

The American College of Medical Genetics and Genomics and the Association for Molecular Pathology have outlined a schema that allows for systematic classification of variant pathogenicity. Although gnomAD is generally accepted as a reliable source of population frequency data and ClinGen has provided guidance on the utility of specific bioinformatic predictors, there is no consensus source for identifying publications relevant to a variant. Multiple tools are available to aid in the identification of relevant variant literature, including manually curated databases and literature search engines. We set out to determine the utility of 4 literature mining tools used for ascertainment to inform the discussion of the use of these tools.

METHODS

Four literature mining tools including the Human Gene Mutation Database, Mastermind, ClinVar, and LitVar 2.0 were used to identify relevant variant literature for 50 RYR1 variants. Sensitivity and precision were determined for each tool.

RESULTS

Sensitivity among the 4 tools ranged from 0.332 to 0.687. Precision ranged from 0.389 to 0.906. No single tool retrieved all relevant publications.

CONCLUSION

At the current time, the use of multiple tools is necessary to completely identify the literature relevant to curate a variant.

Mutational scanning of CRX classifies clinical variants and reveals biochemical properties of the transcriptional effector domain

Cone-Rod Homeobox, encoded by CRX, is a transcription factor (TF) essential for the terminal differentiation and maintenance of mammalian photoreceptors. Structurally, CRX comprises an ordered DNA-binding homeodomain and an intrinsically disordered transcriptional effector domain. Although a handful of human variants in CRX have been shown to cause several different degenerative retinopathies with varying cone and rod predominance, as with most human disease genes the vast majority of observed CRX genetic variants are uncharacterized variants of uncertain significance (VUS). We performed a deep mutational scan (DMS) of nearly all possible single amino acid substitution variants in CRX, using an engineered cell-based transcriptional reporter assay. We measured the ability of each CRX missense variant to transactivate a synthetic fluorescent reporter construct in a pooled fluorescence-activated cell sorting assay and compared the activation strength of each variant to that of wild-type CRX to compute an activity score, identifying thousands of variants with altered transcriptional activity. We calculated a statistical confidence for each activity score derived from multiple independent measurements of each variant marked by unique sequence barcodes, curating a high-confidence list of nearly 2,000 variants with significantly altered transcriptional activity compared to wild-type CRX. We evaluated the performance of the DMS assay as a clinical variant classification tool using gold-standard classified human variants from ClinVar, and determined that activity scores could be used to identify pathogenic variants with high specificity. That this performance could be achieved using a synthetic reporter assay in a foreign cell type, even for a highly cell type-specific TF like CRX, suggests that this approach shows promise for DMS of other TFs that function in cell types that are not easily accessible. Per-position average activity scores closely aligned to a predicted structure of the ordered homeodomain and demonstrated position-specific residue requirements. The intrinsically disordered transcriptional effector domain, by contrast, displayed a qualitatively different pattern of substitution effects, following compositional constraints without specific residue position requirements in the peptide chain. The observed compositional constraints of the effector domain were consistent with the acidic exposure model of transcriptional activation. Together, the results of the CRX DMS identify molecular features of the CRX effector domain and demonstrate clinical utility for variant classification.

Whole genome sequencing in paediatric channelopathy and cardiomyopathy

Background

Precision medicine in paediatric cardiac channelopathy and cardiomyopathy has a rapid advancement over the past years. Compared to conventional gene panel and exome-based testing, whole genome sequencing (WGS) offers additional coverage at the promoter, intronic regions and the mitochondrial genome. However, the data on use of WGS to evaluate the genetic cause of these cardiovascular conditions in children and adolescents are limited.

Methods

In a tertiary paediatric cardiology center, we recruited all patients diagnosed with cardiac channelopathy and cardiomyopathy between the ages of 0 and 18 years old, who had negative genetic findings with prior gene panel or exome-based testing. After genetic counselling, blood samples were collected from the subjects and both their parents for WGS analysis.

Results

A total of 31 patients (11 cardiac channelopathy and 20 cardiomyopathy) were recruited. Four intronic splice-site variants were identified in three cardiomyopathy patients, which were not identified in previous whole exome sequencing. These included a pathogenic variant in TAFAZZIN:c.284+5G>A (Barth syndrome), a variant of unknown significance (VUS) in MYBPC3:c.1224-80G>A and 2 compound heterozygous LP variants in LZTR1 (LZTR1:c.1943-256C>T and LZTR1:c1261-3C>G) in a patient with clinical features of RASopathy. There was an additional diagnostic yield of 1.94% using WGS for identification of intronic variants, on top of conventional gene testing.

Conclusion

WGS plays a role in identifying additional intronic splice-site variants in paediatric patients with isolated cardiomyopathy. With the demonstrated low extra yield of WGS albeit its ability to provide potential clinically important information, WGS should be considered in selected paediatric cases of cardiac channelopathy and cardiomyopathy in a cost-effective manner.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

Functional and phenotypic consequences of an unusual inversion in MSH2

Lynch syndrome is an autosomal dominant disorder that usually results from a pathogenic germline variant in one of four genes (MSH2, MSH6, MLH1, PMS2) involved in DNA mismatch repair. Carriers of such variants are at risk of developing numerous cancers during adulthood. Here we report on a family suspected of having Lynch syndrome due to a history of endometrial adenocarcinoma, ovarian clear cell carcinoma, and adenocarcinoma of the duodenum in whom we identified a germline 29 nucleotide in-frame inversion in exon 3 of MSH2. We further show that this variant is almost completely absent at the protein level, and that the associated cancers have complete loss of MSH2 and MSH6 expression by immunohistochemistry. Functional investigation of this inversion in a laboratory setting revealed a resultant abnormal protein function. Thus, we have identified an unusual, small germline inversion in a mismatch repair gene that does not lead to a premature stop codon yet appears likely to be causal for the observed cancers.

Improved classification and pathogenicity assessment by comprehensive functional studies in a large data set of KCNQ2 variants

AIMS

The paucity of functional annotations on hundreds of KCNQ2 variants impedes the diagnosis and treatment of KCNQ2-related disorders. The aims of this work were to determine the functional properties of 331 clinical KCNQ2 variants, interpreted the pathogenicity of 331 variants using functional data，and explored the association between homomeric channel functions and phenotypes.

MAIN METHODS

We collected 145 KCNQ2 variants from 232 epilepsy patients and 186 KCNQ2 missense variants from the ClinVar database. Whole-cell patch-clamp recording was used to classify the function of 331 variants. Subsequently, we proposed 24 criteria for the pathogenicity interpretation of KCNQ2 variants and used them to assess pathogenicity of 331 variants. Finally, we analyzed the clinical phenotypes of patients carrying these variants, and explored the correlations between functional mechanisms and phenotypes.

KEY FINDINGS

In the homozygous state, 287 were classified as loss-of-function and 14 as gain-of-function. In the more clinically relative heterozygous state, 200 variants exhibited functional impairment, 121 of which showed dominant-negative effects on wild-type KCNQ2 subunits. After introducing functional data as strong-level evidence to interpret pathogenicity, over half of variants (169/331) were reclassified and 254 were classified as pathogenic/likely pathogenic. Moreover, dominant-negative effect and haploinsufficiency were identified as primary mechanisms in DEE/ID and SeLNE, respectively. The degree of impairment of channel function correlated with the phenotype severity.

SIGNIFICANCE

Our study reveals the possible cause of KCNQ2-related disorders at the molecular level, provides compelling evidence for clinical classification of KCNQ2 variants, and expands the knowledge of correlations between functional mechanisms and phenotypes.

Case report: Rare oral manifestations in Cowden syndrome with PTEN mutation

Background

Cowden syndrome (CS) is a rare genetic disorder associated with PTEN gene mutations. It is characterized by macrocephaly, specific mucocutaneous features, and a predisposition to benign and malignant tumors. Cases of CS primarily presenting with oral clinical manifestations are relatively uncommon.

Methods/Results

We report the case of a 41-year-old male proband who presented with bilateral commissural and lingual externally projecting symmetric lesions for over two years. The proband also exhibited other features, including macrocephaly, communication difficulties, and obesity. Similar oral clinical manifestations were observed in family members. Whole exome sequencing analysis revealed PTEN gene mutations associated with CS in both the proband and his younger brother. This case serves as a reminder to be aware of the diverse presentations of CS in oral clinical practice and highlights the importance of genetic testing for guiding diagnosis and treatment.

Conclusion

There are few reported cases of CS primarily presenting with oral lesions. This finding contributes to further understanding of certain aspects of the pathogenesis of CS and enhances awareness of CS cases primarily exhibiting oral clinical manifestations.