Publish in this journal
Journal Information
Vol. 72. Issue 4.
Pages 324-332 (April 2019)
Share
Share
Download PDF
More article options
Visits
Not available
Vol. 72. Issue 4.
Pages 324-332 (April 2019)
Original article
DOI: 10.1016/j.rec.2018.03.012
Full text access
Digenic Heterozigosity in SCN5A and CACNA1C Explains the Variable Expressivity of the Long QT Phenotype in a Spanish Family
La expresividad variable del síndrome de QT largo de una familia española se explica por la heterocigosis digénica en SCN5A y CACNA1C
Visits
...
Paloma Nieto-Marína,1, Juan Jiménez-Jáimezb,1, David Tinaqueroa, Silvia Alfayatea, Raquel G. Utrillaa, María del Mar Rodríguez Vázquez del Reyc, Francesca Perinc, Geòrgia Sarquella-Brugadad, Lorenzo Monserrate, Josep Brugadad, Luis Tercedorb, Juan Tamargoa, Eva Delpóna,2,
Corresponding author
edelpon@med.ucm.es

Corresponding author: Departamento de Farmacología, Facultad de Medicina, Universidad Complutense de Madrid, Plaza de Ramón y Cajal s/n, 28040 Madrid, Spain.
, Ricardo Caballeroa,2
a Departamento de Farmacología, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
b Unidad de Arritmias, Hospital Universitario Virgen de las Nieves, Instituto de Investigación Biosanitaria de Granada, Granada, Spain
c Departamento de Cardiología Pediátrica, Hospital Universitario Virgen de las Nieves, Granada, Spain
d Departamento de Cardiología Pediátrica, Hospital San Joan de Déu, Hospitalet de Llobregat, Barcelona, Spain
e Departamento de Cardiología, Health in Code, A Coruña, Spain
This item has received
...
Visits
(Daily data update)
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (7)
Show moreShow less
Tables (1)
Table. Biophysical Properties of WT and Mutant Cav1.2 Channels
Additional material (1)
Abstract
Introduction and objectives

A known long QT syndrome-related mutation in Nav1.5 cardiac channels (p.R1644H) was found in 4 members of a Spanish family but only 1 of them showed prolongation of the QT interval. In the other 3 relatives, a novel missense mutation in Cav1.2 cardiac channels was found (p.S1961N). Here, we functionally analyzed p.S1961N Cav1.2 channels to elucidate whether this mutation regulates the expressivity of the long QT syndrome phenotype in this family.

Methods

L-type calcium current (ICaL) recordings were performed by using the whole-cell patch-clamp technique in Chinese hamster ovary cells transiently transfected with native and/or p.S1961N Cav1.2 channels.

Results

Expression of p.S1961N channels significantly decreased ICaL density. Using Ba as a charge carrier to suppress the Ca-dependent inactivation of Cav1.2 channels, we demonstrated that the mutation significantly accelerates the voltage-dependent inactivation of Cav1.2 channels decreasing the inactivation time constant. As a consequence, the total charge flowing through p.S1961N Cav1.2 channels significantly decreased. The effects of the p.S1961N Cav1.2 and p.R1644H Nav1.5 mutations alone or their combination on the action potential (AP) morphology were simulated using a validated model of the human ventricular AP. The p.S1961N Cav1.2 mutation shortens the AP duration and abrogates the prolongation induced by p.R1644H Nav1.5 channels.

Conclusions

The p.S1961N mutation in Cav1.2 channels decreased the ICaL, an effect which might shorten ventricular AP. The presence of the loss-of-function Cav1.2 mutation could functionally compensate the prolonging effects produced by the Nav1.5 gain-of-function mutation.

Keywords:
Long QT syndrome
CACNA1C
SCN5A
Cav1.2
Nav1.5
Patch-clamp
Channelopathies
INaL
ICaL
Abbreviations:
APD
IBa
ICaL
LQTS
VDI
Resumen
Introducción y objetivos

En 4 miembros de una familia española se identificó una mutación en los canales cardiacos Nav1.5 (p.R1644H) descrita ya y relacionada con el síndrome de QT largo con anterioridad. Sin embargo, solo 1 de los portadores presentaba el intervalo QT prolongado. En los otros 3 individuos se identificó una nueva mutación con cambio de sentido en los canales cardiacos Cav1.2 (p.S1961N). En este trabajo se analizaron las características funcionales de los canales p.S1961N Cav1.2 para averiguar si dicha mutación regula la expresividad del síndrome de QT largo en esta familia.

Métodos

La corriente de calcio tipo L (ICaL) se registró mediante la técnica de patch-clamp en células de ovario de hámster chino transfectadas transitoriamente con los canales cardiacos humanos en su forma nativa o mutada.

Resultados

La expresión de canales p.S1961N disminuye significativamente la densidad de la ICaL. Al sustituir el ion calcio por bario para suprimir la inactivación dependiente del calcio de los canales Cav1.2, se demostró que la mutación acelera significativamente la inactivación dependiente del voltaje de los canales Cav1.2 y disminuye la constante de tiempo de inactivación. Como consecuencia, la carga total que atraviesa los canales p.S1961N Cav1.2 disminuye significativamente. Los efectos que las mutaciones p.S1961N Cav1.2 y p.R1644H Nav1.5, por separado o en combinación, producen sobre las características de los potenciales de acción (PA) se simularon mediante un modelo matemático de PA ventriculares humanos. Los resultados demuestran que la mutación p.S1961N Cav1.2 abrevia la duración del PA y suprime la prolongación inducida por la mutación p.R1644H de los canales Nav1.5.

