1. Introduction
Cardiovascular disease (CVD) remains the leading cause of mortality and morbidity in Europe [1]. Plasmatic low-density lipoprotein (LDLc) levels is the main causal factor of these vascular events. Statins are the cornerstone of atherosclerotic cardiovascular disease (ASCVD) prevention worldwide. Based on the evidence of the clinical trials, the European Society of Cardiology (ESC)/European Atherosclerosis Society (EAS) guidelines have established a goal of 50% reduction in LDLc concentration in high and very high cardiovascular risk patients to reduce vascular events [1,2]. Most of the trials have demonstrated a remarkable reduction of 22% of major vascular events after five years of statin treatment [3]. Even so, almost 80% of treated patients do not achieve recommended LDLc levels. This low adherence to statin therapy is due to several causes, including cultural reasons, costs, side effects and lack of effectiveness information [4].
Although statins are well tolerated, the major cause of non-adherence to treatment is its adverse effects. The European Medicine Agency defines statin intolerance as the inability to tolerate at least two or more statins at the lowest approved daily dose due to the development of side effects, that began or increase during statin therapy and stop when statin treatment was discontinued [5]. Among these side effects, statin-associated muscle symptoms (SAMS) are the most prevalent reason for the non-adherence and/or discontinuation of statin treatment [6]. Four SAMS forms have been described, namely myalgia, myopathy, myositis and rhabdomyolysis. Rhabdomyolysis is the most severe form of statin-induced muscle damage. Nevertheless, myalgia—muscle complaints without creatine kinase elevation or major functional loss—is usually self-reported in the range of 7–29% [7]. The mechanisms involved in the pathogenesis of SAMS are still unknown. The lack of a gold standard test or biomarker to diagnose this adverse effect, the inconsistent definition of this reality in clinical practice, and the impact on the cardiovascular event rates after dropping this therapy demonstrate the importance of identifying a novel and accessible biomarker that distinguishes this population [8].
MicroRNAs (miRNAs) are small non-coding RNAs identified as important regulators of genes involved in several biological processes [9,10,11]. miRNAs regulate gene expression by degrading messenger RNAs (mRNAs), repressing protein synthesis, or interacting with long non-coding RNAs. Their properties make them the most widely studied extracellular RNAs as diagnostic and therapeutic-tailored markers in the cardiovascular field [12,13,14]. Several studies have described the relationship between statins and miRNA profiles [15,16,17,18], but this is the first time miRNAs have been considered as an alternative diagnostic tool related to statin intolerance detection. Therefore, we aim to determine the role of miRNA as a peripheral biomarker to identify the statin intolerant (SI) population and thus personalize dyslipidemia treatment in these patients at cardiovascular risk.
2. Results
2.1. Clinical Parameters between SI and NSI Patients
The anthropometrical and clinical features of the SI and NSI cohorts are shown in Table 1.
This study included 84 patients, 39 SI patients and 45 NSI patients as the control group. Considering the groups, the male population were larger than the female population in NSI, but SI female (61.5%) were significantly more frequent (p = 0.01); in contrast, a lower percentage of females (31%) was observed in the NSI group. There were no significant differences in terms of age. The presence of ASCVD, the years of dyslipidemia (DLP) (p = 0.003, respectively), and high blood pressure (p = 0.006) were significantly different between cohorts. It is important to note that non-HDLc plasmatic concentration (p < 0.001) displayed significant differences between SI and NSI groups. When considering the medication intake, the use of angiotensin agonist receptors (p = 0.004), diuretic (p = 0.003), beta-blockers (p < 0.001), alpha-blockers (p = 0.03) and aspirin (p < 0.001) was significantly different between SI and NSI cohorts. Regarding the lipid-lowering therapy, the NSI cohort was mostly on Rosuvastatin and, to a lesser degree, on Atorvastatin, among others, such as 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase inhibitors or statins. The SI group was treated using alternative lipid-lowering drugs, but nine patients refused to use any other therapeutic option.
2.2. Plasma miRNA Profile in Patients with SI
To assess whether the plasma miRNAs were differentially expressed between the groups, we first conducted screen profiling of 179 circulating miRNAs commonly found in human plasma in ten SI and ten NSI age-matched patients (Figure 1). The established criteria for the selection of miRNA candidates were high expression levels (median Cq < 35 and detected in at least 80% of all samples) and statistical significance (p < 0.05). A total of 10 miRNAs, let-7c-5p, let-7d-5p, let-7f-5p, miR-128-3p, miR-186-5p, miR-30e-3p, miR-376a-3p, miR-376c-3p, miR-543 and miR-574-3p, were significantly expressed according to the selection criteria and selected for further analysis.
2.3. miRNAs Validation Study and Their Correlation with Clinical Parameters
We validated these ten miRNAs in 45 SI and 39 NSI individuals to confirm the diagnostic robustness to discriminate between NSI and SI patients. Our results showed that let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p were significantly upregulated in plasma from SI cohort versus NSI (Figure 2).
We investigated the association between the differentially expressed miRNAs and some significant clinical parameters in the SI group. No correlation was found but DLP and non-HDLc levels (Table 2).
2.4. Circulating miRNA as a Biological Marker of SI
We assess the ability of the differentially expressed circulating miRNAs to distinguish between SI and NSI patients using the AUC-ROC. Figure 3A shows that miR-376c-3p achieved the highest AUC with a value of 0.736 (95% CI: 0.627–0.845; p < 0.001), indicating a moderate performance. Let-7c-5p, let-7d-5p, let-7f-5p and miR-376a-3p display AUC values of 0.652, 0.627, 0.688 and 0.682, respectively. Afterwards, we consider the diagnostic potential of the 5-miRNA set to differentiate between SI and NSI individuals. The AUC for the combination value of our miRNA panel (let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p) was 0.936 (95% CI: 0.887–0.985; p < 0.001) (Figure 3B), improving the diagnostic ability. The sensitivity, specificity and accuracy of each miRNA and the 5-miRNA panel are shown in Table 3.
2.5. Combination of miRNAs, Years of Dislipidemia and Non-HDLc to Categorize SI Patients
We investigated the potential value of the miRNA candidates to discriminate between the SI vs. the NSI cohort and their association with some clinical parameters. Only years of DLP and non-HDLc levels were significantly higher in SI patients (Figure 4A,B). We assess the potential of these two clinical parameters to distinguish between SI and NSI groups. The ROC curves of these factors, DLP and non-HDLc, showed AUC values of 0.700 and 0.807, respectively, indicating a moderate performance to discriminate SI from NSI patients, whereas the combination of these two parameters achieved an AUC value of 0.844 (Table 4). We developed a multivariate model to increase the diagnostic power to discriminate between SI vs. NSI combining differentially expressed miRNAs, DLP and/or non-HDLc plasmatic levels. The diagnostic performance of the 5-miRNA panel was only slightly improved with DLP as, although the AUC value was similar, this model achieved an accuracy of 84.81% against the 82.72% shown by the 5-miRNA panel (Figure 4C,D). In contrast, the diagnostic ability of the 5-miRNA set and the non-HDLc concentration was lower than that of the 5-miRNA panel by itself, showing an AUC value of 0.855. Interestingly, the combination of the 3-miRNA panel composed of let-7f, miR-376a-3p and miR-376c-3p, and DLP plus non-HDLc reached the highest diagnostic performance with an AUC value of 0.954 and an accuracy of 89.47% (Figure 4D and Table 4). We next use the Weka data mining tool to evaluate the performance of our model using 10-fold cross-validation. We ran all classifiers in Weka and, the most successful algorithm was Ada Boost M1 achieving a mean training accuracy, sensitivity, specificity, MCC and AUC of 86%, 83.67%, 88.57%, 0.714 and 0.891, respectively, keeping higher diagnostic value than the 5-miRNA panel model (Figure 5).
2.6. Kyoto Encyclopedia of Genes and Genomes Pathway and Gene Ontology Enrichment Analysis
We next investigated the biological significance of the five differentially expressed miRNAs using GO enrichment and KEGG pathway analysis. We used the miRNet online software (
3. Discussion
We investigated the circulating miRNA signature of a high and very high cardiovascular risk population treated with statins. We report, for the first time, a fingerprint of miRNAs (let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p) to discriminate between SI and NSI patients. Moreover, when we add some clinical information, this circulating miRNA set increases its accuracy remarkably.
Statins are by far the most important therapy used for the treatment and prevention of cardiovascular diseases. Statins have widely proven their profits in primary and secondary cardiovascular prevention [1,2,3]. Although usually well tolerated, statin dropout rates and non-adherence are high after the first year of therapy, 43% and 24%, respectively [21]. The main reason is SAMS which leads to neglect of the treatment and dropout [22] and dramatically increases cardiovascular morbidity and mortality [5,23].
The identification of SI patients is crucial to avoiding cardiovascular risk and reaching the aims established by the guidelines [2]. These patients need to be focused differently in order to prevent the adverse impact of statins and atherosclerosis over the years. However, the current clinical diagnosis of SI has limitations. Clinically, there are no clues to suspect who will present this condition. Myotoxicity shows very heterogeneous manners, and the phenotype that includes myalgia is the most common. SAMS affects symmetrical groups of muscles such as pelvic and shoulder girdles or limbs; muscular weakness, cramps, and, more rarely, increased creatinine kinase plasmatic levels may be considered an adverse effect [24]. Some studies have proposed various scores to assess myalgia [25,26,27]. Several pathological paths have been suggested to explain statins’ toxicity as the mitochondrial dysfunction through Co-enzyme Q10 deficiency, the reduction in beta-oxidation of fatty acids, the promotion of muscular apoptosis, drug–drug interaction, exercise, or certain genetic polymorphisms [22]. The lack of knowledge of the underlying pathological mechanism that leads to SAMS impedes any preventive measures. Thus, the difficulty reaching a clinical diagnosis, the lack of circulating markers or non-invasive tests to detect this side-effect, and the dropout of the preventive therapy leads to a crucial need for a tool that solves this problem. Some groups have related miRNAs with the effectiveness of different statins, and their pleiotropic effect on the endothelial cells [28,29,30,31]. To our knowledge, the use of circulating miRNAs as a diagnostic tool for SI detection remains unexplored.
In the present study, five plasmatic miRNAs, (let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p) were identified as being enriched in SI patients as compared with NSI. Moreover, together, these miRNAs showed a promising potential to discriminate between SI and NSI cohorts. Thus, let-7c-5p is involved in muscle attenuation through the TFG-β factor; its intracellular effector Smad3 is known to inhibit myogenesis and stimulate adipogenesis and myofiber lipid accumulation [32,33]. Let-7f-5p, regulated by the FoxO transcription factors family, is overexpressed in myasthenia gravis in muscle-specific tyrosine kinase antibody seropositive [34]. Myasthenia gravis is a neuromuscular autoimmune disease caused by antibodies which attack receptors at the neuromuscular junction. MiR-376a-3p, proposed as a biomarker of coronary artery disease, acts by nuclear receptor-interacting protein 1 (NRIP1) leading to the cellular apoptosis. Hence, miR-376a-3p might be related to the effect of statins reducing the risk of coronary artery disease, however, this hypothesis needs to be investigated [35]. Related to statin toxicity, others circulating miRNAs have been described as miR-499 and miR-145 [36]. In 2016, Min et al., using combined in vivo and in vitro modelling, demonstrated a release of miR-499-5p linked to muscle injury during exercise when muscular contraction induced by carbachol is combined with statins intake. Circulating levels of miR-1, miR-133a and miR-206 were significantly upregulated in patients in statin therapy who practiced endurance sports and showed muscular adverse effects [37]. Fu et al. demonstrated the statin-induced injury on the skeletal muscle through the overexpression of miR-1a by regulating the mitogen-activated kinase 1 (MKK1) pathway [38]. Some murine models have been used to demonstrate that moderate and gradual exercise may have a beneficial role on statin intake and its impact on skeletal muscle [39,40].
In the biological process analysis, we identified that these miRNAs were highly enriched with genes in the TGF-β pathway, suggesting that they are plausible candidates. Recently, has been reported that TGF-β through Smad3 transcriptional repression leads to the inhibition of myogenesis differentiation [41,42]. Myogenic differentiation and function may need several cellular signal and growth factors, but also an appropriate environment that may be impair by statins.
The genes predicted to be impacted by these five miRNAs are shown in Figure 6. All miRNAs described, except let-7f-5p, shared POTEG, BEND4, POTEM and CASTOR that are genes linked with muscular homeostasis and apoptotic processes. In this sense, CASTOR2, among other genes, has been described as a target of the proteolytic pathway of FoxO driving to muscular atrophy [43]. TGFβR1, common to all miRNAs except miR-376a-3p, has been reported to affect the muscular growth and differentiation [44]. TGFβ signaling plays a critical role in regulating muscle growth and atrophy, in both inherited and acquired myopathies [45]. It is important to note that miRNA target prediction analysis identified insulin-like growth factor (IGF1) as a target gene shared by these five miRNAs. Several studies have reported that statins regulate IGF1-R signaling at different levels [46,47,48]. AKT is a key-regulator of the synthesis and degradation muscle growth and regeneration through AKT/mTOR and FoxO transcription factors [49]. IGF1 is related to statin myotoxicity by inhibiting the IGF1/AKT pathway leading to an increased myofibrillar proteins degradation, impaired proteins synthesis. Thus, statins inhibit AKT phosphorylation activates caspases and PAPR inducing apoptosis and, most importantly, are closely associated with the inhibition of cholesterol synthesis at the level of hydroxymethylglutaryl-CoA (HMG)-CoA reductase. These findings have been demonstrated in both in vitro and in vivo models, enhancing the importance of these pathways in statin-induced myotoxicity [49]. FoxO are a family of transcription factors that are targeted by these miRNAs. They are implicated in the protein degradation and apoptosis. FoxO are related to muscular atrophy and caquexia with the use of certain chemotherapy treatments through the IGF1/AKT/FoxO pathways [50]. Finally, GO enrichment analysis revealed significant enrichment of mitochondrial myopathy-related diseases. In agreement with this, numerous studies have shown that mitochondria play an important role in statin-induced myopathies. Moreover, it has been suggested that statins may have major effects on mitochondrial function, and some of their adverse effects might be mediated through mitochondrial pathways [51,52,53]. Hence, our data suggest that these circulating miRNAs might be involved in the mechanisms of action of statins on mitochondrial function.
In this study we demonstrate that a panel of five miRNAs have a better performance in predicting patients with SI than individual miRNAs. To the best of our knowledge, these five circulating miRNAs have not been reported to be associated with SI. Moreover, we present a multiparametric model that increases the discriminative power of this diagnostic tool in clinical practice. We worked with several clinical parameters and a combination of miRNAs. Thus, our miRNAs set, let-7f-5p, miR-376a-3p and miR-376c-3p, with two clinical parameters such as DLP and non-HDLc concentration showed a high-yield diagnostic accuracy with an AUC of 0.954 and an accuracy of 89.47%. A 10-fold cross-validation evaluation to estimate the prediction performance resulted in high diagnostic accuracy (86%) with an AUC of 0.891 supporting the robustness of the multiparametric model.
Non-HDLc, a reliable marker for cardiovascular risk and residual risk, encompasses all lipoproteins that contain cholesterol and apo B, being closely related to atherogenic dyslipidemia. Thus, in any clinical context associated with insulin resistance, such as diabetes mellitus 2, metabolic syndrome or visceral obesity, it is suggested to use non-HDLc or apo B as a therapeutic aim rather than classic LDLc [54]. This highlights the importance of this non-HDLc for better control of dyslipidemia [2] and an update of the risk calculation indices has recently been published SCORE2 [55] and SCORE2-OP [56]. This update introduces different variations on the original SCORE index assessing not only cardiovascular mortality but the risk of developing a cardiovascular event. In this sense, non-HDLc has been enhanced as a major cardiovascular risk factor with the same impact as LDLc or global cholesterol. Thus, non-HDLc levels are equivalent in terms of risk to high blood pressure, diabetes mellitus or chronic kidney disease. On the other hand, DLP play a key role in ASCVD. The time of exposition to a reduced LDLc concentration is independent of the reduced mechanism used; the clinical benefit to more decreased LDLc levels, whether genetically or pharmacologically, is determinant; the lower and the earlier the better [57].
We expose a circulating miRNAs fingerprint as an accessible and helpful tool for SI diagnosis. In addition, we propose a multiparametric model, a panel of three miRNAs and two clinical parameters, that detects patients with this adverse event with the highest diagnostic value reported to date. The bioinformatics analysis leads to the biological and molecular processes present in this entity. This proposal permits a therapy-tailored strategy to avoid the dropout of a crucial therapy and offer an alternative to treat patients with high cardiovascular risk.
Our current study has several limitations. Firstly, our sample was prospectively recruited from the outpatient clinic. The size of the study sample, recruited from various Lipid Units and Cardiology Departments, did not allow us to obtain a robust multivariate logistic regression model. Furthermore, a larger sample size is needed to validate these data. As a consequence, these results should be extended and replicated to a larger population before the novel biomarkers can be routinely applied in clinical practice. Finally, even though databases registered the expression of the parental genes, we have no confirmation about the direct secretion of these circulating miRNAs into the extracellular space. Hence, the association of miRNAs with SI and all the interactions are putative. Therefore, larger confirmative studies on circulating miRNAs are needed to generate a better, reproducible prognostic signature for SI.
4. Materials and Methods
4.1. Study Population and Design
This is a multicenter prospective case-control study. Patients were recruited from three field centers (Puerta del Mar University Hospital, Cádiz; Virgen del Rocio University Hospital, Sevilla; and Reina Sofia University Hospital of Córdoba, Spain). Inclusion criteria were patients over 18 years old with a high and very high cardiovascular risk in treatment with statins referred to the Lipid Unit of these centers. We have focused on the most frequent symptom, SAMS [5]. A total of 84 consecutive subjects (Figure 7) were included and divided into two different cohorts: (i) SI patients (n = 39) and, (ii) patients non-intolerants to statins (NSI) (n = 45). SI group was defined based on the European Medicine Agency criteria [5]. Before including these patients in the study, we applied the 2019 ESC/EAS Lipid Management Guidelines, which establish the requirement to offer the patient alternative dosing such as every other day or twice a week with Atorvastatin, Pitavastatin or Rosuvastatin. All patients included dropped out of the alternative proposed due to intolerance to statins. All patients recruited were high or very high cardiovascular risk based on the classification of the cardiovascular risk grade of the ESC/EAS guidelines (Table 5) [2]. ASCVD can be documented as either clinical or unequivocal based on imaging. In the clinical branch, ASCVD includes previous acute coronary syndrome (myocardial infarction or unstable angina), stable angina, coronary revascularization procedures including percutaneous coronary intervention, coronary artery bypass graft or arterial revascularization procedures, ischemic stroke, and peripheral arterial disease. Conditions considered to explicitly document ASCVD on imaging include findings known to be predictive of clinical events, such as significant plaque on coronary angiography or coronary tomography scan (a multivessel coronary disease with two major epicardial arteries having >50% stenosis) or on carotid ultrasound.
Exclusion criteria were patients over 80 years old or with marked frailty, intensive physical activity, heavy alcohol consumption and patients who presented chronic inflammatory, autoimmune or neoplasia disease or other pharmacological therapies that interfere with the metabolic pathways of the different statins. Subjects with high training exercise workload, previous muscular symptoms, myopathy disease, high levels of creatine kinase levels or hypothyroidism were excluded from the study.
Detailed anthropometric, clinical and pharmacological information was obtained from each subject including cardiovascular risk factors and previous vascular events. All patients performed a blood test including thyroid, hepatic, kidney, lipid profiles and creatine kinase at the inclusion moment. Our institution’s ethics committee (Comité de Ética de la Investigation de Cádiz) approved the study protocol. The study was performed in full compliance with the Helsinki II Declaration. All participants provided written informed consent.
4.2. Blood Collection
Ten milliliters of peripheral blood were collected in K2-ethylenediaminetetraacetic acid tubes after 10 h overnight fasting and were immediately centrifuged (1500× g, 15 min, 4 °C). The blood was processed within 4 h after isolation. The upper layer containing plasma was divided into aliquots and stored at −80 °C until further analysis.
4.3. RNA Isolation
Total RNA was extracted from 200 µL of plasma using the miRNeasy Serum/Plasma Advanced Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. A total of 200 µL of plasma aliquots were used to obtain enough final volume. During the extraction, 3.5 μL of miRNeasy Serum/Plasma Spike-In Control (1.6 × 108 copies/μL of the C. elegans miR-39 miRNA mimic) were added to each sample as an internal control. RNA was eluted in 20 μL of RNase-free water.
4.4. MiRNA Real-Time Reverse Transcriptase-Polymerase Chain Reaction
To evaluate miRNA expression levels and profiles of SI and NSI groups, an analysis of 179 miRNA species known to be present in human serum was performed using miRCURY LNA RT Kit (Qiagen, Hilden, Germany) for reverse transcription, miRCURY LNA SYBR Green PCR Kit (Qiagen) for quantitative real-time (qRT-PCR) amplification and Human Serum/Plasma Focus, miRCURY LNA miRNA Focus PCR panels (Quiagen, Hilden, Germany) according to the manufacturer’s instructions. Raw cycle threshold (Cq) values were inter-plate calibrated using UniSp3. Cqs above 35 cycles were censored at the minimum level observed for each miRNA. As housekeeping miRNA for data normalization, we selected the most stably expressed pair of miRNAs (miR-148a-3p and let-7b-5p) as determined by the Normfinder algorithm [58]. Accordingly, each single miRNA was normalized as ∆Cq = mean (CqmiR-148a-3p and let-7b-5p)—Cq miRNA. MiRNA levels were log-transformed before being used in the statistical analyses.
4.5. Validation of miRNA Profiles
For the validation study, each selected miRNA candidate was quantified using miRCURY LNA miRNA Custom PCR Panels (Qiagen, Hilden, Germany). RNA was reverse transcribed using the miRCURY LNA RT Kit (Qiagen, Hilden, Germany). qRT-PCR was performed using miRCURY LNA SYBR Green PCR Kit (Qiagen, Hilden, Germany), as previously described [59] and amplification curves were evaluated with CFX Manager™ software (BioRad). The specificity of the amplification was corroborated by melting curve analysis.
4.6. miRNA-Gene Network Analysis
The miRNAs obtained were tested using the miRNet database (
4.7. KEGG and GO Term Enrichment Analysis
The miRNet online software (
4.8. Statistical Analysis
Continuous variables are shown as mean ± standard deviation. Categorical variables are expressed as frequency and percentage of patients (%). Outliers were identified through the Rout method, using a Q = 1% [64]. The normal distribution of each variable was verified with the Shapiro–Wilk test. Intergroup comparisons of miRNAs levels were performed using non-parametric Mann–Whitney and Kruskal–Wallis rank tests for continuous variables. An analysis of differences between groups was performed using the analysis of variance. ROC curves that characterize the diagnostic performance of candidate miRNAs and logistic regression models were plotted to determine the area under the curve (AUC) and the specificity and sensitivity of the optimal cut-offs. ROC curves were generated by plotting sensitivity against 1-specificity. Data were presented as the AUC and 95% confidence intervals. The changes in p-values of their variables were evaluated by the Wald test and a likelihood ratio. Training was performed with 10-fold cross-validation using Waikato environment for knowledge analysis (Weka) data mining tool. The statistical software package R was used for all analyses (Team RC. R: A Language and Environment for Statistical Computing.
5. Conclusions
In conclusion, we identified for the first time a signature of five miRNAs differentially expressed in the SI population. Our multiparametric model also combined circulating miRNAs and clinical variables, and emerges as a potential clinical tool to identify patients with muscular toxic effects as a result of statins. In daily practice, this panel with accessible clinical parameters may improve the accuracy and may lead to a tailored-treatment strategy improving the cardiovascular risk and outcome of these patients.
R.T. and A.M. conceived the experiments; A.M., O.M., F.F., A.G.-E., J.S.R.R., O.C., A.P.-S. and E.A.-V. recruited the subjects; F.B. and A.P.-S. conducted the experiments, and analysed the results; R.T., A.M., F.B. and A.P.-S. wrote the manuscript. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.
The study was conducted in accordance with the Declaration of Helsinki and approved by the Andalusian Ethic Committee Board (protocol code FIB-ATO 2018-01; date of approval 3 May 2010).
Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients for paper publication.
Data transparency is guaranteed. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. We used various software for functional enrichment and statistical analysis. All of them are cited in our manuscript.
The authors declare no conflict of interest.
The study protocol was approved by the Andalusian Biomedical Research Ethics committee. The study was performed in full compliance with the Declaration of Helsinki.
ACEI | Angiotensin-converting enzyme inhibitors |
ARB | Angiotensin II receptors blockers |
ASCVD | Atherosclerotic cardiovascular diseases |
AUC | Area under the ROC curve |
CCB | Calcium channel blockers |
CI | Confidence interval |
CPK | Creatine kinase |
CVC | Cardiovascular disease |
DLP | Years of dyslipidemia |
ESC | European Society of Cardiology |
EAS | European Atherosclerosis Society |
FoxO | Forkhead members of the O class family |
GO | Gene Ontology |
HDLc | High-density lipoprotein cholesterol |
KEGG | Kyoto Encyclopedia of Genes and Genomes |
LDLc | Plasmatic low-density lipid |
miRNAs | MicroRNAs |
Non-HDLc | Non-HDL cholesterol |
NSI | Non-statin intolerant |
SAMS | Statin-associated muscle symptoms |
SD | Standard deviation |
SI | Statin intolerant |
STRING | Search Tool for the Retrieval of Interacting Genes |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. Color heatmap based on raw miRNA expression values where each column represents a patient, and each row represents a miRNA. The color scale illustrates the relative expression level of miRNAs (red and yellow represent low expression and blue and purple represent high expression). MiRNA expression levels were normalized to miR-148a-3p and let-7b-5p. MiRNA: microRNA; NSI: non-statin intolerant; SI: statin intolerant.
Figure 2. Boxplots of miRNA expression levels in NSI and SI cohorts. The analysis was carried out using qPCR. Data are presented in log2. Data represent the mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.005. Error bars represent SDs. NSI, non-statin intolerant; SI: statin intolerant.
Figure 3. ROC curves for evaluating the predictive performance of differentially expressed miRNAs to discriminate between SI vs. NSI. (A) ROC curves for let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p. (B) The ROC curve of the 5-miRNA panel combination value of let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p. AUC: area under the curve; miRNA: microRNA; NSI: non-statin intolerant; SI: statin intolerant.
Figure 4. ROC curves for evaluating the predictive performance of clinical factors with differentially expressed miRNAs. (A) Box plot of DLP in NSI (n = 45) and SI (n = 39) subjects. (B) Box plot of non-HDLc levels in NSI (n = 45) and SI (n = 39) subjects. (C) ROC curves for each clinical parameter, DLP and non-HDLc, and for the association of DLP plus non-HDLc. (D) The ROC curve of the combined value of the 3-miRNA panel (let-7f-5p, miR-376a-3p and miR-376c-3p), DLP and non-HDLc plasmatic concentration. DLP: years of dyslipidemia SI: statin intolerant; NSI: non-statin intolerant; non-HDLc: non-high-density lipoprotein cholesterol. * p < 0.05; *** p < 0.005.
Figure 5. ROC curve of 10-fold cross validation test for 3-miRNA panel + DLP + non-HDLc model (Ada Boost M1).
Figure 6. KEGG and GO analysis of differentially expressed miRNAs. (A) miRNA-gene network for of let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p. (B) Venn diagram showing overlap of gene targets of let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p. (C) GO and KEGG functional enrichment analysis of let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p.
Figure 7. The flowchart of the study design. This figure illustrates the experimental workflow of the study including screening and validation. MiRNA: microRNA; NSI: non-statin intolerant; SI: statin intolerant.
Baseline demographics, clinical characteristics, and treatment of NSI and SI population. Data are presented as mean ± SD for continuous variables and as percentage for categorical variables. The difference between NSI and SI patients was evaluated with unpaired Student t test a, Pearson Chi-square b and Wilcoxon test c. * Therapy mostly in association. NS, no significant; ACEI: angiotensin-converting enzyme inhibitors; ARB: angiotensin II receptor blockers; ASCVD: atherosclerotic cardiovascular disease; CCB: calcium channel blockers; CPK: creatine kinase; DLP: years of dyslipidemia; MDRD-4: glomerular filtration rates; non-HDLc: non-high-density lipoprotein cholesterol; NSI: non-statin intolerant; OAD: oral antidiabetic drugs; PSCK-9: Monoclonal anti-proprotein convertase subtilisin/kexin type 9; SI: statin intolerant.
Variables | NSI | SI | p Value |
---|---|---|---|
n | 45 | 39 | |
Demographics | |||
Age (years) a | 66.6 ± 11.5 | 63.6 ± 10.9 | NS |
Sex (female, %) b | 31 | 61.5 | 0.01 |
DLP (years) c | 5.6 ± 6.3 | 9.1 ± 7.5 | 0.003 |
High blood pressure (%) b | 73 | 41 | 0.006 |
Diseases | |||
Diabetes Mellitus (%) b | 38 | 20.5 | NS |
ASCVD (%) b | 71 | 13 | <0.001 |
Chronic kidney disease (%) b | 18 | 12.8 | NS |
Analytical profile | |||
Basal blood glucose a | 118.5 ± 36.2 | 110 ± 44.6 | NS |
Non-HDLc (mg/dL) a | 104.9 ± 32.7 | 169.9 ± 63.5 | <0.001 |
Triglycerides (mg/dL) a | 151.7 ± 97.4 | 164 ± 93 | NS |
MDRD-4 (mL/min) a | 76.5 ± 25.6 | 82 ± 31.6 | NS |
Transaminase GOT (U/L) a | 23.1 ± 13 | 23.5 ± 8 | NS |
Transaminase GPT (U/L) a | 27.3 ± 23 | 24.6 ± 18.9 | NS |
CPK (U/L) a | 82.8 ± 40.4 | 165.9 ± 14.7 | NS |
Medication | |||
ACEI (%) b | 11 | 2.5 | NS |
ARB (%) b | 67 | 33.3 | 0.004 |
OAD (%) b | 33 | 18 | NS |
Insulin (%) b | 11 | 10.5 | NS |
Diuretic (%) b | 51 | 18 | 0.003 |
CCB (%) b | 31 | 18 | NS |
Beta-blockers (%) b | 62 | 15.4 | <0.001 |
Alpha-blockers (%) b | 22 | 5 | 0.03 |
Aspirin (%) b | 64 | 18 | <0.001 |
Atorvastatin 40 mg (%) | 20 | - | |
Atorvastatin 80 mg (%) | 15.5 | - | |
Rosuvastatin 10 mg (%) | 17.7 | - | |
Rosuvastatin 20 mg (%) | 35.5 | - | |
Pitavastatin 4mg (%) | 6.6 | - | |
Simvastatin 40 mg (%) | 4.4 | - | |
PCSK9 inhibitors: Evolucumab (%) b | 13 | 7.6 | NS |
PCSK9 inhibitors Alirocumab (%) b | 11 | 2.5 | NS |
Fenofibrate (%) b | 18 | 7.6 | NS |
Omega-3 (%) | - | 5.1 | |
Colesevelam (%) | - | 2.5 | |
Colestiramine (%) | - | 12.8 | |
Armolipid plus * (%) | - | 41.0 | |
Ezetimibe 10 mg * (%) b | 42 | 38 | NS |
Acenocumarol (%) b | 7 | 7.6 | NS |
Correlation between the clinical parameter (DLP and non-HDLc) and individual microRNAs in SI subjects. DLP: years of dyslipidemia; non-HDLc: non-high-density lipoprotein cholesterol SI: statin intolerant. Coefficient significant at p < 0.05.
microRNAs | SI Cohort | |||
---|---|---|---|---|
DLP (Years) | Non-HDLc (mg/dL) | |||
Pearson r | p | Pearson r | p | |
Let-7c-5p | −0.164 | 0.281 | 0.177 | 0.244 |
Let-7d-5p | −0.079 | 0.603 | 0.226 | 0.131 |
Let-7f-5p | −0.136 | 0.368 | 0.168 | 0.266 |
miR-376a-3p | −0.173 | 0.250 | 0.103 | 0.498 |
miR-376c-3p | −0.265 | 0.079 | 0.176 | 0.248 |
Assessment of the potential diagnostic value of differentially expressed miRNAs and the 5-miRNA panel (let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p) as biomarkers to categorize statin intolerant patients. AUC: area under the curve; CI: confidence interval; miRNA: microRNA.
miRNA | AUC (95% CI) | Sensitivity % | Specificity % | Accuracy % | p Value |
---|---|---|---|---|---|
Let-7c-5p | 0.652 (0.535 to 0.770) | 61.70 | 55.56 | 59.04 | 0.017 |
Let-7d-5p | 0.627 (0.507 to 0.747) | 52.63 | 58.70 | 55.95 | 0.046 |
Let-7f-5p | 0.688 (0.573 to 0.803) | 60.53 | 64.44 | 62.65 | 0.003 |
miR-376a-3p | 0.682 (0.563 to 0.800) | 68.89 | 64.10 | 66.67 | 0.004 |
miR-376c-3p | 0.736 (0.627 to 0.845) | 70.45 | 64.10 | 67.47 | <0.001 |
5-miRNA panel | 0.936 (0.887 to 0.985) | 81.25 | 84.85 | 82.72 | <0.001 |
Evaluation of the potential of clinical parameters (DLP and non-HDLc) and the multivariate models as SI biomarkers. AUC, area under the curve; CI, confidence interval. DLP: years of dislipidemia; miRNA: microRNA; non-HDLc: non-high-density lipoprotein cholesterol.
Multiparametric Model | AUC (95% CI) | Sensitivity % | Specificity % | Accuracy % | p Value |
---|---|---|---|---|---|
DLP (years) | 0.700 (0.587 to 0.814) | 57.89 | 60.00 | 58.54 | 0.017 |
Non-HDLc (mg/dL) | 0.807 (0.703 to 0.911) | 77.08 | 84.38 | 80.00 | <0.001 |
DLP + non-HDLc | 0.844 (0.751 to 0.937) | 79.55 | 85.29 | 82.05 | <0.001 |
5-miRNA panel + DLP | 0.940 (0.892 to 0.989) | 85.71 | 83.78 | 84.81 | <0.001 |
5-miRNA panel + non-HDLc | 0.889 (0.814 to 0.964) | 85.00 | 81.08 | 83.12 | <0.001 |
3-miRNA panel + DLP + non-HDLc | 0.954 (0.911 to 0.998) | 89.74 | 89.19 | 89.47 | <0.001 |
Cardiovascular risk criteria based on the 2019 ESC/EAS guidelines used to the recruitment of patients [
Very High CVD Risk | High CVD Risk |
---|---|
Presence of ASCVD clinically/imaging. | Total cholesterol over 310 mg/dL, |
DM patients with target organ damage or |
DM patients without target organ damage. |
Severe CKD |
Moderate CKD |
A calculated SCORE ≥ 10% for 10 years risk of fatal CVD. | A calculated SCORE ≥ 5% and < 10% for 10 years’ risk of fatal CVD. |
FH with a ASCVD or with another major risk factor. | FH without any other major risk factor. |
References
1. Catapano, A.L.; Graham, I.; De Backer, G.; Wiklund, O.; Chapman, M.J.; Drexel, H.; Hoes, A.W.; Jennings, C.S.; Landmesser, U.; Pedersen, T.R. et al. 2016 ESC/EAS Guidelines for the Management of Dyslipidaemias. Eur. Heart J.; 2016; 37, pp. 2999-3058. [DOI: https://dx.doi.org/10.1093/eurheartj/ehw272] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27567407]
2. Mach, F.; Baigent, C.; Catapano, A.L.; Koskinas, K.C.; Casula, M.; Badimon, L.; Chapman, M.J.; De Backer, G.G.; Delgado, V.; Ference, B.A. et al. ESC Scientific Document Group (2020). 2019 ESC/EAS Guidelines for the management of dyslipidaemias: Lipid modification to reduce cardiovascular risk. Eur. Heart J.; 2020; 41, pp. 111-188. [DOI: https://dx.doi.org/10.1093/eurheartj/ehz455] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31504418]
3. Cholesterol Treatment Trialists’ (CTT) Collaboration Baigent, C.; Blackwell, L.; Emberson, J.; Holland, L.E.; Reith, C.; Bhala, N.; Peto, R.; Barnes, E.H.; Keech, A. et al. Efficacy and safety of more intensive lowering of LDL cholesterol: A meta-analysis of data from 170,000 participants in 26 randomised trials. Lancet; 2010; 376, pp. 1670-1681. [DOI: https://dx.doi.org/10.1016/S0140-6736(10)61350-5]
4. Banach, M.; Stulc, T.; Dent, R.; Toth, P.P. Statin non-adherence and residual cardiovascular risk: There is need for substantial improvement. Int. J. Cardiol.; 2016; 225, pp. 184-196. [DOI: https://dx.doi.org/10.1016/j.ijcard.2016.09.075] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27728862]
5. Rosenson, R.S.; Baker, S.K.; Jacobson, T.A.; Kopecky, S.L.; Parker, B.A. The National Lipid Association’s Muscle Safety Expert Panel. An assessment by the Statin Muscle Safety Task Force: 2014 update. J. Clin. Lipidol.; 2014; 8, pp. S58-S71. [DOI: https://dx.doi.org/10.1016/j.jacl.2014.03.004]
6. Toth, P.P.; Granowitz, C.; Hull, M.; Anderson, A.; Philip, S. Long-term statin persistence is poor among high-risk patients with dyslipidemia: A real-world administrative claims analysis. Lipids Health Dis.; 2019; 18, 175. [DOI: https://dx.doi.org/10.1186/s12944-019-1099-z]
7. Stroes, E.S.; Thompson, P.D.; Corsini, A.; Vladutiu, G.D.; Raal, F.J.; Ray, K.K.; Roden, M.; Stein, E.; Tokgözoğlu, L.; Nordestgaard, B.G. et al. European Atherosclerosis Society Consensus Panel (2015). Statin-associated muscle symptoms: Impact on statin therapy-European Atherosclerosis Society Consensus Panel Statement on Assessment, Aetiology and Management. Eur. Heart J.; 2015; 36, pp. 1012-1022. [DOI: https://dx.doi.org/10.1093/eurheartj/ehv043]
8. Zhang, H.; Plutzky, J.; Shubina, M.; Turchin, A. Continued Statin Prescriptions after Adverse Reactions and Patient Outcomes: A Cohort Study. Ann. Intern. Med.; 2017; 167, pp. 221-227. [DOI: https://dx.doi.org/10.7326/M16-0838]
9. Henning, R.J. Cardiovascular Exosomes and MicroRNAs in Cardiovascular Physiology and Pathophysiology. J. Cardiovasc. Transl. Res.; 2021; 14, pp. 195-212. [DOI: https://dx.doi.org/10.1007/s12265-020-10040-5]
10. Mellis, D.; Caporali, A. MicroRNA-based therapeutics in cardiovascular disease: Screening and delivery to the target. Biochem. Soc. Trans.; 2018; 46, pp. 11-21. [DOI: https://dx.doi.org/10.1042/BST20170037]
11. Kalayinia, S.; Arjmand, F.; Maleki, M.; Malakootian, M.; Singh, C.P. MicroRNAs: Roles in cardiovascular development and disease. Cardiovasc. Pathol.; 2021; 50, 107296. [DOI: https://dx.doi.org/10.1016/j.carpath.2020.107296] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33022373]
12. Costa, M.C.; Kurc, S.; Drożdż, A.; Cortez-Dias, N.; Enguita, F.J. The circulating non-coding RNA landscape for biomarker research: Lessons and prospects from cardiovascular diseases. Acta Pharmacol. Sin.; 2018; 39, pp. 1085-1099. [DOI: https://dx.doi.org/10.1038/aps.2018.35]
13. Lu, D.; Thum, T. RNA-based diagnostic and therapeutic strategies for cardiovascular disease. Nat. Rev. Cardiol.; 2019; 16, pp. 661-674. [DOI: https://dx.doi.org/10.1038/s41569-019-0218-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31186539]
14. Huang, C.K.; Kafert-Kasting, S.; Thum, T. Preclinical and Clinical Development of Noncoding RNA Therapeutics for Cardiovascular Disease. Circ. Res.; 2020; 126, pp. 663-678. [DOI: https://dx.doi.org/10.1161/CIRCRESAHA.119.315856]
15. Sodi, R.; Eastwood, J.; Caslake, M.; Packard, C.J.; Denby, L. Relationship between circulating microRNA-30c with total- and LDL-cholesterol, their circulatory transportation and effect of statins. Clin. Chim. Acta.; 2017; 466, pp. 13-19. [DOI: https://dx.doi.org/10.1016/j.cca.2016.12.031]
16. Zambrano, T.; Hirata, R.; Hirata, M.H.; Cerda, Á.; Salazar, L.A. Statins differentially modulate microRNAs expression in peripheral cells of hyperlipidemic subjects: A pilot study. Eur. J. Pharm. Sci.; 2018; 117, pp. 55-61. [DOI: https://dx.doi.org/10.1016/j.ejps.2018.02.007]
17. Stormo, C.; Kringen, M.K.; Lyle, R.; Olstad, O.K.; Sachse, D.; Berg, J.P.; Piehler, A.P. RNA-sequencing analysis of HepG2 cells treated with atorvastatin. PLoS ONE; 2014; 9, e105836. [DOI: https://dx.doi.org/10.1371/journal.pone.0105836]
18. Cerda, A.; Bortolin, R.H.; Manriquez, V.; Salazar, L.; Zambrano, T.; Fajardo, C.M.; Hirata, M.H.; Hirata, R. Effect of statins on lipid metabolism-related microRNA expression in HepG2 cells. Pharmacol. Rep.; 2021; 73, pp. 868-880. [DOI: https://dx.doi.org/10.1007/s43440-021-00241-3]
19. Lee, S.; Dong, H.H. FoxO integration of insulin signaling with glucose and lipid metabolism. J. Endocrinol.; 2017; 233, pp. R67-R79. [DOI: https://dx.doi.org/10.1530/JOE-17-0002]
20. Mallinson, J.E.; Constantin-Teodosiu, D.; Sidaway, J.; Westwood, F.R.; Greenhaff, P.L. Blunted Akt/FOXO signalling and activation of genes controlling atrophy and fuel use in statin myopathy. J. Physiol.; 2009; 587, pp. 219-230. [DOI: https://dx.doi.org/10.1113/jphysiol.2008.164699]
21. Naderi, S.H.; Bestwick, J.P.; Wald, D.S. Adherence to drugs that prevent cardiovascular disease: Meta-analysis on 376,162 patients. Am. J. Med.; 2012; 125, pp. 882-887.e1. [DOI: https://dx.doi.org/10.1016/j.amjmed.2011.12.013] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22748400]
22. Vinci, P.; Panizon, E.; Tosoni, L.M.; Cerrato, C.; Pellicori, F.; Mearelli, F.; Biasinutto, C.; Fiotti, N.; Di Girolamo, F.G.; Biolo, G. Statin-Associated Myopathy: Emphasis on Mechanisms and Targeted Therapy. Int. J. Mol. Sci.; 2021; 22, 11687. [DOI: https://dx.doi.org/10.3390/ijms222111687] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34769118]
23. Turner, R.M.; Yin, P.; Hanson, A.; FitzGerald, R.; Morris, A.P.; Stables, R.H.; Jorgensen, A.L.; Pirmohamed, M.J. Investigating the prevalence, predictors, and prognosis of suboptimal statin use early after a non-ST elevation acute coronary syndrome. Clin. Lipidol.; 2017; 11, pp. 204-214. [DOI: https://dx.doi.org/10.1016/j.jacl.2016.12.007]
24. Alonso, R.; Cuevas, A.; Cafferata, A. Diagnosis and Management of Statin Intolerance. J. Atheroscler. Thromb.; 2019; 26, pp. 207-215. [DOI: https://dx.doi.org/10.5551/jat.RV17030]
25. Parker, B.A.; Capizzi, J.A.; Grimaldi, A.S.; Clarkson, P.M.; Cole, S.M.; Keadle, J.; Chipkin, S.; Pescatello, L.S.; Simpson, K.; White, C.M. et al. Effect of statins on skeletal muscle function. Circulation; 2013; 127, pp. 96-103. [DOI: https://dx.doi.org/10.1161/CIRCULATIONAHA.112.136101]
26. Colantonio, L.D.; Kent, S.T.; Huang, L.; Chen, L.; Monda, K.L.; Serban, M.C.; Manthripragada, A.; Kilgore, M.L.; Rosenson, R.S.; Muntner, P. Algorithms to identify statin intolerance in Medicare administrative claim data. Cardiovasc. Drugs Ther.; 2016; 30, pp. 525-533. [DOI: https://dx.doi.org/10.1007/s10557-016-6680-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27389413]
27. Mancini, G.B.; Baker, S.; Bergeron, J.; Fitchett, D.; Frohlich, J.; Genest, J.; Gupta, M.; Hegele, R.A.; Ng, D.; Pearson, G.J. et al. Diagnosis, Prevention, and Management of Statin Adverse Effects and Intolerance: Canadian Consensus Working Group Update. Can. J. Cardiol.; 2016; 32, pp. S35-S65. [DOI: https://dx.doi.org/10.1016/j.cjca.2016.01.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27342697]
28. Leal, K.; Saavedra, K.; Rebolledo, C.; Salazar, L.A. MicroRNAs hsa-miR-618 and hsa-miR-297 Might Modulate the Pleiotropic Effects Exerted by Statins in Endothelial Cells through the Inhibition of ROCK2 Kinase: In-silico Approach. Front. Cardiovasc. Med.; 2021; 8, 704175. [DOI: https://dx.doi.org/10.3389/fcvm.2021.704175]
29. Raucci, A.; Macrì, F.; Castiglione, S.; Badi, I.; Vinci, M.C.; Zuccolo, E. MicroRNA-34a: The bad guy in age-related vascular diseases. Cell. Mol. Life Sci.; 2021; 78, pp. 7355-7378. [DOI: https://dx.doi.org/10.1007/s00018-021-03979-4]
30. Hirota, T.; Fujita, Y.; Ieiri, I. An updated review of pharmacokinetic drug interactions and pharmacogenetics of statins. Expert Opin. Drug Metab. Toxicol.; 2020; 16, pp. 809-822. [DOI: https://dx.doi.org/10.1080/17425255.2020.1801634]
31. Mohammadzadeh, N.; Montecucco, F.; Carbone, F.; Xu, S.; Al-Rasadi, K.; Sahebkar, A. Statins: Epidrugs with effects on endothelial health?. Eur. J. Clin. Investig.; 2020; 50, e13388. [DOI: https://dx.doi.org/10.1111/eci.13388] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32854143]
32. Mitchell, C.J.; D’Souza, R.F.; Schierding, W.; Zeng, N.; Ramzan, F.; O’Sullivan, J.M.; Poppitt, S.D.; Cameron-Smith, D. Identification of human skeletal muscle miRNA related to strength by high-throughput sequencing. Physiol. Genom.; 2018; 50, pp. 416-424. [DOI: https://dx.doi.org/10.1152/physiolgenomics.00112.2017]
33. Liu, D.; Black, B.L.; Derynck, R. TGF-beta inhibits muscle differentiation through functional repression of myogenic transcription factors by Smad3. Genes. Dev.; 2001; 15, pp. 2950-2966. [DOI: https://dx.doi.org/10.1101/gad.925901] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11711431]
34. Fiorillo, A.A.; Heier, C.R.; Huang, Y.F.; Tully, C.B.; Punga, T.; Punga, A.R. Estrogen Receptor, Inflammatory, and FOXO Transcription Factors Regulate Expression of Myasthenia Gravis-Associated Circulating microRNAs. Front. Immunol.; 2020; 11, 151. [DOI: https://dx.doi.org/10.3389/fimmu.2020.00151] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32153563]
35. Du, L.; Xu, Z.; Wang, X.; Liu, F. Integrated bioinformatics analysis identifies microRNA-376a-3p as a new microRNA biomarker in patient with coronary artery disease. Am. J. Transl. Res.; 2020; 12, pp. 633-648. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32194911]
36. Min, P.K.; Park, J.; Isaacs, S.; Taylor, B.A.; Thompson, P.D.; Troyanos, C.; D’Hemecourt, P.; Dyer, S.; Chan, S.Y.; Baggish, A.L. Influence of statins on distinct circulating microRNAs during prolonged aerobic exercise. J. Appl. Physiol.; 2016; 120, pp. 711-720. [DOI: https://dx.doi.org/10.1152/japplphysiol.00654.2015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26472872]
37. Siracusa, J.; Koulmann, N.; Sourdrille, A.; Chapus, C.; Verret, C.; Bourdon, S.; Goriot, M.E.; Banzet, S. Phenotype-Specific Response of Circulating miRNAs Provides New Biomarkers of Slow or Fast Muscle Damage. Front. Physiol.; 2018; 9, 684. [DOI: https://dx.doi.org/10.3389/fphys.2018.00684]
38. Fu, C.N.; Song, J.W.; Song, Z.P.; Wang, Q.W.; Bai, W.W.; Guo, T.; Li, P.; Liu, C.; Wang, S.X.; Dong, B. Excessive expression of miR-1a by statin causes skeletal injury through targeting mitogen-activated protein kinase kinase kinase 1. Aging; 2021; 13, pp. 11470-11490. [DOI: https://dx.doi.org/10.18632/aging.202839]
39. Bouitbir, J.; Daussin, F.; Charles, A.L.; Rasseneur, L.; Dufour, S.; Richard, R.; Piquard, F.; Geny, B.; Zoll, J. Mitochondria of trained skeletal muscle are protected from deleterious effects of statins. Muscle. Nerve.; 2012; 46, pp. 367-373. [DOI: https://dx.doi.org/10.1002/mus.23309]
40. Lotteau, S.; Ivarsson, N.; Yang, Z.; Restagno, D.; Colyer, J.; Hopkins, P.; Weightman, A.; Himori, K.; Yamada, T.; Bruton, J. et al. A Mechanism for Statin-Induced Susceptibility to Myopathy. JACC Basic Transl. Sci.; 2019; 4, pp. 509-523. [DOI: https://dx.doi.org/10.1016/j.jacbts.2019.03.012]
41. Liu, L.; Cao, P.; Zhang, L.; Qi, M.; Wang, L.; Li, Z.; Shao, G.; Ding, L.; Zhao, X.; Zhao, X. et al. Comparisons of adipogenesis- and lipid metabolism-related gene expression levels in muscle, adipose tissue and liver from Wagyu-cross and Holstein steers. PLoS ONE; 2021; 16, e0247559. [DOI: https://dx.doi.org/10.1371/journal.pone.0247559]
42. Bai, W.; Ren, M.; Cheng, W.; Lu, X.; Liu, D.; Wang, B. Qindan Capsule Attenuates Myocardial Hypertrophy and Fibrosis in Pressure Overload-Induced Mice Involving mTOR and TGF-β1/Smad Signaling Pathway Inhibition. Evid. Based Complementary Altern. Med.; 2021; 2021, 5577875. [DOI: https://dx.doi.org/10.1155/2021/5577875] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34007292]
43. Oyabu, M.; Takigawa, K.; Mizutani, S.; Hatazawa, Y.; Fujita, M.; Ohira, Y.; Sugimoto, T.; Suzuki, O.; Tsuchiya, K.; Suganami, T. et al. FOXO1 cooperates with C/EBPδ and ATF4 to regulate skeletal muscle atrophy transcriptional program during fasting. FASEB J.; 2022; 36, e22152. [DOI: https://dx.doi.org/10.1096/fj.202101385RR] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35061305]
44. Hillege, M.; Shi, A.; Galli, R.A.; Wu, G.; Bertolino, P.; Hoogaars, W.; Jaspers, R.T. Lack of Tgfbr1 and Acvr1b synergistically stimulates myofibre hypertrophy and accelerates muscle regeneration. Elife; 2022; 11, e77610. [DOI: https://dx.doi.org/10.7554/eLife.77610] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35323108]
45. Burks, T.N.; Cohn, R.D. Role of TGF-beta signalling in inherited and acquired myopathies. Skelet. Muscle; 2011; 1, 19. [DOI: https://dx.doi.org/10.1186/2044-5040-1-19]
46. Forbes, K.; Shah, V.K.; Siddals, K.; Gibson, J.M.; Aplin, J.; Westwood, M. Statins inhibit insulin-like growth factor action in first trimester placenta by altering insulin-like growth factor 1 receptor glycosylation. Mol. Hum. Reprod.; 2015; 21, pp. 105-114. [DOI: https://dx.doi.org/10.1093/molehr/gau093]
47. Jang, H.J.; Hong, E.M.; Park, S.W.; Byun, H.W.; Koh, D.H.; Choi, M.H.; Kae, S.H.; Lee, J. Statin induces apoptosis of human colon cancer cells and downregulation of insulin-like growth factor 1 receptor via proapoptotic ERK activation. Oncol. Lett.; 2016; 12, pp. 250-256. [DOI: https://dx.doi.org/10.3892/ol.2016.4569]
48. Nordstrand, A.; Lundholm, M.; Larsson, A.; Lerner, U.H.; Widmark, A.; Wikström, P. Inhibition of the insulin-like growth factor-1 receptor enhances effects of simvastatin on prostate cancer cells in co-culture with bone. Cancer Microenviron.; 2013; 6, pp. 231-240. [DOI: https://dx.doi.org/10.1007/s12307-013-0129-z]
49. Bonifacio, A.; Sanvee, G.M.; Bouitbir, J.; Krähenbühl, S. The AKT/mTOR signaling pathway plays a key role in statin-induced myotoxicity. Biochim. Biophys. Acta; 2015; 1853, pp. 1841-1849. [DOI: https://dx.doi.org/10.1016/j.bbamcr.2015.04.010]
50. Zhang, H.; Chi, M.; Chen, L.; Sun, X.; Wan, L.; Yang, Q.; Guo, C. Linalool Prevents Cisplatin Induced Muscle Atrophy by Regulating IGF-1/Akt/FoxO Pathway. Front. Pharmacol.; 2020; 11, 598166. [DOI: https://dx.doi.org/10.3389/fphar.2020.598166]
51. Apostolopoulou, M.; Corsini, A.; Roden, M. The role of mitochondria in statin-induced myopathy. Eur. J. Clin. Investig.; 2015; 45, pp. 745-754. [DOI: https://dx.doi.org/10.1111/eci.12461] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25991405]
52. Auer, J.; Sinzinger, H.; Franklin, B.; Berent, R. Muscle- and skeletal-related side-effects of statins: Tip of the iceberg?. Eur. J. Prev. Cardiol.; 2016; 23, pp. 88-110. [DOI: https://dx.doi.org/10.1177/2047487314550804] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25230981]
53. Mollazadeh, H.; Tavana, E.; Fanni, G.; Bo, S.; Banach, M.; Pirro, M.; von Haehling, S.; Jamialahmadi, T.; Sahebkar, A. Effects of statins on mitochondrial pathways. J. Cachexia Sarcopenia Muscle; 2021; 12, pp. 237-251. [DOI: https://dx.doi.org/10.1002/jcsm.12654] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33511728]
54. Verbeek, R.; Hovingh, G.K.; Boekholdt, S.M. Non-high-density lipoprotein cholesterol: Current status as cardiovascular marker. Curr. Opin. Lipidol.; 2015; 26, pp. 502-510. [DOI: https://dx.doi.org/10.1097/MOL.0000000000000237]
55. SCORE2 working group and ESC Cardiovascular risk collaboration. SCORE2 risk prediction algorithms: New models to estimate 10-year risk of cardiovascular disease in Europe. Eur. Heart J.; 2021; 42, pp. 2439-2454. [DOI: https://dx.doi.org/10.1093/eurheartj/ehab309]
56. SCORE2-OP working group and ESC Cardiovascular risk collaboration. SCORE2-OP risk prediction algorithms: Estimating incident cardiovascular event risk in older persons in four geographical risk regions. Eur. Heart J.; 2021; 42, pp. 2455-2467. [DOI: https://dx.doi.org/10.1093/eurheartj/ehab312]
57. Ference, B.A.; Majeed, F.; Penumetcha, R.; Flack, J.M.; Brook, R.D. Effect of naturally random allocation to lower low-density lipoprotein cholesterol on the risk of coronary heart disease mediated by polymorphisms in NPC1L1, HMGCR, or both: A 2 × 2 factorial Mendelian randomization study. J. Am. Coll. Cardiol.; 2015; 65, pp. 1552-1561. [DOI: https://dx.doi.org/10.1016/j.jacc.2015.02.020]
58. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res.; 2004; 64, pp. 5245-5250. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-04-0496]
59. Belmonte, T.; Mangas, A.; Calderon-Dominguez, M.; Quezada-Feijoo, M.; Ramos, M.; Campuzano, O.; Gomez, S.; Peña, M.L.; Cubillos-Arango, A.M.; Dominguez, F. et al. Peripheral microRNA panels to guide the diagnosis of familial cardiomyopathy. Transl. Res.; 2020; 218, pp. 1-15. [DOI: https://dx.doi.org/10.1016/j.trsl.2020.01.003]
60. Mi, H.; Muruganujan, A.; Ebert, D.; Huang, X.; Thomas, P.D. PANTHER version 14: More genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res.; 2019; 47, pp. D419-D426. [DOI: https://dx.doi.org/10.1093/nar/gky1038]
61. Clum, A.; Huntemann, M.; Bushnell, B.; Foster, B.; Foster, B.; Roux, S.; Hajek, P.P.; Varghese, N.; Mukherjee, S.; Reddy, T. et al. DOE JGI Metagenome Workflow. Msystems; 2021; 6, pp. e00804-e00820. [DOI: https://dx.doi.org/10.1128/mSystems.00804-20] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34006627]
62. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P. et al. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res.; 2019; 47, pp. D607-D613. [DOI: https://dx.doi.org/10.1093/nar/gky1131] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30476243]
63. Liao, Y.; Wang, J.; Jaehnig, E.J.; Shi, Z.; Zhang, B. WebGestalt 2019: Gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res.; 2019; 47, pp. W199-W205. [DOI: https://dx.doi.org/10.1093/nar/gkz401] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31114916]
64. Motulsky, H.J.; Brown, R.E. Detecting outliers when fitting data with nonlinear regression—A new method based on robust nonlinear regression and the false discovery rate. BMC Bioinform.; 2006; 7, 123. [DOI: https://dx.doi.org/10.1186/1471-2105-7-123] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16526949]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Atherosclerotic cardiovascular diseases (ASCVD) are the leading cause of morbidity and mortality in Western societies. Statins are the first-choice therapy for dislipidemias and are considered the cornerstone of ASCVD. Statin-associated muscle symptoms are the main reason for dropout of this treatment. There is an urgent need to identify new biomarkers with discriminative precision for diagnosing intolerance to statins (SI) in patients. MicroRNAs (miRNAs) have emerged as evolutionarily conserved molecules that serve as reliable biomarkers and regulators of multiple cellular events in cardiovascular diseases. In the current study, we evaluated plasma miRNAs as potential biomarkers to discriminate between the SI vs. non-statin intolerant (NSI) population. It is a multicenter, prospective, case-control study. A total of 179 differentially expressed circulating miRNAs were screened in two cardiovascular risk patient cohorts (high and very high risk): (i) NSI (n = 10); (ii) SI (n = 10). Ten miRNAs were identified as being overexpressed in plasma and validated in the plasma of NSI (n = 45) and SI (n = 39). Let-7c-5p, let-7d-5p, let-7f-5p, miR-376a-3p and miR-376c-3p were overexpressed in the plasma of SI patients. The receiver operating characteristic curve analysis supported the discriminative potential of the diagnosis. We propose a three-miRNA predictive fingerprint (let-7f, miR-376a-3p and miR-376c-3p) and several clinical variables (non-HDLc and years of dyslipidemia) for SI discrimination; this model achieves sensitivity, specificity and area under the receiver operating characteristic curve (AUC) of 83.67%, 88.57 and 89.10, respectively. In clinical practice, this set of miRNAs combined with clinical variables may discriminate between SI vs. NSI subjects. This multiparametric model may arise as a potential diagnostic biomarker with clinical value.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details


1 Research Unit, Biomedical Research and Innovation Institute of Cadiz (INiBICA), Puerta del Mar University Hospital, 11009 Cádiz, Spain;
2 Cardiology Service, Hospital Josep Trueta, University of Girona, 17007 Girona, Spain;
3 Research Unit, Biomedical Research and Innovation Institute of Cadiz (INiBICA), Puerta del Mar University Hospital, 11009 Cádiz, Spain;
4 UCERV, UCAMI, Servicio de Medicina Interna, Hospital Virgen del Rocío, 41013 Seville, Spain;
5 Lipid and Atherosclerosis Unit, IMIBIC/Hospital Universitario Reina Sofía/Universidad de Córdoba, 14004 Córdoba, Spain;
6 Cardiology Service, Hospital Josep Trueta, University of Girona, 17007 Girona, Spain;
7 Lipid and Atherosclerotic Unit, Internal Medicine Department, Puerta del Mar University Hospital, 11009 Cadiz, Spain;