Effect of diabetes on exosomal miRNA profile in patients with obesity

Abstract

Introduction

Obesity is a risk factor for type 2 diabetes mellitus (T2DM) and cardiovascular disease. T2DM increases the risk of cardiovascular-related death. We investigated changes in circulating exosomal microRNA (miRNA) profiles in patients with DM with obesity compared with patients without DM with obesity.

Research design and methods

This prospective study involved 29 patients with obesity (patients without DM=16, patients with DM=13) and healthy volunteers (HVs) (n=18). We measured circulating levels of exosomal miRNAs by next-generation sequencing and compared miRNA levels across the three groups.

Results

The expression levels of 25 miRNAs (upregulated=14, downregulated=11) differed between patients with obesity with DM and patients with obesity without DM. Compared with HV, patients with DM with obesity had 53 dysregulated miRNAs. Additionally, moving stepwise from HV to patients with obesity without DM to patients with obesity with DM, there was a consistent increase in expression levels of miR-23a-5p and miR-6087 and a consistent decrease in expressions levels of miR-6751-3p.

Conclusions

Our data show that the exosomal miRNAs is altered by dysregulated glucose metabolism in patients with obesity. This circulating exosomal miRNA signature in patients with obesity with or without DM is a potential biomarker and therapeutic target in these patients.

Correspondence to Professor Soon Hyo Kwon; [email protected] ; Professor Seongho Ryu; [email protected]

Significance of this study

What is already known about this subject?

• MicroRNA (miRNA) changes in exsosome contribute to the development of diabetes or obesity.

What are the new findings?

• We found the similarity in the exosomal miRNA profiles of patients with obesity with or without diabetes.

• Our study demonstrated the specific circulating exosomal miRNAs of patients with obesity with diabetes compared with patients with obesity without diabetes and healthy volunteers.

• Our study showed miRNAs that change sequentially in patients with obesity without diabetes and patients with obesity with diabetes compared with healthy volunteers.

How might these results change the focus of research or clinical practice?

• This circulating exosomal miRNAs signature of patients with obesity and with or without diabetes is a potential biomarker and therapeutic target of obesity and obesity-related complications.

Introduction

Modern sedentary lifestyles and excessive calorie intake have increased the risk of developing metabolic diseases.¹ Obesity is a major component of metabolic syndrome, which includes a cluster of metabolic disturbances such as hypertension, dyslipidemia, hyperglycemia, insulin resistance, obesity and inflammation.² Obesity poses a substantial health risk, as it has been linked to several common diseases, such as type 2 diabetes mellitus (T2DM), cardiovascular disease, stroke, arthritis and several types of cancer.^{3 4} The upsurge in obesity rates appears to be the primary factor affecting the recent increase in the prevalence of T2DM.⁵ Although T2DM is associated with a strongly increased cardiovascular disease risk, the process of irreversible subclinical vascular damage begins in the prediabetes stage.⁶

T2DM creates a burden for patients and is a risk factor for many complications. Diabetes mellitus increases the risk of cardiovascular disease by threefold to fourfold in women and twofold to threefold in men after adjusting for other risk factors.⁷ There is a positive correlation between the volume of visceral fat tissue in patients with T2DM and arterial inflammation.⁸ Another study demonstrated that metabolic syndrome is associated with increased arterial stiffness in patients with T2DM, and abdominal obesity is associated increased arterial stiffness.⁹ Obesity and diabetes are risk factors for cardiovascular disease, and cardiovascular risk increases when both factors coexist.^8–10 Therefore, preventing T2DM in obesity is critical issue for preventing cardiovascular complications. T2DM is a progressive disease with a course that typically depends on β-cell dysfunction and worsening of insulin resistance.¹¹ However, the progression of diabetes from insulin resistance and the factors that play an important role in the process are not well-understood.

Recent studies have shown that exosomes are involved in multiple pathological processes, such as insulin resistance, lipid toxicity, dyslipidemia, endocrine disorders, hypercoagulable state and chronic inflammation.^12–14 Most cells in the human body can secrete exosomes containing mRNA, microRNA (miRNA) and other noncoding RNAs¹⁴; exosomal miRNAs can be stabilized in circulation via protection from the exosome, and then are transferred into target cells to inhibit or activate the expression of certain genes.^{15 16} The mechanisms by which obesity causes diabetes and vascular complications have been described. miRNAs affect adipose tissue metabolism (adipogenesis, release and action of signals derived from adipose tissue) and insulin secretion and activity, and an imbalance in these miRNAs can lead to the development of obesity and related complications.¹⁷ Indeed, dysregulation of miRNAs is a common feature in various human diseases including atherosclerosis and obese/diabetes-associated vascular complications.¹⁸ The miRNAs carried by exosomes in the body fluids may act as efficient biomarkers for diabetic complications, such as diabetic nephropathy.^{14 19} Several recent studies have evaluated miRNAs in diabetes and the prediabetes state²⁰; however, most studies focused on serum cell free circulating miRNAs using microarray or real-time quantitative PCR (qPCR) in prediabetic conditions and overt diabetes. Additionally, these studies did not focus on obesity.

Our previous work investigated the difference of circulating exosomal miRNA profiles between patients with obesity without T2DM and healthy volunteers (HVs).¹² However, the effects of persistent high glucose such as in DM on circulating exosomal miRNAs in patients with obesity has not been evaluated. Therefore, we compared the miRNA levels through small RNA sequencing to identify differentially expressed miRNAs to highlight potential pathways involved in the progression of T2DM in patients with obesity, as well as to provide new targets for treating obesity-associated diabetes.

Research design and methods

Study participants

We prospectively recruited HVs and subjects with obesity at Soonchunhyang University Hospital (Seoul, South Korea). All participant’s demographic data and anthropometric parameters were obtained at the time of enrollment. Obesity was defined as body mass index (BMI) ≥30.0 kg/m² according to the Steering Committee of the Regional Office for the Western Pacific Region of WHO. BMI was calculated by dividing the initial body weight (kg) by height (m) squared. DM was confirmed based on WHO criteria: fasting glucose ≥126 mg/dL, serum glucose ≥200 mg/dL after a 2-hour oral glucose tolerance test, hemoglobin A1c level ≥6.5%²¹ and diabetes medication history. Blood samples were obtained at least 8 hours after fasting, and the estimated glomerular filtration rate was calculated using the CKD Epidemiology Collaboration equation.²²

Circulating exosomal RNA isolation and assessment

We isolated exosomes from serum using the ExoQuick isolation agent (System Biosciences, Palo Alto, California, USA), according to the manufacturer’s guidelines. Serum samples were centrifuged at 3000 ×g for 15 min; subsequently, supernatants were mixed with ExoQuick reagent and incubated for 30 min at 4°C. After centrifugation, the supernatants were aspirated. The remaining pellets contained exosomes. We resuspended the pellets in 100–200 µL of sterile phosphate-buffered saline and RNA was extracted with the miRNeasy Mini Kit (Qiagen, Hilden, Germany). Exosome suspensions (200 µL) were mixed with QIAzol lysis buffer (1 mL) (Qiagen), and the mixtures were processed according to the manufacturer’s guidelines. The RNA was eluted in RNase-free water (20 µL). Purified RNA was analyzed for exosomal RNA size distribution using an Agilent 2100 Bioanalyzer with an RNA Pico Chip and Small RNA Chip (Agilent Technologies, Santa Clara, California, USA).

Complementary DNA library preparation and small RNA sequencing

RNA sequencing was conducted as described previously.¹² The samples were processed to produce exosomal RNA (10 ng) as input for each library. Small RNA libraries were constructed using the SMARTer smRNA-Seq Kit for Illumina (Takara Bio, Shiga, Japan) according to the manufacturer’s guidelines. We generated sequencing libraries by polyadenylation, complementary DNA synthesis and PCR amplification.

The libraries were gel-purified and validated based on their size, purity and concentration, using an Agilent Bioanalyzer. They were quantified by qPCR according to the qPCR Quantification Protocol Guide (KAPA Library Quantification Kits for Illumina Sequencing Platforms). We assessed the quality of libraries using a TapeStation D1000 ScreenTape (Agilent Technologies). Equimolar amounts of libraries were pooled and sequenced on an Illumina HiSeq 2500 instrument (Illumina, San Diego, California, USA) to generate 101 base reads. To perform image decomposition and quality value calculations, we used the modules in the Illumina pipeline. We conducted all procedures for next-generation sequencing analysis through Macrogen (Seoul, Korea).

Analysis of RNA sequencing data

Known and novel miRNAs were identified using the miRDeep2 software algorithm after sequence alignment. We retrieved the Homo sapiens reference genome release hg19 from the UCSC Genome Browser, which was indexed using Bowtie (1.1.2), a program for aligning experimental and reference sequences before sequencing alignment. The reads were then aligned to the mature and precursor H. sapiens miRNAs, obtained from miRBase 21.

Proportion of miRNAs and other RNAs

Uniquely clustered reads were sequentially aligned to the reference genome using miRBase 21 and non-coding RNA database Rfam 9.1 to identify known miRNAs and other types of RNAs, respectively.

Analysis of miRNA expression levels

The raw data (reads for each miRNA) were normalized by the relative log expression using DESeq2. For preprocessing, miRNAs absent from >50% of the samples were excluded, leaving mature miRNAs for analysis. We added one to the normalized read count of the filtered miRNAs and facilitated log2 transformation in order to draw a correlation plot. For each miRNA, the baseMean and log fold-change were calculated between case and control. For comparison of the two groups, we conducted a statistical hypothesis test by using the negative binomial Wald test in DESeq2, and defined the miRNAs that differentially expressed between the two groups as having a |fold change|≥2 and false discovery rate-adjusted p value <0.05. We also performed a hierarchical clustering analysis, using complete linkage and Euclidean distance as measures of similarity, to display the expression patterns of differentially expressed miRNAs that satisfied the criteria of fold-change|≥2 and false discovery rate-adjusted p value of<0.05. All data analyses and visualization of differentially expressed genes were conducted using R V.3.3.1 (www.r-project.org).

miRNA target prediction and functional annotation analysis

