ORIGINAL RESEARCH article
Micro-RNAs 518d-3p and 618 Are Upregulated in Individuals With Type 1 Diabetes With Multiple Microvascular Complications
- 1Laboratório de Carboidratos e Radioimunoensaio, LIM-18, Hospital das Clinicas HCFMUSP, Faculdade de Medicina da Universidade de São Paulo, São Paulo, Brazil
- 2Divisão de Oftalmologia do Hospital das Clinicas, HCFMUSP, Faculdade de Medicina da Universidade de São Paulo, São Paulo, Brazil
- 3Laboratório de Lipides, LIM-10, Hospital das Clinicas HCFMUSP, Faculdade de Medicina da Universidade de São Paulo, São Paulo, Brazil
- 4Programa de Pós-Graduação em Medicina, Universidade Nove de Julho (UNINOVE), São Paulo, Brazil
- 5Departamento de Fisiologia e Biofisica, Instituto de Ciencias Biomedicas, Universidade de São Paulo, São Paulo, Brazil
- 6Divisão de Endocrinologia do Hospital das Clinicas, HCFMUSP, Faculdade de Medicina da Universidade de São Paulo, São Paulo, Brazil
Objective: To compare the serum micro-RNAs (miRNAs) profile of individuals with type 1 diabetes without microvascular complications vs. those with multiple severe microvascular complications, in order to identify epigenetically modulated pathways in these two groups of individuals.
Research Design and Methods: A total of 10 subjects were selected among individuals followed in the Diabetes Outpatient Clinic and sorted according to the absence or presence of all microvascular complications. Samples from these participants were used for evaluation of serum miRNA expression profile employing a qRT-PCR assay with hydrolysis probes based on the Taqman Low Density Arrays (TLDA) system. The top six most differentially expressed miRNAs between the aforementioned groups were validated by qRT-PCR in additional 47 type 1 diabetes individuals sorted according to the absence or presence of all microvascular complications and matched for age, sex, degree of metabolic control, diabetes duration, and age at diagnosis.
Results: Twenty one out of three hundred and seventy seven miRNAs were upregulated in the group of individuals with all microvascular complications vs. the group without complications. The following miRs were validated: 518-3p, 34a-5p, 126-5p, 425-5p, 618, and 139-5p and logistic regression analyses showed that miRNA-518-3p and miRNA-618 were positively associated with multiple microvascular complications after adjustment for age, sex, diabetes duration, HbA1c and use of statin, angiotensin-converting enzyme inhibitors and amlodipine.
Conclusions: In this cohort of type 1 diabetes individuals, serum miR-518d-3p and miR-618 were upregulated in those with diabetes kidney disease, diabetes retinopathy, peripheral neuropathy, and cardiovascular autonomic neuropathy in comparison to individuals with no microvascular complications.
Diabetes kidney disease (DKD), diabetic retinopathy (DR), peripheral neuropathy (PN), and cardiovascular autonomic neuropathy (CAN) have hyperglycemia as the main etiopathogenic factor for their development (1). Several studies have already shown that periods of hyperglycemia result in permanent abnormalities in the target tissues of complications, a phenomenon known as “metabolic memory” (2, 3). Post-translational modifications of histones, cytosine methylation, and the action of micro-RNAs (miRNAs) explain deleterious effects of hyperglycemia even after improved glycemic control (4).
miRNAs constitute a class of small RNA with about 19 to 25 non-coding nucleotides that is the main incorporated component of the RNA-induced silencing complex (RISC), which represses translation or degrades target mRNAs (5). Besides acting as post-transcriptional regulators of gene expression, miRNAs can indirectly modulate DNA methylation and even stimulate transcription after binding to promoter regions of target-genes (5).
miRNAs are present in the bloodstream, possibly originated from passive leakage of damaged/dead cells or from active cell secretion, either in exosomes or linked to proteins (6). Secreted miRNAs may play a role in intercellular communication, so that donor cells can affect gene expression of distant or adjacent target-cells. Thus, beyond acting as hormones, serum miRNAs can be functional biomarkers for pathological conditions (7). Several studies have shown that a significant portion of the circulating miRNAs is not in microvesicles but it is associated with the Argonauta 2 (AGO2) protein, the main protein of the RISC complex, and this association explains the stability of circulating miRNAs (8, 9).
Most studies with miRNAs in the context of diabetes have been performed in animal models or in cell culture, but there are already studies addressing the association of serum miRNAs with one or more microvascular complication (10–14). However, there are no studies comparing serum miRNAs between type 1 diabetes individuals without vs. with all microvascular complications. Thus, in the present study we compared the serum miRNA profile of individuals with type 1 diabetes without microvascular complications vs. those with multiple severe microvascular complications: DKD, DR, PN, and CAN, in order to identify epigenetically modulated pathways in these two groups of individuals.
Materials and Methods
A total of 10 individuals were selected among type 1 diabetes individuals followed in the Diabetes Outpatient Clinic of Hospital das Clinicas da Faculdade de Medicina da Universidade de Sao Paulo and sorted according to the absence or presence of microvascular complications: (1) n = 5 individuals without DKD (estimated glomerular filtration rate [eGFR] > 90 mL/min/1.73 m2 and urinary albumin excretion <30 mg/g creatinine), without DR, without PN, and without CAN; and (2) n = 5 individuals with DKD (eGFR < 60 mL/min /1.73 m2 and urinary albumin excretion ≥300 mg/g creatinine), with moderate or severe DR, with PN, and with CAN. The two groups were matched as far as possible for age, sex, and degree of metabolic control (Table 1). Samples from these participants were used for the evaluation of serum miRNA expression profile.
Table 1. Characteristics of type 1 diabetes individuals evaluated in the micro-RNA expression profile.
The top six most differentially expressed miRNAs between the aforementioned groups were further validated in additional 47 type 1 diabetes individuals sorted according to the absence or presence of microvascular complications, using the same criteria described above: (1) n = 20 individuals without DKD, DR, PN and, CAN; and (2) n = 27 individuals with DKD, moderate or severe DR, PN, and CAN. The two groups were matched for age, sex, degree of metabolic control, diabetes duration, and age at diagnosis. The lower body mass index (BMI) and the higher frequency of arterial hypertension in the group with multiple microvascular complications reflect the presence of advanced DKD (Table 2).
Table 2. Characteristics of type 1 diabetes individuals evaluated in the validation of differentially expressed micro-RNAs.
The study was conducted in compliance with the Declaration of Helsinki, in accordance with the Institutional Ethics Committees. After signing informed consent, the participants were evaluated for clinical and biochemical characteristics and for the status of DKD, DR, PN, and CAN as previously described (10).
Serum miRNA Expression Profile
Serum miRNA profile was performed in four steps. Firstly, peripheral blood from the 10 type 1 diabetes individuals was collected in a yellow cap BD tube (Becton Dickinson, Franklin Lakes, USA) specific for serum separation. Tubes were centrifuged at 4°C for 20 min at 1,500 g and serum was separated, aliquoted, and frozen at −80°C. All 10 individuals' sera were analyzed in NanoDrop™ 2000 Spectrophotometer (NanoDrop, Rockland, USA) for the degree of hemolysis at 414 nm (absorbance peak of free hemoglobin) and all presented reading <0.2, being ready for the next step (15, 16), the extraction of total miRNA by the miRNeasy Mini commercial kit (Qiagen, Hilden, Germany).
Before the extraction, 3.5 μL of miRNeasy Serum/Plasma SpikeIn Control (Qiagen) were added to 200 μL of each serum sample, corresponding to 1.6 ×108 copies/μL of miR-39, which was used as an exogenous positive control.
The amount and quality of the extracted miRNAs were evaluated by spectrophotometer ND-1000. A quantitative real-time polymerase chain reaction (qRT-PCR) after reverse transcription was performed for miR-39 using Taqman® miR-39 assay (ThermoFisher Scientific, Carlsbad, USA) and Taqman® Master Mix Universal II at UNG (ThermoFisher Scientific) and confirmed the maintenance and integrity of the extracted miRNAs.
For synthesis and pre-amplification of complementary DNA (cDNA) from total miRNA, the Megaplex™ RT Primers, Human Pool A v2.1 (Thermo Fisher Scientific) and the Thermal cycler Veriti 96 WellsCycler (Thermo Fisher Scientific) were used according to the manufacturer's instructions.
Taqman® Human MicroRNA Array A kit (ThermoFisher Scientific) and Taqman® Universal Master Mix II at UNG (ThermoFisher Scientific) were used for evaluation of serum miRNAs profile. The Taqman® Human MicroRNA Array A kit is a qRT-PCR assay with hydrolysis probes based on the Taqman Low Density Arrays (TLDA) system containing Taqman probes sufficient for expression analysis of 377 miRNAs, in addition to four endogenous controls, totaling 381 miRNAs. The experiment was performed with Viia 7 (Thermo Fisher Scientific).
Validation of Differentially Expressed miRNAs
The top six most differentially expressed miRNAs in the serum profile were validated in serum samples from 47 type 1 diabetes individuals by qRT-PCR. Extraction and evaluation of the miRNAs were performed as described above. Synthesis and pre-amplification of cDNA was performed with TaqMan™ Advanced miRNA cDNA Synthesis kit (Thermo Fisher Scientific) and Thermocycler Veriti 96-Well Thermal Cycler (Thermo Fisher Scientific) according to the manufacturer's instructions.
Taqman® Fast Advanced Master Mix (Thermo Fisher Scientific) and TaqMan® Advanced miRNA Assays (Thermo Fisher Scientific) were used for real-time PCR reaction in StepOne Plus (Thermo Fisher Scientific).
miRNAs content was calculated as relative expression units (2−ΔCt) normalized by the exogenous cel-miR-39 spike-in control and log10 transformed for statistical analysis.
Results are expressed as median and interquartile range. Clinical and biochemical differences between groups were assessed by Pearson's chi-squared test and Wilcoxon/Kruskal-Wallis test. A volcano plot on the results of two-way ANOVA test comparing serum miRNA expression profile between individuals without vs. with complications was generated, with the X axis representing the fold change and the Y axis representing the ANOVA P-value. Cutoff lines were added at−1.5 to 1.5 on the X axis and 0.05 on the Y axis (Figure 1). For the six differentially expressed miRNAs validated, likelihood ratio tests were employed to test models with and without the following covariates: sex, age, diabetes duration, HbA1c, and medicines whose frequency of use was significantly different between the groups with and without complications: statins, angiotensin-converting enzyme inhibitors (ACEI) and amlodipine; a P < 0.05 was considered significant. Statistics were performed with the JMP software (SAS Institute Inc., Cary, USA).
Figure 1. Volcano plot on the results of two-way ANOVA test comparing serum miRNA expression profile between individuals without vs. with complications, with the X axis representing the fold change and the Y axis representing the P-value. The group without complications was used as reference; each dot represents one micro-RNA.
Serum miRNA Expression Profile
Of the 377 target miRNAs evaluated, 193 miRNAs were expressed in the serum of individuals with type 1 diabetes. A total of 21 miRNAs presented a significantly higher expression in the group of individuals with all microvascular complications vs. the group without complications. Table 3 shows these 21 miRNAs ranked in order of decreasing significance.
Table 3. Differentially expressed serum micro-RNAs between type 1 diabetes individuals without vs. with multiple microvascular complications.
Validation of Differentially Expressed miRNAs
Out of the 21 differentially expressed miRNAs, the following miRs were validated by qRT-PCR: 518-3p, 34a-5p, 126-5p, 425-5p, 618, and 139-5p. Logistic regression analyses showed that miRNA-518-3p and miRNA-618 were positively associated with multiple microvascular complications in the unadjusted model (Model 1), after adjustment for age and sex (Model 2), for age, sex, diabetes duration, and HbA1c (Model 3) and for sex, age, diabetes duration, HbA1c and use of statin, ACEI and amlodipine (Model 4) (Table 4).
In silico analyses in the miRDB dataset did not reveal common target-genes between these two miRNAs; there are 31 target-genes for miR-518d-3p and 225 target-genes for miR-618.
The present study identified two miRNAs upregulated in type 1 diabetes individuals with all chronic microvascular complications in relation to individuals with no microvascular complications, miR-518d-3p and miR-618, widening the spectrum of epigenetic changes in the setting of diabetes complications.
miR-518 was identified as a miRNA potentially targeting peroxisome proliferator-activated receptor alpha (PPARα) (17). A functional assay confirmed that miR-518d directly targets PPARα by binding to the 3'-UTR of its mRNA in a study that also suggested that increased miR-518d expression in placenta is implicated in gestational diabetes (18).
PPARα is a nuclear receptor that plays a critical role in lipid homeostasis and inflammation (18). It has been shown that PPARα is decreased in the retina of experimental models of diabetes and that the PPARα knockout worsened DR (19). On the other hand, PPARα has anti-apoptotic and antioxidant effects in the ischemic retina (20), and there are studies showing a lower incidence of DR in individuals exposed to fibrates, which are PPARα agonists (21–23).
Thus, we hypothesized that miR-518d-3p is involved in the PPARα downregulation reported in DR and that the upregulation of miR-518d-3p observed in the group presenting all microvascular complications may reflect moderate to severe DR of individuals included in this group. An additional finding that suggests the involvement of miR-518d-3p in DR is its increased expression in human retinal pigment cells exposed to oxidative stress (24). It is interesting to mention that another miRNA, miR-21, was found to be overexpressed in the retina of a type 2 diabetes rodent model and at least partially responsible for PPARα downregulation (25).
It is also possible that upregulation of miR-518d-3p be a marker of generalized endothelial dysfunction in these individuals since the expression of this miRNA increased in human endothelial cells exposed to acrolein, an unsaturated aldehyde able to induce oxidative stress and inflammation (26), conditions also triggered by hyperglycemia (27).
Regarding miRNA-618, its deregulation has been associated with several neoplasias (28). It has also been identified as a miRNA induced by hypoxia in human primary pulmonary artery smooth muscle cells (29) and upregulated in the blood of individuals during the acute phase of Kawasaki disease, a vasculitis in which transforming growth-factor β (TGFβ) is implicated (30). Among the predicted targets for miR-618 are the following seven genes belonging to TGFβ pathway, TGFB2, THBS1, NOG, KLF10, PPP2CA, PPP2CB, and ID2, according to miRDB database. Given the involvement of TGFβ in diabetes complications (31), further studies are necessary to elucidate whether miR-618 upregulation is related to deregulation of the TGFβ pathway in the setting of diabetes complications.
None of the 21 miRNAs upregulated in the group with all microvascular complications was reported to be associated with hyperglycemia in a previous study of circulating miRNAs in individuals with type 1 diabetes. In addition, the expression of the six validated microRNAs did not differ significantly when individuals were sorted by HbA1c tertiles (data not shown).
Three of the 21 miRNAs upregulated in the group with all microvascular complications (miR-139-5p, miR-106a, and miR-574-3p) were previously identified as upregulated in type 1 diabetes individuals with one or more complications (cardiovascular diseases and/or DR and/or albuminuria) in comparison to individuals with no complications (11), while four miRNAs upregulated in the present study were found downregulated in that study (miR-92a, miR-126, miR-483-5p, and miR-150) (11). These findings suggest that the serum miRNA profile may be sensitive to the presence or absence of specific chronic complications.
This study must be interpreted in light of limitations of cross-sectional studies and requires confirmation in larger cohorts. Among the study's strengths are the use of individual samples (as opposed to pooled serum samples) in the screening phase, the systematic evaluation of all microvascular complications (including CAN) and the homogeneity of the two evaluated groups.
In conclusion, in this cohort of type 1 diabetes individuals, serum miR-518d-3p and miR-618 are upregulated in those with diabetes kidney disease, diabetes retinopathy, peripheral neuropathy and cardiovascular autonomic neuropathy in comparison to individuals with no microvascular complications.
The study was conducted in compliance with the Declaration of Helsinki, in accordance with the Comitê de Ética e Pesquisa (CAPPesq) do Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo (#1.433.880). All participants provided written informed consent prior to sample collection.
DS-B collected and analyzed clinical and biochemical data, performed miRNA experiments, wrote, reviewed, and edited the manuscript. AS performed miRNA experiments and reviewed the manuscript. GG collected clinical and biochemical data, performed miRNA extraction, and reviewed the manuscript. SA performed neuropathy evaluation and reviewed the manuscript. RP and CM performed retinopathy evaluation and reviewed the manuscript. TP collected clinical and biochemical data and reviewed the manuscript. MP, UM, and MdS participated in study design and reviewed the manuscript. MQ collected clinical data and reviewed the manuscript. MC-G designed the study, wrote, reviewed and edited the manuscript.
This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) to MP, UM, and MC-G (2016/15603-0). MC-G, MP, and UM are recipients of a fellowship from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
3. Urbano F, Filippello A, Di Pino A, Barbagallo D, Di Mauro S, Pappalardo A, et al. Altered expression of uncoupling protein 2 in GLP-1-producing cells after chronic high glucose exposure: implications for the pathogenesis of diabetes mellitus. Am J Physiol Cell Physiol. (2016) 310:C558–67. doi: 10.1152/ajpcell.00148.2015
4. Prattichizzo F, Giuliani A, Ceka A, Rippo MR, Bonfigli AR, Testa R, et al. Epigenetic mechanisms of endothelial dysfunction in type 2 diabetes. Clin Epigenet. (2015) 7:56. doi: 10.1186/s13148-015-0090-4
7. Schneider A, Victoria B, Lopez YN, Suchorska W, Barczak W, Sobecka A, et al. Tissue and serum microRNA profile of oral squamous cell carcinoma patients. Sci Rep. (2018) 8:675. doi: 10.1038/s41598-017-18945-z
8. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci USA. (2011) 108:5003–8. doi: 10.1073/pnas.1019055108
10. Monteiro MB, Santos-Bezerra DP, Thieme K, Admoni SN, Perez RV, Machado CG, et al. Thioredoxin interacting protein expression in the urinary sediment associates with renal function decline in type 1 diabetes. Free Radic Res. (2016) 50:101–10. doi: 10.3109/10715762.2015.1109083
11. Barutta F, Bruno G, Matullo G, Chaturvedi N, Grimaldi S, Schalkwijk C, et al. MicroRNA-126 and micro-/macrovascular complications of type 1 diabetes in the EURODIAB Prospective Complications Study. Acta Diabetol. (2017) 54:133–9. doi: 10.1007/s00592-016-0915-4
14. Bijkerk R, Duijs JM, Khairoun M, Ter Horst CJ, van der Pol P, Mallat MJ, et al. Circulating microRNAs associate with diabetic nephropathy and systemic microvascular damage and normalize after simultaneous pancreas-kidney transplantation. Am J Transplant. (2015) 15:1081–90. doi: 10.1111/ajt.13072
16. Kirschner MB, Kao SC, Edelman JJ, Armstrong NJ, Vallely MP, van Zandwijk N, et al. Haemolysis during sample preparation alters microRNA content of plasma. PLoS ONE. (2011) 6:e24145. doi: 10.1371/journal.pone.0024145
18. Zhao C, Zhang T, Shi Z, Ding H, Ling X. MicroRNA-518d regulates PPARα protein expression in the placentas of females with gestational diabetes mellitus. Mol Med Rep. (2014) 9:2085–90. doi: 10.3892/mmr.2014.2058
19. Hu Y, Chen Y, Ding L, He X, Takahashi Y, Gao Y, et al. Pathogenic role of diabetes-induced PPAR-α down-regulation in microvascular dysfunction. Proc Natl Acad Sci USA. (2013) 110:15401–6. doi: 10.1073/pnas.1307211110
20. Moran E, Ding L, Wang Z, Cheng R, Chen Q, Moore R, et al. Protective and antioxidant effects of PPARα in the ischemic retina. Invest Ophthalmol Vis Sci. (2014) 55:4568–76. doi: 10.1167/iovs.13-13127
21. Keech AC, Mitchell P, Summanen PA, O'Day J, Davis TM, Moffitt MS, et al. Effect of fenofibrate on the need for laser treatment for diabetic retinopathy (FIELD study): a randomised controlled trial. Lancet. (2007) 370:1687–97. doi: 10.1016/S0140-6736(07)61607-9
22. Chew EY, Ambrosius WT, Davis MD, Danis RP, Gangaputra S, Greven CM, et al. Effects of medical therapies on retinopathy progression in type 2 diabetes. N Engl J Med. (2010) 363:233–44. doi: 10.1056/NEJMoa1001288
23. Morgan CL, Owens DR, Aubonnet P, Carr ES, Jenkins-Jones S, Poole CD, et al. Primary prevention of diabetic retinopathy with fibrates: a retrospective, matched cohort study. BMJ Open. (2013) 3:e004025. doi: 10.1136/bmjopen-2013-004025
25. Chen Q, Qiu F, Zhou K, Matlock HG, Takahashi Y, Rajala RVS, et al. Pathogenic role of microRNA-21 in diabetic retinopathy through downregulation of PPARα. Diabetes. (2017) 66:1671–82. doi: 10.2337/db16-1246
26. Lee SE, Yang H, Son GW, Park HR, Cho J-J, Ahn H-J, et al. Identification and characterization of MicroRNAs in acrolein-stimulated endothelial cells: implications for vascular disease. BioChip J. (2015) 9, 144–155. doi: 10.1007/s13206-015-9303-3
29. Wu C, So J, Davis-Dusenbery BN, Qi HH, Bloch DB, Shi Y, et al. Hypoxia potentiates microRNA-mediated gene silencing through posttranslational modification of Argonaute2. Mol Cell Biol. (2011) 31:4760–74. doi: 10.1128/MCB.05776-11
30. Shimizu C, Kim J, Stepanowsky P, Trinh C, Lau HD, Akers JC, et al. Differential expression of miR-145 in children with Kawasaki disease. PLoS ONE. (2013) 8:e58159. doi: 10.1371/journal.pone.0058159
Keywords: microvascular diabetes complications, epigenetics, micro-RNAs, diabetes kidney disease, diabetic retinopathy, peripheral neuropathy, cardiovascular autonomic neuropathy
Citation: Santos-Bezerra DP, Santos AS, Guimarães GC, Admoni SN, Perez RV, Machado CG, Pelaes TS, Passarelli M, Machado UF, Queiroz MS, da Silva MER and Correa-Giannella ML (2019) Micro-RNAs 518d-3p and 618 Are Upregulated in Individuals With Type 1 Diabetes With Multiple Microvascular Complications. Front. Endocrinol. 10:385. doi: 10.3389/fendo.2019.00385
Received: 29 October 2018; Accepted: 29 May 2019;
Published: 12 June 2019.
Edited by:Gaetano Santulli, Columbia University, United States
Reviewed by:Antonio Longo, University of Catania, Italy
Davide Barbagallo, University of Catania, Italy
Copyright © 2019 Santos-Bezerra, Santos, Guimarães, Admoni, Perez, Machado, Pelaes, Passarelli, Machado, Queiroz, da Silva and Correa-Giannella. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maria Lucia Correa-Giannella, firstname.lastname@example.org