Next Article in Journal / Special Issue
Cell Type-Selective Expression of Circular RNAs in Human Pancreatic Islets
Previous Article in Journal
Herpes Simplex Virus 1 Deregulation of Host MicroRNAs
Previous Article in Special Issue
Influence of Disease Duration on Circulating Levels of miRNAs in Children and Adolescents with New Onset Type 1 Diabetes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
ncRNA
Volume 4
Issue 4
10.3390/ncrna4040037
ncrna-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Serum Levels of miR-148a and miR-21-5p Are Increased in Type 1 Diabetic Patients and Correlated with Markers of Bone Strength and Metabolism

by
Giuseppina E. Grieco
1,2,
Dorica Cataldo
3,
Elena Ceccarelli
1,
Laura Nigi
1,2,
Giovanna Catalano
1,
Noemi Brusco
1,2,
Francesca Mancarella
1,2,
Giuliana Ventriglia
1,2,
Cecilia Fondelli
3,
Elisa Guarino
3,
Isabella Crisci
3,
Guido Sebastiani
1,2,* and
Francesco Dotta
1,2,3
1
Department of Medicine, Surgery and Neurosciences, University of Siena, 53100 Siena, Italy
2
Fondazione Umberto Di Mario ONLUS c/o Toscana Life Science, 53100 Siena, Italy
3
UOC Diabetologia, Azienda Ospedaliera Universitaria Senese, 53100 Siena, Italy
*
Author to whom correspondence should be addressed.
Non-Coding RNA 2018, 4(4), 37; https://0-doi-org.brum.beds.ac.uk/10.3390/ncrna4040037
Submission received: 18 October 2018 / Revised: 16 November 2018 / Accepted: 22 November 2018 / Published: 27 November 2018
(This article belongs to the Special Issue Non-Coding RNA and Diabetes)
Browse Figures
Versions Notes

Abstract

:
Type 1 diabetes (T1D) is characterized by bone loss and altered bone remodeling, resulting into reduction of bone mineral density (BMD) and increased risk of fractures. Identification of specific biomarkers and/or causative factors of diabetic bone fragility is of fundamental importance for an early detection of such alterations and to envisage appropriate therapeutic interventions. MicroRNAs (miRNAs) are small non-coding RNAs which negatively regulate genes expression. Of note, miRNAs can be secreted in biological fluids through their association with different cellular components and, in such context, they may represent both candidate biomarkers and/or mediators of bone metabolism alterations. Here, we aimed at identifying miRNAs differentially expressed in serum of T1D patients and potentially involved in bone loss in type 1 diabetes. We selected six miRNAs previously associated with T1D and bone metabolism: miR-21; miR-24; miR-27a; miR-148a; miR-214; and miR-375. Selected miRNAs were analyzed in sera of 15 T1D patients (age: 33.57 ± 8.17; BMI: 21.4 ± 1.65) and 14 non-diabetic subjects (age: 31.7 ± 8.2; BMI: 24.6 ± 4.34). Calcium, osteocalcin, parathormone (PTH), bone ALkaline Phoshatase (bALP), and Vitamin D (VitD) as well as main parameters of bone health were measured in each patient. We observed an increased expression of miR-148a (p = 0.012) and miR-21-5p (p = 0.034) in sera of T1D patients vs. non-diabetic subjects. The correlation analysis between miRNAs expression and the main parameters of bone metabolism, showed a correlation between miR-148a and Bone Mineral Density (BMD) total body (TB) values (p = 0.042) and PTH circulating levels (p = 0.033) and the association of miR-21-5p to Bone Mineral Content-Femur (BMC-FEM). Finally, miR-148a and miR-21-5p target genes prediction analysis revealed several factors involved in bone development and remodeling, such as MAFB, WNT1, TGFB2, STAT3, or PDCD4, and the co-modulation of common pathways involved in bone homeostasis thus potentially assigning a role to both miR-148a and miR-21-5p in bone metabolism alterations. In conclusion, these results lead us to hypothesize a potential role for miR-148a and miR-21-5p in bone remodeling, thus representing potential biomarkers of bone fragility in T1D.

1. Introduction

Risk of bone fractures is significantly higher in Type 1 diabetes (T1D) respect to the non-diabetic general population [1]. Indeed, prevalence of hip fractures in T1D patients have been reported to be up to six-fold higher than in the general population and 2.5-fold higher respect to type 2 diabetic (T2D) patients, thus suggesting the presence of specific and critical bone remodeling alterations in autoimmune diabetes [1,2].
One of the main parameters used to measure bone health and structure is represented by bone mineral density (BMD). Interestingly, while BMD was found significantly decreased in T1D patients [3,4], it is often normal or even increased in T2D patients [5,6], indicating the complex interplay among mechanisms characterizing T1D, T2D, and bone alterations.
The prominent increase of bone fractures and bone fragility in T1D is characterized by multiple determinants. For example, pro-inflammatory cytokines, which are major mediators of several diabetic complications, interfere with osteoblast differentiation and bone matrix collagen deposition [7], and contribute to increased bone microvascular dysfunctions [8], thus leading to an altered bone remodeling. In the same way, other common determinants, such as adipokines and bone marrow adiposity [9], hyperglycemia [10], and/or insulin deficiency [11,12], lead to a lower bone turnover, thus increasing the risk of fractures in diabetic patients.
An additional class of molecules potentially involved in the development of many diabetic complications is represented by microRNAs (miRNA). miRNAs are endogenous non-coding small RNAs, 19–24 nucleotides long, which negatively regulate gene expression through their specific binding to 3’UTR region of messenger RNA [13]. miRNAs have been reported to regulate many cellular processes through the modulation of entire signaling pathways. Consequently, alterations of several miRNAs, reported in different diseases, can surely contribute to define a disease phenotype. As a matter of fact, specific miRNA alterations have been associated to T1D both in pancreatic islets [14] and in immune cells [15,16], thus characterizing diabetes outcome [15,17]. Additionally, miRNAs have been found altered in several tissues affected by typical secondary diabetic dysfunctions (e.g., vascular endothelium-miR-126) [18]. Moreover, miRNAs participate in several mechanisms involved in bone development and function, modeling and remodeling, and characterize the phenotype of bone cells, including osteocytes, osteoblasts, and osteoclasts [19]. Of note, miRNAs have been demonstrated to be involved also in bone diseases [20].
Recent studies demonstrated that miRNAs are not confined within cells of origin but can be also secreted, thus being found in different biological fluids such as serum or plasma [21]. In serum/plasma they can be found associated with different cellular components such as protein-bound miRNAs (Ago2 or HDL) or within extracellular vesicles (EVs) (e.g., exosomes) [21]. Although the significance of specific miRNA secretion in biological fluids has not been clearly defined yet, several studies demonstrated that: (i) blood circulating miRNA expression profiles are altered in a variety of diseases [22]; (ii) alterations of multiple circulating miRNAs can define a tissue-specific dysfunction thus being potentially used as disease biomarkers [23]; and (iii) circulating miRNAs are not only byproducts of cell loss but have been reported to be active molecules which mediate communication among different tissues and cells [24]. For these reasons, miRNAs may represent both optimal biomarkers of specific tissue dysfunction and potential molecules orchestrating distant molecular mechanisms [25]. Therefore, the identification of an alteration of a specific set of circulating miRNAs may pave the way for the discovery of new biomarkers and for the identification of new therapeutic targets [26].
Here, we hypothesize that a specific set of circulating miRNAs, found altered in serum/plasma of T1D patients, may represent also optimal biomarkers of bone fragility and/or mediators of bone tissue dysfunctions, given their previously reported role both in autoimmune diabetes and in bone homeostasis, in the regulation of osteoclastogenesis, osteoblastogenesis, and their modulation.

2. Results

2.1. A Specific Set of miRNAs Is Altered in T1D and Involved in Bone Metabolism

An extended literature search was initially performed in order to identify a specific set of candidate circulating miRNAs whose expression has been previously associated with T1D and which have also been proposed as modulators of bone homeostasis, being involved in bone development, modeling and/or remodeling. Then, we looked for their previously reported role in bone cells/metabolism. We selected six miRNAs, all significantly dysregulated in serum or plasma of T1D patients, namely: miR-148a [27,28,29,30], miR-21-5p [28,29], miR-214 [31], miR-375 [28,29,32], miR-27a [27], and miR-24 [27,28,29,31,32]. All have a clear and demonstrated role in osteoblastogenesis or osteoclastogenesis, thus being involved in function of osteoblasts or osteoclasts, or found altered in serum/plasma of patients with bone diseases (e.g., osteoporosis) (Table 1) [33,34,35,36,37,38,39,40]. Above mentioned miRNAs were taken into consideration for further expression analysis in sera collected from a selected cohort of T1D patients and control subjects who were previously characterized in terms of parameters of bone homeostasis and bone strength.