To evaluate the functions of the differentially expressed miRNAs, miRNA target prediction was performed using the miRSystem database. The miRSystem used seven algorithms for predicting miRNA targets (namely, DIANA-microT, miRanda, miRBridge, PicTar, PITA, RNA22 and TargetScan) and two experimentally validated databases of miRNA target genes (TarBase and miRecords). The gene list of predicted targets (prediabetes only or diabetes only miRNA) submitted to the Database for Annotation, Visualization, and Integrated Discovery and the putative targets were annotated using four pathway databases, including Kyoto Encyclopedia of Genes and Genomes, Biocarta, Pathway Interaction Database and Reactome.

Statistical analyses

Continuous variables are expressed as the mean±SD, and categorical variables are expressed as a number and percentage. Variables without normal distribution are expressed as the median and IQR. Continuous variables of HVs and patients with obesity were analyzed by the Mann-Whitney U test. Categorical and non-parametric variables were compared using a χ² test and Kruskal-Wallis test, as appropriate. Statistical analyses were performed using SPSS V.23.0 (IBM, Armonk, New York, USA).

Results

Baseline characteristics

A total of 47 participants (18 HVs, 16 patients with obesity without DM, and 13 patients with obesity with DM), who voluntarily provided informed consent, were enrolled. Baseline characteristics of the participants are shown in table 1. All measurements except diastolic blood pressure, blood urea nitrogen, creatinine, urinary protein and urinary albumin differed between HVs and patients with obesity without DM. The patients with obesity with DM and the HVs were older compared with the patients with obesity without DM. Gender, blood pressure, weight, BMI, waist:hip ratio, serum creatinine, cystatin C, lipid profiles and C-peptide did not differ between the groups with obesity. No participant had hypertension, with the exception of 10 patients (76.9%) with DM and obesity. No other participants had cerebrovascular disease or cancer.

Table 1

Baseline characteristics of study participants

Variables	HVs (n=18)	Non-diabetes (n=16)	Diabetes (n=13)	P value (HVs vs non-diabetes)
Gender (male n, %)	5 (27.8%)	7 (43.8%)	4 (30.8%)	0.331
Age (years)	38.61±7.96	31.38±8.76*	42.77±13.03	0.014
Weight (kg)	57.97±5.87	115.48±24.91	109.84±32.08†	<0.001
BMI (kg/m2)	21.66±1.43	39.93±4.68	37.95±6.43†	<0.001
Systolic BP (mm Hg)	106.11±8.49	127.50±12.38	129.85±9.72†	<0.001
Serum creatinine (mg/dL)	0.76±0.12	0.73±0.14	0.88±0.39	0.369
eGFR by CKD-EPI (mL/min)	107.14±10.89	115.65±11.03*	95.14±30.44	0.031
Total cholesterol (mg/dL)	183.89±24.15	222.19±46.45	191.85±37.36	0.007
HDL cholesterol (mg/dL)	65.33±14.57	47.31±9.76	50.23±13.86†	<0.001
LDL cholesterol (mg/dL)	112.77±23.02	152.38±44.40	125.08±40.85	0.004
Triglyceride (mg/dL)	105.44±82.86	208.56±120.93	196.46±72.50†	0.002
Glucose (mg/dL)	97.00±8.51	108.00±16.47*	161.46±42.32†	0.02
HbA1c (mg/dL)		5.70±0.28	7.90±1.68*
C-peptide (ng/mL)		4.73±1.95	3.84±1.52
Urinary albumin (mg/day)	3.67±3.31†	11.40±15.98	102.60±210.65*	0.037

*P<0.05 (non-diabetes vs diabetes).

†P<0.05 (HVs vs diabetes).

BMI, body mass index; BP, blood pressure; CKD-EPI, chronic kidney disease-improved prediction equations; eGFR, estimated glomerular filtration rate; HbA1c, hemoglobin A1c; HDL, high-density lipoprotein; HVs, healthy volunteers; LDL, low-density lipoprotein.

Profiles of exosomal miRNAs among participants

According to the multidimensional scaling graph, the patients with obesity with DM and patients with obesity without DM had many overlapping areas (online supplementary figure S1). After small RNA sequencing, we identified 26 upregulated and 28 downregulated miRNAs in patients with obesity without DM compared with HVs (figure 1A,B). We also found 24 upregulated and 29 downregulated miRNAs in patients with obesity with DM compared with those of HVs (figure 2A,B). Comparison of the patients with obesity without DM and patients with obesity with DM revealed 14 upregulated miRNAs and 11 downregulated miRNAs in patients with obesity with DM (figure 3A,B). To identify unique miRNAs in patients with obesity without DM or patients with obesity with DM, we selected miRNAs upregulated only in one group compared with HVs (table 2 and online supplementary figures S2 and S3).

Figure 1. Exosomal microRNAs (miRNAs) from patients with obesity without DM differ from those from healthy volunteers (HVs). (A) Heat map showing z-score of exosomal miRNAs from HVs (n=18) and patients with obesity without DM (n=16) with 26 upregulated (yellow) and 28 downregulated (blue) miRNAs. (B) Fold change (FC) of top 10 miRNA expression in patients with obesity without DM compared with HVs. DM, diabetes mellitus.

