| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Integrative Physiology |
From the Departments of Biomedical Engineering (C.M.R., W.L., J.H., J.C., S.A.W., I.R.E.) and Medicine (J.C., S.A.W.), Washington University, St Louis, Mo; Department of Internal Medicine (F.R.), University of Cincinnati College of Medicine, Ohio; and Brown Foundation Institute of Molecular Medicine (R.L., A.J.M.), The University of Texas Health Science Center at Houston, Texas Heart Institute at St. Lukes Episcopal Hospital, Tex.
Correspondence to Igor R. Efimov, Department of Biomedical Engineering, Washington University, St Louis, MO 63130. E-mail igor{at}wustl.edu
| Abstract |
|---|
|
|
|---|
Key Words: arrhythmia vulnerability hypertrophic cardiomyopathy transgenic rabbit
| Introduction |
|---|
|
|
|---|
HCM is characterized by cardiac hypertrophy, myocyte disarray, and interstitial fibrosis.4 Myocyte disarray is considered the pathological hallmark of HCM, whereas myocyte hypertrophy and interstitial fibrosis are common in many cardiovascular conditions. HCM is caused by mutations in sarcomeric proteins (except for phenocopy), with the MYH7 gene encoding ß-myosin heavy chain (MyHC) being the most frequent cause of HCM.5
We have generated a transgenic (TG) rabbit model by cardiac-restricted expression of the mutant ß-MyHC-Q403, known to cause HCM in humans.6 The ß-MyHC-Q403 rabbits recapitulate the phenotype of human HCM and exhibit cardiac hypertrophy, myocyte disarray, interstitial fibrosis, and SCD.7,8 The ß-MyHC-Q403 rabbits afford the opportunity to delineate mechanisms of cardiac arrhythmias in HCM, as the rabbit has long been an experimental model of choice for studies of basic cardiac electrophysiology because of its ability to sustain complex reentrant arrhythmias that faithfully reproduce those in larger species.9,10 The rabbit heart has also been found to be most similar to the human heart with respect to both the effective size and wave patterns during ventricular fibrillation, both of which are principal parameters for the investigation of arrhythmia mechanisms.11
To elucidate the underlying mechanisms contributing to arrhythmia vulnerability and SCD in HCM, we used multiple imaging modalities to assess electrophysiological, structural, and molecular characteristics from the hearts of ß-MyHC-Q403 rabbits and compared them with wild-type (WT) controls.
| Materials and Methods |
|---|
|
|
|---|
TG Rabbits
Generation and initial phenotypic characterization of the ß-MyHC-Q403 TG rabbits have been published.7,8,12 All protocols were approved by the appropriate Institutional Animal Care and Use Committees. We studied 6 TG and 6 WT littermate control rabbits, matched for age, sex, and body weight. These characteristics as well as selected echocardiographic parameters are shown in Table 1.
|
Echocardiography
Echocardiographic and tissue Doppler images were obtained and analyzed as described previously.12,13
Optical Mapping Experiments
The optical mapping system and Langendorff-perfused rabbit heart have been described previously.14 The entire anterior epicardial surface was mapped with a 16x16 photo-diode array (PDA) with a field-of-view of approximately 2.7x2.7 cm (see Figure 2). The excitation–contraction uncoupler 2,3-butanedione monoxime (BDM) (15 mmol/L; Fisher Scientific, Fair Lawn, NJ) was added to the perfusate to suppress motion artifacts in the optical recordings.
Experimental Protocol
The pacing site was located on the anterior left ventricular (LV) epicardium midway between apex and base. The pacing current was at least 4 times the pacing threshold. An S1-S2 protocol was performed with 20 S1s delivered at a cycle length of 300 ms followed by an S2 of decreasing duration until the S2 pacing stimulus no longer produced an action potential. The last S2 interval to capture was deemed the effective refractory period. The S1-S2 protocol was performed on all 12 hearts; however, a full action potential duration (APD) restitution protocol was only performed for 6 hearts (3 WT, 3 TG, 1 male/2 female for each group). The full APD restitution protocol consisted of the same S1-S2 protocol but included longer S2 intervals up to the interval at which the sinus rate interrupted the S2 interval.
A protocol to determine the vulnerable grid15 was then performed. After a train of 20 paced beats at a basic cycle length of 300 ms, a truncated exponential monophasic shock (10-ms duration, 80% tilt) was delivered at a specified coupling interval (CI) that was randomized within the vulnerable grid.15 Shock strengths (SS) ranged from
1 to 18 V/cm. The CI from the last paced beat ranged from 40 to 220 ms with 20-ms steps. The vulnerable grid was completed for both shock polarities. The applied electric field is referred to as RV– when the electrode near the RV was used as the cathode. The opposite electrode configuration is referred to as LV–.
A sustained arrhythmia was defined as lasting more than 6 shock-induced extra beats, consistent with Fabritz et al.15 Any sustained arrhythmias lasting more than several minutes were terminated with a rescue shock. The upper (ULV) and lower (LLV) limit of vulnerability were defined as the highest and lowest SS at which a sustained arrhythmia was induced, respectively. The vulnerable window (VW) was defined as the ULV minus the LLV. The longest and shortest extent of the vulnerable period (VP) were defined as the longest and shortest CIs that produced a sustained arrhythmia, and the VP was defined as the longest minus the shortest extent of VP. After completion of the grid, several grid points were again measured to confirm stability of the grid. These points consisted of a SS slightly larger than the ULV and shocks applied at CIs slightly longer and shorter than the longest and shortest VPs to confirm that sustained arrhythmias were not produced.
Diffusion-Weighted Magnetic Resonance Imaging
Immediately after optical mapping experiments, hearts were perfused with 3.7% formaldehyde, diffusion tensor MRI (DTMRI) was performed, and the myofiber helix angle (
) was calculated for an LV short-axis slice midway between apex and base as described previously.16,17
Immunohistochemistry
After MR imaging, hearts were transferred to 20% sucrose for 2 two days and the LV lateral wall was embedded in Tissue-Tek OCT compound, frozen, and cryosectioned parallel to the epicardium. Immunohistochemistry was performed with a commercially available anti–connexin (Cx)43 antibody (C6219; Sigma) used at a dilution of 1:400. Alexa fluor 555 goat anti-rabbit IgG (A-21428; Molecular Probes) was used for the secondary antibody at a dilution of 1:1000. Confocal imaging was performed with a x40 lens on a Nikon C1/80i confocal microscope.
Data Analysis
Optical Mapping Data Analysis
All signal analysis software was custom developed using MATLAB (The Mathworks Inc, Natick, Mass). This software allowed for automated calculation of activation, repolarization, APD, and calibrated transmembrane voltage in PDA channels having a signal-to-noise ratio of
5 dB. For calibration, a resting potential of –85 mV and an action potential amplitude of 100 mV were assumed.
APD80 was automatically calculated in all channels that met the specified signal-to-noise ratio criteria (
5 dB). Conduction delay in each channel was calculated as the time difference between the pacing stimulus and activation. The percentage APD (%APD) at which a shock was applied was calculated for each channel according to the following equation14: equation
|
|
The %APD was then averaged across all included PDA channels for each vulnerable grid point. The vulnerable grid, the upper and lower limits of VP, and the VP are all expressed as %APD.
Cx43 Quantification
Cx43 density was quantified from immunohistochemical images as previously described.18 Four randomly selected images were acquired from each of 2 sister sections at each tissue location (epi-, mid-, or endocardium) for a total of 8 images per location and 24 images per heart. The Cx43 density at each location was calculated as the average density of the 8 images.
Statistical Analysis
All data are presented as means±SD. To compare group data between TG and WT hearts, either a t test for parameters with equal variance or Kruskal–Wallis test for those with unequal variances were performed. Where appropriate, a 2-way ANOVA and Fishers protected t test were also performed. A value of P<0.05 was considered statistically significant. All nonlinear fitting was performed using the "nlinfit" function in MATLAB. The maximum number of iterations was set to 1000. The termination tolerance on the residual sum of squares and the estimated coefficients were both set to 1e–8.
| Results |
|---|
|
|
|---|
Vulnerable Grid
Large differences in ULV were observed between WT and TG hearts (Figure 1 and Tables 2 and 3
). Significant statistical differences were found for RV– shocks (WT, 7.75±2.06 V/cm; versus TG, 13.80±2.05 V/cm; P=0.003), LV– shocks (WT, 7.00±2.74 V/cm; versus TG, 12.80±2.28 V/cm; P=0.010), and when RV– and LV– shocks were pooled (WT, 7.38±2.28 V/cm; versus TG, 13.30±2.11 V/cm; P=3.2e–5). No statistical differences in ULV were observed between RV– and LV– shocks for either group or when WT and TG data were pooled. No statistical differences in LLV were found between groups. However, the VW was found to be significantly different between WT and TG hearts (WT, 6.25±2.62 V/cm; versus TG, 12.08±2.29 V/cm; P=0.002). Analysis of all parameters relating to the phase of shock application including the longest and shortest extent of VP as well as the VP did not reveal any statistical differences between WT and TG hearts (Table 2). Because of instability of the vulnerable grid, caused by ischemia or prolonged experiment duration, data from 2 TG hearts and 1 WT heart were not included in the analysis of vulnerable grid parameters. Therefore, n=5 for TG hearts (2 male/3 female) and n=4 for WT hearts (1 male/3 female) for grid parameters only.
|
|
|
Patterns of epicardial shock-induced virtual electrode polarization (VEP) and postshock activation were similar for WT and TG hearts. Figure 2 shows the response to an RV– 8 V/cm shock applied at 79% and 86% of the APD to a WT and TG heart, respectively. Figure 2A and 2E shows the anterior views of each heart with the PDA field of view indicated with a red square. From these images, obvious differences in heart size can be observed, with apparent hypertrophy in the TG heart compared with the WT heart. However, the shock-induced VEPs, as observed on the epicardial surface of the heart and shown in Figure 2B and 2F, are almost identical between WT and TG hearts. The postshock spread of activation is also very similar with activation spreading quickly from the RV to the LV then slowly reentering at the base of the RV in both hearts (Figure 2C and 2G). Noteworthy, this results in several shock-induced extra beats in the WT heart (Figure 2D), whereas in the TG heart, a sustained arrhythmia occurs (Figure 2H).
|
In addition, the maximum VEP (VEPmax) and minimum VEP (VEPmin) were calculated for all animals at shock end for an 8 V/cm shock applied at approximately 80% of APD80. The results are reported in Table I in the online data supplement. The magnitudes of polarization as well as VEPmax–VEPmin were comparable between groups. Although the %APD at which shocks were applied varied slightly between animals, there was also no significant difference in this parameter.
APD Restitution
Neither a linear nor biexponential fit of the APD restitution curve revealed significant differences in stability of electrical activation between TG and WT hearts (Figure 3; see the online data supplement for details). The slope of the linear fit was compared and found to be similar between groups (TG, 0.472±0.182; versus WT, 0.424±0.262; P=0.728). For the biexponential fit, the fitting parameter
1 was compared. The TG group was found to have a slightly faster but not statistically significant
1 compared with the WT group (0.281±0.068 versus 0.342±0.209 ms; P=0.827).
|
Several other parameters related to the APD were also investigated, including the mean APD80, the dispersion of APD80 (online data supplement, Equation 2), and the effective refractory period. No statistical differences were observed between groups for any of these calculations (supplemental Table II).
Conduction Velocity
Conduction velocities and the anisotropic ratio of conduction velocity (ratio=Cl/Ct) were calculated for all 12 hearts for basic pacing at a cycle length of 300 ms as well as for an S2 of 180 ms. No statistical differences or interactions were found between WT and TG hearts or between cycle lengths of 180 and 300 ms for any of the velocity parameters when a 2-way ANOVA was performed, indicating that the higher ULV observed in the TG hearts cannot be directly attributed to either slower or more anisotropic conduction (Figure 4; supplemental Table III).
|
Myocardial Fiber Rotation
Myocardial fiber rotation across the LV wall (helix angle:
) was quantified using DTMRI. A larger transmural fiber angle shift was found in TG hearts compared with WT. The total fiber angle shift for the entire wall thickness (Figure 5) was found to be 95.6±10.9° in TG hearts compared with 79.2±7.8° in WT hearts (P=0.039).
|
Immunohistochemistry
Disruption of either the total amount or distribution of the gap junction protein Cx43 has been shown to promote conduction abnormalities and facilitate arrhythmias.19 Therefore, we tested the hypothesis that the transmural distribution of Cx43 may be altered in the TG rabbits. Transmural expression of Cx43 was significantly altered in TG hearts compared with WT hearts (Figure 6). WT hearts showed a general trend of increased expression of endocardial Cx43 compared with the mid- and epicardium. Mid-myocardial Cx43 density was found to be significantly decreased compared with the endocardial density for WT hearts (mid-myocardial, 2.68±0.77%; versus endocardial, 5.31±2.46%; P=0.034). In TG hearts, however, the trend was somewhat reversed, with increased (but not statistically significant) mid-myocardial Cx43 density compared with the endocardium (mid-myocardial, 5.46±2.44%; versus endocardial, 3.92±2.48%; P=0.301). When WT and TG hearts were compared, a significant difference in Cx43 density was found in the mid-myocardium with TG hearts having an
2-fold increase in density compared with WT hearts (5.46±2.44% versus 2.68±0.77%; P=0.024). All Cx43 data are reported in supplemental Table IV.
|
| Discussion |
|---|
|
|
|---|
Both axes of the vulnerability grid (CI and SS) provide a measure of arrhythmogenicity, as the VP is closely related to the dispersion of repolarization15 and the ULV strongly correlates with the defibrillation threshold (DFT).20 An increase in either of these parameters indicates an increased susceptibility to shock-induced ventricular arrhythmias. In this study, similar values of VP were measured, indicating that WT and TG hearts have similar dispersion of repolarization. However, we found that the TG hearts had a much higher ULV and VW (almost 2-fold) compared with WT hearts; thus the TG hearts exhibit an overall increased susceptibility to postshock arrhythmias.
Interestingly, the epicardial postshock activation patterns and maps of VEP were remarkably similar in both groups (Figure 2), with the shock-induced phase singularity mechanism21 responsible for the initial reentrant activation observed on the anterior epicardial surface. The magnitude of VEP produced by the shock was also similar between groups, indicating that the driving force for postshock activation was also similar. However, the major difference between WT and TG hearts appears to be that in TG hearts, shocks below the ULV produced reentrant activation which resulted in a sustained (>6 extra beats up to several minutes) ventricular arrhythmia. However, shocks of the same magnitude and CI that produced sustained arrhythmias in TG hearts produced the same reentrant activation in WT hearts but often resulted in only several extra beats. Thus, it is possible that the mechanisms of arrhythmia "initiation" are not different in TG hearts compared with control, as the same virtual electrode-induced phase singularity mechanism21 is responsible. Rather, it is the mechanisms of arrhythmia "maintenance" that appear to result in a higher ULV.
The increased ULV also indicates that the DFT for TG hearts should be considerably higher20 than WT. Although the DFT may be increased in any hypertrophic state,22 a distinct increase in DFT has yet to be demonstrated in patients with HCM except for a small subset of patients with advanced hypertrophy.23 Because implantable cardioverter defibrillators are now being routinely recommended for both primary and secondary prevention of SCD in HCM,2,3 the possibility of increased DFTs in patients with moderate hypertrophy needs to be considered.
DTMRI revealed a modest yet statistically significant increase in transmural fiber rotation for a mid-ventricle short axis slice in TG hearts compared with WT. The increased transmural fiber rotation in TG hearts may play a role in converting nonsustained ventricular arrhythmias into sustained arrhythmias, thus promoting arrhythmia maintenance. Several numerical studies support this theory24–27 and have shown that an increased amount of transmural fiber rotation can result in scroll wave filament instability, promoting filament bending and scroll wave breakup. However, there is some disagreement regarding the role of a steep APD restitution curve in this phenomenon. Using the Luo–Rudy action potential model, Qu et al26 found that increased fiber rotation led to filament bending and scroll wave breakup only in the settings of a sufficiently steep APD restitution curve. Rappel27 also found that scroll wave breakup did not occur in the Luo–Rudy model but did occur when the Beeler–Reuter model was used. Using a simplified cardiac model, Fenton and Karma24,25 found that scroll wave breakup is a function of fiber rotation rate and the scroll wave trajectory and not of the slope of the APD restitution curve. Regardless of the exact role of the APD restitution curve, these numerical studies indicate that transmural fiber rotation may play a role in arrhythmia maintenance, thus supporting the experimental results obtained in the present study.
Abnormal expression of Cx43 can promote conduction abnormalities and may provide the substrate for arrhythmias.19 In the present study, a statistically significant increase in mid-myocardial Cx43 density was found in TG hearts compared with WT. This altered distribution may affect transmural conduction and facilitate arrhythmia maintenance without affecting epicardial conduction patterns recorded with optical mapping techniques because of the limited depth penetration of optical mapping.
Although total expression of Cx43 is known to be downregulated in the ventricle of patients with hypertrophic, dilated, and ischemic cardiomyopathies,28,29 we observed an increase in mid-myocardial Cx43 in the TG rabbits. Similarly, an increase in total Cx43 was observed in TG rabbits with a mutation in cardiac troponin I that causes HCM in humans.30 These findings may reflect different stages of remodeling or intrinsic differences that take place in familial hypertrophic cardiomyopathies compared with those that are in response to injury.
Most importantly, the increased Cx43 density found in the mid-myocardium of the TG hearts may provide an interesting link between the possibly altered mechanical stress experienced in the mid-myocardial layers of the hypertrophic hearts and the transcription and translation of Cx43. Finite element modeling studies have indicated that muscle fiber stress and strain across the LV wall are very sensitive to the transmural distribution of helix angle,31 changing up to 50%, with only slight variations in helix angle distribution.32 Our DTMRI results show differences in both the distribution and total amount of transmural fiber angle rotation; therefore, the LV wall mechanics may vary greatly between WT and TG hearts. In addition, it has been reported that pulsatile stretch leads to a dramatic upregulation of Cx43 in cultured myocytes.33 Thus, we conclude that altered ventricular mechanics may result in altered Cx43 expression. However, support for this hypothesis must first come in the form of measuring or predicting mechanical strain and/or stress in the hypertrophic ventricle via tagged MRI or numerical simulations.
In summary, we found a 2-fold increase in ULV and VW in the ß-MyHC-Q403 TG rabbits, suggesting an increase in both arrhythmia propensity and DFT. The increased ULV was associated with increased transmural fiber rotation and increased mid-myocardial Cx43 expression, yet the traditional electrophysiological parameters normally associated with arrhythmia were not present. Thus, structural, not electrophysiological, remodeling may underlie the increased arrhythmia propensity in HCM. This finding may have important implications regarding pharmacologic suppression of arrhythmias in HCM, as treatments aimed toward reducing or reversing the structural remodeling may be more effective than traditional antiarrhythmic drugs. Additionally, our findings indicate that an elevated DFT may be present even in the absence of extreme hypertrophy. This important finding should be taken into account as an ever-increasing number of HCM patients are receiving implantable cardioverter defibrillators.
Study Limitations
The limitations of the isolated rabbit heart have been described previously21 and include the use of BDM as an excitation–contraction uncoupler. It is well known that BDM can flatten the restitution curve and shorten APD, causing an offset in the quantitative results obtained in this study. However, because BDM was used in both groups of rabbits, it is very likely that the qualitative observations will still exist if this study were repeated without BDM. To optically map as much of the ventricular epicardium as possible, a tradeoff was made with regard to spatial resolution. Therefore, the fluorescent optical signals at each PDA location are a summary of approximately 1.7x1.7 mm2 of cardiac tissue. This summary effect may mask subtle differences between the WT and TG activation maps, conduction velocities, and VEP maps. In addition, the optical mapping technique cannot penetrate the epicardial surface beyond
1 mm. Therefore, all fluorescence recordings are from the epicardium only. Thus, the change in mid-myocardial Cx43 expression found in this study is not reflected in epicardial conduction velocity measurements. Finally, immunohistochemistry and MRI were performed on the same hearts as those used for the optical mapping study. During the mapping portion of the study, the isolated heart preparation is subject to protein degradation and edema, which can affect immunostaining and structural imaging results. However, all WT and TG hearts were subject to the same optical mapping protocol.
| Acknowledgments |
|---|
This work was supported by NIH grants R01-HL67322 and R01-HL074283 (to I.R.E.) and R01-HL68884 (to A.J.M.); an American Heart Association Predoctoral Fellowship (C.M.R.); Clinician–Scientist Award in Translation Research from Burroughs Wellcome Fund no. 1005907 (to A.J.M.); and the Greater Houston Community Foundation (TexGen) (A.J.M.).
Disclosures
None.
| Footnotes |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. G. Akar The Perfect Storm: Defective Calcium Cycling in Insulated Fibers With Reduced Repolarization Reserve Circ. Res., November 9, 2007; 101(10): 968 - 970. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Research Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |