Prognostic Value of Multiplexed Assays of Variant Effect and Automated Patch-clamping for KCNH2-LQTS Risk Stratification.

Background: Long QT syndrome (LQTS) is a lethal arrhythmia condition, frequently caused by rare loss-of-function variants in the cardiac potassium channel encoded by KCNH2. Variant-based risk stratification is complicated by heterogenous clinical data, incomplete penetrance, and low-throughput functional data. Objective: To test the utility of variant-specific features, including high-throughput functional data, to predict cardiac events among KCNH2 variant heterozygotes. Methods: We quantified cell-surface trafficking of 18,323 variants in KCNH2 and recorded potassium current densities for 506 KCNH2 variants. Next, we deeply phenotyped 1150 KCNH2 missense variant patients, including ECG features, cardiac event history (528 total cardiac events), and mortality. We then assessed variant functional, in silico, structural, and LQTS penetrance data to stratify event-free survival for cardiac events in the study cohort. Results: Variant-specific current density (HR 0.28 [0.13-0.60]) and estimates of LQTS penetrance incorporating MAVE data (HR 3.16 [1.59-6.27]) were independently predictive of severe cardiac events when controlling for patient-specific features. Risk prediction models incorporating these data significantly improved prediction of 20 year cardiac events (AUC 0.79 [0.75-0.82]) over patient-only covariates (QTc and sex) (AUC 0.73 [0.70-0.77]). Conclusion: We show that high-throughput functional data, and other variant-specific features, meaningfully contribute to both diagnosis and prognosis of a clinically actionable monogenic disease.

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.

Multicenter clinical and functional evidence reclassifies a recurrent noncanonical filamin C splice-altering variant.

BACKGROUND: Truncating variants in filamin C (FLNC) can cause arrhythmogenic cardiomyopathy (ACM) through haploinsufficiency. Noncanonical splice-altering variants may contribute to this phenotype. OBJECTIVE: The purpose of this study was to investigate the clinical and functional consequences of a recurrent FLNC intronic variant of uncertain significance (VUS), c.970-4A>G. METHODS: Clinical data in 9 variant heterozygotes from 4 kindreds were obtained from 5 tertiary health care centers. We used in silico predictors and functional studies with peripheral blood and patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Isolated RNA was studied by reverse transcription polymerase chain reaction. iPSC-CMs were further characterized at baseline and after nonsense-mediated decay (NMD) inhibition, using quantitative polymerase chain reaction (qPCR), RNA-sequencing, and cellular electrophysiology. American College of Medical Genetics and Genomics (ACMG) criteria were used to adjudicate variant pathogenicity. RESULTS: Variant heterozygotes displayed a spectrum of disease phenotypes, spanning from mild ventricular dysfunction with palpitations to severe ventricular arrhythmias requiring device shocks or progressive cardiomyopathy requiring heart transplantation. Consistent with in silico predictors, the c.970-4A>G FLNC variant activated a cryptic splice acceptor site, introducing a 3-bp insertion containing a premature termination codon. NMD inhibition upregulated aberrantly spliced transcripts by qPCR and RNA-sequencing. Patch clamp studies revealed irregular spontaneous action potentials, increased action potential duration, and increased sodium late current in proband-derived iPSC-CMs. These findings fulfilled multiple ACMG criteria for pathogenicity. CONCLUSION: Clinical, in silico, and functional evidence support the prediction that the intronic c.970-4A>G VUS disrupts splicing and drives ACM, enabling reclassification from VUS to pathogenic.

Continuous Bayesian variant interpretation accounts for incomplete penetrance among Mendelian cardiac channelopathies.

PURPOSE: The congenital Long QT Syndrome (LQTS) and Brugada Syndrome (BrS) are Mendelian autosomal dominant diseases that frequently precipitate fatal cardiac arrhythmias. Incomplete penetrance is a barrier to clinical management of heterozygotes harboring variants in the major implicated disease genes KCNQ1, KCNH2, and SCN5A. We apply and evaluate a Bayesian penetrance estimation strategy that accounts for this phenomenon. METHODS: We generated Bayesian penetrance models for KCNQ1-LQT1 and SCN5A-LQT3 using variant-specific features and clinical data from the literature, international arrhythmia genetic centers, and population controls. We analyzed the distribution of posterior penetrance estimates across 4 genotype-phenotype relationships and compared continuous estimates with ClinVar annotations. Posterior estimates were mapped onto protein structure. RESULTS: Bayesian penetrance estimates of KCNQ1-LQT1 and SCN5A-LQT3 are empirically equivalent to 10 and 5 clinically phenotype heterozygotes, respectively. Posterior penetrance estimates were bimodal for KCNQ1-LQT1 and KCNH2-LQT2, with a higher fraction of missense variants with high penetrance among KCNQ1 variants. There was a wide distribution of variant penetrance estimates among identical ClinVar categories. Structural mapping revealed heterogeneity among "hot spot" regions and featured high penetrance estimates for KCNQ1 variants in contact with calmodulin and the S6 domain. CONCLUSIONS: Bayesian penetrance estimates provide a continuous framework for variant interpretation.