2.2. CirculatingLevels of miR-148a and miR-21-5p Are Increased in Sera of T1D Patients vs. Control Subjects

In order to analyze the expression of the selected miRNAs reported in Table 1, we took into consideration a specific cohort of T1D patients and non-diabetic volunteers, already characterized for a specific panel of bone metabolic markers and bone physical parameters. T1D patients and healthy subjects were further selected based on age and gender, in order to perfectly match the two groups (Table 2 and Table S1). Of note, as previously reported by our group [41] and in line with several studies published by others [42,43], T1D patients showed a significant reduction of Bone Mineral Density-Femur (BMD-FEM), Bone Mineral Content-Femur (BMC FEM), T score FEM and Z score FEM as well as Bone Mineral Density-Neck (BMD-N), T score N and Z score N vs. non-diabetic subjects, confirming the alterations of bone metabolism in T1D. Additionally, we observed a significant reduction of 25-hydroxy vitamin-D (Vit D) and a concomitant increase in parathormone (PTH) in T1D patients vs. non-diabetic individuals (Table 2). Collectively, these results highlight a reduction of bone strength in T1D patients.
Serum samples from the same cohort of T1D patients (n = 15) and non-diabetic controls (n = 14) were subjected to expression analysis of the 6 selected miRNAs, namely miR-148a, miR-21-5p, miR-214, miR-375, miR-27a, and miR-24, using RT-Stem Loop Real Time PCR. Among selected miRNAs, miR-148a and miR-21-5p expression levels resulted significantly increased in sera of T1D patients vs. non-diabetic controls, while no significant difference was observed for other miRNAs (Figure 1).

2.3. Circulating miR-148a and miR-21-5p Levels Are Associated with BMD-Total Body (TB), PTH, and BMC-FEM

In an attempt to get insight into the potential role of selected circulating microRNAs in bone metabolism, we correlated their expression levels to the main parameters of bone homeostasis and metabolism measured in T1D patients and non-diabetic subjects. Initially, we took into consideration miR-148a and miR-21-5p due to their peculiar differential expression in sera of T1D patients. No significant correlations were identified for miR-148a and miR-21-5p respect to the main clinical T1D patient characteristics (age, BMI, disease duration and autoantibodies positivity) (Figure S1). Conversely, miR-148a correlation analyses with bone metabolic parameters (Figure S2 and Table S2), showed a significant negative correlation between serum miR-148a and BMD-TB (Figure 2a), thus suggesting a potential association of miR-148a with the regulation of bone strength (Figure 2a). Concomitantly, miR-148a was also positively correlated with circulating PTH (Figure 2b), reinforcing the potential link between miR-148a and bone metabolism. Additionally, a significant correlation between miR-21-5p and BMC-Fem was observed as well (Figure 2c), thus highlighting the potential role of miR-21-5p in the modulation of bone strength, in line with miR-148a.
Finally, the correlation analysis between selected miRNAs and parameters of bone metabolism revealed a positive significant association of miR-375 expression levels with PTH (Supplementary Figure S2).

2.4. Bioinformatic Analysis Revealed that miR-148a and miR-21-5p Target Several Genes Involved in Bone Metabolism and Bone Remodeling

Given the observed correlation between miR-148a and miR-21-5p with bone strength parameters and PTH levels, we aimed at evaluating a potential mechanism in which this miRNAs could be involved in the modulation of bone modeling and/or function. To this purpose we looked for predicted and experimentally validated miR-148a and miR-21-5p target genes which could directly or indirectly modulate bone metabolism. Using two manually curated databases (Tarbase and miRTarbase) exclusively containing experimentally validated target genes, and revising previous published studies, we detected several miR-148a and miR-21-5p target genes reported to be involved in bone functions and remodeling (Table 3). Among these, we identified MAFB [44] and WNT1 [45,46] which have been reported to be pivotal modulators of bone remodeling through the regulation of osteoclasts [44] and osteoblasts [47] differentiation, respectively, thus suggesting a potential role for miR-148a in the regulation of bone resorption and deposition in physiologic conditions and in T1D. Additionally, miR-21-5p has been reported to target several genes involved in bone remodeling; among them, we identified STAT3 ascribed as a pivotal modulator of bone metabolism and remodeling [48,49] (Table 3).
Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) analysis, performed through the use of Diana tool miRPath v3.0 software and aimed at identifying miR-148a and miR-21-5p common regulated pathways, identified 1372 validated target genes for hsa-miR-21-5p and 1400 validated target genes for hsa-miR-148a. Most importantly, 24 validated target genes of miR-21-5p and 28 validated target genes of miR-148a-3p are commonly involved in FoxO signaling pathway (p = 0.00087). Moreover, 15 validated target genes of miR-148a-3p are also involved in TGF-β signaling pathway (p = 0.0051) (Figure 3 and Supplementary Table S3). Of note, both FoxO [50,51,52,53,54] and TGF-β [55,56,57,58] signaling pathways have been largely reported to be involved in bone metabolism, resorption and remodeling thus highlighting the putative role of miR-148a and miR-21-5p in the modulation of bone metabolism through the regulation of common pathways.

3. Discussion

