Augmented currents of an HCN2 variant in patients with febrile seizure syndromes

The genetic architecture of common epilepsies is largely unknown. HCNs are excellent epilepsy candidate genes because of their fundamental neurophysiological roles. Screening in subjects with febrile seizures and genetic epilepsy with febrile seizures plus revealed that 2.4% carried a common triple proline deletion (delPPP) in HCN2 that was seen in only 0.2% of blood bank controls. Currents generated by mutant HCN2 channels were ∼35% larger than those of controls; an effect revealed using automated electrophysiology and an appropriately powered sample size. This is the first association of HCN2 and familial epilepsy, demonstrating gain of function of HCN2 current as a potential contributor to polygenic epilepsy. ANN NEUROL 2010;67:542–546

Hyperpolarization-activated cyclic nucleotide-gated ion channels (HCN) conduct I h important for neuronal pacemaker function.][9][10][11] HCN1 and HCN2 variants have been identified, but functional analyses have failed to determine a statistical difference in channel properties. 12The functional changes associated with gene variants that contribute to polygenic disease are, by definition, subtle.Both biological and experimental variability contribute to the variance seen in functional analyses, with the probability of failing to detect a real difference in 2 populations rising dramatically as variation increases.To separate such populations, the sample size needs to increase. 13Here, we used a candidate gene approach to search for HCN1 and HCN2 variants in epilepsy subjects and a medium-throughput automated electrophysiological assay to test channel function with appropriate sample sizes.Currents generated by mutant HCN2 channels were ~35% larger than controls, an increase that could enhance neuronal excitability. 8,9

Patient Collection
Diagnostic criteria for IGE and FS followed that of the Commission on Classification and Terminology of the International League Against Epilepsy, 14 and for GEFS+ followed that of Scheffer and Berkovic. 3Patients were Australian subjects of Caucasian origin and were screened for variants in a randomized fashion where clinical characteristics were blinded.Controls were randomly drawn anonymous Australian blood donors primarily of Caucasian origin.

DNA Preparation and Mutation Analysis
DNA was extracted from peripheral blood using the QIAamp DNA Blood Maxi Kit (Qiagen, Valencia, Germany).Hexosaminidase-labeled intronic primers flanking each exon were used to polymerase chain reaction (PCR)-amplify products between 250bp to 320bp, and the products were analyzed by single-strand conformation analysis (SSCA) on a realtime gel system using a Gel-Scan 2000 DNA fragment analyzer (Corbett Research, Mortlake, Australia).Products showing band shifts were sequenced using the BigDye Terminator Cycle Sequencing Ready Reaction kit (PE Applied Biosystems, Foster City, CA, v2.0), and the sequences were analyzed on an Applied Biosystems ABI Prism 3700 DNA Analyzer.
The HCN1 open reading frame (accession number NM_021072) contains 8 exons, which were divided into 13 amplicons for SSCA.The human HCN2 open reading frame (accession number NM_001194) contains 8 exons, which were divided into 14 amplicons for SSCA.We used single nucleotide polymorphisms (SNPs) to distinguish the genomic sequence spanning the HCN2 gene on chromosome 19 from the highly related chromosome 15 sequence.HCN1 is located on chromosome 5. Primer sequences used for PCR and SSCA can be found in Supporting Information.

HCN Mutagenesis and In Vitro Transcription
See Supporting Information.

HCN 2-Electrode Voltage Clamp Analysis and Statistical Design
Oocytes from Xenopus laevis were prepared as previously described. 15Thirty-five nanoliters of cRNA-encoding the wild-type (WT) and delPPP HCN2 subunits (65ng/μl; stocks confirmed spectrophotometrically and by gel analysis) was injected into stage 5/6 oocytes using the Roboocyte (Multi Channel Systems, Reutlingen, Germany) and stored for 2 days prior to experimentation.For voltage clamp recordings, oocytes were impaled with 2 glass electrodes containing 1.5M potassium acetate (I) and 0.5M KCl (V) and clamped at a holding potential of −30mV.A current-voltage (I-V) relationship was generated by incrementing voltage in 10mV steps from −140mV to 0mV for 15 seconds with a 2-second test potential of −140mV at the end of the pulse.Oocytes were perfused with a bath solution (mM): 96 KCl, 2 NaCl, 2 MgCl 2 , and 10 HEPES (pH 7.5 using KOH).To obtain normalized I-V relationships, peak tail current amplitudes were recorded at −140mV and were divided by the largest peak tail current.This normalized current was plotted against voltage and fit with a Boltzmann curve (GraphPad Prism, GraphPad Software, La Jolla, CA [average fits] and AxoGraph X, AxoGraph Scientific, Sydney, Australia [individual fits).Cyclic adenosine monophosphate (cAMP) modulation was investigated by incubating oocytes in 15μM forskolin (Sigma, Castle Hill, Australia) for a period of 7 minutes.Kinetics of activation was determined by measuring the time to half maximal current using a custom analysis program run in MatLab.Power analysis was performed using the G* Power calculator (www.psycho.uni-duesseldorf.de/aap/projects/gpower).Statistical comparisons were made using an unpaired t test (GraphPad Prism).

Variation of HCN1 and HCN2 in Epilepsy Patients
No major-effect sequence variation was detected in human HCN1 (Table 1).Analysis of HCN2 revealed a number of synonymous SNPs (Table 1) and a variant, c. 2156-2164delCGCCGCCGC, p.719-721PPP, predicted to lead to the deletion of 3 consecutive proline residues (delPPP) in the HCN2 protein.delPPP was found to be heterozygous in 3/65 unrelated patients (allele frequency = 2.3%) with GEFS+ (OMIM #604233), 3/61 unrelated patients (allele frequency = 2.5%) with FS, 3/772 blood bank controls (allele frequency = 0.2%), and 0/72 patients with classical IGE.The 3 FS patients had simple FS with a mean onset of 2 years.Frequency varied in each case, with 1, 3, or multiple events.The 3 GEFS+ patients had myoclonic-astatic epilepsy, FS+ or FS (this patient also carried the SCN1B[C121W] mutation 16 ).Clinical characteristics of patients who were negative for delPPP are described in Supporting Information.

Functional Analysis of delPPP HCN2 Channels
For current magnitude analysis, the inherent variance meant that sample sizes of 180 per group were required to detect differences of 25% with a power of 0.95 to reduce type 2 errors.Although all available data were used, smaller sample sizes were sufficient for all other parameters, because current magnitude had the highest standard deviation.
To achieve these high sample sizes, a medium-throughput electrophysiological assay was used.WT and delPPP HCN2 currents were elicited at hyperpolarizing potentials (Fig) .A significant ~35% increase in tail current was observed for the delPPP HCN2 channel compared with WT (5.32 ± 0.29μA vs 3.89 ± 0.23μA; n = 189 and 181, p = 0.001).Western blot analysis showed that the total level of WT and delPPP HCN2 channel expression were the same (Supporting Information).Normalized current-voltage relationships were also constructed from tail currents.Half maximal activation voltage (V 1/2 ), calculated by fitting a Boltzmann equation to individual activation curves, was identical in the WT and delPPP HCN2 channels (−69.9 ± 0.4mV vs −70.4 ± 0.4mV; p = 0.4).However, the slopes of the Boltzmann curve were significantly different for WT and delPPP HCN2 channels (4.70 ± 0.07 vs 5.18 ± 0.11, p = 0.0003).Activation kinetics were not significantly different across a range of activating voltages (p > 0.35).

cAMP Gating in delPPP HCN2 Channels
Because the delPPP variant is close to the cyclic nucleotide-binding domain, we investigated whether cAMP gating was altered by comparing channel sensitivity to forskolin, an activator of adenylate cyclase.V 1/2 standard deviation was approximately 11% of the mean, with power calculation suggesting that a minimum 10% difference could be detected at power 0.95 with a sample size of 50 per group.This sampling regime would be sensitive to a minimum 2.5% change in slope factor.There was no difference in V 1/2 in forskolin for the WT and delPPP HCN2 channels (−64.2 ± 1.1mV vs −64.4 ± 1.0mV, n = 40 and 55, p = 0.6, see Fig 1F).The slopes of the activation curve in forskolin were also not significantly different (4.05 ± 0.08 vs 4.23 ± 0.13, p = 0.28).

Discussion
Our comparison of HCN1 and HCN2 variation in epilepsy patients and controls revealed little genetic heterogeneity overall, suggesting an intolerance of sequence changes.However, a deletion of 3 consecutive prolines (delPPP) in HCN2 was identified in patients with GEFS+ and FS.The delPPP occurs in a 7-proline repeat in humans and is conserved as a 6-7 repeat in cows, chimpanzees, dogs, and mice, suggesting functional importance.The occurrence of the delPPP allele in controls is not unexpected, because susceptibility alleles are predicted to be present at low frequency in the general population. 2 We were able to detect the difference in the current between WT and delPPP channels by taking advantage of an automated electrophysiology system to record large numbers of oocytes.This enabled statistical detection of a current magnitude increase for the mutant channel.The need for large sample sizes is likely to become more common as an ever-increasing number of putative susceptibility alleles are identified.
The HCN2 delPPP variant occurs with highest frequency in patients with GEFS+ and FS, which both have febrile seizures at presentation.Increases in I h occur following induced febrile seizures in animal models and are proposed to contribute to hippocampal hyperexcitability.8,9 As we did not observe HCN2 delPPP in patients with IGE, which do not present with FS, this suggests that the variant may be a specific susceptibility allele for FS.The impact of changes in HCN function on neuronal excitability is multifaceted, with I h contributing to both resting membrane potential and input resistance. 17Increases in I h , as predicted here, will depolarize membrane potential, taking the neuron closer to the firing potential, and in this way be considered proexcitatory.A recent study that included computer simulation modeling supports this view. 9It is important to note that reduced HCN channel function is also thought to increase neuronal network excitability. 5,6Our kinetic analysis also isolated a very subtle difference in the slope of activation of the delPPP HCN2 channel in comparison to WT, the functional significance of which is unclear.
Simple expression systems lack sufficient complexity to reveal changes that may be neuron specific (eg, subcellular expression).Further, investigations of changes inchaperone protein interactions (eg, TRIP8b 18 ) and a host of other second messengers known to alter HCN function (eg, p38 mitogen-activated protein kinase 19 ) need to be considered.Ultimately, in vivo modeling that allows replication of true complex genetics with the introduction of multiple susceptibility alleles will determine the behavioral impact of the HCN2 delPPP variant. 20

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.

FIGURE .
FIGURE.Electrophysiological characterization of triple proline deletion (delPPP) HCN2 channels.(A) Steady state (top) and tail (bottom) currents from oocytes expressing wild-type (WT) and (B) delPPP HCN2 channels.(C) Average tail current-voltage relationship of WT and delPPP HCN2.*p < 0.05 (D) Activation curve constructed from average normalized current relationships constructed from tail currents fit with the Boltzmann equation.(E) Time to half maximal activation for WT and delPPP HCN2 channels across a range of voltages.(F) Modulation of the voltage dependence of activation of WT and delPPP HCN2 channels by forskolin (forsk), data from D replotted for comparison.