A massively parallel assay accurately discriminates between functionally normal and abnormal variants in a hotspot domain of KCNH2.

Many genes, including KCNH2, contain "hotspot" domains associated with a high density of variants associated with disease. This has led to the suggestion that variant location can be used as evidence supporting classification of clinical variants. However, it is not known what proportion of all potential variants in hotspot domains cause loss of function. Here, we have used a massively parallel trafficking assay to characterize all single-nucleotide variants in exon 2 of KCNH2, a known hotspot for variants that cause long QT syndrome type 2 and an increased risk of sudden cardiac death. Forty-two percent of KCNH2 exon 2 variants caused at least 50% reduction in protein trafficking, and 65% of these trafficking-defective variants exerted a dominant-negative effect when co-expressed with a WT KCNH2 allele as assessed using a calibrated patch-clamp electrophysiology assay. The massively parallel trafficking assay was more accurate (AUC of 0.94) than bioinformatic prediction tools (REVEL and CardioBoost, AUC of 0.81) in discriminating between functionally normal and abnormal variants. Interestingly, over half of variants in exon 2 were found to be functionally normal, suggesting a nuanced interpretation of variants in this "hotspot" domain is necessary. Our massively parallel trafficking assay can provide this information prospectively.

A calibrated functional patch-clamp assay to enhance clinical variant interpretation in KCNH2-related long QT syndrome.

Modern sequencing technologies have revolutionized our detection of gene variants. However, in most genes, including KCNH2, the majority of missense variants are currently classified as variants of uncertain significance (VUSs). The aim of this study was to investigate the utility of an automated patch-clamp assay for aiding clinical variant classification in KCNH2. The assay was designed according to recommendations proposed by the Clinical Genome Sequence Variant Interpretation Working Group. Thirty-one variants (17 pathogenic/likely pathogenic, 14 benign/likely benign) were classified internally as variant controls. They were heterozygously expressed in Flp-In HEK293 cells for assessing the effects of variants on current density and channel gating in order to determine the sensitivity and specificity of the assay. All 17 pathogenic variant controls had reduced current density, and 13 of 14 benign variant controls had normal current density, which enabled determination of normal and abnormal ranges for applying evidence of moderate or supporting strength for VUS reclassification. Inclusion of functional assay evidence enabled us to reclassify 6 out of 44 KCNH2 VUSs as likely pathogenic. The high-throughput patch-clamp assay can provide moderate-strength evidence for clinical interpretation of clinical KCNH2 variants and demonstrates the value of developing automated patch-clamp assays for functional characterization of ion channel gene variants.

Veratridine Can Bind to a Site at the Mouth of the Channel Pore at Human Cardiac Sodium Channel NaV1.5.

The cardiac sodium ion channel (NaV1.5) is a protein with four domains (DI-DIV), each with six transmembrane segments. Its opening and subsequent inactivation results in the brief rapid influx of Na+ ions resulting in the depolarization of cardiomyocytes. The neurotoxin veratridine (VTD) inhibits NaV1.5 inactivation resulting in longer channel opening times, and potentially fatal action potential prolongation. VTD is predicted to bind at the channel pore, but alternative binding sites have not been ruled out. To determine the binding site of VTD on NaV1.5, we perform docking calculations and high-throughput electrophysiology experiments in the present study. The docking calculations identified two distinct binding regions. The first site was in the pore, close to the binding site of NaV1.4 and NaV1.5 blocking drugs in experimental structures. The second site was at the "mouth" of the pore at the cytosolic side, partly solvent-exposed. Mutations at this site (L409, E417, and I1466) had large effects on VTD binding, while residues deeper in the pore had no effect, consistent with VTD binding at the mouth site. Overall, our results suggest a VTD binding site close to the cytoplasmic mouth of the channel pore. Binding at this alternative site might indicate an allosteric inactivation mechanism for VTD at NaV1.5.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study.

