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

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.

Variable Penetrance and Expressivity of a Rare Pore Loss-of-Function Mutation (p.L889V) of Nav1.5 Channels in Three Spanish Families

A novel rare mutation in the pore region of Nav1.5 channels (p.L889V) has been found in three unrelated Spanish families that produces quite diverse phenotypic manifestations (Brugada syndrome, conduction disease, dilated cardiomyopathy, sinus node dysfunction, etc.) with variable penetrance among families. We clinically characterized the carriers and recorded the Na+ current (INa) generated by p.L889V and native (WT) Nav1.5 channels, alone or in combination, to obtain further insight into the genotypic-phenotypic relationships in patients carrying SCN5A mutations and in the molecular determinants of the Nav1.5 channel function. The variant produced a strong dominant negative effect (DNE) since the peak INa generated by p.L889V channels expressed in Chinese hamster ovary cells, either alone (-69.4 ± 9.0 pA/pF) or in combination with WT (-62.2 ± 14.6 pA/pF), was significantly (n ≥ 17, p < 0.05) reduced compared to that generated by WT channels alone (-199.1 ± 44.1 pA/pF). The mutation shifted the voltage dependence of channel activation and inactivation to depolarized potentials, did not modify the density of the late component of INa, slightly decreased the peak window current, accelerated the recovery from fast and slow inactivation, and slowed the induction kinetics of slow inactivation, decreasing the fraction of channels entering this inactivated state. The membrane expression of p.L889V channels was low, and in silico molecular experiments demonstrated profound alterations in the disposition of the pore region of the mutated channels. Despite the mutation producing a marked DNE and reduction in the INa and being located in a critical domain of the channel, its penetrance and expressivity are quite variable among the carriers. Our results reinforce the argument that the incomplete penetrance and phenotypic variability of SCN5A loss-of-function mutations are the result of a combination of multiple factors, making it difficult to predict their expressivity in the carriers despite the combination of clinical, genetic, and functional studies.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants

Background

Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies.

Methods

We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs.

Results

ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing.

Conclusions

ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

The Action Potential Clamp Technique as a Tool for Risk Stratification of Sinus Bradycardia Due to Loss-of-Function Mutations in HCN4: An In Silico Exploration Based on In Vitro and In Vivo Data

These days, in vitro functional analysis of gene variants is becoming increasingly important for risk stratification of cardiac ion channelopathies. So far, such risk stratification has been applied to SCN5A, KCNQ1, and KCNH2 gene variants associated with Brugada syndrome and long QT syndrome types 1 and 2, respectively, but risk stratification of HCN4 gene variants related to sick sinus syndrome has not yet been performed. HCN4 is the gene responsible for the hyperpolarization-activated 'funny' current If, which is an important modulator of the spontaneous diastolic depolarization underlying the sinus node pacemaker activity. In the present study, we carried out a risk classification assay on those loss-of-function mutations in HCN4 for which in vivo as well as in vitro data have been published. We used the in vitro data to compute the charge carried by If (Qf) during the diastolic depolarization phase of a prerecorded human sinus node action potential waveform and assessed the extent to which this Qf predicts (1) the beating rate of the comprehensive Fabbri-Severi model of a human sinus node cell with mutation-induced changes in If and (2) the heart rate observed in patients carrying the associated mutation in HCN4. The beating rate of the model cell showed a very strong correlation with Qf from the simulated action potential clamp experiments (R2 = 0.95 under vagal tone). The clinically observed minimum or resting heart rates showed a strong correlation with Qf (R2 = 0.73 and R2 = 0.71, respectively). While a translational perspective remains to be seen, we conclude that action potential clamp on transfected cells, without the need for further voltage clamp experiments and data analysis to determine individual biophysical parameters of If, is a promising tool for risk stratification of sinus bradycardia due to loss-of-function mutations in HCN4. In combination with an If blocker, this tool may also prove useful when applied to human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) obtained from mutation carriers and non-carriers.

Functional Characterisation of the Rare SCN5A p.E1225K Variant, Segregating in a Brugada Syndrome Familial Case, in Human Cardiomyocytes from Pluripotent Stem Cells

