Identification of cardiomyopathy associated circulating miRNA biomarkers in patients with muscular dystrophy using a complementary cardiovascular magnetic resonance and plasma profiling approach
Journal of Cardiovascular Magnetic Resonance volume 18, Article number: 25 (2016)
Duchenne and Becker muscular dystrophy (DMD and BMD) are X-chromosomal recessive neuromuscular disorders that are caused by mutations in the dystrophin gene and characterized by cardiac involvement. Circulating microRNAs (miRNAs) have been proposed as diagnostic biomarkers for various cardiovascular diseases. However, circulating miRNAs reflecting the presence and/or disease severity of cardiac involvement in DMD/BMD patients have not been described so far.
Sixty-three male patients with known MD and 26 age-matched healthy male controls were prospectively enrolled. All MD patients and controls underwent comprehensive cardiovascular magnetic resonance (CMR) studies as well as venous blood sampling on the same day.
An impaired left ventricular (LV) systolic function (defined as LV-EF <55 %) was detected in 29 (46 %) and presence of late gadolinium enhancement (LGE) indicative of myocardial fibrosis in 48 (76 %) MD patients with an exclusively non-ischemic pattern. Whereas no significant differences were observed for the 27 selected circulating miRNAs in MD patients with abnormal CMR findings (comprising structural and/or functional impairments) compared to those with completely normal CMR studies, a significant up-regulation of three miRNAs was observed in LGE-positive MD patients compared to LGE-negative ones: miR-222 (1.8-fold, p = 0.035), miR-26a (2.1-fold, p = 0.03) and miR-378a-5p (2.4-fold, p = 0.026). A signature of these three miRNAs (miR-26a, miR-222 and miR-378a-5p) resulted in an area under the curve (AUC) value of 0.74 for the diagnosis of LGE-positive MD patients. In a multivariable model, three independent predictors for LGE presence were identified comprising not only clinical and laboratory markers (LV-EF: OR 0.47, 95 % CI 0.24-0.89, p = 0.021 and elevated hs-Trop: OR 2559, 95 % CI 2.97-22.04*105, p = 0.023) but also the circulating miR-222 (OR 938, 95 % CI 938.46, 3.56-24.73*104, p = 0.016).
Up-regulation of circulating miRNAs miR-222, miR-26a and miR-378a-5p indicates the presence of myocardial scars in MD patients. Plasma miR-222 appears to be a promising novel biomarker reflecting structural – but not functional – cardiac alterations in MD patients.
Duchenne and Becker muscular dystrophy (DMD and BMD) are X-chromosomal recessive neuromuscular disorders that are caused by mutations in the dystrophin gene that subsequently lead to either total absence or structural impairment of the dystrophin protein. Since dystrophin is a central protein in the cell membrane of skeletal as well as cardiac muscle cells, MD patients do not only suffer from skeletal muscle weakness and wasting but also from progressive cardiomyopathy [1, 2]. Since therapeutic treatment options for respiratory failure have tremendously improved within the last years, cardiac disease with a characteristic, non-ischemic pattern of left ventricular (LV) myocardial fibrosis leading to non-ischemic dilated cardiomyopathy, heart failure symptoms and ventricular arrhythmias has become a major cause of morbidity and mortality in MD patients [1, 3–5].
Within the last years, cardiovascular magnetic resonance (CMR) has evolved as an excellent tool for the early and sensitive diagnosis of cardiac involvement in MD patients, since this imaging modality does not only allow an accurate functional assessment of the human heart - but also enables non-invasive myocardial fibrosis detection based on techniques such as late gadolinium enhancement (LGE)-imaging or T1-mapping [6–8]. However, previous CMR-based studies in MD patients have also suggested that cardiac serum markers such as troponin or brain-natriuretic peptide (BNP) that are helpful to diagnose and monitor cardiac disease in many ischemic as well as non-ischemic cardiomyopathies, are of limited clinical value in MD patients [9, 10]. Hence, the identification of novel serum biomarkers for diagnosis and monitoring of cardiomyopathy in MD patients (e.g. as a gatekeeper for subsequent CMR studies) is a clinically important, however, still unsolved challenge.
Recent advances in molecular diagnostics have opened novel avenues in defining the molecular basis and underlying pathophysiology of inherited diseases. In particular, microRNAs (miRNAs) have been identified as a novel class of both biomarkers and targets for therapy . miRNAs are small (19–25 nucleotides), non-coding molecules and can modify gene expression by regulating mRNA stability or translation and thereby modify essential cellular functions . Recent studies have shown that miRNAs are not only located inside cells (where they exert their function) but can also be found in different fluids such as blood . Moreover, there is a growing body of evidence that circulating miRNAs can be used as diagnostic as well as prognostic biomarkers for different cardiovascular diseases [14, 15]. However, circulating miRNAs reflecting the presence and/or disease severity of cardiac involvement in DMD/BMD patients have not been described so far.
Hence, the aim of the current study was to identify specific circulating miRNAs in the plasma of DMD/BMD patients that would allow a non-invasive and accurate diagnosis of cardiac disease in these patients.
As part of an ongoing prospective study, 63 male patients with known MD were prospectively enrolled between 2008 and 2015. A diagnosis of DMD (N = 12) or BMD (N = 51) has previously been made in specialized neurology centers based on clinical data, skeletal muscle pathology with dystrophin analyses and/or genetic testing [10, 16–18]. The clinical degree of skeletal myopathy was clinically assessed as follows: 0 = no clinical signs of myopathy; 1 = able to walk, unable to run; 2 = unable to walk; 3 = unable to use hands.
In addition, 26 age-matched healthy male controls were enrolled between 2011 and 2014 and represented the control group. All MD patients and controls underwent comprehensive CMR studies as well as venous blood sampling on the same day. The study protocol complies with the Declaration of Helsinki and was approved by the local ethics committee (Ethik-Kommission Landesärztekammer Baden-Württemberg, Stuttgart, Germany). Informed consent was obtained from the patients prior to study inclusion.
Blood sampling in MD patients and controls
EDTA blood samples were collected on the same day of the CMR study. EDTA plasma was harvested by centrifugation of EDTA Blood Collection Tubes (Sarstedt, Germany) for 10 min at 4000 rpm. Aliquots of plasma supernatant were stored in cryotubes (Sarstedt, Germany) at −20 °C until use. Both in MD patients and controls, laboratory determinations for cardiac biomarkers - high sensitive troponin I (hs-Trop) and brain natriuretic-peptide (NT-proBNP) - were performed using standard methods and considered elevated when serum levels exceeded the upper laboratory reference limit. In addition, creatine kinase (CK) levels were determined.
miRNA extraction and quantification
RNA was extracted from 400 μl of plasma using mirVana miRNA isolation kit (Life technologies, USA) following the manufacturer’s protocol, eluted in 75 μl elution solution and stored at −80 °C.
miRNA selection for quantification was based on literature data related to either cardiovascular diseases and/or DMD/BMD (Additional file 1: Table S1) were reverse transcribed using TaqMan MicroRNA Reverse Transcription Kit (Life technologies, USA). Individual stem-loop reverse transcription primers included in the predeveloped TaqMan miRNA assay (Life technologies, USA) were pooled at a final dilution of 0.05x for each primer. The final RT reaction volume of 7.5 μl contained 0.15 μl 100 mM dNTP, 1.5 μl multiscribe reverse transcriptase (50 U/μl), 0.75 μl 10 × RT buffer, 0.095 μl RNase inhibitor (20 U/μl), 3 μl primer pool and 2 μl of total RNA. The reaction was performed following conditions of manufacturer.
To improve sensitivity of miRNA quantification, a pre-amplification reaction was performed. TaqMan miRNA assays included in the TaqMan miRNA assay (Life technologies, USA) were pooled at a final dilution of 0.2x for each assay. Pre-amplification reaction was done at 10 μl final volume containing 5 μl TaqMan PreAmp Master Mix (2X), 1.5 μl of assay pool, 2.5 μl of nuclease-free water and 1 μl of cDNA. The pre-amplification PCR was run according to the manufacturer’s protocol, the pre-amplification PCR product was diluted 1:5 with suspension buffer (Teknova AS, Norway) and stored at −20 °C until need.
The miRNA expression levels were quantified by real-time PCR using TaqMan® Universal Master Mix II (no UNG) and TaqMan miRNA assays (Life technologies, USA) on a real-time PCR BioMark system (Fluidigm Corporation, USA) following the manufacturer’s protocol. Relative levels of miRNA expression were calculated by normalization to expression levels of miR-16 and thereafter multiplied by 103 in order to increase readability in the respective tables. The following miRNAs had to be excluded from final analysis due to failing measurements: miR-1, miR-31, miR-34c, miR-95, miR-133a, miR-208a, miR-208b, miR-499a-3p, miR-499a-5p and miR-539.
CMR data acquisition
ECG-gated CMR studies were performed on a 1.5-T scanner (Aera, Siemens Medical Solutions, Erlangen, Germany) using commercially available cardiac software, electrocardiographic triggering, and cardiac-dedicated surface coils. Cine-imaging was performed using a steady-state-free-precession (SSFP) sequence in three long-axis slices (four-, three- and two-chamber) and a stack of short-axis slices completely covering the LV. LGE-imaging was performed using a T1-weighted inversion recovery gradient-echo sequence 10–15 min after intravenous contrast administration (0.15 mmol/kg Magnevist®) in the same imaging planes as the cine-images.
CMR data analysis
CMR analysis was performed off-line by two experienced readers blinded to gender and clinical characteristics. Ventricular volumes, ejection fraction and LV mass were derived by contouring endo- and epicardial borders on the short-axis cine images and indexed to body surface area. LGE presence and pattern were first visually assessed on the short-and long-axis images by using the 16-segment AHA model . Second, LGE extent was planimetered on the short-axis contrast images with the use of ImageJ software (National Institutes of Health, Bethesda, Md, USA) and an image intensity level ≥3 SD above the mean of remote myocardium was used to define LGE indicative of damaged myocardium as described previously (National Institutes of Health, Bethesda, Md, USA) and expressed as percentage of total LV mass .
An abnormal CMR was defined by at least one of the following findings: i) LV ejection fraction (LV-EF) less than 55 %, ii) RV ejection fraction (RV-EF) less than 45 %, iii) presence of LGE in at least one myocardial segment (AHA segmentation), and was considered as sign of cardiac involvement.
Continuous variables are expressed as mean ± SD. Skewed variables are expressed as median and interquartile range (IQR). Categorical variables are expressed as frequency with percentage. t-Student test was used for comparison of normally distributed variables, while Mann-Whitney U test was used for comparison of non-normally distributed variables. Non-parametric Kruskal–Wallis test with Bonferroni post-hoc correction was used in case of multiple comparisons of non-normally distributed variables. The Chi-square test with Yate’s correction was used to compare non-continuous variables expressed as proportions. Parametric Pearson or non-parametric Spearman correlations were used as corresponded for correlation analysis. In order to find independent predictors for abnormal CMR findings, i.e. LGE presence, a univariable regression analysis was first performed. Second, the parameters with significant p-values were introduced into the multivariable regression analysis. Extremely skewed miRNAs (skewness statistic < −2 or >2) were Log10 transformed before introduced in the regression analysis. Finally, receiver operating characteristic curves (ROC) were analyzed to assess specificity and sensitivity of single plasma miRNAs as well as their combination using multiple logistic regression analysis. Statistical analysis was performed using SPSS software for Windows (version 20, SPSS, Chicago Illinois, US). A p-value ≤ 0.05 was considered statistically significant.
The study group consisted of 63 male MD patients comprising 12 (19 %) patients with DMD and 51 (81 %) with BMD with a median age of 31 ± 15 yrs (Table 1). DMD patients were younger compared to BMD ones (18 ± 5 yrs vs. 33 ± 15 yrs, p < 0.0001). The control group comprised 26 healthy male volunteers aged 36 ± 13 yrs (p = 0.09). There was no patient with a history of coronary artery disease (CAD) and/or ischemic cardiomyopathy and/or valvular disease. The respective cardiovascular risk profile is illustrated in Table 1. Serum measurements revealed elevated hs-Trop values in 30 (48 %) MD patients (p < 0.0001 vs. controls) whereas NT-proBNP levels were increased in 3 (5 %) MD patients only (p = 0.55 vs. controls).
Major CMR findings
The detailed results of the CMR studies performed in MD patients and controls are given in Table 2. An impaired LV systolic function (defined as LV-EF <55 %) was detected in 29 (46 %) MD patients (p < 0.0001 vs. controls) and a reduced RV systolic function (defined as RV-EF <45 %) in 10 (16 %) MD patients (p = 0.031 vs. controls). Presence of LGE indicative of myocardial fibrosis was detected in 48 (76 %) MD patients (p < 0.0001 vs. controls) with a median LGE extent of 10 % (IQR 5 %–24 %) of LV mass and an exclusively non-ischemic pattern (Fig. 1a-b). There was a strong negative correlation between LGE extent and LV-EF (Spearman’s rho = -0.753, p < 0.0001). Altogether, any pathological CMR finding was observed in 49 (78 %) MD patients whereas all subjects in the control group demonstrated normal CMR findings (p < 0.0001).
No significant differences were seen between DMD and BMD neither in frequencies of LV and/or RV systolic impairment (42 % vs. 47 %, p = 1.00 for LV; 25 % vs. 14 %, p = 0.39 for RV) nor in LGE presence (92 % vs. 73 %, p = 0.26) or LGE extent. Consequently, there was also no significant difference in the prevalence of an abnormal CMR study between these groups (92 % vs. 75 %, p = 0.27).
miRNA findings in MD patients vs. controls
Plasma miRNA results are shown in Table 3. A significant up-regulation in MD patients compared to controls was found for seven out of 17 plasma miRNAs included in the analysis: miR-206 (91-fold increase, p < 0.0001), miR-20b (1.5-fold increase, p = 0.048), miR-222 (24-fold, p < 0.0001), miR-26a (6-fold, p < 0.0001), miR-342 (33-fold, p < 0.0001), miR-378a-3p (minimum 3000-fold; almost undetectable in controls, p < 0.0001), miR-378a-5p (48-fold, p < 0.0001). Additionally, a significant down-regulation was found for the three miRNAs: miR-221, miR-29a and miR-29c, all almost undetectable in the plasma of MD patients.
Among MD patients, none of the miRNAs differed significantly in patients with elevated vs. normal hs-Trop and NT-proBNP, respectively. Moreover, we did not detect any significant differences in miRNA expression in MD patients with mild myopathy (clinical degree 0-1 as defined in our Methods section) compared to those with advanced/severe myopathy (clinical degree of 2–3) (Additional file 2: Table S2).
When we assessed the relationship between plasma miRNAs and age, significant but moderate correlations were found only for two of them: miR-26a (Spearman’s rho = +0.312, p = 0.013) and miR-378a-3p (Spearman’s rho = −0.299, p = 0.017) in the MD group. In contrast, no significant relationship was found in controls regarding plasma miRNA levels and age.
Assessment of CMR and miRNA findings considering all CMR results
When MD patients with any abnormal CMR findings (n = 49; 78 %) were compared to those without (n = 14; 22 %), there was no significant difference regarding age or MD type (Tables 4 and 5). An elevated hs-Trop plasma level was more frequently observed in patients with abnormal CMR findings compared to those without. However, only 57 % (n = 28) of those patients with pathological CMR findings also demonstrated elevation of hs-Trop. Moreover, there were no significant differences in NT-proBNP levels between MD patients with abnormal CMR findings and those without.
As shown in Table 6, there were no significant differences in any of the measured plasma miRNAs in MD patients with any abnormal CMR finding (comprising structural and/or functional impairments) compared to those with completely normal CMR studies.
Assessment of CMR and miRNA findings according to “functional” CMR results
Fifty-one percent (n = 32) of the MD patients showed a reduced LV and/or RV systolic function in the CMR study as measured by the respective EF. In this subgroup of MD patients with impaired systolic function, both LGE prevalence and extent were significantly increased compared to those MD patients with normal biventricular systolic function (n = 31; 49 %) as depicted in Tables 7 and 8.
However, there were again no significant differences in any of the measured plasma miRNAs in MD patients with an impaired systolic function compared to those without (Table 9). The only significant correlation between a miRNA and LV-EF was a moderate negative one for miR-29c (Spearman’s rho = −0.323, p = 0.010).
Assessment of CMR and miRNA findings according to “structural” CMR results
In the next step, MD patients with presence of structural abnormalities - defined as presence of LGE (n = 48; 76 %) - were compared to those without any LGE (n = 15; 24 %). LGE-positive MD patients showed more frequently a reduced LV-EF (28 (58 %) vs. 1 (7 %); p = 0.001) and an elevated hs-Trop level (28 (58 %) vs. 2 (13 %); p = 0.003) compared to LGE-negative ones (Tables 10 and 11).
As illustrated in Table 12, a significant up-regulation of three miRNAs was observed in LGE-positive MD patients compared to LGE-negative ones: miR-222 (1.8-fold, p = 0.035), miR-26a (2.1-fold, p = 0.03) and miR-378a-5p (2.4-fold, p = 0.026). There was no miRNA that was significantly down-regulated in the LGE-positive group compared to the LGE-negative one.
When the miRNA results of LGE-positive and LGE-negative patients were compared to controls, there were significant differences when all three groups were considered (Fig. 2a-c). However, when post hoc analysis for multiple group comparisons was performed, differences between MD patients with and without LGE were only borderline significant for any of the miRNAs – as expected in case of multiple testing in a small-sized study group. While in this post hoc comparison circulating miR-222 and miR-378a-5p levels were significantly higher both in LGE-positive and LGE-negative patients compared to controls (p = 0.003 and p < 0.0001 for miR-222; p = 0.017 and p < 0.0001 for miR-378a-5p), a significant difference was noted only between LGE-positive patients and controls for miR-26a (p < 0.0001) – but not for LGE-negative ones vs. controls.
Assessment of possible predictors of LGE presence
In order to further characterize the association between myocardial scarring as detected by LGE-CMR and the three up-regulated plasma miRNAs, we first performed univariable logistic regression analyses for a series of potential predictors of LGE occurrence - including miR-222, miR-26a and miR-378a-5p (Table 13). In this analysis, a significant association with LGE presence was found for the following parameters: a) CMR parameters such as LV-EF, LV-ESVi and LV-mass, b) the serum marker elevated hs-Trop and c) the miRNAs miR-222 and miR-26a.
In the next step, we performed a multivariable regression analysis focusing on four variables of most interest among those being statistically significant variables from the aforementioned univariable analysis (LV-EF, elevated hs-Trop, miR-222 and miR-26a). In this multivariable model, three independent predictors for LGE presence were found: a) LV-EF (OR 0.47, 95 % CI 0.24-0.89, p = 0.021), b) an elevated hs-Trop (OR 2559, 95 % CI 2.97-22.04*105, p = 0.023) and c) serum miR-222 levels (OR 938, 95 % CI 3.56-24.73*104, p = 0.016).
Plasma miRNAs for the identification of MD patients with LGE presence
As shown in Fig. 3a-c, the individual receiver operating characteristics curves (ROC) for the three significantly up-regulated miRNAs in MD patients with vs. without LGE revealed areas under the curve (AUC) close to 0.70. Sensitivities, specificities and overall accuracies for these miRNAs were as follows: a) miR-26a: 65, 73 and 67 % (cut-off value of 241.54*10-3), b) miR-222: 62, 80 and 65 % (cut-off value of 1833.18*10-3) and c) miR-378a-5p: 69, 60 and 67 % (cut-off value of 18.84*10-3).
Furthermore, by combining these three miRNAs (miR-26a, miR-222 and miR-378a-5p, Fig. 3d) as potential diagnostic signature for cardiomyopathy in MD patients, an improved AUC value of 0.74 was reached with a sensitivity of 70 %, a specificity of 73 % and an overall accuracy of 70 %. In comparison, ROC analysis for the conventional serum marker hs-Trop regarding the identification of LGE-positive MD patients revealed an AUC value of 0.75 with a sensitivity of 63 %, a specificity of 87 % and an overall accuracy of 65 % (Fig. 3e). Moreover, there was no significant correlation between hs-Trop and each of the three miRNAs (miR-26a, miR-222 and miR-378a-5p).
To the best of our knowledge, this is the first study that assessed the diagnostic value of circulating miRNAs for the detection and prediction of functional as well as structural cardiac impairments (that were assessed by comprehensive CMR studies) in MD patients. The present results indicate that some circulating miRNAs are differently expressed in the plasma of MD patients compared to healthy volunteers. However, for the first time we also identified three circulating miRNAs (miR-26a, miR-222 and miR-378a-5p) that were significantly up-regulated in MD patients with proof of myocardial scarring (based on LGE-CMR) compared to those without – independent of the patients’ age and skeletal muscle status. Interestingly, the association between the plasma profile of these three miRNAs and the presence of cardiac disease was only observed in case of structural changes (proof of myocardial scarring) – but not in case of exclusive functional abnormalities.
Pathophysiology of cardiac involvement in MD
As outlined previously, the detailed molecular pathomechanism leading to cardiac disease in MD patients is still to be elucidated, although the underlying genetic dystrophin defect can be identified easily by appropriate mutation screening . Based on preclinical studies in dystrophin deficient animal models, early alterations in cardiomyocyte metabolism and signal transduction were suggested . In addition, excessive intracellular calcium signalling and reactive oxygen species (ROS) generation with breakdown of the mitochondrial membrane potential were described and may constitute the link between the initial sarcolemmal injury due to dystrophin deficiency and mitochondrial dysfunctions. Of note, De Arcangelis et al. showed in an animal model that miRNAs are involved in the regulation of the dystrophin-glycoprotein complex components in dystrophic muscles of the mdx mice. Especially, miR-222 expression was increased in dystrophic muscle and led to a decrease of the b1-syntrophin expression . Taken together, the fragility of the cell membrane that is caused by deficient sarcolemmal dystrophin may predispose not only skeletal muscle cells but also cardiomyocytes to metabolic dysfunctions, which in turn may be enhanced by mechanical stress.
In recent years, CMR imaging has gained wide acceptance for non-invasive evaluation of ischaemic as well as non-ischaemic cardiomyopathies since multi-parametric CMR does not only allow the assessment of functional parameters (such as systolic function) but also enables a non-invasive assessment of myocardial scarring based on LGE-imaging or T1-mapping . Within the last years, the distribution pattern of myocardial scarring in MD patients (as well as female MD carriers) was assessed based on such comprehensive CMR studies and it could be consistently shown that the posterolateral LV wall segments represent the first and most extensive sites of myocardial fibrosis in MD patients – as well as their female MD carrier relatives [17, 23–25] (Fig. 1). Moreover, we could recently show that a “transmural” pattern of myocardial fibrosis independently predicts the occurrence of adverse cardiac events in DMD/BMD patients .
Conventional serum markers of cardiac disease in MD patients
Today, measurement of serum BNP levels is routinely performed in patients presenting with heart failure symptoms and/or suspected cardiomyopathy (e.g. dilated cardiomyopathy) since normal BNP values allow to elegantly rule out congestive heart failure whereas elevated BNP values (in patients with normal renal function) point to an underlying cardiac disease and necessitate further cardiac examinations such as echocardiography and/or CMR . Interestingly, in the present study (comprising 63 MD patients), elevated Nt-proBNP levels were measured in three (5 %) MD patients only – although any abnormal cardiac findings were detected in 49 (78 %) and a reduced LV-EF was observed in 29 (46 %) of them by CMR. Hence, measurement of serum BNP is not very sensitive for diagnosing or ruling out cardiac disease in MD patients. As discussed previously, one explanation for this observation may be that serum BNP levels are supposed to reflect rapid changes of pressure gradients in the atria and ventricles and one may argue that the underlying pathophysiology in MD patients is mainly driven by continuous myocardial cell death (due to dystrophin absence or fragility) resulting secondarily in a continuous and adapted ventricular enlargement with progressive decrease of systolic function, however, without quick changes in atrial/ventricular filling pressures [10, 26, 27].
Furthermore, measurement of serum troponin is routinely performed nowadays not only in patients with acute chest pain syndromes but also in those with (suspected) non-ischemic cardiomyopathy for further detailed classification and risk stratification of such patients [28, 29]. In particular, an elevated level of high-sensitive troponin is believed to be an accurate serum marker of myocardial cell death in various acute settings. In the present study, an elevated serum level of hs-Trop was detected in 30 (48 %) MD patients only - although (as aforementioned) any abnormal cardiac findings were detected in 49 (78 %) and proof of myocardial scarring/damage in 48 (76 %) of the study MD patients by CMR. Hence, measurement of serum hs-Trop in MD patients does not allow to detect all of those MD patients who have at least some structural cardiac abnormalities - as can be depicted by LGE-CMR. Moreover, even hs-Trop is not 100 % specific for cardiac disease and may also originate from cell death of skeletal muscle myocytes – particularly in patients with neuromuscular disorders and advanced skeletal myopathy. Taken together, a normal hs-Trop level in a MD patient does not always rule out cardiac disease and in contrast, an elevated hs-Trop level does not necessarily indicate the presence of cardiac abnormalities. Therefore, novel (more sensitive as well as more specific) markers are wanted for early and accurate diagnosis of cardiac involvement in MD patients.
Circulating miRNAs as novel biomarkers of cardiac disease in MD patients
Recent advances in molecular diagnostics have shown that circulating miRNAs can be used as diagnostic as well as prognostic biomarkers for different cardiovascular diseases [14, 15]. Therefore, we tested the hypothesis that the pattern of plasma miRNA expression would also allow the diagnosis and prediction of cardiac disease in MD patients. In a first step, we searched for those miRNAs that were known to be associated with either cardiovascular diseases and/or DMD/BMD based on the available literature data [30–38]. Compared to our control group, we identified seven significantly up-regulated plasma miRNAs (miR-206, miR-20b, miR-222, miR-26a, miR-342, miR-378a-3p, miR-378a-5p) and three down-regulated miRNAs (miR-221, miR-29a and miR-29c) in MD patients (Table 3). Since we did not detect any significant differences in miRNA expression in MD patients with mild myopathy (clinical degree 0-1 as defined in our Methods section) compared to those with advanced/severe myopathy (clinical degree of 2-3), the aforementioned seven up-regulated and three down-regulated miRNAs do not seem to reflect the skeletal myopathy status of our MD patients.
In the next step, we looked at those MD patients who had any abnormal CMR findings compared to those with completely normal CMR results. Interestingly, there were no significant differences in any of the measured plasma miRNAs (Table 6). Then, we focused first on “functional” cardiac parameters and thereafter on “structural” ones. Again, we did not detect significant differences in any of the measured plasma miRNAs in MD patients when those with an impaired systolic function were compared to those without (Table 9). However, patients with proof of myocardial scarring (LGE-positive) compared to those without (LGE-negative) showed a significant up-regulation of three miRNAs (miR-222, miR-26a and miR-378a-5p, Table 12). Importantly, there was no significant difference in the degree of skeletal myopathy or serum CK in LGE-positive MD patients compared to LGE-negative ones. Hence, the elevation of these three miRNAs in LGE-positive MD patients was not (only) caused by skeletal myopathy, but rather reflected cardiac disease.
Additional univariable analyses revealed a significant and substantial association between the expression level of two miRNAs (miR-222 and miR-26a) and the presence of myocardial scarring. Furthermore, multivariable regression analysis revealed three independent predictors for the presence of LGE – one being a miRNA (miR-222) in addition to a clinical (LV-EF) and a laboratory marker (elevated hs-Trop). Taken together, an up-regulation of three circulating miRNAs miR-222, miR-26a and miR-378a-5p is a useful signal for the presence of cardiac (structural) disease in MD patients – with miR-222 being the strongest and most important biomarker.
Diagnostic value of plasma miRNAs compared to hs-Trop
In the present study, serum hs-Trop was a significant predictor of myocardial scarring – similar to miR-222. For further evaluation of the diagnostic value of the detected miRNAs in comparison to hs-Trop, we performed ROC analyses regarding the identification of LGE-positive MD patients. Area under the curve (AUC) values for the individual miRNAs (miR-222, miR-26a and miR-378a-5p) were close to 0.70. However, combining these three miRNAs resulted in an improved AUC value of 0.74 – compared to 0.75 in case of hs-Trop. However, the respective sensitivity, specificity and overall accuracy was 70, 73 and 70 % for the combination of these three miRNAs compared to 63, 87 and 65 % in case of hs-Trop. Hence, a combined approach based on the measurement of the three plasma miRNAs miR-222, miR-26a and miR-378a-5p seems to be at least as sensitive and accurate for the non-invasive diagnosis of cardiac disease in MD patients as the conventional hs-Trop measurement. Moreover, hs-Trop elevations may also have extra-cardiac causes (e.g. renal failure or myositis) – which so far were not described for the three identified circulating miRNAs. In future studies, the potential superior diagnostic as well as prognostic value of a miRNA-based approach for diagnosis of cardiomyopathy compared to established measurements (e.g. hs-Trop) needs to be evaluated in larger study groups.
Comparison of present miRNA results to previous literature data
In the last years, several studies have been published addressing both pre-clinical animal models and humans as well as different biological materials such as skeletal muscle, heart muscle and serum/plasma. Previously published studies that (amongst others) addressed the three miRNAs miR-222, miR-26a and miR-378a-5p are summarized in the Additional file 3: Table S3. For example, miRNA analyses in the skeletal muscle of mdx mice and DMD patients revealed an up-regulation of miRNAs associated with muscle regeneration (e.g. miR-206) and inflammation (e.g. miR-222) and a down-regulation of miRNAs associated with muscle degeneration (e.g. miR-29c) .
Recently, Jeanson-Leh et al. analysed circulating miRNA profiles in golden retriever muscular dystrophy dogs and also evaluated the association of miRNA to cardiac disease . Importantly, evaluation of cardiac disease was based only on functional analyses using echocardiography. Although a dysregulation of numerous miRNAs (such as miR-1, miR-95, miR-133, miR-208a/b, miR-206, miR-378 and miR-499) was detected in the serum of these dogs, there was no correlation between any miRNAs and cardiac functional parameters in this model. In agreement with this study, we also detected a dysregulation of some circulating miRNAs (e.g. miR-206 and miR-378) in the plasma of MD patients. Notably, we elucidated only a significant association between miRNA and cardiac involvement in MD patients when we looked specifically at structural – but not functional - abnormalities using myocardial scarring based on LGE-CMR. Obviously, miRNA profiling may improve the understanding of the pathophysiology of a certain disease if in depth phenotypic characterisation is considered.
Since the number of study patients was limited, replication of our data in an independent cohort is mandatory to demonstrate clinical utility. Moreover, in the present study we focused on selected miRNAs with evidence from previous studies. Thus, it cannot be excluded that additional miRNAs might even improve the predictive ability of the identified miRNA signature e.g. by comprehensive microarray profiling for miRNA expression. Finally, this was a pilot study and the results need verification in a larger cohort, with a multi-center approach and with consideration of clinical outcomes to define the respective role of miRNAs.
Up-regulation of circulating miRNAs miR-222, miR-26a and miR-378a-5p indicates the presence of myocardial scars in MD patients. Plasma miR-222 appears to be a promising novel biomarker reflecting structural – but not functional – cardiac alterations in MD patients.
Becker muscular dystrophy
cardiovascular magnetic resonance
Duchenne muscular dystrophy
left ventricular end-diastolic volume
left ventricular ejection fraction
left ventricular end-systolic volume
brain natriuretic peptide
right ventricular end-diastolic volume
Yilmaz A, Sechtem U. Cardiac involvement in muscular dystrophy: advances in diagnosis and therapy. Heart. 2012;98(5):420–9.
Verhaert D, Richards K, Rafael-Fortney JA, Raman SV. Cardiac involvement in patients with muscular dystrophies: magnetic resonance imaging phenotype and genotypic considerations. Circ Cardiovasc Imaging. 2011;4(1):67–76.
Finsterer J, Stollberger C. The heart in muscular dystrophies. Z Kardiol. 2005;94(11):791–2.
Hermans MC, Pinto YM, Merkies IS, de Die-Smulders CE, Crijns HJ, Faber CG. Hereditary muscular dystrophies and the heart. Neuromuscul Disord. 2010;20(8):479–92.
Mavrogeni S, Papavasiliou A, Skouteli E, Magoutas A, Dangas G. Cardiovascular magnetic resonance imaging evaluation of two families with Becker muscular dystrophy. Neuromuscul Disord. 2010;20(11):717–9.
Florian A, Ludwig A, Rosch S, Yildiz H, Sechtem U, Yilmaz A. Myocardial fibrosis imaging based on T1-mapping and extracellular volume fraction (ECV) measurement in muscular dystrophy patients: diagnostic value compared with conventional late gadolinium enhancement (LGE) imaging. Eur Heart J Cardiovasc Imaging. 2014;15(9):1004–12.
Ikeda S, Kong SW, Lu J, et al. Altered microRNA expression in human heart disease. Physiol Genomics. 2007;31(3):367–73.
Matkovich SJ, Van Booven DJ, Youker KA, et al. Reciprocal regulation of myocardial microRNAs and messenger RNA in human cardiomyopathy and reversal of the microRNA signature by biomechanical support. Circulation. 2009;119(9):1263–71.
McMurray JJ, Adamopoulos S, Anker SD, et al. ESC guidelines for the diagnosis and treatment of acute and chronic heart failure 2012: The Task Force for the Diagnosis and Treatment of Acute and Chronic Heart Failure 2012 of the European Society of Cardiology. Developed in collaboration with the Heart Failure Association (HFA) of the ESC. Eur J Heart Fail. 2012;14(8):803–69.
Yilmaz A, Gdynia HJ, Baccouche H, et al. Cardiac involvement in patients with Becker muscular dystrophy: new diagnostic and pathophysiological insights by a CMR approach. J Cardiovasc Magn Reson. 2008;10:50.
Quiat D, Olson EN. MicroRNAs in cardiovascular disease: from pathogenesis to prevention and treatment. J Clin Invest. 2013;123(1):11–8.
Zampetaki A, Mayr M. MicroRNAs in vascular and metabolic disease. Circ Res. 2012;110(3):508–22.
Chen X, Ba Y, Ma L, et al. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18(10):997–1006.
Gupta SK, Bang C, Thum T. Circulating microRNAs as biomarkers and potential paracrine mediators of cardiovascular disease. Circ Cardiovasc Genet. 2010;3(5):484–8.
Wang GK, Zhu JQ, Zhang JT, et al. Circulating microRNA: a novel potential biomarker for early diagnosis of acute myocardial infarction in humans. Eur Heart J. 2010;31(6):659–66.
Cerqueira MD, Weissman NJ, Dilsizian V, et al. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Int J Cardiovasc Imaging. 2002;18(1):539–42.
Florian A, Rosch S, Bietenbeck M, et al. Cardiac involvement in female Duchenne and Becker muscular dystrophy carriers in comparison to their first-degree male relatives: a comparative cardiovascular magnetic resonance study. Eur Heart J Cardiovasc Imaging. 2015;25.
Yilmaz A, Gdynia HJ, Mahrholdt H, Sechtem U. Cardiovascular magnetic resonance reveals similar damage to the heart of patients with Becker and limb-girdle muscular dystrophy but no cardiac symptoms. J Magn Reson Imaging. 2009;30(4):876–7.
Yan AT, Shayne AJ, Brown KA, et al. Characterization of the peri-infarct zone by contrast-enhanced cardiac magnetic resonance imaging is a powerful predictor of post-myocardial infarction mortality. Circulation. 2006;114(1):32–9.
Khairallah M, Khairallah R, Young ME, Dyck JR, Petrof BJ, Des RC. Metabolic and signaling alterations in dystrophin-deficient hearts precede overt cardiomyopathy. J Mol Cell Cardiol. 2007;43(2):119–29.
De AV, Serra F, Cogoni C, Vivarelli E, Monaco L, Naro F. beta1-syntrophin modulation by miR-222 in mdx mice. PLoS One. 2010;5(8). doi:10.1371/journal.pone.0012098.
Ugander M, Oki AJ, Hsu LY, et al. Extracellular volume imaging by magnetic resonance imaging provides insights into overt and sub-clinical myocardial pathology. Eur Heart J. 2012;33(10):1268–78.
Florian A, Ludwig A, Engelen M, et al. Left ventricular systolic function and the pattern of late-gadolinium-enhancement independently and additively predict adverse cardiac events in muscular dystrophy patients. J Cardiovasc Magn Reson. 2014;16:81.
Silva MC, Meira ZM, Gurgel GJ, et al. Myocardial delayed enhancement by magnetic resonance imaging in patients with muscular dystrophy. J Am Coll Cardiol. 2007;49(18):1874–9.
Florian A, Ludwig A, Ong P, et al. Cause of cardiac disease in a female carrier of Duchenne muscular dystrophy: myocarditis versus genetic cardiomyopathy without skeletal myopathy? Circulation. 2014;129(18):e482–4.
Gabriel RS, Kerr AJ, Sharma V, Zeng IS, Stewart RA. B-type natriuretic peptide and left ventricular dysfunction on exercise echocardiography in patients with chronic aortic regurgitation. Heart. 2008;94(7):897–902.
Sutton TM, Stewart RA, Gerber IL, et al. Plasma natriuretic peptide levels increase with symptoms and severity of mitral regurgitation. J Am Coll Cardiol. 2003;41(12):2280–7.
Baccouche H, Mahrholdt H, Meinhardt G, et al. Diagnostic synergy of non-invasive cardiovascular magnetic resonance and invasive endomyocardial biopsy in troponin-positive patients without coronary artery disease. Eur Heart J. 2009;30(23):2869–79.
Hamm CW, Bassand JP, Agewall S, et al. ESC guidelines for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation. The Task Force for the management of acute coronary syndromes (ACS) in patients presenting without persistent ST-segment elevation of the European Society of Cardiology (ESC). G Ital Cardiol (Rome). 2012;13(3):171–228.
Eisenberg I, Eran A, Nishino I, et al. Distinctive patterns of microRNA expression in primary muscular disorders. Proc Natl Acad Sci U S A. 2007;104(43):17016–21.
Greco S, De SM, Colussi C, et al. Common micro-RNA signature in skeletal muscle damage and regeneration induced by Duchenne muscular dystrophy and acute ischemia. FASEB J. 2009;23(10):3335–46.
Hu J, Kong M, Ye Y, Hong S, Cheng L, Jiang L. Serum miR-206 and other muscle-specific microRNAs as non-invasive biomarkers for Duchenne muscular dystrophy. J Neurochem. 2014;129(5):877–83.
Jeanson-Leh L, Lameth J, Krimi S, et al. Serum profiling identifies novel muscle miRNA and cardiomyopathy-related miRNA biomarkers in Golden Retriever muscular dystrophy dogs and Duchenne muscular dystrophy patients. Am J Pathol. 2014;184(11):2885–98.
Orenes-Pinero E, Montoro-Garcia S, Patel JV, Valdes M, Marin F, Lip GY. Role of microRNAs in cardiac remodelling: new insights and future perspectives. Int J Cardiol. 2013;167(5):1651–9.
Roberts TC, Blomberg KE, McClorey G, et al. Expression analysis in multiple muscle groups and serum reveals complexity in the microRNA transcriptome of the mdx mouse with implications for therapy. Mol Ther Nucleic Acids. 2012;1, e39.
Thum T. Noncoding RNAs and myocardial fibrosis. Nat Rev Cardiol. 2014;11(11):655–63.
Vignier N, Amor F, Fogel P, et al. Distinctive serum miRNA profile in mouse models of striated muscular pathologies. PLoS ONE. 2013;8(2), e55281.
Wei C, Kim IK, Kumar S, et al. NF-kappaB mediated miR-26a regulation in cardiac fibrosis. J Cell Physiol. 2013;228(7):1433–42.
We would like to acknowledge the expert technical assistance of Ursula Waldherr and Andrea Jarmuth.
This work was supported by the Robert-Bosch-Foundation (Stuttgart, Germany) and the Deutsche Forschungsgemeinschaft (Bonn, Germany) grant DFG SCH 858/1-2.
The authors declare that they have no competing interests.
SB carried out the molecular genetic studies and performed the statistical analysis; AF performed the CMR studies, performed the statistical analysis and helped to draft the manuscript; AP performed the CMR studies; SR participated in performing the CMR studies and collected the blood specimens; JW helped to draft the manuscript; US participated in the design of the study and helped to draft the manuscript; MS participated in the design of the study and helped to draft the manuscript; ES participated in the design and coordination of the study and helped to draft the manuscript; AY conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Becker, S., Florian, A., Patrascu, A. et al. Identification of cardiomyopathy associated circulating miRNA biomarkers in patients with muscular dystrophy using a complementary cardiovascular magnetic resonance and plasma profiling approach. J Cardiovasc Magn Reson 18, 25 (2016). https://doi.org/10.1186/s12968-016-0244-3