BACKGROUND: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. METHODS: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. RESULTS: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. CONCLUSIONS: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier; NCT03394859.

Estimating the Posttest Probability of Long QT Syndrome Diagnosis for Rare KCNH2 Variants.

BACKGROUND: The proliferation of genetic profiling has revealed many associations between genetic variations and disease. However, large-scale phenotyping efforts in largely healthy populations, coupled with DNA sequencing, suggest variants currently annotated as pathogenic are more common in healthy populations than previously thought. In addition, novel and rare variants are frequently observed in genes associated with disease both in healthy individuals and those under suspicion of disease. This raises the question of whether these variants can be useful predictors of disease. To answer this question, we assessed the degree to which the presence of a variant in the cardiac potassium channel gene KCNH2 was diagnostically predictive for the autosomal dominant long QT syndrome. METHODS: We estimated the probability of a long QT diagnosis given the presence of each KCNH2 variant using Bayesian methods that incorporated variant features such as changes in variant function, protein structure, and in silico predictions. We call this estimate the posttest probability of disease. Our method was applied to over 4000 individuals heterozygous for 871 missense or in-frame insertion/deletion variants in KCNH2 and validated against a separate international cohort of 933 individuals heterozygous for 266 missense or in-frame insertion/deletion variants. RESULTS: Our method was well-calibrated for the observed fraction of heterozygotes diagnosed with long QT syndrome. Heuristically, we found that the innate diagnostic information one learns about a variant from 3-dimensional variant location, in vitro functional data, and in silico predictors is equivalent to the diagnostic information one learns about that same variant by clinically phenotyping 10 heterozygotes. Most importantly, these data can be obtained in the absence of any clinical observations. CONCLUSIONS: We show how variant-specific features can inform a prior probability of disease for rare variants even in the absence of clinically phenotyped heterozygotes.

High-throughput discovery of trafficking-deficient variants in the cardiac potassium channel KV11.1.

BACKGROUND: KCHN2 encodes the KV11.1 potassium channel responsible for IKr, a major repolarization current during the cardiomyocyte action potential. Variants in KCNH2 that lead to decreased IKr have been associated with long QT syndrome type 2 (LQT2). The mechanism of LQT2 is most often induced loss of KV11.1 trafficking to the cell surface. Accurately discriminating between variants with normal and abnormal trafficking would aid in understanding the deleterious nature of these variants; however, the volume of reported nonsynonymous KCNH2 variants precludes the use of conventional methods for functional study. OBJECTIVE: The purpose of this study was to report a high-throughput, multiplexed screening method for KCNH2 genetic variants capable of measuring the cell surface abundance of hundreds of missense variants in the resulting KV11.1 channel. METHODS: We developed a method to quantitate KV11.1 variant trafficking on a pilot region of 11 residues in the S5 helix. RESULTS: We generated trafficking scores for 220 of 231 missense variants in the pilot region. For 5 of 5 variants, high-throughput trafficking scores validated when tested in single variant flow cytometry and confocal microscopy experiments. We further explored these results with planar patch electrophysiology and found that loss-of-trafficking variants do not produce IKr. Conversely, but expectedly, some variants that traffic normally were still functionally compromised. CONCLUSION: We describe a new method for detecting KV11.1 trafficking-deficient variants in a multiplexed assay. This new method accurately generated trafficking data for variants in KV11.1 and is extendable both to all residues in KV11.1 and to other cell surface proteins.

A Bayesian method to estimate variant-induced disease penetrance.

A major challenge emerging in genomic medicine is how to assess best disease risk from rare or novel variants found in disease-related genes. The expanding volume of data generated by very large phenotyping efforts coupled to DNA sequence data presents an opportunity to reinterpret genetic liability of disease risk. Here we propose a framework to estimate the probability of disease given the presence of a genetic variant conditioned on features of that variant. We refer to this as the penetrance, the fraction of all variant heterozygotes that will present with disease. We demonstrate this methodology using a well-established disease-gene pair, the cardiac sodium channel gene SCN5A and the heart arrhythmia Brugada syndrome. From a review of 756 publications, we developed a pattern mixture algorithm, based on a Bayesian Beta-Binomial model, to generate SCN5A penetrance probabilities for the Brugada syndrome conditioned on variant-specific attributes. These probabilities are determined from variant-specific features (e.g. function, structural context, and sequence conservation) and from observations of affected and unaffected heterozygotes. Variant functional perturbation and structural context prove most predictive of Brugada syndrome penetrance.