Brugada syndrome (BrS) is an inherited autosomal dominant cardiac channelopathy. Pathogenic rare mutations in the SCN5A gene, encoding the alpha-subunit of the voltage-dependent cardiac Na⁺ channel protein (Nav1.5), are identified in 20% of BrS patients, affecting the correct function of the channel. To date, even though hundreds of SCN5A variants have been associated with BrS, the underlying pathogenic mechanisms are still unclear in most cases. Therefore, the functional characterization of the SCN5A BrS rare variants still represents a major hurdle and is fundamental to confirming their pathogenic effect. Human cardiomyocytes (CMs) differentiated from pluripotent stem cells (PSCs) have been extensively demonstrated to be reliable platforms for investigating cardiac diseases, being able to recapitulate specific traits of disease, including arrhythmic events and conduction abnormalities. Based on this, in this study, we performed a functional analysis of the BrS familial rare variant NM_198056.2:c.3673G>A (NP_932173.1:p.Glu1225Lys), which has been never functionally characterized before in a cardiac-relevant context, as the human cardiomyocyte. Using a specific lentiviral vector encoding a GFP-tagged SCN5A gene carrying the specific c.3673G>A variant and CMs differentiated from control PSCs (PSC-CMs), we demonstrated an impairment of the mutated Nav1.5, thus suggesting the pathogenicity of the rare BrS detected variant. More broadly, our work supports the application of PSC-CMs for the assessment of the pathogenicity of gene variants, the identification of which is increasing exponentially due to the advances in next-generation sequencing methods and their massive use in genetic testing.

Automated Patch-Clamp and Induced Pluripotent Stem Cell-Derived Cardiomyocytes: A Synergistic Approach in the Study of Brugada Syndrome

The development of high-throughput automated patch-clamp technology is a recent breakthrough in the field of Brugada syndrome research. Brugada syndrome is a heart disorder marked by abnormal electrocardiographic readings and an elevated risk of sudden cardiac death due to arrhythmias. Various experimental models, developed either in animals, cell lines, human tissue or computational simulation, play a crucial role in advancing our understanding of this condition, and developing effective treatments. In the perspective of the pathophysiological role of ion channels and their pharmacology, automated patch-clamp involves a robotic system that enables the simultaneous recording of electrical activity from multiple single cells at once, greatly improving the speed and efficiency of data collection. By combining this approach with the use of patient-derived cardiomyocytes, researchers are gaining a more comprehensive view of the underlying mechanisms of heart disease. This has led to the development of more effective treatments for those affected by cardiovascular conditions.

Dominant negative effects of SCN5A missense variants

PURPOSE

Up to 30% of patients with Brugada syndrome (BrS) carry loss-of-function (LoF) variants in the cardiac sodium channel gene SCN5A encoding for the protein Na_V1.5. Recent studies suggested that Na_V1.5 can dimerize, and some variants exert dominant negative effects. In this study, we sought to explore the generality of missense variant Na_V1.5 dominant negative effects and their clinical severity.

METHODS

We identified 35 LoF variants (<10% of wild type [WT] peak current) and 15 partial LoF variants (10%-50% of WT peak current) that we assessed for dominant negative effects. SCN5A variants were studied in HEK293T cells, alone or in heterozygous coexpression with WT SCN5A using automated patch clamp. To assess the clinical risk, we compared the prevalence of dominant negative vs putative haploinsufficient (frameshift, splice, or nonsense) variants in a BrS consortium and the Genome Aggregation Database population database.

RESULTS

In heterozygous expression with WT, 32 of 35 LoF and 6 of 15 partial LoF variants showed reduction to <75% of WT-alone peak current, showing a dominant negative effect. Individuals with dominant negative LoF variants had an elevated disease burden compared with the individuals with putative haploinsufficient variants (2.7-fold enrichment in BrS cases, P = .019).

CONCLUSION

Most SCN5A missense LoF variants exert a dominant negative effect. This class of variant confers an especially high burden of BrS.

The Genetics of Brugada Syndrome

