Biophysical classification of a CACNA1D de novo mutation as a high-risk mutation for a severe neurodevelopmental disorder
Molecular Autism volume 11, Article number: 4 (2020)
There is increasing evidence that de novo CACNA1D missense mutations inducing increased Cav1.3 L-type Ca2+-channel-function confer a high risk for neurodevelopmental disorders (autism spectrum disorder with and without neurological and endocrine symptoms). Electrophysiological studies demonstrating the presence or absence of typical gain-of-function gating changes could therefore serve as a tool to distinguish likely disease-causing from non-pathogenic de novo CACNA1D variants in affected individuals. We tested this hypothesis for mutation S652L, which has previously been reported in twins with a severe neurodevelopmental disorder in the Deciphering Developmental Disorder Study, but has not been classified as a novel disease mutation.
For functional characterization, wild-type and mutant Cav1.3 channel complexes were expressed in tsA-201 cells and tested for typical gain-of-function gating changes using the whole-cell patch-clamp technique.
Mutation S652L significantly shifted the voltage-dependence of activation and steady-state inactivation to more negative potentials (~ 13–17 mV) and increased window currents at subthreshold voltages. Moreover, it slowed tail currents and increased Ca2+-levels during action potential-like stimulations, characteristic for gain-of-function changes. To provide evidence that only gain-of-function variants confer high disease risk, we also studied missense variant S652W reported in apparently healthy individuals. S652W shifted activation and inactivation to more positive voltages, compatible with a loss-of-function phenotype. Mutation S652L increased the sensitivity of Cav1.3 for inhibition by the dihydropyridine L-type Ca2+-channel blocker isradipine by 3–4-fold.
Conclusions and limitations
Our data provide evidence that gain-of-function CACNA1D mutations, such as S652L, but not loss-of-function mutations, such as S652W, cause high risk for neurodevelopmental disorders including autism. This adds CACNA1D to the list of novel disease genes identified in the Deciphering Developmental Disorder Study. Although our study does not provide insight into the cellular mechanisms of pathological Cav1.3 signaling in neurons, we provide a unifying mechanism of gain-of-function CACNA1D mutations as a predictor for disease risk, which may allow the establishment of a more reliable diagnosis of affected individuals. Moreover, the increased sensitivity of S652L to isradipine encourages a therapeutic trial in the two affected individuals. This can address the important question to which extent symptoms are responsive to therapy with Ca2+-channel blockers.
In electrically excitable cells, Ca2+ inward current through voltage-gated Ca2+-channels (VGCCs; Cav) contributes to membrane depolarization and plays a key role in coupling electrical activity to intracellular Ca2+-dependent signaling processes (for review, see [1,2,3]). Therefore, VGCCs control essential physiological processes, such as hormone secretion, muscle contraction, sinoatrial node pacemaking, and sensory function. In the brain, they trigger neurotransmitter release, shape neuronal excitability, and couple excitation to gene expression associated with synaptic plasticity and different types of learning and memory [1, 3,4,5]. VGCC activity is fine-tuned to the specific requirements of cellular functions. This is accomplished by the functional heterogeneity and diverse subcellular targeting of ten pore-forming α1-subunit isoforms , several accessory β- and α2δ-subunits as well as by alternative splicing and post-translational modification (for a review, see ).
L-type Ca2+ channels (LTCCs; Cav1) form one of the three main families of VGCCs. From the four members (Cav1.1–Cav1.4), Cav1.2 and Cav1.3 are expressed in most electrically excitable cells, including a wide variety of brain regions . Both are located post-synaptically at somatodendritic locations and serve a key role in activity-dependent gene transcription [1, 3, 5]. Genetic variants in both channels have been associated with neuropsychiatric disease risk. Multiple common intronic single nucleotide polymorphisms (SNPs) in Cav1.2 (CACNA1C gene) LTCCs have consistently been associated with bipolar disorder and schizophrenia, although the consequences of these polymorphisms for Cav1.2 function remain unknown (for review, see [7,8,9]). Very rare de novo CACNA1C missense mutations cause Timothy Syndrome, a severe disease with lethal arrhythmias, facial dysmorphism, syndactyly and autism spectrum disorder (ASD) in surviving patients [10,11,12]. Together these findings have triggered new interest in clinical trials to repurpose LTCC blockers (“Ca2+-antagonists”), licensed as antihypertensive drugs since decades, also for the treatment of mood disorders .
We [14,15,16] and others [17,18,19,20] have recently provided accumulating evidence that de novo missense mutations in the pore-forming α1-subunit of Cav1.3 LTCCs (CACNA1D) confer high risk for neurodevelopmental disorders in humans. Symptoms range from ASD with (mutations A749G, Q547H [14, 20];) and without (G407R ) intellectual disability to more severely affected patients with seizures, muscle hypotonia, and global developmental delay (V401L ). Some patients also exhibit additional endocrine symptoms (primary aldosteronism or hyperinsulinism; G403D, I750M [17, 18]) due to the expression of Cav1.3 in adrenal zona glomerulosa cells and pancreatic β-cells (for a review, see ). Our studies revealed that all these genetic variants are not present in healthy parents or unaffected siblings and are absent in 141,456 reference genomes of controls without pediatric disease (gnomAD database, ). This is consistent with high penetrance and strongly supports a likely causative role of these mutations. Moreover, electrophysiological analysis of six mutations (from seven of these patients) after expression in HEK-293 cells revealed a highly consistent pattern of functional changes: they all induce gating changes that can enhance Cav1.3 Ca2+-current through these channels, in particular at subthreshold voltages. This gain-of-function is evident from a drastic slowing of channel inactivation and/or by facilitation of channel opening at more negative voltages . Therefore, the demonstration of such typical gain-of-function gating changes in functional studies may allow to distinguish likely pathogenic from non-pathogenic CACNA1D missense variants and help in the genetic diagnosis of individuals with neurodevelopmental disorders. This appears necessary because several genetic studies failed to classify CACNA1D missense variants as high-risk mutations and CACNA1D as a high-risk gene for neurodevelopmental disorders, including ASD [14, 15, 22, 23]. For example, gain-of-function CACNA1D mutation G407R in a patient with ASD has been identified, but has not been classified as high-risk mutation. However, functional analysis revealed typical gain-of-function changes, which strongly support its pathogenic potential .
In contrast to de novo gene-disrupting mutations (nonsense, splice site, frameshift), which cause a protein loss-of-function, the prediction of the pathogenic potential of missense variants is more difficult because in most cases their functional consequences cannot be predicted by bioinformatics tools. While our data argue for a high disease risk due to Cav1.3 gain-of-function, heterozygous de novo CACNA1D variants resulting in a loss of Cav1.3 activity are unlikely to cause human disease. This is strongly supported by previous findings both in knockout mice (for a review, see ) and Cav1.3-deficient humans with sinoatrial node dysfunction and deafness (SANDD; OMIM #614896 [24, 25]), in which functional loss of one or both CACNA1D alleles did not lead to a central nervous system (CNS) disease phenotype. This complicates the classification of new CACNA1D variants as high-risk mutations in genetic studies.
Here, we provide further convincing evidence for the high disease risk of gain-of-function de novo CACNA1D mutations for neurodevelopmental disorders. This is shown for mutation S652L, which has previously been identified in the Deciphering Developmental Disorders study, in a cohort of individuals with a severe developmental disorder of unknown cause ). However, in this study, it has not been classified as novel disease mutation with compelling evidence for pathogenicity and therefore CACNA1D has not been included as one of the 12 novel high-risk genes. Moreover, we demonstrate that a rare variant at the same position, S652W, induces a gating defect compatible with a loss-of-function, which explains its presence in apparently healthy individuals. Our data should raise awareness for the pathogenic potential of CACNA1D mutations, especially in patients without additional congenital endocrine symptoms as diagnostic features. De novo CACNA1D missense mutations may be underdiagnosed in clinical practice.
Complementary DNA constructs
Human wild-type (WT) Cav1.3 α1-subunits contained either exons 8a and 42 (WTL; long C-terminal splice variant; Genbank accession number: EU363339) or exons 8a and 43S (WTS; short C-terminal splice variant ;). All constructs were previously cloned into a pGFPminus vector containing a CMV promoter, an ampicillin resistance gene and no GFP tag as described [26, 27].
Cloning of S652 constructs S652LL, S652LS, and S652WL: To introduce mutations S652L or S652W into various Cav1.3 splice variants SOE PCR was used. Briefly, nt 1685-4059 of WTL or WTS were PCR amplified with overlapping primers (primer pair 1 and 2) introducing the point mutations C>T (Ser>Leu) or C>G (Ser>Trp) at position nt 1967 in separate PCR reactions (PCR a and b) using WTL or WTS as templates. The two separate PCR products were then used as templates for the final PCR reaction (PCR c) with primer pair 3. This fragment was then AauI/HindIII digested and cloned into respective sites of WTL or WTS yielding hCav1.3L S652L (S652LL) or hCav1.3L S652W (S652WL) and hCav1.3S S652L (S652LS). The following primer pairs were used for SOE PCR of S652L or S652W constructs (purchased from Eurofins MWG Operon, Ebersberg, Germany): primer pair 1, AauI (BrsGI) fwd: 5′-CCAACAAAGTCCTCTTGGCTCTGTTC-3′, S652L SOE rev: 5′-GATAATGAAGAGAAAAAGCAGAAGCAACAGCAAAGCGATGGACTTCATGGAGTTTAATAAG -3′ or S652W SOE rev: 5′-GATAATGAAGAGAAAAAGCAGAAGCAACAGCCAAGCGATGGACTTCATGGAGTTTAATAAG -3′ (314 bp); primer pair 2: S652L SOE fwd: 5′-CTTATTAAACTCCATGAAGTCCATCGCTTTGCTGTTGCTTCTGCTTTTTCTCTTCATTATC-3′ or S652W SOE fwd: 5′-CTTATTAAACTCCATGAAGTCCATCGCTTGGCTGTTGCTTCTGCTTTTTCTCTTCATTATC-3′, HindIII rev: 5′-ATAGATGAAGAACAGCATGGCTATGAGG-3′ (2122 bp); primer pair 3: AauI (BrsGI) fwd, HindIII rev (2375 bp).
PCR reaction mix for PCR a, b and c contained 5 μl 10x Pfu buffer with 20 mM MgSO4 (Cat# EP0571; Thermo Fisher Scientific,Waltham, MA, USA), 2 mM dNTP mix (Cat# R0241; Thermo Fisher Scientific, Waltham, MA, USA), 2.5 μl DMSO, 500 ng DNA template (PCR a and b) or 0.5 μl of PCR products a and b (PCR c), 10 pMol/μl forward primer, 10 pMol/μl reverse primer, 0.5 μl Pfu polymerase (native) (2.5 units/μl; Cat# EP0571; Thermo Fisher Scientific, Waltham, MA, USA) and nuclease-free water to a final volume of 50 μl. Following PCR program was performed: initial denaturation at 95 °C for 3 min, then 35 cycles of 30 s denaturation at 95 °C, 30 s annealing at 50 °C, variable extension time at 72 °C depending on fragment size (for Pfu polymerase: 2 min/kb) followed by a final elongation step of 72 °C for 7 min. The integrity of all cloned constructs was confirmed by restriction site mapping and sequencing (Eurofins MWG Operon, Ebersberg, Germany).
Cell culture and transfection
For whole-cell patch-clamp recordings, tsA-201 cells (a human embryonic kidney (HEK)-293 subclone stably expressing SV40 temperature-sensitive T-antigen, ECACC, 96121229) were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Cat# D6546; Merck KGaA, Darmstadt, Germany) containing 4500 mg/l L-glucose, 10% fetal bovine serum (FBS; Cat# 10270106; Thermo Fisher Scientific, Waltham, MA, USA), 2 mM L-glutamine (Cat# 25030032; Thermo Fisher Scientific, Waltham, MA, USA), 10 units/ml penicillin G (Cat# P-3032; Merck KGaA, Darmstadt, Germany), 10 μg/ml streptomycin (Cat# S-6501; Merck KGaA, Darmstadt, Germany) and maintained at 37 °C in a humidified incubator with 5% CO2. Cells were grown to ~ 80% confluency and split using 0.05% trypsin for cell dissociation. Cells were transiently transfected using the Ca2+-phosphate precipitation method always including EGFP (1.5 μg) as a transfection marker. For recordings of WTL vs S652LL or S652WL tsA-201 cells were transiently transfected with human ɑ1 (3 μg), rat β3 (2 μg; Genbank accession number NM_012828), and rabbit ɑ2δ-1 (2.5 μg, Genbank accession number NM_001082276) subunits whereas for recordings of WTS vs S652LS HEK-293 cells stably expressing β3 and ɑ2δ-1 were used and required only LTCC ɑ1 (3 μg) transient transfection [26, 28, 29]. HEK-293 cells stably expressing β3 and ɑ2δ-1 were periodically treated with selection agents for each subunit (β3, 500 μg/ml geneticin (Cat# 10131027; Thermo Fisher Scientific, Waltham, MA, USA); ɑ2δ-1, 10 μg/ml blasticidin S HCl (Cat# A1113903; Thermo Fisher Scientific, Waltham, MA, USA)). All data were obtained from > 3 independent transfections. On the following day, cells were trypsinized (0.05% trypsin) and plated onto poly-l-lysine-(Cat# P-2636; Merck KGaA, Darmstadt, Germany) precoated 35-mm culture dishes. Cells were kept at 30 °C and 5% CO2 and were subjected to electrophysiological experiments 20–72 h after transfection.
Electrophysiological recordings in tsA-201 cells
For whole-cell patch-clamp experiments, patch pipettes were pulled in a micropipette puller (Sutter Instrument, Novato, CA, USA) using borosilicate glass capillaries (borosilicate glass; Cat# 64-0792, Warner Instruments, Hamden, CT, USA) and fire-polished using a MF-830 microforge (Narishige Co, Tokyo, Japan). Pipettes with a resistance of 1.5–3 MΩ were backfilled with internal solution containing (in mM): 135 CsCl, 10 Cs-EGTA, 1 MgCl2, 10 HEPES, 4 Na2ATP adjusted to pH 7.4 with CsOH. The bath solution contained (in mM): 15 CaCl2 or 15 BaCl2, 150 Choline-Cl, 1 MgCl2, 10 HEPES, adjusted to pH 7.3 with CsOH. Whole-cell patch-clamp recordings were performed at room temperature (20-23 °C) using an Axopatch 200B Amplifier (Molecular Devices, San José, CA, USA). Data were digitized (Digidata, 1322A digitizer, Molecular Devices, San José, CA, USA) at 50 kHz, low-pass filtered at 1–5 kHz and analyzed using pClamp 10.2 software (Molecular Devices, San José, CA, USA). Series resistance was compensated by 60–90% and all voltages were corrected for a liquid junction potential of − 9.3 mV . Currents were leak subtracted either offline using a 50-ms hyperpolarizing voltage step from − 89 to − 99 mV or using an online P/4 protocol. Current-voltage (I-V) relationships were measured by applying 50 ms depolarizing square pulses to various test potentials (Δ 5 mV increments) starting from a holding potential (HP) of − 89 mV. I-V curves were fitted to the equation I = Gmax (V − Vrev)/(1 + exp [− (V − V0.5)/k]) where I is the peak current, Gmax is the maximum conductance, V is the test potential, Vrev is the extrapolated reversal potential, V0.5 is the half-maximal activation voltage, and k is the slope factor. The voltage dependence of activation was obtained from the I-V relationship by calculating the conductance (G = I/V − Vrev) followed by normalization (G/Gmax) and plotting as a function of voltage. The G-V curve was fitted using the following Boltzmann relationship: G = Gmax/(1 + exp[− (V – V0.5)/k]. The steady-state inactivation was determined by calculating the ratio between current amplitudes of a control versus a test pulse (I/Icontrol; both 20 ms to Vmax) separated by a 5-s conditioning step to various potentials (10 mV increments; 30 s intersweep interval; HP: − 89 mV) and plotting as a function of voltage. Steady-state inactivation curves were fitted using a modified Boltzmann equation: G = (1 – Gmax)/(1 + exp [(V – V0.5,inact)/kinact] + Gmax where V0.5inact is the half-maximal inactivation voltage and kinact is the inactivation slope factor. Channel open probability was estimated by dividing the peak ionic tail current (Itail) by the integrated “ON” gating charge (QON) at the potential where there is no ionic inward or outward current (Vrev). ON-gating currents were filtered at 5 kHz and digitized at 50 kHz. Window currents were obtained by multiplying the steady-state inactivation at a given voltage (Fig. 1c,d, Fig. 6b) with the corresponding current densities (pA/pF) at the given potentials of the I-V relationships (Fig. 1a,b, Fig. 6a). Pulses to physiologically relevant potentials i.e. − 20 mV which corresponds to ∼ − 35–38 mV in physiological Ca2+-concentrations were applied for 5-s. Normalized inward Ca2+-currents (ICa) were multiplied with the corresponding conductance at − 20 mV (G: WTL, 0.1253; S652LS, 0.4330; WTS, 0.2647; S652LS, 0.6325) and plotted as a function of time. The percentage of inactivation during a 5-s long depolarizing pulse from a HP of − 89 mV to the potential of maximal inward current (Vmax) was determined after 50, 100, 250, 500, 1000, and 5000 ms with Ca2+ or Ba2+ as a charge carrier. Ca2+-dependent inactivation (CDI) was determined over a broad voltage range by analyzing the fraction of remaining currents at the end of 250 ms depolarizations (expressed as fraction of the peak current amplitude, r250) to different test potentials (∆ 10-mV increments) with Ca2+ or Ba2+ as a charge carrier resulting in typical U-shaped dependence of voltage. The difference between Ca2+- and Ba2+-curves indicates the strength of CDI and is defined as parameter f for each voltage step and plotted as a function of voltage. The fractional Ca2+-dependent component of inactivation after 250 ms was calculated as CDI = 1 − rCa/rBa, where rCa and rBa is the fraction of current remaining at a given time point of inactivation, respectively. Persistent currents were determined after 5-s long depolarizations to different potentials expressed as fractional persistent current (%) normalized to the peak current amplitude measured by a 20-ms pre-pulse to the Vmax in the same sweep. Tail currents were obtained from a transition from + 80 mV to − 60 mV or − 40 mV and normalized tail currents were fitted to a bi-exponential equation. In general, experiments with currents < 100 pA and > 1000 pA (range was defined prospectively) were excluded from analysis to avoid errors in the measurements of V0.5,act, which tends to become more hyperpolarized at larger current amplitudes. For pharmacological experiments, cells were depolarized using a 100-ms square pulse to the Vmax of each individual cell (0.1 Hz; HP: − 89 mV). Cells were perfused using an air pressure-driven perfusion system (BPS-8 Valve Control System, ALA Scientific Instruments) with external bath solution (15 mM Ca2+) in the presence or absence of isradipine (Fisher scientific, 50-850-70001) with a flow rate of 0.5 ml/min. Isradipine stocks were prepared in DMSO and freshly diluted 1:1000 in bath solution to the final concentration prior to the experiment. On each recording day, individual control recordings with bath solution only were performed using the same tubes subsequently used for isradipine experiments. Drug application was started after at least three constant control sweeps during perfusion with bath solution. Drug effects were corrected for linear current decay (“run-down”) measured in control cells.
HEK-293 cells stably expressing β3 and α2δ-1 were transfected with WTL or S652LL Cav1.3 α1-subunits (3 μg) together with EGFP (1.5 μg) for visualization of transfected cells. Cells were patched 24–72 h after transfection at room temperature (21–23 °C) with internal solution containing (in mM): 114 CsMeSO3, 5 CsCl, 1 MgCl2, 4 Na2ATP, 10 HEPES and 0.5 Cs-EGTA (pH 7.3) and 0.2 mM Fluo-4 pentapotassium salt (cat#: F14200, Thermo Fisher). During the recording, cells were kept in Tyrode’s solution containing (in mM): 135 NaCl, 5.4 KCl, 1.8 CaCl2, 0.33 MgCl2, 0.33 NaH2PO4, 5 HEPES and 5 glucose (pH 7.4). Fluorescence from transfected cells was detected using Photon Technology International (PTI) photomultipliers and manufacturer’s software. Excitation was performed at 488 nm and fluorescent emission from each sample was recorded at 520 nm. Action potential waveform (APW)-like stimulus trains were applied at a frequency of 10 Hz and 300 sweeps per run. The APW protocol was elicited from – 80-mV HP composed of 3 voltage ramps: step: − 80 to − 60 mV for 2.5 ms, 1st ramp: − 60 to + 20 mV in 1 ms, 2nd ramp: + 20 to - 70 mV in 1.5 ms, 3rd ramp: − 70 to − 60 mV in 5 ms (afterhyperpolarization), step: − 60 mV for 90 ms. Fluorescent signals were normalized to baseline fluorescence (F0) and current density (pA/pF) determined in a ramp before the start of the train. Ca2+-charge of WTL and S652LL was obtained by integrating the area of Itail normalized to maximum ICa determined in a ramp before the start of the train.
Protein preparation and immunoblot analysis of HEK-293 cells
HEK-293 cells stably expressing β3 and ɑ2δ-1 were transiently transfected with WTL/S and mutant α1-subunits using JetPrime® transfection reagent (VWR International, Radnor, PA, USA) according to the manufacturer’s protocol. Membrane preparations were performed 48 h after transfection by first washing cells with phosphate-buffered saline (in mM: 137 NaCl, 2.7 KCl, 8 Na2HPO4, 1.5 KH2PO4). Subsequently, cells were harvested and resuspended in 2 ml lysis buffer (10 mM Tris-HCl, 1 μg/ml aprotinin, 0.1 mg/ml trypsin inhibitor, 1 μM pepstatin A, 0.5 mM benzamidine, 0.2 mM phenylmethylsulfonylfluoride, 2 mM iodacetamide, 1 μl/ml leupeptin, pH 7.4) and lysed on ice for 15 min. After resuspension, the mixture was homogenized by passing it through a 27-gauge cannula and centrifuged for 20 min at 726×g to remove cell debris. Membranes were collected by subjecting the resulting supernatant to an ultracentrifugation step at 110,561×g for 30 min. The pellet was resuspended in 200 μl of lysis buffer and stored at − 80 °C. Protein concentrations were measured via Bradford assay. Proteins were mixed with 4× NuPAGETM LDS sample buffer (Cat# NP0008; Thermo Fisher Scientific, Waltham, MA, USA) and incubated at 70 °C for 10 min before loading onto NuPAGETM 3–8% Tris-acetate protein gels (Thermo Fisher Scientific, Waltham, MA, USA) together with a PageRulerTM Plus prestained protein ladder (10-250 kDa; Cat# 26619; Thermo Fisher Scientific, Waltham, MA, USA). Proteins were separated using a 20× NuPAGETM Tris-acetate SDS running buffer (Cat# LA0041; Thermo Fisher Scientific, Waltham, MA, USA) at 50 V for 40 min and 100 V for 1 h. Proteins were transferred onto a polyvinyliden fluoride membrane (PVDF, Immobilon-P membrane, Cat# IPVH00010; Merck KGaA, Darmstadt, Germany) in transfer buffer (20x NuPAGETM Tris-acetate SDS running buffer, 20% Methanol (v/v)) at 100 V for 90 min. After the transfer, the membrane was blocked for 2 h at room temperature in blocking buffer (5% milk powder and 0.2% gelatine in 20 mM Tris, 150 mM NaCl, 0.5% Triton X-100, 0.1% Tween 20) and incubated with the primary antibody in blocking buffer at 4 °C overnight. α1-subunits were detected using anti-Cav1.3 antibody (rabbit-polyclonal; Cat# ACC-005, Lot# AN2150; Alomone labs, Jerusalem, Israel; diluted 1:100,000). Anti-ɑ-tubulin (mouse-monoclonal (DM1A); Cat# CP06, Lot# D00143511; Merck KGaA, Darmstadt, Germany; diluted 1:100,000) was used as loading control. After washing with washing buffer (20 mM Tris, 150 mM NaCl, 0.5% Triton X-100, 0.1% Tween 20) the membrane was incubated with the secondary antibody in blocking buffer for 2 h at room temperature followed by another wash step. Peroxidase conjugated goat anti-rabbit IgG (whole molecule; Cat# A0545, Lot# SLBC3108; Merck KGaA, Darmstadt, Germany; diluted 1:15,000) and peroxidase conjugated goat anti-mouse IgG (H+L; Cat# 31430, Lot# SC245915; Thermo Fisher Scientific, Waltham, MA, USA; diluted 1:5000) were used as secondary antibodies. For detection, the SuperSignalTM West Femto maximum sensitivity substrate (Cat# 34096; Thermo Fisher Scientific, Waltham, MA, USA) was used. Signals were quantified with a FUSION FX7 Imager (Vilber Lourmat Deutschland GmbH, Eberhardzell, Germany) and analysis of band intensity using ImageJ 1.46 (National Institute of Health). Integrated densities of mutant and WT signals were normalized to the loading control.
The absolute expression level of α1-subunit protein varies between membrane preparations of individual transfections. Therefore, transfections, membrane preparations, and Western blot analysis were always carried out in parallel for WT and mutant channels.
We predicted the structure of the WT Cav1.3 α1-subunit and two mutants by developing a homology model based on the cryo-electron microscopy (EM) structure of the Cav1.1 α1-subunit in the closed (potentially inactivated) state . The high sequence conservation of the Cav1.1 α1-subunit to the Cav1.3 α1-subunit (~ 75% similarity and ~ 62% identity measured by MOE) allowed us to predict a reliable structure model.
Homology modeling has been performed using Rosetta and MOE (Molecular Operating Environment, version 2018.08, Molecular Computing Group Inc., Montreal, Canada). Additionally, ab initio Rosetta was used to generate structures for loops that were not resolved in the original Cav1.1 α1-subunit template. The structures for the mutants were derived from the WT model by replacing the mutated residue and carrying out a local energy minimization using MOE. The C-terminal and N-terminal parts of each domain were capped with acetylamide (ACE) and N-methylamide to avoid perturbations by free charged functional groups. The structure model was embedded in a plasma membrane consisting of POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and cholesterol in a 3:1 ratio, using the CHARMM-GUI Membrane Builder. Water molecules and 0.15 M KCl were included in the simulation box. Energy minimizations of WT and mutant structures in the membrane environment were performed. The topology was generated with the LEaP tool of the AmberTools18, using force fields for proteins and lipids, ff14SBonlysc and Lipid14, respectively. The WT and mutant structures were gradually heated from 0 to 300 K in two steps, keeping the lipids fixed, and then equilibrated over 1 ns. Then molecular dynamics simulations were performed for 10 ns, with time steps of 2 fs, at 300 K and in anisotropic pressure scaling conditions. Van der Waals and short-range electrostatic interactions were cut off at 10 Å, whereas long-range electrostatics were calculated by the Particle Mesh Ewald (PME) method. MOE was used to visualize the key interactions and point out differences in the WT and mutant structures.
The Deciphering Developmental Disorders Study  has UK Research Ethics Committee approval (10/H0305/83, granted by the Cambridge South REC, and GEN/284/12 granted by the Republic of Ireland REC). Parental informed consent has been obtained for updated clinical information of the probands.
Data were analyzed using Clampfit 10.2 (Axon Instruments) and Sigma Plot 11 (Systat Software, Chicago, IL). For statistical analysis Graph Pad Prism 5.01 software (GraphPad Software, La Jolla, CA) was used. Significance of group differences between two groups was determined using unpaired Student’s t test for normally distributed data or Mann-Whitney U test for non-normally distributed data. Significance of group differences between three and more groups was determined using one-way analyses of variance (ANOVA) or two-way ANOVA for normally distributed data (with Bonferroni post-test as indicated). All data are represented as mean ± SEM. Significance level was set to α-error lower than p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***). All original datasets are available from the corresponding author on reasonable request.
S652 L, a novel CACNA1D de novo mutation in monozygotic twins with a severe neurodevelopmental disorder and ASD
The novel CACNA1D variant has been reported in the Deciphering Developmental Disorders Study . In a large unbiased genotype-driven approach this study analyzed 1,133 children with severe, undiagnosed developmental disorders and their parents using a combination of exome sequencing and array-based detection of chromosomal rearrangements . The CACNA1D variant (chr3: 53757881 C>T, human reference genome hg19) was predicted to cause a p.Ser652Leu (S652L) mutation (reference sequence NM_001128839). However, it was not considered a novel disease mutation in this study. The patients (DECIPHER database individuals #262954 and #262956; decipher.sanger.ac.uk) are male monozygotic twins, thirteen years old at present, and both harbor one copy of the mutation. Their clinical phenotype has been updated recently. It manifests as a severe neurodevelopmental phenotype with delayed speech and language development and a global developmental delay. Both individuals show self-injurious behaviors and have been diagnosed with ASD two years ago. One patient is tall (5 ft 5 in.) and heavy (57 kg) for his age, has undescended testes and suffered from seizures at the age of two without recurrence. In addition, he shows challenging behavior with attention deficit hyperactivity disorder-like symptoms. No abnormal blood pressure has so far been reported and both are currently not treated with any medication. In addition, facial dysmorphism characterized by epicanthus, abnormality of the nose, microtia, a small vermillion border, and widely spaced teeth was noted. Prediction of the possible impact of this amino acid substitution on protein structure and function using bioinformatics prediction tools PolyPhen2, SIFT and MutationTaster indicated a probably damaging (score: 1.00), deleterious (100%) or disease-causing (probability: 0.999) effect, respectively. This variant is not reported in the gnomAD database (a reference database that lists exomes and genomes from a total of 141,456 unrelated individuals harboring mutations without pediatric disease ).
A frameshift mutation (c.1934_1935insG (p.Glu646GlyfsTer)) in the KIF22 gene (chr16: 29816479 T>TG, human reference genome hg19) was also identified in both patients. This gene is highly expressed in bones, cartilage, skin, ligaments and joint capsules . Mutations in the KIF22 gene have so far not been associated with neurodevelopmental disorders but result in a syndrome called spondyloepimetaphyseal dysplasia with joint laxity (SEMDJL; OMIM #603213), with malformations of the spine, skeletal dysplasia and malalignment of limbs but no intellectual impairment. Complete knock-out of Kif22 in mice results in premature intrauterine death but surviving Kif22−/− embryos develop into healthy adult mice . Moreover, protein loss-of-function variants (stop gained, frameshift mutations) are reported in more than 40 controls in the gnomAD database. Therefore, the KIF22 mutation is unlikely to explain the severe neurodevelopmental phenotype of the two patients. In contrast, CACNA1D mutation S652L has also been reported as a somatic aldosterone-producing-adenoma (APA) mutation  in a patient suffering from resistant hypertension, providing additional indirect evidence for a pathogenic role of this novel germline CACNA1D variant.
Mutation S652L changes the voltage-dependence of activation and inactivation
Based on recent findings of unique gating changes induced by pathogenic CACNA1D de novo mutations [14,15,16,17,18, 20, 34], we therefore hypothesized that S652 L could also explain the neurodevelopmental phenotype in both patients. For a detailed biophysical characterization in tsA-201 cells we introduced this mutation into the biophysically distinct C-terminally long (WTL, S652LL) and short (WTS, S652LS) Cav1.3 splice variants . Both splice variants are abundantly expressed in the brain and differ considerably with respect to their biophysical properties, with higher voltage-sensitivity and faster Ca2+-dependent inactivation of WTS [3, 16]. Mutant α1-subunit proteins were expressed as intact proteins with the expected molecular mass (Additional file 1: Figure S1).
Mutation S652L induced pronounced gating changes. It significantly shifted the voltage-dependence of activation (Fig. 1a, b) and steady-state inactivation to more negative potentials indicating a phenotype that can support a channel gain-of-function in both splice variants (Fig. 1c, d; for parameters see Table 1) by promoting Ca2+-inward currents (ICa) at negative voltages. As a consequence the mutation induced a higher window current at subthreshold potentials (− 50 and − 40 mV) compared to WT in the short Cav1.3 splice variant (Fig. 1e, f). These variants comprise about half of the Cav1.3 α1-subunits in the brain . At − 50 mV significant window current was only measurable in S652 LS but not in WTS, and it was two times larger in the mutant at − 40 mV. Whereas current amplitudes were larger at negative voltages, ICa above the potential of maximal inward current (Vmax) were significantly smaller (Fig. 1a, b). Since gating currents, which are a measure for the number of active channels at the cell surface, were not different in the mutant channels (QON [mean ± SEM; pA*ms)]: WTL, 158.9 ± 26.3, n = 23; S652 LL, 140.3 ± 25.1, n = 21; Mann-Whitney test), reduced current density is likely due to a decreased open probability (PO). This is further supported by a significantly decreased slope of the Itail/QON relationship for S652 LL channels (Itail/QON [linear regression slopes, mean ± SEM; ms−1]: WTL: − 7.22 ± 0.916, r2 = 0.72, n = 26; S652 LL: − 4.24 ± 0.657, r2 = 0.72, n = 25; slopes are significantly different: F = 6.43, p = 0.015, F test; Additional file 2: Figure S2).
Note that short Cav1.3 splice variants have very small, non-measurable ON-gating currents  and were therefore not further analyzed here.
Mutation S652L accelerates voltage-dependent inactivation but has opposing effects on Ca2+-dependent inactivation
Since either acceleration or slowing of the inactivation time course of Cav1.3 is also a hallmark of pathogenic CACNA1D mutations, we studied voltage- (VDI) and Ca2+-dependent inactivation (CDI). Mutation S652L significantly accelerated inactivation kinetics (Fig. 2) during 5-s depolarizations to Vmax with both Ba2+ (which reports VDI) and Ca2+ (which, in addition, induces CDI) as charge carriers in both C-terminally long (Fig. 2a) and short (Fig. 2b; for statistics, see Table 2) splice variants. It also significantly reduced ICa during long-lasting depolarizations as shown as the percentage of remaining current after 5-s depolarizations to various test potentials in both Cav1.3 splice variants (Fig. 2c, d). By comparing the difference of fractional inactivation between ICa and inward Ba2+-current (IBa) after 250-ms depolarizations to different test potentials (f value, see Fig. 3) mutational effects on the voltage-dependence of CDI could also be determined. Whereas maximal CDI was unchanged in the long Cav1.3 splice variant (Fig. 3a, c), it was significantly reduced in the short variant (Fig. 3b, d). Therefore, faster inactivation of ICa must be due to the acceleration of VDI, despite being partially compensated by reduced CDI in WTS. Assuming that CDI and VDI are independent processes, we also calculated the fractional Ca2+-dependent component of inactivation as previously described  for WTS vs S652LS from the data shown in Table 2. After 250 ms of inactivation there was no difference between WTS (0.778 ± 0.027) and S652LS (0.805 ± 0.021; p = 0.40, n = 21; unpaired Student’s t test) and this was also true for all other time points. This further confirms that S652L can promote Cav1.3 inactivation largely by affecting VDI.
Mutation S652L increases Ca2+-signaling during stimulation protocols simulating neuronal firing patterns
In order to predict the consequences of mutation S652L on Ca2+-influx during neuronal activity patterns, we simulated its activity during sustained upstate potentials and during action potential firing.
At low potentials, Cav1.3 channels can contribute to the formation of plateau potentials due to their known negative activation range (for a review, see ). This was described in medium spiny neurons when transient upstate potentials were induced by glutamatergic excitatory input . To quantify changes of ICa amplitudes induced by the mutation at sustained subthreshold depolarizations we mimicked plateau potentials by prolonged depolarizations to − 20 mV with 15 mM Ca2+ as charge carrier. This voltage would correspond to ~− 35–38 mV at physiological Ca2+-concentrations  (Fig. 4a, b). The resulting ICa traces were normalized to maximal ICa at Vmax in individual cells and corrected for splicing- and mutation-dependent differences in V0.5 by multiplying with the corresponding conductance measured for the different WT and mutant constructs at − 20 mV (derived from activation curves in Fig. 1c, d). As shown in Fig. 4a, b, mutated channels caused a marked and highly significant increase of ICa over the first 300 ms of the depolarization compatible with enhanced channel activity suitable to support upstate potentials during this time period (for details, see legend to Fig. 4).
In order to predict the consequences of mutation S652L on Ca2+-influx during action potential firing we simultaneously measured ICa and cytosolic Ca2+-responses upon stimulation of HEK-293 cells transfected with WTL and S652LL with 10-Hz-trains of action potential-like waveforms (APW) (Fig. 5a; for details, see legend) [26, 38]. These stimuli resulted in typical ICa transients (Fig. 5a) with maximal ICa reached during the repolarization phase of the APW [29, 38, 39]. Peak ICa slowly decreased during trains, an effect that was significantly enhanced by the mutation (Fig. 5a, b). However, the mutation enhanced the increase in intracellular Ca2+-levels measured simultaneously during this stimulation protocol (Fig. 5c). To explain this discrepancy, we measured the total Ca2+-charge during the 30-s train. This was significantly higher in S652L mutated channels even after 15 s of stimulation (Fig. 5d; for details, see legend). It was not due to higher S652L channel expression, because we normalized data to current density for each cell. As shown in Fig. 5e, ICa amplitude during APW repolarization was significantly larger for S652L as compared to WTL (peak 1st AP [pA/pF]: WTL, − 12.43 ± 1.95, n = 19; S652LL, − 33.64 ± 3.13***, n = 21; unpaired Student’s t test, ***p < 0.001). This could be explained by the more negative activation voltage range as well as a pronounced slowing of ICa deactivation. Deactivation of tail currents following repolarizations from + 80 mV to − 60 mV or − 40 mV was significantly slower in S652LL as compared to WTL. This was primarily caused by a decrease in the contribution of the fast component and an increase of the slow component of the bi-exponential deactivation process (Fig. 5f; for statistics, see Table 3). Thus, higher Ca2+-levels and Ca2+-charge can be attributed to the slower deactivation kinetics and higher current amplitudes induced by mutation S652L during action potential like firing.
Mutation S652W produces a loss of Cav1.3 channel function
Our data predict that only CACNA1D mutations that are able to enhance Cav1.3 channel activity can confer high risk for neurodevelopmental symptoms. Accordingly, like any of the other previously described pathogenic variants, S652L has not been reported in the genomes of 141,456 control individuals free from pediatric disease (gnomAD database ;). In contrast, the pathogenic potential should be low or absent from missense mutations causing gating defects favoring reduced function, as outlined above, in mice [40, 41] and humans [24, 25]. Interestingly, the gnomAD database also reports the rare variant p.Ser672Trp (chr3: 53757881 C>G, human reference genome hg19) variant (S652W), located in the same position as S652L, in three healthy unrelated individuals. PolyPhen2 and SIFT predict a probably damaging (score: 1.00) and deleterious (100%) effect on protein function, respectively, but these algorithms cannot predict gating changes. This provided us with a unique opportunity to further test our above hypothesis by studying the biophysical properties of S652W. As shown in Fig. 6b, this mutation has opposite effects on the voltage-dependence of gating compared to S652L. Both steady-state activation and inactivation were significantly shifted to more positive voltages by 4–5 mV (for statistics see Table 4).
Consequently, window currents were not increased and even tended to be shifted to more positive voltages also compatible with a loss-of-function at threshold voltages (Fig. 6c). S652W also failed to enhance channel function by other mechanisms: it did neither slow inactivation kinetics (Fig. 6d; for statistics, see Table 5) nor change the fraction of non-inactivating current as evident from the steady-state inactivation analysis (Fig. 6b).
Molecular modeling of Cav1.3 WT, S652L, and S652W α1-subunits
On the structural level, the loss of a newly discovered inter-domain hydrogen bond connecting the S4–S5 linkers of repeats II and I could explain the gating differences between the two variants (Fig. 7). Our Cav1.3 channel homology model, based on the cryo-EM structure of the Cav1.1 α1-subunit , localizes S652 at the C-terminal end of the S4–S5 linker in channel repeat II (Fig. 7, upper, left). The S4–S5 linkers in each repeat are known to form contacts with the cytoplasmic end of their corresponding S6 helices, which together form the inner mouth of the channel (activation gate ;). This allows the S4–S5 linkers to transmit voltage-sensor movements to the activation gate. Our model predicts that S652 forms a hydrogen bond with S256 in the S4–S5 linkers of the neighboring repeat I (Fig. 7a). Therefore, this hydrogen bond connects the S4–S5 linkers in two different repeats and could be crucial for controlling the voltage-sensitivity of channel gating. This is supported by our finding that in the S652L mutant, this inter-domain interaction is not stabilized by any hydrogen bonds. Instead, leucine forms much weaker hydrophobic contacts with the residues V259, V260, and L261 located in the S4–S5 linkers of repeat I (Fig. 7b). In contrast, the aromatic side chain of the tryptophan in the S652 W mutant is capable of forming an inter-domain pi-H interaction with S256 in IS4–S5 as well as an intra-domain hydrogen bond with its own IIS4-S5 backbone (K648) (Fig. 7c), which should also allow a rigidifying effect similar as in the WT channel.
Together with our functional data, this reveals the importance of an inter-domain hydrogen bond for normal electromechanical coupling in Cav1.3 channels, which has not been described before. Weakening this interaction by substituting S652 with leucine causes a dramatic change in channel gating. Since we modeled the mutation with the activation gate in a closed channel conformation, this hydrogen bond likely stabilizes the channel in a closed state. Its weakening would favor the transition to and/or the stabilization in the open state, which can explain the strong shift of the voltage-dependence of activation to more negative potentials. Stabilization of the open state can also explain the slower transition from the open to the closed state evident as slowing of deactivation in the S652L mutant channel at a given voltage.
Mutation S652L increases the sensitivity of Cav1.3 channels for inhibition by the dihydropyridine LTCC blocker isradipine
The fact that pathogenicity is associated with enhanced channel function but that reduced Cav1.3 channel function in the brain is not associated with CNS symptoms (see above), make available LTCC blockers a potential therapeutic option for the symptomatic treatment of individuals affected by gain-of-function mutations.
These drugs, such as the dihydropyridines (DHPs) nifedipine, felodipine, or isradipine, are safely used since decades for the treatment of arterial hypertension and angina. Since DHPs preferentially bind to channels in an inactivated state [42, 43], mutations affecting the gating properties, such as S652L may change the sensitivity of the channel for inhibition by DHPs. In order to quantify mutation-induced changes on DHP sensitivity, we assessed inhibition of the C-terminally long WT and S652L mutant channels by the DHP isradipine using a standard square pulse protocol (100 ms to Vmax, 0.1 Hz, HP: − 89 mV). As illustrated in Fig. 8, S652L-mutated channels required significantly lower isradipine concentrations for channel inhibition with a 3–4-fold decrease of their half maximal inhibitory concentration (IC50, mean (95% coincidence interval); WTL: 60.3 (52.0 – 70.0) nM, S652LL: 18.1 (15.3 – 21.5) nM; Fig. 8). This finding suggests that DHPs with good brain penetrance, such as isradipine [29, 44] may preferentially inhibit S652L-mutated Cav1.3 channels. Based on their good clinical safety profile this encourages therapeutic trials with DHPs in the affected individuals.
In this study, we provide compelling evidence for the CACNA1D S652L variant as a high-risk and likely disease-causing mutation in two individuals of the Deciphering Developmental Disorders cohort (decipher.sanger.ac.uk) of children with severe, undiagnosed developmental disorders . This evidence builds on a detailed biophysical characterization that demonstrates gating changes able to also induce enhanced channel activity as is typical for six other de novo germline missense mutations in patients with ASD with and without other neurodevelopmental symptoms [14,15,16,17,18, 20]. Our data adds the CACNA1D gene to the other 12 developmental-disorder-linked genes identified in the Deciphering Developmental Disorders study and thus further increases its diagnostic yield. Moreover, since S652L has also been found in APAs as a somatic mutation we also confirm S652L as a disease-causing mutation in patients with primary aldosteronism . We also demonstrate that although heterozygous missense mutations in the same position can be damaging, they can have opposite effects on channel function and, in the case of Cav1.3 α1-subunits, may comprise low or no pathogenic risk when the resulting gating changes does not support a gain-of-function. This complicates predictions of the disease-causing role of rare de novo CACNA1D variants in genetic studies and emphasizes the need for functional analysis as described in this report. We clearly show this for S652W. Whereas neither S652L nor one of the other pathogenic CACNA1D mutations is reported in the gnomAD database, the heterozygous S652W variant is reported in three neurologically apparently healthy individuals in this database . In contrast to S652L, this variant tended to shift the window current to more positive rather than negative voltages, which is due to the positive shift of the voltage-dependence of activation and inactivation. This nicely fits our prediction that only CACNA1D mutations which can also support channel gain-of-function confer high risk for neurodevelopmental disorders. It is also in line with previous studies in mice (for a review, see ) and humans [24, 25]. These found that heterozygous loss-of-function of Cav1.3 is unlikely to cause symptomatic neurodevelopmental disorders and even homozygous loss of Cav1.3 function has not been reported to cause neuropsychiatric behavioral changes [3, 24, 25, 41]. Therefore, a unifying feature of all pathogenic mutations described so far is the potential to induce gating changes that can enhance Cav1.3 function during neuronal activity. The potential for enhanced channel function may, however, vary between different neurons. For example, in neurons firing from more depolarized membrane potentials, the negative shift in steady-state inactivation may also reduce the availability of Cav1.3 channels.
S652L adds to another six de novo missense mutations functionally characterized so far by us [14,15,16] and others [17, 18, 20, 38] in a total of seven patients with neurodevelopmental disorders. Functionally they fall into two major classes. Channel gain-of-function is either predominantly induced by stabilizing a large non-inactivating current component (type 1: in particular, G403D, G407R) or by inducing a strong shift of activation voltage to more negative potentials (type 2: A749G, I750M and V401L) [14,15,16,17], as also observed for S652L. Both types of mutations enhance intracellular Ca2+-load when expressed in HEK-293 (A749G, ; S652L, this paper) or GLT muscle cells (G407R, ). Pronounced negative shifts of activation and/or pronounced slowing of channel inactivation can therefore be taken as "diagnostic" feature for the pathogenicity of CACNA1D de novo mutations in patients with neurodevelopmental disease.
Importantly, our data also strongly suggest that pathogenicity can also be assumed if the same variant has also been reported as a somatic mutation in APAs. As shown here, this is the case for S652L  and likewise for G403D, I750M and V401L [17, 45].
We have also observed a small but significant shift of Vrev (Fig. 1, Table 1) to more negative voltages. This could indicate a potential change in ion selectivity by the mutation. Interestingly, we have also detected similar shifts by other Cav1.3 α1-subunit gain-of-function mutations . If confirmed by single-channel recordings, this may also contribute to the mutation-induced pathological signaling changes.
Very recently identified de novo mutation V259A has been reported in another severely affected individual, a 1-year-old male with seizures, global developmental delay and primary aldosteronism . So far this mutation has not been functionally characterized. Although its pathogenic potential appears to be high based on two different APA mutations reported in the same position (V259D, V259G) ), final proof requires functional analysis as described here.
Finally, another important and clinically highly relevant result of our study was the observation that Cav1.3 channels harboring the S652L-mutation require lower concentrations of the DHP isradipine for inhibition. This can be explained by the known voltage-dependence of DHP action due to their preferred binding to inactivated channel states [29, 42]. S652L induces a pronounced negative shift of the voltage-dependence of inactivation thus increasing the availability of inactivated channels. This preclinical finding is a strong motivation to test if repurposing of already available DHPs could ameliorate symptoms in affected individuals. Since also other mutations inducing strong shifts in steady-state inactivation are likely to increase DHP sensitivity (e.g. A749G, ), this treatment approach, if successful, could also be offered to individuals with other CACNA1D mutations.
Although our data strongly support CACNA1D as a high-risk gene for neurodevelopmental disorders and emphasize the need of functional analysis to distinguish likely pathogenic (able to increase Cav1.3 activity) from non-pathogenic de novo mutations (unable to increase Cav1.3 activity), our studies do not provide insight into altered signaling cascades downstream of Cav1.3 channels. This will require introduction of one or more of these human mutations into the mouse Cacna1d gene for electrophysiological and biochemical studies in native cells. Such animal models will also allow to address the important question, if currently available LTCC blocker, such as isradipine, felodipine of nimodipine, can normalize cellular function and, perhaps, even behavioral phenotypes in these mice. Although LTCC blockers may normalize the mutation-induced increase in channel function after diagnosis, it is possible that the mutation may have already caused permanent developmental deficiencies resistant to drug treatment. Therefore the clinical potential of this therapy needs to be tested in small clinical trials in affected individuals.
Taken together our data have important implications for genetic diagnostics. We provide evidence that CACNA1D is a neurodevelopmental disorder-linked gene. Although initially considered to cause high risk only for ASD with or without intellectual disability , the increasing number of affected individuals, including S652L, now strongly indicates that the majority presents with a more severe phenotype. This can involve seizures, intellectual disability and, due to the role of Cav1.3 for aldosterone and insulin secretion [40, 46]), also with (often transient, [17, 18]) endocrine symptoms. Our findings with S652W (the loss-of-function mutation) emphasize that, in the case of CACNA1D, the amino acid position itself does not allow predictions about the disease risk of a variant, even if bioinformatics prediction tools provide high scores for protein damage. A high probability for pathogenicity can also be assumed if a variant identical to the germline mutation has also been found in at least two different individuals as a somatic mutation in an APA or an aldosterone-producing cell-cluster [34, 47, 48]. Our report should raise awareness for the pathogenic potential of CACNA1D mutations, especially in patients without additional congenital endocrine symptoms as diagnostic features. At present, de novo CACNA1D missense mutations may be underdiagnosed in clinical practice.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Action potential-like waveform
Autism spectrum disorder
Central nervous system
- G :
Human embryonic kidney
- I Ba :
- I Ca :
- P O :
- Q ON :
“ON” gating charge
- V 0.5 :
Half-maximal activation voltage
- V max :
Potential of maximal inward current
- V rev :
Nanou E, Catterall WA. Calcium Channels, Synaptic Plasticity, and Neuropsychiatric Disease. Neuron. 2018;98(3):466–81.
Dolphin AC. Voltage-gated calcium channels and their auxiliary subunits: physiology and pathophysiology and pharmacology. J Physiol. 2016;594(19):5369–90.
Zamponi GW, Striessnig J, Koschak A, Dolphin AC. The Physiology, Pathology, and Pharmacology of Voltage-Gated Calcium Channels and Their Future Therapeutic Potential. Pharmacol Rev. 2015;67(4):821–70.
Striessnig J, Pinggera A, Kaur G, Bock G, Tuluc P. L-type Ca2+ channels in heart and brain. Wiley Interdiscip Rev Membr Transp Signal. 2014;3(2):15–38.
Cohen SM, Li B, Tsien RW, Ma H. Evolutionary and functional perspectives on signaling from neuronal surface to nucleus. Biochem Biophys Res Commun. 2015;460(1):88–99.
Alexander SP, Striessnig J, Kelly E, Marrion NV, Peters JA, Faccenda E, et al. THE CONCISE GUIDE TO PHARMACOLOGY 2017/18: Voltage-gated ion channels. Br J Pharmacol. 2017;174(Suppl 1):S160–s94.
Heyes S, Pratt WS, Rees E, Dahimene S, Ferron L, Owen MJ, et al. Genetic disruption of voltage-gated calcium channels in psychiatric and neurological disorders. Prog Neurobiol. 2015;134:36–54.
Moon AL, Haan N, Wilkinson LS, Thomas KL, Hall J. CACNA1C: Association With Psychiatric Disorders, Behavior, and Neurogenesis. Schizophr Bull. 2018;44(5):958–65.
Kabir ZD, Martinez-Rivera A, Rajadhyaksha AM. From Gene to Behavior: L-Type Calcium Channel Mechanisms Underlying Neuropsychiatric Symptoms. Neurotherapeutics. 2017;14(3):588–613.
Splawski I, Timothy KW, Sharpe LM, Decher N, Kumar P, Bloise R, et al. Cav1.2 calcium channel dysfunction causes a multisystem disorder including arrhythmia and autism. Cell. 2004;119(1):19–31.
Splawski I, Timothy KW, Decher N, Kumar P, Sachse FB, Beggs AH, et al. Severe arrhythmia disorder caused by cardiac L-type calcium channel mutations. Proc Natl Acad Sci U S A. 2005;102(23):8089–96 discussion 6-8.
Gillis J, Burashnikov E, Antzelevitch C, Blaser S, Gross G, Turner L, et al. Long QT, syndactyly, joint contractures, stroke and novel CACNA1C mutation: expanding the spectrum of Timothy syndrome. Am J Med Genet A. 2012;158a(1):182–7.
Atkinson LZ, Colbourne L, Smith A, Harmer CH, Nobre AC, Rendell J, et al. The Oxford study of Calcium channel Antagonism, Cognition, Mood instability and Sleep (OxCaMS): study protocol for a randomised controlled, experimental medicine study. Trials. 2019;20(1):120.
Pinggera A, Lieb A, Benedetti B, Lampert M, Monteleone S, Liedl KR, et al. CACNA1D de novo mutations in autism spectrum disorders activate Cav1.3 L-type calcium channels. Biol Psychiatry. 2015;77(9):816–22.
Pinggera A, Mackenroth L, Rump A, Schallner J, Beleggia F, Wollnik B, et al. New gain-of-function mutation shows CACNA1D as recurrently mutated gene in autism spectrum disorders and epilepsy. Hum Mol Genet. 2017;26(15):2923–32.
Pinggera A, Negro G, Tuluc P, Brown MJ, Lieb A, Striessnig J. Gating defects of disease-causing de novo mutations in Cav1.3 Ca2+ channels. Channels (Austin). 2018;12(1):388–402.
Scholl UI, Goh G, Stolting G, de Oliveira RC, Choi M, Overton JD, et al. Somatic and germline CACNA1D calcium channel mutations in aldosterone-producing adenomas and primary aldosteronism. Nat Genet. 2013;45(9):1050–4.
Flanagan SE, Vairo F, Johnson MB, Caswell R, Laver TW, Lango Allen H, et al. A CACNA1D mutation in a patient with persistent hyperinsulinaemic hypoglycaemia, heart defects, and severe hypotonia. Pediatr Diabetes. 2017;18(4):320–3.
Semenova NA, Ryzhkova OR, Strokova TV, Taran NN. The third case report a patient with primary aldosteronism, seizures, and neurologic abnormalities (PASNA) syndrome de novo variant mutations in the CACNA1D gene. Zh Nevrol Psikhiatr Im S S Korsakova. 2018;118(12):49–52.
Garza-Lopez E, Lopez JA, Hagen J, Sheffer R, Meiner V, Lee A. Role of a conserved glutamine in the function of voltage-gated Ca2+ channels revealed by a mutation in human CACNA1D. J Biol Chem. 2018;293(37):14444–54.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv. 2019:531210.
Iossifov I, Ronemus M, Levy D, Wang Z, Hakker I, Rosenbaum J, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74(2):285–99.
The Deciphering Developmental Disorders S. Large-scale discovery of novel genetic causes of developmental disorders. Nature. 2015;519(7542):223–8.
Baig SM, Koschak A, Lieb A, Gebhart M, Dafinger C, Nurnberg G, et al. Loss of Cav1.3 (CACNA1D) function in a human channelopathy with bradycardia and congenital deafness. Nat Neurosci. 2011;14(1):77–84.
Liaqat K, Schrauwen I, Raza SI, Lee K, Hussain S, Chakchouk I, et al. Identification of CACNA1D variants associated with sinoatrial node dysfunction and deafness in additional Pakistani families reveals a clinical significance. J Hum Genet. 2019;64(2):153–60.
Bock G, Gebhart M, Scharinger A, Jangsangthong W, Busquet P, Poggiani C, et al. Functional properties of a newly identified C-terminal splice variant of Cav1.3 L-type Ca2+ channels. J Biol Chem. 2011;286(49):42736–48.
Koschak A, Reimer D, Huber I, Grabner M, Glossmann H, Engel J, et al. alpha 1D (Cav1.3) subunits can form L-type Ca2+ channels activating at negative voltages. J Biol Chem. 2001;276(25):22100–6.
Lieb A, Ortner N, Striessnig J. C-terminal modulatory domain controls coupling of voltage-sensing to pore opening in Cav1.3 L-type Ca2+ channels. Biophys J. 2014;106(7):1467–75.
Ortner NJ, Bock G, Dougalis A, Kharitonova M, Duda J, Hess S, et al. Lower Affinity of Isradipine for L-Type Ca2+ Channels during Substantia Nigra Dopamine Neuron-Like Activity: Implications for Neuroprotection in Parkinson's Disease. J Neurosci. 2017;37(28):6761–77.
Wu J, Yan Z, Li Z, Qian X, Lu S, Dong M, et al. Structure of the voltage-gated calcium channel Cav1.1 at 3.6 A resolution. Nature. 2016;537(7619):191–6.
Min BJ, Kim N, Chung T, Kim OH, Nishimura G, Chung CY, et al. Whole-exome sequencing identifies mutations of KIF22 in spondyloepimetaphyseal dysplasia with joint laxity, leptodactylic type. Am J Hum Genet. 2011;89(6):760–6.
Ohsugi M, Adachi K, Horai R, Kakuta S, Sudo K, Kotaki H, et al. Kid-mediated chromosome compaction ensures proper nuclear envelope formation. Cell. 2008;132(5):771–82.
Fernandes-Rosa FL, Williams TA, Riester A, Steichen O, Beuschlein F, Boulkroun S, et al. Genetic spectrum and clinical correlates of somatic mutations in aldosterone-producing adenoma. Hypertension. 2014;64(2):354–61.
Pinggera A, Striessnig J. Cav 1.3 (CACNA1D) L-type Ca2+ channel dysfunction in CNS disorders. J Physiol. 2016;594(20):5839–49.
Scharinger A, Eckrich S, Vandael DH, Schonig K, Koschak A, Hecker D, et al. Cell-type-specific tuning of Cav1.3 Ca2+-channels by a C-terminal automodulatory domain. Front Cell Neurosci. 2015;9:309.
Barrett CF, Tsien RW. The Timothy syndrome mutation differentially affects voltage- and calcium-dependent inactivation of Cav1.2 L-type calcium channels. Proc Natl Acad Sci U S A. 2008;105(6):2157–62.
Olson PA, Tkatch T, Hernandez-Lopez S, Ulrich S, Ilijic E, Mugnaini E, et al. G-protein-coupled receptor modulation of striatal Cav1.3 L-type Ca2+ channels is dependent on a Shank-binding domain. J Neurosci. 2005;25(5):1050–62.
Limpitikul WB, Dick IE, Ben-Johny M, Yue DT. An autism-associated mutation in Cav1.3 channels has opposing effects on voltage- and Ca2+-dependent regulation. Sci Rep. 2016;6:27235.
Helton TD, Xu W, Lipscombe D. Neuronal L-type calcium channels open quickly and are inhibited slowly. J Neurosci. 2005;25(44):10247–51.
Sinnegger-Brauns MJ, Hetzenauer A, Huber IG, Renstrom E, Wietzorrek G, Berjukov S, et al. Isoform-specific regulation of mood behavior and pancreatic beta cell and cardiovascular function by L-type Ca2+ channels. J Clin Invest. 2004;113(10):1430–9.
Platzer J, Engel J, Schrott-Fischer A, Stephan K, Bova S, Chen H, et al. Congenital deafness and sinoatrial node dysfunction in mice lacking class D L-type Ca2+ channels. Cell. 2000;102(1):89–97.
Liao P, Yu D, Li G, Yong TF, Soon JL, Chua YL, et al. A smooth muscle Cav1.2 calcium channel splice variant underlies hyperpolarized window current and enhanced state-dependent inhibition by nifedipine. J Biol Chem. 2007;282(48):35133–42.
Ortner NJ, Striessnig J. L-type calcium channels as drug targets in CNS disorders. Channels (Austin, Tex). 2016;10(1):7–13.
Uchida S, Yamada S, Nagai K, Deguchi Y, Kimura R. Brain pharmacokinetics and in vivo receptor binding of 1,4-dihydropyridine calcium channel antagonists. Life Sci. 1997;61(21):2083–90.
Akerstrom T, Willenberg HS, Cupisti K, Ip J, Backman S, Moser A, et al. Novel somatic mutations and distinct molecular signature in aldosterone-producing adenomas. Endocr Relat Cancer. 2015;22(5):735–44.
Azizan EA, Poulsen H, Tuluc P, Zhou J, Clausen MV, Lieb A, et al. Somatic mutations in ATP1A1 and CACNA1D underlie a common subtype of adrenal hypertension. Nat Genet. 2013;45(9):1055–60.
Omata K, Anand SK, Hovelson DH, Liu C-J, Yamazaki Y, Nakamura Y, et al. Aldosterone-Producing Cell Clusters Frequently Harbor Somatic Mutations and Accumulate With Age in Normal Adrenals. J Endocrine Soc. 2017;1(7):787–99.
Nanba K, Omata K, Gomez-Sanchez CE, Stratakis CA, Demidowich AP, Suzuki M, et al. Genetic Characteristics of Aldosterone-Producing Adenomas in Blacks. Hypertension. 2019;73(4):885–92.
We thank Jennifer Müller, Gospava Stojanovic, Bettina Tschugg for expert technical assistance. The Deciphering Developmental Disorders study presents independent research commissioned by the Health Innovation Challenge Fund [grant number HICF-1009-003], a parallel funding partnership between Wellcome and the Department of Health, and the Wellcome Sanger Institute [grant number WT098051]. The views expressed in this publication are those of the author(s) and not necessarily those of Wellcome or the Department of Health. The research team acknowledges the support of the National Institute for Health Research, through the Comprehensive Clinical Research Network.
This work was supported by grants of the Austrian Science Fund (FWF, P27809, W1101, DOC 30/CavX).
Ethics approval and consent to participate
The Deciphering Developmental Disorders Study  has UK Research Ethics Committee approval (10/H0305/83, granted by the Cambridge South REC, and GEN/284/12, granted by the Republic of Ireland REC). Parental informed consent has been obtained for updated clinical information of the probands.
Consent for publication
Parental consent was obtained for anonymized data being included in publications.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Expression of WT and S652L Cav1.3 α1-subunits in HEK-293 cells. (A) Expression of C-terminally long and short WT and S652L α1-subunits by Western blot analysis. One representative Western blot of >3 experiments from >2 membrane preparations of transfected HEK-293 cells stably expressing β3 and α2δ-1 is shown. The apparent molecular mass of the full length forms of the long and short Cav1.3 α1-subunit splice variants obtained under our experimental conditions is indicated. (B) Quantification of relative total protein expression levels was carried out by integrating densities of WT and mutant signals and normalization to the loading control α-tubulin (α-tub). Statistics: unpaired student´s t-test compared to WT (S652 LL: 118.26 ± 23.43, n = 4; S652 LS: 129.62 ± 8.94, n = 3). Data are presented as mean ± SEM.
Relationship of QON-gating charge movement and integrated tail current amplitude (Itail) as an indirect estimate of open probability. (A) Channel open probability was estimated from the slope of the QON-Itail relationships measured at Vrev. Slopes were obtained by linear regression: Itail/QON [mean ± SEM; ms-1]: WTL: -7.22 ± 0.916, r2 = 0.72, n = 26; S652LL: -4.24 ± 0.657, r2 = 0.72, n = 25; slopes are significantly different: F = 6.43, p = 0.015, F test). (B) Representative traces obtained by depolarization to the reversal potential from a HP of -89 mV. Data were collected from more than three independent transfections.
About this article
Cite this article
Hofer, N.T., Tuluc, P., Ortner, N.J. et al. Biophysical classification of a CACNA1D de novo mutation as a high-risk mutation for a severe neurodevelopmental disorder. Molecular Autism 11, 4 (2020). https://doi.org/10.1186/s13229-019-0310-4