High-Throughput Reclassification of SCN5A Variants.

Partial or complete loss-of-function variants in SCN5A are the most common genetic cause of the arrhythmia disorder Brugada syndrome (BrS1). However, the pathogenicity of SCN5A variants is often unknown or disputed; 80% of the 1,390 SCN5A missense variants observed in at least one individual to date are variants of uncertain significance (VUSs). The designation of VUS is a barrier to the use of sequence data in clinical care. We selected 83 variants: 10 previously studied control variants, 10 suspected benign variants, and 63 suspected Brugada syndrome-associated variants, selected on the basis of their frequency in the general population and in individuals with Brugada syndrome. We used high-throughput automated patch clamping to study the function of the 83 variants, with the goal of reclassifying variants with functional data. The ten previously studied controls had functional properties concordant with published manual patch clamp data. All 10 suspected benign variants had wild-type-like function. 22 suspected BrS variants had loss of channel function (<10% normalized peak current) and 22 variants had partial loss of function (10%-50% normalized peak current). The previously unstudied variants were initially classified as likely benign (n = 2), likely pathogenic (n = 10), or VUSs (n = 61). After the patch clamp studies, 16 variants were benign/likely benign, 45 were pathogenic/likely pathogenic, and only 12 were still VUSs. Structural modeling identified likely mechanisms for loss of function including altered thermostability and disruptions to alpha helices, disulfide bonds, or the permeation pore. High-throughput patch clamping enabled reclassification of the majority of tested VUSs in SCN5A.

Deep Mutational Scan of an SCN5A Voltage Sensor.

BACKGROUND: Variants in ion channel genes have classically been studied in low throughput by patch clamping. Deep mutational scanning is a complementary approach that can simultaneously assess function of thousands of variants. METHODS: We have developed and validated a method to perform a deep mutational scan of variants in SCN5A, which encodes the major voltage-gated sodium channel in the heart. We created a library of nearly all possible variants in a 36 base region of SCN5A in the S4 voltage sensor of domain IV and stably integrated the library into HEK293T cells. RESULTS: In preliminary experiments, challenge with 3 drugs (veratridine, brevetoxin, and ouabain) could discriminate wild-type channels from gain- and loss-of-function pathogenic variants. High-throughput sequencing of the pre- and postdrug challenge pools was used to count the prevalence of each variant and identify variants with abnormal function. The deep mutational scan scores identified 40 putative gain-of-function and 33 putative loss-of-function variants. For 8 of 9 variants, patch clamping data were consistent with the scores. CONCLUSIONS: These experiments demonstrate the accuracy of a high-throughput in vitro scan of SCN5A variant function, which can be used to identify deleterious variants in SCN5A and other ion channel genes.

Patient-independent human induced pluripotent stem cell model: A new tool for rapid determination of genetic variant pathogenicity in long QT syndrome.

BACKGROUND: Commercial genetic testing for long QT syndrome (LQTS) has rapidly expanded, but the inability to accurately predict whether a rare variant is pathogenic has limited its clinical benefit. Novel missense variants are routinely reported as variant of unknown significance (VUS) and cannot be used to screen family members at risk for sudden cardiac death. Better approaches to determine the pathogenicity of VUS are needed. OBJECTIVE: The purpose of this study was to rapidly determine the pathogenicity of a CACNA1C variant reported by commercial genetic testing as a VUS using a patient-independent human induced pluripotent stem cell (hiPSC) model. METHODS: Using CRISPR/Cas9 genome editing, CACNA1C-p.N639T was introduced into a previously established hiPSC from an unrelated healthy volunteer, thereby generating a patient-independent hiPSC model. Three independent heterozygous N639T hiPSC lines were generated and differentiated into cardiomyocytes (CM). Electrophysiological properties of N639T hiPSC-CM were compared to those of isogenic and population control hiPSC-CM by measuring the extracellular field potential (EFP) of 96-well hiPSC-CM monolayers and by patch clamp. RESULTS: Significant EFP prolongation was observed only in optically stimulated but not in spontaneously beating N639T hiPSC-CM. Patch-clamp studies revealed that N639T prolonged the ventricular action potential by slowing voltage-dependent inactivation of CaV1.2 currents. Heterologous expression studies confirmed the effect of N639T on CaV1.2 inactivation. CONCLUSION: The patient-independent hiPSC model enabled rapid generation of functional data to support reclassification of a CACNA1C VUS to likely pathogenic, thereby establishing a novel LQTS type 8 mutation. Furthermore, our results indicate the importance of controlling beating rates to evaluate the functional significance of LQTS VUS in high-throughput hiPSC-CM assays.

Protein structure aids predicting functional perturbation of missense variants in SCN5A and KCNQ1.

Rare variants in the cardiac potassium channel KV7.1 (KCNQ1) and sodium channel NaV1.5 (SCN5A) are implicated in genetic disorders of heart rhythm, including congenital long QT and Brugada syndromes (LQTS, BrS), but also occur in reference populations. We previously reported two sets of NaV1.5 (n = 356) and KV7.1 (n = 144) variants with in vitro characterized channel currents gathered from the literature. Here we investigated the ability to predict commonly reported NaV1.5 and KV7.1 variant functional perturbations by leveraging diverse features including variant classifiers PROVEAN, PolyPhen-2, and SIFT; evolutionary rate and BLAST position specific scoring matrices (PSSM); and structure-based features including "functional densities" which is a measure of the density of pathogenic variants near the residue of interest. Structure-based functional densities were the most significant features for predicting NaV1.5 peak current (adj. R2 = 0.27) and KV7.1 + KCNE1 half-maximal voltage of activation (adj. R2 = 0.29). Additionally, use of structure-based functional density values improves loss-of-function classification of SCN5A variants with an ROC-AUC of 0.78 compared with other predictive classifiers (AUC = 0.69; two-sided DeLong test p = .01). These results suggest structural data can inform predictions of the effect of uncharacterized SCN5A and KCNQ1 variants to provide a deeper understanding of their burden on carriers.

Multiple mechanisms underlie increased cardiac late sodium current.

BACKGROUND: We recently reported a quantitative relationship between the degree of functional perturbation reported in the literature for 356 variants in the cardiac sodium channel gene SCN5A and the penetrance of resulting arrhythmia phenotypes. In the course of that work, we identified multiple SCN5A variants, including R1193Q, that are common in populations but are reported in human embryonic kidney (HEK) cells to generate large late sodium current (INa-L). OBJECTIVE: The purpose of this study was to compare the functional properties of R1193Q with those of the well-studied type 3 long QT syndrome mutation ΔKPQ. METHODS: We compared functional properties of SCN5A R1193Q with those of ΔKPQ in Chinese hamster ovary (CHO) cells at baseline and after exposure to intracellular phosphatidylinositol (3,4,5)-trisphosphate (PIP3), which inhibits INa-L generated by decreased Phosphoinositide 3-kinase (PI3K) activity. We also used CRISPR/Cas9 editing to generate R1193Q in human-induced pluripotent stem cells differentiated to cardiomyocytes (hiPSC-CMs). RESULTS: Both R1193Q and ΔKPQ generated robust INa-L in CHO cells. PIP3 abrogated the late current phenotype in R1193Q cells but had no effect on ΔKPQ. Homozygous R1193Q hiPSC-CMs displayed increased INa-L and long action potentials with frequent triggered beats, which were reversed with the addition of PIP3. CONCLUSION: The consistency between the late current produced in HEK cells, CHO cells, and hiPSC-CMs suggests that the late current is a feature of the SCN5A R1193Q variant in human cardiomyocytes but that the mechanism by which the late current is produced is distinct and indirect, as compared with the more highly penetrant ΔKPQ. These data suggest that observing a late current in an in vitro setting does not necessarily translate to highly pathogenic type 3 long QT syndrome phenotype but depends on the underlying mechanism.

SCN5A (NaV1.5) Variant Functional Perturbation and Clinical Presentation: Variants of a Certain Significance.

BACKGROUND: Accurately predicting the impact of rare nonsynonymous variants on disease risk is an important goal in precision medicine. Variants in the cardiac sodium channel SCN5A (protein NaV1.5; voltage-dependent cardiac Na+ channel) are associated with multiple arrhythmia disorders, including Brugada syndrome and long QT syndrome. Rare SCN5A variants also occur in ≈1% of unaffected individuals. We hypothesized that in vitro electrophysiological functional parameters explain a statistically significant portion of the variability in disease penetrance. METHODS: From a comprehensive literature review, we quantified the number of carriers presenting with and without disease for 1712 reported SCN5A variants. For 356 variants, data were also available for 5 NaV1.5 electrophysiological parameters: peak current, late/persistent current, steady-state V1/2 of activation and inactivation, and recovery from inactivation. RESULTS: We found that peak and late current significantly associate with Brugada syndrome (P<0.001; ρ=-0.44; Spearman rank test) and long QT syndrome disease penetrance (P<0.001; ρ=0.37). Steady-state V1/2 activation and recovery from inactivation associate significantly with Brugada syndrome and long QT syndrome penetrance, respectively. Continuous estimates of disease penetrance align with the current American College of Medical Genetics classification paradigm. CONCLUSIONS: NaV1.5 in vitro electrophysiological parameters are correlated with Brugada syndrome and long QT syndrome disease risk. Our data emphasize the value of in vitro electrophysiological characterization and incorporating counts of affected and unaffected carriers to aid variant classification. This quantitative analysis of the electrophysiological literature should aid the interpretation of NaV1.5 variant electrophysiological abnormalities and help improve NaV1.5 variant classification.