Brugada syndrome is a heritable channelopathy characterized by a peculiar electrocardiogram (ECG) pattern and increased risk of cardiac arrhythmias and sudden death. The arrhythmias originate because of an imbalance between the repolarizing and depolarizing currents that modulate the cardiac action potential. Even if an overt structural cardiomyopathy is not typical of Brugada syndrome, fibrosis and structural changes in the right ventricle contribute to a conduction slowing, which ultimately facilitates ventricular arrhythmias. Currently, Mendelian autosomal dominant transmission is detected in less than 25% of all clinical confirmed cases. Although 23 genes have been associated with the condition, only SCN5A, encoding the cardiac sodium channel, is considered clinically actionable and disease causing. The limited monogenic inheritance has pointed toward new perspectives on the possible complex genetic architecture of the disease, involving polygenic inheritance and a polygenic risk score that can influence penetrance and risk stratification. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Causative Variants for Inherited Cardiac Conditions in a Southeast Asian Population Cohort

Background:

Variable penetrance and late-onset phenotypes are key challenges for classifying causal as well as incidental findings in inherited cardiac conditions. Allele frequencies of variants in ancestry-specific populations, along with clinical variant analysis and interpretation, are critical to determine their true significance.

Methods:

Here, we carefully reviewed and classified variants in genes associated with inherited cardiac conditions based on a population whole-genome sequencing cohort of 4810 Singaporeans representing Southeast Asian ancestries.

Results:

Eighty-nine (1.85%) individuals carried either pathogenic or likely pathogenic variants across 25 genes. Forty-six (51.7%) had variants in causal genes for familial hyperlipidemia, but there were also recurrent variants in SCN5A and MYBPC3 , causal genes for inherited arrhythmia and cardiomyopathy, which, despite previous reports, we determined to lack criteria for pathogenicity.

Conclusions:

Our findings highlight the incidence of disease-related variants in inherited cardiac conditions and emphasize the value of large-scale sequencing in specific ancestries. Follow-up detailed phenotyping and analysis of pedigrees are crucial because assigning pathogenicity will significantly affect clinical management for individuals and their family members.

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.

Genomic and Non-Genomic Regulatory Mechanisms of the Cardiac Sodium Channel in Cardiac Arrhythmias

Nav1.5 is the predominant cardiac sodium channel subtype, encoded by the SCN5A gene, which is involved in the initiation and conduction of action potentials throughout the heart. Along its biosynthesis process, Nav1.5 undergoes strict genomic and non-genomic regulatory and quality control steps that allow only newly synthesized channels to reach their final membrane destination and carry out their electrophysiological role. These regulatory pathways are ensured by distinct interacting proteins that accompany the nascent Nav1.5 protein along with different subcellular organelles. Defects on a large number of these pathways have a tremendous impact on Nav1.5 functionality and are thus intimately linked to cardiac arrhythmias. In the present review, we provide current state-of-the-art information on the molecular events that regulate SCN5A/Nav1.5 and the cardiac channelopathies associated with defects in these pathways.

Functionally validated SCN5A variants allow interpretation of pathogenicity and prediction of lethal events in Brugada syndrome

AIMS

The prognostic value of genetic variants for predicting lethal arrhythmic events (LAEs) in Brugada syndrome (BrS) remains controversial. We investigated whether the functional curation of SCN5A variations improves prognostic predictability.

METHODS AND RESULTS

Using a heterologous expression system and whole-cell patch clamping, we functionally characterized 22 variants of unknown significance (VUSs) among 55 SCN5A mutations previously curated using in silico prediction algorithms in the Japanese BrS registry (n = 415). According to the loss-of-function (LOF) properties, SCN5A mutation carriers (n = 60) were divided into two groups: LOF-SCN5A mutations and non-LOF SCN5A variations. Functionally proven LOF-SCN5A mutation carriers (n = 45) showed significantly severer electrocardiographic conduction abnormalities and worse prognosis associated with earlier manifestations of LAEs (7.9%/year) than in silico algorithm-predicted SCN5A carriers (5.1%/year) or all BrS probands (2.5%/year). Notably, non-LOF SCN5A variation carriers (n = 15) exhibited no LAEs during the follow-up period. Multivariate analysis demonstrated that only LOF-SCN5A mutations and a history of aborted cardiac arrest were significant predictors of LAEs. Gene-based association studies using whole-exome sequencing data on another independent SCN5A mutation-negative BrS cohort (n = 288) showed no significant enrichment of rare variants in 16 985 genes including 22 non-SCN5A BrS-associated genes as compared with controls (n = 372). Furthermore, rare variations of non-SCN5A BrS-associated genes did not affect LAE-free survival curves.

CONCLUSION

In vitro functional validation is key to classifying the pathogenicity of SCN5A VUSs and for risk stratification of genetic predictors of LAEs. Functionally proven LOF-SCN5A mutations are genetic burdens of sudden death in BrS, but evidence for other BrS-associated genes is elusive.

Brugada syndrome masked by complete left bundle branch block: A clinical and functional study of its association with the p.1449Y>H SCN5A variant

SCN5A gene variants are associated with both Brugada syndrome and conduction disturbances, sometimes expressing an overlapping phenotype. Functional consequences of SCN5A variants assessed by patch-clamp electrophysiology are particularly beneficial for correct pathogenic classification and are related to disease penetrance and severity. Here, we identify a novel SCN5A loss of function variant, p.1449Y>H, which presented with high penetrance and complete left bundle branch block, totally masking the typical findings on the electrocardiogram. We highlight the possibility of this overlap combination that makes impossible an electrocardiographic diagnosis and, through a functional analysis, associate the p.1449Y>H variant to SCN5A pathogenicity.

Functional cross-talk between phosphorylation and disease-causing mutations in the cardiac sodium channel Nav1.5

The cardiac sodium channel (Na_v1.5) is crucial for generating a regular heartbeat. It is thus not surprising that Na_v1.5 mutations have been linked to life-threatening arrhythmias. Interestingly, Na_v1.5 activity can also be altered by posttranslational modifications, such as tyrosine phosphorylation. Our combination of protein engineering and molecular modeling has revealed that the detrimental effect of a long QT3 patient mutation is only exposed when a proximal tyrosine is phosphorylated. This suggests a dynamic cross-talk between the genetic mutation and a neighboring phosphorylation, a phenomenon that could be important in other classes of proteins. Additionally, we show that phosphorylation can affect the channel’s sensitivity toward clinically relevant drugs, a finding that may prove important when devising patient-specific treatment plans.

The voltage-gated sodium channel Na_v1.5 initiates the cardiac action potential. Alterations of its activation and inactivation properties due to mutations can cause severe, life-threatening arrhythmias. Yet despite intensive research efforts, many functional aspects of this cardiac channel remain poorly understood. For instance, Na_v1.5 undergoes extensive posttranslational modification in vivo, but the functional significance of these modifications is largely unexplored, especially under pathological conditions. This is because most conventional approaches are unable to insert metabolically stable posttranslational modification mimics, thus preventing a precise elucidation of the contribution by these modifications to channel function. Here, we overcome this limitation by using protein semisynthesis of Na_v1.5 in live cells and carry out complementary molecular dynamics simulations. We introduce metabolically stable phosphorylation mimics on both wild-type (WT) and two pathogenic long-QT mutant channel backgrounds and decipher functional and pharmacological effects with unique precision. We elucidate the mechanism by which phosphorylation of Y1495 impairs steady-state inactivation in WT Na_v1.5. Surprisingly, we find that while the Q1476R patient mutation does not affect inactivation on its own, it enhances the impairment of steady-state inactivation caused by phosphorylation of Y1495 through enhanced unbinding of the inactivation particle. We also show that both phosphorylation and patient mutations can impact Na_v1.5 sensitivity toward the clinically used antiarrhythmic drugs quinidine and ranolazine, but not flecainide. The data highlight that functional effects of Na_v1.5 phosphorylation can be dramatically amplified by patient mutations. Our work is thus likely to have implications for the interpretation of mutational phenotypes and the design of future drug regimens.

Specific decreasing of Na+ channel expression on the lateral membrane of cardiomyocytes causes fatal arrhythmias in Brugada syndrome

