1. Introduction
Type 1 Diabetes Mellitus (T1DM) is a chronic autoimmune disease characterized by the degeneration of beta cells, resulting in insulin deficiency [
1]. T1DM incidence is strongly associated with Vitamin D deficiency, and 1,25(OH)2D3 can potentially prevent islet cell death and enhance insulin production. Several studies suggest that low 1,25(OH)2D3 levels may affect beta cell regulation [
2]. Numerous epidemiological studies have reported a high prevalence of Vitamin D deficiency among Saudi children with T1DM, indicating a strong correlation between these two factors [
3]
The Vitamin D metabolism pathway is controlled by various genes, including
CYP2R1,
CYP27B1 (responsible for activation), 24-hydroxylase (
CYP24A1) (responsible for inactivation), VDR (responsible for action), and GC. The concentration of 1,25(OH)2D3, the physiologically active form of Vitamin D, is tightly regulated by both 1-hydroxylase and the catabolic enzyme 24-hydroxylase (encoded by the
CYP24A1 gene).
CYP24A1 catalyzes the hydroxylation reaction leading to the degradation of 1,25(OH)2D3, resulting in the excretion of calcitroic acid and other metabolites in the bile.
CYP24A1 expression is induced by both 1,25(OH)2D3 and 25(OH)D3, making it one of the most highly inducible genes in humans, capable of increasing its transcription by 20,000-fold [
4].
Vitamin D deficiency might enhance autoimmune responses in the context of T1DM. 1,25(OH)2D3 can exert its immunomodulatory effects through genomic responses and its ability to alter gene transcription, as most human cells contain Vitamin D receptors (VDR) [
5]. Numerous immune cells express
VDRs and
CYP27B1 enzymes, the synthesis of which is influenced by various immune-specific signals [
6]. In autoimmune diseases, the Vitamin D metabolite plays a crucial role in downregulating all mechanisms associated with adaptive immunity, inducing immunological tolerance, and promoting anti-inflammatory activity [
7]. Therefore, understanding the molecular mechanism responsible for Vitamin D deficiency in T1DM patients is paramount.
miRNAs are small, non-coding RNA molecules that post-transcriptionally regulate gene expression, often by preventing or triggering the degradation of mRNAs [
8]. miRNAs that target genes involved in Vitamin D metabolism can impact the circulating 1,25(OH)2D3 levels in individuals with T1DM. Several miRNAs have been predicted to target
CYP24A1, but only has-miR-125b-5p has been experimentally validated [
9].Studies conducted on ovarian granulosa and breast cancer cells have shown that overexpression or antisense knockdown of has-miR-125b-5p suppresses and enhances the expression of
CYP24A1 protein, respectively [
9]. Besides its role in Vitamin D catabolism, has-miR-125b-5p also targets the
VDR. Reducing has-miR-125b-5p expression enhances the response to 1,25(OH)2D3 in melanoma cell lines [
10]. Conversely, increased has-miR-125b-5p levels suppress endogenous
VDR protein levels in MCF-7 breast cancer cells, contributing to resistance to 1,25(OH)2D [
11]. While other miRNAs may potentially target VDR, only one, miR-326, has been validated in the peripheral blood lymphocytes of individuals with T1DM [
12]. Considering the proposed role of 1,25(OH)2D3 in protecting against autoimmunity [
13], it is reasonable to speculate that miR-326 may inhibit the immunomodulatory effects of 1,25(OH)2D3 in preventing inflammatory and autoimmune disorders. Moreover, studies have shown that miR-125a-5p is upregulated in Treg cells isolated from the lymph nodes draining the pancreas of T1DM patients. This upregulation potentially contributes to the reduced expression of its target gene, C-C Chemokine Receptor type-2 (
CCR2)[
14].
The Vitamin D -activating enzymes,
CYP27B1 is predicted to be targeted by multiple miRNAs, although only one has been experimentally validated. Several studies have shown that miR-21 inhibits
CYP27B1 expression in monocytes infected with Mycobacterium leprae (
M. Leprae), thus inhibiting downstream antibacterial responses induced by Vitamin D intracrine signaling [
15]. Another miRNA predicted to target
CYP27B1 is has-miR-216b-5p. T1DM is associated with increased miR-216a expression in pancreatic islets, possibly as a compensatory mechanism [
16]. Both miR-377 and miR-216a have been identified as early biomarkers of nephropathy in children with T1DM. Their correlation with carotid intimal thickness (
CIMT) provides insights into the subclinical atherosclerotic processes in diabetic nephropathy [
17].
Accumulating evidence suggests that miRNAs are involved in T1DM pathogenesis through multiple mechanisms, including the regulation of immune cell differentiation, development, activation, and the disruption of immune system equilibrium. Overexpression of miR-34a in diabetic mice reduces B lymphopoiesis capacity, disturbs pancreatic islet defense, and increases sensitivity to damage [
18]. Differential miRNA expression also influences the production of specific T lymphocytes. Moreover, miR-26 and miR-101 have been associated with the differentiation of cells toward the T helper 1 (Th1) phenotype. Additionally, miR-21, miR-93, miR-326, and miR-31 are believed to alter T cell functions and play a role in T1DM autoimmunity [
19,
20].
Several miRNAs have been associated with T1DM in peripheral blood mononuclear cells (PBMCs). Patients with recently diagnosed T1DM exhibited significant downregulation of miR-21a and miR-93, which target NF-KB signaling to regulate apoptosis and inflammation [
21]. Furthermore, miR-326 is overexpressed in PBMCs from T1DM patients [
22,
23], suggesting the involvement of miRNAs in T1DM autoimmunity as they target significant immune modulators—VDR and erythroblastosis virus E26 oncogenic homolog 1 (ETS-1). Additionally, miRNA signatures in PBMCs were correlated with autoantibodies in T1DM patients, with increased miR-326 levels correlated with antibodies against glutamic acid decarboxylase (GAD) and tyrosine phosphatase-like protein (IA2) [
22,
23], and reduced miR-146a levels correlated with antibodies against GAD [
21]. Detecting miRNAs in PBMCs offers the advantage of using them as biomarkers for monitoring disease progression.
The cause of T1DM is not fully understood, but is believed to be the development of autoantibodies and autoreactive Th1 and cytotoxic T lymphocytes (CTLs), which cause the immune system to destroy insulin-producing pancreatic cells [
1,
24]. Vitamin D deficiency appears to contribute to increased activation of B cells and autoantibody production, and long-term supplementation with Vitamin D led to an increase in T-regulator cells in individuals with SLE [
25,
26].
Here, we hypothesize that irregular miRNA expression in PBMCs from patients with T1DM might influence the downstream target genes involved in Vitamin D metabolism. This could potentially play a significant role in the development of T1DM and 1,25(OH)2D3 insufficiency. Additionally, circulating miRNAs and Vitamin D-related genes may affect the pathogenesis of T1DM by modulating autoimmune responses.
In this study, we aimed to ascertain whether there are differential expressions of Vitamin D metabolism genes between patients with T1DM and healthy control. Additionally, we sought to identify any potential correlations between circulating serum 1,25(OH)2D3 concentrations and the expression levels of these metabolic genes in individuals with T1DM. We further investigated the potential influence of circulating miRNAs on the expression of vitamin D metabolism genes within the PBMCs of those with T1DM. Notably, alterations in vitamin D metabolism genes and associated miRNAs could hold significance in relation to the autoimmune aspect of T1DM.
3. Discussion
Vitamin D is pivotal in modulating autoimmune diseases, including diabetes. Its deficiency has been implicated in the pathogenesis of T1DM. In this context, we analyzed whether serum 1,25(OH)2D3 concentrations and Vitamin D metabolic genes' expression correlate in T1DM patients. Additionally, we investigated if miRNAs might modulate the expression of these Vitamin D metabolism genes in PBMCs of T1DM, given that miRNAs targeting genes in the Vitamin D metabolism pathway could influence circulating 1,25(OH)2D3 levels in T1DM. To the best of our knowledge, this is the first study linking circulating miRNA expression with alterations in Vitamin D metabolism genes within PBMCs of T1DM, considering Vitamin D status and ongoing islet autoimmunity.
There is a robust association between T1DM incidence and Vitamin D deficiency [
27]. Many epidemiological studies have underscored the prevalence of Vitamin D deficiency among Saudi children diagnosed with T1DM, alluding to a potent connection between the two [
3]. Our research findings revealed a heightened 1,25(OH)2D3 deficiency among T1DM patients compared to non-diabetic controls, with 53.33% of T1DM patients exhibiting reduced 1,25(OH)2D3 levels. This is congruent with prior studies highlighting the pervasive nature of 1,25(OH)2D3 deficiency among T1DM patients [
28]. The relationship between Vitamin D deficiency and T1DM onset—whether causative or consequential—remains elusive.
Our investigations into a potential link between Vitamin D deficiency in T1DM and Vitamin D metabolism genes led us to compare gene expressions in T1DM patients and healthy individuals. We noted a significant decline in
CYP27B1 mRNA levels among T1DM patients, with no observable changes for
CYP24A1. Our results are coherent with previous findings showing attenuated
CYP27B1 mRNA levels in T1DM patients relative to controls. The implications of the
CYP27B1 gene on T1DM pathogenesis are evident as it modulates mRNA expression and influences 1,25(OH)2D3 serum levels, potentially via the −1260 C/A polymorphism [
29]. However, our data did not showcase a direct correlation between
CYP27B1 mRNA levels and serum 1,25(OH)2D3. There is speculation that a correlation might exist at the protein level. Further research is needed to uncover these connections, suggesting the possibility of significant breakthroughs ahead.
Given that
CYP7B1 expression didn't correlate with 1,25(OH)2D3 status in T1DM, we explored other factors potentially influencing its expression. Prior research suggests epigenetic modifications, such as miRNAs, might impact
CYP27B1 gene expression. Given the significant role of miRNAs in T1DM pathogenesis and
CYP27B1 expression regulation in other tissues [
30], we probed the potential of miRNA in modifying
CYP27B1 expression. Bioinformatic predictions posited that has-miRNA-216b-5p and has-miRNA-21-5p could target the 3' UTR of
CYP27B1, while has-miRNA-125b-5p might target
CYP24A1's 3' UTR. The expressions of influencing miRNAs, governed by factors like BMI and age, were consistent across both groups in our study. Notably, our findings indicated an upregulation of has-miR-216b-5p and a downregulation of has-miR-21-5p, with no alterations for has-miR-125b-5p. Prior research corroborates the differential expression of has-miR-21-5p and has-miR-216b-5p in T1DM [
16,
21,
31,
32]. Interestingly, while our results showcased differential expression of specific miRNAs, consistent with previous studies, they did not directly correlate with
CYP27B1 in PBMCs of T1DM. This raises the tantalizing possibility of other underlying mechanisms, perhaps genetic polymorphism or insulin resistance, playing a role in modulating
CYP27B1 expression.
While the crucial role of miRNAs in preserving immune equilibrium is apparent, the perturbed expression of specific miRNAs and its subsequent contribution to autoimmune diseases introduces new avenues of exploration [
33]. The accumulating evidence pointing towards miRNAs' involvement in T1DM pathogenesis—particularly in domains like immune responses and beta-cell metabolism—poses thought-provoking questions [
34]. Might other unidentified miRNAs be just as critical, if not more so, in influencing T1DM onset? The prevailing research, emphasizing miRNAs' role in immune stability and the negative impact of their anomalous expression on autoimmunity, raises queries about potential therapeutic interventions harnessing miRNA regulation [
33]–[
35]. With multiple studies documenting shifts in miRNA expression linked to T1DM pathology, one cannot help but speculate if reversing these alterations could mitigate disease progression [
21] [
22,
23]. Remarkably, the diminished expression of miR-21 and miR-93 in T1DM patient PBMCs, vital to inflammatory and apoptotic pathways, hints at potential therapeutic targets [
21]. Given the consistent reduction of has-miR-21-5p in our study and others, can enhancing its expression offer therapeutic benefits? Further comprehensive studies might unravel the potential opportunities within miRNA regulation in T1DM.
Several pivotal antibodies can be detected in T1DM patients' years before the onset of the disease, notably antibodies against IA-2, IA-2b, and GAD. Intriguingly, research suggests that miRNAs may play a crucial role in the biosynthesis of specific auto-antibodies, as a cluster of 32 miRNAs has been found to modulate the expression pattern of these T1DM-related auto-antibodies [
36]. However, our findings reveal no correlation between miRNAs and auto-antibodies.
Subsequently, we sought to understand the potential role of reduced
CYP27B1 expression on the presence of autoantibodies in patients with T1DM. This investigation entailed correlating
CYP27B1 mRNA levels with specific autoantibodies, notably IA2 and GAD65. We discovered a negative association between
CYP27B1 mRNA levels and IA-2 autoantibodies, suggesting a compelling possibility: the
CYP27B1 gene might play a pivotal role in the autoimmune cascade characteristic of T1DM. Moreover, previous studies have revealed associations between polymorphisms in codon 416 of the Vitamin D-binding protein gene and autoimmune markers of T1DM, particularly IA2 antibodies [
37]. Given these findings, it raises the question: Could the intricate interplay of genetic factors, miRNAs, and Vitamin D metabolism be fundamental to the complex etiology of T1DM? Future research endeavors are poised to unravel this intricate network of interactions and illuminate the underlying mechanisms of the disease.
The clinical implications of this study underscore the profound significance of 1,25(OH)2D3 levels and related metabolic gene expressions such as CYP27B1 in the early detection and monitoring of the disease. Furthermore, altering miRNA expression in T1DM may provide promising targets for future pharmaceutical interventions, enhancing the overall management and treatment avenues. Additionally, as the study reveals correlations between CYP27B1 expressions and T1DM, it emphasizes the potential of personalized medicine, wherein treatments could be tailored to the unique genetic profiles of patients. This could revolutionize how T1DM is managed, shifting from generic treatments to more individual-centric ones. Clinicians might be encouraged to adopt a more holistic approach to managing T1DM, considering an array of factors beyond just the traditional metrics of glucose and insulin.
The study presented here offers several commendable strengths, making it a valuable contribution to the current scientific understanding. Firstly, its pioneering nature cannot be understated; this research appears to be the initial foray into connecting circulating miRNA expression with alterations in the Vitamin D metabolism genes, specifically in PBMCs of T1DM, thereby potentially bridging an existing knowledge gap. This comprehensive study approach stands out, as it delves into the surface-level investigation of serum 1,25(OH)2D3 levels and ventures deeper to uncover the underlying genetic and epigenetic mechanisms that might influence these levels in T1DM patients. Furthermore, the depth of analysis is worth noting. The research did not merely stop at identifying genes with altered expressions but extended its exploration to establish correlations with other physiological markers, such as autoantibodies. Considering the burgeoning interest in the role of Vitamin D in T1DM and the increasing significance attached to miRNAs in various disease processes, this study is timely and relevant.
However, like all scientific endeavors, this study is not devoid of limitations. A significant constraint is that the gene expression was gauged only at the mRNA level. While this offers a snapshot, it might not always mirror the protein levels or biological activity, which are critical for a complete understanding. The sample size, being limited, raises concerns about the robustness and generalizability of the findings. Moreover, while the study adeptly uncovers various correlations, the age-old causality conundrum remains. It is yet to be definitively ascertained whether Vitamin D deficiency is a consequence or a precursor to T1DM. Not to mention, there might be lurking variables.
In light of the above, it becomes evident that while this study paves the way, there is a long journey ahead. More extensive cohort studies, augmented with in vitro or animal model experiments, must validate and expand on these preliminary findings. Evaluations at the protein level, coupled with functional assays, might be the key to unlocking a more profound mechanistic understanding of the processes at play.
4. Materials and Methods
4.1. Ethics and Consent
This research project underwent a rigorous review and received approval from the Institutional Review Board at the General Directorate of Health Affairs (GDHA) in Madinah (IRB Number: 276). The GDHA, an official government ethics committee, issues ethical approvals for human research on behalf of the Ministry of Health (MOH). All studies involving human participants were conducted in compliance with the ethical guidelines set forth in the Declaration of Helsinki. Every participant provided written informed consent to participate in the study.
4.2. Study Design, Recruitment of T1DM Patients, Healthy Control, Inclusion and Exclusion Criteria
In this case-control study, patients with T1DM were evaluated for symptoms of diabetes and a casual plasma glucose concentration of 7.0 mmol/L or a 2-hour post-load glucose concentration of 11.1 mmol/L during an oral glucose tolerance test, based on the criteria set by the American Diabetes Society. Age-matched control children were selected from the Maternity and Children's Hospital.
Inclusion criteria for cases included an age group of 0 to 16 years and diagnosed cases of T1DM with classical symptoms (polyuria, polydipsia, and polyphagia) along with a random plasma glucose of ≥ 200 mg/dl. Patients with pancreatic disease, hepatic disease, renal disease, bone diseases, malignancy, and any history of drug use such as calcium and Vitamin D were excluded from the cases.
For healthy controls, individuals with blood sugar within normal limits, normal growth and puberty range, and no endocrine abnormalities or autoimmune conditions were included. Subjects with pancreatic disease, hepatic disease, renal disease, bone disease, malignancy, and any history of drug use such as calcium and Vitamin D were excluded from the control group.
Data, including demographic data, medical history, comorbidities, and medications, were collected for each patient from the medical records at the Maternity and Children's Hospital.
4.3. Isolation of Peripheral Blood Mononuclear Cells (PBMCs)
Blood samples (5 ml) were collected from both T1DM patients and healthy controls into tubes containing EDTA and were diluted 1:1 in PBS. Following this, density gradient separation of the diluted samples was carried out using Hypaque-Ficoll medium (Innotrain, Germany). PBMC samples were isolated through centrifugation of the Hypaque-Ficoll gradient, and fresh PBMCs were immediately utilized for RNA extraction.
4.4. Total RNA and miRNA Extraction from PBMCs
Total RNA, including miRNA, was extracted from the purified PBMCs using the Qiagen miRNeasy Mini Kit (Qiagen, Valencia, CA, USA), following the manufacturer's protocol. The purity and concentration of RNA were determined by obtaining OD260/280 readings using a dual-beam UV spectrophotometer (Eppendorf AG, Hamburg, Germany).
4.5. Reverse Transcription and Quantitative Real-Time PCR
Complementary DNA synthesis from miRNA samples was carried out using the miScript II RT Kit and HiSpec Buffer (Qiagen), following the manufacturer’s instructions. Quantitative Real-Time PCR (qRT-PCR) was conducted using the miScript SYBR Green PCR Kit (Qiagen). The 25 μl reaction included 1 μl of cDNA, forward and reverse primers at optimized concentrations, and RNase-free water to make up the volume. The PCR was set up in a 96-well plate with triplicates, under a UV-irradiated hood on an ABI 7700HT PCR machine (Applied Biosystems 7500 Fast Real-Time PCR System) (Applied Biosystems, USA). Minus RT and no-template controls were included for each assay. The cycling conditions were: 50ºC for 2 min, 95ºC for 10 min followed by 40 cycles of 95ºC for 15 sec and 60ºC for 60 sec, and a melt curve stage (95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s). All primer sequences for the selected miRNAs and Vitamin D-related genes for real-time RT-PCR analyses are listed in
Table S2.
Results from the real-time PCR were computed using the equation [2ˆ(-∆Ct)]. The expression levels of CYP27B1 and CYP24A1 mRNA were normalized to the housekeeping gene, acidic ribosomal phosphoprotein P0 (RPLP0). For miRNAs, hsa-miR-U6B was used as a housekeeping gene to normalize all Ct values. Ct values were determined using the SDS software v.2.1, with manual baseline settings applied to fix the same threshold for both target miRNAs and reference genes
4.6. Serum Sample Collection and 1,25(OH)2D3 Level Measurement
A 3 ml blood sample was centrifuged at 1500g for 10 minutes to obtain serum samples from both patients and controls. These serum samples were then frozen at -80°C until required for analysis. Serum 1,25(OH)2D3 levels were measured using a Vitamin D Kit (IDS, UK) in a consistent laboratory setting. A concentration of 20 ng/ml or higher was considered a normal level of 1,25(OH)2D3. Concentrations between 10 and 20 ng/ml were classified as 1,25(OH)2D3 insufficiency, while levels below 10 ng/ml were categorized as 1,25(OH)2D3 deficiency.
4.7. Serological Analysis
Anti-GAD65 and anti-IA2 antibodies were assessed using enzyme immunoassay (ELISA) with Medizym commercial kits (Berlin, Germany). The detection of antibodies was performed semi-quantitatively, referencing a value of 5 IU/ml for GAD65 and 10 IU/ml for IA2. HbA1c levels were determined using a commercially available automatic system (DCA 2000, Bayer Diagnostics, Tarrytown, NY, USA).
4.8. Target Gene Identification and Bioinformatic Analysis
A comprehensive literature search was undertaken to identify a particular set of miRNAs whose expression has been previously linked with T1DM and are thought to influence the Vitamin D signaling pathway. Our study concentrated on three miRNAs, namely hsa-miR-125b-5p [
31,
32,
38], has-miR-216b-5p [
16,
17],and has-miR-21-5p [
39,
40], which have been significantly dysregulated in serum, plasma, and PBMCs of T1DM patients, as identified using bioinformatics tools (TargetScan Human, available at
http://www.targetscan.org). Based on the extant literature encompassing in vitro studies, we selected miRNAs potentially interacting with the 3′UTR sequences of
CYP24A1 and
CYP27B1 (refer to Supplementary table S1). Among several miRNAs that might interact with the
CYP24A1 3′UTR sequence, we chose hsa-miR-125b-5p. For the
CYP27B1 3′UTR, hsa-miR-216b-5p and hsa-miR-21 were selected. The involvement of these miRNAs in regulating
CYP27B1 expression was substantiated by in vitro studies.
4.9. Statistical Analysis
Statistical analyses were conducted using GraphPad Prism 6.0 software (GraphPad Software, Inc. CA, USA). Clinical parameters are represented as mean ± standard deviation (SD). To compare statistical significance between the two groups (T1DM vs. healthy control), the unpaired t-test or Mann–Whitney U test was employed. Expression of miRNA or Vitamin D-related genes is denoted as mean ± standard error. To identify correlations between miRNA expression and mRNA of Vitamin D-related genes in all participants, Pearson or Spearman correlation testing was utilized. A two-tailed p-value (p) < 0.05 was deemed statistically significant for all statistical tests.