Conclusiones

La mutación p.S1961N en los canales Cav1.2 disminuye la ICaL, un efecto que podría abreviar la duración de los PA ventriculares humanos. La presencia de esta mutación que disminuye la función de los canales Cav1.2 compensa funcionalmente los efectos producidos por la mutación de los canales Nav1.5 que aumenta su función y prolonga la duración de los PA.

Palabras clave:
Síndrome de QT largo
CACNA1C
SCN5A
Cav1.2
Nav1.5
Patch-clamp
Canalopatías
INaL
ICaL
Full Text
INTRODUCTION

The long QT syndrome (LQTS) is an inherited primary arrhythmogenic syndrome characterized by a prolonged QT interval on the electrocardiogram (ECG) due to delayed ventricular repolarization and is associated with syncope, seizures, ventricular tachycardia (mainly torsade de pointes), and a high risk of sudden cardiac death.1,2 Ventricular action potential (AP) duration (APD) is the result of a balance between depolarizing (Na and Ca) and repolarizing (K) currents. Thus, both an increase in the depolarizing or a decrease in the repolarizing currents could give rise to the different types of LQTS.1,2

LQTS type 3 is due to gain-of-function mutations in the SCN5A gene, which encodes the α-subunit of the sodium channels (Nav1.5) responsible for the cardiac sodium current (INa).3 These mutations usually modify Nav1.5 channel kinetics by several biophysical mechanisms ultimately leading to an increase in the late component of the INa (INaL) that allows a persistent inward current during the plateau phase of the ventricular AP that lengthens the APD.3

CACNA1C, CACNB2, and CACNA2D1 genes encode the pore forming α-subunit (Cav1.2) and the β2 and α2δ ancillary subunits, respectively, that form the cardiac Ca channel which generates the L-type Ca current (ICaL).4 Loss-of-function mutations in CACNA1C and CACNB2 genes shorten the APD and have been associated with Brugada syndrome and short QT syndrome.5,6 Conversely, gain-of-function mutations in the CACNA1C gene lead to prolongation of the ventricular APD and LQTS type 8.7

In this study, we identified a previously characterized gain-of-function mutation in SCN5A (encoding p.R1644H Nav1.5)8 in 4 relatives of a Spanish family in which only 1 of the carriers had a LQTS type 3 phenotype. Genetic analysis of the other 3 relatives who had normal corrected QT (QTc) values demonstrated that they also carry a novel missense mutation in the CACNA1C gene encoding for p.S1961N Cav1.2 channels. Therefore, these family members harbour 2 variants in different genes (digenic heterozygosity). Thus, this study aimed to characterize the electrophysiological properties of p.S1961N Cav1.2 channels to elucidate whether the presence of this variant can ameliorate the APD prolongation and LQTS produced by the p.R1644H Nav1.5 mutant.

METHODSGenetic Analysis

Clinical evaluation, including 12-lead ECG, of the family members was performed at the Hospital Universitario Virgen de las Nieves (Granada, Spain). The study was approved by the local ethics committee and conforms to the principles outlined in the Declaration of Helsinki. Each participant gave written informed consent.9–11

The whole codifying sequence and the flanking intronic regions of KCNQ1, KCNH2, SCN5A, KCNJ2, KCNJ8, CACNA1C, AKAP9, KCNE1, and KCNE2 genes were sequenced using an Illumina 1500 Hiseq next-generation sequencing platform. Potentially pathogenic variants were confirmed using Sanger sequencing.9–11 The pathogenicity of the variants was predicted according to the current recommendations of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (see details in the “Genetic analysis” section of the supplementary material).12 Other exonic nonsynonymous variants found, not predicted as pathogenic, are described in the Table of the supplementary material.

Mutagenesis and Cell Transfection

The p.S1961N substitution in Cav1.2 (NP_000710.5) was introduced using the QuikChange Site-Directed Mutagenesis kit (Stratagene, United States) and confirmed by direct deoxyribonucleic acid (DNA) sequencing.5,9,11,13 Chinese hamster ovary cells were transiently transfected with the cDNA encoding native (WT) or mutated Cav1.2 channels and the α2δ and β subunits (1:1.7:4 ratio).

Patch-clamping

Currents were recorded using the whole-cell patch-clamp configuration following previously described methods.9,13 Uncompensated access resistance, cell capacitance, and peak ICaL amplitude generated by WT Cav1.2 channels were 1.8 ± 0.4 M?, 16.2 ± 1.2 pF, and -521 ± 60 pA (n = 25), respectively. Therefore, no significant voltage errors (< 5mV) due to series resistance were expected with the micropipettes used.

Mathematical Modelling of Ventricular Action Potential

To simulate the shapes of human ventricular AP, we employed the O’Hara-Rudy model previously validated and used for similar purposes.14 The model was run at different driving frequencies (0.1-3Hz) under baseline conditions (WT) or by incorporating the changes in INaL and ICaL produced by p.R1644H Nav1.5 and p.S1961N Cav1.2 mutations, respectively.

Statistical Analysis

Results are expressed as mean ± standard error of the mean. Paired or unpaired t test or 1-way ANOVA followed by the Newman-Keuls test were used to assess statistical significance where appropriate. To take into account repeated sample assessments, data were analyzed with multilevel mixed-effects models. A value of P < .05 was considered statistically significant. An expanded Methods section is available in the supplementary material.

RESULTSCase Description and Genetic Analysis