Figure 2. Exosomal microRNAs (miRNAs) from patients with obesity with DM differ from those from patients with obesity without DM. (A) Heat map showing z-score of exosomal miRNAs from patients with obesity without DM (n=16) and patients with obesity with DM (n=13) with 14 upregulated (yellow) and 11 downregulated (blue) miRNAs. (B) Fold change (FC) of top 10 miRNA expression in patients with obesity with DM compared with patients with obesity without DM. DM, diabetes mellitus.

Figure 3. Exosomal microRNAs (miRNAs) from patients with obesity with DM differ from those from healthy volunteers (HVs). (A) Heat map showing z-score of exosomal miRNAs from HVs (n=18) and patients with obesity with DM (n=13) with 24 upregulated (yellow) and 29 downregulated (blue) miRNAs. (B) Fold change (FC) of top 10 miRNA expression in patients with obesity with DM compared with HVs. DM, diabetes mellitus.

Table 2

Upregulated exosomal miRNAs only in patients without diabetes/patients with diabetes compared with HVs

	Non-diabetes/HVs			Diabetes/HVs
hsa-let-7c-5p	2.1	<0.001	hsa-miR-150–3 p	2.16	<0.05
hsa-let-7d-3p	2.13	<0.001	hsa-let-7b-5p	2.18	<0.001
hsa-miR-4488	2.36	<0.05	hsa-miR-664a-5p	2.23	<0.05
hsa-miR-130b-5p	2.41	<0.05	hsa-miR-320b	2.28	<0.01
hsa-miR-29b-3p	2.43	<0.01	hsa-miR-423–5 p	2.32	<0.001
hsa-miR-151a-3p	2.43	<0.001	hsa-miR-483–5 p	2.52	<0.01
hsa-miR-556–3 p	2.52	<0.05	hsa-miR-320a	2.52	<0.01
hsa-miR-21–5 p	2.57	<0.001	hsa-miR-3960	2.57	<0.01
hsa-let-7b-3p	2.74	<0.01	hsa-miR-424–3 p	3.25	<0.01
hsa-miR-192–5 p	2.74	<0.001
hsa-miR-183–5 p	2.85	<0.01
hsa-miR-378g	2.9	<0.05
hsa-miR-126–5 p	3.03	<0.001
hsa-let-7f-1–3 p	4.23	<0.001
hsa-let-7a-3p	4.72	<0.001

HVs, healthy volunteers; miRNAs, microRNAs.

Exosomal miRNAs associated with T2DM development in obesity

We identified the specific miRNAs that continuously changed from HVs to patients with obesity without and with diabetes. Expression of miR-23a-5p in patients with obesity without DM increased by 2.12-fold compared with HVs; the value in the patients with obesity with DM was 2.64-fold higher than that in patients with obesity without DM. Additionally, miR-6087 levels in patients with obesity without DM were 2.64-fold higher than those in HVs, whereas these values in patients with obesity with DM increased 2.03-fold compared with those patients with obesity without DM. In contrast, miR-6751-3p levels in patients with obesity without DM were 4.03-fold lower than those in HVs, whereas these values in patients with obesity with DM decreased 2.86-fold compared with those in patients with obesity without DM.

Prediction of targets and functional annotation of obesity with or without diabetes-specific miRNAs

We identified 15 obesity-specific miRNAs and 9 diabetes-specific miRNAs with significantly different expression compared with HVs (table 2). Targets of these core miRNAs were predicted using miRSystem. Obesity-specific miRNAs and diabetes-specific miRNAs were predicted to target 1543 and 2405 genes, respectively. The possible pathways targeted by these miRNAs are presented in online supplementary figures S4 and S5.

Discussion

In this study, we found that the miRNA profiles in exosomes of patients with obesity with DM were similar to those in patients with obesity without DM compared with HVs. These results suggest that obesity has an overwhelming effect on exosomal miRNAs and that these changes compared with individuals without obesity contribute to metabolic complications. We also identified the miRNAs that differed according to presence of diabetes, which may play an important role in the pathological mechanism by which patients with insulin resistance develop diabetes.

Recent in vivo studies demonstrated that exosomes can transfer mature miRNAs between organs, resulting in functional changes in the receiving cells and affecting whole-body insulin sensitivity.^{13 23} The transfer of obesity-associated exosomal miRNAs induced glucose intolerance, adipose inflammation and hepatic steatosis in lean mice.¹³ These results support a role for obesity-related exosomal miRNAs in metabolic complications. The similarity of exosomal miRNAs in the patients with obesity without DM and patients with obesity with DM in our study also provide insight into the pathogenesis of vascular complications in obesity. miR-122–5 p, miR-193b-5P, miR-193a-5P and miR-1290 were the most upregulated miRNAs in the patients with obesity compared with HVs, whereas miR-1273a and miR-20a-5p were the most downregulated miRNAs. These miRNAs may contribute to the pathogenesis of vascular complications in the patients with obesity and to insulin resistance, thus contributing to the development of diabetes. For example, miR-122–5 p is known to be associated with the prognosis of acute myocardial infarction²⁴ and is a fatty liver and cardiovascular risk factor.²⁵ One study reported that miR-20a-5p is significantly decreased in gestational diabetes and is associated with glucose homeostasis, including insulin signaling, phosphatidyl-inositol-3-kinase (PI3K)-protein kinase B (Akt), and forkhead box protein O (FoxO).²⁶ Both of miR-193a-5P and miR-193b-5P have been reported to contribute to tumor suppressor in various cancers.²⁷ miR-1290 has also been studied mainly in cancer, and most studies showed that overexpression of miR-1290 increases cancer cell proliferation and invasion.²⁸ Reduction of miR-1273a is also known to increase the activity of circular PIP5K1A associated with progression of colon cancer.²⁹ These results support that obesity is associated with cancer development via miRNAs and that these miRNAs should be validated in a cancer cohort.