To our knowledge, this is the first study which associates circulating miRNAs expression to bone metabolism alterations in type 1 diabetes (T1D). Multiple studies suggested that circulating miRNAs could be novel potential disease biomarkers [22], also critically involved in the regulation of molecular mechanisms at distant sites from their cells of origin [25]. Such properties make miRNAs ideal predictive, diagnostic and follow-up biomarkers, as well as putative therapeutic targets.
In T1D, several studies identified multiple alterations of circulating miRNAs, even though the results were in some cases discordant [17,26]. Of note, circulating miRNAs have been found associated with several life-threatening T1D complications, thus being suggested as potential therapeutic targets or as precious biomarkers for an early detection respect to the clinical onset of such complications [18,59,60].
We and others have previously demonstrated that T1D patients show an increased risk of fractures [61] due to peculiar alterations of bone metabolism/strength and remodeling [41]. Indeed, such patients were characterized by lower BMD and BMC, higher PTH levels and by an additional plethora of alterations of bone metabolism which highlighted a lower bone turnover, suggesting, in turn, an alteration of bone modeling and remodeling in autoimmune diabetes. However, although BMD represents a good indicator of risk of fractures, it could underestimate it, thus making risk assessment a challenging analysis. To this regard, the introduction of novel biomarkers, such as miRNAs, may help in risk assessment and, potentially, in an accurate prediction and evaluation of bone fragility. Moreover, among several potential causes of bone fragility that have been suggested, at present a role for circulating miRNAs in T1D as mediators of bone metabolic alterations has not yet been proposed.
As a matter of fact, in T1D as well as in bone diseases, miRNAs have been shown to be differentially expressed in tissues (e.g., pancreatic islet cells, immune cells, and bone cells) as well as in serum/plasma [20,29], thus underlining the importance of evaluating circulating miRNAs both as biomarkers and/or as causative/contributing factors to a specific disease.
In the present study, a set of miRNAs selected based on their reported involvement both in T1D pathogenesis and in bone metabolism/bone remodeling (miR-148a, miR-21-5p, miR-214, miR-375, miR-27a, and miR-24) were evaluated in serum samples of 15 T1D patients and of 14 age-matched control subjects; results showed and confirmed the hyperexpression of miR-148a and miR-21-5p in sera of T1D patients vs. non-diabetic individuals. miR-148a and miR-21-5p were previously observed as differentially expressed in T1D serum/plasma; in the present study such differential expression is confirmed and additionally associated to bone fragility and bone metabolism alterations in T1D, being inversely correlated with BMD-TB (miR-148a) and BMC-FEM (miR-21-5p) and positively associated with PTH levels (miR-148a).
Particularly, hyperexpression of miR-148a and miR-21-5p in serum/plasma from T1D patients were previously observed in several different studies which evaluated circulating miRNAs expression profiles in T1D patients cohort of similar composition (number of T1D patients and non-diabetic controls and clinical characteristics) with respect to the present study [27,28,29,30,32]. Of note, in two independent studies, miR-148a was also reported to be increased in plasma of osteoporotic (OP) postmenopausal women [33] and in sera from OP patients with reported hip fractures [34], thus suggesting a link between miR-148a, T1D and bone remodeling. Similarly, miR-21-5p was previously associated to bone remodeling [62] and its circulating expression signature was suggested as a potential biomarker of bone diseases [63]. The analysis of clinical characteristics of T1D patients led us to exclude a role for some major confounders (e.g., diabetic complications, other inflammatory diseases, etc.—see Table S1) in dictating miR-148a and miR-21-5p expression levels, since no differential expression were detected as associated with these clinical characteristics (Supplementary Table S1).
Although such results are in line to what was previously reported and demonstrated for miR-148a and miR-21-5p expression in T1D and in bone diseases, our data should be further confirmed in additional patients and controls thus aiming at reinforcing such results. Similarly, the correlations between miR-148a and miR-21-5p with parameters of bone metabolism should be performed and confirmed by taking into consideration an extended cohort of patients. In fact, although the associations of such miRNAs to bone alterations interestingly suggest their potential involvement in bone homeostasis, our cohort did not reach the statistical power to unequivocally correlate them exclusively in T1D patients or in non-diabetic control cohort, thus representing a limit of our study. On the other hand, our results indicate that miR-148a and miR-21-5p are significantly associated to parameters of bone metabolism in T1D patients and non-diabetic controls as a potential global effect of a physiological regulatory mechanism which could be exacerbated and altered in T1D patients due to hyperexpression of these circulating miRNAs. Of note, although the correlations between miR-148a/miR-21-5p and parameters of bone metabolism and strength did not reach statistical significance by analyzing only T1D group, such correlations showed an interesting tendency (data not shown) which should be further examined in future studies.
An additional potential issue regarding circulating miRNAs is the origin of their secretion. Indeed, although at present we cannot decipher the cells of origin of circulating miR-148a or miR-21-5p, a previously published study revealed that an increased expression of miR-148a occurred in CD14+ osteoclasts precursor cells derived from lupus patients [44]. Such observation suggests that circulating miR-148a alteration might represent a common phenomenon occurring in autoimmune diseases which potentially affects bone strength. As a matter of fact, in the same study, it was demonstrated that miR-148a negatively regulates MAFB, which, in turn, negatively modulates Receptor Activator of Nuclear factor Kappa-Β Ligand (RANKL) induced osteoclastogenesis, thus resulting in an increased rate of osteoclasts differentiation and, consequently, elevated bone resorption. Collectively, these data suggest that the increased circulating levels of miR-148a may derive from CD14+ osteoclast precursor cells, thus potentially explaining the inverse correlation between miR-148a and BMD, being the latter associated with an increased rate of bone resorption. Additionally, it has been observed that miR-148a expression is physiologically reduced during the differentiation of bone marrow derived mesenchymal stem cells (MSCs) towards osteogenic lineage [64] while its upregulation promoted adipogenic differentiation (by repressing WNT1) [46]; moreover, an increased expression of miR-148a was observed in exosomes derived from MSC cultures during osteogenic differentiation [65], showing miR-148a secretion and modulation during bone modeling. Altogether, these data suggest that in-situ miR-148a hyperexpression may be detrimental to bone modeling and remodeling and that miR-148a may be secreted through exosomes, thus being potentially detectable in biological fluids.
As for miR-21-5p, it has been previously extensively demonstrated its role in the modulation of osteoclastogenesis [66]. Of note, a recent report demonstrated that circulating miR-21-5p secreted through exosomes may control osteoclastogenesis [67], thus underlining the potential role of this circulating miRNA in the modulation of bone remodeling or as promising biomarker [68].
Collectively, these evidences suggest that miR-148a and miR-21-5p regulate common pathways involved in bone metabolism and remodeling.
In addition, we demonstrate that the expression levels of miR-148a are positively correlated with PTH. It has been shown that increased levels of miR-148a in rat parathyroid cultures induce an abnormal PTH secretion; importantly, a specific reduction of such miRNA, by using miR-148a antagomiRs, affected PTH secretion in vivo, thus suggesting an association between miR-148a levels and PTH serum concentration [69]. We can hypothesize that the increased circulating expression of miR-148a may have a broad range of effects, thus potentially modifying the secretion of PTH by targeting parathyroid cells as well. Although miR-148a has been reported to be an in vivo modulator of PTH secretion, we cannot exclude that the observed correlation between this miRNA and PTH may also derive from an indirect effect of PTH itself on bone metabolism, which may in turn alter miR-148a expression and secretion. Of note, we also identified miR-375 as positively correlated with PTH circulating levels, thus suggesting a role for miR-375 in the regulation of PTH in T1D patients and non-diabetic controls, even though its expression levels did not differ between the two groups analyzed.
In conclusion, in this study we confirmed the hyperexpression of circulating miR-148a and miR-21-5p in sera from T1D patients; additionally we reported the association of miR-148a hyperexpression to several markers of bone fragility. Further, we showed and discussed that miR-148a and miR-21-5p targets pivotal genes involved in bone remodeling pathways, thus raising the hypothesis of miR-148a/miR-21-5p as contributing factors to bone fragility or as a potential biomarkers of bone metabolic alterations in T1D.

4. Materials and Methods

4.1. Ethical Statement

All procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation (institutional and national) and with the Helsinki Declaration of 1975, as revised in 2008. The research protocol was submitted for consideration, comment, guidance and approval to the concerned research ethics committee, independently on the researcher, the sponsor and any other undue influence and is duly qualified. Every law and regulation of the country in which the research has been performed has been taken into consideration. The study was approved by the local Institutional Review Board (Comitato Etico Regionale per la Sperimentazione Clinica della Regione Toscana—Prot: 11807_2017—16/10/2017), and written informed consent was obtained from all participants as previously reported [41].

4.2. Study Population

A total of n = 15 patients with type 1 diabetes (T1D) (age range 24–49 years, mean disease duration 15.71 ± 11.33 years) referred to the Diabetes Unit of our department were included in the study. All patients had normal serum creatinine levels and no major comorbidities impairing normal daily activity. All patients were positive for at least one autoantibody at the onset of T1D. All patients were on treatment with insulin. Main clinical characteristics are reported in Table 2 and extended Table S1 which includes the clinical details of T1D patients regarding diabetic complications (2/15 T1D patients reported retinopathy; neither micro- nor macro-albuminuria was detected), current medications, other inflammatory diseases, and smoking habits alongside with miRNAs expression levels for each patient analyzed.
Age- and sex-matched controls (n = 14) were recruited from healthy volunteers. All control subjects had normal fasting glucose levels and HbA1c, 5/14 reported smoking habits, no chronic inflammatory diseases or current medications were reported.
Subjects with Paget’s disease of bone, primary hyperparathyroidism, congestive heart failure, recent myocardial infarction, multiple myeloma, or other neoplasia were excluded from the study. Moreover, subjects were also excluded if they received treatment with anti-resorptive or anabolic compounds for osteoporosis and current corticosteroid therapy, or any other treatment known to affect bone metabolism.
General and clinical characteristics of patients and controls are reported in Table 2.

4.3. Samples Collections and Clinical Analysis

Blood samples were collected in the morning after an overnight fast, in serum gel tube (BD Vacutainer) and stored at room temperature for maximum 2h. Blood samples were then centrifuged 1800×g for 15min at 4 °C and then stored at −80 °C.
Serum calcium concentrations (Ca2+) (corrected for albumin concentration), were measured using standard automated laboratory techniques. Circulating levels of intact osteocalcin [DiaSorin Diagnostics (Saluggia, Italy); interassay CV 7.1%, normal range 1.8–6.6 ng/mL] and of bone-specific alkaline phosphatase (bALP) (Beckman Coulter, Fullerton, CA, USA; with an interassay CV of_7.9%, normal range 9–21_g/L) were measured in serum samples as markers of bone turnover. Circulating PTH (DiaSorin, Stillwater, MN, USA; interassay CV _7.3%; normal range 10–60 pg/mL) and 25OH Vitamin D (DiaSorin Diagnostics; sensitivity, 1.5 ng/mL; interassay CV _11%; normal range for vitamin D sufficiency, _30 ng/mL) were evaluated by RIA.
At the time of blood sampling, areal BMD of the lumbar spine and of the proximal femur was determined by a dual-energy X-ray absorptiometry device (Lunar Prodigy; GE Healthcare, Waukesha, WI, USA).

4.4. RNA Extraction