A Mechanism of Calmodulin Modulation of the Human Cardiac Sodium Channel.

The function of the human cardiac sodium channel (NaV1.5) is modulated by the Ca2+ sensor calmodulin (CaM), but the underlying mechanism(s) are controversial and poorly defined. CaM has been reported to bind in a Ca2+-dependent manner to two sites in the intracellular loop that is critical for inactivation of NaV1.5 (inactivation gate [IG]). The affinity of CaM for the complete IG was significantly stronger than that of fragments that lacked both complete binding sites. Structural analysis by nuclear magnetic resonance, crystallographic, and scattering approaches revealed that CaM simultaneously engages both IG sites using an extended configuration. Patch-clamp recordings for wild-type and mutant channels with an impaired CaM-IG interaction revealed CaM binding to the IG promotes recovery from inactivation while impeding the kinetics of inactivation. Models of full-length NaV1.5 suggest that CaM binding to the IG directly modulates channel function by destabilizing the inactivated state, which would promote resetting of the IG after channels close.

Exploiting ion channel structure to assess rare variant pathogenicity.

BACKGROUND: A 27-year-old woman was seen for long QT syndrome. She was found to be a carrier of 2 variants, KCNQ1 Val162Met and KCNH2 Ser55Leu, and both were classified as "pathogenic" by a diagnostic laboratory, in part because of sequence proximity to other known pathogenic variants. OBJECTIVE: The purpose of this study was to assess the relationship between both the KCNQ1 and KCNH2 variants and clinical significance using protein structure, in vitro functional assays, and familial segregation. METHODS: We used co-segregation analysis of family, patch clamp in vitro electrophysiology, and structural analysis using recently released cryo-electron microscopy structures of both channels. RESULTS: The structural analysis indicates that KCNQ1 Val162Met is oriented away from functionally important regions while Ser55Leu is positioned at domains critical for KCNH2 fast inactivation. Clinical phenotyping and electrophysiology study further support the conclusion that KCNH2 Ser55Leu is correctly classified as pathogenic but KCNQ1 Val162Met is benign. CONCLUSION: Proximity in sequence space does not always translate accurately to proximity in 3-dimensional space. Emerging structural methods will add value to pathogenicity prediction.

Predicting the Functional Impact of KCNQ1 Variants of Unknown Significance.

BACKGROUND: An emerging standard-of-care for long-QT syndrome uses clinical genetic testing to identify genetic variants of the KCNQ1 potassium channel. However, interpreting results from genetic testing is confounded by the presence of variants of unknown significance for which there is inadequate evidence of pathogenicity. METHODS AND RESULTS: In this study, we curated from the literature a high-quality set of 107 functionally characterized KCNQ1 variants. Based on this data set, we completed a detailed quantitative analysis on the sequence conservation patterns of subdomains of KCNQ1 and the distribution of pathogenic variants therein. We found that conserved subdomains generally are critical for channel function and are enriched with dysfunctional variants. Using this experimentally validated data set, we trained a neural network, designated Q1VarPred, specifically for predicting the functional impact of KCNQ1 variants of unknown significance. The estimated predictive performance of Q1VarPred in terms of Matthew's correlation coefficient and area under the receiver operating characteristic curve were 0.581 and 0.884, respectively, superior to the performance of 8 previous methods tested in parallel. Q1VarPred is publicly available as a web server at http://meilerlab.org/q1varpred. CONCLUSIONS: Although a plethora of tools are available for making pathogenicity predictions over a genome-wide scale, previous tools fail to perform in a robust manner when applied to KCNQ1. The contrasting and favorable results for Q1VarPred suggest a promising approach, where a machine-learning algorithm is tailored to a specific protein target and trained with a functionally validated data set to calibrate informatics tools.