Abstract
Cardiovascular mortality is high in ESRD, partly driven by sudden cardiac death and recurrent heart failure due to uremic cardiomyopathy. We investigated whether speckle-tracking echocardiography is superior to routine echocardiography in early detection of uremic cardiomyopathy in animal models and whether it predicts cardiovascular mortality in patients undergoing dialysis. Using speckle-tracking echocardiography in two rat models of uremic cardiomyopathy soon (4–6 weeks) after induction of kidney disease, we observed that global radial and circumferential strain parameters decreased significantly in both models compared with controls, whereas standard echocardiographic readouts, including fractional shortening and cardiac output, remained unchanged. Furthermore, strain parameters showed better correlations with histologic hallmarks of uremic cardiomyopathy. We then assessed echocardiographic and clinical characteristics in 171 dialysis patients. During the 2.5-year follow-up period, ejection fraction and various strain parameters were significant risk factors for cardiovascular mortality (primary end point) in a multivariate Cox model (ejection fraction hazard ratio [HR], 0.97 [95% confidence interval (95% CI), 0.95 to 0.99; P=0.012]; peak global longitudinal strain HR, 1.17 [95% CI, 1.07 to 1.28; P<0.001]; peak systolic and late diastolic longitudinal strain rates HRs, 4.7 [95% CI, 1.23 to 17.64; P=0.023] and 0.25 [95% CI, 0.08 to 0.79; P=0.02], respectively). Multivariate Cox regression analysis revealed circumferential early diastolic strain rate, among others, as an independent risk factor for all-cause mortality (secondary end point; HR, 0.43; 95% CI, 0.25 to 0.74; P=0.002). Together, these data support speckle tracking as a postprocessing echocardiographic technique to detect uremic cardiomyopathy and predict cardiovascular mortality in ESRD.
Patients with ESRD exhibit a dramatically increased risk for cardiovascular morbidity and mortality.1–3 Interestingly, this fact is not mainly driven by the traditional atherosclerotic risk factors of the general population.4 In contrast to the general population, where myocardial infarction and stroke are the predominant causes of cardiovascular death, sudden cardiac death (SCD) and recurrent heart failure are the most common causes of cardiovascular mortality in this patient group.5–7 This might be caused by uremic cardiomyopathy (UC), a common complication in patients with ESRD that is characterized by cardiac fibrosis, capillary rarefaction, left ventricular hypertrophy, and both systolic and diastolic dysfunction.8,9 In fact, two of these characteristics—cardiac fibrosis and capillary rarefaction—trigger electric instability of the myocardium and thus may be directly related to SCD in patients with ESRD.5 Therefore, early detection of UC in patients with ESRD might identify patients at risk for SCD and might have clinical impact on assessment and therapeutic decision-making in those patients.
Two-dimensional speckle-tracking echocardiography (STE) imaging is a relatively novel echocardiographic approach to assess regional left ventricular function by tracking acoustic markers (speckles) within the myocardium from frame to frame in B-mode images.10,11 Many studies used this technique in comparison with conventional echocardiographic measures. Most conventional echocardiographic measures determine myocardial motion; the motion of any part of the myocardium is influenced by translational and tethering effects. STE analysis determines deformation and subtracts the effect of tethering and translational motion of the whole heart and is therefore superior to other echocardiographic measures regarding early changes of myocardial function. This was shown even in the setting of special patient cohorts.12–15 Because it does not depend on the Doppler angle, as in tissue Doppler imaging, it allows assessment of strain and strain rate analyses in different myocardial axes (i.e., longitudinal, circumferential, and radial). STE has been used to gain further insight into the pathophysiology of cardiac ischemia and infarction, and several primary diseases of the myocardium, such as hypertrophic or diabetic cardiomyopathy.16–18 It has recently been appreciated that STE can also be used to identify areas of myocardial fibrosis,19–24 one of the hallmarks of UC. This prompted us to evaluate in an animal study whether STE can detect early changes of left ventricular function and correlates with histologically determined severity of myocardial fibrosis in rat models of UC. Furthermore, in a clinical study we determined whether STE can predict cardiovascular and all-cause mortality in patients with ESRD.
Results
Animal Study
Uremic Cardiomyopathy Was Characterized by Cardiac Hypertrophy and Interstitial Myocardial Fibrosis in Both Rat Models of CKD
To determine whether echocardiographically measured strain and strain rate values can detect early myocardial changes and myocardial fibrosis, we performed an animal study using STE combined with histological quantification of fibrosis in two rat models of CKD with uremic cardiomyopathy (i.e., adenine-nephropathy and 5/6 nephrectomy; Figure 1). After induction of CKD, serum creatinine, urea, and phosphorus were significantly increased, whereas creatinine clearance was significantly decreased in both models of UC compared with the control group (Table 1). CKD was accompanied by significant myocardial hypertrophy (significantly increased cardiomyocyte cross-sectional area [Figure 2, A and B]) and increased heart weight (Figure 2C). Both heart weight and cardiomyocyte cross-sectional area showed a significant linear regression with serum phosphate (Figure 2, D and E). Quantitative real-time PCR of whole-heart mRNA expression revealed increased expression of β-myosin heavy chain and brain natriuretic peptide in both rat models of UC at the early and late time points (Figure 2, F and G).
Experimental scheme of the animal study. High-phosphorus diet was given after subtotal nephrectomy (5/6 Nx).
Humoral data from rats stratified for experimental groups and time points
The severity of cardiac hypertrophy correlates with serum phosphate levels. (A and B) Wheat germ agglutinin (WGA) staining (A) and quantification of cardiomyocyte cross-sectional area (B) indicates hypertrophy of cardiomyocytes in both models of UC: early (i.e., 4–6 weeks after induction of CKD) and late (8–10 weeks after induction of CKD). (C) Significantly increased heart weight in both models of UC at early and late time points. (D and E) Parameters of left ventricular hypertrophy show a significant linear regression with serum phosphate levels. (F and G) Quantitative real-time PCR analysis revealed significant cardiac upregulation of β-myosin heavy chain (MHC) (D) and brain natriuretic peptide (BNP) in both models of UC. *P<0.05; **P<0.01; ***P<0.001 (one-way ANOVA with post hoc Scheffe test). All scale bars=50 μm.
Quantification of Masson trichrome–stained sections revealed a significant increase of interstitial extracellular matrix in both models of UC (Figure 3, A and B). The interstitial myocardial fibrosis was further characterized by increased expression of extracellular matrix proteins collagen 1α1, collagen 3α1, and fibronectin (Figure 3, C–F, Supplemental Figure 1) and upregulation of TGF-β mRNA transcripts (Figure 3G). Serum phosphate levels correlated significantly with the degree of interstitial myocardial fibrosis (Figure 3H).
Uremic cardiomyopathy is characterized by interstitial myocardial fibrosis with upregulation of TGF-β expression. (A and B) Trichrome staining of hearts (A) and quantification (B) demonstrates significantly increased interstitial myocardial fibrosis in both models of UC at early and late time points. (C and D) Staining (C) and quantitative real-time PCR (D) for collagen 1α1 revealed increased expression in both rat models of UC. (E and F) Other extracellular matrix proteins, such as collagen 3α1 (E) and fibronectin (F), also showed increased expression of mRNA and protein levels (immunostaining: Supplemental Figure 1). (G) Quantitative real-time PCR showed increase expression of TGF-β mRNA transcript in both rat models of UC at early and late time points. (H) The severity of interstitial myocardial fibrosis did correlate significantly with serum phosphate levels. *P<0.05; **P<0.01; ***P<0.001 (one-way ANOVA with post hoc Scheffe test). All scale bars=50 μm.
STE Detected Early Impairment of Left Ventricular Systolic and Diastolic Deformation in Experimental Uremic Cardiomyopathy
Myocardial hypertrophy in both models of UC was confirmed by significantly increased thickness of the interventricular septum and the posterior wall as assessed by echocardiography (Table 2). Whereas standard echocardiographic measures, such as fractional shortening (FS), left ventricular stroke volume, and cardiac output (CO) were altered only at late time points (Figure 4, A and B, Table 2), left ventricular deformation parameters assessed by STE as peak global and systolic strain radial and circumferential (SR and SC peak G and peak S, respectively) were already significantly reduced at early time points (Figure 4, A and B, Table 2). Furthermore, strain rate parameters (i.e., myocardial deformation per time unit) were significantly reduced at early and late time points in the radial orientation at both the systolic peak (SrR peak systolic) and the early diastolic peak (SrR diastolic peak E) (Figure 4, C and D). Thus, STE detected an early global (systolic and diastolic) impairment of left ventricular deformation in UC, which was not reflected in reduced conventional echocardiographic readouts of left ventricular function.
Echocardiographic standard and speckle-tracking deformation parameters of rats stratified for experimental groups
Uremic cardiomyopathy leads to early impairment of left ventricular systolic and diastolic deformation. (A) At the early time points (i.e., 4–6 weeks after induction of CKD), both peak global radial and circumferential strain (SR peak G, SC peak G) were already significantly reduced, whereas standard-echocardiographic parameters, such as FS and Duplex ultrasonographically determined CO, remained unchanged. Representative parasternal short-axis views with color-coded radial and circumferential strain show a dramatic reduction in predominantly radial strain. (B) At late time points (8–10 weeks after induction of CKD), in addition to the myocardial deformation parameters (SR peak G and SC peak G), conventional echocardiographic readouts of left ventricular function FS and CO show a significant decrease. (C–E) Strain rate parameters indicating the rate by which the myocardial deformation occurs (deformation/ time unit) were significantly reduced at early and late time points in the radial orientation at both the systolic peak (SrR peak systolic [C]) and the early diastolic peak (SrR diastolic peak E [E]) in both rat models of UC. In contrast, the peak systolic rate of myocardial deformation in the circumferential direction (SrC peak systolic) was significantly reduced only in the 5/6 nephrectomy (Nx) model of UC. *P<0.05; **P<0.01; ***P<0.001 versus control group.
Myocardial Deformation Parameters Correlated Significantly with the Degree of Interstitial Myocardial Fibrosis and Myocardial Hypertrophy in Rats
We further asked whether the myocardial deformation parameters assessed by STE correlate with characteristics of UC. Several of the determined deformation parameters showed a higher correlation with the degree of interstitial fibrosis, the myocardial weight, and the cardiomyocyte cross-sectional area than the commonly used parameters of left ventricular function, such as FS or CO (Figure 5, Supplemental Table 1). Deformation parameters of the systolic and global left ventricular function as peak global radial and circumferential strain (Pearson correlation coefficient (PCC)=0.701 [P<0.001] and 0.678 [P<0.001], respectively), as well as the peak systolic radial and circumferential strain rate (PCC=0.613 [P<0.001] and 0.611 [P<0.001], respectively), showed highly significant correlations with interstitial myocardial fibrosis (Figure 5, Supplemental Table 1). We additionally observed a significant correlation of these deformation parameters with factors of left ventricular hypertrophy as cardiomyocyte cross-sectional area and heart weight (Figure 5, Supplemental Table 1).
Global circumferential and radial strain and hallmarks of UC are significantly correlated. Speckle-tracking parameters correlate highly significant with hallmarks of UC, such as interstitial fibrosis, myocardial weight, and cross-sectional cardiomyocyte diameter. CM, cardiomyocyte; SC peak G, strain circumferential peak global; SR peak G, strain radial peak global.
Importantly, the early diastolic peak SrR, as a diastolic left ventricular deformation parameter, also showed a close correlation with interstitial myocardial fibrosis (PCC=0.622; P<0.001) and with cardiomyocyte cross-sectional area and heart weight (PCC=0.7 [P<0.001] and 0.56 [P<0.001], respectively). This finding indicates that the systolic myocardial deformation intensity and the rate by which this deformation occurred during systole and diastole was significantly reduced in relation to increased cardiac hypertrophy and interstitial fibrosis.
Clinical Study
Patients
Of 985 patients with ESRD who underwent echocardiography in 2006, 2007, and 2008; 206 patients were identified as eligible with two-dimensional strain analysis accessible echocardiography. Thirty-five of these patients were excluded because of inappropriate image quality for strain analysis (Figure 6A). Thus, the study population comprised 171 patients with ESRD (111 men [64%]) undergoing dialysis (93% hemodialysis and 7% peritoneal dialysis; mean dialysis vintage±SD, 39±55 months; mean age, 64±14 years). Thirty-one of these patients had previously undergone kidney transplantation and had returned to dialysis 4.6±5.3 years before echocardiography.
Less negative strain longitudinal peak global is associated with a higher cardiovascular mortality in patients with ESRD. (A) Flow chart of the clinical study. (B) Kaplan–Meier analysis of cumulative cardiovascular survival for patients with a peak global longitudinal strain stratified at the median (−11.83), indicating that a less negative peak global longitudinal strain (lower myocardial deformation) was associated with a higher cardiovascular mortality (primary end point). Statistical comparisons between survival curves were made by log-rank test. 2D, two-dimensional; SL peak G, strain longitudinal peak global.
Follow-Up
During the 2.1±0.9-year follow-up period, 42 patients (25%) reached the primary end point (cardiovascular mortality), on average after 311±302 days. The total mortality (secondary end point) during follow-up was 44% (i.e., 75 patients died, after a mean of 302±285 days). In addition to the primary end point, causes of death were septicemia in 19 patients (11%), cancer in 7 (4%), and unknown in 7 (4%). Seventeen patients (10%) received a kidney transplant during follow-up, on average after 442±274 days. All 171 patients were included in the statistical analyses.
Baseline characteristics of survivors, nonsurvivors, and patients who reached the primary end point are outlined in Table 3. Compared with healthy persons (n=50), patients with ESRD showed a significant reduction of various left ventricular strain values (Supplemental Table 2).
Baseline demographic and clinical characteristics stratified for survivors, nonsurvivors, and patients who died of cardiovascular disease during follow-up
In ESRD patients with diabetes mellitus (type 1 or 2), peak global longitudinal strain and early diastolic strain rate longitudinal were reduced by 10% and 17%, respectively, compared with ESRD patients without diabetes mellitus (−12.51±3.23 versus −11.26±3.3 [P=0.02] and 0.93±0.32 versus 0.77±0.32 [P=0.002], respectively) (Supplemental Table 3). Of note, less negative peak global longitudinal strain reflects reduced left ventricular deformation. We did not observe significantly altered strain parameters in hypertensive patients with ESRD compared with ESRD patients without hypertension (Supplemental Table 4).
Strain longitudinal peak global was correlated with age (r=0.19; P=0.01); however, there was not a significant correlation with body mass index (r=0.13; P=0.11) (Supplemental Table 5).
Intra- and Interobserver Agreement
For intraobserver agreement, the Lin coefficient was 0.98 (95% confidence interval [95% CI], 0.97 to 0.99) for peak global longitudinal strain and 0.95 (95% CI, 0.90 to 0.990) for peak systolic longitudinal strain. For interobserver agreement, the Lin coefficient was 0.96 (95% CI, 0.94 to 0.98) and 0.90 (95% CI, 0.86 to 0.97), respectively. The Lin coefficients for the other strain parameters showed similar values.
Strain Parameters Predict Cardiovascular and All-Cause Mortality
In a Kaplan–Meier survival analysis, hemodialysis patients were categorized as being above (more negative) or below (less negative) the median (−11.83) of longitudinal peak global strain. This analysis demonstrated that less negative longitudinal strain (peak global) was associated with an increased cardiovascular mortality (log-rank test: P<0.001) (Figure 6B).
Univariate Cox regression analysis revealed that apart from ejection fraction, several strain variables significantly differed between the patients who reached the primary or secondary endpoint and those who did not (Table 4).
Univariate Cox analysis
Multivariate analysis (Table 5) for cardiovascular death or all-cause mortality did not change this observation. Various strain parameters were strong independent risk factors for cardiovascular and all-cause mortality (Table 5). Adding parameters such as protein, albumin, phosphate, potassium, C-reactive protein, or hemoglobin as potential confounders did not substantially change the results.
Multivariate Cox analysis
Crude effects of the clinical variables age, hypertension, diabetes, and body mass index, plus the parameters ejection fraction and longitudinal peak global strain, are shown in Supplemental Table 6. We included all these parameters in a multivariate analysis where only age and longitudinal peak global strain remained independent predictors of cardiovascular mortality (Supplemental Table 7).
Discussion
Sudden cardiac death and chronic heart failure are the major causes of cardiovascular death in CKD and ESRD.5–7 UC is the underlying pathophysiologic process and is characterized by myocardial hypertrophy and interstitial myocardial fibrosis.9,25 Our study points toward a novel echocardiographic approach, STE, to detect UC and to predict cardiovascular and all-cause mortality in ESRD.
The pathogenesis of UC is incompletely understood. Increased levels of circulating factors, such as inflammatory cytokines, endogenous cardiotonic steroids, vasoactive hormones, and uremic toxins, might be involved.26–31 Activation of the renin-angiotensin system is thought to be an important contributor to myocardial remodeling because in patients with ESRD, the severity of left ventricular hypertrophy correlates with plasma aldosterone concentration and administration of angiotensin-converting enzyme inhibitors improved structural changes as interstitial fibrosis in experimental and clinical studies.32–34 Anemia, increased afterload due to aortic stiffness and hypertension, increased oxidative stress, and insulin resistance are other important contributors to the development of UC.25,35,36 Another factor is hyperphosphatemia; in rats with subtotal nephrectomy, a high-phosphorus diet increased the severity of interstitial fibrosis and microvascular disease.37 In patients free of known cardiovascular disease, dietary phosphorus was associated with a greater left ventricular mass.38 The possible link between high phosphorus levels and cardiac hypertrophy might be increasing serum levels of fibroblast growth factor-23, which can induce left ventricular hypertrophy.39
We observed uremia with hyperphosphatemia and cardiac hypertrophy with interstitial fibrosis in both rat models of UC. The hyperphosphatemia was achieved by dietary phosphorus supplementation in the model of subtotal nephrectomy (5/6 nephrectomy), whereas the adenine nephropathy model is known for its severe hyperparathyroidism.40,41 Serum phosphate levels correlated significantly with the grade of cardiac hypertrophy as heart weight and cardiomyocyte cross-sectional area and the severity of myocardial fibrosis.
TGF-β is a known key player in solid organ fibrosis42 and has been implicated in the development of UC.31,43 We observed a significant upregulation of TGF-β in both rat models of UC. Additionally, UC was accompanied by re-expression of the fetally expressed β-myosin heavy chain and increased cardiac expression of brain natriuretic peptide, both of which have been linked to cardiac fibrosis, hypertrophy, and UC.31,44,45
Our first major finding was that, in rat models, UC led to early changes of myocardial systolic and diastolic deformation detectable by STE, whereas classic echocardiographic parameters of left ventricular function, such as FS or CO, changed only at a later time point, when the rats exhibited greater myocardial fibrosis and hypertrophy. Martin et al. recently reported early changes in early diastolic circumferential strain rate assessed by STE in rats with mild renal insufficiency after uninephrectomy.43 This is in line with our finding of early changes in peak diastolic strain rate, but we observed predominant changes in peak global and systolic strain in our rat models of UC. This discrepancy might be explained by our more severe model of UC with more advanced CKD and higher degrees of interstitial fibrosis resulting in a left ventricular systolic impairment in addition to the reduced diastolic deformation.
Our second important finding was that speckle-tracking parameters correlated significantly with the grade of myocardial fibrosis and hypertrophy in both rat models of UC.
In line with our results, Pen et al. described a significant correlation of both radial and circumferential strain to interstitial myocardial fibrosis and myocyte cross-sectional area in mouse models of cardiac dysfunction (transverse aortic constriction and profound isoflurane anesthesia).24 However, these authors observed a higher correlation for circumferential strain with interstitial fibrosis than for radial strain.24 In a rat model of doxorubicin cardiomyopathy, Migrino et al. observed an early reduction of radial strain after 8 weeks; the fractional shortening changed later, after 12 weeks.46 The correlation of global and systolic strain and strain rate parameters with the severity of cardiac fibrosis might be explained by the fact that fibrosis represents a loss of functional cardiomyocytes that are replaced with scar tissue, resulting in impaired systolic left ventricular deformation.24 Moreover, interstitial myocardial fibrosis might lead to increased left ventricular stiffness, thereby explaining the correlation with diastolic strain parameters in particular early diastolic strain rate radial. Early diastolic strain rate has the advantage of reflecting the diastolic performance of all myocardial segments by a regional measurement.47 Recently, a significant correlation between left ventricular filling pressure and early diastolic strain rate was shown.48 Thus, this parameter may describe a sensitive and physiologic approach to integrating myocardial relaxation characteristics and hemodynamics. Furthermore, researchers reported that early diastolic strain rate is superior to tissue Doppler imaging–derived analysis and is independently associated with the outcome in acute myocardial infarction.47 In line with this, we found that early diastolic strain rate is an independent risk factor for all-cause mortality in patients with ESRD. Interestingly, this parameter is impaired in segments at risk despite near-normal systolic deformation.49
Late gadolinium enhancement imaging via cardiac magnetic resonance imaging is the gold standard method for the noninvasive evaluation of cardiac fibrosis.50 Fibrosis severity detected via late gadolinium enhancement correlates significantly with STE-assessed strain parameters in patients with hypertrophic cardiomyopathy and Fabry disease.19,20,23 However, in patients with higher-stage CKD and those with ESRD, gadolinium is contraindicated because of the risk of nephrogenic systemic fibrosis.51 The earlier detection of myocardial fibrosis in UC using STE compared with standard echocardiographic readouts can thus possibly be used for interventional trials in patients with CKD and ESRD.
Our third major finding was that speckle-tracking parameters, such as peak global longitudinal strain, peak systolic longitudinal strain rate, and peak global circumferential and radial strain, were independent risk factors for cardiovascular mortality. In 201 patients with heart failure (number of patients with CKD not stated), Cho et al. reported that global circumferential strain is an independent and better predictor of cardiac events (readmission for heart failure or cardiac death) than left ventricular ejection fraction.52
Several studies described a reduction of left ventricular strain parameters in patients with CKD or ESRD despite preserved left ventricular ejection fraction.53–55 Edwards et al. reported that patients with early-stage CKD (CKD stages 2–3) have impaired systolic strain and strain rate values (assessed by tissue Doppler imaging), whereas the conventional echocardiographic parameters of left ventricular function as ejection fraction remain unchanged compared with healthy persons without CKD.56 Rahkit et al. used tissue Doppler imaging in patients with CKD and did show that diastolic tissue velocity is an independent predictor of cardiovascular events.57 Recently, Liu et al. reported that global longitudinal strain is an independent predictor of all-cause mortality in patients with ESRD who have preserved left ventricular function.58 Although radial, circumferential, and longitudinal fibers are predominantly in different layers of the myocardium, we observed reduced myocardial deformation in all these orientations when we compared our ESRD patient cohort to healthy persons. As heart muscle is incompressible, strain parameters are interrelated in three dimensions; thus, by the time the chamber contracts and shortens in systole (circumferential and longitudinal strain), the wall thickens (radial strain).59
As a postprocessing computer algorithm of routine grayscale images, speckle-tracking analysis can easily been integrated into clinical standard procedures and may help to detect UC and define patients with ESRD at high risk for cardiovascular death. An early detection of UC might be the basis for interventional trials and could help to initiate secondary prevention strategies, such as improved chronic heart failure therapy and implantation of cardioverter-defibrillators, which might lead to a reduction of cardiovascular mortality in ESRD. In fact, the current rates of implantable cardioverter-defibrillators in patients with ESRD are very low despite these patients’ high risk for SCD60; thus, new strategies for risk assessment are urgently needed.
The study had some limitations. The radial and circumferential strain and strain rate values in this study reflect the myocardial contractility in the six midventricular segments because we determined strain values only from midventricular short-axis views in both rats and humans. Thus, we cannot exclude the possibility that the strain and strain rate parameters of the whole left ventricle, including the apex and basal levels, would be different. Because left ventricular contractility may differ between rats and humans,61 we cannot entirely extrapolate our observed animal echocardiography data to dialysis patients.
In conclusion, STE detects early interstitial fibrosis and changes in left ventricular contractility in rat models of UC and predicts cardiovascular and all-cause mortality in patients with ESRD.
Concise Methods
Animal Study
Local government authorities approved all animal experiments reported in this article, and all animal experiments adhered to the National Institutes of Health’s Guide for the Care and Use of Laboratory Animals. A total of 37 male Wistar rats (Charles River Laboratories, Borken, Germany) weighing 290–320 g were studied. We used two different models of UC:
1. We performed 5/6 nephrectomy in a one-step surgery: Rats were anesthetized (xylazine/ketamine intraperitoneal), and a 6-cm incision was made in the linea alba; the right kidney was exposed, unencapsulated, and removed after ligation of its pedicle. Thereafter, the left kidney artery was exposed. Two of the normally three segment arteries were ligated (7-0 suture; Monofil, Ethicon) after visualization of the renal perfusion of these arteries via short forceps clamping/unclamping to make sure that only one third of the kidney remains perfused. To aggravate the myocardial fibrosis, rats were fed a high-phosphorus (1.2%) diet (Altromin, Lage, Germany).
2. Adenine nephropathy was induced by feeding rats a diet containing 0.75% adenine and 2.5% protein (Altromin) over 4 weeks, as described previously.62 Because decrease in renal function in the adenine model takes about 2 weeks to reach a level of CKD similar to that in the 5/6 nephrectomy model directly after surgery, we compared rats in the adenine nephropathy model 6 and 10 weeks after start of the diet with the 5/6 nephrectomy rats 4 and 8 weeks after surgery, respectively (Figure 1). The 10-week adenine group was fed an additional week of the adenine diet after 6 weeks.
Echocardiography was performed at 6 and 10 weeks (i.e., 4 and 8 weeks after 5/6 nephrectomy surgery) (Figure 1). Serum samples were collected weekly from the tail vein. After echocardiography the rats were euthanized by exsanguination.
Echocardiography and Two-Dimensional Strain Analysis in Rats
After rats were anesthetized (80 mg/kg ketamine delivered intraperitoneally) and shaved, echocardiographic data were acquired using a Vivid I echocardiography machine (GE Healthcare, Horton, Norway) and an 11.5-MHz 10S-RS sector transducer. Image depth was 2–2.5 cm, and the frame rate was 190–225 frames per second. Images were stored and analyzed offline using EchoPAC-PC software (GE Healthcare). Left ventricular size and systolic and diastolic diameter of the interventricular and the posterior wall were measured using M-Mode in the parasternal short axis view at the level of the papillary muscles. FS was calculated as follows:
Aortic annulus diameter was measured during systole from two-dimensional images. Diameter of pulmonary annulus and velocity time integral of continuous wave Doppler images were obtained from the short-axis view. Stroke volume was estimated as described previously.63 CO was calculated by multiplication of stroke volume with heart rate.
Speckle tracking analysis was performed using the Q-analysis software within EchoPAC-PC software (GE Healthcare) as described previously.64 Left ventricular radial and circumferential peak systolic strain and strain rate were assessed in the midventricular parasternal short axis view using papillary muscles as landmarks.
Quantitative Real-Time PCR
Quantitative real-time PCR was carried out as previously described.65 TaqMan primers and probes were designed from sequences in the GenBank database using the Primer Express software (Applied Biosystems, Foster City, CA). Primers are listed in Supplemental Table 8.
Miscellaneous Measurements
Serum and urine biochemistries were performed by standard laboratory procedures using an automated analyzer.
Histomorphological Analysis
For histologic and immunohistochemical analyses, tissue sections were fixed in 3.7% formaldehyde for 24 hours, paraffin embedded, cut with a rotating microtome at 3-µm thickness (Leica Microsystems), and stained according to routine histology protocols. Immunohistochemical analysis was performed using primary antibodies specific for fibronectin (rabbit polyclonal, 1:50; EMD Millipore, Billerica, MA), collagen I and III (both goat polyclonal, 1:100; Southern Biotec, Birmingham, AL). Antigen demarcation was achieved by microwave pretreatment in citrate buffer (pH, 6.0). After blocking (Blocking Kit; Vector Laboratories, Burlingame, CA) and 3% H2O2 treatment, sections were incubated with the primary antibody. Biotinylated horse monoclonal anti-goat and horse monoclonal anti-rabbit antibody (Vector Laboratories) were used as secondary antibodies. Detection was realized using the Vectastin ABC Kit (Vector Laboratories). Slides were counterstained with hematoxylin.
Histologic Measurements
To visualize cardiomyocyte borders, fixed tissue was incubated in wheat germ agglutinin conjugated to Alexa Fluor 488 (Invitrogen) at 1 mg/ml in PBS containing 10 mM sodium azide. ImageJ software (National Institutes of Health, Bethesda, MD) was used to quantify cross-sectional cell surface area of 30 cells per field in each of four different myocardial regions (lateral, anterior, septal, and posterior) in a short-axis midventricular section. Quantification of myocardial fibrosis in Masson trichrome–stained sections, collagen I and III– and fibronectin-positive surface area, were determined in midventricular short-axis sections at 200× magnification (10 pictures per rat) using the number of stained pixels per total pixels in Adobe Photoshop CS5 (Adobe Systems, Inc., San Jose, CA).
Clinical Study
Patients
Patients with ESRD (undergoing hemodialysis or peritoneal dialysis treatment with a mean dialysis vintage of 38 months [determinable in 153 of 171 patients]) who had echocardiography in the RWTH Aachen University Hospital Department of Cardiology between January 1, 2006, and December 31, 2008, were screened for study eligibility. Patients with echocardiography that was no suitable for speckle-tracking analysis (i.e., echocardiography was not performed on the adequate machine or picture quality was too low for strain analysis) were excluded. Additionally, to compare the echocardiographic speckle-tracking values of patients with ESRD with strain values of healthy persons, we included 50 age-matched participants (31 men [62%]); mean age, 63±12 years; serum creatinine level <106 μmol/L) who had a normal left ventricular ejection fraction (>55%, biplane Simpson method) and, apart from a small percentage with hypertension (20%), were free of any known cardiovascular disease or other cardiovascular risk factors.
The study was approved by the ethical committee of the RWTH Aachen University Hospital (EK167/11) and carried out according to the Declaration of Helsinki. Follow-up data were obtained from patients with ESRD. Data on medical history, medications, and biochemical data were obtained by retrospective chart review.
Two-Dimensional Echocardiography and Speckle Tracking Analysis
The present study is based on echocardiographic two-dimensional loops that were all prospectively recorded in the echocardiography unit of the Cardiology Department in our hospital using a digital ultrasound scanner (Vivid 7; GE Healthcare) equipped with a 2.5-MHz transducer. Left ventricular apical two-, three-, and four-chamber views and left ventricular parasternal short-axis view at midventricle were acquired and processed. Left ventricular ejection fraction was measured offline on a personal computer with EchoPAC-PC software (version 110.1.3; GE Healthcare) by manual tracing of end-systolic and end-diastolic endocardial borders in apical two- and four-chamber views, applying the Simpson biplane method. The frame rate for the investigations was between 50 and 90 frames per second.
Speckle-tracking analysis was performed retrospectively and offline with the above-mentioned software package; we used the first cardiac cycle of each acquired loop as described elsewhere.16,66 Peak systolic and global longitudinal strain and peak systolic and diastolic longitudinal strain rate were obtained for each segment in an 18-segment left ventricular model (6 segments at each apical view). Peak systolic and global circumferential and radial strain and peak systolic and diastolic circumferential and radial strain rates were assessed in 6 segments of the parasternal short-axis view at the mid-ventricular level. The software package automatically grades the tracking quality of each segment on a scale ranging from 1.0 (optimal) to 3.0 (inadequate). Only segments with a tracking quality of ≤2.0 were included in the analysis. The investigator was blinded to the patients’ outcome data. All hemodialysis patients underwent echocardiography on a nondialysis day.
Outcomes
Cardiovascular death was defined as the primary end point. Death by any cause was the secondary end point. Follow-up was performed by review of the patients’ hospital chart or telephone interview with the patient or the patient’s primary healthcare providers and through the government death registry. Death occurring outside the hospital for which no other cause was specified was regarded as SCD and was included in the definition of cardiovascular death if the patient had known cardiovascular disease (i.e., prior myocardial infarction or chronic heart failure) and was not diagnosed with cancer.
Statistical Analyses
All analyses were done using SPSS software, version 18 (SPSS, Inc., Chicago, IL) using Windows 7. Continuous factors were summarized by means±SDs. Categorical data were presented by relative frequencies and percentages. A t test was used to compare differences between two groups for continuous variables, and the Fisher exact test was used to compare differences between two groups for categorical variables.
The Kaplan–Meier method was used to estimate cumulative survival. To identify the prognostic factors of all-cause and cardiovascular death, comparisons between survival curves were made by log-rank test. Cox regression was used (the “enter” method) to determine the effect of strain parameters on all-cause and cardiovascular mortality. The multivariable adjusted Cox regression analysis accounted for age, sex, diabetes, and presence of cardiovascular disease at the start of the study as possible confounders (i.e., known coronary artery disease, coronary stent, prior ST-segment elevation myocardial infarction or non–ST-segment elevation myocardial infarction, prior stroke, ischemic cardiomyopathy).
All test results are reported as hazard ratios, corresponding 95% CIs, and P values. In the animal study, for multiple group comparison, ANOVA with post hoc Scheffe correction was applied. Analyses were performed using the above-mentioned SPSS software or GraphPad Prism (5.0c, GraphPad Software, Inc., La Jolla, CA). Statistical significance was defined as P<0.05.
To define interobserver and intraobserver variability in analysis of myocardial deformation parameters, the same observer and a second independent observer repeated the analysis of 10 consecutive patients 1 week apart using the same two-dimensional echocardiographic loop and the same cardiac cycle. The Lin coefficient was calculated as an aggregate measure for agreement with a maximum range between −1 and +1. For continuous data, it represents an analog of the weighted κ coefficient determined for ordinal data.67
Disclosures
None.
Footnotes
M.B. and G.S. contributed equally to this work.
Published online ahead of print. Publication date available at www.jasn.org.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2013070734/-/DCSupplemental.
- Copyright © 2014 by the American Society of Nephrology

















