| |||||||||||||||||||||||||||||||||||||||||
We found that sinus bradycardia in members of a large family was associated with a mutation in the gene coding for the pacemaker HCN4 ion channel. Pacemaker channels of the sinoatrial node generate spontaneous activity and mediate cyclic AMP (cAMP)dependent autonomic modulation of the heart rate. The mutation associated with bradycardia is located near the cAMP-binding site; functional analysis found that mutant channels respond normally to cAMP but are activated at more negative voltages than are wild-type channels. These changes, which mimic those of mild vagal stimulation, slow the heart rate by decreasing the inward diastolic current. Thus, diminished function of pacemaker channels is linked to familial bradycardia.
The pacemaker ("funny") current (If) of sinoatrial-node myocytes determines the slope of the diastolic depolarization of pacemaker cells and thus has a key role in the generation and autonomic regulation of sinus rhythm and rate.2 Because of their specific involvement in pacemaking, f-channels are the target for the pharmacologic control of heart rate. Indeed, substances acting by specific f-channel blockade slow the heart rate without side effects, such as the negative inotropic effects typical of beta-blockers or calcium antagonists.3 This property is particularly useful in the treatment of coronary heart disease.4
The f-channels are encoded by the hyperpolarization-activated, cyclic nucleotidegated (HCN) channel gene family.5,6 Of the four known HCN subunits (1 through 4), HCN4 is the most highly expressed in the mammalian sinoatrial node.7 We therefore screened a panel of persons with asymptomatic or symptomatic bradycardia (including sick-sinus syndrome and atrioventricular block) for mutations in the hHCN4 gene. The hHCN4 gene is located on chromosome 15 (15q2425)8; the genes coding for other channels contributing to the electrical activity of sinoatrial-node cells, including the L- and T-type calcium channels and delayed-rectifying potassium channels, are all located on different chromosomes.9
We report here on a family with a hereditary form of asymptomatic bradycardia associated with a mutation in the pacemaker-channel
-subunit HCN4. Functional analysis revealed that the mutant channels are activated at voltages more negative than those at which wild-type channels are activated; as a result, they supply less current during diastolic depolarization, which in turn results in slowing of the heart rate. The mutation mimics the effect of moderate vagal stimulation.
Methods
Our research protocol was reviewed and approved by the institutional review board for the Department of Biomolecular Sciences and Biotechnology, University of Milan.
Sequence Analysis and Mutagenesis
All subjects gave written informed consent before undergoing genetic analysis. Screening for mutations was performed on genomic DNA samples extracted from whole blood, saliva, or both (Puragene, Gentra Systems). The primers were designed to amplify DNA fragments of 200 to 350 bp in order to screen all of the coding portion of hHCN4. The polymerase-chain-reaction (PCR) products were analyzed by single-strand conformation polymorphism (SSCP) analysis. We used the primers 5'ATGCCTCATCCTGAGTCCTG3' (F) and 5'CTCACCAATGCGGTCCAG3' (R) to amplify exon 7 by Pfu polymerase (Stratagene). Mutagenesis was identified by DNA sequencing (MWG Biotech). A control group of 373 healthy subjects was used to exclude DNA polymorphisms. The hHCN4 complementary DNA was inserted in the eukaryotic expression vector pcDNA 1.1 (Clontech Laboratories), and the point mutation (2016 C
A) was generated by a PCR overlap method. The oligonucleotide used to incorporate the mutation into hHCN4 was 5'ACAGCCAGAGTGAGGGCC3'.
Electrophysiological Methods
Wild-type and mutant human (h) HCN4 channel complementary DNA was transfected for transient functional expression into HEK293 cells with a plasmid containing green fluorescent protein, as described previously.10 From one to five days after transfection, the cells were dispersed by trypsin and plated on 35-mm plastic petri dishes. A dish was placed under the stage of an inverted microscope, and GFP-expressing cells were selected for patch-clamp analysis at room temperature (25° to 26°C). The cells were initially superfused with Tyrode's solution containing 140 mM sodium chloride, 5.4 mM potassium chloride, 1.8 mM calcium chloride, 5.5 mM D-glucose, and 5 mM HEPESsodium hydroxide buffer (pH 7.4). The pipettes used in the whole-cell patch-clamp experiment contained 10 mM sodium chloride, 130 mM potassium chloride, 1 mM egtazic acid (EGTA), 0.5 mM magnesium chloride, 2 mM ATP (sodium salt), 0.1 mM guanosine triphosphate (sodium salt), 5 mM phosphocreatine, and 5 mM HEPESpotassium hydroxide buffer (pH 7.2). The control extracellular solution in whole-cell experiments contained 110 mM sodium chloride, 30 mM potassium chloride, 1.8 mM calcium chloride, 0.5 mM magnesium chloride, and 5 mM HEPESsodium hydroxide buffer (pH 7.4); 1 mM barium chloride, 2 mM manganese chloride, 100 µM nickel chloride, and 20 µM nifedipine were added to improve dissection of the pacemaker current. The pipettes used in inside-out macropatch11 experiments contained 70 mM sodium chloride, 70 mM potassium chloride, 1.8 mM calcium chloride, 1 mM magnesium chloride, 1 mM barium chloride, 2 mM manganese chloride, and 5 mM HEPESsodium hydroxide buffer (pH 7.4); a solution containing 130 mM potassium aspartate, 10 mM sodium chloride, 2 mM calcium chloride, 5 mM EGTApotassium hydroxide, and 10 mM HEPESpotassium hydroxide buffer (pH 7.2; pCa [the negative log of the concentration of Ca2+ ions] 7) perfused the intracellular sides of the patches.
For the investigation of mutant channels, HEK293 cultures were split into two halves, one to express control channels and one to express mutant channels, and the two halves were treated by identical procedures; experiments comparing the properties of mutant and control channels were always performed on cells matched according to the day of culture.
The activation curves for hHCN4 currents recorded under whole-cell conditions were obtained by standard activation and deactivation protocols and analyzed by the Boltzmann equation, y=1÷{1+exp[(VV1/2)÷s]}, where y is the fractional activation, V is the voltage in millivolts, V1/2 is the half-activation voltage in millivolts, and s is the inverse slope factor in millivolts. Mean activation curves were obtained by fitting individual curves from each cell to the Boltzmann equation and averaging half-activation voltages and inverse slope factors. The activation curves in inside-out macropatches, the time constants of current activation and deactivation, and the shifts induced by cAMP were calculated as previously reported.12 The doseresponse curves of cAMP-induced shifts in inside-out macropatches were analyzed by the Hill equation, as follows: S÷Smax = 1÷ [1+(k1/2÷cAMP concentration)h], where S is the shift, k1/2 is the half-maximal concentration, and h is the Hill factor. Each patch was exposed to cAMP only once. The holding potential was 35 mV in all experiments.
Measurements of Heart Rate
Heart-rate measurements were performed during the daytime (between 10 a.m. and 6 p.m.) with the patient at rest. The rates for children 10 years of age or younger were normalized to adult rates by linear scaling according to published data (Table A1.46 in Macfarlane and Lawrie13). The mean adult rates used for normalization were calculated by averaging the rates for the age groups 18 through 29 years, 30 through 39 years, 40 through 49 years, and 50 years or more from Table A1.1 in Macfarlane and Lawrie.13 These rates (±SE) were 71.0±0.7 beats per minute for men and 76.7±0.2 beats per minute for women.
Statistical Analysis
Data were compared by the independent Student's t-test or two-way analysis of variance, and P values of 0.05 or less were considered to indicate statistical significance. Genetic-linkage analysis and lod-score calculations were performed with the MLINK program of the LINKAGE software package.14
Results
We screened 52 persons with bradycardia for mutations anywhere in the coding region of the pacemaker-channel gene hHCN4. This procedure led us initially to identify a missense mutation in exon 7, S672R, in one member of an Italian family (the proband) with asymptomatic sinus bradycardia (heart rate, 43 beats per minute). We then collected and examined DNA from a total of 27 members of the same family. PCRSSCP analysis showed that the S672R mutation cosegregates with the bradycardic phenotype, indicating an autosomal dominant pattern (Figure 1). The heart rate varied from 43 to 60 beats per minute in persons with the mutated gene (mean value for 15 persons, 52.2±1.4 beats per minute) and from 64 to 81 beats per minute in those with the wild-type gene (mean value for 12 persons, 73.2±1.6 beats per minute; P<0.001 by Student's t-test) (Figure 1 and Figure 2B). The presence of the S672R mutation was confirmed by electropherograms (Figure 2A). The same mutation was absent in 746 control chromosomes from unrelated persons of Italian origin.
|
|
=0, the range of
is 0 to 0.5, step 0.01), indicating tight linkage. Changing the penetrance over the range from 0.7 to 1.0 generated maximum lod scores (at
=0) in the range of 5.07 to 5.66. The HCN channels are six transmembrane-domain channels that are dually activated by voltage hyperpolarization and cAMP,2,5 and the native cardiac pacemaker f-channels are composed mainly of hHCN4 subunits.7 The cAMP activates the channels by direct binding11 and partial removal of an intrinsic inhibitory mechanism,15,16 resulting in a depolarizing shift of the open probability curve12); this shift underlies the modulation of the heart rate by autonomic neurotransmitters.17,18,19
Since the S672R mutation is located within the cyclic nucleotidebinding domain (CNBD) of hHCN4 (Figure 2C), it might affect cAMP binding and thus interfere with the autonomic modulation of heart rate. The fact that S672 is highly conserved throughout the animal kingdom (Figure 2C) supports the view that it has an essential functional role.
To evaluate the functional effect of the S672R mutation, we transfected wild-type and mutated hHCN4 complementary DNA into HEK293 cells (Figure 3A). The mutant S672R channels were expressed in HEK293 cells as efficiently as were wild-type channels. However, the S672R mutant channels were activated at more negative voltages than were wild-type channels. The half-activation voltages, as calculated from the Boltzmann equation, were 76.1±1.3 mV for the wild-type channels and 84.5±1.5 mV for the mutant channels (P=0.002 by Student's t-test), with a negative shift of 8.4 mV; inverse slope factors did not change significantly (10.2±0.8 and 11.2±0.8 mV for wild-type and mutant channels, respectively; P=0.41) (Figure 3A, left-hand panel). Moreover, the mutant channels were deactivated faster than were the wild-type channels, as shown by plotting mean activation and deactivation time-constant curves for wild-type and mutant channels (Figure 3A, right-hand panel). The deactivation time constants were significantly faster in mutant than in wild-type channels, according to two-way analysis of variance (P<0.001). These properties of the mutant channels are compatible with a reduced contribution of the pacemaker current to diastolic depolarization.
|
Since the S672R mutation does not interfere with cAMP-induced channel activation, the negative shift of the activation curve in Figure 3A probably reflects a constitutive new biophysical property of the channel caused by the mutation. This was confirmed by measuring activation curves in inside-out macropatches in the absence of cAMP: the V1/2 of the mutant channels was 94.4 mV, 9.1 mV more negative than that of the wild-type channels (103.5 mV), and s was 8.7 mV in both cases (Figure 3B, right-hand panel). Thus, S672R channels are constitutively activated at voltages about 8 to 9 mV more negative than are wild-type channels.
Finally, in order to produce channels heteromeric for wild-type and mutant
-subunits, such as those occurring in heterozygotes, we cotransfected identical amounts of wild-type and S672R complementary DNA into HEK293 cells and measured the resulting currents (Figure 3C). The properties of heteromeric wild-typemutant channels were intermediate between those of wild-type and homomeric S672R channels. The mean activation curve had a V1/2 of 81.0±1.6 mV, with a shift of 4.9 mV relative to the control (P=0.042 by Student's t test); s was 11.2±0.7mV (P=0.38) (Figure 3C, left-hand panel); the deactivation time constants were significantly different from those of wild-type channels according to two-way analysis of variance (P=0.008) (Figure 3C, right-hand panel).
Shifting of the If activation curve underlies physiologic frequency changes in pacemaker cells.2 Moderate vagal activity slows the heart rate by a muscarinic-induced negative shift of the If activation curve,18 a mechanism contributing to the basal slow rate associated with vagal tone19); this effect is accompanied by a marked acceleration of channel deactivation18 and is fully explained by a reduction of intracellular cAMP.11 The effects of the heterozygous S672R mutation on hHCN4 channels (a shift of 4.9 mV in the activation curve and acceleration of deactivation) are therefore analogous to those of a low dose (10 to 30 nM) of acetylcholine.19 This led us to speculate that the bradycardic rate slowing in persons heterozygous for the S672R mutation might be similar to the slowing induced in freely beating pacemaker myocytes by 10 to 30 nM acetylcholine. Published data indicate that the slowing in pacemaker rate caused by acetylcholine at these concentrations varies between 15 and 41 percent.19 Indeed, in the family investigated, the bradycardic mutant rates were on average 29 percent slower than nonbradycardic rates (Figure 1), a result compatible with the above concept.
Discussion
We report that familial sinus bradycardia is associated with a mutation in the HCN4 pacemaker-channel
-subunit. HCN4 is the major HCN isoform contributing to native f-channels in the sinoatrial node, the natural cardiac pacemaker region.7 Defective HCN4 channels have been found to be associated with disorders of cardiac rhythm, but previous reports either were based on a single patient21 or described a complex array of rhythm disturbances without clarifying which specific effect was mechanistically related to the mutation.22
The f-channels control pacemaker activity by generating the diastolic depolarization phase of the action potential, and they mediate the chronotropic action of autonomic neurotransmitters.2 They are activated by cAMP11); by increasing cAMP levels,
-adrenergic stimulation increases diastolic If and consequently steepens diastolic depolarization, thus causing an acceleration of the heart rate17); the opposite mechanism operates when muscarinic stimulation by acetylcholine decreases cAMP levels and If and slows the heart rate.18,19
The hHCN4 mutation associated with bradycardia, S672R, is located in the CNBD (Figure 2C), a region composed of several residues that affect cAMP binding.23 We found, however, that the S672R mutation does not affect cAMP-induced channel activation, but it does modify channel kinetics by shifting the current activation range to hyperpolarized voltages and slowing current deactivation; these changes mimic those caused by a low concentration (10 to 30 nM) of acetylcholine19 and lead to a reduced flow of inward current during diastolic depolarization and hence to slowing of the heart rate.
Finally, it is known that several mutations may alter the kinetics and cAMP dependence of HCN channels24,25); furthermore, HCN kinetics are modified by
-subunits.6 Thus, the mutation reported here may represent a specific case of a broader mechanism for sinus bradycardia based on constitutive inhibition of pacemaker channels.
Supported by grants from the Ministero dell'Istruzione, dell'Università e della Ricerca (FISR 02.00652.ST97, Cofin 2003052919, and FIRB RBNE01LNX7) to Dr. DiFrancesco.
No potential conflict of interest relevant to this article was reported.
We are indebted to Professor D. Noble (University of Oxford) for comments; to Drs. C. Viscomi, B. Terragni, A. Moroni, G. Consalez, and L. Gianfranceschi for assistance; to D. Lopez and G. Suffia for technical help; to Dr. B.U. Kaupp (University of Jülich) for supplying hHCN4; and to Dr. S. Priori (University of Pavia) for giving us access to control DNA.
Source Information
From the Department of Biomolecular Sciences and Biotechnology, Laboratory of Molecular Physiology and Neurobiology, University of Milan, Milan (R.M., M.B., D.D.); and the Department of Cardiology, Merate Hospital, Merate, Italy (T.G.-R.).
Address reprint requests to Dr. DiFrancesco at the Department of Biomolecular Sciences and Biotechnology, Laboratory of Molecular Physiology and Neurobiology, University of Milan, Via Celoria 26, 20133 Milan, Italy, or at dario.difrancesco{at}unimi.it.
References
| |||||||||||||||||||||||||||||||||||||||||
This article has been cited by other articles:
HOME | SUBSCRIBE | SEARCH | CURRENT ISSUE | PAST ISSUES | COLLECTIONS | PRIVACY | TERMS OF USE | HELP | beta.nejm.org Comments and questions? Please contact us. The New England Journal of Medicine is owned, published, and copyrighted © 2009 Massachusetts Medical Society. All rights reserved. |