1. Introduction
Cardiovascular disease (CVD), which refers to a wide range of diseases that affect the heart and blood vessels, including ischemic heart disease, stroke, and heart failure, is one of the leading causes of morbidity and mortality worldwide, accounting for more than 17 million deaths each year [
1,
2]. This condition appears to be highly prevalent in both developed and developing countries, and the rising prevalence of risk factors such as obesity and diabetes, along with population aging, is expected to increase the global burden of CVDs in the coming decades, bringing a massive burden on public health worldwide.
The pathogenic mechanisms of CVDs are complex and multifactorial, involving a combination of genetic, environmental, and lifestyle factors, and although recent research has been successful in uncovering novel risk factors and the underlying genetic and molecular mechanisms, seeking novel molecular players involved in the development of CVD is crucial to improve early diagnosis in high-risk patients and develop new potential therapeutic targets.
MicroRNAs (miRNAs) are endogenous, conserved, single-stranded non-coding RNAs of 21-25 nucleotides in length that regulate gene expression post-transcriptionally by base pairing to complementary sequences usually located to the 3′-untranslated region (UTR) of target messenger RNAs (mRNAs), leading to their degradation or translational repression [
3]. They are important regulatory elements able to fine-tune almost all biological pathways and are essential for maintaining cellular homoeostasis. Consequently, their dysregulation plays a relevant role in the onset and progression of diverse diseases [
4,
5].
In the cardiovascular system, miRNAs regulate a wide range of biological processes relevant to cardiovascular health, such as inflammation, oxidative stress, fibrosis, angiogenesis, cardiac cell contractility and growth [
6,
7], thus acting as promoters or suppressors of disease processes. Altered expression of several miRNAs has indeed been associated with several CVDs, holding promise as biomarkers for early disease detection [
8,
9,
10,
11,
12,
13].
Recently, a small group of miRNAs selectively enriched or expressed in cardiac and/or skeletal muscle (also called myomiRs), which includes miRs-1, -133a, -133b, -206, -208a, -208b, -499a, and -499b, has been proved to play important activities in the functional maintenance of the cardiovascular system as regulators of several processes such as proliferation, differentiation, regeneration, and metabolism [
14,
15,
16,
17,
18]. Evidence for altered myomiRs expression patterns in both human and animal models of cardiovascular related diseases has been reported and summarized in some recent reviews [
19,
20] and references therein]. To name a few, decreased expression of both miR-1 and miR-133 was found in mouse and human models of cardiac hypertrophy [
21] and reported in patients with myocardial infarction [
22]. In the same study by Boštjančič and colleagues [
22] miR-208 was found upregulated, in agreement with high expression levels in acute myocardial infarction patients found by Han and colleagues [
23]. Likewise, it has been reported that circulating miR-499 was substantially elevated in patients with acute myocardial infarction as compared with healthy control groups [
24]. In one report, Kumar and colleagues [
25] showed that miR-133b was downregulated in plasma samples of coronary artery disease patients.
In the present study the plasma expression patterns of miRNA-1, miRNA-133a, miRNA-133b, miRNA-206, miRNA-208b, and miRNA-499 were investigated in patients with different types of CVD and non-CVD controls, with the aim to further highlight the importance of these miRNAs as contributors to CVD pathogenesis and to explore whether they could be used as biomarkers for discrimination among different forms of CVD.
2. Results
As stated in materials and methods, subjects were classified as having cardiovascular disease (CVD) whether they had one or more of the following physician-confirmed diagnoses: ischemic cardiomyopathy, atrial fibrillation, heart failure and stroke. The demographic and clinical characteristics of the study cohort stratified according to absence/presence of CVD (CVD−/CVD+) are shown in
Table 1. Compared to CVD− controls, CVD+ patients were older (p = 0.001), while the ratio of males to females was not significantly different between the two groups (p > 0.05).
Risk factors for CVD such as diastolic pressure, serum total cholesterol, low-density lipoprotein (LDL) cholesterol levels were significantly lower in CVD patients than in control subjects (P<0.01); this is probably due to the use by patients of pharmacological treatments for blood pressure and lipid control. Conversely, levels of uric acid, an emergent biomarker of CVD, as well as marker of kidney dysfunction (urea and creatinine), resulted to be higher in CVD patients compared to controls. Finally, as expected, CVD patients were characterized by a higher mortality rate during the study follow-up (49 % vs 20%).
qRT–PCR was conducted to quantify the levels of miR-1, miR-133a, miR-133b, miR-206, miR-208a, miR-208b, and miR-499 in the plasma of healthy controls and CVD patients. For data analysis, we eliminated miR-1, miR-208b, and miR-499 since the expression levels were below the threshold of detection be quantified reliably and analyzed for differential expression.
First, we explored the association between these miRNAs and the variables in
Table 1 in the overall study population. A slightly positive correlation between and age (r = 0.14, p = 0.042) and a very weak negative correlation with the level of albumin and HDL cholesterol (r = -0.15, p = 0.037 and r = -0.17, p = 0.02, respectively) were observed for miR-133a. No other significant correlations between miRNA expression and laboratory data were observed.
We then examined whether the levels of circulating miR-133a, miR-133b, and miR-206 were associated with the diagnosis of CVD. Logistic regression analysis adjusted for age, sex and for statin exposure (
Table 2) showed that unitary increase in the expression of miR-133b and miR-206 was associated with a decreased risk of CVD (OR= 0.66; 95% CI= 0.44-0.97; p=0.038 and OR= 0.57; 95% CI= 0.36-0.85; p=0.009 for miR-133b and miR-206 respectively).
Since all subjects were followed-up for a median duration of approximately three years for their overall survival [28 healthy controls (27%) and 40 CVD patients (49%) died during this period], we also evaluate if the detected effects of the analyzed miRNAs on CVD risk might affect survival of the subjects. MiRNAs levels were divided by quartiles and the highest quartile was compared to others. The Kaplan-Meier survival curves are plotted and shown in
Figure S1. For all analyzed miRNAs, p-values for both log-rank test and Cox proportional hazards regression analysis adjusted for age and sex indicated no statistically significant differences in all-cause death between “high-expression” and “low/middle-expression” groups.
Next, we performed a sub-analysis to evaluate whether a characteristic miRNA expression profile is linked to a specific cardiovascular condition. The results, reported in
Table 3, indicated that, compared with the control group, unitary increase in plasma levels of miR-133a were significantly associated with decreased odds of heart failure by 53% (OR = 0.47, C.I. 0.24-0.95; p=0.031); a similar trend, although not significant at nominal level, was observed for miR-206 (OR= 0.62, C.I.0.37-1.06; p=0.074). Furthermore, increase in serum levels of miR-133b was significantly associated with decreased risk of stroke of about 37% (OR = 0.63, C.I. 0.39-0.99; p=0.048).
2.1. In Silico Functional Characterization
We performed miRNA-target enrichment analysis by using MIENTURNET tool setting the three miRNAs as the inputs and using experimentally validated miRNA-target interactions (miRTarBase). We retrieved one hundred thirty-nine genes significantly associated with these miRNA markers, see
Supplementary Table S1. As this table reports, a total of fifty-seven mRNAs were predicted to be targeted by miR-133a, forty-seven by miR-133b, and thirty-five by miR-206. Fifty-six mRNAs were in common between miR-133a and miR-133b, while four targets were shared between miR-133a and miR-206, and two between miR-133b and miR-206. Three mRNAs were targeted by all three miRNAs. The top-ranked targets of these miRNAs and the number of interactions is indicated in
Figure 1a, while the significant interaction network is shown in
Figure 1b. Functional enrichment analysis using the database REACTOME by MIENTURNET revealed the main pathways involved with the three miRNAs (
Figure 1c). Among these pathways, target genes of miR-133a and miR133b were involved in potassium channels, EGFR signaling, and Jak-stat signaling pathways. The latter pathway was also targeted by miR-206; other targets of this miRNA are involved in the dissolution of fibrin clot and mainly in signaling by the receptor tyrosine kinase MET.
3. Discussion
In recent years, interest has grown in determining the role of miRNAs in cardiovascular biology and in identifying circulating miRNA signatures related to the onset and progression of CVD, a leading cause of worldwide morbidity and mortality.
In the present study, we focused our attention on a subset of miRNAs, known as myomiRs, that are highly enriched in both skeletal and cardiac muscles and that have recently garnered attention for their involvement in most of the CVD-related processes, such as cardiac hypertrophy, cardiac fibrosis, apoptosis, inflammation (see references in the Introduction section). Six myomiRs were selected, but only three of them were detectable and analyzed in the plasma samples, namely miR-133a, miR-133b, and miR-206. The study showed that in the whole sample of patients with CVD (which included ischemic cardiomyopathy, atrial fibrillation, heart failure or stroke) the expression levels of both miR-133b and miR-206 were downregulated compared to those in healthy controls, suggesting that low levels of these miRNAs can have a detrimental effect on cardiovascular risk. When analyzing the levels of these miRNAs in subjects categorized according to specific CVDs, we observed a significant downregulation of miR-133b in patients with stroke and a significant downregulation of miR-133a in patients with heart failure, with a trend in the same direction for miR-206.
Taken all together, these findings suggest that specific molecular targets and cellular mechanisms involved in the action of these miRNAs likely act in response to different cardiac damages to obtain distinct expression patterns in distinct types of CVDs. This view is supported by findings of several studies from different cardiovascular conditions, which highlight the complex and context-dependent roles of these miRNAs in CVD development. In line with our findings, studies have found that patients with stroke have lower levels of miR-133b in their blood compared to healthy individuals, and that increasing miR-133b levels could promote recovery after a stroke [
26,
27]. Also, it has been reported that downregulation of miR-133b might contribute to the development of ventricular fibrillation in patients with myocardial infarction (MI) [
28], as well as to the risk of coronary artery disease (CAD) [
25]. By contrast, its upregulation has been reported in patients with ST-segment elevation MI and in a mouse model of MI [
29]. Also, association between increased concentration of miR-133b and increased odds of lipid-rich coronary plaque has been reported in patients with stable CAD [
30]. Furthermore, consistent with our observations, de Gonzalo-Calvo et al [
31] found that patients with heart failure have lower levels of miR-133a compared to healthy individuals. Consistently, it has been reported that synthetic miR-133a has the potential to prevent heart failure in a model of pressure overload [
32], and that a miR-133a mimic and miR-133a overexpression significantly decreased heart fibrosis in rats with chronic heart failure [
33]. On the other hand, the expression of this miRNA, which plays a role in the pathogenesis of most cardiovascular disease, is either upregulated or downregulated due to the pathological changes in particular cell types or specific CVD conditions. For instance, it is downregulated in patients with MI but, by contrast, its expression level is elevated significantly in patients with acute myocardial infarction (AMI) (for a review see [
34]). Regarding miR-206, while it was found to be downregulated when the whole group of patients with CVD was considered, only a slight decrease in its expression levels in the group of patients with heart failure was observed. Although not statistically significant, this fits well with the broad role of this miRNA in regulating heart muscle growth (hypertrophy), which is known to contribute to heart failure [
35,
36]. Several studies in model organisms evaluated miR-206 expression levels in different cardiac lesions. In rat heart with acute MI, it has been found that miR-206 expression is significantly decreased in infarcted myocardial areas and hypoxia-induced cardiomyocytes compared to non-infarcted areas; further investigation revealed that loss-of-function of miR-206 increased cardiomyocyte apoptosis [
37]. On the other hand, Shan et al [
38] found miR-206 expression significantly increased in a rat model of MI. Also, studies have shown that miR-206 expression is upregulated in cardiac hypertrophy [
35], and that its overexpression in adult cardiomyocytes causes arrhythmia [
39].
From all the above discussion, it appears that miRs-133a, -133b, and -206 are likely to exhibit differential expression patterns depending on the type of CVD and the clinical and physiological context. This raises the possibility that they can serve as diagnostic and/or prognostic markers not only for overall CVD but also for specific CVD types, pathophysiological processes and/or clinical outcomes. This possibility only comes through a better understanding of the regulation of these miRNAs, as well as the elucidation of their relevant targets, downstream pathways and biological functions in the context of the cardiovascular system, that presently have not been fully clarified.
The miR-133a forms a cluster with miR-1-2 on chromosome 18 (18q11.2), which is transcribed as a bicistronic transcript regulated by many myocytes’ differentiation factors, including myocardin, a nuclear protein that promotes activation of cardiac gene expression by associating with the transcription factor SRF (serum response factor) [
40]. MiR-133b and miR-206 are clustered on chromosome 6 (6p12.2) and transcribed as either a single bicistronic lncRNA, or under certain conditions, from separate promoters to yield separate primary transcripts [
41]. The primary transcript encoding miR-133b/miR-206 is induced by Wnt/β-catenin signals [
42]. This signaling pathway regulates the stability and nuclear localization of β-catenin and plays an important role cardiac function and cardiac tissue homeostasis [
43,
44]. Its activation is related to pathological stages after MI, including inflammation, angiogenesis, and fibrosis [
45,
46].
The target enrichment analysis showed that mRNAs of eleven genes were the most significant interacting targets, some of which deserve to be mentioned for their possible association with the cardiac function and pathophysiological state. Among the mRNAs targeted by both miR-133a and miR-133b were PRDM16, a zing-finger transcription factor whose cardiac-specific deletion results in cardiac hypertrophy, fibrosis, mitochondrial dysfunction, and heart failure [
47], the hyperpolarization-activated cyclic nucleotide–gated 4 (HCN4) channel, a transmembrane protein expressed in the cardiovascular system reported to play an important role in heart rate regulation [
48], and the glutathione transferase P1 (GSTP1), an enzyme involved in detoxification processes, that has gained attention as a biomarker for myocardial stress, particularly in the context of heart failure [
49,
50]. Among the mRNAs targeted by miR-133a and miR-206 were Annexin A2 (AnxA2), a Ca2+-dependent phospholipid-binding protein important for heart health that has been linked to mechanisms involved in CVD, such as angiogenesis, thrombosis and lipid metabolism [
51,
52], and the nuclear receptors NR4A2 that has been reported as associated to impairment of cardiac function and heart failure [
53,
54]. The results of the functional enrichment analysis showed that validated targets of both miR-133a and miR133b were enriched in biological processes related to regulation of cardiovascular function such as the potassium channels, membrane-spanning proteins that regulate cardiac repolarization and the shape and duration of the cardiac action potential and that are actively involved in the mechanism of several heart diseases [
55,
56,
57] and processes mediated by the EGFR signaling, an important signaling hub associated with vascular physiology whose relevance in cardiovascular health and disease has been shown [
58,
59,
60]. Target genes of all the three miRNAs were enriched in the Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway. This pathway, regulated by cytokines, triggers multiple signals involved in development, homeostasis and inflammation, participates plays a central role in transducing stress and growth signals in various myocardial injuries such as myocardial apoptosis, hypertrophy, fibrosis and ischemia [
46,
61,
62]. For miR-206 there was also a significant enrichment for gene sets related to the signaling pathways activated by the receptor tyrosine kinase MET, the HGF (hepatocyte growth factor) receptor. MET-dependent signals control heart homeostasis and prevents oxidative stress in normal cardiomyocytes and plays also important cardioprotective roles after heart injuries, such as ischemic injury, by promoting pro-survival effect such as protection from cardiomyocytes apoptosis and autophagy [
63,
64]. Furthermore, one of the enriched pathways was associated to dissolution of fibrin clot, the final step of the coagulation cascade. Impaired plasma fibrin clot properties characterize patients with acute and chronic manifestations of CVD [
65].
This study has many strengths: first, it provides valuable insights into the role of myomiRs in the development of cardiovascular diseases and highlights their potential as biomarkers, contributing to a growing body of evidence that could lead to improved diagnostic and therapeutic strategies for CVD. Second, in silico functional enrichment analysis add robustness to study results, highlighting targets and pathways that are likely involved in the pathophysiological mechanisms in CVD risk; third, this study directly addressed the clinical relevance of circulating miRNA in a real-world population of older nursing home residents, making study findings potentially highly relevant to clinical practice.
The study has also some limitations: the small sample size of the whole cohort and of subgroups of patients with different types of CVDs limit the generalizability of study findings; furthermore, the study cohort consisted of older nursing home residents enrolled in the Calabria region of Southern Italy; while this is relevant for studying cardiovascular disease in the aging population, the results may not be applicable to other age groups and geographic populations, as genetic and environmental factors may influence myomiRs expression.
4. Materials and Methods
4.1. Study Population
The cohort of this retrospective study included 209 individuals, which consisted of 90 patients (31% were males) clinically diagnosed with CVD including ischemic cardiomyopathy (52%), atrial fibrillation (29%), heart failure (18%) and stroke (31%) and a control group of 119 healthy controls (28% were males). The mean age of the study population was 85 ± 6.3 years for patients with CVD and 82 ± 7.8 years for the healthy controls.
All participants were recruited from nursing homes located in the provinces of Crotone and Cosenza, Calabria region (southern Italy), as part of a larger study aimed at assessing and monitoring the quality of aging in the entire region. Subjects were eligible to participate in the study if they were of Calabrian ancestry.
4.2. Clinical and Laboratory Assessments
Each participant underwent a complete clinical and geriatric assessment designed to evaluate cognitive status, functional abilities, physical health, and social aspects. A structured questionnaire administered by a trained interviewer was used to collect data. Demographic (age, sex) and anthropometric (height, weight, BMI) variables were available for all participants; clinical variables included systolic and diastolic blood pressure (SBP and DBP, respectively), history of hypertension and diabetes, and presence of CVD; the diagnosis of CVD was based on patients’ detailed medical histories, clinical symptoms, imaging studies and physical or laboratory findings as determined by a board-certified cardiologist in accordance with the international clinical guidelines. Additionally, a peripheral blood sample was drawn from each participant for laboratory examinations, including traditional biomarkers and myomiR, assessed via standardized protocols. Traditional biomarkers included fasting lipid levels (total cholesterol, triglycerides, and high- and low-density lipoprotein cholesterol (HDL-c, and LDL-c, respectively), fasting plasma glucose, HbA1c, urea, creatinine, electrolytes, uric acid, total proteins, albumin, C-reactive protein (CRP), iron, ferritin, and total bilirubin.
4.3. Plasma RNA Extraction
The peripheral blood specimens from all participants were collected in an EDTA-containing tube and centrifuged at 12,000 rpm for 15 min. Afterward, the plasma was stored in a freezer at −80°C until analysis. The MiRNeasy Serum/Plasma Kit (Qiagen, Germany) was used to extract the total RNA from 200 µL of the plasma. In brief, QIAzol lysis buffer (1000µl) were added to the sample to facilitate lysis. Upon addition of 200 µL of chloroform and 3.5 μL of Arabidopsis thaliana miR-159a (assay ID 000338) as a spike-in control, the lysate was separated into an organic and aqueous phase by centrifugation, with RNA partitioning to the upper aqueous phase which was separated and mixed with 1.5 volumes of 100% ethanol and purified on the RNeasy MinElute spin column provided in the kit. RNA was eluted in 14 µL of RNase-free water. The RNA was quantified on a Qubit 2.0 Fluorometer (Life Technologies, Milan, Italy). The yield was around 30–50 ng/mL each sample.
4.4. Quantitative Real-Time Polymerase Chain Reaction
5 μL of RNA was reverse transcribed to cDNA using TaqMan® microRNA Reverse Transcription Kit (Life Technologies) and stem–loop specific RT primers for each of the miRNA examined: hsa-miR-1 (ID 002222); hsa-miR-133a (ID 002246); hsa-miR133b (ID 002247); hsa-miR-206 (ID 000510); hsa-miR-208b (ID 002290); hsa-miR-499 (ID 001045). RNU6B (U6) (ID 001973) was used as an endogenous control to normalize the miRNA qPCR data.
After incubation, quantitative real-time PCR reactions were performed using the TaqMan MicroRNA Assay kit (Applied Biosystems). Each PCR reaction contained 1.33 μl of the RT reaction product, 10 μl of TaqMan 2x Universal PCR Master Mix without uracil-N-glycoslyase (UNG), and 1 μl of 20x TaqMan MicroRNA Assay reagent in a total volume of 20 μl. Real-time reaction was carried out at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 60 s., using the QuantStudio3™ Real-Time PCR System (Applied Biosystems, Milan, Italy) with automatic baseline setting. All reactions, including the no-template controls, were run in triplicate. The relative expression levels of each miRNA in comparison with the normalizer were calculated using the comparative threshold (Ct) method 2−ΔCt [
66]. Expression levels (2
−ΔCt) were log-transformed to better fit a normal distribution.
4.5. Pathway Enrichment Analysis of the Predicted Targets
The Mienturnet (MicroRNA ENrichment TURned NETwork) web tool (
http://userver.bio.uniroma1.it/apps/mienturnet/) [
67] was used to access the miRTarBase 9.0 database to search for experimentally validated target genes of the myomiRs analyzed in our study, and to generate a graphical representation of the network of miRNA-target interactions identified by the enrichment analysis. Networks were constructed considering miRNA-mRNA interactions with strong functional level of evidence (such as Western blot and luciferase assay), a threshold of 1 for the minimum number of miR–target interactions and FDR < 0.05. The Reactome database was used for pathway enrichment analysis to search pathways related to miRNA targets with strong empirical evidence and FDR < 0.01. Mienturnet tool also allows the graph of the enriched plot, where each dot indicates one miRNA whose targets are participate in the corresponding pathway. The dots are colored according to the p-adjusted value and sized according to the ratio of genes to pathways.
4.6. Statistical Analysis
The descriptive data were presented as means (SD) for continuous variables and percentages for categorical variables. The Shapiro–Wilk test was used to evaluate the distribution of the variables. Demographic and clinical variables were compared between groups using independent sample t-test for continuous variables and Fisher's exact test for categorical variables. Correlation analyses were performed using Spearman’s correlation coefficient to determine the magnitude of association between plasma miRNA levels and biochemical markers. A binary logistic regression analysis was used to evaluate the association between presence of CVD and the variability of the assessed plasma miRNA levels. Multivariate analysis was performed after adjustment for age, sex, and assumption of statin to exclude their effect. Survival data for all-cause of mortality were available for patients, gathered during an average of three years. For each miRNA, Kaplan-Meier survival curves of subjects in the highest quartile were compared to those of remaining quartiles by using the log-rank test. Cox proportional hazard regression analyses were also performed considering age and gender as confounder variables. P < 0.05 was considered statistically significant. All statistical data were analyzed with R (v4.3.3).