Next Article in Journal
A Deep Learning Approach for Atrial Fibrillation Classification Using Multi-Feature Time Series Data from ECG and PPG
Previous Article in Journal
Detecting Diabetic Ketoacidosis with Infection: Combating a Life-Threatening Emergency with Practical Diagnostic Tools
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Euglycemia to Recent Onset of Type 2 Diabetes Mellitus: A Proof-of-Concept Study on Circulating microRNA Profiling Reveals Distinct, and Early microRNA Signatures

by
Marta Greco
1,
Maria Mirabelli
1,
Alessandro Salatino
1,
Francesca Accattato
1,
Vincenzo Aiello
2,
Francesco S. Brunetti
1,
Eusebio Chiefari
1,
Salvatore A. Pullano
1,
Antonino S. Fiorillo
1,
Daniela P. Foti
3,* and
Antonio Brunetti
1,*
1
Department of Health Sciences, University “Magna Græcia” of Catanzaro, 88100 Catanzaro, Italy
2
Department of Precision Medicine, Vanvitelli University, 80133 Naples, Italy
3
Department of Experimental and Clinical Medicine, University “Magna Græcia” of Catanzaro, 88100 Catanzaro, Italy
*
Authors to whom correspondence should be addressed.
Diagnostics 2023, 13(14), 2443; https://doi.org/10.3390/diagnostics13142443
Submission received: 13 June 2023 / Revised: 19 July 2023 / Accepted: 19 July 2023 / Published: 21 July 2023
(This article belongs to the Section Pathology and Molecular Diagnostics)

Abstract

:
Background and aim—Alterations in circulating microRNA (miRNA) expression patterns are thought to be involved in the early stages of prediabetes, as well as in the progression to overt type 2 diabetes mellitus (T2D) and its vascular complications. However, most research findings are conflicting, in part due to differences in miRNA extraction and normalization methods, and in part due to differences in the study populations and their selection. This cross-sectional study seeks to find new potentially useful biomarkers to predict and/or diagnose T2D by investigating the differential expression patterns of circulating miRNAs in the serum of patients with impaired fasting glucose (IFG) and new-onset T2D, with respect to euglycemic controls, using a high-throughput 384-well array and real-time PCR. Methods—Thirty subjects, aged 45–65 years, classified into three matched groups (of 10 participants each) according to their glycometabolic status, namely (1) healthy euglycemic controls, (2) patients with IFG and (3) patients with new-onset, uncomplicated T2D (<2 years since diagnosis) were enrolled. Circulating miRNAs were extracted from blood serum and profiled through real-time PCR on a commercial 384 well-array, containing spotted forward primers for 372 miRNAs. Data analysis was performed by using the online data analysis software GeneGlobe and normalized by the global Ct mean method. Results—Of the 372 analyzed miRNAs, 33 showed a considerably different expression in IFG and new-onset T2D compared to healthy euglycemic controls, with 2 of them down-regulated and 31 up-regulated. Stringent analysis conditions, using a differential fold regulation threshold ≥ 10, revealed that nine miRNAs (hsa-miR-3610, hsa-miR-3200-5p, hsa-miR-4651, hsa-miR-3135b, hsa-miR-1281, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p and hsa-let-7a-5p) showed a specific increase in new-onset T2D patients compared to IFG patients, suggesting their possible role as early biomarkers of progression from prediabetes to T2D. Moreover, by conventional fold regulation thresholds of ±2, hsa-miR-146a-5p was down-regulated and miR-1225-3p up-regulated in new-onset T2D patients only. Whereas hsa-miR-146a-5p has a well-known role in glucose metabolism, insulin resistance and T2D complications, no association between hsa-miR-1225-3p and T2D has been previously reported. Bioinformatic and computational analysis predict a role of hsa-miR-1225-3p in the pathogenesis of T2D through the interaction with MAP3K1 and HMGA1. Conclusions—The outcomes of this study could aid in the identification and characterization of circulating miRNAs as potential novel biomarkers for the early diagnosis of T2D and may serve as a proof-of-concept for future mechanistic investigations.

1. Introduction

MicroRNAs (miRNAs) are a class of short noncoding single-strand RNAs that are generated from genomic DNA. They are characterized by tissue-specific expression patterns and are involved in the regulation of fundamental biological processes, such as tissue development, cell proliferation and differentiation, metabolism, apoptosis and aging, displaying a role in the pathophysiology of diseases [1,2,3,4,5,6,7]. From a biofunctional perspective, miRNAs typically exert a negative effect on gene expression by targeting sites in the 3′ untranslated region (3′UTR) of specific mRNAs, resulting in destabilization, translational repression or mRNA cleavage. However, miRNAs up-regulating gene expression have also been reported [8,9]. To date, more than 1000 miRNAs have been described in humans [10,11], with many of them identifiable in the blood stream [12,13,14,15]. As single-strand RNA molecules, miRNAs are susceptible to degradation by extracellular RNAses, unless they are protected by the encapsulation into lipid-vesicles (i.e., exosomes, microvesicles and apoptotic bodies) [15]. This process preserves the stability of circulating miRNAs, allowing these molecules to act as mediators of intercellular communication [16]. Since their initial discovery, circulating miRNAs have shown specific signatures in response to pathological stimuli, with the potential to be used as biomarkers for diagnosing and/or monitoring diseases [16,17,18]. The interest in circulating miRNAs has increased intensively in recent years, while efficient methods for recovering miRNAs from biological fluids and for quantitative detection, such as real-time PCR and microarray, have become commercially available [15]. With the aid of technical improvements in miRNA profiling, new miRNA molecules have been associated with the dysregulation of many human diseases, including type 2 diabetes (T2D) [16,17,18,19]. T2D can evolve silently over a long period of time, from normal glucose metabolism to impaired fasting glucose (IFG) or prediabetes, and then to overt hyperglycemia and microvascular and macrovascular complications, therefore the scientific community has now shifted its focus on finding strategies for early diagnosis and prevention of T2D, rather than disease control [20,21,22]. Intervening early and specifically with patients who are in the early stages of the disease can help to reduce the high prevalence of T2D and associated risks for adverse clinical outcomes such as blindness, kidney failure, cardiovascular events, stroke and lower limb amputation. This, in turn, would considerably reduce the overall impact of T2D on healthcare [22]. Even though it is described that tissues targeted by insulin release miRNAs into the extracellular space [23,24], and these are involved in the onset of microvascular complications of T2D [15,19,25], there are only few studies investigating the circulating miRNAs signature of the initial stages of disease, when the alterations in glucose metabolism are still mild and clinically asymptomatic [26,27]. Here, we evaluated the expression of circulating miRNAs in age- and sex-matched patients with or without IFG and recent onset of T2D, in order to determine the differences between groups and identify miRNA signatures specifically associated with disease risk. To ensure the accuracy and reliability of our high-throughput miRNA profiling data analysis, we implemented strict criteria for patient enrollment, took careful measures to eliminate any preanalytical biases that could affect miRNA extraction and purification (such as preventing contamination from intracellular miRNAs resulting from hemolysis), and applied a stringent mean-centric method to normalize miRNA expression. This strategy could serve in the identification and characterization of potential novel biomarkers for the early diagnosis of T2D and may represent a proof-of-concept for future mechanistic investigations.

2. Materials and Methods

2.1. Study Design

For this proof-of-concept, cross-sectional study aimed at discovering novel circulating biomarkers in T2D, 30 unrelated, well-selected, individuals, in the age range of 45–65 years, were enrolled at the Endocrinology Unit “Renato Dulbecco” University Hospital of Catanzaro, Italy. The study population was divided into three age- and sex-matched groups, of 10 patients each, according to their glycometabolic status: (1) euglycemic individuals (control group); (2) patients with impaired fasting glucose (IFG group); (3) patients with a recent (<2 years) diagnosis of T2D (new-onset T2D group), without any micro and macrovascular complications. Diagnosis of IFG and T2D was performed based on 2019 American Diabetes Association Criteria [28]. Inflammation (C-Reactive Protein (CRP) > 6 mg/L; Erythrocyte Sedimentation rate (ESR) > 30 mm/h), obesity (BMI > 30 kg/m2), poor glycemic control (glycated hemoglobin, HbA1c > 8%, 64 mmol/mol), smoking, alcohol and drug use, including antidiabetic pharmacological medications with possible impact on miRNAs (i.e., metformin [29]), as well as intercurrent illnesses and moderate-to-severe hypertriglyceridemia (triglycerides > 175 mg/dL, 2.0 mmol/L) were assumed as exclusion criteria. Additionally, in consideration of the possible alterations in miRNA expression throughout the various stages of chronic kidney disease [18,30], to mitigate any potential bias in the detection of circulating miRNAs that are specifically linked to the progression from prediabetes to T2D, individuals with indications of renal impairment (such as an estimated glomerular filtration rate < 60 mL/min/1.73 m2 or the presence of microalbuminuria), regardless of whether it was related to diabetic nephropathy [31], were excluded from the study.

2.2. Anthropometric and Laboratory Determinations