Reduced cardiac sodium (Na+) channel current (INa) resulting from the loss-of-function of Na+ channel is a major cause of lethal arrhythmias in Brugada syndrome (BrS). Inspired by previous experimental studies which showed that in heart diseases INa was reduced along with expression changes in Na+ channel within myocytes, we hypothesized that the local decrease in INa caused by the alteration in Na+ channel expression in myocytes leads to the occurrence of phase-2 reentry, the major triggering mechanism of lethal arrhythmias in BrS. We constructed in silico human ventricular myocardial strand and ring models, and examined whether the Na+ channel expression changes in each myocyte cause the phase-2 reentry in BrS. Reducing Na+ channel expression in the lateral membrane of each myocyte caused not only the notch-and-dome but also loss-of-dome type action potentials and slowed conduction, both of which are typically observed in BrS patients. Furthermore, the selective reduction in Na+ channels on the lateral membrane of each myocyte together with spatial tissue heterogeneity of Na+ channel expression caused the phase-2 reentry and phase-2 reentry-mediated reentrant arrhythmias. Our data suggest that the BrS phenotype is strongly influenced by expression abnormalities as well as genetic abnormalities of Na+ channels.

Relationship between sodium channel function and clinical phenotype in SCN5A variants associated with Brugada syndrome

The identification of a pathogenic SCN5A variant confers an increased risk of conduction defects and ventricular arrhythmias (VA) in Brugada syndrome (BrS). However, specific aspects of sodium channel function that influence clinical phenotype have not been defined. A systematic literature search identified SCN5A variants associated with BrS. Sodium current (I_Na) functional parameters (peak current, decay, steady‐state activation and inactivation, and recovery from inactivation) and clinical features (conduction abnormalities [CA], spontaneous VA or family history of sudden cardiac death [SCD], and spontaneous BrS electrocardiogram [ECG]) were extracted. A total of 561 SCN5A variants associated with BrS were identified, for which data on channel function and clinical phenotype were available in 142. In the primary analysis, no relationship was found between any aspect of channel function and CA, VA/SCD, or spontaneous BrS ECG pattern. Sensitivity analyses including only variants graded pathogenic or likely pathogenic suggested that reduction in peak current and positive shift in steady‐state activation were weakly associated with CA and VA/SCD, although sensitivity and specificity remained low. The relationship between in vitro assessment of channel function and BrS clinical phenotype is weak. The assessment of channel function does not enhance risk stratification. Caution is needed when extrapolating functional testing to the likelihood of variant pathogenicity.

Functional evaluation of gene mutations in Long QT Syndrome: strength of evidence from in vitro assays for deciphering variants of uncertain significance

Background

Genetic screening is now commonplace for patients suspected of having inherited cardiac conditions. Variants of uncertain significance (VUS) in disease-associated genes pose problems for the diagnostician and reliable methods for evaluating VUS function are required. Although function is difficult to interrogate for some genes, heritable channelopathies have established mechanisms that should be amenable to well-validated evaluation techniques.

The cellular electrophysiology techniques of ‘voltage-’ and ‘patch-’ clamp have a long history of successful use and have been central to identifying both the roles of genes involved in different forms of congenital Long QT Syndrome (LQTS) and the mechanisms by which mutations lead to aberrant ion channel function underlying clinical phenotypes. This is particularly evident for KCNQ1, KCNH2 and SCN5A, mutations in which underlie > 90% of genotyped LQTS cases (the LQT1-LQT3 subtypes). Recent studies utilizing high throughput (HT) planar patch-clamp recording have shown it to discriminate effectively between rare benign and pathological variants, studied through heterologous expression of recombinant channels. In combination with biochemical methods for evaluating channel trafficking and supported by biophysical modelling, patch clamp also provides detailed mechanistic insight into the functional consequences of identified mutations. Whilst potentially powerful, patient-specific stem-cell derived cardiomyocytes and genetically modified animal models are currently not well-suited to high throughput VUS study.

Conclusion

The widely adopted 2015 American College of Medical Genetics (ACMG) and Association for Molecular Pathology (AMP) guidelines for the interpretation of sequence variants include the PS3 criterion for consideration of evidence from well-established in vitro or in vivo assays. The wealth of information on underlying mechanisms of LQT1-LQT3 and recent HT patch clamp data support consideration of patch clamp data together (for LQT1 and LQT2) with information from biochemical trafficking assays as meeting the PS3 criterion of well established assays, able to provide ‘strong’ evidence for functional pathogenicity of identified VUS.

Update on Genetic Basis of Brugada Syndrome: Monogenic, Polygenic or Oligogenic?