We determined the exosomal levels of 15 obesity-specific miRNAs. One-third of these miRNAs belongs to the let-7 family (let-7b-3p, let-7f-1-3p, let-7d-3p, let-7a-3p and let-7c-5p), which is downregulated in the development of diabetes and vascular complications.³⁰ Although miR-126-5p is known to have a protective effect on atherosclerosis, miR-192-5p has been shown to protect the kidney from damage caused by hypertension.^{31 32} Urinary exosomal miR-21-5p is increased in diabetic nephropathy, and is associated with kidney damage caused by diabetes.³³ In addition, miR-29b-3p was found to be specific to islet beta-cell stress.³⁴ This is the first study to report that the other six miRNAs (miR-378g, miR-556-3p, miR-151a-3p, miR-130b-5p, miR-4488 and miR-183-5p) are associated with metabolic syndrome or diabetes. The changes in these miRNAs indicate that vascular complications are already occurring in patients with obesity without DM. Further studies are needed to determine whether the newly discovered miRNAs play a role in the promotion or inhibition of the progression of metabolic syndrome to diabetes.

Additionally, there were nine diabetes-specific miRNAs, which are thought to play a key role in the development of diabetes and vascular complication. Some miRNAs were reported to be associated with hyperglycemia or vascular disease. In a study comparing patients with T2DM with healthy individuals, miR-3960 was one of five miRNAs that was significantly different in serum.³⁵ One study showed that N-acetyl-seryl-aspartyl-lysyl-proline (AcSDKP) associated with the fibroblast growth factor receptor 1 signaling pathway maintained endothelial mitochondrial homeostasis by inducing miR-let-7b-5p expression.³⁶ Additionally, overexpression of miR-320a is a key risk factor for atherosclerosis³⁷ and another study showed that overexpression of miR-320b abolishes enhanced mesenchymal stem cell-medicated cardioprotection.³⁸ Five other diabetes-specific miRNAs have not been reported to be directly related to obesity or diabetes.

Two miRNAs were observed to be sequentially increased in patients with obesity without DM and patients with obesity with DM compared with HVs; the only sequentially reduced miRNA was considered responsible for the mechanism of diabetes development. One of these miRNAs, miR-23a-5p, has been found to be associated with atherosclerosis plaque progression and vulnerability by promoting foam cell formation. This miRNA attenuates cholesterol efflux to promote foam cell formation by downregulating the expression of ATP-binding cassette transporters (ABCA1 and ABCG1).³⁹ The other, miR-6087, is known to inhibit endoglin expression, a type of transforming growth factor-β receptor on endothelial cells. It is known to play an important role in angiogenesis,⁴⁰ and is associated with expression of the tumor suppressor gene p53.⁴¹ Recent studies showed that p53 gene polymorphisms affect life expectancy and lifestyle-related diseases, such as T2DM and obesity, indicating a relationship between p53 function and metabolic disorders.⁴² No studies have examined the correlation between miR-6751-3p and diabetes, and the function of miR-6751-3p is unknown. Diabetes does not develop in all patients with obesity, but the continuously changing miRNAs identified in our study appear to play an important role in diabetes development in patients with obesity. Further studies are needed to clarify the functions of miRNA in diabetes.

Pathway analysis related to miRNAs revealed common biological pathways in both the patients with obesity with DM and patients with obesity without DM. As expected, diabetes-specific miRNAs significantly affect the pathways associated with insulin secretion, insulin regulation and energy metabolism. In the patients with obesity without DM, miRNAs significantly affected pathways associated with insulin secretion, inflammatory cytokine, interferon, muscle atrophy and vascular smooth muscle atrophy. The insulin receptor, WNT, AMP-activated protein kinase, transforming growth factor-β signaling and mitogen-activated protein kinase (MAPK) signaling pathways were among the top canonical pathways, identified in previous studies of metabolic syndrome.^{12 43} Pathway analysis in our current study was similar to that in our previous study.^{12 19} Representative examples of important pathways in our study include the MAPK signaling and integrin pathways. MAPK signaling pathway, axon guidance signaling pathway, neurotrophin signaling pathway and notch signaling pathway are associated with insulin resistance.⁴⁴ MAPK signaling is known to play an important role in diabetic nephropathy.¹⁹ In addition, the integrin-vascular endothelial growth factor axis is a key factor in pathological angiogenesis, including in tumor metastasis, tissue remodeling and diabetic complications.¹⁹ Further studies are required to evaluate the association between these miRNAs and the biological pathways identified, as well as the mechanisms and critical points in the development of diabetes from metabolic syndrome.

There were some limitations to our study. First, a relatively small number of participants was enrolled in this study. Further studies with a larger cohort are warranted to confirm our findings. Second, some demographic factors and clinical parameters differed between groups, which may have affected the circulating exosomal miRNA profiles. Third, although we conducted exosomal RNA sequencing using validated commercial kits, the procedures used to isolate the exosomes are not standardized, and systematic bias in library preparation may lead to false-positive identification of miRNAs.

Conclusions

We identified specific miRNAs and associated biological pathways in patients with obesity without DM and patients with obesity with DM. Moreover, specific exosomal miRNAs of patients with obesity without DM and patients with obesity with DM may provide clues to the mechanism of hyperglycemic complications. Adipokines, inflammation and reactive oxygen species in patients with obesity may affect the cells in the body, resulting in exosomal miRNA changes, leading to developing diabetes and vascular complications; however, further research is needed to confirm this. In conclusion, these miRNAs can be considered as candidates for treating metabolic disease and diabetes.

Footnote

HyosK and Y-UB contributed equally.

Contributors HyosK contributed to data analysis, statistical analysis, literature search and writing the manuscript. Y-UB involved in data analysis, conduction of the experiments, statistical analysis and cowriting the manuscript. HL, HyouK, JSJ, HN, DCH and DWB involved in data analysis, searching literature and editing the manuscript. SaHK and HKP involved in study design, data collection and editing the manuscript. SR and SoHK contributed to study design, submitting the grant, conduction of the experiments, data collection and analysis, literature search and editing the manuscript.

Funding This study was supported by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Korean government (MSIT) (NRF-2019M3E5D3073092) and the Soonchunhyang University Research Fund (20191124).

Competing interests None declared.

Patient consent for publication Not required.

Ethics approval This study was conducted in accordance with the Declaration of Helsinki, and the study protocol was approved by the institutional review board of Soonchunhyang University Seoul Hospital (IRB No: 2015-11-020).

Provenance and peer review Not commissioned; externally peer reviewed.

Data availability statement Data are available on reasonable request.

References

¹ Lee M-W, Lee M, Oh K-J. Adipose tissue-derived signatures for obesity and type 2 diabetes: adipokines, Batokines and microRNAs. J Clin Med 2019; 8. doi: doi:10.3390/jcm8060854 31208019

² Grundy SM, Brewer HB, Cleeman JI, et al. Definition of metabolic syndrome: report of the National heart, lung, and blood Institute/American heart association conference on scientific issues related to definition. Circulation 2004; 109: 433–8. doi:10.1161/01.CIR.0000111245.75752.C6 http://www.ncbi.nlm.nih.gov/pubmed/14744958

³ Pescador N, Pérez-Barba M, Ibarra JM, et al. Serum circulating microRNA profiling for identification of potential type 2 diabetes and obesity biomarkers. PLoS One 2013; 8: e77251. doi:10.1371/journal.pone.0077251 http://www.ncbi.nlm.nih.gov/pubmed/24204780

⁴ Kwon SH. Changes in kidney function markers after bariatric surgery in morbidly obese patients. Kidney Res Clin Pract 2020; 39: 115–7. doi:10.23876/j.krcp.20.058 http://www.ncbi.nlm.nih.gov/pubmed/32522893

⁵ Sonmez A, Yumuk V, Haymana C, et al. Impact of obesity on the metabolic control of type 2 diabetes: results of the Turkish nationwide survey of glycemic and other metabolic parameters of patients with diabetes mellitus (TEMD obesity study). Obes Facts 2019; 12: 167–78. doi:10.1159/000496624 http://www.ncbi.nlm.nih.gov/pubmed/30893706

⁶ Pickup JC. Inflammation and activated innate immunity in the pathogenesis of type 2 diabetes. Diabetes Care 2004; 27: 813–23. doi:10.2337/diacare.27.3.813 http://www.ncbi.nlm.nih.gov/pubmed/14988310

⁷ Norhammar A, Schenck-Gustafsson K. Type 2 diabetes and cardiovascular disease in women. Diabetologia 2013; 56: 1–9. doi:10.1007/s00125-012-2694-y http://www.ncbi.nlm.nih.gov/pubmed/22945305

⁸ Reijrink M, de Boer SA, Spoor DS, et al. Visceral adipose tissue volume is associated with premature atherosclerosis in early type 2 diabetes mellitus independent of traditional risk factors. Atherosclerosis 2019; 290: 87–93. doi:10.1016/j.atherosclerosis.2019.09.016 http://www.ncbi.nlm.nih.gov/pubmed/31604171

⁹ Levisianou D, Melidonis A, Adamopoulou E, et al. Impact of the metabolic syndrome and its components combinations on arterial stiffness in type 2 diabetic men. Int Angiol 2009; 28: 490–5. http://www.ncbi.nlm.nih.gov/pubmed/20087288

¹⁰ Daousi C, Casson IF, Gill GV, et al. Prevalence of obesity in type 2 diabetes in secondary care: association with cardiovascular risk factors. Postgrad Med J 2006; 82: 280–4. doi:10.1136/pmj.2005.039032 http://www.ncbi.nlm.nih.gov/pubmed/16597817

¹¹ Fonseca VA. Defining and characterizing the progression of type 2 diabetes. Diabetes Care 2009; 32 Suppl 2: S151–6. doi:10.2337/dc09-S301 http://www.ncbi.nlm.nih.gov/pubmed/19875543