Serum samples were thawed on ice and then further centrifuged 3000× g for 5 min at 4 °C in order to completely remove contaminant cells and cell debris. A total of 50 μL of serum from each patient was diluted in 350 μL of nuclease-free water to avoid protein aggregates. Then, RNA was extracted from diluted serum by using miRNeasy miRNA extraction kit (Qiagen, Hilden, Germany), by adding 1.2 mL of Trizol LS (Lifetechnologies, Carlsbad, CA, USA) as lysis buffer and finally eluted in 30 μL of nuclease-free water.

4.5. Single Assay qRT Real-Time PCR

The expression of miRNAs miR-21-5p, miR-24-3p, miR-214, miR-375, miR-27a and miR-148a was analyzed in all 29 serum samples through single assay qRT Real-Time PCR using TaqMan miRNA assay primers (Life technologies, CA, USA). RNA was reverse-transcribed employing Custom RT primers pool and preamplified using Custom Preamp primers pool. Briefly, 5 µL each RT or TM primer was diluted in a total volume of 500 µL Tris-EDTA 1X and used for RT or preamplification reaction. Then, 3 μL of extracted RNA were added to 6 μL of custom primers pool, 0.30 μL 100 mM dNTPs, 3 μL of 50 U/μL Multiscribe RT, 1.50 μL 10× RT Buffer, 0.19 μL 20 U/μL RNase Inhibitor and 1.01 μL H2O. The reaction product was incubated at 16 °C for 30 min, 42 °C for 30 min and then at 85 °C for 5 min. Afterwards, the synthesized cDNA was preamplified using Custom Preamp primer pool: 2.5 μL of cDNA from each sample were added to 12.5 μL 2× TaqMan Preamp Master Mix, 3.75 μL 10× Custom Preamp primers and 6.75 μL H2O. The reaction was incubated at 95 °C for 10 min, at 55 °C for 2 min and at 72 °C for 2 min, then for 12 cycles at 95 °C for 15 s and 60°C for 4 min and, finally, at 99 °C for 10 min. In each well, 5 μL of preamplified cDNA (diluted 1:8 in Tris-Edta 0.1X) were added to 15 μL reaction mix composed of 10 μL TaqMan Universal Master Mix, 1 μL of TaqMan miRNA expression assay and 4 μL of nuclease-free H2O. The reaction was incubated at 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and at 60°C for 1min.
Data analysis was performed using 2−ΔCt method; samples with resulting raw cycle-threshold (Ct) >35.0 were considered as not detected/not expressed. Normalization analysis was performed by using miRNAs hsa-miR-451 to exclude hemolyzed samples alongside with 414 nm spectrophotometer absorbance [70,71,72], hsa-miR-191 as endogenous control previously reported by other studies [73,74,75] and exogenous ath-miR-159a as technical control of RNA extraction reproducibility [76,77].

4.6. Bioinformatic Analysis