The proband (III.1) (arrow in Figure 1A) had a syncopal episode when he was 28 years old after a single dose of ciprofloxacin, a QT prolonging drug.15 This suggested that he might have a genetic predisposition to drug-induced QT prolongation. Conversely, his clinical evaluation revealed a normal ECG with a Bazzett QTc of 420ms (Figure 1B). However, both his mother (II.2) and grandmother (I.1), who had never experienced syncopal episodes, suddenly died at rest when they were 49 years old (clinical or genetic data were not available). Therefore, his family history of sudden cardiac death encouraged the evaluation of the proband's relatives. The proband's sister (III.2) showed an inferior atrial rhythm at 60 beats per minute, normal PR and QRS durations, and a slightly widened T wave (Figure 1A,1B). However, her QT properly shortened during exercise and the QTc was also normal (430ms). One of the proband's nephews (IV.1), a 13-year-old boy, had an ECG with a late T wave after a long flat ST-segment, which resulted in a markedly prolonged QTc (558ms) (Figure 1B). Conversely, IV.2 showed a QTc of 420ms with normal T wave (Figure 1B). The rest of the relatives studied (white symbols in Figure 1A) have been and are currently asymptomatic, with completely normal ECGs.

Figure 1.

A: pedigree of the family (females [○] and males [□]). The arrow indicates the proband. Diagonal lines indicate deceased patients. + and – indicate individuals with and without the p.R1644H Nav1.5 (blue) and p.S1961N Cav1.2 (red) variants, respectively. Black symbols indicate long QT syndrome-affected individuals, while grey symbols depict individuals without genetic testing. White symbols represent studied individuals not affected by any clinical phenotype. B: electrocardiogram of III.1, III.2, IV.1, and IV.2. SCD, sudden cardiac death.

(0.44MB).

Considering the clinical data, a genetic analysis by means of next-generation sequencing (see Methods) was offered to both the proband, who initially refused it, and his sister, who accepted. Genotyping of the proband's sister identified 2 mutations (digenic heterozigosity). These mutations were thereafter sequenced, using the Sanger method, in the proband (who finally accepted) and the rest of the relatives studied. The first mutation, predicted as pathogenic (see the “Genetic analysis” section of the supplementary material), appeared in SCN5A (NC_000003.11:g.38592932C>T according to GRCh37) leading to the arginine-to-histidine substitution of the residue 1644 (p.R1644H) in Nav1.5 channels. This mutation was previously described in several families with LQTS type 3 and has never been related to any other arrhythmic phenotype, including the Brugada syndrome.8,16–18 Functional studies revealed that the p.R1644H mutation produces a marked increase in INaL (gain-of-function), which explains the APD and QT prolongation.8,16 The p.R1644H Nav1.5 mutation was also found in the proband and nephews, but not in any other of the relatives tested (Figure 1A). In the proband's sister, the proband, and 1 of his nephews (Figure 1A) a mutation in the CACNA1C gene was also identified (NC_000012.11:g.2797710G>A) (Figure 2B) that leads to a serine-to-asparagine substitution in the Cav1.2 channel at position 1961 (p.S1961N) (Figure 2A). The serine residue is highly preserved at the equivalent position among different species (Figure 2C) and the mutation has not been annotated previously in any public database. Importantly, the p.S1961N Cav1.2 mutation is also carried by II.4 and III.4 but both individuals had normal ECGs and have never experienced any arrhythmic event.

Figure 2.

A: scheme of the L-type cardiac calcium channel. B: DNA sequence chromatograms of IV.1 and IV.2 depicting the heterozygous variation (NC_000003.11:g.38592932C>T) of the SCN5A gene leading to the missense mutation p.R1644H in both individuals (top) and a WT sequence (left) and the heterozygous mutation (NC_000012.11:g.2797710G>A) (right) of the CACNA1C gene, which results in the p.S1961N missense mutation (bottom). C: sequence alignment of the regions surrounding S1961 residue in Cav1.2 channels in several species. The box highlights the conservation of this residue. “*” identical residues. “:” and “.” indicate conservation between groups of strongly and weakly similar properties, respectively. DNA, deoxyribonucleic acid; WT, native form of the channel.

(0.37MB).

Therefore, family members who showed digenic heterozygosity do not exhibit LQTS type 3 phenotype, while IV.1, who is the only one who carries only the SCN5A mutation, does. Thus, we decided to functionally characterize the p.S1961N Cav1.2 mutation to elucidate whether the presence of this variant can ameliorate the APD prolongation and LQTS type 3 produced by the p.R1644H Nav1.5 mutant.

Functional Characterization of p.S1961N Cav1.2 Channels

