ABSTRACT
-
Background
- Diabetes is a chronic disease with several long-term complications. Several glucose-lowering drugs are used to treat type 2 diabetes mellitus (T2DM), e.g., glimepiride and liraglutide, in which both having different modes of action. Circulating microRNAs (miRNAs) are suggested as potential biomarkers that are associated with the disease development and the effects of the treatment. In the current study we evaluated the effect of glimepiride, liraglutide on the expression of the circulating miRNAs.
-
Methods
- The present study is a post hoc trial from a previously randomized control trial comparing liraglutide versus glimepiride both in combination with metformin in subjects with T2DM, and subclinical heart failure. miRNAs were determined in the subjects’ serum samples with next generation sequencing. Expression patterns of the circulating miRNAs were analyzed using bioinformatic univariate and multivariate analyses (clinical trial registration: NCT01425580).
-
Results
- Univariate analyses show that treatment with glimepiride altered expression of three miRNAs in patient serum, miR-206, miR-182-5p, and miR-766-3p. Both miR-182-5p and miR-766-3p were also picked up among the top contributing miRNAs with penalized regularised logistic regressions (Lasso). The highest-ranked miRNAs with respect to Lasso coefficients were miR-3960, miR-31-5p, miR-3613-3p, and miR-378a-3p. Liraglutide treatment did not significantly influence levels of circulating miRNAs.
-
Conclusion
- Present study indicates that glucose-lowering drugs differently affect the expression of circulating miRNAs in serum in individuals with T2DM. More studies are required to investigate possible mechanisms by which glimepiride is affecting the expression of circulating miRNAs.
-
Keywords: Circulating microRNA; Diabetes mellitus, type 2; Glimepiride; Liraglutide
GRAPHICAL ABSTRACT
INTRODUCTION
- Diabetes is a prevalent and costly chronic disease associated with several long-term complications. Worldwide, the prevalence of diabetes has increased dramatically over the past two decades and is expected to affect more than 500 million adults by 2030, with most having type 2 diabetes mellitus (T2DM) [1]. Myocardial infarction is the most common cause of death in individuals with diabetes [2,3]. Although factors that are known to reduce the risk of myocardial infarction [4,5] include the use of lipid-lowering and antihypertensive medications and better glycemic control [6-8] an excess risk of death still exists.
- There are several glucose-lowering agents available for the treatment of T2DM, including glimepiride, a second-generation sulfonylurea, e.g., glimepiride, improving the relative insulin secretory deficit found in T2DM, and glucagon-like peptide 1 (GLP-1) receptor antagonists, e.g., liraglutide, which acts via insulin release and reduction in glucagon secretion in the presence of elevated glucose concentrations.
- MicroRNAs (miRNAs) consisting of short (21 to 22 nucleotides), single-stranded, non-coding RNAs are prevalent in human plasma and other biological fluids, and their expression has been proven to be disease-specific [9-13]. miRNAs have been identified as key regulators of complex biological processes linked to multiple cardiovascular pathologies, including left ventricular hypertrophy, ischemic heart disease, heart failure, hypertension, and arrhythmias. Compared with conventional biomarkers, miRNAs offer advantages with high stability [14]. In addition, miRNAs have been identified as regulators of the β-cells’ mass and function, and thus putatively involved in the pathogenesis of T2DM [13,15,16]. The aim of the present study is to evaluate the effect of two different classes of glucoselowering drugs, liraglutide and glimepiride, on the expression pattern of circulating miRNAs in patients with T2DM.
METHODS
- Cohort
- The present study is a post hoc trial from the main study [17-19]; which was an assessor-blinded and active-controlled parallel-group trial comparing liraglutide 1.8 mg subcutaneously versus glimepiride 4 mg orally both once daily and in combination with metformin during 18 weeks in subjects with T2DM and subclinical heart failure.
- The trial protocol was approved (number D-nr 2011/233-31/1) by regional ethics committees at both participating centers and was conducted in accordance with International Conference on Harmonization-Good Clinical Practice guidelines. Participants gave written informed consent before being enrolled. The clinical trial registration information can be found at ClinicalTrials.gov with the identifier NCT01425580.
- The definition of subclinical heart failure was fulfillment of one of the following criteria: left ventricular ejection fraction ≤50%, or evidence of systolic dysfunction as shown by reduced relaxation velocity, in the left ventricular long axis, of 20% compared to a normal population. Key exclusion criteria were heart failure classified according to the New York Heart Association (NYHA) classification III–IV; history of atrial fibrillation or flutter; presence of acute myocarditis or significant valvulopathies [19]. The primary endpoint was treatment change in either longitudinal functional reserve index (LFRI) and (diastolic functional reserve index [DFRI] and, systolic functional reserve index [SFRI], respectively). Briefly, LFRI was constructed by collecting e´ and s´ at rest (e´rest, s´rest), and during bicycle exercise stress test (e´exercise, s´exercise) using the formula: LFRIdiastolic=Δe´×[1−(1/e´rest)], and LFRIsystolic=Δs´×[1−(1/s´rest)], respectively. Where Δe´ is the difference between e´exercise and e´rest and Δs´ is the difference between s´exercise and s´rest (Note that LFRI has no magnitude).
- The baseline characteristics of the cohort are shown in Table 1.
- Serum samples collection
- Plasma samples were collected at baseline and after 18 weeks of treatment with liraglutide or glimepiride and the plasma samples were after spinning and aliquoting into the Eppendorf vessels and stored at –80°C.
- miRNA isolation and quality control
- Serum-ethylenediaminetetraacetic acid (EDTA) samples were preserved at –80°C until dry ice shipment to the Qiagen Genomic Services lab (Qiagen, Hilden, Germany) for further miRNA purification and analysis. According to Qiagen miRNA was purified from the serum-EDTA samples using miRNeasy Kit (Qiagen) according to the manufacturer’s instructions. Purified miRNA quality was checked by quantitative polymerase chain reaction (qPCR) to evaluate if the expression levels of selected miRNAs in our samples were within the expected range. In this analysis, the expression levels of miRNAs hsa-miR-103a-3p, hsa-miR-191-5p, hsa-miR-451a, hsa-miR23a-3p, and hsa-miR-30c-5p was done since those miRNAs are stably expressed in biofluids like serum and plasma. To assess hemolysis qPCR quantification of two miRNAs was used: miRNA-451 which is expressed in red blood cells and miRNA23a which is relatively stable in serum and plasma and not affected by hemolysis (miRNA-23a).
- Library preparation and sequencing
- The library preparation was done using the QIAseq miRNA Library Kit (Qiagen). A total of 5 μL total RNA was converted into miRNA next generation sequencing libraries. Adapters containing unique molecular identifiers were ligated to the RNA which was subsequently converted to cDNA. The cDNA was amplified using PCR and purified. Library preparation quality check was performed using Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Based on the quality of the inserts and the concentration measurements the libraries were pooled in equimolar ratios. The library pool(s) were quantified using qPCR and sequenced on a NextSeq550 sequencing instrument according to the manufacturer’s instructions. Raw data were demultiplexed and FASTQ files for each sample were generated using the bcl2fastq software (Illumina, San Diego, CA, USA). FASTQ data were checked using the FastQC tool (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Sequencing and processed miRNA data are available under gene accession no. GSE223538 at the public functional genomics data repository Gene Expression Omnibus (GEO).
- The raw sequencing data and its corresponding meta data for the miRNA results, are in the process of being uploaded for evaluation by the public functional genomics data repository GEO with minimum information about a microarray experiment (MIAME)-compliant data submission.
- Statistical and bioinformatical analysis
- For the statistical analysis of the cohort, continuous data are summarized as mean±standard deviation and categorical data are presented as percentages. A two-sided P value of less than 0.05 was considered statistically significant. Data were analyzed using Stata 14.2 (StataCorp., College Station, TX, USA).
- Qiagen Genomic Service generated the raw counts of miRNAs. miRNAs with a mean total read count ≥50 were used as input into the DESeq2 Bioconductor package in R (R Foundation for Statistical Computing, Vienna, Austria) for statistical analysis. Differential expression tests were conducted with DESeq function to estimate size factors, dispersions and for fitting a patient paired binomial generalized linear models for each miRNA and Wald test for significance testing (design: patient+gender+cardiovascular disease background+treatment). miRNAs were considered as differentially expressed if the absolute value of fold change was >1.5 and the Benjamini-Hochberg adjusted P value was <0.05.
- Linear mixed effects models evaluating effect of differentially expressed miRNA and treatment changes in SFRI and DFRI were fitted for fixed effects of gender, cardiovascular disease background and treatment, with individual patients as random effect (lme4 and lmerTest package in R) [19].
- Unsupervised dimensionality reduction of miRNA expressions was performed with principal component analysis (PCA) and with Uniform Manifold Approximation and Projection (UMAP; embed package in R) algorithms in R. Unsupervised Gaussian finite mixture model clustering and evaluation were performed on miRNA expressions (mclust package in R). Supervised UMAP and partial least squares (PLS) in R were used for supervised dimensionality reduction of miRNA expressions for dimensional separation of treatment groups.
- Evaluation of multivariate effects of miRNA expressions related to treatments was performed with a nested cross-validation machine learning workflow assessing miRNAs as predictors and classifiers of treatments. miRNAs with a total variance <0.2 and missing values >30% were excluded as predictors. Data were partitioned into training (75%) and testing (25%) sets. Each set was z-normalized, and knn-imputation of missing values was performed on the training set. Thereafter, penalized regularised logistic regressions (Lasso; glmnet package in R), Random forests (ranger package in R), XGBoosted trees (xgboost package in R), and neural network, performed on both miRNAs and PCA-derived principal components as main variables (TensorFlow and Keras integration with R) were in parallel hyperparameter tuned. Hyperparameters for Lasso (penalty); Random Forest (number of variables included in each random tree and minimum n for split); XGBoosted trees (number of variables included in each tree, tree depth, loss reduction, learning rate, and minimum n for split) and neural network (number of hidden layers and penalty) were tuned with a Latin hypercube search approach. Tuning was performed with internal validation on 25 bootstraps of the training data with evaluation through the receiver operating characteristic (ROC) area under the curve (AUC). Optimal hyperparameters were used in the final models and performance was tested on hold-out testing data. This workflow procedure was repeated iteratively 100 times as randomized nested cross-validation with random partitioning of samples to training and testing set at each iteration to estimate the robustness of predictions with regards to random effects of partitioning. Model variable importance was assessed with coefficients (Lasso) and permutation (Random forest, XGBoosted trees) with the variable importance plots (VIP) package in R. Model performance was compared using random effects and Bayesian statistics with Stan [20] through specifying and fitting the models via the rstanarm package in R for posterior distributions of mean AUC and difference in AUC.
- The CyTargetLinker plug-in tool (version 4.1.0) in Cytoscape (version 3.9.1) was used to predict miRNA target genes (linkset miRTarBase release 8.0). The stringApp (version 1.7) plug-in tool in Cytoscape was used to retrieve and filter for Gene Ontology (GO) annotations for the targeted genes.
RESULTS
- One hundred and five patients were screened for eligibility, 62 met the criteria and were randomized to liraglutide or glimepiride. Basal characteristics of the study population are shown in Table 1. Treatment groups were well-balanced except for triacylglycerols and glycosylated hemoglobin which was higher in subjects randomized to liraglutide treatment. Treatment change of secondary outcomes after 18-week treatment showed no significant changes in clinical or echocardiographic parameters apart from body weight (–3.7 kg vs. –0.2 kg), waist circumference (–3.1 cm vs. –0.8 cm), and heart rate (6.3 bpm vs. –2.3 bpm) in liraglutide and glimepiride respectively, as shown in Table 2.
- In total 2,632 different miRNAs were detected among the patient’s serum samples. Of these 683 miRNAs were expressed in at least 25% of the samples and 320 in at least 90% of the samples. Two hundred and three miRNAs with a total mean read count ≥50 were further evaluated for univariate differential expression and multivariate patterns in machine learning associated with treatment with liraglutide or glimepiride in individuals with T2DM. Additionally, we evaluated if the manifestation of cardiovascular disease or their gender has an impact on the miRNA expression.
- Unsupervised discrimination of global variation of miRNA expression in patient serum before and after treatment
- Unsupervised dimensional reduction with PCA and UMAP shows no apparent capture of global variation in principal components one to four associated with treatment groups (Fig. 1A and B). This was also true for gender and cardiovascular disease background (Supplementary Fig. 1). In addition, unsupervised Gaussian finite mixture model-based clustering shows no apparent clusters related to the treatment groups (Fig. 1C), an indication that treatment did not induce a large variation in miRNA expressions in blood.
- Differential univariate expression of miRNA in patient serum after treatments
- Univariate analysis with patient paired design showed that glimepiride treatment significantly altered the expression of three circulating serum miRNAs: miR-206, miR-182-5p, and miR-766-3p (Fig. 2, Supplementary Table 1). These changes in miRNA expression were not dependent on either the sex of the patients or their cardiovascular disease background (data not shown). Treatment with liraglutide didn’t cause any significant changes in miRNAs in serum (Supplementary Table 2).
- Linear mixed effect models show an increased performance with miR-206 and miR-182-5p on DFRI when controlling for gender, cardiovascular disease background and treatment. When added miR-206 and miR-182-5p significantly increased model R2 from 0.73 to 0.75 (P=0.028) and from 0.73 to 0.74 (P=0.019) respectively. However, only miR-182-5p, and not miR-206, show a significant model effect size (0.629, P=0.036; and 0.556, P=0.053, respectively). Analysis on SFRI shows no significant changes correlated to the miR-206, miR-182-5p, and miR-766-3p (data not shown).
- Supervised discrimination of global variation of miRNA expression in patient serum before and after treatment
- For evaluation of possible multivariate effects, supervised dimensional reduction with PLS and UMAP was performed and show partial separation between some of the treatment groups, especially with glimepiride treatment (Supplementary Fig. 2). This multivariate separation was further evaluated with supervised machine learning to find signatures for before and after treatment.
- Multivariate expression patterns of miRNA in serum for discrimination of before and after treatment
- The different machine learning algorithms used to model multivariate patterns between before and after treatment show a moderate predictive power of miRNAs to separate samples before and after glimepiride treatment. This was not observed for the liraglutide treatment (data not shown). Evaluation of model performances with posterior distributions of AUC (Fig. 3A) showed that the Lasso models had the best performance with an average iteration mean AUC of 0.68 (95% cofidence interval, 0.61 to 0.75) (Fig. 3B). The variable importance for the Lasso models was evaluated by the averaged iterated variable coefficients of all included miRNAs. The top 40 contributing miRNAs are presented in Fig. 3C. miR-182-5p and miR-766-3p, but not miR-206, were among the top 40 contributing miRNAs. However, they were not among the highest-ranked miRNAs. The two miRNAs with the highest positive coefficients were miR-3960 and miR-31-5p, and the highest negative coefficients were miR-3613-3p and miR-378a-3p. These were in the univariate analyses either not statistically significant or did not reach the cut-off for fold change, reflecting their dependence of multivariate co-expression of other miRNAs for combined predictive strength.
- Functional annotation and target interaction of differentially expressed miRNA
- The miRNAs miR-206, miR-182-5p, and miR-766-3p that are seen to be differentially expressed with glimepiride have many biological targets that can be predicted with miRTarBase target analysis. Their involvement in processes and pathways associated with insulin signalling and signalling pathways in diabetes was given focus with target predictions of the miRNAs towards targets belonging to gene ontologies and pathways of relevance. miRTarBase predicted targets connected to insulin signalling/responses are shown in Fig. 4A and connected to advanced glycation endproduct/receptor for advanced glycation end product (AGE-RAGE) signalling pathway in diabetic complications are shown in Fig. 4B.
DISCUSSION
- In this post hoc randomized control trial, we evaluated the effect of two different classes of glucose-lowering drugs, liraglutide and glimepiride, on the expression pattern of circulating miRNAs in people with T2DM before and after 3-month of treatment. The univariate results of our study show that glimepiride significantly changes serum levels of three miRNAs: miR-206, miR-182-5p, and miR-766-3p. Although a combined multivariate expression of miRNAs reaches a moderate level of machine learning predictive power, neither treatment appears to extensively affect the expression of circulating miRNAs in the present study.
- miR-206 belongs to the myomiR family and is probably one of the most studied miRNAs today. Profiling of miRNA expression in different tissues showed that miR-206 was found to be primarily expressed in the skeletal muscle but also in the heart [21]. In skeletal muscle, this miRNA positively regulates muscle growth and differentiation and showed to be specifically induced during myogenesis [22]. In cardiomyocytes, an increase of the miR-206 expression was shown to negatively regulate apoptosis after ischemic injury [23]. It was also shown that an intracellular increase in miR-1/miR206 in vitro was able to promote high glucose-induced cardiac cell apoptosis [24].
- The activity of miR-206 has also been shown in other tissues, such as pancreatic islets, brain, intestine, and liver where miR-206 in mice was strongly induced by a high-fat diet, and involved in the repression of the activity of the rate-limiting glucose sensor–glucokinase [25]. This study showed that miR-206 was involved in the regulation of an important glycolytic gene glucokinase, in the murine islets, thus indicating involvement in glucose homeostasis. However, this relationship between miR-206 and glucokinase has not been confirmed in humans [26].
- miRNAs are acting intracellularly but their presence has been shown in different biological fluids including plasma and serum and is suggested as biomarkers for different physiological and pathological conditions [27-29]. Levels of circulating miR-206 were shown to be affected by several conditions such as muscle damage induced by physical exercise [30], muscular dystrophy [31], or amyotrophic lateral sclerosis [32] as well as cancer [33,34]. miR-206 has been suggested as a biomarker of disease development, especially in cancer patients. Also, cardiovascular diseases are shown to affect levels of miR-206. Decreased levels of circulating miR-206 have been associated with increased blood pressure in patients with left heart disease, LHD [35]. miR-206 together with let-7c-5p, miR-765, miR483-5p, and miR-31-5p were shown to be elevated in the plasma of patients with coronary heart disease [36]. Albeit, in several studies reporting the changes of miRNA levels in the serum/plasma levels in individuals with diabetes (reviewed in [37]), the levels of circulating miR-206 were not shown to be affected by this condition. Using both in vivo (mice) and in vitro (human hepatocytes) studies, Wu et al. [38] have shown miR-206 to be an effective inhibitor for gluconeogenesis and lipogenesis in studied models thus having therapeutic potential in the treatment of hyperglycemia and nonalcoholic fatty liver disease. The effect of mir-206 was shown to be mediated by protein tyrosine phosphatase non-receptor type 1-insulin receptor (PTPN1-INSR)/insulin receptor substrate 1 (IRS1) axis.
- Studies have also established links between miR-182-5p and miR-766-3p and cardiovascular diseases. Elevated levels of miR-182-5p together with miR-199a-5p and miR-5187-5p were found to be associated with the presence of the unprotected left main coronary artery disease (uLMCAD) in patients with symptoms of coronary artery disease [39]. In a study aimed to investigate the diagnostic potential of miR-182-5p together with miR30a-5p, those miRNAs were shown to be associated with dysglycemia. Expression levels of miR-182-p were found to increase in individuals with prediabetes or screen-detected diabetes individuals; however, in individuals with T2DM expression of that miRNA was decreasing. This suggests that levels of that miRNA were decreasing with the duration of the disease [40,41]. This is in line with conclusions made by others [42] that also showed the miR-182 to be a potential modulator of the transcription factor Forkhead box protein O1 (FOXO1), known to regulate hepatic gluconeogenesis and glycogenolysis in response to insulin. miR-182 is a negative regulator of the FOXO1 protein and was suggested to reduce gluconeogenesis in a compensatory way during prediabetes but increase during T2DM.
- miR-766-3p has previously been shown to be involved in cancer, be connected to aging, and have anti-inflammatory properties by reducing nuclear factor-κB activity [43-45]. T2DM individuals were having a lower expression of miR-766-3p compared to healthy subjects [46], but levels of that miRNA were positively affected by physical training [47]. In the study by Cao et al. [48], it was shown that miR-766-3p has its putative binding site in the sequence of long non-coding RNA LINC-P21 and that the levels of LINC-P21 and miR-766-3p expression were correlated in a reverse ratio in a pancreatic β-cell line INS-1. The study also concludes that LINC-P21 inhibits the proliferation of β-cells and decreases insulin secretion by competitively binding miR-766-3p and thus up-regulating another target of miR-766-3p—a nuclear receptor nuclear receptor subfamily 3 group C member 2 (NR3C2) gene.
- Findings on drug-induced changes in plasma/serum miRNA levels are scarce. To our knowledge, only a few studies have evaluated the effect of glucose-lowering drugs on circulating miRNAs. Demirsoy et al. [49] in a targeted miRNA study shows that in subjects with T2DM, glucose-lowering treatment with metformin downregulated 13 different miRNAs in plasma samples (hsa-let-7e-5p, hsa-let-7f-5p, hsa-miR-21-5p, hsamiR-24-3p, hsa-miR-26b-5p, hsa-miR-126-5p, hsa-miR-129-5p, hsa-miR-130b-3p, hsa-miR-146a-5p, has-miR-148a-3p, hsa-miR-152-3p, hsa-miR-194-5p, hsa-miR-99a-5p). In addition, metformin has also been shown to significantly change the expression of miR-192, miR-140-5p, and miR-222 [50]. The results of the two studies are not confirmative and the analysis depth of miRNAs differ significantly between them. In the present study glimepiride and liraglutide was an add-on to metformin treatment, which may have demasked the effects of the metformin on circulating miRNAs. Also, in the present study, with a more in-depth characterisation of circulating miRNAs, none of the above miRNAs were significantly affected by glimepiride or liraglutide, in combination with metformin. However, differences in treatment’s impact on the expression of circulated miRNA might depend on the differences in mechanism action between the different glucose-lowering drugs. Metformin decreases blood glucose through inhibition of mitochondrial complex I activity and decrease of adenosine triphosphate (ATP) production in hepatocytes, leading to activation of adenosine monophosphate (AMP) activated protein kinase-activated protein kinase and inhibition of fructose-1,6-bisphosphatase enzyme. Glimepiride blocks ATP-sensitive potassium channels (KATP channels) by binding to sulfonylurea receptors 1 and 2 on pancreatic β-cells. Blockage of (KATP channels) causes depolarisation of the cell membrane and release of insulin granules. In contrast, liraglutide increases intracellular production of cyclic AMP which in turn stimulates the glucose-dependent release of insulin and glucose-dependent inhibition of glucagon. Despite knowledge about the mechanisms of action of those drugs it is still difficult to explain why glimepiride changes the expression of the discussed above miRNA. Interestingly, expression levels of miR-186-5p were previously shown to be insensitive to the treatment with metformin and/or insulin [41].
- The current study was conducted on a relatively small number of patients which, with respect to high interindividual variations in circulating miRNAs, might introduce some uncertainty in the results with regards to subtle changes after treatment. In addition, a longer treatment period with glucose-lowering drugs may also result in more pronounced differential expressions of the circulating serum miRNA. Nevertheless, the in-depth characterisation of circulating miRNAs and the interesting results obtained in the current study prompts the necessity of further studies to evaluate mechanisms that involve miRNAs in patients’ response to glimepiride and liraglutide.
- In conclusion, the present study indicates that glucose-lowering drugs differently affect the expression of circulating miRNAs in serum in individuals with T2DM. More studies are required to investigate possible mechanisms by which glimepiride is affecting the expression of circulating miRNAs. Since the miRNAs, studied in the present study, are suggested as potential biomarkers in patients with different diseases, further studies on the effect of glimepiride on plasma levels of those biomarkers are much needed.
SUPPLEMENTARY MATERIALS
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2022.0342.
Supplementary Fig. 1.
Distribution of variation in microRNA (miRNA) expression in samples with dimensionality reduction of
miRNA expressions with respect to (A) gender and (B) cardiovascular disease background. Figures are overlaid with colour representing the groups (male or female and with cardiovascular disease background [CVD] or without [NO_CVD]).
dmj-2022-0342-Supplementary-Fig-1.pdf
Supplementary Fig. 2.
Distribution of variation in microRNA (miRNA) expression in samples with respect to treatments. Supervised dimensionality reduction of miRNA expressions with (A) partial least squares (PLS) and with (B) supervised uniform manifold approximation and projection (UMAP), with first four dimensions shown. Figures are overlaid with colour representing the treatment groups (before treatment with liraglutide [L0] or glimepiride [G0] and after treatment with liraglutide [L1] or glimepiride [G1]).
dmj-2022-0342-Supplementary-Fig-2.pdf
NOTES
-
CONFLICTS OF INTEREST
Nikolai N. Scherbak and Robert Kruse has nothing to declare. Johan Jendle has received grants from Medtronic and Novo Nordisk; and lecture/other fees from Abbott, AstraZeneca, Boehringer Ingelheim, Eli Lilly, Nordic Infucare, Novo Nordisk, and Sanofi. Thomas Nyström as received unrestricted grants from AstraZeneca and NovoNordisk and has served on national advisory boards of Amgen, Novo Nordisk, Sanofi-Aventis, Eli Lilly, MSD, and Boehringer Ingelheim.
-
AUTHOR CONTRIBUTIONS
Conception or design: N.S., T.N., J.J.
Acquisition, analysis, or interpretation of data: N.S., R.K., T.N., J.J.
Drafting the work or revising: N.S., R.K., J.J.
Final approval of the manuscript: N.S., R.K., T.N., J.J.
-
FUNDING
This study was investigator-initiated and investigator-designed. The investigators received an unrestricted grant from Novo Nordisk A/S, but the company was not involved in data collection, study management, analysis, or interpretation of data. Nor was the company involved in the decisions regarding the submission of the manuscript.
Acknowledgements- None
Fig. 1.Distribution of variation in microRNA (miRNA) expression in samples with respect to treatments. Unsupervised dimensionality reduction of miRNA expressions with (A) principal component analysis (PCA) and with (B) uniform manifold approximation and projection (UMAP), with first four dimensions shown, and (C) Gaussian finite mixture model-based clustering of miRNA expressions into four clusters. Figures are overlaid with colour representing the treatment groups (before treatment with liraglutide [L0] or glimepiride [G0] and after treatment with liraglutide [L1] or glimepiride [G1]).
Fig. 2.Univariate differential expression in microRNAs (miRNAs) before and after glimepiride treatment. (A) Volcano plot of miRNA log fold changes (FCs) and statistical significance show three differentially expressed miRNAs (red dots) further depicted in (B, C, D) box-violin plots with patient paired individual points overlaid. Statistical analysis performed with DESeq fitting patient paired binomial generalized linear models and Wald test for significance testing. Biological and statistical significance was set at a FC >1.5 and a Benjamini-Hochberg adjusted P<0.05, respectively. NS, not significant.
Fig. 3.Machine learning multivariate classification of microRNA (miRNA) expression before and after glimepiride treatment. (A) Posterior distributions of receiver operating characteristic (ROC) and area under the curve (AUC) from penalized regularised logistic regressions (Lasso), Random forests, XGBoosted trees, and Neural network (performed on both miRNAs and principal component analysis [PCA]-derived principal components). (B) Show the ROC-curve and (C) show the averaged iterated variable coefficients of the top 40 contributing miRNAs for the Lasso models. LR, penalized regularised logistic regressions (Lasso); NN, Neural network; NN_PCA, Neural network on PCA-derived principal components; RF, random forest; XG, XGBoosted trees; ST, stack of LR, XG and RF.
Fig. 4.The three microRNAs (miRNAs) differentially expressed with glimepiride treatment; miR-206, miR-182-5p, and miR-766-3p, and their mRNA-associated targets with a focus on targets belonging to (A) the gene ontologies/WikiPathways “Cellular response to insulin stimulus” (GO:0032869), “Response to insulin” (GO:0032868), “Insulin signalling” (WP481), and “Insulin receptor signalling pathway” (GO:0008286), and (B) targets belonging to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway “advanced glycation endproduct/receptor for advanced glycation end product (AGE-RAGE) signalling pathway in diabetic complications” (hsa04933). Node color of miRNAs reflect their fold change and edges reflects miRTarBase prediction towards gene products (grey nodes).
Table 1.Baseline characteristics of the study population
Characteristic |
Liraglutide (n=33) |
Glimepiride (n=29) |
P value |
Age, yr |
60.8±7.6 |
63.0±6.8 |
0.240a
|
Male sex |
24 (72.7) |
21 (72.4) |
1.000b
|
Diabetes duration, yr |
5 (1–10) |
1 (3–7) |
0.368c
|
Smoking |
3 (9.1) |
4 (13.8) |
0.852b
|
BMI, kg/m2
|
30.5±4.4 |
29.0±3.2 |
0.152a
|
Body weight, kg |
91.8±15.9 |
89.0±9.9 |
0.411a
|
Waist circumference, cm |
109.0±13.0 |
106.3±9.7 |
0.366a
|
Mean systolic BP, mm Hg |
131.9±14.0 |
129.3±10.9 |
0.414a
|
Mean diastolic BP, mmHg |
76.7±7.9 |
77.2±7.9 |
0.838a
|
eGFR, mL/min/1.72 m2
|
88.3±15.0 |
87.4±13.1 |
0.799a
|
Complications |
|
|
|
|
Hypertension |
29 (87.9) |
21 (72.4) |
0.224b
|
|
Hyperlipidemia |
25 (75.8) |
23 (79.3) |
0.980b
|
|
Coronary artery disease |
10 (30.3) |
11 (37.9) |
0.714b
|
|
Stroke |
1 (3) |
2 (6.9) |
0.902b
|
|
Proliferative retinopathy |
1 (3) |
1 (3.4) |
1.000b
|
Treatment |
|
|
|
|
Antiplatelet therapy |
11 (33.3) |
12 (41.4) |
0.696b
|
|
Anticoagulant treatment |
3 (9.1) |
1 (3.4) |
0.714b
|
|
ACE inhibitors/ARB blockers |
25 (75.8) |
20 (69) |
0.754b
|
|
Beta-blockers |
14 (42.4) |
13 (44.8) |
1.000b
|
|
Calcium inhibitors |
13 (39.4) |
10 (34.5) |
0.894b
|
|
Diuretics |
11 (33.3) |
6 (20.7) |
0.408b
|
|
Statins |
22 (66.7) |
24 (82.8) |
0.248b
|
Biochemical parameters |
|
|
|
|
HbA1c, mmol/mol |
54 (50–60) |
50 (49–54) |
0.036c
|
|
Triglycerides, mmol/L |
2.0 (1.4–2.6) |
1.5 (1.0–2.2) |
0.029c
|
|
Total cholesterol, mmol/L |
4.4 (4.0–6.0) |
4.5 (3.7–4.8) |
0.370c
|
|
LDL-cholesterol, mmol/L |
2.8±1.2 |
2.5±1.0 |
0.440a
|
Echocardiographic parameters |
|
|
|
|
LVEF at rest, % |
53.6±11.5 |
53.7±9.2 |
0.994 |
|
GLS at rest, % |
−15.3±4.2 |
−16.5±3.6 |
0.319 |
Table 2.Treatment change of secondary outcomes during 18-week treatment
Clinical characteristic |
Liraglutide (n=31) |
Glimepiride (n=28) |
95% CI |
P value |
Body weight, kg |
–3.7±3.5 |
–0.2±4.3 |
–5.5 to –1.4 |
0.001 |
Waist circumference, cm |
–3.1±2.8 |
–0.8±4.4 |
–4.2 to –0.4 |
0.019 |
Mean systolic BP, mm Hg |
–6.6±12.7 |
–1.9±12.0 |
–11.2 to 1.7 |
0.146 |
Mean diastolic BP, mm Hg |
–0.3±8.0 |
–0.9±7.9 |
–3.5 to 4.8 |
0.761 |
HR at rest, bpm |
6.3±9.6 |
–2.3±9.4 |
3.0 to 14.2 |
0.003 |
HR during exercise, bpm |
1.3±16.6 |
–4.9±11.3 |
–2.3 to 14.7 |
0.149 |
HbA1c, mmol/mol |
–11.2±7.7 |
–9.0±10.6 |
–7.0 to 2.6 |
0.366 |
Triglyceride, mmol/L |
–0.2±0.4 |
–0.1±0.8 |
–0.4 to 0.2 |
0.505 |
LDL, mmol/L |
–0.2±0.4 |
–0.3±0.9 |
–0.3 to 0.4 |
0.644 |
HDL, mmol/L |
–0.1±0.2 |
0.0±0.2 |
0.0 to –0.1 |
0.471 |
CRP, mmol/L |
0.5±6.9 |
1.45±8.9 |
–5.1 to 3.2 |
0.650 |
NT-proBNP, mg/L |
–53.6±505.4 |
6.6±57.1 |
–252.8 to 132.3 |
0.533 |
Longitudinal functional reserve index |
26 |
22 |
|
|
DFRI |
−0.18 |
–0.53 |
–0.28 to 2.59 |
0.190 |
SFRI |
–0.10 |
–0.18 |
–1.0 to 1.7 |
0.540 |
REFERENCES
- 1. NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in diabetes since 1980: a pooled analysis of 751 population-based studies with 4.4 million participants. Lancet 2016;387:1513-30.PubMedPMC
- 2. Holman RR, Paul SK, Bethel MA, Matthews DR, Neil HA. 10-Year follow-up of intensive glucose control in type 2 diabetes. N Engl J Med 2008;359:1577-89.ArticlePubMed
- 3. Rao Kondapally Seshasai S, Kaptoge S, Thompson A, Di Angelantonio E, Gao P, Sarwar N, et al. Diabetes mellitus, fasting glucose, and risk of cause-specific death. N Engl J Med 2011;364:829-41.ArticlePubMedPMC
- 4. UK Prospective Diabetes Study Group. Tight blood pressure control and risk of macrovascular and microvascular complications in type 2 diabetes: UKPDS 38. BMJ 1998;317:703-13.ArticlePubMedPMC
- 5. Colhoun HM, Betteridge DJ, Durrington PN, Hitman GA, Neil HA, Livingstone SJ, et al. Primary prevention of cardiovascular disease with atorvastatin in type 2 diabetes in the Collaborative Atorvastatin Diabetes Study (CARDS): multicentre randomised placebo-controlled trial. Lancet 2004;364:685-96.ArticlePubMed
- 6. Grant RW, Pirraglia PA, Meigs JB, Singer DE. Trends in complexity of diabetes care in the United States from 1991 to 2000. Arch Intern Med 2004;164:1134-9.ArticlePubMed
- 7. American Diabetes Association. Standards of medical care in diabetes: 2014. Diabetes Care 2014;37 Suppl 1:S14-80.ArticlePubMedPDF
- 8. Swedish National Diabetes Register (NDR) Centre of Registers Västra Götaland. Nation-wide results 1996-2020. Gothenburg: NDR; 2020 Available from: https://www.ndr.nu/pdfs/NationWideResults_1996-2020.pdf.
- 9. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A 2008;105:10513-8.ArticlePubMedPMC
- 10. Zampetaki A, Kiechl S, Drozdov I, Willeit P, Mayr U, Prokopi M, et al. Plasma microRNA profiling reveals loss of endothelial miR-126 and other microRNAs in type 2 diabetes. Circ Res 2010;107:810-7.ArticlePubMed
- 11. De Smaele E, Ferretti E, Gulino A. MicroRNAs as biomarkers for CNS cancer and other disorders. Brain Res 2010;1338:100-11.ArticlePubMed
- 12. Ciccacci C, Morganti R, Di Fusco D, D’Amato C, Cacciotti L, Greco C, et al. Common polymorphisms in MIR146a, MIR128a and MIR27a genes contribute to neuropathy susceptibility in type 2 diabetes. Acta Diabetol 2014;51:663-71.ArticlePubMedPDF
- 13. Barbagallo D, Piro S, Condorelli AG, Mascali LG, Urbano F, Parrinello N, et al. miR-296-3p, miR-298-5p and their downstream networks are causally involved in the higher resistance of mammalian pancreatic α cells to cytokine-induced apoptosis as compared to β cells. BMC Genomics 2013;14:62.ArticlePubMedPMCPDF
- 14. Romaine SP, Tomaszewski M, Condorelli G, Samani NJ. MicroRNAs in cardiovascular disease: an introduction for clinicians. Heart 2015;101:921-8.ArticlePubMedPMC
- 15. Guay C, Regazzi R. New emerging tasks for microRNAs in the control of β-cell activities. Biochim Biophys Acta 2016;1861(12 Pt B):2121-9.ArticlePubMed
- 16. Sebastiani G, Po A, Miele E, Ventriglia G, Ceccarelli E, Bugliani M, et al. MicroRNA-124a is hyperexpressed in type 2 diabetic human pancreatic islets and negatively regulates insulin secretion. Acta Diabetol 2015;52:523-30.ArticlePubMedPDF
- 17. Jendle J, Fang X, Cao Y, Bojo L, Nilsson BK, Hedberg F, et al. Effects on repetitive 24-hour ambulatory blood pressure in subjects with type II diabetes randomized to liraglutide or glimepiride treatment both in combination with metformin: a randomized open parallel-group study. J Am Soc Hypertens 2018;12:346-55.ArticlePubMed
- 18. Nystrom T, Santos-Pardo I, Fang X, Cao Y, Hedberg F, Jendle J. Heart rate variability in type 2 diabetic subjects randomized to liraglutide or glimepiride treatment, both in combination with metformin: a randomized, open, parallel-group study. Endocrinol Diabetes Metab 2019;2:e00058.PubMedPMCPDF
- 19. Nystrom T, Santos-Pardo I, Hedberg F, Wardell J, Witt N, Cao Y, et al. Effects on subclinical heart failure in type 2 diabetic subjects on liraglutide treatment vs. glimepiride both in combination with metformin: a randomized open parallel-group study. Front Endocrinol (Lausanne) 2017;8:325.PubMedPMC
- 20. McElreath R. Statistical rethinking: a Bayesian course with examples in R and Stan. 2nd ed. London: Chapman and Hall/CRC; 2020. p. 594.
- 21. Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V. Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol 2004;5:R13.ArticlePubMedPMC
- 22. Rao PK, Kumar RM, Farkhondeh M, Baskerville S, Lodish HF. Myogenic factors that regulate expression of muscle-specific microRNAs. Proc Natl Acad Sci U S A 2006;103:8721-6.ArticlePubMedPMC
- 23. Yan Y, Dang H, Zhang X, Wang X, Liu X. The protective role of MiR-206 in regulating cardiomyocytes apoptosis induced by ischemic injury by targeting PTP1B. Biosci Rep 2020;40:BSR20191000.ArticlePubMedPMCPDF
- 24. Shan ZX, Lin QX, Deng CY, Zhu JN, Mai LP, Liu JL, et al. miR-1/miR-206 regulate Hsp60 expression contributing to glucose-mediated apoptosis in cardiomyocytes. FEBS Lett 2010;584:3592-600.ArticlePubMedPDF
- 25. Vinod M, Patankar JV, Sachdev V, Frank S, Graier WF, Kratky D, et al. MiR-206 is expressed in pancreatic islets and regulates glucokinase activity. Am J Physiol Endocrinol Metab 2016;311:E175-85.ArticlePubMedPMC
- 26. Joglekar MV, Wong WK, Maynard CL, Umrani MR, Martin D, Loudovaris T, et al. Expression of miR-206 in human islets and its role in glucokinase regulation. Am J Physiol Endocrinol Metab 2018;315:E634-7.ArticlePubMed
- 27. Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, et al. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res 2008;18:997-1006.ArticlePubMedPDF
- 28. Raffort J, Hinault C, Dumortier O, Van Obberghen E. Circulating microRNAs and diabetes: potential applications in medical practice. Diabetologia 2015;58:1978-92.ArticlePubMedPDF
- 29. Peters LJ, Biessen EA, Hohl M, Weber C, van der Vorst EP, Santovito D. Small things matter: relevance of microRNAs in cardiovascular disease. Front Physiol 2020;11:793.ArticlePubMedPMC
- 30. Gomes CP, Oliveira GP, Madrid B, Almeida JA, Franco OL, Pereira RW. Circulating miR-1, miR-133a, and miR-206 levels are increased after a half-marathon run. Biomarkers 2014;19:585-9.ArticlePubMed
- 31. Li X, Li Y, Zhao L, Zhang D, Yao X, Zhang H, et al. Circulating muscle-specific miRNAs in Duchenne muscular dystrophy patients. Mol Ther Nucleic Acids 2014;3:e177.ArticlePubMedPMC
- 32. Toivonen JM, Manzano R, Olivan S, Zaragoza P, Garcia-Redondo A, Osta R. MicroRNA-206: a potential circulating biomarker candidate for amyotrophic lateral sclerosis. PLoS One 2014;9:e89065.ArticlePubMedPMC
- 33. Heinemann FG, Tolkach Y, Deng M, Schmidt D, Perner S, Kristiansen G, et al. Serum miR-122-5p and miR-206 expression: non-invasive prognostic biomarkers for renal cell carcinoma. Clin Epigenetics 2018;10:11.ArticlePubMedPMCPDF
- 34. Liu X, Zheng W, Zhang X, Dong M, Sun G. The diagnostic and prognostic value of serum miR-206 in colorectal cancer. Int J Clin Exp Pathol 2017;10:7528-33.PubMedPMC
- 35. Jin P, Gu W, Lai Y, Zheng W, Zhou Q, Wu X. The circulating microRNA-206 level predicts the severity of pulmonary hypertension in patients with left heart diseases. Cell Physiol Biochem 2017;41:2150-60.ArticlePubMedPDF
- 36. Infante T, Forte E, Punzo B, Cademartiri F, Cavaliere C, Soricelli A, et al. Correlation of circulating miR-765, miR-93-5p, and miR-433-3p to obstructive coronary heart disease evaluated by cardiac computed tomography. Am J Cardiol 2019;124:176-82.ArticlePubMed
- 37. Kim M, Zhang X. The profiling and role of miRNAs in diabetes mellitus. J Diabetes Clin Res 2019;1:5-23.PubMedPMC
- 38. Wu H, Zhang T, Pan F, Steer CJ, Li Z, Chen X, et al. MicroRNA-206 prevents hepatosteatosis and hyperglycemia by facilitating insulin signaling and impairing lipogenesis. J Hepatol 2017;66:816-24.ArticlePubMedPMC
- 39. Zhu L, Chen T, Ye W, Wang JY, Zhou JP, Li ZY, et al. Circulating miR-182-5p and miR-5187-5p as biomarkers for the diagnosis of unprotected left main coronary artery disease. J Thorac Dis 2019;11:1799-808.ArticlePubMedPMC
- 40. Weale CJ, Matshazi DM, Davids SF, Raghubeer S, Erasmus RT, Kengne AP, et al. Circulating miR-30a-5p and miR-182-5p in prediabetes and screen-detected diabetes mellitus. Diabetes Metab Syndr Obes 2020;13:5037-47.ArticlePubMedPMCPDF
- 41. Weale CJ, Matshazi DM, Davids SF, Raghubeer S, Erasmus RT, Kengne AP, et al. Expression profiles of circulating microRNAs in South African type 2 diabetic individuals on treatment. Front Genet 2021;12:702410.ArticlePubMedPMC
- 42. Karolina DS, Armugam A, Tavintharan S, Wong MT, Lim SC, Sum CF, et al. MicroRNA 144 impairs insulin signaling by inhibiting the expression of insulin receptor substrate 1 in type 2 diabetes mellitus. PLoS One 2011;6:e22839.ArticlePubMedPMC
- 43. Liu S, Lin Z, Zheng Z, Rao W, Lin Y, Chen H, et al. Serum exosomal microRNA-766-3p expression is associated with poor prognosis of esophageal squamous cell carcinoma. Cancer Sci 2020;111:3881-92.ArticlePubMedPMCPDF
- 44. Owczarz M, Polosak J, Domaszewska-Szostek A, Kolodziej P, Kurylowicz A, Puzianowska-Kuznicka M. Age-related epigenetic drift deregulates SIRT6 expression and affects its downstream genes in human peripheral blood mononuclear cells. Epigenetics 2020;15:1336-47.ArticlePubMedPMC
- 45. Hayakawa K, Kawasaki M, Hirai T, Yoshida Y, Tsushima H, Fujishiro M, et al. MicroRNA-766-3p contributes to anti-inflammatory responses through the indirect inhibition of NF-κB signaling. Int J Mol Sci 2019;20:809.ArticlePubMedPMC
- 46. Yang ZM, Chen LH, Hong M, Chen YY, Yang XR, Tang SM, et al. Serum microRNA profiling and bioinformatics analysis of patients with type 2 diabetes mellitus in a Chinese population. Mol Med Rep 2017;15:2143-53.ArticlePubMedPMC
- 47. Olioso D, Dauriz M, Bacchi E, Negri C, Santi L, Bonora E, et al. Effects of aerobic and resistance training on circulating microRNA expression profile in subjects with type 2 diabetes. J Clin Endocrinol Metab 2019;104:1119-30.ArticlePubMed
- 48. Cao Z, Yao F, Lang Y, Feng X. Elevated circulating LINC-P21 serves as a diagnostic biomarker of type 2 diabetes mellitus and regulates pancreatic β-cell function by sponging miR-766-3p to upregulate NR3C2. Exp Clin Endocrinol Diabetes 2022;130:156-64.ArticlePubMed
- 49. Demirsoy lH, Ertural DY, Balci S, Cınkır U, Sezer K, Tamer L, et al. Profiles of circulating MiRNAs following metformin treatment in patients with type 2 diabetes. J Med Biochem 2018;37:499-506.ArticlePubMedPMC
- 50. Ortega FJ, Mercader JM, Moreno-Navarrete JM, Rovira O, Guerra E, Esteve E, et al. Profiling of circulating microRNAs reveals common microRNAs linked to type 2 diabetes that change with insulin sensitization. Diabetes Care 2014;37:1375-83.ArticlePubMedPDF
Citations
Citations to this article as recorded by
- Glimepiride Compared to Liraglutide Increases Plasma Levels of miR-206, miR-182-5p, and miR-766-3p in Type 2 Diabetes Mellitus: A Randomized Controlled Trial (Diabetes Metab J 2023;47:668-81)
Nikolai N. Scherbak, Robert Kruse, Thomas Nyström, Johan Jendle
Diabetes & Metabolism Journal.2023; 47(6): 882. CrossRef - Glimepiride Compared to Liraglutide Increases Plasma Levels of miR-206, miR-182-5p, and miR-766-3p in Type 2 Diabetes Mellitus: A Randomized Controlled Trial (Diabetes Metab J 2023;47:668-81)
Da Young Lee
Diabetes & Metabolism Journal.2023; 47(6): 879. CrossRef