Diana tool mirPath v.3 online software (http://snf-515788.vm.okeanos.grnet.gr/) was used to identify validated target genes of miR-148a and miR-21-5p. Furthermore, through the use of the same software, KEGG pathway and Gene Ontology (GO) analysis has been performed to address a putative role for these microRNAs in molecular pathways with a specific involvement in bone metabolism and remodeling.

4.7. Statistics

Data have been reported as means or median ± SD; sample values distribution analysis was performed for all parameters by using D’Agostino & Pearson normality test. Non-parametric Mann-Whitney U Test or Student’s T-test were adopted for p-value evaluation according to values distribution; a p value ≤0.05 was considered as statistically significant.
MicroRNAs expression reported as normalized 2−ΔcT values were compared between T1D and the control group using Mann-Whitney U test (p value ≤ 0.05). Logistic regression and Spearman correlation analysis were performed to evaluate the association between miR-148a expression and measured bone parameters.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2311-553X/4/4/37/s1: Table S1: extended clinical characteristics of T1D patients, Table S2: R and P values of correlations between microRNAs expression and main bone parameters, Table S3: KEGG pathway analysis for validated target genes of miR-148a and miR-21-5p, Figure S1: correlation analysis between miR-148a and miR-21-5p expression and clinical characteristics of patients, Figure S2: correlation analysis between main parameters of bone metabolism and six selected miRNAs expression.

Author Contributions

G.E.G. and G.C. performed the experiments and contributed to the experimental design and manuscript drafting; E.C. evaluated bone metabolic parameters of T1D patients and contributed to the experimental design and manuscript drafting; D.C., L.N., C.F., E.G., and I.C. contributed to the recruitment of T1D patients and non-diabetic controls. L.N., F.M., G.V., and N.B. contributed to the experimental design and to the scientific discussion. G.S. and F.D. performed the experimental design, supervised the work and drafted the manuscript.

Funding

This work has also received funding from Italian Ministry of Research (grant no. 2015373Z39_007) and by Fondazione Roma.

Acknowledgments

We acknowledge the precious secretarial help of Maddalena Prencipe.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Janghorbani, M.; Feskanich, D.; Willett, W.C.; Hu, F. Prospective study of diabetes and risk of hip fracture: The Nurses’ Health Study. Diabetes Care 2006, 29, 1573–1578. [Google Scholar] [CrossRef] [PubMed]
  2. Weber, D.R.; Haynes, K.; Leonard, M.B.; Willi, S.M.; Denburg, M.R. Type 1 diabetes is associated with an increased risk of fracture across the life span: A population-based cohort study using The Health Improvement Network (THIN). Diabetes Care 2015, 38, 1913–1920. [Google Scholar] [CrossRef] [PubMed]
  3. Eller-Vainicher, C.; Zhukouskaya, V.V.; Tolkachev, Y.V.; Koritko, S.S.; Cairoli, E.; Grossi, E.; Beck-Peccoz, P.; Chiodini, I.; Shepelkevich, A.P. Low bone mineral density and its predictors in type 1 diabetic patients evaluated by the classic statistics and artificial neural network analysis. Diabetes Care 2011, 34, 2186–2191. [Google Scholar] [CrossRef] [PubMed]
  4. Joshi, A.; Varthakavi, P.; Chadha, M.; Bhagwat, N. A study of bone mineral density and its determinants in type 1 diabetes mellitus. J. Osteoporos. 2013, 2013, 397814. [Google Scholar] [CrossRef] [PubMed]
  5. Bonds, D.E.; Larson, J.C.; Schwartz, A.V.; Strotmeyer, E.S.; Robbins, J.; Rodriguez, B.L.; Johnson, K.C.; Margolis, K.L. Risk of fracture in women with type 2 diabetes: The Women’s Health Initiative Observational Study. J. Clin. Endocrinol. Metab. 2006, 91, 3404–3410. [Google Scholar] [CrossRef]
  6. Ma, L.; Oei, L.; Jiang, L.; Estrada, K.; Chen, H.; Wang, Z.; Yu, Q.; Zillikens, M.C.; Gao, X.; Rivadeneira, F. Association between bone mineral density and type 2 diabetes mellitus: A meta-analysis of observational studies. Eur. J. Epidemiol. 2012, 27, 319–332. [Google Scholar] [CrossRef]
  7. Gilbert, L.; He, X.; Farmer, P.; Boden, S.; Kozlowski, M.; Rubin, J.; Nanes, M.S. Inhibition of osteoblast differentiation by tumor necrosis factor-α. Endocrinology 2000, 141, 3956–3964. [Google Scholar] [CrossRef]
  8. De Waard, E.A.C.; de Jong, J.J.A.; Koster, A.; Savelberg, H.H.C.M.; van Geel, T.A.; Houben, A.J.H.M.; Schram, M.T.; Dagnelie, P.C.; van der Kallen, C.J.; Sep, S.J.S.; et al. The association between diabetes status, HbA1c, diabetes duration, microvascular disease, and bone quality of the distal radius and tibia as measured with high-resolution peripheral quantitative computed tomography—The Maastricht Study. In Osteoporosis International; Springer: Berlin, Germany, 2018. [Google Scholar]
  9. Veldhuis-Vlug, A.G.; Rosen, C.J. Mechanisms of marrow adiposity and its implications for skeletal health. Metab. Clin. Exp. 2017, 67, 106–114. [Google Scholar] [CrossRef] [Green Version]
  10. Balint, E.; Szabo, P.; Marshall, C.F.; Sprague, S.M. Glucose-induced inhibition of in vitro bone mineralization. Bone 2001, 28, 21–28. [Google Scholar] [CrossRef]
  11. Yang, J.; Zhang, X.; Wang, W.; Liu, J. Insulin stimulates osteoblast proliferation and differentiation through ERK and PI3K in MG-63 cells. Cell Biochem. Funct. 2010, 28, 334–341. [Google Scholar] [CrossRef]
  12. Gandhi, A.; Beam, H.A.; O’Connor, J.P.; Parsons, J.R.; Lin, S.S. The effects of local insulin delivery on diabetic fracture healing. Bone 2005, 37, 482–490. [Google Scholar] [CrossRef] [PubMed]
  13. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  14. Roggli, E.; Gattesco, S.; Caille, D.; Briet, C.; Boitard, C.; Meda, P.; Regazzi, R. Changes in microRNA expression contribute to pancreatic β-cell dysfunction in prediabetic NOD mice. Diabetes 2012, 61, 1742–1751. [Google Scholar] [CrossRef] [PubMed]
  15. Sebastiani, G.; Ventriglia, G.; Stabilini, A.; Socci, C.; Morsiani, C.; Laurenzi, A.; Nigi, L.; Formichi, C.; Mfarrej, B.; Petrelli, A.; et al. Regulatory T-cells from pancreatic lymphnodes of patients with type-1 diabetes express increased levels of microRNA miR-125a-5p that limits CCR2 expression. Sci. Rep. 2017, 7, 6897. [Google Scholar] [CrossRef] [PubMed]
  16. Sebastiani, G.; Grieco, F.A.; Spagnuolo, I.; Galleri, L.; Cataldo, D.; Dotta, F. Increased expression of microRNA miR-326 in type 1 diabetic patients with ongoing islet autoimmunity. Diabetes Metab. Res. Rev. 2011, 27, 862–866. [Google Scholar] [CrossRef] [PubMed]
  17. Dotta, F.; Ventriglia, G.; Snowhite, I.V.; Pugliese, A. MicroRNAs: Markers of β-cell stress and autoimmunity. Curr. Opin. Endocrinol. Diabetes Obes. 2018, 25, 237–245. [Google Scholar] [CrossRef]
  18. Barutta, F.; Bellini, S.; Mastrocola, R.; Bruno, G.; Gruden, G. Microrna and microvascular complications of diabetes. Int. J. Endocrinol. 2018, 2018, 6890501. [Google Scholar] [CrossRef]
  19. Chen, Z.; Zhang, Y.; Liang, C.; Chen, L.; Zhang, G.; Qian, A. Mechanosensitive miRNAs and Bone Formation. Int. J. Mol. Sci. 2017, 18, 1684. [Google Scholar] [CrossRef]
  20. Gennari, L.; Bianciardi, S.; Merlotti, D. MicroRNAs in bone diseases. Osteoporos. Int. 2017, 28, 1191–1213. [Google Scholar] [CrossRef]
  21. Turchinovich, A.; Weiz, L.; Langheinz, A.; Burwinkel, B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011, 39, 7223–7233. [Google Scholar] [CrossRef] [Green Version]
  22. Weiland, M.; Gao, X.-H.; Zhou, L.; Mi, Q.-S. Small RNAs have a large impact: Circulating microRNAs as biomarkers for human diseases. RNA Biol. 2012, 9, 850–859. [Google Scholar] [CrossRef] [PubMed]
  23. Waters, P.S.; McDermott, A.M.; Wall, D.; Heneghan, H.M.; Miller, N.; Newell, J.; Kerin, M.J.; Dwyer, R.M. Relationship between circulating and tissue microRNAs in a murine model of breast cancer. PLoS ONE 2012, 7, e50459. [Google Scholar] [CrossRef] [PubMed]
  24. Guay, C.; Regazzi, R. Exosomes as new players in metabolic organ cross-talk. Diabetes Obes. Metab. 2017, 19 (Suppl. 1), 137–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Thomou, T.; Mori, M.A.; Dreyfuss, J.M.; Konishi, M.; Sakaguchi, M.; Wolfrum, C.; Rao, T.N.; Winnay, J.N.; Garcia-Martin, R.; Grinspoon, S.K.; et al. Adipose-derived circulating miRNAs regulate gene expression in other tissues. Nature 2017, 542, 450–455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sebastiani, G.; Nigi, L.; Grieco, G.E.; Mancarella, F.; Ventriglia, G.; Dotta, F. Circulating microRNAs and diabetes mellitus: A novel tool for disease prediction, diagnosis, and staging? J. Endocrinol. Investig. 2017, 40, 591–610. [Google Scholar] [CrossRef] [PubMed]
  27. Nielsen, L.B.; Wang, C.; Sørensen, K.; Bang-Berthelsen, C.H.; Hansen, L.; Andersen, M.-L.M.; Hougaard, P.; Juul, A.; Zhang, C.-Y.; Pociot, F.; Mortensen, H.B. Circulating levels of microRNA from children with newly diagnosed type 1 diabetes and healthy controls: Evidence that miR-25 associates to residual β-cell function and glycaemic control during disease progression. Exp. Diabetes Res. 2012, 2012, 896362. [Google Scholar] [PubMed]
  28. Seyhan, A.A.; Nunez Lopez, Y.O.; Xie, H.; Yi, F.; Mathews, C.; Pasarica, M.; Pratley, R.E. Pancreas-enriched miRNAs are altered in the circulation of subjects with diabetes: A pilot cross-sectional study. Sci. Rep. 2016, 6, 31479. [Google Scholar] [CrossRef]
  29. Assmann, T.S.; Recamonde-Mendoza, M.; De Souza, B.M.; Crispim, D. MicroRNA expression profiles and type 1 diabetes mellitus: Systematic review and bioinformatic analysis. Endocr. Connect. 2017, 6, 773–790. [Google Scholar] [CrossRef]
  30. Åkerman, L.; Casas, R.; Ludvigsson, J.; Tavira, B.; Skoglund, C. Serum miRNA levels are related to glucose homeostasis and islet autoantibodies in children with high risk for type 1 diabetes. PLoS ONE 2018, 13, e0191067. [Google Scholar] [CrossRef]
  31. Erener, S.; Marwaha, A.; Tan, R.; Panagiotopoulos, C.; Kieffer, T.J. Profiling of circulating microRNAs in children with recent onset of type 1 diabetes. JCI Insight 2017, 2, e89656. [Google Scholar] [CrossRef]
  32. Samandari, N.; Mirza, A.H.; Nielsen, L.B.; Kaur, S.; Hougaard, P.; Fredheim, S.; Mortensen, H.B.; Pociot, F. Circulating microRNA levels predict residual β cell function and glycaemic control in children with type 1 diabetes mellitus. Diabetologia 2017, 60, 354–363. [Google Scholar] [CrossRef] [PubMed]
  33. Bedene, A.; Mencej Bedrač, S.; Ješe, L.; Marc, J.; Vrtačnik, P.; Preželj, J.; Kocjan, T.; Kranjc, T.; Ostanek, B. miR-148a the epigenetic regulator of bone homeostasis is increased in plasma of osteoporotic postmenopausal women. Wien. Klin. Wochenschr. 2016, 128, 519–526. [Google Scholar] [CrossRef] [PubMed]
  34. Seeliger, C.; Karpinski, K.; Haug, A.T.; Vester, H.; Schmitt, A.; Bauer, J.S.; van Griensven, M. Five freely circulating miRNAs and bone tissue miRNAs are associated with osteoporotic fractures. J. Bone Miner. Res. 2014, 29, 1718–1728. [Google Scholar] [CrossRef]
  35. Li, H.; Wang, Z.; Fu, Q.; Zhang, J. Plasma miRNA levels correlate with sensitivity to bone mineral density in postmenopausal osteoporosis patients. Biomarkers 2014, 19, 553–556. [Google Scholar] [CrossRef] [PubMed]
  36. Yang, N.; Wang, G.; Hu, C.; Shi, Y.; Liao, L.; Shi, S.; Cai, Y.; Cheng, S.; Wang, X.; Liu, Y.; et al. Tumor necrosis factor α suppresses the mesenchymal stem cell osteogenesis promoter miR-21 in estrogen deficiency-induced osteoporosis. J. Bone Miner. Res. 2013, 28, 559–573. [Google Scholar] [CrossRef] [Green Version]
  37. Zhao, C.; Sun, W.; Zhang, P.; Ling, S.; Li, Y.; Zhao, D.; Peng, J.; Wang, A.; Li, Q.; Song, J.; et al. miR-214 promotes osteoclastogenesis by targeting Pten/PI3k/Akt pathway. RNA Biol. 2015, 12, 343–353. [Google Scholar] [CrossRef] [Green Version]
  38. Du, F.; Wu, H.; Zhou, Z.; Liu, Y.U. microRNA-375 inhibits osteogenic differentiation by targeting runt-related transcription factor. Exp. Ther. Med. 2015, 10, 207–212. [Google Scholar] [CrossRef]
  39. You, L.; Pan, L.; Chen, L.; Gu, W.; Chen, J. MiR-27a is Essential for the Shift from Osteogenic Differentiation to Adipogenic Differentiation of Mesenchymal Stem Cells in Postmenopausal Osteoporosis. Cell. Physiol. Biochem. 2016, 39, 253–265. [Google Scholar] [CrossRef]
  40. Zeng, H.-C.; Bae, Y.; Dawson, B.C.; Chen, Y.; Bertin, T.; Munivez, E.; Campeau, P.M.; Tao, J.; Chen, R.; Lee, B.H. MicroRNA miR-23a cluster promotes osteocyte differentiation by regulating TGF-β signalling in osteoblasts. Nat. Commun. 2017, 8, 15000. [Google Scholar] [CrossRef]
  41. Gennari, L.; Merlotti, D.; Valenti, R.; Ceccarelli, E.; Ruvio, M.; Pietrini, M.G.; Capodarca, C.; Franci, M.B.; Campagna, M.S.; Calabrò, A.; et al. Circulating sclerostin levels and bone turnover in type 1 and type 2 diabetes. J. Clin. Endocrinol. Metab. 2012, 97, 1737–1744. [Google Scholar] [CrossRef]
  42. Shah, V.N.; Harrall, K.K.; Shah, C.S.; Gallo, T.L.; Joshee, P.; Snell-Bergeon, J.K.; Kohrt, W.M. Bone mineral density at femoral neck and lumbar spine in adults with type 1 diabetes: A meta-analysis and review of the literature. Osteoporos. Int. 2017, 28, 2601–2610. [Google Scholar] [CrossRef] [PubMed]
  43. Kaur, H.; Joshee, P.; Franquemont, S.; Baumgartner, A.; Thurston, J.; Pyle, L.; Nadeau, K.J.; Shah, V.N. Bone mineral content and bone density is lower in adolescents with type 1 diabetes: A brief report from the RESISTANT and EMERALD studies. J. Diabetes Complicat. 2018, 32, 931–933. [Google Scholar] [CrossRef] [PubMed]
  44. Cheng, P.; Chen, C.; He, H.-B.; Hu, R.; Zhou, H.-D.; Xie, H.; Zhu, W.; Dai, R.-C.; Wu, X.-P.; Liao, E.-Y.; et al. miR-148a regulates osteoclastogenesis by targeting V-maf musculoaponeurotic fibrosarcoma oncogene homolog B. J. Bone Miner. Res. 2013, 28, 1180–1190. [Google Scholar] [CrossRef]
  45. Mu, J.; Zhu, D.; Shen, Z.; Ning, S.; Liu, Y.; Chen, J.; Li, Y.; Li, Z. The repressive effect of miR-148a on Wnt/β-catenin signaling involved in Glabridin-induced anti-angiogenesis in human breast cancer cells. BMC Cancer 2017, 17, 307. [Google Scholar] [CrossRef]
  46. Shi, C.; Zhang, M.; Tong, M.; Yang, L.; Pang, L.; Chen, L.; Xu, G.; Chi, X.; Hong, Q.; Ni, Y.; et al. miR-148a is Associated with Obesity and Modulates Adipocyte Differentiation of Mesenchymal Stem Cells through Wnt Signaling. Sci. Rep. 2015, 5, 9930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Joeng, K.S.; Lee, Y.-C.; Lim, J.; Chen, Y.; Jiang, M.-M.; Munivez, E.; Ambrose, C.; Lee, B.H. Osteocyte-specific WNT1 regulates osteoblast function during bone homeostasis. J. Clin. Investig. 2017, 127, 2678–2688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Li, J. JAK-STAT and bone metabolism. Jakstat 2013, 2, e23930. [Google Scholar] [CrossRef] [PubMed]
  49. Zhou, H.; Newnum, A.B.; Martin, J.R.; Li, P.; Nelson, M.T.; Moh, A.; Fu, X.-Y.; Yokota, H.; Li, J. Osteoblast/osteocyte-specific inactivation of Stat3 decreases load-driven bone formation and accumulates reactive oxygen species. Bone 2011, 49, 404–411. [Google Scholar] [CrossRef]
  50. Rached, M.-T.; Kode, A.; Xu, L.; Yoshikawa, Y.; Paik, J.-H.; Depinho, R.A.; Kousteni, S. FoxO1 is a positive regulator of bone formation by favoring protein synthesis and resistance to oxidative stress in osteoblasts. Cell Metab. 2010, 11, 147–160. [Google Scholar] [CrossRef]
  51. Dou, C.; Li, N.; Ding, N.; Liu, C.; Yang, X.; Kang, F.; Cao, Z.; Quan, H.; Hou, T.; Xu, J.; et al. HDAC2 Regulates FOXO1 during RANKL-Induced Osteoclastogenesis. Am. J. Physiol. Cell Physiol. 2016, 310, C780–C787. [Google Scholar] [CrossRef]
  52. Iyer, S.; Ambrogini, E.; Bartell, S.M.; Han, L.; Roberson, P.K.; de Cabo, R.; Jilka, R.L.; Weinstein, R.S.; O’Brien, C.A.; Manolagas, S.C.; et al. FOXOs attenuate bone formation by suppressing Wnt signaling. J. Clin. Investig. 2013, 123, 3409–3419. [Google Scholar] [CrossRef] [Green Version]
  53. Kode, A.; Mosialou, I.; Silva, B.C.; Joshi, S.; Ferron, M.; Rached, M.T.; Kousteni, S. FoxO1 protein cooperates with ATF4 protein in osteoblasts to control glucose homeostasis. J. Biol. Chem. 2012, 287, 8757–8768. [Google Scholar] [CrossRef]
  54. Wang, Y.; Dong, G.; Jeon, H.H.; Elazizi, M.; La, L.B.; Hameedaldeen, A.; Xiao, E.; Tian, C.; Alsadun, S.; Choi, Y.; et al. FOXO1 mediates RANKL-induced osteoclast formation and activity. J. Immunol. 2015, 194, 2878–2887. [Google Scholar] [CrossRef] [PubMed]
  55. Meng, X.; Vander Ark, A.; Lee, P.; Hostetter, G.; Bhowmick, N.A.; Matrisian, L.M.; Williams, B.O.; Miranti, C.K.; Li, X. Myeloid-specific TGF-β signaling in bone promotes basic-FGF and breast cancer bone metastasis. Oncogene 2016, 35, 2370–2378. [Google Scholar] [CrossRef]
  56. Meng, X.; Vander Ark, A.; Daft, P.; Woodford, E.; Wang, J.; Madaj, Z.; Li, X. Loss of TGF-β signaling in osteoblasts increases basic-FGF and promotes prostate cancer bone metastasis. Cancer Lett. 2018, 418, 109–118. [Google Scholar] [CrossRef]
  57. Wu, M.; Chen, G.; Li, Y.-P. TGF-β and BMP signaling in osteoblast, skeletal development, and bone formation, homeostasis and disease. Bone Res. 2016, 4, 16009. [Google Scholar] [CrossRef]
  58. Lamora, A.; Talbot, J.; Mullard, M.; Brounais-Le Royer, B.; Redini, F.; Verrecchia, F. TGF-β Signaling in Bone Remodeling and Osteosarcoma Progression. J. Clin. Med. 2016, 5, 96. [Google Scholar] [CrossRef]
  59. Regazzi, R. MicroRNAs as therapeutic targets for the treatment of diabetes mellitus and its complications. Expert Opin. Ther. Targets 2018, 22, 153–160. [Google Scholar] [CrossRef]
  60. Wang, C.; Wan, S.; Yang, T.; Niu, D.; Zhang, A.; Yang, C.; Cai, J.; Wu, J.; Song, J.; Zhang, C.-Y.; et al. Increased serum microRNAs are closely associated with the presence of microvascular complications in type 2 diabetes mellitus. Sci. Rep. 2016, 6, 20032. [Google Scholar] [CrossRef] [Green Version]
  61. Weber, D.R.; Schwartz, G. Epidemiology of skeletal health in type 1 diabetes. Curr. Osteoporos. Rep. 2016, 14, 327–336. [Google Scholar] [CrossRef]
  62. Sugatani, T.; Vacher, J.; Hruska, K.A. A microRNA expression signature of osteoclastogenesis. Blood 2011, 117, 3648–3657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Hackl, M.; Heilmeier, U.; Weilner, S.; Grillari, J. Circulating microRNAs as novel biomarkers for bone diseases—Complex signatures for multifactorial diseases? Mol. Cell. Endocrinol. 2016, 432, 83–95. [Google Scholar] [CrossRef] [PubMed]
  64. Gao, J.; Yang, T.; Han, J.; Yan, K.; Qiu, X.; Zhou, Y.; Fan, Q.; Ma, B. MicroRNA expression during osteogenic differentiation of human multipotent mesenchymal stromal cells from bone marrow. J. Cell. Biochem. 2011, 112, 1844–1856. [Google Scholar] [CrossRef] [PubMed]
  65. Xu, J.-F.; Yang, G.-H.; Pan, X.-H.; Zhang, S.-J.; Zhao, C.; Qiu, B.-S.; Gu, H.-F.; Hong, J.-F.; Cao, L.; Chen, Y.; et al. Altered microRNA expression profile in exosomes during osteogenic differentiation of human bone marrow-derived mesenchymal stem cells. PLoS ONE 2014, 9, e114627. [Google Scholar] [CrossRef] [PubMed]
  66. Hu, C.-H.; Sui, B.-D.; Du, F.-Y.; Shuai, Y.; Zheng, C.-X.; Zhao, P.; Yu, X.-R.; Jin, Y. miR-21 deficiency inhibits osteoclast function and prevents bone loss in mice. Sci. Rep. 2017, 7, 43191. [Google Scholar] [CrossRef] [Green Version]
  67. Xu, Z.; Liu, X.; Wang, H.; Li, J.; Dai, L.; Li, J.; Dong, C. Lung adenocarcinoma cell-derived exosomal miR-21 facilitates osteoclastogenesis. Gene 2018, 666, 116–122. [Google Scholar] [CrossRef] [PubMed]
  68. Eguchi, T.; Watanabe, K.; Hara, E.S.; Ono, M.; Kuboki, T.; Calderwood, S.K. OstemiR: A novel panel of microRNA biomarkers in osteoblastic and osteocytic differentiation from mesencymal stem cells. PLoS ONE 2013, 8, e58796. [Google Scholar] [CrossRef]
  69. Shilo, V.; Mor-Yosef Levi, I.; Abel, R.; Mihailović, A.; Wasserman, G.; Naveh-Many, T.; Ben-Dov, I.Z. Let-7 and MicroRNA-148 Regulate Parathyroid Hormone Levels in Secondary Hyperparathyroidism. J. Am. Soc. Nephrol. 2017, 28, 2353–2363. [Google Scholar] [CrossRef]
  70. Pizzamiglio, S.; Zanutto, S.; Ciniselli, C.M.; Belfiore, A.; Bottelli, S.; Gariboldi, M.; Verderio, P. A methodological procedure for evaluating the impact of hemolysis on circulating microRNAs. Oncol. Lett. 2017, 13, 315–320. [Google Scholar] [CrossRef]
  71. Kirschner, M.B.; Edelman, J.J.B.; Kao, S.C.-H.; Vallely, M.P.; van Zandwijk, N.; Reid, G. The Impact of Hemolysis on Cell-Free microRNA Biomarkers. Front. Genet. 2013, 4, 94. [Google Scholar] [CrossRef]
  72. Murray, M.J.; Watson, H.L.; Ward, D.M.; Bailey, S.; Ferraresso, M.; Nicholson, J.C.; Gnanapragasam, V.J.; Thomas, B.; Scarpini, C.G.; Coleman, N. “Future-proofing” blood processing for measurement of circulating microRNAs in samples from biobanks and prospective clinical trials. Cancer Epidemiol. Biomarkers Prev. 2017, 27, 208–218. [Google Scholar] [CrossRef]
  73. Murray, M.J.; Bell, E.; Raby, K.L.; Rijlaarsdam, M.A.; Gillis, A.J.M.; Looijenga, L.H.J.; Brown, H.; Destenaves, B.; Nicholson, J.C.; Coleman, N. A pipeline to quantify serum and cerebrospinal fluid microRNAs for diagnosis and detection of relapse in paediatric malignant germ-cell tumours. Br. J. Cancer 2015, 114, 151–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Wang, J.; Zhu, Y.; Jin, F.; Tang, L.; He, Z.; He, Z. Differential expression of circulating microRNAs in blood and haematoma samples from patients with intracerebral haemorrhage. J. Int. Med. Res. 2016, 44, 419–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Rice, J.; Roberts, H.; Rai, S.N.; Galandiuk, S. Housekeeping genes for studies of plasma microRNA: A need for more precise standardization. Surgery 2015, 158, 1345–1351. [Google Scholar] [CrossRef]
  76. Marabita, F.; de Candia, P.; Torri, A.; Tegnér, J.; Abrignani, S.; Rossi, R.L. Normalization of circulating microRNA expression data obtained by quantitative real-time RT-PCR. Brief Bioinform. 2016, 17, 204–212. [Google Scholar] [CrossRef] [PubMed]
  77. Gevaert, A.B.; Witvrouwen, I.; Vrints, C.J.; Heidbuchel, H.; Van Craenenbroeck, E.M.; Van Laere, S.J.; Van Craenenbroeck, A.H. MicroRNA profiling in plasma samples using qPCR arrays: Recommendations for correct analysis and interpretation. PLoS ONE 2018, 13, e0193173. [Google Scholar] [CrossRef] [PubMed]
Ncrna 04 00037 g001 550
Figure 1. The expression of hsa-miR-148a and miR-21-5p is increased in the serum of patients with T1D. Single assay RT-qPCR validation of miR-21-5p (a), miR-24 (b), miR-27a (c), miR-148a (d), miR-214 (e), and miR-375 (f) in n = 14 non-diabetic and n = 15 T1D patients. Data are reported as mean ± SD of normalized 2−ΔCT values. Statistics using Mann–Whitney U test, p < 0.05.
Figure 1. The expression of hsa-miR-148a and miR-21-5p is increased in the serum of patients with T1D. Single assay RT-qPCR validation of miR-21-5p (a), miR-24 (b), miR-27a (c), miR-148a (d), miR-214 (e), and miR-375 (f) in n = 14 non-diabetic and n = 15 T1D patients. Data are reported as mean ± SD of normalized 2−ΔCT values. Statistics using Mann–Whitney U test, p < 0.05.
Ncrna 04 00037 g001
Ncrna 04 00037 g002 550
Figure 2. The expression of hsa-miR-148a and miR-21-5p levels are correlated with bone metabolism parameters in T1D patients and non-diabetic control subjects. Correlation analysis between miR-148a serum expression levels, reported as normalized 2−ΔCT values, and BMD total body (TB) reported as g/cm2 (a) and circulating levels of parathormone (PTH) reported as pg/mL (b). Correlation analysis between miR-21-5p serum expression levels and BMC-FEM reported as g/cm2 (c). Spearman R test was performed to evaluate r-values and p-values (p < 0.05).
Figure 2. The expression of hsa-miR-148a and miR-21-5p levels are correlated with bone metabolism parameters in T1D patients and non-diabetic control subjects. Correlation analysis between miR-148a serum expression levels, reported as normalized 2−ΔCT values, and BMD total body (TB) reported as g/cm2 (a) and circulating levels of parathormone (PTH) reported as pg/mL (b). Correlation analysis between miR-21-5p serum expression levels and BMC-FEM reported as g/cm2 (c). Spearman R test was performed to evaluate r-values and p-values (p < 0.05).
Ncrna 04 00037 g002
Ncrna 04 00037 g003 550
Figure 3. miR-148a and miR-21-5p regulate common pathways involved in bone metabolism and remodeling. Hierarchical Clustering Heatmap bioinformatic analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of hsa-miR-148a-3p and hsa-miR-21-5p target genes, showed FoxO and TGF-β signaling pathway as main pathways involved in bone metabolism and remodeling. Of note, FoxO signaling pathway is commonly identified both in miR-148a-3p and miR-21-5p analysis. Color key indicates Log p values from the less significant pathway (light yellow) to the most significant pathway (dark red).
Figure 3. miR-148a and miR-21-5p regulate common pathways involved in bone metabolism and remodeling. Hierarchical Clustering Heatmap bioinformatic analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of hsa-miR-148a-3p and hsa-miR-21-5p target genes, showed FoxO and TGF-β signaling pathway as main pathways involved in bone metabolism and remodeling. Of note, FoxO signaling pathway is commonly identified both in miR-148a-3p and miR-21-5p analysis. Color key indicates Log p values from the less significant pathway (light yellow) to the most significant pathway (dark red).
Ncrna 04 00037 g003
Table
Table 1. Selected miRNAs analyzed in the study.
Table 1. Selected miRNAs analyzed in the study.
Selected miRNAsEvidences of Circulating miRNAs Differentially Expressed in T1D PatientsEvidences of miRNAs Involved in Bone Remodeling
hsa-miR-148aNielsen et al. 2012 [27]
Seyan et al.2016 [28]
Assman et al. 2017 [29]
Akerman 2018 [30]		Bedene et al. 2016 [33]
Seeliger et al. 2015 [34]
hsa-miR-21-5pSeyan et al.2016 [28]
Assman et al. 2017 [29]		Li et al. 2014 [35]
Seeliger et al. 2015 [34]
Yang et al. 2013 [36]
hsa-miR-214Erener et al. 2017 [31]Zhao et al. 2015 [37]
hsa-miR-375Seyan et al.2016 [28]Du et al. 2015 [38]
Samandari et al.2016 [32]
Assman et al. 2017 [29]
hsa-miR-27aNielsen et al.2012 [27]You et al. 2016 [39]
Zeng et al. 2017 [40]
hsa-miR-24Seyan et al.2016 [28]Seeliger et al. 2015 [34]
Samandari et al.2016 [32]
Assman et al. 2017 [29]
Nielsen et al. 2012 [27]
Erener et al. 2017 [31]
Table
Table 2. Main clinical characteristics of non-diabetic subjects and type 1 diabetic patients
Table 2. Main clinical characteristics of non-diabetic subjects and type 1 diabetic patients
Non-Diabetic SubjectsType 1 Diabetic Patients
Subjects (n)1415
Male/Female (n/n)10/410/5
Age (years)28 ± 7.9332 ± 7.9
Disease duration (years)/15.71 ± 11.33
HbA1c (%)N/A7.36 ± 0.80593
Chronic complications (n/n)N/A2/15
Smoking habits (n/n)N/A3/15
Other therapies (n/n)N/A6/15
Chronic inflammatory diseases (n/n)N/A4/15
BMI #24 ± 4.3421 ± 1.73 #
Ca (mg/dL)9.37 ± 0.509.54 ± 0.46
Osteocalcin (ng/mL)4.24 ± 1.824.50 ± 2.21
bALP (μg/L)10.5 ± 7.3310.2 ± 4.72
25OHD (ng/mL) #27.9 ± 10.3314.3 ± 9.83 #
PTH (pg/mL) #18 ± 9.3626 ± 10.23 #
BMD TB (g/cm2)1.23 ± 0.111.15 ± 0.09
T score TB0.7 ± 1.18−0.2 ± 1.15
Z score TB0.43 ± 1.100.08 ± 1.12
BMC TB (g/cm2)2.82 ± 0.462.61 ± 0.38
BMD FEM (g/cm2) #1.156 ± 0.170.967 ± 0.13 #
T score FEM #0.8 ± 1.18−0.8 ± 1.10 #
Z score FEM #0.6 ± 1.04−0.8 ± 1.11 #
BMC FEM (g/cm2) *37.86 ± 8.1332.04 ± 4.82*
BMD N (g/cm2) *1.09 ± 0.180.94 ± 0.12*
T score N *0.3 ± 1.27−0.80 ± 1.08*
Z score N #0.2 ± 1.13−0.8 ± 1.10 #
BMC N (g/cm2)5.34 ± 1.354.40 ± 1.04
BMD L (g/cm2)1.20 ± 0.151.14 ± 0.17
T score L−0.04 ± 1.24−0.57 ± 1.48
Z score L−0.14 ± 1.26−0.42 ± 1.59
Note: BMI= Body Mass Index; Ca = Calcium; bALP = Bone Alkaline Phosphatase; 25 OHD = 25 OH Vitamin D; PTH = Parathormone; BMD = Bone Mineral Density; BMC= Bone Mineral Content; TB = Total Body; FEM= Femur; N = Neck; L = Lumbar. Samples distribution analysis was performed for all parameters by using D’Agostino & Pearson normality test. BMI, bALP, 25 OHD, PTH, T score TB, BMD FEM, T score FEM, Z score FEM and Z score N are indicated as median ± SD (# Mann-Whitney U test = p ≤ 0.05); age, Ca, osteocalcin, BMD TB, Z score TB, BMC TB, BMC FEM, BMD N, T score N, BMC N, BMD L, T score L, Z score L are indicated as mean ± SD. (* Unpaired t test = p ≤ 0.05).
Table
Table 3. Validated target genes of hsa-miR-148a and hsa-miR-21-5p.
Table 3. Validated target genes of hsa-miR-148a and hsa-miR-21-5p.
Target GeneGene NameFunctionTargeting miRNA
MAFBV-mafavianmusculoaponeurotic fibrosarcoma oncogenehomolog BNegatively regulates RANKL-induced osteoclast differentiationmiR-148a
WNT1Wingless-type MMTV integration site family, member 1Canonical pathwayWnt/β catenin: stimulates osteoblast differentiation and proliferationmiR-148a
BCL2L11BCL2-like 11 (apoptosis facilitator)Induces apoptosismiR-148a
PTENPhosphatase and tensin homologInteracts with Canonical pathwayWnt/β cateninmiR-148a
TGFB2Transforming growth factor, β 2Stimulates osteoblast differentiation and proliferation and bone formationmiR-148a
IGF1Insulin-like growth factor 1 (somatomedin C)Stimulates bone mineralization and positively regulates osteoblast differentiationmiR-148a
GADD45AGrowth arrest and DNA-damage-inducible, αPositively regulates apoptosismiR-148a
STAT3Signal transducer and activator of transcription 3Specifically inactivates osteoblast/osteocytesmiR-21-5p
TGFBR1Transforming Growth Factor β Receptor 1Regulates bone resorptionmiR-21-5p
IGF1RInsulin-like growth factor 1 receptorInvolved in bone remodelingmiR-21-5p
BCL6B-Cell lymphoma 6Regulates osteoclast differentiationmiR-21-5p
TGFB2Transforming growth factor-β 2Involved in osteoblast differentiation and bone remodelingmiR-21-5p
PDCD4Programmed Cell Death 4Involved in osteoclastogenesis through cFos regulationmiR-21-5p
SOD2Superoxide dismutase 2Controls osteoclastogenesis and maintains osteoblast differentiationmiR-21-5p
FOXO3Forkhead box O-3Its overexpression increases bone resorptionmiR-21-5p
USP7Ubiquitin-specific protease 7Involved in osteogenic differentiationmiR-21-5p
PTENPhosphatase and tensin homologInvolved in adiponectin-induced osteogenesismiR-21-5p
IL-10Interleukin 10Inhibits bone resorptionmiR-21-5p

Share and Cite

MDPI and ACS Style

Grieco, G.E.; Cataldo, D.; Ceccarelli, E.; Nigi, L.; Catalano, G.; Brusco, N.; Mancarella, F.; Ventriglia, G.; Fondelli, C.; Guarino, E.; et al. Serum Levels of miR-148a and miR-21-5p Are Increased in Type 1 Diabetic Patients and Correlated with Markers of Bone Strength and Metabolism. Non-Coding RNA 2018, 4, 37. https://0-doi-org.brum.beds.ac.uk/10.3390/ncrna4040037

AMA Style

Grieco GE, Cataldo D, Ceccarelli E, Nigi L, Catalano G, Brusco N, Mancarella F, Ventriglia G, Fondelli C, Guarino E, et al. Serum Levels of miR-148a and miR-21-5p Are Increased in Type 1 Diabetic Patients and Correlated with Markers of Bone Strength and Metabolism. Non-Coding RNA. 2018; 4(4):37. https://0-doi-org.brum.beds.ac.uk/10.3390/ncrna4040037

Chicago/Turabian Style

Grieco, Giuseppina E., Dorica Cataldo, Elena Ceccarelli, Laura Nigi, Giovanna Catalano, Noemi Brusco, Francesca Mancarella, Giuliana Ventriglia, Cecilia Fondelli, Elisa Guarino, and et al. 2018. "Serum Levels of miR-148a and miR-21-5p Are Increased in Type 1 Diabetic Patients and Correlated with Markers of Bone Strength and Metabolism" Non-Coding RNA 4, no. 4: 37. https://0-doi-org.brum.beds.ac.uk/10.3390/ncrna4040037

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