Brugada syndrome is a rare inherited arrhythmogenic disease leading to ventricular fibrillation and high risk of sudden death. In 1998, this syndrome was linked with a genetic variant with an autosomal dominant pattern of inheritance. To date, rare variants identified in more than 40 genes have been potentially associated with this disease. Variants in regulatory regions, combinations of common variants and other genetic alterations are also proposed as potential origins of Brugada syndrome, suggesting a polygenic or oligogenic inheritance pattern. However, most of these genetic alterations remain of questionable causality; indeed, rare pathogenic variants in the SCN5A gene are the only established cause of Brugada syndrome. Comprehensive analysis of all reported genetic alterations identified the origin of disease in no more than 40% of diagnosed cases. Therefore, identifying the cause of this rare arrhythmogenic disease in the many families without a genetic diagnosis is a major current challenge in Brugada syndrome. Additional challenges are interpretation/classification of variants and translation of genetic data into clinical practice. Further studies focused on unraveling the pathophysiological mechanisms underlying the disease are needed. Here we provide an update on the genetic basis of Brugada syndrome.

Predicting functional effects of missense variants in voltage-gated sodium and calcium channels

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes.

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.

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.

The evolution of gene-guided management of inherited arrhythmia syndromes: Peering beyond monogenic paradigms towards comprehensive genomic risk scores

Inherited arrhythmia syndromes have traditionally been viewed as monogenic forms of disease whose pathophysiology is driven by a single highly penetrant rare genetic variant. Although an accurate depiction of a proportion of genetic variants, the variable penetrance frequently noted in genotype positive families and the presence of sporadic genotype negative cases have long highlighted a more nuanced truth being operative. Coupled with our more recent recognition that many rare variants implicated in inherited arrhythmia syndromes possess unexpectedly high allele frequencies within the general population, these observations have contributed to the realization that a spectrum of pathogenicity exists among clinically relevant genetic variants. Notably, variable mutation pathogenicity and corresponding variable degrees of penetrance emphasize a limitation of contemporary guidelines, which attempt to dichotomize genetic variants as pathogenic or benign. Recognition of the existence of low and intermediate penetrant variants insufficient to be causative for disease in isolation has served to emphasize the importance of additional genetic, clinical, and environmental factors in the pathogenesis of rare inherited arrhythmia syndromes. Despite being rare, it has also become increasingly evident that common genetic variants play critical roles in both heritable channelopathies and cardiomyopathies and in aggregate may even be the primary drivers in certain instances, such as genotype negative Brugada syndrome. Our growing realization that the genetic substrates of inherited arrhythmia syndromes have intricacies that extend beyond traditionally perceived monogenic paradigms has highlighted a potential value of leveraging more comprehensive genomic risk scores for predicting disease development and arrhythmic risk.

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.

Late sodium current blocker GS967 inhibits persistent currents induced by familial hemiplegic migraine type 3 mutations of the SCN1A gene

Background

Familial hemiplegic migraine (FHM) is a group of genetic migraine, associated with hemiparesis and aura. Three causative different genes have been identified, all of which are involved in membrane ion transport. Among these, SCN1A encodes the voltage-gated Na⁺ channel Nav1.1, and FHM caused by mutations of SCN1A is named FHM3. For 7 of the 12 known FHM3-causing SCNA1 mutations functional consequences have been investigated, and even if gain of function effect seems to be a predominant phenotype, for several mutations conflicting results have been obtained and the available data do not reveal a univocal FHM3 pathomechanism.

Methods

To obtain a more complete picture, here, we characterized by patch clamp approach the remaining 5 mutations (Q1489H, I1498M, F1499 L, M1500 V, F1661 L) in heterologous expression systems.

Results

With the exception of I1498M, all mutants exhibited the same current density as WT and exhibited a shift of the steady state inactivation to more positive voltages, an accelerated recovery from inactivation, and an increase of the persistent current, revealing that most FHM3 mutations induce a gain of function. We also determined the effect of GS967, a late Na⁺ current blocker, on the above mentioned mutants as well as on previously characterized ones (L1649Q, L1670 W, F1774S). GS967 inhibited persistent currents of all SCNA1 FMH3-related mutants and dramatically slowed the recovery from fast inactivation of WT and mutants, consistent with the hypothesis that GS967 specifically binds to and thereby stabilizes the fast inactivated state.