Figure 3A shows ICaL traces recorded at potentials between -50 and +70mV in Chinese hamster ovary cells expressing human cardiac WT or p.S1961N Cav1.2 channels along with the α2δ and β subunits. Peak density of the currents generated by p.S1961N Cav1.2 channels (-26.4 ± 1.4 pA/pF; P < .05) was approximately 25% lower than that generated by WT channels (-35.5 ± 3.0 pA/pF) (Figure 3A, 3B). In another group of experiments, cells were cotransfected (0.5:0.5 ratio) with the cDNA encoding WT and p.S1961N Cav1.2 channels considering the heterozygous condition of the mutation carriers. Under these conditions, peak current density was not statistically different to that generated by WT and p.S1961N channels separately (Figure 3B; P > .05). Figure 3C shows the voltage dependence of the activation of WT, p.S1961N, and WT+p.S1961N Cav1.2 channels obtained by plotting the normalized conductance (estimated by using equation #1; see section patch-clamping of the supplementary material) as a function of the membrane potential. The fit of a Boltzmann function to the data, yielded the midpoint (Vh) and the slope (k) of the curves, which almost overlap (Figure 3C). Thus, Vh and k values of the 3 curves were not significantly different (Table). The fit of the equation #2 (see section patch-clamping of the supplementary material) to the current density-voltage data in each experiment yielded the Erev which averaged 74.9 ± 2.7mV for WT channels. The mutation, either expressed alone or in combination with WT channels, did not modify the Erev (P > .05) (Table).

Figure 3.

A: currents generated by WT or p.S1961N Cav1.2 channels. B,C: mean density-voltage (B) and conductance-voltage curves (C) for ICaL generated by WT, p.S1961N or WT+p.S1961N Cav1.2 channels. Each point represents the mean ± standard error of the mean of (n) experiments. In B, *P < .05 vs Cav1.2 WT. ICaL, L-type Ca current. WT, native form of the channel.

(0.23MB).
Table.

Biophysical Properties of WT and Mutant Cav1.2 Channels

      ActivationInactivation
Cav1.2 construct (ICaL)  Charge at +20 mV (pC/pF)  Erev(mV)  Vh(mV)  k  τ (ms)  Vh(mV)  k  τfast(ms)  τslow(ms) 
WT (n = 8)  -1.7 ± 0.2  74.9 ± 2.7  7.3 ± 1.4  6.4 ± 0.3  1.2 ± 0.03  -19.1 ± 1.2  7.0 ± 0.9  22.4 ± 0.3  89.4 ± 5.1 
p.S1961N (n = 8)  -1.0 ± 0.1a  73.2 ± 1.3  5.5 ± 1.3  6.8 ± 0.2  1.1 ± 0.1  -20.6 ± 2.2  7.8 ± 0.5  14.7 ± 1.1a  78.3 ± 6.0 
WT+p.S1961N (n = 9)  -1.4 ± 0.1a,b  76.1 ± 0.9  6.1 ± 1.6  6.9 ± 0.2  1.2 ± 0.1  -21.1 ± 2.3  7.9 ± 0.6  18.1 ± 1.0a,b  90.3 ± 3.9 
Cav1.2 construct (IBa)  Charge at +10 mV (pC/pF)  Erev(mV)  Vh(mV)  k  τ (ms)  Vh(mV)  k  τ (ms)   
WT (n = 13)  -4.5 ± 1.1  55.7 ± 0.9  -3.4 ± 1.3  4.9 ± 0.4  1.8 ± 0.1  -27.9 ± 4.4  8.7 ± 1.6  67.5 ± 3.3   
p.S1961N (n = 10)  -1.6 ± 0.3a  54.8 ± 1.7  -0.8 ± 1.3  5.1 ± 0.3  1.5 ± 0.2  -30.0 ± 4.6  7.9 ± 1.0  43.1 ± 3.9a   

Erev, reversal potential; IBa, current generated using Ba as a charge carrier; ICaL: L-type Ca current; τ, time constants of activation (ICaL)/inactivation (IBa) obtained by fitting a monoexponential function to the current traces recorded at +20mV. τfast and τslow, fast and slow time constants of inactivation (ICaL) obtained by fitting a biexponential function to the current traces recorded at +20mV; Vh and k, midpoint and slope values of the activation and inactivation curves; WT, native form of the channels.

Values are the mean ± standard error of the mean of (n) experiments in each group.

a

P < .05 vs WT.

b

P < .05 vs p.S1961N.

Figure 4A shows superimposed ICaL traces elicited when applying 500 ms-pulses from -80 to +20mV in cells expressing WT or p.S1961N Cav1.2 channels (currents generated by WT+p.S1961N channels were omitted for clarity). Activation kinetics of Cav1.2 channels was not modified by the p.S1961N mutation alone or cotransfected with WT channels (Table; P > .05). Conversely, the mutation significantly accelerated the time course of current inactivation (Figure 4A), decreasing the fast time constant of inactivation (τfast) calculated by the biexponential fit to the current decline (Figure 4B and Table). Interestingly, τfast of the current generated by WT+p.S1961N channels was intermediate (Figure 4B and Table). Therefore, to accurately analyze the effects produced by the p.S1961N mutation on the total Ca influx we measured the charge density generated by WT, p.S1961N, and WT+p.S1961N Cav1.2 channels at each membrane potential (Figure 4C). The density was calculated by normalizing the charge crossing the membrane, estimated from the integral of current traces, to the cell capacitance (methods of the supplementary material). Due to the combined effects of the mutation on the peak current density and on the inactivation kinetics, the charge was significantly reduced by the p.S1961N mutation either alone or together with WT (Figure 4C).

Figure 4.

A: superimposed current traces generated at +20mV by WT or p.S1961N Cav1.2 channels. B: fast and slow time constants of inactivation obtained at +20mV in cells expressing WT, p.S1961N or WT+p.S1961N Cav1.2 channels. C: charge estimated as the integral of the current traces recorded as a function of the membrane potential in the 3 groups. D: Top: currents generated by WT Cav1.2 channels using the protocol shown. Bottom: inactivation curves for WT, p.S1961N or WT+p.S1961N Cav1.2 channels. Each bar/point represents the mean ± standard error of the mean of ≥ 8 experiments. *P < .05 vs Cav1.2 WT. WT, native form of the channel.

(0.25MB).

Figure 4D (top) shows ICaL traces generated by WT Cav1.2 channels with the protocol used to assess the voltage dependence of inactivation and construct the inactivation curves (methods of the supplementary material) (Figure 4D bottom). Neither the Vh nor the k of the inactivation curves were modified by the p.S1961N mutation either alone or in combination with WT channels (Figure 4D and Table).

The p.S1961N Mutation Accelerated the Cav1.2 Voltage-dependent Inactivation

Inactivation of Cav1.2 channels consists of 2 independent mechanisms, the Ca-dependent (CDI) and the voltage-dependent inactivation (VDI). We tested the effects of p.S1961N mutation on the VDI by using Ba as charge carrier since Ba currents (IBa) do not exhibit CDI.19Figure 5A shows IBa traces generated in 2 Chinese hamster ovary cells transfected with WT or p.S1961N Cav1.2 channels, respectively. Inactivation decay of IBa defined by a monoexponential function was markedly accelerated by the mutation (Figure 5B and Table). Furthermore, the mutation significantly decreased the peak IBa and Ba charge densities, the decrease in the charge being greater than that of the peak current (Figure 5C and Figure 5D). Finally, the voltage dependence of IBa inactivation was not modified by the mutation (Figure 5E and Table). All these results strongly suggested that p.S1961N mutation accelerates the VDI of Cav1.2 channels.

Figure 5.

A: IBa traces generated by WT or p.S1961N Cav1.2 channels. B-E: inactivation time constants obtained at +10mV (B), mean current density- (C), charge-density voltage (D), and inactivation (E) curves for IBa generated by WT or p.S1961N Cav1.2 channels. Each point represents the mean ± standard error of the mean of 10 experiments. In B-D, *P < .05 vs Cav1.2 WT. WT, native form of the channel.

(0.27MB).
Simulation of the Effects of p.S1961N Mutation on the Cardiac Action Potential

The effects of the p.S1961N Cav1.2 and p.R1644H Nav1.5 mutations, alone or combined, on the AP characteristics were simulated using a validated model of human endocardial ventricular myocytes (Figure 6).14 The model was run in baseline conditions (WT) or by introducing the kinetic and conductance changes produced by the p.S1961N Cav1.2 mutation. To simulate the effects of the p.R1644H mutation the INaL conductance was increased ≈3-fold, as previously described.8Figure 6A shows the simulated AP traces obtained at 1Hz. p.R1644H Nav1.5 lengthened the APD20, APD50, and APD90 by 17.4%, 17.1%, and 15.2%, respectively. Due to the ICaL reduction, p.S1961N Cav1.2 shortened APD20, APD50, and APD90 by 21.7%, 14.3%, and 10.8%, respectively (Figure 6A). Interestingly, the AP traces corresponding to p.R1644H+p.S1961N and WT were almost undistinguishable, and only a small prolongation of the APD90 (4.3%) was observed. Figure 6B depicts the APD90 change (prolongation or shortening) induced by each mutation in cells driven at different frequencies (0.1-3Hz). Prolongation induced by p.R1644H Nav1.5 was maximum at low frequencies and decreased progressively at faster rates. The shortening induced by p.S1961N Cav1.2 was similar at frequencies between 0.1 and 2Hz (≈10%), decreasing slightly at 3Hz. In the presence of p.S1961N, the prolongation induced by p.R1644H was reduced, being completely suppressed at high frequencies (2 and 3Hz).

Figure 6.

A: action potential traces generated by the O’Hara-Rudy model in endocardial ventricular myocytes at 1Hz in baseline conditions (WT), or in the presence of the p.R1644H Nav1.5 and p.S1961N Cav1.2 mutations alone or in combination. B: percentage of APD90 prolongation (positive values) or shortening (negative values) induced by each mutation alone or in combination in endocardial cells driven at 0.1-3Hz. APD, action potential duration; WT, native form of the channel.

(0.14MB).
Therapeutic Implications of the Functional Results

Expert consensus recommendations on LQTS therapeutic interventions consider that INaL inhibitors, such as flecainide, can be useful as add-on therapy for the treatment of LQTS type 3 patients.20 Accordingly, we tested the safety and efficacy of flecainide by infusing the drug (2mg/kg) under a controlled situation, before the implementation of a treatment of the IV.1 patient with this drug. Figure 7A shows that in the gain-of-function SCN5A mutant carrier, flecainide shortened the QT interval (QTc decreased from 558 to 470ms) without inducing any arrhythmic response. Therefore, in this child, oral flecainide was added-on the maximal tolerated dose of atenolol (50mg od) for arrhythmic event prophylaxis.

Figure 7.

Electrocardiogram recordings (V1-V6) of IV.1 (top) and IV.2 (bottom) before (baseline) and after 10minutes of intravenous infusion of flecainide at 2mg/kg.

(0.17MB).

Conversely, loss-of-function CACNA1C mutations have been related to Brugada syndrome5 which is definitively diagnosed when a type I ST-segment elevation in precordial leads is observed after intravenous administration of flecainide. Since our functional analysis demonstrated that the p.S1961N Cav1.2 mutation reduces the ICaL density, we wondered whether the presence of the CACNA1C mutation would lead to a Brugada syndrome-like pattern in the presence of flecainide.11 Therefore, we performed a flecainide test in IV.2, who carries both the SCN5A and the loss-of-function CACNA1C mutation. Interestingly, infusion for 10minutes with flecainide (also 2mg/kg) in the patient IV.2 produced a QRS widening and a Brugada syndrome-like pattern, without modifying the QTc (Figure 7).

DISCUSSION

Here we functionally describe a novel missense mutation located at the distal C-terminal domain of Cav1.2 (p.S1961N). This mutation was found in a proband and 4 of his relatives belonging to a Spanish family with history of syncope and sudden cardiac death. The results demonstrated that the p.S1961N mutation decreased the ICaL density, which probably shortens the duration of the plateau phase of the cardiac AP.

The oldest nephew of the proband had LQTS type 3 and was found to carry the p.R1644H Nav1.5 mutation. However, his mother, his brother, and his uncle (the proband) do carry the p.R1644H Nav1.5 mutation but did not show LQTS features. This could be explained by the frequent variable LQTS expressivity among carriers of a mutation, which, in some families can be explained by the presence of multiple mutations in the same or different genes (compound or digenic heterozygosity, respectively) or by single nucleotide polymorphisms (SNPs).9,10,21 Indeed, the functional compensation between mutations or polymorphisms within the SCN5A gene has been previously described.10,22,23 Such is the case of the common Nav1.5 p.H558R polymorphism when it is presented in the same (cis) or in a different allele (trans) from the mutation.10,21,23 Furthermore, it has been shown that 2 SCN5A mutations in trans (compound heterozygosity) can interact, ameliorating their deleterious effects reciprocally.10,24 However, no additional polymorphisms or mutations were found in the SCN5A gene in any of the family members studied here.

Some SCN5A mutations can give rise to a wide spectrum of disease phenotypes among the different carriers, thus producing “overlapping syndromes”.11,25 In fact, the substitution of Arg by Cys at the same position (p.R1644C) of the Nav1.5 channels has been identified in families with either LQTS type 326 or Brugada syndrome.27 However, this is not the case in our family, in which besides the LQTS type 3 phenotype of patient IV.1, no other electrocardiographic features suggesting other arrhythmic phenotype were identified. This is consistent with data previously reported regarding other families carrying the p.R1644H gain-of-function mutation, which has exclusively been associated with LQTS type 3.17,18

p.S1961N Cav1.2 Mutation Reduces L-type Ca Current Density

At physiologically relevant membrane potentials, the peak ICaL density generated by p.S1961N Cav1.2 channels decreased by ∼25% and the total charge crossing the membrane by ∼40%. Conversely, when cells were simultaneously transfected with WT and p.S1961N Cav1.2 channels (0.5:0.5 ratio) to “reproduce” the heterozygous condition of the participants, the current decrease was lessened. These results agree with the mild-normal phenotype of the p.S1961N mutation carriers who exhibited normal ECG and QT intervals, either when they had it concurrently with the Nav1.5 mutation (III.1, III.2, IV.2) or when not (II.4 and III.4). Importantly, the hallmark of the mutation was the acceleration of the time course of Cav1.2 channel inactivation without any alteration of the voltage dependence of either activation or inactivation.

As mentioned, inactivation of Cav1.2 channels consists of CDI and VDI components.19,28 CDI is mainly controlled by Ca/calmodulin binding to an IQ motif and is greatly inhibited when Ba is used as charge carrier.19 Molecular determinants of VDI, which include, among others, the C-terminal domain of Cav1.2 channels, are multiple and complex.28 Previous results demonstrated that mutations in the calcineurin (CaN)-binding site at the Cav1.2 channel (residues 1913-1941; NP_000710.5) markedly increase the VDI rate.29 Furthermore, it was suggested that the CaN-binding site is the VDI regulatory motif of the Ca cardiac channel,30 which partially overlaps with the protein phosphatase 2A (PP2A)-binding site (residues 1928-1970).31 According to this, the p.S1961 residue lies in the C-terminal region of the PP2A-binding site. Our results demonstrated that the p.S1961N mutation mainly increases the VDI rate of Cav1.2 channels, and thus, it could be speculated that the mutation enhances the PP2A binding and/or function. Indeed, some studies have suggested that PP2A-mediated dephosphorylation of Cav1.2 channels inhibits their function by antagonizing the cAMP-dependent protein kinase-mediated ICaL increase.31 Finally, the novel loss-of-function mutation in CACNA1C gene that we have described might, by removing the serine at position 1961, eliminate a putative phosphorylation site. Indeed, this serine is predicted as a putative phosphorylation site for phosphoinositide-dependent kinase 1 (PDPK1) by bioinformatic predictions (Group-based prediction System 3.0), although it has not been demonstrated experimentally. Therefore, further studies are needed to determine the actual mechanism by which the mutation causes an acceleration of the VDI of Cav1.2 channels.

Clinical Implications of the Functional Findings in this Family

The results of the mathematical model confirmed that the INaL increase produced by the p.R1644H Nav1.5 mutation could be responsible for the prolongation of the APD and QT interval duration observed in IV.1. Conversely, the slight decrease in the Ca influx during the plateau phase of the human cardiac AP produced by the p.S1961N Cav1.2 mutation could compensate the increase in the Na influx produced by the p.R1644H Nav1.5 mutation. Consequently, those family members that harbour both mutations would exhibit “almost normal” ventricular APs and ECG. However, we must keep in mind that the proband, who carries both mutations, had a syncopal episode when he was under treatment with ciprofloxacin, an antibiotic that increases the QT interval.15 Furthermore, both the mother and the grandmother of the proband, who seem to be “obligate carriers” of both mutations (the proband's father does not have any of the mutations), experienced sudden cardiac death. Of note, in patient IV.2, who also carries both mutations, flecainide infusion produced a noticeable QRS widening and a Brugada syndrome-like pattern, without any modification of the QT duration. Therefore, it can be hypothesized that the “functional compensation” between the gain-of-function and loss-of-function mutations of the Na and Ca channels, respectively, seems to be limited to conditions in which the genetic background of the patients does not concur with other factors that either prolong the ventricular repolarization (such as drugs, bradycardia, and electrolyte disturbances) or decrease cardiac excitability and/or shorten the duration of the plateau phase of the ventricular AP (such as flecainide). Therefore, it seems that in the patients carrying both mutations, either torsades de pointes- or Brugada syndrome-associated arrhythmias could be generated as a function of different proarrhythmic factors.

CONCLUSIONS

These results further support the contention that variable penetrance and phenotypic expressivity of arrhythmogenic inherited syndromes are partially attributable to genetic factors. In this context, functional analysis could help to select phenotype-based therapy in each family member.

FUNDING

This work was supported by Fondos Europeos de Desarrollo Regional; Ministerio de Economía y Competitividad [SAF2014-58769-P; SAF2017-88116-P]; Comunidad Autónoma de Madrid [B2017/BMD-3738]; Instituto de Salud Carlos III [PI16/00398]; The ERA-Net for Research on Rare Diseases (AC14/00029), Mutua Madrileña and Banco Bilbao Vizcaya Argentaria Foundations, and the Spanish Society of Cardiology Grants.

CONFLICTS OF INTEREST

None declared.

WHAT IS KNOWN ABOUT THE TOPIC?

At least 15 genes have been linked so far to LQTS, an inherited primary arrhythmogenic syndrome characterized by a prolonged QT interval that increases the risk of sudden cardiac death due to ventricular fibrillation. LQTS shows marked phenotypic expressivity, which has been mainly attributed to demographic factors and greatly hampers risk stratification in patients and therapeutic decisions. However, in some families, the variable expressivity is due to genetic factors and can be explained by the presence of multiple mutations in the same or different genes (compound or digenic heterozygosity, respectively) or by SNPs.

WHAT DOES THIS STUDY ADD?

We functionally analyzed a mutation in the CACNA1C gene encoding p.S1691N Cav1.2 channels found in digenic heterozygosity in a Spanish family. The results demonstrated that this mutation decreased the Ca influx during the plateau phase of the AP, an effect that functionally counterbalances the increase in the Na influx produced by the long QT-associated SCN5A mutation encoding p.R1644H Nav1.5 channels. As a consequence, only the family member carrying the SCN5A mutation exclusively exhibits LQTS type 3. Therefore, phenotypic expression of LQTS could be modulated by genetic factors and the functional analysis of concurrent mutations could help to guide a personalized therapeutic approach.

.

Acknowledgements

We thank Sandra Sacristán, Lorena Ondo, and Paloma Vaquero for their invaluable technical assistance.

References
[1]
Y. Nakano, W. Shimizu.
Genetics of long-QT syndrome.
J Hum Genet., 61 (2016), pp. 51-55
[2]
M.S. Bohnen, G. Peng, S.H. Robey, et al.
Molecular pathophysiology of congenital Long QT Syndrome.
Physiol Rev., 97 (2017), pp. 89-134
[3]
E. Savio-Galimberti, M. Argenziano, C. Antzelevitch.
Cardiac Arrhythmias Related to Sodium Channel Dysfunction.
Handb Exp Pharmacol., (2017),
Accessed 22 Feb 2018
[4]
J.S. Rougier, H. Abriel.
Cardiac voltage-gated calcium channel macromolecular complexes.
Biochim Biophys Acta., 1863 (2016), pp. 1806-1812
[5]
E. Burashnikov, R. Pfeiffer, H. Barajas-Martínez, et al.
Mutations in the cardiac L-type calcium channel associated with inherited J-wave syndromes and sudden cardiac death.
Heart Rhythm., 7 (2010), pp. 1872-1882
[6]
D.M. Béziau, J. Barc, T. O’Hara, et al.
Complex Brugada syndrome inheritance in a family harbouring compound SCN5A and CACNA1C mutations.
Basic Res Cardiol., 109 (2014), pp. 446
[7]
M.J. Betzenhauser, G.S. Pitt, C. Antzelevitch.
Calcium channel mutations in cardiac arrhythmia syndromes.
Curr Mol Pharmacol., 8 (2015), pp. 133-142
[8]
R. Dumaine, Q. Wang, M.T. Keating, et al.
Multiple mechanisms of Na+ channel--linked long-QT syndrome.
Circ Res., 78 (1996), pp. 916-924
[9]
R. Caballero, R.G. Utrilla, I. Amorós, et al.
Tbx20 controls the expression of the KCNH2 gene and of hERG channels.
Proc Natl Acad Sci USA., 114 (2017), pp. E416-E425
[10]
L. Núñez, A. Barana, I. Amorós, et al.
p.D1690N Nav1.5 rescues p.G1748D mutation gating defects in a compound heterozygous Brugada syndrome patient.
Heart Rhythm., 10 (2013), pp. 264-272
[11]
P. Dolz-Gaitón, M. Núñez, L. Núñez, et al.
Functional characterization of a novel frameshift mutation in the C-terminus of the Nav1.5 channel underlying a Brugada syndrome with variable expression in a Spanish family.
PLoS One., 8 (2013), pp. 1-15
[12]
S. Richards, N. Aziz, S. Bale, et al.
Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology.
Genet Med., 17 (2015), pp. 405-424
[13]
M. Matamoros, M. Pérez-Hernández, G. Guerrero-Serna, et al.
Nav1.5 N-terminal domain binding to α1-syntrophin increases membrane density of human Kir2.1, Kir2.2 and Nav1.5 channels.
Cardiovasc Res., 110 (2016), pp. 279-290
[14]
T. O’Hara, L. Virág, A. Varró, Y. Rudy.
Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation.
PLoS Comput Biol., 7 (2011), pp. e1002061
[15]
J. Tamargo.
Drug-induced torsade de pointes: from molecular biology to bedside.
Jpn J Pharmacol., 83 (2000), pp. 1-19
[16]
D.W. Wang, K. Yazawa, A.L. George, P.B. Bennett.
Characterization of human cardiac Na+ channel mutations in the congenital long QT syndrome.
Proc Natl Acad Sci USA., 93 (1996), pp. 13200-13205
[17]
I. Splawski, J. Shen, K.W. Timothy, et al.
Spectrum of mutations in long-QT syndrome genes. KVLQT1, HERG, SCN5A, KCNE1, and KCNE2.
Circulation., 102 (2000), pp. 1178-1185
[18]
S. Kapa, D.J. Tester, B.A. Salisbury, et al.
Genetic testing for long-QT syndrome: distinguishing pathogenic mutations from benign variants.
Circulation., 120 (2009), pp. 1752-1760
[19]
B.Z. Peterson, C.D. DeMaria, J.P. Adelman, D.T. Yue.
Calmodulin is the Ca2+ sensor for Ca-dependent inactivation of L-type calcium channels.
Neuron., 22 (1999), pp. 549-558
[20]
S.G. Priori, A.A. Wilde, M. Horie, et al.
Executive summary: HRS/EHRA/APHRS expert consensus statement on the diagnosis and management of patients with inherited primary arrhythmia syndromes.
Heart Rhythm., 10 (2013), pp. e85-e108
[21]
C. Napolitano, V. Novelli, M.D. Francis, S.G. Priori.
Genetic modulators of the phenotype in the long QT syndrome: state of the art and clinical impact.
Curr Opin Genet Dev., 33 (2015), pp. 17-24
[22]
S. Poelzing, C. Forleo, M. Samodell, et al.
SCN5A polymorphism restores trafficking of a Brugada syndrome mutation on a separate gene.
Circulation., 114 (2006), pp. 368-376
[23]
P.C. Viswanathan, D.W. Benson, J.R. Balser.
A common SCN5A polymorphism modulates the biophysical effects of an SCN5A mutation.
J Clin Invest., 111 (2003), pp. 341-346
[24]
J. Clatot, A. Ziyadeh-Isleem, S. Maugenre, et al.
Dominant-negative effect of SCN5A N-terminal mutations through the interaction of Na(v)1.5 α-subunits.
Cardiovasc Res., 96 (2012), pp. 53-63
[25]
C.A. Remme, A.A. Wilde, C.R. Bezzina.
Cardiac sodium channel overlap syndromes: different faces of SCN5A mutations.
Trends Cardiovasc Med., 18 (2008), pp. 78-87
[26]
C. Napolitano, S.G. Priori, P.J. Schwartz, et al.
Genetic testing in the long QT syndrome: development and validation of an efficient approach to genotyping in clinical practice.
JAMA., 294 (2005), pp. 2975-2980
[27]
A. Frustaci, S.G. Priori, M. Pieroni, et al.
Cardiac histological substrate in patients with clinical phenotype of Brugada syndrome.
Circulation., 112 (2005), pp. 3680-3687
[28]
F. Hofmann, V. Flockerzi, S. Kahl, J.W. Wegener.
L-type Cav1.2 calcium channels: from in vitro findings to in vivo function.
Physiol Rev., 94 (2014), pp. 303-326
[29]
S. Tandan, Y. Wang, T.T. Wang, et al.
Physical and functional interaction between calcineurin and the cardiac L-type Ca2+ channel.
Circ Res., 105 (2009), pp. 51-60
[30]
M. Cohen-Kutner, Y. Yahalom, M. Trus, D. Atlas.
Calcineurin Controls Voltage-Dependent-Inactivation (VDI) of the Normal and Timothy Cardiac Channels.
Sci Rep., 2 (2012), pp. 366
[31]
D.D. Hall, J.A. Feekes, Arachchige, et al.
Binding of protein phosphatase 2A to the L-type calcium channel Cav1.2 next to Ser1928, its main PKA site, is critical for Ser1928 dephosphorylation.
Biochemistry., 45 (2006), pp. 3448-3459

Both authors contributed equally to this work.

Both authors are co-senior authors.

Copyright © 2018. Sociedad Española de Cardiología
Idiomas
Revista Española de Cardiología (English Edition)

Subscribe to our newsletter

Article options
Tools
Supplemental materials
es en

¿Es usted profesional sanitario apto para prescribir o dispensar medicamentos?

Are you a health professional able to prescribe or dispense drugs?

es en
Política de cookies Cookies policy
Utilizamos cookies propias y de terceros para mejorar nuestros servicios y mostrarle publicidad relacionada con sus preferencias mediante el análisis de sus hábitos de navegación. Si continua navegando, consideramos que acepta su uso. Puede cambiar la configuración u obtener más información aquí. To improve our services and products, we use "cookies" (own or third parties authorized) to show advertising related to client preferences through the analyses of navigation customer behavior. Continuing navigation will be considered as acceptance of this use. You can change the settings or obtain more information by clicking here.