¹² Bae Y-U, Kim Y, Lee H, et al. Bariatric surgery alters microRNA content of circulating exosomes in patients with obesity. Obesity 2019; 27: 264–71. doi:10.1002/oby.22379 http://www.ncbi.nlm.nih.gov/pubmed/30624857

¹³ Castaño C, Kalko S, Novials A, et al. Obesity-Associated exosomal miRNAs modulate glucose and lipid metabolism in mice. Proc Natl Acad Sci U S A 2018; 115: 12158–63. doi:10.1073/pnas.1808855115 http://www.ncbi.nlm.nih.gov/pubmed/30429322

¹⁴ Yao Z-Y, Chen W-B, Shao S-S, et al. Role of exosome-associated microRNA in diagnostic and therapeutic applications to metabolic disorders. J Zhejiang Univ Sci B 2018; 19: 183–98. doi:10.1631/jzus.B1600490 http://www.ncbi.nlm.nih.gov/pubmed/29504312

¹⁵ de Jong OG, Verhaar MC, Chen Y, et al. Cellular stress conditions are reflected in the protein and RNA content of endothelial cell-derived exosomes. J Extracell Vesicles 2012; 1. doi: doi:10.3402/jev.v1i0.18396 24009886

¹⁶ Xiao M, Li J, Li W, et al. Micrornas activate gene transcription epigenetically as an enhancer trigger. RNA Biol 2017; 14: 1326–34. doi:10.1080/15476286.2015.1112487 http://www.ncbi.nlm.nih.gov/pubmed/26853707

¹⁷ Landrier J-F, Derghal A, Mounien L. Micrornas in obesity and related metabolic disorders. Cells 2019; 8. doi: doi:10.3390/cells8080859 31404962

¹⁸ Araldi E, Suárez Y. Micrornas as regulators of endothelial cell functions in cardiometabolic diseases. Biochim Biophys Acta 2016; 1861: 2094–103. doi:10.1016/j.bbalip.2016.01.013 http://www.ncbi.nlm.nih.gov/pubmed/26825686

¹⁹ Kim H, Bae Y-U, Jeon JS, et al. The circulating exosomal microRNAs related to albuminuria in patients with diabetic nephropathy. J Transl Med 2019; 17: 236. doi:10.1186/s12967-019-1983-3 http://www.ncbi.nlm.nih.gov/pubmed/31331349

²⁰ Jiménez-Lucena R, Camargo A, Alcalá-Diaz JF, et al. A plasma circulating miRNAs profile predicts type 2 diabetes mellitus and prediabetes: from the CORDIOPREV study. Exp Mol Med 2018; 50: 1–12. doi:10.1038/s12276-018-0194-y http://www.ncbi.nlm.nih.gov/pubmed/30598522

²¹ WHO. Diabetes. Available: https://www.who.int/nmh/publications/fact_sheet_diabetes_en.pdf?ua=1

²² Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med 2009; 150: 604–12. doi:10.7326/0003-4819-150-9-200905050-00006 http://www.ncbi.nlm.nih.gov/pubmed/19414839

²³ Ying W, Riopel M, Bandyopadhyay G, et al. Adipose tissue macrophage-derived exosomal miRNAs can modulate in vivo and in vitro insulin sensitivity. Cell 2017; 171: 372–84. doi:10.1016/j.cell.2017.08.035 http://www.ncbi.nlm.nih.gov/pubmed/28942920

²⁴ Wang Y, Chang W, Zhang Y, et al. Circulating miR-22-5p and miR-122-5p are promising novel biomarkers for diagnosis of acute myocardial infarction. J Cell Physiol 2019; 234: 4778–86. doi:10.1002/jcp.27274 http://www.ncbi.nlm.nih.gov/pubmed/30256407

²⁵ Raitoharju E, Seppälä I, Lyytikäinen L-P, et al. Blood hsa-miR-122-5p and hsa-miR-885-5p levels associate with fatty liver and related lipoprotein metabolism-The young Finns study. Sci Rep 2016; 6: 38262. doi:10.1038/srep38262 http://www.ncbi.nlm.nih.gov/pubmed/27917915

²⁶ Pheiffer C, Dias S, Rheeder P, et al. Decreased expression of circulating miR-20a-5p in South African women with gestational diabetes mellitus. Mol Diagn Ther 2018; 22: 345–52. doi:10.1007/s40291-018-0325-0 http://www.ncbi.nlm.nih.gov/pubmed/29556924

²⁷ Khordadmehr M, Shahbazi R, Sadreddini S, et al. miR-193: a new weapon against cancer. J Cell Physiol 2019; 234: 16861–72. doi:10.1002/jcp.28368 http://www.ncbi.nlm.nih.gov/pubmed/30779342

²⁸ Zhang WC, Chin TM, Yang H, et al. Tumour-Initiating cell-specific miR-1246 and miR-1290 expression converge to promote non-small cell lung cancer progression. Nat Commun 2016; 7: 11702. doi:10.1038/ncomms11702 http://www.ncbi.nlm.nih.gov/pubmed/27325363

²⁹ Zhang Q, Zhang C, Ma J-X, et al. Circular RNA PIP5K1A promotes colon cancer development through inhibiting miR-1273a. World J Gastroenterol 2019; 25: 5300–9. doi:10.3748/wjg.v25.i35.5300 http://www.ncbi.nlm.nih.gov/pubmed/31558874

³⁰ Brennan E, Wang B, McClelland A, et al. Protective effect of let-7 miRNA family in regulating inflammation in diabetes-associated atherosclerosis. Diabetes 2017; 66: 2266–77. doi:10.2337/db16-1405 http://www.ncbi.nlm.nih.gov/pubmed/28487436

³¹ Baker MA, Wang F, Liu Y, et al. MiR-192-5p in the kidney protects against the development of hypertension. Hypertension 2019; 73: 399–406. doi:10.1161/HYPERTENSIONAHA.118.11875 http://www.ncbi.nlm.nih.gov/pubmed/30595117

³² Schober A, Nazari-Jahantigh M, Wei Y, et al. MicroRNA-126-5p promotes endothelial proliferation and limits atherosclerosis by suppressing DLK1. Nat Med 2014; 20: 368–76. doi:10.1038/nm.3487 http://www.ncbi.nlm.nih.gov/pubmed/24584117

³³ Zang J, Maxwell AP, Simpson DA, et al. Differential expression of urinary exosomal microRNAs miR-21-5p and miR-30b-5p in individuals with diabetic kidney disease. Sci Rep 2019; 9: 10900. doi:10.1038/s41598-019-47504-x http://www.ncbi.nlm.nih.gov/pubmed/31358876

³⁴ Saravanan PB, Vasu S, Yoshimatsu G, et al. Differential expression and release of exosomal miRNAs by human islets under inflammatory and hypoxic stress. Diabetologia 2019; 62: 1901–14. doi:10.1007/s00125-019-4950-x http://www.ncbi.nlm.nih.gov/pubmed/31372667

³⁵ Ding L, Ai D, Wu R, et al. Identification of the differential expression of serum microRNA in type 2 diabetes. Biosci Biotechnol Biochem 2016; 80: 461–5. doi:10.1080/09168451.2015.1107460 http://www.ncbi.nlm.nih.gov/pubmed/26554942

³⁶ Hu Q, Li J, Nitta K, et al. FGFR1 is essential for N-acetyl-seryl-aspartyl-lysyl-proline regulation of mitochondrial dynamics by upregulating microRNA let-7b-5p. Biochem Biophys Res Commun 2018; 495: 2214–20. doi:10.1016/j.bbrc.2017.12.089 http://www.ncbi.nlm.nih.gov/pubmed/29269295

³⁷ Chen C, Wang Y, Yang S, et al. MiR-320a contributes to atherogenesis by augmenting multiple risk factors and down-regulating SRF. J Cell Mol Med 2015; 19: 970–85. doi:10.1111/jcmm.12483 http://www.ncbi.nlm.nih.gov/pubmed/25728840

³⁸ Xu B, Luo Y, Liu Y, et al. Platelet-Derived growth factor-BB enhances MSC-mediated cardioprotection via suppression of miR-320 expression. Am J Physiol Heart Circ Physiol 2015; 308: H980–9. doi:10.1152/ajpheart.00737.2014 http://www.ncbi.nlm.nih.gov/pubmed/25724494

³⁹ Yang S, Ye Z-M, Chen S, et al. MicroRNA-23a-5p promotes atherosclerotic plaque progression and vulnerability by repressing ATP-binding cassette transporter A1/G1 in macrophages. J Mol Cell Cardiol 2018; 123: 139–49. doi:10.1016/j.yjmcc.2018.09.004 http://www.ncbi.nlm.nih.gov/pubmed/30227118

⁴⁰ Yoo JK, Kim J, Choi S-J, et al. Discovery and characterization of novel microRNAs during endothelial differentiation of human embryonic stem cells. Stem Cells Dev 2012; 21: 2049–57. doi:10.1089/scd.2011.0500 http://www.ncbi.nlm.nih.gov/pubmed/22142236

⁴¹ Zhang Y, Kang R, Liu W, et al. Identification and analysis of p53-mediated competing endogenous RNA network in human hepatocellular carcinoma. Int J Biol Sci 2017; 13: 1213–21. doi:10.7150/ijbs.21502 http://www.ncbi.nlm.nih.gov/pubmed/29104512

⁴² Hashimoto N, Nagano H, Tanaka T. The role of tumor suppressor p53 in metabolism and energy regulation, and its implication in cancer and lifestyle-related diseases. Endocr J 2019; 66: 485–96. doi:10.1507/endocrj.EJ18-0565 http://www.ncbi.nlm.nih.gov/pubmed/31105124

⁴³ Ramzan F, D'Souza RF, Durainayagam BR, et al. Circulatory miRNA biomarkers of metabolic syndrome. Acta Diabetol 2020; 57: 203–14. doi:10.1007/s00592-019-01406-6 http://www.ncbi.nlm.nih.gov/pubmed/31435783

⁴⁴ Chen G-Q, Lian W-J, Wang G-M, et al. Altered microRNA expression in skeletal muscle results from high-fat diet-induced insulin resistance in mice. Mol Med Rep 2012; 5: 1362–8. doi:10.3892/mmr.2012.824 http://www.ncbi.nlm.nih.gov/pubmed/22395471