Simulation of neuronal firing showed that enhanced persistent currents cause an increase of ionic fluxes during action potential repolarization and consequent accumulation of K⁺ and/or exhaustion of neuronal energy resources. In silico application of GS967 largely reduced net ionic currents in neurons without impairing excitability.

Conclusion

In conclusion, late Na⁺ current blockers appear a promising specific pharmacological treatment of FHM3.

Genetic interpretation and clinical translation of minor genes related to Brugada syndrome

Brugada syndrome (BrS) is an inherited arrhythmogenic disease associated with sudden cardiac death. The main gene is SCN5A. Additional variants in 42 other genes have been reported as deleterious, although these variants have not yet received comprehensive pathogenic analysis. Our aim was to clarify the role of all currently reported variants in minor genes associated with BrS. We performed a comprehensive analysis according to the American College of Medical Genetics and Genomics guidelines of published clinical and basic data on all genes (other than SCN5A) related to BrS. Our results identified 133 rare variants potentially associated with BrS. After applying current recommendations, only six variants (4.51%) show a conclusive pathogenic role. All definitively pathogenic variants were located in four genes encoding sodium channels or related proteins: SLMAP, SEMA3A, SCNN1A, and SCN2B. In total, 33.83% of variants in 19 additional genes were potentially pathogenic. Beyond SCN5A, we conclude definitive pathogenic variants associated with BrS in four minor genes. The current list of genes associated with BrS, therefore, should include SCN5A, SLMAP, SEMA3A, SCNN1A, and SCN2B. Comprehensive genetic interpretation and careful clinical translation should be done for all variants currently classified as potentially deleterious for BrS.

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

Rare variants in the cardiac potassium channel K_V7.1 (KCNQ1) and sodium channel Na_V1.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 Na_V1.5 (n = 356) and K_V7.1 (n = 144) variants with in vitro characterized channel currents gathered from the literature. Here we investigated the ability to predict commonly reported Na_V1.5 and K_V7.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 Na_V1.5 peak current (adj. R² = 0.27) and K_V7.1 + KCNE1 half-maximal voltage of activation (adj. R² = 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 (I_Na-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 (PIP₃), which inhibits I_Na-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 I_Na-L in CHO cells. PIP₃ abrogated the late current phenotype in R1193Q cells but had no effect on ΔKPQ. Homozygous R1193Q hiPSC-CMs displayed increased I_Na-L and long action potentials with frequent triggered beats, which were reversed with the addition of PIP₃.

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.

Systematic re-evaluation of SCN5A variants associated with Brugada syndrome

BACKGROUND

A large number of SCN5A variants have been reported to underlie Brugada syndrome (BrS). However, the evidence supporting individual variants is highly heterogeneous.

OBJECTIVE

We systematically re-evaluated all SCN5A variants reported in BrS using the 2015 American college of medical genetics and genomics and the association for molecular pathology (ACMG-AMP) guidelines.

METHODS

A PubMed/Embase search was performed to identify all reported SCN5A variants in BrS. Standardized bioinformatic re-analysis (SIFT, PolyPhen, Mutation Taster, Mutation assessor, FATHMM, GERP, PhyloP, and SiPhy) and re-evaluation of frequency in the gnomAD database were performed. Fourteen ACMG-AMP rules were deemed applicable for SCN5A variant analysis.

RESULTS

Four hundred and eighty unique SCN5A variants were identified, the majority of which 425 (88%) were coding variants. One hundred and fifty-six of 425 (37%) variants were classified as pathogenic/likely pathogenic. Two hundred and fifty-eight (60%) were classified as variants of uncertain significance, while a further 11 (3%) were classified as benign/likely benign. When considering the subset of variants that were considered "null" variants separately, 95% fulfilled criteria for pathogenicity/likely pathogenicity. In contrast, only 17% of missense variants fulfilled criteria for pathogenicity/likely pathogenicity. Importantly, however, only 25% of missense variants had available functional data, which was a major score driver for pathogenic classification.

CONCLUSION

Based on contemporary ACMG-AMP guidelines, only a minority of SCN5A variants implicated in BrS fulfill the criteria for pathogenicity or likely pathogenicity.