After 12 h overnight fast, all patients were checked for blood pressure and anthropometric determinations (body weight and height). Venous blood was collected to measure biochemical metabolic parameters, such as fasting glucose, insulin, lipid profile (total, HDL- and LDL-cholesterol and triglycerides) and HbA1c. ESR and CRP were determined to exclude ongoing acute inflammatory processes. In addition, for miRNA testing, venous blood (~10 mL) was collected from each patient in a serum separator tube (Vacutainer R SSTTM, BD, San Jose, CA, USA) and centrifuged within 1 h at 2000× g for 15 min at 4 °C. Serum fractions were rapidly aliquoted into cryovial tubes and stored at −80 °C for subsequent miRNA determinations. Measurements of fasting glucose and lipid profile were performed on Cobas 6000 (Roche, Basel, Switzerland). Serum CRP concentrations were determined by immunonephelometric method (Siemens, Alpharetta, GA, USA) using BNII nephelometer (Siemens, Alpharetta, GA, USA). ESR was determined by capillary photometry technology using TEST 1 analyzer (ALIFAX, Padova, Italy). HbA1c was measured by high-performance liquid chromatography using Premier Hb9210™ HbA1c Analyzer (A. Menarini, Firenze, Italy). Serum insulin concentration was measured by the ADVIA Centaur Immunoassay system (Siemens, Alpharetta, GA, USA), using a chemiluminescent immunoassay. All the aforementioned tests were carried out in accordance with the manufacturer’s instructions. Quality control was assessed on a daily basis for all determinations, and the results met the precision targets indicated by the manufacturers.

2.3. Circulating miRNAs Profiling

Stored serum aliquots were thawed on ice and centrifuged at 16,000× g for 5 min at 4 °C to remove insoluble material. To exclude hemolysis, all samples were subjected to spectrophotometric measurement at 414 nm before continuing with further processing, and, as confirmation test, delta Ct difference between miR-23a-3p and miR-451a had to be less than 5. Circulating miRNAs were then isolated from serum samples (200 µL) of individual patients (n = 10 per group) as enriched low-molecular-weight RNA fractions (LMW RNA) using a commercial column-based system (miRNeasy Serum/Plasma Kit®, QIAGEN, Hilden, Germany), without specifically isolating miRNA-containing lipid-vesicles. As recommended by the manufacturer, synthetic C. elegans miR-39 was used as exogenous control and spiked-in to monitor LMW RNA extraction efficiency. The quantity and purity of the extracted LMW RNA (A260/280 ratio = 1.9–2.1 and A260/230 ratio = 2.0–2.2) were determined by a spectrometric method (Nanodrop 1000, Thermo Scientific, Waltham, MA, USA). Afterwards, equal amounts (100 ng) of the individual LMW RNA fractions were pooled for each study group. The pooled LMW RNA fraction (500 ng) was used for a subsequent one-step tailed universal reverse transcription reaction (reaction volume of 100 µL), carried out by miScript II RT Kit (Qiagen) using a thermocycler system (GeneAmp 2700 Thermal Cycler, Applied Biosystems, Waltham, MA, USA). The cDNA products obtained from reverse transcription were used as templates for miRNA profiling by real-time PCR on a high-throughput 384 well-array, containing spotted forward primers for 372 miRNAs with known detectable expression in the human serum (miScript miRNA PCR Array Human Serum & Plasma 384HC®, QIAGEN [19] and SYBR Green chemistry (SYBR Green PCR Kit, QIAGEN). Reactions were carried out on QuantStudio™ 12K Flex Real-Time PCR System (Applied Biosystems). miRNAs expression analysis was conducted by a free online software (https://geneglobe.qiagen.com, accessed on 1 February 2023), using 6 reference endogenous controls (small nucleolar RNA: SNORD61, SNORD68, SNORD72, SNORD95, SNORD96A and RNU6B/RNU6-2) and spike-in C. elegans miR-39 as normalizers. MiRNAs whose amplification had reported raw Ct values > 35 were excluded.

2.4. Global Ct Mean Normalization and Bioinformatic Analysis

For this large-scale miRNA profiling analysis and biomarker identification, the free online data analysis software GeneGlobe was employed (https://geneglobe.qiagen.com, accessed on 1 February 2023). The Global Ct Mean, which uses the calculated mean expression level of all miRNAs in a given sample, including the reference endogenous controls and spike-in miR-39, was automatically implemented as a strongly reliable method of miRNA expression normalization. This normalization strategy assumes that the mean expression level of all miRNAs in a sample, from either the healthy control group or the patient group (IFG or new-onset T2D), is constant when using the same total RNA input [32,33]. By convention, miRNAs with a differential fold regulation > 2 or <−2 between the two groups were considered significant in GeneGlobe analysis, as previously described [34].
To predict the target genes and metabolic phenotypes associated with aberrant circulating miRNAs, a combination of bioinformatic analysis via internet databases and computational tools, including miRbase (http://www.mirbase.org, accessed on 1 February 2023), mirDB (http://Mirdb.org, accessed on 1 February 2023), mirDIP 4.1-integrative database (http://ophid.utoronto.ca/mirDIP, accessed on 1 February 2023), TargetScan (http://www.targetscan.org/vert_71, accessed on 1 February 2023) was used.

2.5. Statistical Analysis

Data were expressed as median and interquartile ranges. The Kruskal–Wallis non-parametric test was used to determine intergroup variations in clinical and biochemical parameters. A p value less than 0.05 was considered statistically significant. Data analysis was carried out using SPSS, software version 20.0 (SPSS Inc., Chicago, IL, USA).

3. Results

3.1. Characteristics of the Study Popupulation

The biochemical and anthropometric characteristics of patients enrolled in this study are summarized in Table 1. The three groups under study, namely control, IFG and new-onset T2D, exhibited similar distribution according to gender and age. Additionally, no variations were noted regarding BMI, lipid profile, liver function tests and inflammatory parameters. The sole discrepancy was observed in fasting glucose and HbA1c levels, which accurately reflected the gradual decline in glucose tolerance characteristic of the transition from euglycemia to prediabetes and T2D. This finding confirmed that the selection of cases (IFG and new-onset T2D patients) and controls was appropriate for the purpose of the study. Although diastolic blood pressure levels exhibited a slight, yet statistically significant difference, it did not have clinical implications as none of the patients were considered hypertensive or needed medications to manage their blood pressure (Table 1).

3.2. Circulating miRNAs Profiling by Real-Time PCR

A high-throughput 384 well array profiling of circulating miRNAs was performed on pooled serum LMW RNA fractions, obtained from the three studies and reversely transcribed into cDNAs. In the control group, 3 miRNAs were found to be highly expressed (Ct ≤ 25) and 35 miRNAs were found to be moderately expressed (Ct between 25 and 30). In the IFG group, 4 miRNAs were found to be highly expressed (Ct ≤ 25) and 78 miRNAs were found to be moderately expressed (Ct ≤ 30). In the new-onset T2D group, 7 miRNAs were found to be highly expressed (Ct ≤ 25) and 88 miRNAs were found to be moderately expressed (Ct ≤ 30), therefore suggesting that a progressive change of expression of circulating miRNA could occur in the transition from euglycemia to recent onset of T2D. The supplementary material (Table S1) provides the mature ID of the highly expressed miRNAs for each study group. Out of all miRNAs detectable by real-time PCR (Ct ≤ 35), 95 miRNAs were found to be over-expressed and 54 under-expressed in the IFG group vs. the control group; 109 miRNAs were found to be over-expressed and 52 under-expressed in the new-onset T2D vs. the control group. When comparing new-onset T2D and IFG patients, 49 miRNAs were found to be over-expressed and 33 under-expressed. The scatter plot analysis of differential miRNA expression is illustrated in Figure 1. The mature ID of the differentially expressed miRNAs and their relative quantification (fold regulation) are detailed in Table 2 (for over-expressed miRNAs) and Table 3 (for under-expressed miRNAs).
A total of 33 circulating miRNAs showed considerably different expression levels between patients with IFG or new-onset T2D and matched euglycemic controls (Figure 2). This was determined by using a stringent analysis method (global Ct mean) and considering only those miRNAs whose expression level was relatively low (fold regulation < 10) in either the control or the test group and is reasonably high in the other group (fold regulation ≥ 10). Some miRNAs, such as hsa-miR-1260a, hsa-miR-4687-5p, hsa-miR-3135b, hsa-let-7b-5p, hsa-let-7a-5p, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p, hsa-miR-26b-5p were differentially expressed in both IFG and new-onset T2D compared to euglycemic control subjects (Figure 2). Notably, 9 miRNAs (hsa-miR-3610, hsa-miR-3200-5p, hsa-miR-4651, hsa-miR-3135b, hsa-miR-1281, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p, hsa-let-7a-5p) displayed a significant increase in fold regulation values, calculated using the euglycemic control subjects as the reference group, specifically in patients with new-onset T2D compared to those with IFG (Figure 3). This finding suggests their potential as early biomarkers for disease progression from prediabetes to T2D. Moreover, from a comparative analysis of relative fold regulation data, using the conventional threshold > 2 for up-regulation and <−2 for down-regulation automatically implemented by the Geneglobe software, it was found that hsa-miR-146a-5p was considerably down-regulated while hsa-miR-1225-3p was up-regulated uniquely in individuals with newly diagnosed T2D in comparison to euglycemic controls (fold regulation in new-onset T2D vs. controls, hsa-miR-1225-3p: 3.6; hsa-miR-146a-5p: −2.6). These changes in expression levels were not observed in individuals with IFG, despite a similar trend in direction (fold regulation in IFG vs. controls, hsa-miR-1225-3p: 1.9; hsa-miR-146a-5p: −1.3). Previous research has reported that hsa-miR-146a-5p plays a role in glucose metabolism, insulin resistance and development of T2D complications [35,36,37], whereas, to our knowledge, no association between hsa-miR-1225-3p and T2D has been formerly described.

4. Discussion

Several studies have suggested that circulating miRNAs could be useful for identifying various diseases and predicting clinical outcomes. However, there is no consensus on a miRNA signature that would enable early identification of T2D in its preclinical, asymptomatic stage and prediction of the development of its related complications [15,23].
This proof-of-concept study sought to identify differentially expressed miRNAs in serum of IFG and new-onset T2D patients compared to a control population of healthy euglycemic subjects, in order to evaluate a circulating miRNA signature of the initial onset of T2D. From comparative analysis of high-throughput real-time PCR profiling data, obtained from the total pool of reversely transcribed serum LMW RNA fractions, nine miRNAs, namely hsa-miR-1260a, hsa-miR-4687-5p, hsa-miR-3135b, hsa-let-7b-5p, hsa-let-7a-5p, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p and hsa-miR-26b-5p, were found over-expressed in both IFG and new-onset T2D patients compared to age- and sex-matched healthy controls. These miRNAs could be used to non-invasively identify patients in the very early stages of the natural history of T2D and strengthen the need for insulin-sensitizing approaches to revert the condition before any complication develops.
Further comparisons revealed that hsa-miR-3610, hsa-miR-3200-5p, hsa-miR-4651, hsa-miR-3135b, hsa-miR-1281, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p and hsa-let-7a-5p undergo a significant, progressive expression change from euglycemic individuals to patients with IFG and new-onset T2D. This suggests that this circulating miRNA signature can be used as biomarker for progression of disease in those with prediabetes and/or initial laboratory evidence of altered glucose control.
Some of these circulating miRNAs have been shown to be involved in the development of complications of T2D, such as hsa-miR-1281, which has been demonstrated by our group, through functional in vitro assays, to have a pathogenic role in diabetic retinopathy, the first microvascular alteration in the natural history of T2D [19]. By targeting HIF1AN mRNA, and de-repressing the hypoxia-inducible factor HIF-1a, hsa-miR-1281 positively regulates the expression of neoangiogenic VEGFA in endothelial cells exposed to a glucose excess environment, therefore contributing to retinal vascular damage [19].
Circulating levels of hsa-miR-146a-5p and hsa-miR-1225-3p were found, respectively, to be lower and higher, only in the group of new-onset T2D patients. Hsa-miR-146a-5p is expressed in white blood cells, and endothelial cells, while it is known to down-regulate the expression of NF-kB and inhibit the NF-kB-induced inflammatory signaling pathways [35,36,37]. A decrease in circulating hsa-miR-146a-5p may play an important role in the onset and development of T2D complications through the overexpression of NF-kB and modulation of systemic inflammatory responses that cause glucose control dysfunction [36,37]. In patients with T2D, changes in the blood expression levels of hsa-miR-146a-5p have been reported; however, there is conflicting information about whether these changes implicate an increase or decrease [38]. To date, there is no report on hsa-miR-1225-3p and its involvement in T2D. However, bioinformatic and computational analysis suggest that this miRNA may be involved in the pathogenesis and complications of T2D through its predicted targets, the mitogen activated protein kinase MAP3K1 and HMGA1, whose transcripts and 3′UTR regions have putative binding sites for hsa-miR-1225-3p.
MAP3K1 activates the ERK and JNK kinase cascade in response to various stress stimuli (i.e., cytokines, free fatty acids and hyperglycemia), and appears a crucial mediator in the transition from obesity to T2D [39]. The activation of JNK in liver, adipose tissue and skeletal muscle has been shown to interfere with the transmission of insulin signaling, resulting in insulin resistance. Furthermore, in pancreatic beta cells, this activation results in decreased compensatory insulin secretion [39]. In addition, MAP3K1 activates CHUK and IKBKB, the central protein kinases of the NF-kB pathway [40,41]. MiRNA hsa-miR-1225-3p is also predicted to target HMGA1, an architectural transcription factor involved in regulating the expression of several glucose-responsive genes. HMGA1 is known to modulate fundamental biological processes, such as cell proliferation, migration and cell metabolism [29,42,43]. Defective variants of the HMGA1 gene have been linked to reduced insulin receptor expression and increased susceptibility to T2D [44,45,46]. It is tempting for us to speculate that, by binding to the 3′UTR of the HMGA1 mRNA and facilitating the decay of transcript, increased levels of hsa-miR-1225-3p could produce a down-regulation of the insulin receptor expression and insulin resistance.
The bioinformatic analysis of human miRNA target predictions offers potential insights into how other miRNAs that are differentially expressed in our study might interfere with the insulin signaling pathway. For example, hsa-miR-1260a might play a role in regulating glucose metabolism and insulin sensitivity, by targeting PPARGC1A and PTPN1. Decreased expression of PPARG coactivator 1 alpha, encoded by the PPARGC1A gene, has been associated with impaired insulin function in skeletal myocytes [47,48,49,50], whereas the tyrosine phosphatase, encoded by the PTPN1 gene, has been evidenced to act as a negative regulator of the insulin signaling pathway, by dephosphorylating multiple phospho-thyrosine residues within the activated insulin receptor [51,52,53,54]. Moreover, the let-7 family of miRNA has been shown to be involved in the regulation of glucose metabolism and insulin sensitivity by acting on targets associated with the insulin/IGF-1 receptor pathway [55].
MiRNA hsa-miR-4687-5p, instead, might be involved in the regulation of TNF-alpha induced NF-kappa-B activation and apoptosis by targeting ZFAND6 (AN1-type zinc finger protein 6). Persistent hyperglycemia leads to the activation of the NF-kB pathway, which has a key role in the pathogenesis of the micro- and macrovascular complications of T2D. This is caused by an increased expression of different pro-inflammatory cytokines, including TNF-alpha [56,57,58]. TNF-alpha is one of the earliest markers of obesity and T2D, and its raised levels represent a major contributing factor to the development and maintenance of insulin resistance [59,60].
Figure 4 shows the hypothetic mechanisms emerged from our findings, by which circulating deregulated miRNAs cause insulin resistance, inflammation and vascular damage, both leading to susceptibility to T2D development.
Although there is a growing number of individual miRNA studies focused on the prediction of T2D and its vascular complications, it is important to acknowledge some limitations commonly found in these types of studies. Specifically, individual miRNA-disease association studies often suffer from insufficient statistical power due to small sample sizes, are susceptible to many kinds of experimental biases, and may lead to inconsistent results. Very recently, meta-analytic approaches operating on public high-throughput transcriptomic datasets in the NCBI Gene Expression Omnibus (GEO) repository have identified hsa-miR-1260a and hsa-miR-146a-5p as two out of nine potential circulating miRNA markers for incident T2D [61]. The findings we have observed in our study are partially congruent with the small signature of miRNAs that has been derived from the aforementioned meta-analysis [61]. However, meta-analytical investigations on circulating miRNAs are unavoidably influenced by the diverse characteristics of available high-throughput miRNA profiling datasets and the unequal quantities of dysregulated miRNAs in T2D patients compared to healthy subjects across different studies. Furthermore, this approach excludes individuals with IFG [61]. A further concern is the considerable diversity in algorithms employed for normalizing miRNAs in the context of disease association studies, encompassing methods such as spike-in normalization, small nucleolar RNA normalization and mean-centric normalization. In light of this, the Global Ct mean normalization has recently undergone validation as the most effective normalization method in high-throughput miRNA profiling studies where the process of miRNA biomarker identification is initiated [31]. The use of diverse normalization techniques in high-throughput miRNA profiling studies may explain the lack of prior reports indicating an association between T2D and the increased expression of circulating hsa-miR-1225-3p [62]. With the broader implementation of the Global Ct mean normalization method in high-throughput miRNA profiling studies, it is conceivable that greater reproducibility will be achieved, particularly in the systematic identification of circulating miRNA markers that are candidate for further validation processes. In this regard, it is noteworthy to mention that in our study, the differential expression of circulating miRNAs and the direction of variation remained consistent even when alternative normalization methods, such as those based on small nucleolar RNAs, were used. Nevertheless, there was a discrepancy in the magnitude of fold regulation observed, as normalizing with the Global Ct mean yielded lower values compared to less stringent normalization methods.
Although the mechanisms underlying changes in circulating miRNA levels in T2D are not fully understood, previous studies have suggested that miRNAs can be actively or passively released by a variety of cells and, by the association with microvesicles or protein complexes, transferred in an active form to recipient cells [15]. Based on the present findings, although the source and target cells of extracellular miRNAs remain uncertain, as well as their specific carriers [63], we can hypothesize that modifications in specific circulating miRNAs, potentially with a pathogenic role, may occur simultaneously during the course of disease progression from normal glucose tolerance to IFG and overt T2D. While the main strength of this work is the very strict selection of the study population, in which both clinical confounders and pharmacological therapies were excluded, and the technical rigor in the procedures applied to blood samples, the major limitation of the study is the lack of individual validation of the identified circulating miRNA signature in a separate cohort of patients with new-onset T2D and prediabetes. The design of this work, which involves the high-throughput profiling of pooled LMW RNA fractions from participants, limits the potential for conducting correlation analysis between circulating miRNAs and conventional clinical and laboratory parameters. Additionally, it hinders the identification of the precise miRNA-regulated network responsible for these disease traits. However, this study was designed as a proof-of-concept study that can now encourage new research on selected miRNAs, including the novel hsa-miR-1225-3p, in the pathogenesis and therapeutics of T2D.

5. Conclusions

Despite intensive research and rapid progress in recent decades, the key mechanisms underlying T2D pathogenesis, and the development of complications are still not fully elucidated. Circulating miRNAs can be used as noninvasive biomarkers for the identification and monitoring of various diseases, including T2D. High-throughput profiling data analysis revealed that some miRNAs could be considered as noninvasive diagnostic biomarkers of the recent onset of T2D (hsa-miR-1260a, hsa-miR-4687-5p, hsa-miR-3135b, hsa-let-7b-5p and hsa-let-7a-5p, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p and hsa-miR-26b-5p) while others (hsa-miR-3610, hsa-miR-3200-5p and hsa-miR-4651, hsa-miR-3135b, hsa-miR-1281, hsa-miR-4301, hsa-miR-195-5p, hsa-miR-523-5p and hsa-let-7a-5p) as biomarkers of disease progression from prediabetes to T2D. For the first time in this study, it is suggested that the up-regulation of hsa-miR-1225-3p may be responsible for insulin resistance in T2D patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/diagnostics13142443/s1, Table S1: Highly expressed miRNAs in euglycemic controls, and IFG and new-onset T2D patients.

Author Contributions

Conceptualization, D.P.F. and A.B.; investigation, M.G., A.S. and F.A.; data curation, M.G., M.M., V.A., F.S.B. and E.C.; funding acquisition, S.A.P. and A.S.F.; writing—original draft preparation, M.G. and M.M.; writing—review and editing, D.P.F. and A.B.; project administration, D.P.F. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Regione Calabria, Sez. Area Centro (protocol code no. 116, date of approval 14 May 2015).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data supporting the reported results are available from the corresponding author upon reasonable request.

Acknowledgments

Abstract of this study has been accepted for presentation at the Association for Molecular Pathology (AMP) Europe Congress 2023 (Milan, 18–20 June 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, Y.; Samal, E.; Srivastava, D. Serum response factor regulates a muscle-specific microRNA that targets Hand2 during cardiogenesis. Nature 2005, 436, 214–220. [Google Scholar] [CrossRef] [PubMed]
  2. Chen, J.F.; Mandel, E.M.; Thomson, J.M.; Wu, Q.; Callis, T.E.; Hammond, S.M.; Conlon, F.L.; Wang, D.Z. The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat. Genet. 2006, 38, 228–233. [Google Scholar] [CrossRef] [PubMed]
  3. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  4. Inukai, S.; Slack, F. MicroRNAs and the genetic network in aging. J. Mol. Biol. 2013, 425, 3601–3608. [Google Scholar] [CrossRef] [Green Version]
  5. Garzon, R.; Croce, C.M. MicroRNAs in normal and malignant hematopoiesis. Curr. Opin. Hematol. 2008, 15, 352–358. [Google Scholar] [CrossRef]
  6. Latronico, M.V.; Condorelli, G. MicroRNAs and cardiac pathology. Nat. Rev. Cardiol. 2009, 6, 419–429. [Google Scholar] [CrossRef]
  7. Dimmeler, S.; Nicotera, P. MicroRNAs in age-related diseases. EMBO Mol. Med. 2013, 5, 180–190. [Google Scholar] [CrossRef] [PubMed]
  8. Vasudevan, S.; Steitz, J.A. AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell 2007, 128, 1105–1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Cordes, K.R.; Sheehy, N.T.; White, M.P.; Berry, E.C.; Morton, S.U.; Muth, A.N.; Lee, T.-H.; Miano, J.M.; Ivey, K.N.; Srivastava, D. miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature 2009, 460, 705–710. [Google Scholar] [CrossRef] [Green Version]
  10. Strachan, T.; Read, A. Human Molecular Genetics; Chapter 9; Organization of the Human Genome: Geneva, Switzerland, 2010. [Google Scholar]
  11. Friedländer, M.R.; Lizano, E.; Houben, A.J.; Bezdan, D.; Báñez-Coronel, M.; Kudla, G.; Mateu-Huertas, E.; Kagerbauer, B.; González, J.; Chen, K.C.; et al. Evidence for the biogenesis of more than 1000 novel human microRNAs. Genome Biol. 2014, 15, R57. [Google Scholar] [CrossRef] [Green Version]
  12. Cortez, M.A.; Bueso-Ramos, C.; Ferdin, J.; Lopez-Berestein, G.; Sood, A.K.; Calin, G.A. MicroRNAs in body fluids—the mix of hormones and biomarkers. Nat. Rev. Clin. Oncol. 2011, 8, 467–477. [Google Scholar] [CrossRef] [Green Version]
  13. Mitchell, P.S.; Parkin, R.K.; Kroh, E.M.; Fritz, B.R.; Wyman, S.K.; Pogosova-Agadjanyan, E.L.; Peterson, A.; Noteboom, J.; O’Briant, K.C.; Allen, A.; et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc. Natl. Acad. Sci. USA 2008, 105, 10513–10518. [Google Scholar] [CrossRef]
  14. Weber, J.A.; Baxter, D.H.; Zhang, S.; Huang, D.Y.; Huang, K.H.; Lee, M.J.; Galas, D.J.; Wang, K. The microRNA spectrum in 12 body fluids. Clin. Chem. 2010, 56, 1733–1741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Greco, M.; Chiefari, E.; Mirabelli, M.; Salatino, A.; Pullano, S.A.; Aversa, A.; Foti, D.P.; Brunetti, A. Insights into the World of MicroRNAs. In Biomarkers in Diabetes. Biomarkers in Disease: Methods, Discoveries and Applications; Patel, V.B., Preedy, V.R., Eds.; Springer: Cham, Switzerland, 2023. [Google Scholar] [CrossRef]
  16. Chen, X.; Liang, H.; Zhang, J.; Zen, K.; Zhang, C.Y. Secreted microRNAs: A new form of intercellular communication. Trends Cell Biol. 2012, 22, 125–132. [Google Scholar] [CrossRef] [PubMed]
  17. Guay, C.; Regazzi, R. Circulating microRNAs as novel biomarkers for diabetes mellitus. Nat. Rev. Endocrinol. 2013, 9, 513–521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Bolignano, D.; Greco, M.; Presta, P.; Duni, A.; Vita, C.; Pappas, E.; Mirabelli, M.; Lakkas, L.; Naka, K.K.; Brunetti, A.; et al. A small circulating miRNAs signature predicts mortality and adverse cardiovascular outcomes in chronic hemodialysis patients. Clin. Kidney J. 2023, 16, 868–878. [Google Scholar] [CrossRef]
  19. Greco, M.; Chiefari, E.; Accattato, F.; Corigliano, D.M.; Arcidiacono, B.; Mirabelli, M.; Liguori, R.; Brunetti, F.S.; Pullano, S.A.; Scorcia, V.; et al. MicroRNA-1281 as a Novel Circulating Biomarker in Patients with Diabetic Retinopathy. Front. Endocrinol. 2020, 11, 528. [Google Scholar] [CrossRef]
  20. Brannick, B.; Wynn, A.; Dagogo-Jack, S. Prediabetes as a toxic environment for the initiation of microvascular and macrovascular complications. Exp. Biol. Med. 2016, 241, 1323–1331. [Google Scholar] [CrossRef] [Green Version]
  21. US Preventive Services Task Force. Screening for Prediabetes and Type 2 Diabetes: US Preventive Services Task Force Recommendation Statement. JAMA 2021, 326, 736–743. [Google Scholar] [CrossRef]
  22. Ortiz-Martínez, M.; González-González, M.; Martagón, A.J.; Hlavinka, V.; Willson, R.C.; Rito-Palomares, M. Recent Developments in Biomarkers for Diagnosis and Screening of Type 2 Diabetes Mellitus. Curr. Diabetes Rep. 2022, 22, 95–115. [Google Scholar] [CrossRef]
  23. Nigi, L.; Grieco, G.E.; Ventriglia, G.; Brusco, N.; Mancarella, F.; Formichi, C.; Dotta, F.; Sebastiani, G. MicroRNAs as Regulators of Insulin Signaling: Research Updates and Potential Therapeutic Perspectives in Type 2 Diabetes. Int. J. Mol. Sci. 2018, 19, 3705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Özcan, S. microRNAs in Pancreatic β-Cell Physiology. Adv. Exp. Med. Biol. 2015, 887, 101–117. [Google Scholar] [CrossRef] [PubMed]
  25. Ismail, A.; El-Mahdy, H.A.; Gomaa Eldeib, M.; Doghish, A.S. miRNAs as cornerstones in diabetic microvascular complications. Mol. Genet. Metab. 2023, 138, 106978. [Google Scholar] [CrossRef]
  26. Párrizas, M.; Brugnara, L.; Esteban, Y.; González-Franquesa, A.; Canivell, S.; Murillo, S.; Gordillo-Bastidas, E.; Cussó, R.; Cadefau, J.A.; García-Roves, P.M.; et al. Circulating miR-192 and miR-193b are markers of prediabetes and are modulated by an exercise intervention. J. Clin. Endocrinol. Metab. 2015, 100, E407-15. [Google Scholar] [CrossRef] [Green Version]
  27. Sidorkiewicz, I.; Niemira, M.; Maliszewska, K.; Erol, A.; Bielska, A.; Szalkowska, A.; Adamska-Patruno, E.; Szczerbinski, L.; Gorska, M.; Kretowski, A. Circulating miRNAs as a Predictive Biomarker of the Progression from Prediabetes to Diabetes: Outcomes of a 5-Year Prospective Observational Study. J. Clin. Med. 2020, 9, 2184. [Google Scholar] [CrossRef] [PubMed]
  28. American Diabetes Association Professional Practice Committee. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes-2022. Diabetes Care 2022, 45 (Suppl. S1), S17–S38. [Google Scholar] [CrossRef]
  29. Salatino, A.; Mirabelli, M.; Chiefari, E.; Greco, M.; Di Vito, A.; Bonapace, G.; Brunetti, F.S.; Crocerossa, F.; Epstein, A.L.; Foti, D.P.; et al. The anticancer effects of Metformin in the male germ tumor SEM-1 cell line are mediated by HMGA1. Front. Endocrinol. 2022, 13, 1051988. [Google Scholar] [CrossRef] [PubMed]
  30. Peters, L.J.F.; Floege, J.; Biessen, E.A.L.; Jankowski, J.; van der Vorst, E.P.C. MicroRNAs in Chronic Kidney Disease: Four Candidates for Clinical Application. Int. J. Mol. Sci. 2020, 21, 6547. [Google Scholar] [CrossRef] [PubMed]
  31. Greco, M.; Chiefari, E.; Mirabelli, M.; Salatino, A.; Tocci, V.; Cianfrone, P.; Foti, D.P.; Brunetti, A. Plasma or Urine Neutrophil Gelatinase-Associated Lipocalin (NGAL): Which Is Better at Detecting Chronic Kidney Damage in Type 2 Diabetes? Endocrines 2022, 3, 175–186. [Google Scholar] [CrossRef]
  32. Kang, K.; Peng, X.; Luo, J.; Gou, D. Identification of circulating miRNA biomarkers based on global quantitative real-time PCR profiling. J. Anim. Sci. Biotechnol. 2012, 3, 4. [Google Scholar] [CrossRef] [Green Version]
  33. Chekka, L.M.S.; Langaee, T.; Johnson, J.A. Comparison of Data Normalization Strategies for Array-Based MicroRNA Profiling Experiments and Identification and Validation of Circulating MicroRNAs as Endogenous Controls in Hypertension. Front. Genet. 2022, 13, 836636. [Google Scholar] [CrossRef]
  34. Marima, R.; Hull, R.; Dlamini, Z.; Penny, C. The profiling, identification, quantification and analysis of differentially expressed genes (DEGs) in response to drug treatment in lung cancer. MethodsX 2021, 8, b101381. [Google Scholar] [CrossRef]
  35. La Sala, L.; Prattichizzo, F.; Ceriello, A. The link between diabetes and atherosclerosis. Eur. J. Prev. Cardiol. 2019, 26 (Suppl. S2), 15–24. [Google Scholar] [CrossRef] [PubMed]
  36. Rasoulinejad, S.A.; Akbari, A.; Nasiri, K. Interaction of miR-146a-5p with oxidative stress and inflammation in complications of type 2 diabetes mellitus in male rats: Anti-oxidant and anti-inflammatory protection strategies in type 2 diabetic retinopathy. Iran. J. Basic Med. Sci. 2021, 24, 1078–1086. [Google Scholar] [CrossRef] [PubMed]
  37. Lo, W.Y.; Wang, S.J.; Wang, H.J. Non-canonical Interaction Between O-Linked N-Acetylglucosamine Transferase and miR-146a-5p Aggravates High Glucose-Induced Endothelial Inflammation. Front. Physiol. 2020, 11, 1091. [Google Scholar] [CrossRef] [PubMed]
  38. Alipoor, B.; Ghaedi, H.; Meshkani, R.; Torkamandi, S.; Saffari, S.; Iranpour, M.; Omrani, M.D. Association of MiR-146a Expression and Type 2 Diabetes Mellitus: A Meta-Analysis. Int. J. Mol. Cell Med. 2017, 6, 156–163. [Google Scholar] [CrossRef] [PubMed]
  39. Yung, J.H.M.; Giacca, A. Role of c-Jun N-terminal Kinase (JNK) in Obesity and Type 2 Diabetes. Cells 2020, 9, 706. [Google Scholar] [CrossRef] [Green Version]
  40. Régnier, C.H.; Song, H.Y.; Gao, X.; Goeddel, D.V.; Cao, Z.; Rothe, M. Identification and characterization of an IkappaB kinase. Cell 1997, 90, 373–383. [Google Scholar] [CrossRef] [Green Version]
  41. Woronicz, J.D.; Gao, X.; Cao, Z.; Rothe, M.; Goeddel, D.V. IkappaB kinase-beta: NF-kappaB activation and complex formation with IkappaB kinase-alpha and NIK. Science 1997, 278, 866–869. [Google Scholar] [CrossRef]
  42. Sgarra, R.; Pegoraro, S.; Ros, G.; Penzo, C.; Chiefari, E.; Foti, D.; Brunetti, A.; Manfioletti, G. High Mobility Group A (HMGA) proteins: Molecular instigators of breast cancer onset and progression. Biochim. Biophys. Acta Rev. Cancer 2018, 1869, 216–229. [Google Scholar] [CrossRef]
  43. Chiefari, E.; Foti, D.P.; Sgarra, R.; Pegoraro, S.; Arcidiacono, B.; Brunetti, F.S.; Greco, M.; Manfioletti, G.; Brunetti, A. Transcriptional Regulation of Glucose Metabolism: The Emerging Role of the HMGA1 Chromatin Factor. Front. Endocrinol. 2018, 9, 357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Brunetti, A.; Brunetti, L.; Foti, D.; Accili, D.; Goldfine, I.D. Human diabetes associated with defects in nuclear regulatory proteins for the insulin receptor gene. J. Clin. Investig. 1996, 97, 258–262. [Google Scholar] [CrossRef] [PubMed]
  45. Chiefari, E.; Nevolo, M.T.; Arcidiacono, B.; Maurizio, E.; Nocera, A.; Iiritano, S.; Sgarra, R.; Possidente, K.; Palmieri, C.; Paonessa, F.; et al. HMGA1 is a novel downstream nuclear target of the insulin receptor signaling pathway. Sci. Rep. 2012, 2, 251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chiefari, E.; Tanyolaç, S.; Paonessa, F.; Pullinger, C.R.; Capula, C.; Iiritano, S.; Mazza, T.; Forlin, M.; Fusco, A.; Durlach, V.; et al. Functional variants of the HMGA1 gene and type 2 diabetes mellitus. JAMA. 2011, 305, 903–912. [Google Scholar] [CrossRef] [Green Version]
  47. Chanprasertyothin, S.; Jongjaroenprasert, W.; Ongphiphadhanakul, B. The Association of Soluble IGF2R and IGF2R Gene Polymorphism with Type 2 Diabetes. J. Diabetes Res. 2015, 2015, 216383. [Google Scholar] [CrossRef] [PubMed]
  48. Chen, Y.L.; Pei, D.; Hung, Y.J.; Lee, C.H.; Hsiao, F.C.; Wu, C.Z.; Lin, J.; Hsu, C.; Chang, J.; Hsieh, C. Associations between genetic variants and the severity of metabolic syndrome in subjects with type 2 diabetes. Genet. Mol. Res. 2015, 14, 2518–2526. [Google Scholar] [CrossRef] [PubMed]
  49. Orlowski, J.; Grinstein, S. Diversity of the mammalian sodium/proton exchanger SLC9 gene family. Pflugers Arch. 2004, 447, 549–565. [Google Scholar] [CrossRef]
  50. Bellomo, E.A.; Meur, G.; Rutter, G.A. Glucose regulates free cytosolic Zn2+ concentration, Slc39 (ZiP), and metallothionein gene expression in primary pancreatic islet β-cells. J. Biol. Chem. 2011, 286, 25778–25789. [Google Scholar] [CrossRef] [Green Version]
  51. Khananshvili, D. The SLC8 gene family of sodium-calcium exchangers (NCX)—Structure, function, and regulation in health and disease. Mol. Asp. Med. 2013, 34, 220–235. [Google Scholar] [CrossRef]
  52. Gurzov, E.N.; Stanley, W.J.; Brodnicki, T.C.; Thomas, H.E. Protein tyrosine phosphatases: Molecular switches in metabolism and diabetes. Trends Endocrinol. Metab. 2015, 26, 30–39. [Google Scholar] [CrossRef]
  53. Nasarre, L.; Juan-Babot, O.; Gastelurrutia, P.; Llucia-Valldeperas, A.; Badimon, L.; Bayes-Genis, A.; Llorente-Cortés, V. Low density lipoprotein receptor-related protein 1 is upregulated in epicardial fat from type 2 diabetes mellitus patients and correlates with glucose and triglyceride plasma levels. Acta Diabetol. 2014, 51, 23–30. [Google Scholar] [CrossRef] [PubMed]
  54. Rosa, B.L.; Vicenta, L.C.; Victor, G.C.; Dolores, L.C.; Fernando, C.; Ricardo, G.H.; Lina, B. Alterations of specific biomarkers of metabolic pathways in vascular tree from patients with Type 2 diabetes. Cardiovasc. Diabetol. 2012, 11, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Chiefari, E.; Mirabelli, M.; La Vignera, S.; Tanyolaç, S.; Foti, D.P.; Aversa, A.; Brunetti, A. Insulin Resistance and Cancer: In Search for a Causal Link. Int. J. Mol. Sci. 2021, 22, 11137. [Google Scholar] [CrossRef]
  56. Patel, S.; Santani, D. Role of NF-kappa B in the pathogenesis of diabetes and its associated complications. Pharmacol. Rep. 2009, 61, 595–603. [Google Scholar] [CrossRef]
  57. Meyerovich, K.; Ortis, F.; Cardozo, A.K. The non-canonical NF-κB pathway and its contribution to β-cell failure in diabetes. J. Mol. Endocrinol. 2018, 61, F1–F6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Donath, M.Y.; Shoelson, S.E. Type 2 diabetes as an inflammatory disease. Nat. Rev. Immunol. 2011, 11, 98–107. [Google Scholar] [CrossRef]
  59. Arcidiacono, B.; Chiefari, E.; Foryst-Ludwig, A.; Currò, G.; Navarra, G.; Brunetti, F.S.; Mirabelli, M.; Corigliano, D.M.; Kintscher, U.; Britti, D.; et al. Obesity-related hypoxia via miR-128 decreases insulin-receptor expression in human and mouse adipose tissue promoting systemic insulin resistance. eBioMedicine 2020, 59, 102912. [Google Scholar] [CrossRef]
  60. Greco, M.; Mirabelli, M.; Tocci, V.; Mamula, Y.; Salatino, A.; Brunetti, F.S.; Dragone, F.; Sicilia, L.; Tripolino, O.; Chiefari, E.; et al. Prothymosin-Alpha, a Novel and Sensitive Biomarker of the Inflammatory and Insulin-Resistant Statuses of Obese Individuals: A Pilot Study Involving Humans. Endocrines 2023, 4, 427–436. [Google Scholar] [CrossRef]
  61. De Silva, K.; Demmer, R.T.; Jönsson, D.; Mousa, A.; Forbes, A.; Enticott, J. A data-driven biocomputing pipeline with meta-analysis on high throughput transcriptomics to identify genome-wide miRNA markers associated with type 2 diabetes. Heliyon 2022, 8, e08886. [Google Scholar] [CrossRef]
  62. Zahari Sham, S.Y.; Ng, C.T.; Azwar, S.; Yip, W.K.; Abdullah, M.; Thevandran, K.; Osman, M.; Seow, H.F. Circulating miRNAs in Type 2 Diabetic Patients with and without Albuminuria in Malaysia. Kidney Blood Press. Res. 2022, 47, 81–93. [Google Scholar] [CrossRef]
  63. Singh, T.; Kaushik, M.; Mishra, L.C.; Behl, C.; Singh, V.; Tuli, H.S. Exosomal miRNAs as novel avenues for breast cancer treatment. Front. Genet. 2023, 14, 1134779. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Scatter plot analysis of differential miRNA expression on a logarithmic scale. Red and green dots represent the number of miRNAs that were found to be over-expressed (fold regulation values > 2) and under-expressed (fold regulation values < −2), respectively, in (A) IFG vs. control group; (B) new-onset T2D vs. control group, (C) new-onset T2D vs. IFG group. Black dots represent a lack of differential expression (unvaried miRNAs).
Figure 1. Scatter plot analysis of differential miRNA expression on a logarithmic scale. Red and green dots represent the number of miRNAs that were found to be over-expressed (fold regulation values > 2) and under-expressed (fold regulation values < −2), respectively, in (A) IFG vs. control group; (B) new-onset T2D vs. control group, (C) new-onset T2D vs. IFG group. Black dots represent a lack of differential expression (unvaried miRNAs).
Diagnostics 13 02443 g001
Figure 2. Global Ct mean normalization of miRNA expression and intergroup comparison. The fold regulation values are plotted on the X axis. The Y axis starts at the arbitrary cut-off value of 10 for up-regulation that is considered to be indicative of a significantly different miRNA expression level between both the IFG and new-onset T2D patient groups and the euglycemic control group. Global Ct mean normalization was performed using the GeneGlobe software (https://geneglobe.qiagen.com).
Figure 2. Global Ct mean normalization of miRNA expression and intergroup comparison. The fold regulation values are plotted on the X axis. The Y axis starts at the arbitrary cut-off value of 10 for up-regulation that is considered to be indicative of a significantly different miRNA expression level between both the IFG and new-onset T2D patient groups and the euglycemic control group. Global Ct mean normalization was performed using the GeneGlobe software (https://geneglobe.qiagen.com).
Diagnostics 13 02443 g002
Figure 3. Progressive increase in specific miRNA expression levels in the evolution from IFG to overt T2D with respect to the euglycemic status. The fold regulation values—calculated using the euglycemic control group as reference for both the IFG and new-onset T2D patient groups—are plotted on the X axis. Global Ct mean normalization was performed using the GeneGlobe software (https://geneglobe.qiagen.com).
Figure 3. Progressive increase in specific miRNA expression levels in the evolution from IFG to overt T2D with respect to the euglycemic status. The fold regulation values—calculated using the euglycemic control group as reference for both the IFG and new-onset T2D patient groups—are plotted on the X axis. Global Ct mean normalization was performed using the GeneGlobe software (https://geneglobe.qiagen.com).
Diagnostics 13 02443 g003
Figure 4. Hypothetical mechanisms by which some of the circulating deregulated miRNAs described in this study may cause insulin resistance (A), inflammation and vascular damage (B), leading to susceptibility to T2D development. Circulating miRNAs and their predicted target genes are indicated with matched colors.
Figure 4. Hypothetical mechanisms by which some of the circulating deregulated miRNAs described in this study may cause insulin resistance (A), inflammation and vascular damage (B), leading to susceptibility to T2D development. Circulating miRNAs and their predicted target genes are indicated with matched colors.
Diagnostics 13 02443 g004
Table 1. General characteristics and laboratory parameters of study participants.
Table 1. General characteristics and laboratory parameters of study participants.
Control
(n = 10)
IFG
(n = 10)
New-Onset T2D
(n = 10)
p Value
Female sex (%)6 (60%)6 (60%)6 (60%)-
Age (yrs)53.50 (47–58)51.50 (47–61)58 (57–65)0.167
BMI (kg/m2)26 (22–28)28 (25–30)29 (27–30)0.060
SBP (mmHg)120 (119–120)120 (120–130)120 (120–125)0.170
DBP (mmHg)80 (78–80)85 (80–85)80 (80–85)0.048
Fasting glucose (mg/dL)88 (83–91)113 (105–118)135 (118–165)<0.001
HbA1c (%)5.10 (4.98–5.38)5.70 (5.43–6.10)6.20 (6.0–6.60)<0.001
Total cholesterol (mg/dL)179 (160–200)204 (178–223)185 (165–210)0.337
LDL-cholesterol (mg/dL)110 (93–125)142 (109–158)113 (103–125)0.133
HDL-cholesterol (mg/dL)57 (44–60)42 (37–56)55 (40–56)0.207
Triglycerides (mg/dL)88 (75–127)103 (67–158)120 (93–144)0.480
CRP (mg/L)3.16 (3.16–3.36)3.45 (3.16–10)4.08 (0.64–5.75)0.659
ESR (mm/h)6.50 (4.75–14.25)9.50 (3.75–28.25)9.0 (5.0–10)0.806
T2D duration (yrs)--2 (0–2)-
Data are expressed as medians and interquartile ranges (IQR). The inter-group variability is determined by Kruskal–Wallis test. Bold values denote statistical significance at p < 0.05. BMI, Body Mass Index; SBP, Systolic Blood Pressure; DBP, Diastolic Blood Pressure; HbA1c, Glycated hemoglobin; ESR, Erythrocyte Sedimentation Rate; CRP, C Reactive Protein.
Table 2. Over-expressed miRNAs in IFG vs. control group, new-onset T2D vs. control group and new-onset T2D vs. IFG group.
Table 2. Over-expressed miRNAs in IFG vs. control group, new-onset T2D vs. control group and new-onset T2D vs. IFG group.
IFG vs. ControlNew-Onset T2D vs. ControlNew-Onset T2D vs. IFG
Mature IDFold RegulationMature IDFold RegulationMature IDFold Regulation
hsa-miR-4301124.8hsa-miR-4301142.74hsa-miR-4732-3p10.29
hsa-let-7b-5p100.66hsa-let-7b-5p90.34hsa-miR-142-5p9.28
hsa-miR-26b-5p38.41hsa-miR-1287-5p50.12hsa-miR-877-3p7.75
hsa-let-7a-5p37.1hsa-let-7a-5p48.41hsa-miR-31416.43
hsa-miR-523-5p32.75hsa-miR-191342.73hsa-miR-2467-3p5.91
hsa-miR-1287-5p26.97hsa-miR-523-5p34.47hsa-miR-7-5p5.52
hsa-miR-191324.14hsa-miR-26b-5p31.06hsa-miR-324-5p5.11
hsa-miR-195-5p23.16hsa-miR-195-5p25.76hsa-miR-942-5p4.87
hsa-miR-100-5p22.37hsa-miR-158721.07hsa-miR-550a-5p4.64
hsa-miR-1260a18.94hsa-miR-3135b18.99hsa-miR-10b-5p4.33
hsa-miR-365b-3p17.43hsa-miR-1260a17.23hsa-miR-361-3p3.87
hsa-miR-4687-5p16.95hsa-miR-1260b16.53hsa-miR-21103.85
hsa-miR-3135b15.49hsa-miR-20b-5p15.97hsa-miR-548o-5p3.79
hsa-let-7f-5p14.66hsa-miR-128114.9hsa-miR-331-3p3.74
hsa-miR-313113.12hsa-miR-313114.49hsa-miR-769-5p3.66
hsa-miR-318312.67hsa-miR-4687-5p13.24hsa-miR-18a-5p3.66
hsa-miR-20b-5p10.8hsa-miR-3200-5p13.06hsa-miR-3191-3p3.64
hsa-miR-1260b10.44hsa-miR-100-5p12.88hsa-miR-140-5p3.61
hsa-miR-223-3p9.4hsa-let-7c-5p11.94hsa-miR-15873.59
hsa-miR-370-3p8.9hsa-miR-365b-3p11.77hsa-miR-130b-3p3.49
hsa-miR-3200-5p8.3hsa-let-7f-5p11.37hsa-miR-99a-5p3.35
hsa-miR-39118.02hsa-miR-671-3p11.29hsa-miR-31853.33
hsa-miR-46517.91hsa-miR-625-3p11.29hsa-miR-5963.3
hsa-miR-12817.85hsa-miR-465110.76hsa-miR-181c-5p3.26
hsa-miR-378g7.85hsa-miR-45059.9hsa-miR-138-1-3p3.17
hsa-miR-30e-3p7.8hsa-miR-144-3p9.76hsa-miR-45053.15
hsa-miR-185-5p7.75hsa-miR-1237-3p9.56hsa-miR-15393.1
hsa-miR-1237-3p7.59hsa-miR-378g9.36hsa-miR-29b-3p3
hsa-let-7c-5p7.23hsa-miR-11839.11hsa-miR-193a-5p2.92
hsa-miR-144-3p7.08hsa-miR-5968.92hsa-miR-376c-3p2.82
hsa-miR-671-3p7.08hsa-miR-1910-5p8.86hsa-miR-31762.56
hsa-miR-101-3p7.03hsa-miR-31418.44hsa-miR-11832.4
hsa-miR-1301-3p6.88hsa-miR-185-5p8.15hsa-miR-501-5p2.38
hsa-miR-1910-5p6.7hsa-miR-36107.93hsa-miR-574-3p2.32
hsa-miR-12906.38hsa-miR-877-3p7.71hsa-miR-1207-5p2.3
hsa-miR-36466.2hsa-miR-1207-5p7.61hsa-miR-28-3p2.27
hsa-miR-36106.16hsa-miR-39117.55hsa-miR-874-3p2.24
hsa-miR-42915.95hsa-miR-99a-5p7.55hsa-miR-92b-3p2.19
hsa-miR-39075.95hsa-miR-31837.4hsa-miR-361-5p2.19
hsa-miR-12035.95hsa-let-7e-5p7.35hsa-miR-3622a-5p2.18
hsa-miR-15875.87hsa-miR-223-3p7.2hsa-miR-339-3p2.16
hsa-miR-4732-5p5.79hsa-miR-373-5p6.31hsa-miR-3752.16
hsa-miR-625-3p5.63hsa-miR-42746.18hsa-miR-43062.13
hsa-miR-106b-5p5.36hsa-miR-370-3p6.05hsa-miR-154-5p2.12
hsa-miR-181c-3p5.25hsa-miR-101-3p5.93hsa-miR-199b-5p2.1
hsa-miR-18a-3p5.08hsa-miR-181c-5p5.69hsa-miR-454-3p2.08
hsa-miR-42745.01hsa-miR-46895.65hsa-miR-19762.02
hsa-miR-219a-1-3p4.8hsa-miR-142-5p5.57hsa-miR-378b2.01
hsa-let-7e-5p4.74hsa-miR-1301-3p5.53hsa-miR-625-3p2.01
hsa-miR-34c-3p4.54hsa-miR-106b-5p5.23
hsa-miR-877-5p4.51hsa-miR-39074.92
hsa-miR-1909-5p4.48hsa-miR-12904.81
hsa-miR-146b-5p4.45hsa-miR-877-5p4.71
hsa-miR-46884.39hsa-miR-1909-5p4.68
hsa-miR-363-3p4.39hsa-miR-4732-5p4.68
hsa-miR-490-3p4.33hsa-miR-320e4.52
hsa-miR-373-5p4.18hsa-miR-138-1-3p4.49
hsa-miR-46893.93hsa-let-7d-5p4.46
hsa-miR-42673.85hsa-miR-10a-5p4.25
hsa-miR-11833.79hsa-miR-12034.25
hsa-miR-629-5p3.77hsa-miR-36464.19
hsa-miR-10a-5p3.69hsa-miR-34c-3p4.05
hsa-miR-628-3p3.66hsa-miR-181c-3p4.02
hsa-miR-345-5p3.64hsa-miR-15393.96
hsa-let-7i-5p3.61hsa-miR-31853.96
hsa-miR-320e3.4hsa-miR-46883.96
hsa-let-7d-5p3.37hsa-miR-3191-3p3.91
hsa-miR-1207-5p3.3hsa-miR-42913.88
hsa-miR-27b-3p3.17hsa-miR-363-3p3.88
hsa-miR-45053.15hsa-miR-219a-1-3p3.83
hsa-miR-45163.12hsa-miR-628-3p3.8
hsa-miR-330-3p3.06hsa-miR-1225-3p3.62
hsa-miR-487b-3p2.94hsa-miR-2467-3p3.55
hsa-miR-425-5p2.91hsa-miR-3753.5
hsa-miR-2276-3p2.89hsa-let-7i-5p3.43
hsa-miR-675-3p2.87hsa-miR-324-3p3.36
hsa-miR-605-5p2.78hsa-miR-27b-3p3.36
hsa-miR-5962.7hsa-miR-45163.33
hsa-miR-19a-3p2.7hsa-miR-92b-3p3.2
hsa-miR-451a2.65hsa-miR-675-3p3.2
hsa-miR-7-1-3p2.54hsa-miR-490-3p3.15
hsa-miR-127-3p2.5hsa-miR-146b-5p3.09
hsa-miR-324-3p2.43hsa-miR-42672.98
hsa-miR-193b-3p2.33hsa-miR-425-3p2.96
hsa-miR-378a-5p2.3hsa-miR-942-5p2.92
hsa-miR-99a-5p2.26hsa-miR-374c-5p2.86
hsa-miR-374c-5p2.19hsa-miR-31762.8
hsa-miR-425-3p2.16hsa-miR-550a-5p2.78
hsa-miR-489-3p2.15hsa-miR-2276-3p2.75
hsa-miR-1247-5p2.12hsa-miR-425-5p2.75
hsa-miR-188-5p2.1hsa-miR-190a-5p2.67
hsa-miR-44222.06hsa-miR-10b-5p2.6
hsa-miR-664a-3p2.06hsa-miR-139-3p2.58
hsa-miR-7-2-3p2hsa-miR-188-5p2.46
hsa-miR-45382hsa-miR-451a2.41
hsa-miR-1247-5p2.39
hsa-miR-664a-3p2.37
hsa-miR-361-3p2.32
hsa-miR-21102.31
hsa-miR-605-5p2.29
hsa-miR-43022.29
hsa-miR-1180-3p2.28
hsa-miR-548o-5p2.28
hsa-miR-193b-3p2.28
hsa-miR-331-3p2.25
hsa-miR-7-1-3p2.25
hsa-miR-769-5p2.2
hsa-miR-45382.18
hsa-miR-140-5p2.17
hsa-miR-19a-3p2.17
hsa-miR-378b2.04
hsa-miR-44542.04
Table 3. Under-expressed miRNAs in IFG vs. control group, new-onset T2D vs. control group and new-onset T2D vs. IFG group.
Table 3. Under-expressed miRNAs in IFG vs. control group, new-onset T2D vs. control group and new-onset T2D vs. IFG group.
IFG vs. ControlNew-Onset T2D vs. ControlNew-Onset T2D vs. IFG
Mature IDFold RegulationMature IDFold RegulationMature IDFold Regulation
hsa-miR-328-3p−52.22hsa-miR-328-3p−30.54hsa-miR-629-5p−6.13
hsa-miR-324-5p−15.1hsa-miR-598-3p−13.96hsa-miR-152-3p−5.53
hsa-miR-154-5p−14.89hsa-miR-152-3p−12.84hsa-miR-1307-3p−5.3
hsa-miR-382-5p−14.59hsa-miR-502-3p−11.74hsa-miR-3651−5.23
hsa-miR-7-5p−11.05hsa-miR-382-5p−11.49hsa-miR-345-5p−5.08
hsa-miR-409-3p−9.49hsa-miR-409-3p−11.49hsa-miR-4323−4.52
hsa-miR-210-3p−8.44hsa-miR-181d-5p−8.07hsa-miR-30e-3p−4.25
hsa-miR-181d-5p−8.26hsa-miR-28-5p−7.91hsa-miR-487b-3p−4.05
hsa-miR-28-5p−8.09hsa-miR-1307-3p−7.58hsa-miR-18a-3p−3.83
hsa-miR-4258−7.71hsa-miR-154−5p-7.03hsa-miR-378a-5p−3.75
hsa-miR-340-5p−6.85hsa-miR-424-3p−6.98hsa-miR-502-3p−3.72
hsa-miR-627-5p−6.8hsa-miR-132-3p−6.79hsa-miR-181a-5p−3.67
hsa-miR-4732-3p−6.35hsa-miR-181a-5p−6.65hsa-miR-378e−3.2
hsa-miR-1976−6.26hsa-miR-30a-3p−6.56hsa-miR-590-5p−3
hsa-miR-151a-5p−5.92hsa-miR-4323−6.25hsa-miR-421−2.84
hsa-miR-34a-5p−5.53hsa-miR-3651−6.03hsa-miR-301a-3p−2.84
hsa-miR-29b-3p−5.53hsa-miR-210-3p−5.91hsa-miR-206−2.74
hsa-miR-132-3p−5.53hsa-miR-627-5p−5.87hsa-let-7a-3p−2.72
hsa-miR-224-5p−5.49hsa-miR-339-5p−5.47hsa-miR-598-3p−2.72
hsa-miR-598-3p−5.12hsa-miR-34a-5p−5.4hsa-miR-374a-5p−2.72
hsa-miR-339-5p−4.98hsa-miR-199a-5p−4.87hsa-miR-26b-3p−2.67
hsa-miR-130b-3p−4.68hsa-miR-4258−4.57hsa-miR-192-5p−2.67
hsa-miR-130a-3p−4.58hsa-miR-378a-3p−4.45hsa-miR-505-3p−2.61
hsa-miR-18a-5p−4.31hsa-miR-340-5p−4.12hsa-miR-199a-5p−2.6
hsa-miR-424-3p−3.99hsa-miR-93-3p−3.87hsa-miR-3923−2.58
hsa-miR-361-5p−3.94hsa-miR-130a-3p−3.69hsa-miR-3653-3p−2.51
hsa-miR-30a-3p−3.88hsa-miR-15a-5p−3.35hsa-miR-15a-5p−2.49
hsa-miR-199b-5p−3.75hsa-miR-151a-5p−3.3hsa-miR-1-3p−2.42
hsa-miR-16-5p−3.22hsa-miR-421−3.26hsa-miR-378a-3p−2.2
hsa-miR-342-3p−3.17hsa-miR-199b-3p−3.19hsa-miR-485-5p−2.15
hsa-miR-93-5p−3.17hsa-miR-1976−3.1hsa-miR-338-3p−2.12
hsa-miR-150-5p−3.15hsa-miR-194-5p−3.06hsa-miR-199b-3p−2.08
hsa-miR-502-3p−3.15hsa-miR-324-5p−2.95hsa-miR-151b−2.05
hsa-miR-194-5p−2.98hsa-miR-93-5p−2.93
hsa-miR-5095−2.98hsa-miR-125a-5p−2.81
hsa-miR-193a-5p−2.9hsa-miR-191-5p−2.78
hsa-miR-125a-5p−2.84hsa-miR-224-5p−2.76
hsa-miR-433-3p−2.65hsa-miR-151a-3p−2.76
hsa-miR-24-3p−2.6hsa-miR-24-3p−2.76
hsa-miR-19b-3p−2.6hsa-miR-192-5p−2.74
hsa-miR-186-5p−2.58hsa-miR-19b-3p−2.68
hsa-miR-130b-5p−2.56hsa-miR-16-5p−2.66
hsa-miR-151a-3p−2.54hsa-miR-146a-5p−2.64
hsa-miR-152-3p−2.32hsa-miR-186-5p−2.61
hsa-miR-484−2.31hsa-miR−433-3p−2.59
hsa-miR-320a−2.26hsa-miR-126-3p−2.55
hsa-miR-93-3p−2.2hsa-miR-505-3p−2.52
hsa-let-7g-5p−2.14hsa-miR-21-5p−2.37
hsa-miR-92a-3p−2.14hsa-miR-151b−2.35
hsa-let-7d-3p−2.12hsa-miR-25-3p−2.16
hsa-miR-25-3p−2.07hsa-miR-342-3p−2.02
hsa-miR-423-5p−2.02hsa-miR-15b-5p−2.02
hsa-miR-191-5p−2.02hsa-miR-7-5p−2
hsa-miR-378a-3p−2.02
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Greco, M.; Mirabelli, M.; Salatino, A.; Accattato, F.; Aiello, V.; Brunetti, F.S.; Chiefari, E.; Pullano, S.A.; Fiorillo, A.S.; Foti, D.P.; et al. From Euglycemia to Recent Onset of Type 2 Diabetes Mellitus: A Proof-of-Concept Study on Circulating microRNA Profiling Reveals Distinct, and Early microRNA Signatures. Diagnostics 2023, 13, 2443. https://doi.org/10.3390/diagnostics13142443

AMA Style

Greco M, Mirabelli M, Salatino A, Accattato F, Aiello V, Brunetti FS, Chiefari E, Pullano SA, Fiorillo AS, Foti DP, et al. From Euglycemia to Recent Onset of Type 2 Diabetes Mellitus: A Proof-of-Concept Study on Circulating microRNA Profiling Reveals Distinct, and Early microRNA Signatures. Diagnostics. 2023; 13(14):2443. https://doi.org/10.3390/diagnostics13142443

Chicago/Turabian Style

Greco, Marta, Maria Mirabelli, Alessandro Salatino, Francesca Accattato, Vincenzo Aiello, Francesco S. Brunetti, Eusebio Chiefari, Salvatore A. Pullano, Antonino S. Fiorillo, Daniela P. Foti, and et al. 2023. "From Euglycemia to Recent Onset of Type 2 Diabetes Mellitus: A Proof-of-Concept Study on Circulating microRNA Profiling Reveals Distinct, and Early microRNA Signatures" Diagnostics 13, no. 14: 2443. https://doi.org/10.3390/diagnostics13142443

APA Style

Greco, M., Mirabelli, M., Salatino, A., Accattato, F., Aiello, V., Brunetti, F. S., Chiefari, E., Pullano, S. A., Fiorillo, A. S., Foti, D. P., & Brunetti, A. (2023). From Euglycemia to Recent Onset of Type 2 Diabetes Mellitus: A Proof-of-Concept Study on Circulating microRNA Profiling Reveals Distinct, and Early microRNA Signatures. Diagnostics, 13(14), 2443. https://doi.org/10.3390/diagnostics13142443

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop