Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 24 May 2024
Sec. Diabetes: Molecular Mechanisms
This article is part of the Research Topic Digital Technology in the Management and Prevention of Diabetes View all 11 articles

Comprehensive machine learning models for predicting therapeutic targets in type 2 diabetes utilizing molecular and biochemical features in rats

Marwa Matboli*Marwa Matboli1*Hiba S. Al-AmodiHiba S. Al-Amodi2Abdelrahman KhaledAbdelrahman Khaled3Radwa Khaled,Radwa Khaled4,5Marian M. S. RoushdyMarian M. S. Roushdy1Marwa AliMarwa Ali1Gouda Ibrahim DiabGouda Ibrahim Diab6Mahmoud Fawzy ElnagarMahmoud Fawzy Elnagar7Rasha A. Elmansy,Rasha A. Elmansy8,9Hagir H. TAhmedHagir H. TAhmed10Enshrah M. E. AhmedEnshrah M. E. Ahmed11Doaa M. A. ElzoghbyDoaa M. A. Elzoghby12Hala F. M.Kamel,Hala F. M.Kamel1,2Mohamed F. FaragMohamed F. Farag13Hind A. ELsawiHind A. ELsawi14Laila M. FaridLaila M. Farid15Mariam B. AbouelkhairMariam B. Abouelkhair15Eman K. HabibEman K. Habib16Heba FikryHeba Fikry17Lobna A. SalehLobna A. Saleh18Ibrahim H. Aboughaleb*Ibrahim H. Aboughaleb6*
  • 1Medical Biochemistry and Molecular Biology Department, Faculty of Medicine, Ain Shams University, Cairo, Egypt
  • 2Biochemistry Department, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
  • 3Bioinformatics Group, Center of Informatics Sciences (CIS), School of Information Technology and Computer Sciences, Nile University, Giza, Egypt
  • 4Biotechnology/Biomolecular Chemistry Department, Faculty of Science, Cairo University, Cairo, Egypt
  • 5Medicinal Biochemistry and Molecular Biology Department, Modern University for Technology and Information, Cairo, Egypt
  • 6Biomedical Engineering Department, Egyptian Armed Forces, Cairo, Egypt
  • 7Zoology Department, Faculty of Science, Ain Shams University, Cairo, Egypt
  • 8Anatomy Unit, Department of Basic Medical Sciences, College of Medicine and Medical Sciences, Qassim University, Buraydah, Saudi Arabia
  • 9Department of Anatomy and Cell Biology, Faculty of Medicine, Ain Shams University, Cairo, Egypt
  • 10Anatomy Unit, Department of Basic Medical Sciences, College of Medicine and Medical Sciences, AlNeelain University, Khartoum, Sudan
  • 11Pathology Unit, Department of Basic Medical Sciences, College of Medicine and Medical Sciences, Gassim University, Buraydah, Saudi Arabia
  • 12Clinical Pathology, Faculty of Medicine, Ain Shams University, Cairo, Egypt
  • 13Medical Physiology Department, Armed Forces College of Medicine, Cairo, Egypt
  • 14Department of Internal Medicine, Badr University in Cairo, Badr, Egypt
  • 15Pathology Department, Faculty of Medicine, Ain Shams University, Cairo, Egypt
  • 16Department of Anatomy and Cell Biology, Faculty of Medicine, Galala University, Attaka, Suez Governorate, Egypt
  • 17Department of Histology, Faculty of Medicine, Ain Shams University, Cairo, Egypt
  • 18Department of Clinical Pharmacology, Faculty of Medicine, Ain Shams University, Cairo, Egypt

Introduction: With the increasing prevalence of type 2 diabetes mellitus (T2DM), there is an urgent need to discover effective therapeutic targets for this complex condition. Coding and non-coding RNAs, with traditional biochemical parameters, have shown promise as viable targets for therapy. Machine learning (ML) techniques have emerged as powerful tools for predicting drug responses.

Method: In this study, we developed an ML-based model to identify the most influential features for drug response in the treatment of type 2 diabetes using three medicinal plant-based drugs (Rosavin, Caffeic acid, and Isorhamnetin), and a probiotics drug (Z-biotic), at different doses. A hundred rats were randomly assigned to ten groups, including a normal group, a streptozotocin-induced diabetic group, and eight treated groups. Serum samples were collected for biochemical analysis, while liver tissues (L) and adipose tissues (A) underwent histopathological examination and molecular biomarker extraction using quantitative PCR. Utilizing five machine learning algorithms, we integrated 32 molecular features and 12 biochemical features to select the most predictive targets for each model and the combined model.

Results and discussion: Our results indicated that high doses of the selected drugs effectively mitigated liver inflammation, reduced insulin resistance, and improved lipid profiles and renal function biomarkers. The machine learning model identified 13 molecular features, 10 biochemical features, and 20 combined features with an accuracy of 80% and AUC (0.894, 0.93, and 0.896), respectively. This study presents an ML model that accurately identifies effective therapeutic targets implicated in the molecular pathways associated with T2DM pathogenesis.

1 Introduction

Globally, the burden of diabetes mellitus (DM) is projected to rise to 1.3 billion people by 2050, making it one of the most widely spread diseases worldwide (1). Type 2 diabetes mellitus (T2DM) is the prevalent form of DM, it is hallmarked by hyperglycemia, insulin resistance, and ultimate decrease in β-cells insulin secretion (2). When it develops further worsening comorbidities emerge including micro- and macrovascular disease, leading to kidney dysfunction, diabetic retinopathy, blindness, heart disease, stroke, and lower limb amputations (3).

Diabetes mellitus is a highly heterogeneous entity (4). To enhance our comprehension of the underlying biological mechanisms and identify individuals at risk, it is crucial to investigate the genetic contributions to diabetes. Such knowledge can ultimately lead to the development of more precise and effective therapeutic approaches. As T2DM progresses, it often necessitates the simultaneous administration of multiple medications that target different pathophysiologic pathways (5). This combined treatment approach aims to regulate blood glucose levels and mitigate the progression of complications.

However, emerging evidence suggests that inflammatory pathways play a pivotal role as common mediators in the natural course of diabetes when influenced by risk factors (6). Interestingly, a previous study discovered that HFD-induced mice increased mitochondrial DNA (mtDNA) release into the cytosol of adipocytes, activating the cGAS-STING pathway and inflammatory response, resulting in chronic inflammation in adipose tissue and insulin resistance (7, 8). Moreover, autophagy plays a crucial role in T2DM pathogenesis as it protects cells from the damaging effects of oxidative stress and endoplasmic reticulum stress, which is essential for the survival and proper functioning of β-cell and insulin sensitivity, however, when the autophagic system in β-cells fails, it can exacerbate β-cell dysfunction, particularly in the presence of insulin resistance, potentially leading to hyperglycemia (9). A growing body of evidence suggests that enhanced autophagy, triggered by insulin resistance, may act as a safeguarding mechanism against the deterioration and increased apoptosis of pancreatic β-cells (10). This underscores the potential significance of autophagy modulation in the pursuit of therapeutic strategies aimed at preserving β-cell function in T2DM and ultimately managing hyperglycemia.

The growing availability of high-throughput technologies in large populations, such as genomics and transcriptomics with linked medical record data supports the development of new computational approaches for drug targeting using molecular biomarkers in addition to the traditional biomarkers. Noncoding RNAs (ncRNAs) including microRNAs (miRNAs) and long noncoding RNAs (LncRNAs), exhibit diverse functions in post-transcriptional gene regulation, epigenetic gene silencing, modulation of insulin secretion, and endoplasmic reticulum stress (11, 12). These ncRNAs have been intricately linked to the development of T2DM (13, 14) For instance, mir-375, a highly expressed miRNA in islet cells, plays a crucial role in insulin secretion and β-cell functioning. Dysregulation of mir-375 has been linked to impaired insulin secretion and β-cell dysfunction (15). On the other hand, lncRNAs have been implicated in the pathogenesis of insulin resistance and the maintenance of glucose homeostasis by regulating inflammatory and lipogenic processes (16, 17).

Despite the presence of chemical anti-diabetic agents, possible adverse effects and limited efficiency could occur. As a result, recent endeavors have explored other treatment options for the rising T2DM prevalence (18). Medicinal plants are of paramount importance in maintaining body health with no side effects compared with synthetic drugs (19). Likewise, several clinical studies have been conducted to examine the impact of probiotics on improving glycemic index, lipid profile, glucose metabolism, and insulin sensitivity (20). These investigations were fuelled by the observation that the gut microbiota of diabetic patients tends to be altered (21). Such changes in the gut microbiota could lead to metabolic endotoxemia, which occurs through the release of lipopolysaccharides and could trigger inflammation and insulin resistance (22).

Several researches has confirmed that inflammation is closely linked to the pathogenesis of T2DM and its complications. Many anti-diabetic drugs are usually prescribed to diabetic patients, to decrease the progression of T2DM through modulation of inflammation. However, those anti-diabetic drugs are often not successful as a result of side effects; such as sulphonylureas & biguanides may cause acute severe hypoglycemia and lactic acidosis. Therefore, researchers are searching for efficient natural therapeutic targets with less or no side effects. Natural products’ derived bioactive molecules have been proven to improve insulin resistance and associated complications through suppression of inflammatory signaling pathways (23). Moreover, the utilization of probiotics as dietary supplements gains popularity because gut microbiota dysbiosis significantly contributes to T2DM (24). Zbiotics as newly implemented enginerred probiotics as a new complementary therapeutic strategy that alleviate oxidative stress and beneficial effect on reducing blood glucose levels, HOMA-IR, and HbA1c.

This study investigates the potency of medicinal plant-based drugs and probiotics in modulating T2DM in streptozotocin-high-fat diet-induced induced rats. Isorhamnetin is a flavonoid found in sea buckthorn, medicinal plants, and ginkgo fruits (25). It possesses various pharmacological effects, including anti-inflammatory, anti-tumor, antioxidant, antibacterial, and antiviral properties (26, 27). It has been found to promote glucose uptake, maintain glucose homeostasis, and improve dyslipidemia in mice with T2DM (28). It can also reduce the expression of inflammatory cytokines and enhance the health of the gut microbiota in T2DM mice (29). Rhodiola rosea L. is abundant in flavonoids, glycosides, coumarins, and organic acid compounds (30). Rosavin, the R. rosea bioactive compound, possesses protective properties against inflammation, reduces blood glucose levels, exhibits antiviral and antitumor effects, and promotes blood circulation activation (31, 32). The caffeic acid extract derived from Artemisia dracunculus L. has been proven to enhance insulin receptor signaling. Semisynthetic compounds derived from caffeic acid induce DNA damage and apoptosis in tumor cells by activating autophagy (33). Furthermore, it plays a protective role in preventing renal damage (34, 35). ZBiotics is an engineered probiotic that involves the use of a genetically modified strain of B. subtilis by incorporating an acetaldehyde dehydrogenase gene. The modified strain converts acetaldehyde derived from ethanol into acetic acid,thereby reducing the potential harm caused by alcohol consumption (36). Animal toxicity studies have indicated ZBiotics’ high level of safety (37). However,its potential as a therapeutic drug for diabetes has not been investigated yet.

In order to examine the clinical trials’ data related to caffeic acid and related compounds in diabetic patients, we searched the largest clinical trial database at ‘https://clinicaltrials.gov’. No search results were obtained with the keywords ‘saffeic acid and diabetes mellitus’. Since propolis found in beehive is a major source of caffeic acid derivatives, therefore we searched with keywords ‘propolis and diabetes mellitus’ on the database. Three studies were found in which propolis was administered orally or applied topically to diabetic patients(ClinicalTrials.gov Identifier: NCT03416127 phase 2 in patients with type 2 DM for 12 weeks, NCT02794506 Phase: 4. in Type 2 DM, periodontitis, and NCT03649243 in diabetic foot ulcer). Similarly, no results found for rosavin in T2DM except one recent clinical trial in diabetic kidney disease(NCT06176599). More than 20 clinical trials were found correlating Rhodiola Rosea the origin of rosavin in many mental disorders, metabolic diseases, coronary diseases. Additionally, one clinical trial(NCT00961909) assessing the efficacy of Artemisia dracunculus plant rich in isorhamantin in T2DM. Moreover, no results found for Zbiotics in T2DM (Supplementary Table S1).

To minimize the consequences of diabetes and improve patient care, researchers have explored various fields, such as machine learning (ML) and artificial intelligence (AI), by applying ML techniques in the field of biology, researchers have significantly improved the precision of prediction models (38). Numerous studies have looked into using ML techniques to predict the occurrence of diabetes (39, 40).

Deberneh et al. created an ML model that can predict the occurrence of type 2 diabetes (T2D). The models categorize input data instances into three conditions: normal (non-diabetic), prediabetes, or diabetes. To construct their prediction model, they identified key features using a data-driven technique that includes an analysis of variance (ANOVA) test and recursive feature elimination methods. Also, they compared the performance of various machine learning models, such as LR, support vector machine (SVM), RF, and XGBoost algorithms (41). Wei S et al. built an ML model for diabetes detection. The study assessed two crucial data processing techniques: Principal Component Analysis (PCA) and Linear Discriminant Analysis (LDA) across various machine learning algorithms. The highest accuracy achieved among the five algorithms tested (Neural Network, Support Vector Machine, Decision Tree, LR, and Naïve Bayes) was 77.86% using 10-fold cross-validation (42).

Elsherbini AM et al. employed machine learning techniques to identify significant genes associated with diabetes and assess their potential as biomarkers for early detection. The analysis highlighted the HLA-DQB1 gene as a promising biomarker to detect diabetes through ML algorithms with decent accuracy (43). Xu et al. introduced a computational model utilizing stochastic gradient boosting. They incorporated six features, encompassing molecular structures, structural similarities, ATC code similarities, protein–protein interaction, chemical-chemical interaction, and disease pathways (44). Moreover, Costello et al. and Jang et al. performed extensive comparative analyses of machine learning methods for drug response prediction in cancer cell lines, recommending using elastic net or ridge regression with input features from all genomic profiling platforms (45, 46).

In this study, we aimed to utilize standard statistical methods and machine learning techniques to leverage the gene expression patterns and their epigenetic regulators in the livers and adipose tissues, along with conventional biochemical parameters to identify predictive features that could be used to assess the response to medicinal plant-based drugs and probiotics in an animal model of T2DM.

2 Materials and methods

2.1 Chemicals

Rosavin was obtained from Aktin Chemicals, Inc (Cat.No. APC-380, China), While isorhamnetin 3-O-acylglucosides, Z-Biotics, Caffeic acid, and sodium citrate buffer were purchased from Sigma Chemical Co, St. Louis, Mo, USA.

2.2 Experimental design

A hundred male Wistar rats (150–170 g) were purchased from the Holding Company for Biological Products and Vaccines based in Giza, Egypt. Rats were randomized into 10 groups (N=10 per group) and kept one week for acclimation with free access to normal rat chow and water under well-controlled condition (20 ± 2°C), 12 h light/dark cycle. The animal procedures followed the guidelines outlined in the National Institutes of Health guide for the care and use of Laboratory Animals (8th edition, 2011). Ethical approval for the experiments was obtained from the Institutional Animal Ethics Committee of Ain Shams University, Faculty of Medicine NO. FWA000017585. The Diabetic rat model was induced by feeding rats with a high-fat diet (HFD) that consists of 58% fat, 17% carbohydrate, and 25% protein and libitum for a total of 12 weeks. After the initial 4 weeks, the rats were administered two low-dose intraperitoneal (i.p.) injections of streptozotocin (STZ; 30 mg/kg) dissolved in citrate buffer (pH 4.5), with a one-week interval between injections. The Normal control rats were administered citrate buffer only. After one week following the last STZ injection, blood samples were collected from the tail vein, and blood glucose levels were measured using a glucometer (Accu-check, Roche Diagnostics, Risch-Rotkreuz, Switzerland). The onset of diabetes was confirmed when the non-fasting blood glucose levels were equal to or higher than 200 mg/dl. Throughout the study, the rats were allowed to continue consuming their respective diets until its completion. Then, T2DM-induced rats were randomly divided into 9 groups, 10 rats each, including; (I) the T2DM control group (N=10): that received intraperitoneal STZ, fed HFD, and 0.9% saline orally, (II) Rosavin-10, and 30 groups: the T2DM induced rats received 10 mg/kg-30 mg/kg of rosavin, respectively dissolved in 0.9% saline for 4 weeks, (III) Z-Biotic 0.5, and Z-Biotic 1 groups: T2DM induced rats received 0.5 mg/kg- 1 mg/kg Z-Biotic, respectively dissolved in DMSO for 3 weeks. (IV) Isorhamnetin-10, and Isorhamnetin-40 groups: the T2DM induced rats received 10 mg/kg-40 mg/kg Isorhamnetin dissolved in DMSO for 3 weeks, (V) Caffeic acid groups: the T2DM induced rats received 10 mg/kg-50 mg/kg caffeic acid, respectively dissolved in cold water. Rats received the medicinal plant-based drugs and the probiotics orally by gastric gavage. (VI) The Normal group (N=10): rats received sodium citrate buffer 1 ml/kg intraperitoneally, the same amount injected in the weight-matched rat in the other groups. At the end of the experiment, blood samples were obtained from the retro-orbital veins of the animals under ether anesthesia. Subsequently, the animals were euthanized through cervical dislocation. Both the right and left gastrocnemius muscles, along with adipose tissue, were collected from all the animal groups. One gastrocnemius muscle and a portion of adipose tissue from each animal were frozen at -80°C for biochemical analysis. The other gastrocnemius muscle and a portion of adipose tissue were fixed in 10% formalin for histopathological evaluation (Figure 1).

Figure 1
www.frontiersin.org

Figure 1 Workflow describing the animal groups subjected to molecular and biochemical analysis for ML-model building identify signatures associated with drug response.T2DM, Type 2 Diabetes Mellitus; (L) extracted from the Liver tissues; (A) extracted from the Adipose tissues; HOMA-IR, Homeostatic Model Assessment of Insulin Resistance; LDL-C, low-density lipoprotein cholesterol; HDL-C, high-density lipoprotein cholesterol; AST, aspartate transaminase; ALT, alanine transaminase; BUN, Blood Urea Nitrogen; ACR, urine albumin to creatinine ratio; MLA, Machine learning Algorithm.

2.3 Biochemical analysis

Prior to euthanization, blood samples were obtained from the optical vein of the rats and then subjected to centrifugation at 2000g for 10 minutes at 4°C. The resulting serum was collected for further analysis. Commercial ELISA kits from RayBiotech, USA, were utilized to measure the levels of cholesterol, triglycerides, high-density lipoprotein (HDL), low-density lipoprotein (LDL), alanine transaminase (ALT), aspartate transaminase (AST), total cholesterol (TC), serum creatinine, Blood Urea Nitrogen (BUN), fasting blood glucose, postprandial blood glucose, and fasting blood insulin following the manufacturer’s instructions. The calculation of the Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) was performed using the equation: fasting insulin (μU/L) multiplied by fasting glucose (nmol/L)/22.5. Urine samples were collected for one day using individual metabolic cages, which took place on the day before the completion of the treatment. The level of albumin in the urine was evaluated using commercially available colorimetric kits from RayBiotech, USA, following the instructions provided by the supplier.

2.4 Histopathological analysis

Liver and adipose tissue specimens underwent dissection and fixation in neutral buffered formalin for 72 hours. Subsequently, the samples underwent a series of processing steps involving ethyl alcohol, Xylene clearance, and embedding in paraplast tissue-embedding media. Employing a rotatory microtome, tissue sections with a thickness of 5μm were acquired and affixed onto glass slides. The staining procedure utilized Hematoxylin and Eosin (H and E), following standard protocols outlined by Culling, C.F.A. 2013. Skilled histologists, operating in a blinded fashion, scrutinized these tissue sections.

2.5 Morphometric analysis

For the analysis of white fat cells and the determination of average cell diameter, a minimum of 6 non-overlapping fields were randomly selected and scanned. Following the approach outlined by Batts and Ludwig (47), individual biopsy specimens were assessed for the grade of inflammation, rated on a scale from 0 to 4 [0: no activity; 1: minimal; 2: mild; 3: moderate; and 4: severe]. All data were acquired through the Leica Application module linked to Leica Microsystems GmbH (Germany).

2.6 Bioinformatics analysis

● In this study, we conducted a comprehensive analysis of genes that are differentially expressed relevant to diabetes. To identify these genes, we performed a search in the Gene Expression Omnibus (GEO) database http://www.ncbi.nlm.nih.gov/geo, accessed on Jan 2024. In this study, we conducted a comprehensive analysis of genes that are differentially expressed relevant to diabetes. To identify these genes, we performed a search in the Gene Expression Omnibus (GEO) database http://www.ncbi.nlm.nih.gov/geo, accessed on Jan 2024. We utilized keywords related to diabetes mellitus type 2 such as ‘type 2 diabetes,’ ‘diabetic nephropathy,’ ‘tissue,’ ‘pancreas’. After obtaining a pool of potential datasets, we proceeded to screen them based on specific inclusion criteria. The first criterion was related to the tissue type, where we focused on tissues relevant to T2DM. The second inclusion criterion was the availability of normal tissues used as controls in the dataset. The inclusion of normal tissues ensures a suitable baseline for comparative analysis with tissues from T2DM patients. Furthermore, we also considered the sample size in the dataset. To ensure statistical robustness and reduce potential bias, we set a minimum threshold of ten or more samples in each dataset.

Specifically, we selected two datasets, GSE20966 and GSE142025, and used the “GEO2R/Limma R” package to screen for highly significant differentially expressed genes (DEGs). We considered genes with a p-value less than 0.05 and a log twofold change (LogFC) value of ≥ 1 or ≤ -1, using the Benjamini and Hochberg method for false discovery rate adjustment. In order to identify the biological pathways of DE-mRNAs, a thorough enrichment analysis was conducted using Enrichr (http://amp.pharm.mssm.edu/Enrichr, Jan 2022) with the Kyoto Encyclopedia of Genes and Genomes (KEGG) selected as the analysis tool. The analysis revealed The DE-mRNAs were related to Pancreatic secretion, ECM-receptor interaction, HNF3B pathway,RANKL regulation of apoptosis and immune response(Supplementary Figure S2). then filtered based on their association with, specific pathways that are of interest in diabetes research, such as insulin resistance, autophagy, cGAS/STING, and NOD-like receptor pathways(Supplementary Figure S2). To identify genes associated with these pathways, we utilized the GeneCards database https://www.genecards.org/, accessed on Jan 2024. Accordingly, we selected nine genes: Z-DNA Binding Protein 1 (ZBP1), Stimulator Of Interferon Response CGAMP Interactor 1 (STING1), RIG-I, Retinoic Acid-Inducible Gene 1 Protein (DDX58), mammalian target of rapamycin (mTOR), Nuclear Factor Kappa B Subunit 1 (NFKB1), insulin-like growth factor-1 (IGF-1), Component Of Inhibitor Of Nuclear Factor Kappa B Kinase Complex (CHUK), RAC-beta serine/threonine-protein kinase (AKT2), and Ret Proto-Oncogene (RET). To explore the interactions between the differentially expressed genes, we utilized the String database https://string-db.org/, accessed on Jan 2024 to construct a protein-protein interaction (PPI) network. Furthermore, to identify miRNAs that target the selected DEGs, we used the mirwalk database http://mirwalk.umm.uni-heidelberg.de/, accessed on Jan 2024. Finally, to predict the interaction of miRNAs-LncRNAs, we used RNA22 https://cm.jefferson.edu/rna22/, accessed on Jan 2024, mirwalk database http://mirwalk.umm.uni-heidelberg.de/, accessed on Jan 2024 and DIANA Tools https://diana.e-ce.uth.gr/lncbasev3/interactions, accessed on Jan 2024. Thus, the following LncRNAs (AC074117.2, RP11–773H22.4, and RP4–605O3.4) and miRNAs (miR-1976, miR-1, miR-611, and miR-3163) were chosen (Supplementary Table S2, S3, Supplementary Figures S3-S8).

2.7 Total RNA extraction and quantitative real-time PCR

Total RNA was isolated from the Liver (L) and Adipose tissue (A) samples using the miRNeasy kit (Qiagen, USA; Cat no. 74104). Next, the RNA quality, integrity, and concentration were measured using the DeNovix DS-11 microvolume spectrophotometer (Wilmington, USA) and stored at -80°C. The obtained RNA samples were reverse transcribed into cDNA by the two-step RT-PCR using the miScript II RT kit (Qiagen, USA; Cat no. 218161). Then the qRT-PCR was performed using Applied Biosystems Tm 7500 system (Foster City, California, United States). Regarding the expression of mRNAs, QuantiTect SYBR Green PCR Kit (Qiagen, Helman Germany; Cat no. 204143) was used. The relative expression of the miRNAs was obtained using the miScript SYBR Green PCR Kit (Qiagen, Helman Germany; Cat no. 218073). The relative expression of the lncRNAs was performed using RT2 SYBR Green ROX qPCR Master mix (Qiagen, Helman Germany; Cat no: 330500). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and U6 were regarded as internal controls for the mRNAs, LncRNAs, and miRNAs, respectively. The Livak method RQ = 2-ΔΔCt method was adopted to analyze the relative expression. For each sample, 2 replicates were set. Ct values of more than 35 were deemed negative. Melting curve analysis validated the amplicons’ specificities for SYBR Green-based PCR amplification. In this study, proper standardization procedures were used to detect any experimental error produced at any stage of the RNA extraction and processing according to MIQE recommendations. The PCR procedure was as follows: an initial activation phase at 95°C for 15 minutes was followed by 40 cycles of PCR at 94°C for 15 s, 55°C for 30 seconds, and 72°C for 30 s.

2.8 Statistical analysis

The Statistical analysis was performed with SPSS 26. Software. Data are presented as mean ± SD and significant differences were compared using a 1-way analysis of variance followed by post hoc Tukey’s test. The Shapiro-Wilk test confirmed that the data in this study followed a normal distribution. Categorical data were expressed as percentages and compared using the chi-square test.

3 Machine learning model

3.1 Data

One of the main objectives of this study is to create a predictive model using machine learning algorithms to thoroughly identify a promising set of features that strongly indicate the drug’s response to T2DM. The dataset in this study is a mouse model data with 100 samples distributed among normal, treated, and diseased groups (Table 1) with molecular features expressed in adipose tissues (A) and liver tissues (L) besides the biochemical features (Table 2). To determine treatment response, two key features were focused on: inflammation grade and fat cell diameter. Samples exhibiting a fat cell diameter exceeding 70 and an inflammation grade of 3 or higher were categorized as “not improved” (binary 0), while those falling below these thresholds were labeled as “improved” (binary 1). The analysis showed the relative ratio is 1:2.44 for the improved and not improved samples as shown in Figure 2.

Table 1
www.frontiersin.org

Table 1 The number of samples per normal, disease model, and treatment groups.

Table 2
www.frontiersin.org

Table 2 Molecular and biochemical features used in ML models.

Figure 2
www.frontiersin.org

Figure 2 Show the number of improved (48) and not improved (29) samples.

The dataset was divided into 70 samples for training and 30 samples for testing. Also, Multi-classifiers approaches were implemented using 5 classifiers including K-Nearest Neighbors (KNN), Light Gradient Boosting Machine (LGBM), Random Forest (RF), Logistic Regression (LR), and Ada Boost Classifier.

3.2 Dataset preprocessing

In the development of machine learning models, data preprocessing stands as a critical step. Herein, thorough measures were taken to ensure the data’s quality by addressing several factors and they were null feature removal, noise reduction, and outlier elimination. Additionally, the column index was reset to accurately assign each column its respective index. Moreover, the creation of a ‘status’ column in binary format (1,0) was undertaken to reflect treatment response whether it was improved or not improved (49).

3.3 Feature selection

The feature selection technique was utilized, this method effectively reduces the complexity and size of the data, enhancing the learning process. Additionally, by selecting only the relevant features, the model becomes faster and more precise, thereby enhancing its predictive capabilities through noise reduction.

3.3.1 Selecting top features with recursive feature elimination cross validation

An initial trial was carried out to determine the selected machine learning features for predicting T2DM drug response. The Recursive Feature Elimination Cross-Validation (RFECV) technique was employed for this purpose. RFECV functions by iteratively eliminating features and evaluating prediction accuracy. The most optimal reduced set of features, which either matches or surpasses the original accuracy of the complete feature set, is then identified as the optimal set of predictive features. Various combinations of features yield different accuracy levels (50). Furthermore, relevant model hyperparameters for the RF model (e.g., maximum tree depth, minimum samples split, and estimators) were fine-tuned to achieve the best-performing model. Optimization involves experimenting with different parameter sets to determine the most effective configuration for our dataset. The top selected features for molecular, biochemical, and both combined were reported.

3.4 Cross validation

In k-fold cross-validation, the training set is partitioned into k distinct and equally-sized subsets. The classifier is trained on each subset using the union of all other subsets. Consequently, the average error rate across all subsets provides an estimate of the classifier’s error rate. Each data point is included in the test set precisely once and appears in the training set k-1 times. Increasing k helps reduce the variance of the resulting estimate. We utilized the k-fold cross-validation method, employing k=3 in our analysis (51).Of note, we utilized Stratified K-Fold cross-validation with a number of folds set to 3 to validate the performance of our machine learning model. This technique was chosen to ensure that each fold preserves the proportion of samples for each class, thereby mitigating potential biases and enhancing the reliability of our results. During the validation process, we partitioned the dataset into three subsets, with each subset containing a representative distribution of the target variable classes. We then iteratively trained the model on two-thirds of the data (training set) and evaluated its performance on the remaining one-third (validation set), ensuring that the model was tested on unseen data in each fold. By adopting Stratified K-Fold cross-validation, we aimed to demonstrate the independence and rigor of our model validation, thereby enhancing the reproducibility and accuracy of our results.

3.5 Models predictions

A multi-class approach using five machine learning algorithms (KNN, RF, LGBM, LR, and Ada Boost) was employed on three separate models as follows: one using only Molecular features, one using only Biochemical features, and one combining both feature groups as shown in Table 3. The top selected features obtained from RFECV were used in each model. This comprehensive strategy aimed to identify a potential set of highly predictive features for drug response. Also, the top-performing classifier applied on the training set for each feature group was selected and then applied to the testing set to evaluate their performance on an unseen dataset.

Table 3
www.frontiersin.org

Table 3 The three predictive models were applied to the five classifiers.

3.6 Machine learning evaluation

The three model’s performance were assessed using the test dataset. The evaluation was based on key performance metrics, such as the area under the curve (AUC). To further analyze the model’s efficacy, a 2 × 2 confusion matrix was constructed using the test dataset, allowing for the calculation of true positive, false positive, true negative, and false negative values.

3.7 Visual analysis for molecular and biochemical features

a graphical implementation has been used to analyze the relationship between each molecular and biochemical features of treated samples against the T2DM group as well, as to the normal group to capture insights for drug response.

3.8 Python packages

This study’s data were processed using Python 3.7 as the programming language. We used many Python-based packages and modules as well to ease the processing pipeline. The ‘pandas’ package (version 1.3.5) and ‘NumPy’ (version 1.20.3) are utilized for data manipulation and analysis. ‘Seaborn’ (version 0.13. 2) enhances data visualization capabilities. ‘Matplotlib.pyplot’ (version 3.5.0) is another data visualization package, offering a versatile toolkit for creating static, interactive, and animated plots. The scikit-learn (version 1.0.2) is extensively used for machine learning tasks. ‘statistics’ provides functions for mathematical statistics.

4 Results

4.1 Effect on diabetic parameters, serum biochemicals, and liver inflammation

Table 4 presents the observed changes in different parameters in the T2DM group compared to the normal group, along with the effects of medicinal plant-based drugs and probiotics drug administration. In the T2DM group, there was a significant increase in serum glucose, insulin, and HOMAIR levels. However, these elevated levels were modulated after the administration of our selected drugs, particularly at higher doses, suggesting that the drugs have the potential to suppress hyperinsulinemia and improve insulin resistance. Furthermore, renal function biomarkers including serum creatinine, BUN, and albumin-to-creatinine ratio (ACR) showed a significant increase in the T2DM group. On the other hand, the administration of different doses of drugs led to a modulation of these biomarkers. The lipid profile analysis revealed a significant increase in total cholesterol, triglycerides, and LDL levels, along with reduced HDL levels in the T2DM rat group. While the treated groups exhibited improvements in the lipid profile levels. Moreover, hepatic damage biomarkers such as AST and ALT were elevated in the T2DM group but decreased in the treated groups. Similarly, liver inflammation grades exhibited an increase in the T2DM group but decreased in the treated groups, particularly at higher doses. The fat cells’ diameter increased significantly in the T2DM group as well as the treated groups except for the Rosavin-30 treated group which showed decreasing in the fat cells’ diameter.

Table 4
www.frontiersin.org

Table 4 Serum biochemical parameters, liver inflammation in the normal, T2DM, and treated rats.

4.2 Histopathological results

The normal group demonstrated a normal architecture. The hepatic lobule is composed of hepatocytes (h) that are arranged into branching cords with sinusoids (S) in between. Hepatocytes appeared as polygonal cells with granular eosinophilic cytoplasm and round nuclei. Central vein (CV) is located at the center of the lobule whereas the T2DM group showed an altered architecture with marked hepatocellular microvesicular steatosis (Arrow) all over hepatic lobule accompanied, widespread of necrotic cells (Arrowhead), inflammatory cells infiltrate (asterisks) and pyknotic cells (Curved arrow). The Rosavin-10 treated group showed altered architecture with mild hepatocellular microvesicular steatosis (Arrow), in addition to the presence of pyknotic cells (Curved arrow), necrotic cells (Arrowhead), and inflammatory cells (asterisks) whereas the Rosavin-30 treated group showed a normal histological feature of the hepatic lobule with wide sinusoids (S). The Z-Biotic 0.5 treated group showed an altered architecture of the liver with moderate hepatocellular microvesicular steatosis (Arrow), and the presence of inflammatory cell infiltrates (asterisks) whereas the Z-Biotic 1 treated group showed an altered architecture of the liver with mild hepatocellular microvesicular steatosis (Arrow) with minimum mononuclear inflammatory cells infiltrates. The Isorhamnetin-10 treated group demonstrated necrotic changes (Arrowhead) in the hepatocytes, several hepatocytes showed pale cytoplasm with vacuolated nuclei. Inflammatory cells (asterisks) are scattered in the hepatic lobule whereas the Isorhamnetin-40 treated group demonstrated a slightly normal histological feature. The Caffeic acid-10 treated group showed an altered architecture with hepatocellular microvesicular steatosis (Arrow), scattered necrotic cells (Arrowhead), and inflammatory cell infiltrates (asterisks) whereas the Caffeic acid-50 treated group demonstrated a slightly normal histological feature (Figure 3). Analysis of adipose tissue revealed a significant increase in fat cell diameter in T2DM animals. However, treatment of the diabetic animals with rosavin, Z-biotic, isorhamnetin, or Caffeic acid exhibited notable decreases in epididymal adipocyte cell size compared to the normal group. Interestingly, the most pronounced ameliorative effect was observed in the high-dose-treated groups (Figure 4).

Figure 3
www.frontiersin.org

Figure 3 Histopathology examination of liver tissue in experimental animals. (A) the Normal group, (B) the T2DM group, (C) Rosavin-10 group, (D) Rosavin-30, (E) Z-Biotic 0.5, (F) Z-Biotic 1, (G) Isorhamnetin-10, (H) Isorhamnetin-40, (I) Caffeic acid-10, (J) Caffeic acid-50. Microvesicular steatosis (Arrow), Inflammatory cells infiltrate (asterisks), Necrotic cells (Arrowhead), Pyknotic hepatocytes (Curved arrow), Sinusoids (S), Central vein (CV), Portal vein (PV), Bile duct (BD). Hepatocytes (h) (10X magnifications).

Figure 4
www.frontiersin.org

Figure 4 Histopathology examination of adipose tissue in the experimental animals. (A) the Normal group, (B) the T2DM group, (C) the Rosavin-10 group, (D) Rosavin-30, (E) Z-Biotic 0.5, (F) Z-Biotic 1, (G) Isorhamnetin-10, (H) Isorhamnetin-40, (I) Caffeic acid-10, (J) Caffeic acid-50. (40X magnifications).

4.3 Effect on molecular targets mRNAs-miRNAs-LncRNAs

The results presented in Table 5 demonstrate significant changes in various biomarkers among the different experimental groups expressed in liver tissues (L) and adipose tissues (A). In comparison to the normal group, the T2DM group exhibited markedly elevated levels of ZBP1, STING1, DDX58, mTOR, NFKB1, IGF1R, CHUK, RET, hsa-miR-1976, hsa-miR-611, and lnc-RP11–773H22.4 (p < 0.05) while AKT2, hsa-miR-1, hsa-miR-3163, lnc-AC074117.2, and lnc-RP4–605O3.4 were decreased. However, after drug administration, particularly at higher doses, this disturbance was significantly modulated.

Table 5
www.frontiersin.org

Table 5 The differential expression of mRNAs-miRNAs-LncRNAs in liver tissues (L) and adipose tissues (A) among different animal groups, including a normal group, T2DM group, and various treatment groups with Rosavin, Z-Biotic, Isorhamnetin, and caffeic acid.

4.5 ML Analysis

4.5.1 Feature selection using RFECV-based Random Forest for T2DM drug response prediction

The feature selection procedure employing RFECV yielded 13 features out of 32 for the molecular model, 10 out of 12 for the biochemical model, and 20 out of 44 for the combined model. These selections, as shown in Table 6, were made while maintaining the same prediction accuracy level as shown in Figure 5. Molecular accuracy achieved 80% while the biochemical and both of them combined achieved 85%.

Table 6
www.frontiersin.org

Table 6 Show the top selected features by RFECV for each model.

Figure 5
www.frontiersin.org

Figure 5 Number and accuracy scores for each feature set and both of them combined.

4.5.2 Model prediction results

Initial results applied on the training set showed the top-performing classifiers for each feature group. Table 7 summarizes the accuracy of the adopted classifiers for each feature group. Notably, the KNN was the top classifier for the molecular model and the combined model while LGBM Classifier performed the best for the biochemical model.

Table 7
www.frontiersin.org

Table 7 The evaluation metric for the top-performing classifier on the training set.

Then, the selected classifiers were applied to the testing set for the T2DM drug response prediction and to evaluate their predictive performance on unseen data. Using this strategy, we ensured that only the most effective classifiers were applied for prediction for each feature group, thereby enhancing the reliability and robustness of our predictive models. Table 8 summarizes the evaluation metric for the testing set. Notably, all classifiers in the three models achieved the same prediction power (0.8).

Table 8
www.frontiersin.org

Table 8 The evaluation metric for the best classifiers on the testing set for each feature group.

4.5.3 Performance evaluation of machine learning models for drug response prediction

The confusion matrix presented in Figure 6 displays the correctness of the prediction of whether the sample is improved (1), or not improved (0) on the test set, for the molecular, biochemical, and combined model. Roc curve displays how accurate the prediction models are (Figure 7).

Figure 6
www.frontiersin.org

Figure 6 Confusion matrix for the best-performing classifiers for each feature group. (A) Molecular, (B) Biochemical, (C) Molecular+Biochemical.

Figure 7
www.frontiersin.org

Figure 7 Roc curve for the best-performing classifiers for each feature group for the prediction of drug response. (A) Molecular, (B) Biochemical, (C) Molecular+Biochemical.

4.5.4 Visual analysis results for molecular and biochemical features

The implemented graphical approach displayed the differences in feature values by showing measures such as mean values and standard deviations between the different treated groups (drugs), against the normal and the T2DM groups (Supplementary Figures S6, S7).

5 Discussion

Diabetes mellitus is a metabolic disorder that impacts various organs and is primarily characterized by a deficiency in insulin production or response. This deficiency leads to impaired glucose tolerance and high blood glucose levels, accompanied by disruptions in the metabolism of fats, carbohydrates, and proteins (52).

Various factors, including modified dietary patterns, metabolic stress, and genetic predisposition, have the potential to activate the innate immune system. This activation leads to insulin resistance, T2DM, and related complications like dyslipidemia, diabetic nephropathy, retinopathy, and atherosclerosis. Among the drugs tested in our study, rosavin has garnered attention due to its antioxidant and anti-inflammatory properties, which help mitigate oxidative stress and chronic inflammation often observed in diabetes (53). Rosavin has been shown to lower total cholesterol, TGs, and LDL levels while increasing HDL, thereby exerting lipid-lowering effects. Insulin resistance precedes the abnormal elevation of blood glucose levels, which is the primary clinical indicator of T2DM. In the prediabetic stage, the body responds to insulin resistance by increasing insulin production to meet the normal insulin requirements. This leads to a state of chronic hyperinsulinemia and, as a result of sustained high blood sugar levels, the failure of pancreatic beta cells. Eventually, this progression leads to the development of T2DM. Prolonged and exaggerated metabolic stress can lead to detrimental inflammatory reactions, resulting in insulin resistance and inflammatory diseases. This chronic inflammatory state eventually gives rise to long-term complications associated with diabetes, including microvascular complications such as diabetic liver disease, diabetic nephropathy, neuropathy, and macrovascular complications like cardiovascular and cerebrovascular diseases. The cGAS-STING pathway plays a crucial role in diabetic complications and has been increasingly reported in diabetic nephropathy and diabetic angiopathy, which is linked to mitochondrial dysfunction caused by lipotoxicity. Evidence suggests that the cGAS-STING pathway is over-activated in diabetes and its complications. This heightened activation of cGAS-STING may serve as a protective mechanism, considering that diabetic patients are more susceptible to infections. Notably, knocking out STING has been shown to reduce insulin resistance induced by a high-fat diet in peripheral tissues and improve overall glucose intolerance. However, it is important to note that STING deficiency also impairs the ability of beta cells to secrete insulin in response to glucose stimulation. Autophagy is crucial in T2DM, serving to sustain cellular energy levels during fasting while eliminating damaged components like organelles, lipids, and misfolded proteins. Moreover, it contributes to pancreatic beta cell function and insulin resistance. Recent findings highlight the role of autophagy in T2DM pathophysiology, particularly in maintaining pancreatic beta cell function. Furthermore, increased autophagy serves as a protective mechanism against oxidative stress in insulin-targeted tissues like the liver, adipose tissue, and skeletal muscle. NOD1 and NOD2 are implicated in the development of diabetes, likely through their interaction with the gut microbiota. Antibiotic-induced changes in the gut microbiota are crucial for enhancing insulin sensitivity. Nucleotide-binding oligomerization domain (NOD)-like receptors, such as NOD1 and NOD2, both recruit receptor-interacting protein kinase 2 (RIPK2), but they exert opposite effects on blood glucose regulation. While NOD1 links signals from bacterial cell walls to metabolic inflammation and insulin resistance, NOD2 can foster immune tolerance, improve insulin sensitivity, and enhance blood glucose control in obesity. Similarly, NLR family pyrin domain-containing (NLRP) inflammasomes can elicit different metabolic outcomes. NLRP1 may protect against obesity and metabolic inflammation, possibly due to its preference for regulating IL-18, whereas NLRP3 tends to promote IL-1β-mediated metabolic inflammation and insulin resistance. Additionally, Rosavin may protect against diabetic complications such as neuropathy and nephropathy by reducing nerve damage, kidney injury, and hepatic damage as indicated by decreased levels of creatinine, BUN, ALT, and AST (30). Specifically, rosavin has demonstrated protective anti-inflammatory effects in various models, including bleomycin-induced pulmonary fibrosis. Rosavin exerts its effects by downregulating the expression of pro-inflammatory molecules such as NF-κB, p65, TGF-β1, and α-SMA while upregulating the expression of nuclear erythroid 2-related factor 2 (Nrf2), a transcription factor involved in antioxidant defense (54) Our findings supported the hypoglycemic (lower serum glucose, insulin, and HOMA-IR) and hypolipidemic (lower serum TG, TC, LDL, and higher HDL levels) effects of rosavin treatment (55). Mao. reported that Rhodiola rosea L. root extracts improved oral glucose tolerance, decreased serum TG, and LDL levels, and increased HDL levels in KKAy mice, a T2D model (56). Liu et al. proved that rosavin has been shown to attenuate LPS-induced activation of the TLR-4/NF-κB signaling pathway in RAW264.7 cells and inhibit the release of inflammatory factors in A549 cells. In a dose-dependent manner, rosavin ameliorated histopathological alterations, reduced the levels of inflammatory factors, and inhibited the TLR-4/NF-κB/MAPK signaling pathway and apoptosis activation. It also significantly reduced the number of inflammatory cells in bronchoalveolar lavage fluid and the expression of NF-κB p65 protein in the lung tissue of a mouse model. Moreover, it reduced the expression of hydroxyproline and malondialdehyde while enhancing the activities of superoxide dismutase and glutathione peroxidase in lung tissue (57).

Zbiotics, a newly engineered probiotic, has not been fully investigated for its impact on the pathogenesis of DM. Therefore, we conducted a study to assess its effectiveness on biochemical and molecular markers associated with DM. Both low and high doses of zbiotics demonstrated hypoglycemic and hypolipidemic effects, along with hepatoprotective and renoprotective effects evidenced by reductions in AST, ALT, creatinine, and BUN levels, particularly noticeable with the high dose. Histopathological analysis revealed mild hepatocellular microvesicular steatosis with minimal inflammatory cell infiltration. Our findings indicate that zbiotics with engineered acetaldehyde dehydrogenase can eloquently explain the results of the present study regarding how ZBiotics® reversed the toll of diabetes on the studied histopathological, biochemical and molecular levels. The concept was based on two major mechanisms that could elicit inflammation that leads over time to diabetic settings and evolves diabetic complications; firstly, the oxidative-stress-induced inflammation, and secondly the disturbance in diabetic gut microbiota and barriers that could lead to the activation of subsequently inflammation.

Isorhamnetin has emerged as a promising therapeutic agent for T2DM by improving gut health and insulin resistance (29). Previous studies have demonstrated its efficacy in lowering fasting blood glucose levels, improving renal function, and ameliorating dyslipidemia in T2DM rats by upregulating autophagy in renal tissues (26). In animal models, Isorhamnetin supplementation has been found to reduce reactive oxygen species levels, inhibit atherosclerotic plaque formation, and mitigate myocardial hypertrophy and fibrosis ((58) (Mechanistically, Isorhamnetin acts through various pathways, including inhibition of the PI3K/AKT signaling pathway, activation of the AMPK/mTOR pathways, and modulation of insulin secretion via phosphorylation of insulin receptor substrate-2 (IRS-2), phosphatidylinositol 3-kinase (PI3K), Akt, and activated pancreatic and duodenal homeobox-1 (PDX-1) (59, 60). Moreover, Yang et al. proved that isorhamnetin exhibited hepatoprotective effects by reducing liver fibrosis through inhibition of HSC activation and ECM formation, and by downregulating the TGF-β1/Smad3 and TGF-β1/p38 MAPK pathways (61). In our previous study, we demonstrated that Isorhamnetin acted as an effective therapy for DM by modulating the insulin resistance signaling pathway and autophagy-related RNA network (28). Building upon our previous findings, our current study expands on the beneficial effects of isorhamnetin. We observed that isorhamnetin improved insulin resistance parameters, reduced elevated glucose levels, and alleviated inflammation infiltration in liver sections. Furthermore, isorhamnetin treatment led to the recovery of elevated levels of ALT and AST, indicating its hepatoprotective effects against chronic injury. isorhamnetin efficiently improves altered lipid metabolism by decreasing TGs, LDL, and TC, while increasing HDL. These lipid-lowering effects contribute to the hepatoprotective role of isorhamnetin, which can be attributed to its anti-inflammatory properties. Additionally, Lu et al. demonstrate that Isorhamnetin affects the P38/PPAR-α pathway, which in turn regulates the expression of apoptosis and autophagy-related proteins (62). This finding is consistent with our study, where Isorhamnetin modulated autophagy-related genes such as STING1, IGF-1, and AKT2.

Caffeic acid demonstrates anti-diabetic effects through multiple mechanisms, including the enhancement of antioxidant enzymes, inhibition of NF-κB signaling pathways, and activation of the NrF2 transcription factor (63). In a study by Xu et al., it was shown that caffeic acid administration at doses of 25 and 35 mg/kg significantly reduced plasma glucose, TG, TC, and LDL levels, while notably increasing HDL, insulin, and antioxidant levels in streptozotocin-induced diabetic Wistar rats after five weeks of treatment (64). The findings from Xu’s study align with our research, where we investigated the effects of caffeic acid treatment at doses of 10 and 50 mg/kg and we observed significant improvements in plasma glucose, lipid profiles, and insulin levels in diabetic rats treated with caffeic acid. Additionally, our results showed hepatoprotective and renoprotective effects of caffeic acid, along with a reduction in inflammatory cell infiltration in the liver and a decrease in epididymal adipocyte cell size, particularly with the higher dose. Furthermore, our previous research indicated that caffeic acid may activate the mitogen-activated protein kinases (MAPK) signaling pathway through the regulation of miR-636, leading to the induction of autophagy and attenuation of diabetic nephropathy (65). Consistent with this, our recent findings demonstrate that caffeic acid treatment significantly downregulates the expression of autophagy-related genes such as IGF1R, NFKB1, and STING1 in adipose and liver tissues, further supporting its beneficial effects in diabetes management. Bhattacharya et al. investigated the effects of caffeic acid on glucose sensitivity, glucose-stimulated insulin secretion (GSIS), and gene expression in INS-1E cells under normoglycemic conditions (NC) and glucotoxic conditions (GC). They found that caffeic acid significantly increased the expression of Insulin-1 (Ins-1), Ins-2, pancreatic and duodenal homeobox 1 (Pdx-1), Akt-1, Akt-2, insulin receptor substrate-1 (Irs1), Bcl2, heat shock protein 90 and 70 (Hsp90 and Hsp70) during NC. Additionally, caffeic acid downregulated acetyl coenzyme A carboxylase 1 (ACC1) without affecting Glucokinase (Gck) and Glucose transporter-2 (Glut-2) expressions in INS-1E cells. However, under GC conditions, caffeic acid did not change the expression of GLUT-2, Gck, Ins2, Beta2, Pdx1, Akt2, Irs1, Bcl2, and Hsp90. Instead, it upregulated Ins1, Akt2, and Hsp70, while downregulating Beta2, Caspase 3 (Casp3), and Bax. caffeic acid also significantly increased glucose sensitivity and GSIS in INS-1E cells and thereby caffeic acid may enhance the survival and function of β-cells during glucotoxic conditions by modulating the expression of these genes (66).

Many human and animal studies investigated prolonged drug exposure on both safety and efficacy outcomes. Lekomtseva et al. investigated the effects of Rhodiola rosea extract on prolonged or chronic fatigue symptoms, 100 subjects were administered 2 × 200 mg of the extract daily over 8 weeks of an open-label clinical trial. Results showed the greatest improvement after just 1 week of treatment, with continued reduction in fatigue symptoms throughout the study, reaching statistically significant improvement by week 8. Importantly, safety assessments indicated favorable outcomes, with most adverse events being mild and unrelated to the study drug (67). Ochoa-Morales et al. conducted a 12-week double-blind, randomized placebo-controlled trial that assessed the efficacy of propolis compared to placebo in controlling glycemic levels in 36 patients with T2DM. Administered twice daily before breakfast and dinner, propolis (300 mg) significantly reduced fasting plasma glucose (FPG) and 2-hour post-load glucose (2-h PG) levels compared to placebo (48). In a 6-month masked, randomized clinical trial conducted by El‐Sharkawy et al. individuals with chronic periodontitis (CP) and T2DM undergoing scaling and root planning (SRP) were given either a placebo or a daily regimen of 400 mg oral propolis. Results showed that the propolis group exhibited significant reductions in hemoglobin A1c (HbA1c) levels by 0.82% and 0.96% units at 3 and 6 months, respectively, along with decreases in fasting plasma glucose (FPG) and serum N€-(carboxymethyl) lysine (CML) levels. Additionally, both groups showed improvements in periodontal parameters after therapy, but the propolis group demonstrated significantly greater reductions in probing depth and gains in clinical attachment level compared to the control group (68). In a study investigating the effects of caffeic acid on diabetic cardiomyopathy, it was found that caffeic acid, along with ellagic acid, demonstrated protective effects in diabetic mice. Various parameters, including lipid profile, coagulability, oxidative stress, and inflammation, were assessed. After 12 weeks, the treated animals showed beneficial effects, including decreased triglyceride levels, increased plasma insulin levels, decreased plasma glucose levels, anti-coagulatory effects, antioxidative effects, and anti-inflammatory properties in the cardiac tissue (69). Moreover, Rodríguez-Rodríguez et al. investigated the metabolic effects of an extract from Opuntia ficus-indica (OFI) for 12 weeks, known for its high isorhamnetin glycoside content, in a mouse model of diet-induced obesity and isolated pancreatic islets. Mice fed a high-fat (HF) diet supplemented with OFI extract exhibited reduced body weight gain and lower levels of circulating total cholesterol, LDL cholesterol, and HDL cholesterol compared to those on the HF diet alone. Furthermore, HF-OFI diet-fed mice showed lower glucose and insulin concentrations but slightly higher insulin levels compared to control mice. These metabolic enhancements were associated with decreased adipocyte size, enhanced hepatic phosphorylation of IRS1 tyr-608 and S6 K thr-389, and reduced hepatic lipid content (70). Jamali-Raeufy et al. studied the effect of Isorhamnetin on diabetic male rats for 12 weeks. Isorhamnetin, administered intraperitoneally at a dose of 10 mg/kg body weight once daily, elicited significant effects on various parameters. Notably, Isorhamnetin treatment led to a marked reduction in pain severity and blood glucose levels, while also promoting a significant increase in body weight compared to the control group. Moreover, Isorhamnetin demonstrated inhibitory effects on astroglial activation, acetyl-cholinesterase activity, oxidative stress markers, apoptosis, and inflammatory markers within diabetic rats (58).

In this study, we incorporated molecular biomarkers (mRNAs-miRNAs-LncRNAs) expressed in livers and adipose tissues of animal models representing normal, T2DM, and treated groups alongside conventional biochemical parameters. Our objective was to utilize these predictive targets to select the most potent candidates for treating T2DM using medicinal plant-based drugs such as Rosavin, isorhamnetin, and Caffeic acid, as well as probiotics like Z-biotic. To achieve this goal, we developed a machine-learning model using 5 different algorithms (KNN, RF, LR, LGBM, and Ada Boost) to recognize the significant features associated with T2DM to achieve a reliable and effective improvement prediction.

As expected, the diabetic rat models exhibited the highest levels of serum insulin resistance index (HOMA-IR), insulin, fasting blood glucose, as well as biomarkers indicating renal function impairment, liver damage, and lipid profiles. Conversely, the normal group, as well as the groups treated with medicinal plants and probiotics, displayed the lowest levels of these markers, more obvious in the highest drug doses, with an opposite relationship observed for HDL (Table 4). As a heterogeneous condition, T2DM is typified by impaired insulin secretion (known as the beta cell secretory defect) and insulin resistance, leading to elevated blood glucose levels. Various disturbances in biochemical parameters are associated with T2DM, including dyslipidemia, which involves elevated fasting and postprandial TG, decreased levels of HDL-C, increased levels of LDL-C, and a prevalence of small, dense LDL-C particles (71). Liver enzymes play a crucial role in regulating metabolism, particularly in maintaining normal blood glucose levels during fasting and after meals. Insulin resistance in the liver leads to increased glycogenolysis and lipolysis, resulting in elevated hepatic glucose production and abnormal triglyceride and fatty acid synthesis. These abnormalities, including elevated levels of AST, ALT, and ALP, are indicative of liver dysfunction and can precede the onset of fasting hyperglycemia. The elevation of these enzymes may be attributed to the direct hepatotoxic effects of excess fatty acids, which can disrupt cell membranes and impair mitochondrial function (72). Additionally, oxidative stress, peroxisomal beta-oxidation, and inflammation mediated by pro-inflammatory cytokines like TNF-α contribute to hepatic injury. Furthermore, the increased activity of ALT, a gluconeogenic enzyme suppressed by insulin, suggests a disruption in insulin signaling rather than solely hepatocyte injury (73, 74). Prolonged high blood glucose levels lead to increased oxidative stress, inflammation, and dysfunction of microvascular endothelial cells, which are common underlying mechanisms in both diabetic nephropathy (DN) and diabetic retinopathy (DR). These conditions often coexist, with elevated blood urea nitrogen (BUN) levels indicating their presence. BUN levels are also linked to the body’s catabolic activity and may reflect reduced blood flow and increased oxidative stress. Additionally, microvascular hypoperfusion and oxidative stress contribute to the development of DR, potentially explaining the association between elevated BUN levels and DR (75). Serum creatinine, primarily metabolized by skeletal muscle, correlates with total skeletal muscle mass. Low serum creatinine levels are considered a risk factor for T2DM and dysglycemia, as they reflect reduced skeletal muscle mass. Skeletal muscle plays a crucial role in glucose uptake, with reduced uptake contributing to insulin resistance and the development of T2DM long before hyperglycemia becomes apparent. The decrease in skeletal muscle mass is thus associated with increased insulin resistance and the risk of developing T2DM (76, 77).

We next determined the histopathological alterations in the studied groups, the livers of the T2DM-induced rats showed significant changes, including fat accumulation, cell death, inflammation, and abnormal cell morphology, in addition to an increase in the adipocyte cell size. Treatment with Rosavin at a higher dose restored the liver structure to normal, while Z-Biotic and isorhamnetin showed moderate improvements. caffeic acid had some positive effects but to a lesser extent. The administration of rosavin, particularly at a high dosage, significantly decreased the severity of liver inflammation grades. This outcome aligns with our previous findings in a rat model of non-alcoholic steatohepatitis (NASH), where rosavin treatment improved liver functions, and lipid profile, and mitigated hepatic inflammation, fibrosis, and cell death (78). Overall, the drugs Rosavin, Z-Biotic, isorhamnetin, and Caffeic acid exhibited varying effects on the destructive liver architecture induced by diabetes, with higher doses generally demonstrating a more restorative effect and this was also obvious in decreasing the adipocyte cell diameter and liver inflammation grades.

Inflammation and autoimmunity play significant roles in the development of diabetes, and hence, targeting the inflammatory response has shown therapeutic benefits (79, 80). The NOD signaling pathway is involved in the inflammation triggered by the cGAS-STING pathway (81, 82). Perturbations in the gut microbiota in diabetic patients can, in turn, activate the cGAS-STING-NOD pathway, leading to inflammation (83). Paradoxically, we found that the probiotics Z-Biotic 1mg significantly modulated most of the inflammatory-enriched mRNAs (DDX58, NFKB1, CHUK, RET) that were implicated in the STING/NOD/IR pathways and made them decrease significantly to be close to the normal group ratios, interestingly, the high dose of caffeic acid behaved similarly.

The cytoplasmic DNA sensor known as cyclic GMP-AMP synthase (cGAS) directly binds to DNA, including mitochondrial DNA (mtDNA) and abnormal bacterial or viral DNA, leading to its activation. Upon activation, cGAS triggers the synthesis of a secondary messenger called 2′,3′-cGAMP, which then binds to and activates the stimulator of the interferon gene (STING). STING, a transmembrane protein primarily found in the endoplasmic reticulum (ER), undergoes translocation to the Golgi and recruits the downstream TANK-binding kinase 1 (TBK1) to form a complex. This complex, in turn, phosphorylates and activates interferon regulatory factor 3 (IRF3) and nuclear factor kappa B (NF-κB), initiating a cascade of signals that activate various innate immune-related genes, including type I interferon (IFN) (84). Recent studies have implicated the cGAS-STING pathway in the development of diabetic cardiomyopathy, a condition characterized by aseptic inflammatory activation. Notably, increased mtDNA in the cytosol and the activation of the cGAS-STING signaling pathway, along with its downstream targets such as IRF3, NF-κB, IL-18, and IL-1β, have been observed in the context of diabetic cardiomyopathy (85). A previous study discovered that activating STING resulted in an induction of the NLRP3 inflammasome and pyroptosis, which culminated in an inflammatory response observed in diabetic mice (86). Qiao et al. found that STING plays a unique role in regulating insulin action in peripheral metabolic tissues and insulin secretion from β-cells (87).DDX58 participates in innate immune responses through the NOD-like receptor signaling pathway, detecting cytoplasmic DNA and triggering the production of interferon 1 and inflammatory cytokines (88, 89). Moreover, the cGAS-STING1 and DDX58-MAVS pathways are connected to the innate immune response (90).In a study conducted by Frietze et al. using a NASH model, it was shown that activation of DDX58 leads to the cleavage of LC3, a marker linked to autophagosome formation, DDX58 interacts with the autophagy receptor protein SQSTM1 to degrade itself selectively after viral stimulation. Activation of DDX58 enhances autophagic responses, aiding in the removal of harmful lipid inclusion bodies associated with inflammation and cell death. Excessive fatty acids hinder DDX58 activity, reducing crucial autophagic responses and worsening lipotoxicity. The study revealed that DDX58 directly influences SQSTM1 gene expression, protein accumulation, and targeted autophagic degradation. Furthermore, sustained overexpression of DDX58 markedly reduces inflammation mediated by JAK-STAT signaling (91).Studies have demonstrated that ZBP1, a cytosolic nucleic acid sensor, plays a pivotal role in coordinating innate immune responses by triggering both NF-κB and interferon regulatory factors (IRFs) signaling pathways (92). Additionally, ZBP1 activation leads to the expression of inflammatory cytokines and induces various forms of inflammatory cell death, including necroptosis, pyroptosis, apoptosis, and PANoptosis, in response to different host-derived nucleic acids (93, 94). Moreover, ZBP1 acts as a regulator of IFN-1-mediated disease progression, sensing mitochondrial DNA instability through the cGAS-STING pathway and sustaining IFN-1 signaling, which is involved in heart failure and cardiac cell remodeling (95, 96).

Stimulation of NOD and TNF receptors in intestinal epithelial cells activates CHUK, which stabilizes ATG16L1, preventing endoplasmic reticulum stress during inflammation (97). It also inhibits Kappa Beta kinases that suppress NFKB1 (98). High glucose promotes proliferation and invasiveness in pancreatic cancer cells by upregulating RET, a proto-oncogene encoding a receptor tyrosine kinase (99). ZBP1 acts as a regulator of IFN-1-mediated disease progression, sensing mitochondrial DNA instability through the cGAS-STING pathway and sustaining IFN-1 signaling, which is involved in heart failure and cardiac cell remodeling (95, 96) A previous study discovered that activating STING resulted in an induction of the NLRP3 inflammasome and pyroptosis, which culminated in an inflammatory response observed in diabetic mice (86). Qiao et al. found that STING plays a unique role in regulating insulin action in peripheral metabolic tissues and insulin secretion from β-cells (87).

Intense evidence supports the involvement of miRNAs in the intricate regulation of glucose homeostasis, making them potential contributors to the development of diabetes (100). These small molecules exert their influence by impacting various aspects of insulin levels, including insulin production, and exocytosis (101) Additionally, lncRNAs can act as miRNA sponges by binding to miRNAs, thus preventing them from interacting with their mRNA targets (102). This regulatory mechanism plays a crucial role in diabetes pathogenesis, influencing beta-cell function, apoptosis, insulin secretion, glucose metabolism, and insulin resistance (103, 104).

We previously reported that miRNAs (has-miR-1976 and has-miR-611) acted as sponges on lncRNAs (AC074117.2 and RP4–605O3.4), decreasing the expression of mRNA (CHUK mRNA). Notably, these findings demonstrated the potential of has-miR(-1976 and 611) miRNAs as distinguishing factors between insulin-resistant and insulin-sensitive patients, and their involvement in modulating the STING/NOD/IR pathways (105). These results were consistent with our observations in animal models. Previous studies have highlighted the role of miR-1976 as a tumor suppressor in non-small cell lung cancer (106). In vitro experiments demonstrated that the suppression of miR-1976 resulted in a significant acceleration of wound healing and cell migration. Furthermore, it stimulated cell proliferation, decreased cell apoptosis, and increased the populations of CD44+/CD24− cells (107).

Our results suggested that RP11–773H22.4 could bind to miR-1, and miR-3163 and serve as a competing endogenous RNA (ceRNA) to upregulate m-TOR, IGF1-R, and downregulate AKT2 expression, thus contributing to increased autophagy and the progression of T2DM, consistent with our previous results (28). Previous studies have indicated that insulin-like growth factor-1 (IGF-1) hinders the process of autophagy by inhibiting AKT. This inhibition is believed to be mediated by AKT’s activation of rapamycin complex 1 (mTORC1), a known inhibitor of autophagy, there for Inhibition of IGF-1R signaling cascade reduces autophagy, impacting autophagosome precursor formation and suggests that targeting the IGF-1R receptor or its downstream pathway may have implications for therapeutic purposes (108). Moreover, Akt2 can inhibit the expression of the cGAS-STING pathway and suppress the inflammatory response in diabetes (109). Autophagy is crucial for maintaining cellular balance during stressful conditions, excessive or uncontrolled autophagy can trigger autophagy-dependent apoptosis, cardiac injury, and impaired function, primarily through autophagic cell death (110112). MiR-1, which has been studied extensively, is often suppressed in various biological samples from patients with T2DM (113). Overexpression of miR-1 has been shown to inhibit cardiac fibrosis and apoptosis by altering the expression levels of Bcl-2, TGFb1, and Bax (114). Additionally, mTOR plays a significant role in cardioprotection, diabetes, cellular metabolism, apoptosis, autophagy, and mitochondrial biogenesis (115117) Studies have demonstrated that inhibiting mTOR with rapamycin (RAPA) can reduce cardiomyocyte apoptosis following myocardial infarction (MI) stress (118). Moreover, in retinoblastoma cancer stem cells, miR-3163 influences cell proliferation, apoptosis, and drug resistance (119). Additionally, miR-3163 has been shown to enhance the sensitivity of hepatocellular carcinoma (HCC) cells to sorafenib by suppressing the cleavage of Notch protein (120).

The present study involved; i, Data preprocessing and feature selection: The dataset was carefully preprocessed to ensure data integrity and consistency with the study objectives. Null values ​​and duplicate entries were carefully handled to maintain the quality of the dataset, which was essential for subsequent analysis. Targeted variables created following data correction methods ‘standardization’ are needed to predict response to type 2 diabetes mellitus (T2DM) treatment Based on ‘fat cell diameter’ and ‘. Inflammation (liver)’ features used to divide the sample into “unimproved” and “improved” categories. Categorical features of the dataset were adjusted to facilitate their incorporation into the modeling pipeline. Robust coding techniques such as one-hot coding were used to ensure compatibility with machine learning algorithms. Feature selection was performed using cross-validation (RFECV) method and recursive feature elimination, which is the cornerstone of predictive modeling. This approach systematically identified factors that contribute to understanding T2DM drug response patterns. Notably, the selected molecules include ZBP1 (A), STING1(L), DDX58 (L), mTOR (L), NFKB1 (A), CHUK (A), RET (L), RET (A), miR-1. (A), miR-611 (L), miR-611 (A), LncRNA-RP11–773H22.4 (L), and LncRNA-RP11–773H22.4 (A). Furthermore, biological parameters such as glucose, insulin, HOMA-IR, total cholesterol, triglycerides, AST, ALT, creatinine, BUN, and ACR were identified as important contributors to the prediction of T2DM treatment response. ii, Training and evaluation model: The trained machine learning models, including K-Nearest Neighbors (KNN) and Light Gradient Boosting Machine (LGBM), were subjected to rigorous evaluation to assess their predictive prowess. Leveraging cross-validation techniques, the models demonstrated commendable performance, achieving an impressive accuracy of 80%. These findings underscore the robustness and efficacy of the selected features in discerning intricate patterns inherent in T2DM treatment response. Moreover, the systematic approach adopted in this study holds promise for enhancing therapeutic strategies and optimizing patient outcomes in clinical settings.

In this study, a feature selection approach based on the RFECV technique was adopted to select the top molecular and biochemical features that had the highest prediction accuracy for T2DM drug response in the dataset. The selected molecular features were ZBP1 (A), STING1 (L), DDX58 (L), mTOR (L), NFKB1 (A), CHUK (A), RET (L), RET (A), miR-1976 (A), miR-611 (L), miR-611 (A), LncRNA-RP11–773H22.4 (L), and LncRNA-RP11–773H22.4 (A). While the biochemical features were Glucose, Insulin, HOMA-IR, Total cholesterol, Triglycerides, AST, ALT, Creatinine, BUN, and ACR. The utilized molecular and biochemical models using these features combined or separately maintained decent prediction performance with top classifiers KNN and LGBM at an accuracy of 80%.

While our study provides valuable insights into the potential therapeutic efficacy of medicinal plant-based drugs and probiotics in treating T2DM, there are several limitations that should be acknowledged. Our study focused on a limited number of medicinal plant-based drugs and a single probiotic, which may not fully represent the spectrum of potential therapeutic interventions for T2DM. Future research incorporating a broader range of pharmacological agents and therapeutic modalities could provide a more comprehensive understanding of effective treatment strategies for this complex condition. Further larger studies addressing the possibility of multiple dose and long-term treatment regimens and their impact on the evaluation of drug safety and efficacy should be explored. Also, Further future validation through independent datasets is strongly needed to ensure the findings’ applicability in human T2DM treatment. Moreover, while our study assessed various molecular and biochemical features associated with T2DM pathogenesis, it is important to recognize that these chosen parameters may not capture the full complexity of the disease process. Additional factors such as genetic predisposition, environmental influences, and lifestyle factors may also contribute to individual variability in drug response and treatment outcomes. Further research is needed to address the limitations and translate these findings into clinically meaningful interventions for patients with T2DM. Researchers stressed the pivotal role of all members of brain, kidneys, pancreatic cells, alpha cells, and the gastrointestinal tract in the development of glucose intolerance besides the four cell types(liver, muscle and adipose tissues) (121). Thus further larger validation in other tissues is strongly recommended to ensure the generalization of the results.

6 Conclusion

We developed a prediction system for identification of potential therapeutic targets using machine learning algorithms with feature selection using mRNAs-miRNAs-LncRNAs implicated in autophagy and STING/NOD/IR pathways that directly correlated with T2DM pathogenesis in addition to biochemical features. Our results demonstrated that the KNN algorithm outperformed other classifiers in both the molecular and combined models, for the biochemical model, the LGBM Classifier exhibited the highest performance with AUC values of 0.9177, 0.9, and 0.9397, respectively. Notably, our machine learning approach successfully identified 20 significant features out of the total 44 features in the combined model with an accuracy of 80%.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics statement

The animal study was approved by Research Ethics Committee of the Faculty of Medicine at Ain Shams University under federal-wide assurance NO. FWA000017585. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

MS: Writing – review & editing, Validation, Methodology, Investigation, Funding acquisition, Conceptualization. HA-A: Writing – review & editing, Investigation. AK: Writing – original draft, Visualization, Data curation. RK: Writing – original draft, Investigation, Formal analysis. MR: Writing – review & editing, Software. MA: Writing – review & editing, Investigation. GD: Writing – original draft, Methodology. ME: Writing – original draft, Visualization, Methodology. RE: Writing – review & editing, Supervision. HA: Writing – review & editing, Investigation. EA: Writing – original draft, Data curation. DE: Writing – review & editing, Investigation. HK: Writing – review & editing, Investigation. MF: Writing – original draft, Writing – review & editing, Methodology. HE: Writing – review & editing, Investigation. LF: Writing – review & editing, Validation. MBA: Writing – review & editing, Validation. EH: Writing – review & editing, Visualization, Resources. HF: Writing – review & editing, Validation. LS: Writing – review & editing, Validation. IA: Writing – review & editing, Writing – original draft, Methodology, Data curation.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study received fund from the Academy of Scientific Research and Technology, Egypt, JESOR call 2019 ID 5090.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1384984/full#supplementary-material

References

1. Watkins DA, Ali MK. Measuring the global burden of diabetes: implications for health policy, practice, and research. Lancet. (2023) 402(10397):163–16. doi: 10.1016/S0140-6736(23)01287-4

PubMed Abstract | CrossRef Full Text | Google Scholar

2. American Diabetes Association. Standards of medical care in diabetes—2022 abridged for primary care providers. Clin Diabetes. (2022) 40:10–38. doi: 10.2337/cd22-as01

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ray KK, Nicholls SJ, Li N, Louie MJ, Brennan D, Lincoff AM, et al. Efficacy and safety of bempedoic acid among patients with and without diabetes: prespecified analysis of the CLEAR Outcomes randomised trial. Lancet Diabetes Endocrinology. (2024) 12:19–28. doi: 10.1016/S2213-8587(23)00316-9

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Antonio-Villa NE, Fernández-Chirino L, Vargas-Vazquez A, Bahena-López JP, Bello-Chavolla OY. Trends in diabetes subgroups and their risk for all-cause, cardiovascular disease and diabetes-specific mortality in the US: A data-driven reproducible machine learning approach. J Endocrine Soc. (2021) 5:A423. doi: 10.1210/jendso/bvab048.863

CrossRef Full Text | Google Scholar

5. Buse JB, Wexler DJ, Tsapas A, Rossing P, Mingrone G, Mathieu C, et al. 2019 update to: management of hyperglycemia in type 2 diabetes, 2018. A consensus report by the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetes Care. (2020) 43:487–93. doi: 10.2337/dci19-0066

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Rayego-Mateos S, Morgado-Pascual JL, Opazo-Ríos L, Guerrero-Hue M, García-Caballero C, Vázquez-Carballo C, et al. Pathogenic pathways and therapeutic approaches targeting inflammation in diabetic nephropathy. Int J Mol Sci. (2020) 21:3798. doi: 10.3390/ijms21113798

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gong J, Gao X, Ge S, Li H, Wang R, Zhao L. The role of cGAS-STING signalling in metabolic diseases: from signalling networks to targeted intervention. Int J Biol Sci. (2024) 20:152. doi: 10.7150/ijbs.84890

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Tsalamandris S, Antonopoulos AS, Oikonomou E, Papamikroulis GA, Vogiatzi G, Papaioannou S, et al. The role of inflammation in diabetes: current concepts and future perspectives. Eur Cardiol review. (2019) 14:50. doi: 10.15420/ecr

CrossRef Full Text | Google Scholar

9. Watada H, Fujitani Y. Minireview: autophagy in pancreatic β-cells and its implication in diabetes. Mol Endocrinology. (2015) 29:338–48. doi: 10.1210/me.2014-1367

CrossRef Full Text | Google Scholar

10. Lim YM, Lim H, Hur KY, Quan W, Lee HY, Cheon H, et al. Systemic autophagy insufficiency compromises adaptation to metabolic stress and facilitates progression from obesity to diabetes. Nat Commun. (2014) 5:4934. doi: 10.1038/ncomms5934

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sakshi S, Jayasuriya R, Ganesan K, Xu B, Ramkumar KM. Role of circRNA-miRNA-mRNA interaction network in diabetes and its associated complications. Mol Therapy-Nucleic Acids. (2021) 26:1291–302. doi: 10.1016/j.omtn.2021.11.007

CrossRef Full Text | Google Scholar

12. Yang Q, Fang D, Chen J, Hu S, Jiang J, Zeng M, et al. LncRNAs associated with oxidative stress in diabetic wound healing: regulatory mechanisms and application prospects. Theranostics. (2023) 13:3655. doi: 10.7150/thno.85823

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Pandey A, Ajgaonkar S, Jadhav N, Saha P, Gurav P, Panda S, et al. Current Insights into miRNA and lncRNA Dysregulation in Diabetes: Signal Transduction, Clinical Trials and Biomarker Discovery. Pharmaceuticals. (2022) 15:1269. doi: 10.3390/ph15101269

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Smit-McBride Z, Morse LS. MicroRNA and diabetic retinopathy—biomarkers and novel therapeutics. Ann Trans Med. (2021) 9. doi: 10.21037/atm

CrossRef Full Text | Google Scholar

15. Eliasson L, Regazzi R. Micro (RNA) management and mismanagement of the islet. J Mol Biol. (2020) 432:1419–28. doi: 10.1016/j.jmb.2019.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Erfanian Omidvar M, Ghaedi H, Kazerouni F, Kalbasi S, Shanaki M, Miraalamy G, et al. Clinical significance of long noncoding RNA VIM-AS1 and CTBP1-AS2 expression in type 2 diabetes. J Cell Biochem. (2019) 120:9315–23. doi: 10.1002/jcb.28206

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Tello-Flores VA, Beltrán-Anaya FO, Ramírez-Vargas MA, Esteban-Casales BE, Navarro-Tito N, Alarcón-Romero LD, et al. Role of long non-coding RNAs and the molecular mechanisms involved in insulin resistance. Int J Mol Sci. (2021) 22:7256. doi: 10.3390/ijms22147256

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shuey MM, Lee KM, Keaton J, Khankari NK, Breeyear JH, Walker VM, et al. A genetically supported drug repurposing pipeline for diabetes treatment using electronic health records. EBioMedicine. (2023) 94. doi: 10.1016/j.ebiom.2023.104674

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kavishankar GB, Lakshmidevi N, Murthy SM, Prakash HS, Niranjana SR. Diabetes and medicinal plants-A review. Int J Pharm BioMed Sci. (2011) 2:65–80.

Google Scholar

20. Akbari V, Hendijani F. Effects of probiotic supplementation in patients with type 2 diabetes: systematic review and meta-analysis. Nutr Rev. (2016) 74:774–84. doi: 10.1093/nutrit/nuw039

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Biassoni R, Di Marco E, Squillario M, Barla A, Piccolo G, Ugolotti E, et al. Gut microbiota in T1DM-onset pediatric patients: machine-learning algorithms to classify microorganisms as disease linked. J Clin Endocrinol Metab. (2020) 105:e3114–26. doi: 10.1210/clinem/dgaa407

CrossRef Full Text | Google Scholar

22. Lv C, Sun Y, Zhang ZY, Aboelela Z, Qiu X, Meng ZX. β-cell dynamics in type 2 diabetes and in dietary and exercise interventions. J Mol Cell Biol. (2022) 14:mjac046. doi: 10.1093/jmcb/mjac046

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gothai S, Ganesan P, Park SY, Fakurazi S, Choi DK, Arulselvan P. Natural phyto-bioactive compounds for the treatment of type 2 diabetes: inflammation as a target. Nutrients. (2016) 8:461. doi: 10.3390/nu8080461

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ayesha IE, Monson NR, Klair N, Patel U, Saxena A, Patel D, et al. Probiotics and their role in the management of type 2 diabetes mellitus (Short-term versus long-term effect): A systematic review and meta-analysis. Cureus. (2023) 15:e46741. doi: 10.7759/cureus.46741

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kalai FZ, Boulaaba M, Ferdousi F, Isoda H. Effects of isorhamnetin on diabetes and its associated complications: A review of in vitro and in vivo studies and a post hoc transcriptome analysis of involved molecular pathways. Int J Mol Sci. (2022) 23:704. doi: 10.3390/ijms23020704

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Qiu S, Sun G, Zhang Y, Li X, Wang R. Involvement of the NF-κB signaling pathway in the renoprotective effects of isorhamnetin in a type 2 diabetic rat model. Biomed Rep. (2016) 4:628–34. doi: 10.3892/br.2016.636

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Khaled R. Biological activities of isorhamnetin: A review. Plantae Scientia. (2020) 3:78–81. doi: 10.32439/ps.v3i5.78-81

CrossRef Full Text | Google Scholar

28. Matboli M, Saad M, Hasanin AH, Saleh LA, Baher W, Bekhet MM, et al. New insight into the role of isorhamnetin as a regulator of insulin signaling pathway in type 2 diabetes mellitus rat model: Molecular and computational approach. Biomedicine Pharmacotherapy. (2021) 135:111176. doi: 10.1016/j.biopha.2020.111176

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Li J, Yuan H, Zhao Z, Li L, Li X, Zhu L, et al. The mitigative effect of isorhamnetin against type 2 diabetes via gut microbiota regulation in mice. Front Nutr. (2022) 9:1070908. doi: 10.3389/fnut.2022.1070908

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Bernatoniene J, Jakstas V, Kopustinskiene DM. Phenolic compounds of Rhodiola rosea L. as the potential alternative therapy in the treatment of chronic diseases. Int J Mol Sci. (2023) 24:12293. doi: 10.3390/ijms241512293

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zheng T, Bian F, Chen L, Wang Q, Jin S. Beneficial effects of Rhodiola and salidroside in diabetes: potential role of AMP-activated protein kinase. Mol Diagnosis Ther. (2019) 23:489–98. doi: 10.1007/s40291-019-00402-4

CrossRef Full Text | Google Scholar

32. Wang ZS, Gao F, Lu FE. Effect of ethanol extract of Rhodiola rosea on the early nephropathy in type 2 diabetic rats. J Huazhong Univ Sci Technol [Medical Sciences]. (2013) 33:375–8. doi: 10.1007/s11596-013-1127-6

CrossRef Full Text | Google Scholar

33. Dib I, El Alaoui-Faris FE. Artemisia campestris L.: review on taxonomical aspects, cytogeography, biological activities and bioactive compounds. Biomedicine Pharmacotherapy. (2019) 109:1884–906. doi: 10.1016/j.biopha.2018.10.149

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Hbika A, Daoudi NE, Bouyanzer A, Bouhrim M, Mohti H, Loukili EH, et al. Artemisia absinthium L. Aqueous and ethyl acetate extracts: Antioxidant effect and potential activity in vitro and in vivo against pancreatic α-amylase and intestinal α-glucosidase. Pharmaceutics. (2022) 14:481. doi: 10.3390/pharmaceutics14030481

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Olennikov DN, Chirikova NK, Kashchenko NI, Nikolaev VM, Kim SW, Vennos C. Bioactive phenolics of the genus Artemisia (Asteraceae): HPLC-DAD-ESI-TQ-MS/MS profile of the Siberian species and their inhibitory potential against α-amylase and α-glucosidase. Front Pharmacol. (2018) 9:756. doi: 10.3389/fphar.2018.00756

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ma J, Lyu Y, Liu X, Jia X, Cui F, Wu X, et al. Engineered probiotics. Microbial Cell factories. (2022) 21:1–9. doi: 10.1186/s12934-022-01799-0

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Appala Naidu B, Kannan K, Santhosh Kumar DP, Oliver JW, Abbott ZD. Lyophilized B. subtilis ZB183 spores: 90-Day Repeat Dose Oral (Gavage) toxicity study in wistar rats. J Toxicol. (2019) 2019. doi: 10.1101/724542

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Zou Q, Qu K, Luo Y, Yin D, Ju Y, Tang H. Predicting diabetes mellitus with machine learning techniques. Front Genet. (2018) 9:515. doi: 10.3389/fgene.2018.00515

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Woldaregay AZ, Årsand E, Botsis T, Albers D, Mamykina L, Hartvigsen G. Data-driven blood glucose pattern classification and anomalies detection: machine-learning applications in type 1 diabetes. J Med Internet Res. (2019) 21:e11030. doi: 10.2196/11030

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Kavakiotis I, Tsave O, Salifoglou A, Maglaveras N, Vlahavas I, Chouvarda I. Machine learning and data mining methods in diabetes research. Comput Struct Biotechnol J. (2017) 15:104–16. doi: 10.1016/j.csbj.2016.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Deberneh HM, Kim I. Prediction of type 2 diabetes based on machine learning algorithm. Int J Environ Res Public Health. (2021) 18:3317. doi: 10.3390/ijerph18063317

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Wei S, Zhao X, Miao C. A comprehensive exploration to the machine learning techniques for diabetes identification. IEEE 4th World Forum on Internet of Thing (2018). p. 291–5.

Google Scholar

43. Elsherbini AM, Alsamman AM, Elsherbiny NM, El-Sherbiny M, Ahmed R, Ebrahim HA, et al. Decoding diabetes biomarkers and related molecular mechanisms by using machine learning, text mining, and gene expression analysis. Int J Environ Res Public Health. (2022) 19:13890. doi: 10.3390/ijerph192113890

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Xu Q, Xiong Y, Dai H, Kumari KM, Xu Q, Ou HY, et al. PDC-SGB: Prediction of effective drug combinations using a stochastic gradient boosting algorithm. J Theor Biol. (2017) 417:1–7. doi: 10.1016/j.jtbi.2017.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Jang IS, Neto EC, Guinney J, Friend SH, Margolin AA. Systematic assessment of analytical methods for drug sensitivity prediction from cancer cell line data. Pac. Symp Biocomput. (2014) 19:63–74.

Google Scholar

46. Costello JC, Heiser LM, Georgii E, Gönen M, Menden MP, Wang NJ, et al. A community effort to assess and improve drug sensitivity prediction algorithms. Nat Biotechnol. (2014) 32:1202–12. doi: 10.1038/nbt.2877

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Batts KP, Ludwig J. An update on terminology and reporting. Am J Surg pathology. (1995) 19:1409–17. doi: 10.1097/00000478-199512000-00007

CrossRef Full Text | Google Scholar

48. Ochoa-Morales PD, González-Ortiz M, Martínez-Abundis E, Pérez-Rubio KG, Patiño-Laguna AD. Anti-hyperglycemic effects of propolis or metformin in type 2 Diabetes Mellitus. Int J Vitamin Nutr Res. (2022) 72:59–139. doi: 10.1024/0300-9831/a000760

CrossRef Full Text | Google Scholar

49. García S, Luengo J, Herrera F. Data preprocessing in data mining. Cham, Switzerland: Springer International Publishing (2015).

Google Scholar

50. Yin Z, Zhang J. Operator functional state classification using least-square support vector machine based recursive feature elimination technique. Comput Methods programs biomedicine. (2014) 113:101–15. doi: 10.1016/j.cmpb.2013.09.007

CrossRef Full Text | Google Scholar

51. Kohavi R. A study of cross-validation and bootstrap for accuracy estimation and model selection. InIjcai. (1995) 14:1137–45. doi: 10.1038/s41581-020-0278-5

CrossRef Full Text | Google Scholar

52. Cole JB, Florez JC. Genetics of diabetes mellitus and diabetes complications. Nat Rev nephrology. (2020) 16:377–90. doi: 10.1038/s41581-020-0278-5

CrossRef Full Text | Google Scholar

53. Li Y, Pham V, Bui M, Song L, Wu C, Walia A, et al. Rhodiola rosea L.: an herb with anti-stress, anti-aging, and immunostimulating properties for cancer chemoprevention. Curr Pharmacol Rep. (2017) 3:384–95. doi: 10.1007/s40495-017-0106-1

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Pomari E, Stefanon B, Colitti M. Effects of two different Rhodiola rosea extracts on primary human visceral adipocytes. Molecules. (2015) 20:8409–28. doi: 10.3390/molecules20058409

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Verpeut JL, Walters AL, Bello NT. Citrus aurantium and Rhodiola rosea in combination reduce visceral white adipose tissue and increase hypothalamic norepinephrine in a rat model of diet-induced obesity. Nutr Res. (2013) 33:503–12. doi: 10.1016/j.nutres.2013.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Mao Y. Hypoglycemic and hypolipidaemic activities of polysaccharides from Rhodiola rosea in KKAy mice. J Food Process preservation. (2017) 41:e13219. doi: 10.1111/jfpp.2017.41.issue-6

CrossRef Full Text | Google Scholar

57. Liu QH, Zhang K, Feng SS, Zhang LJ, Li SY, Wang HY, et al. Rosavin alleviates LPS-induced acute lung injure by modulating the TLR-4/NF-κB/MAPK singnaling pathways. Int J Mol Sci. (2024) 25:1875. doi: 10.3390/ijms25031875

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Jamali-Raeufy N, Baluchnejadmojarad T, Roghani M. Isorhamnetin exerts neuroprotective effects in STZ-induced diabetic rats via attenuation of oxidative stress, inflammation and apoptosis. J Chem Neuroanatomy. (2019) 102:101709. doi: 10.1016/j.jchemneu.2019.101709

CrossRef Full Text | Google Scholar

59. Zhang XX, Jiang CH, Liu Y, Lou DX, Huang YP, Gao M, et al. Cyclocarya paliurus triterpenic acids fraction attenuates kidney injury via AMPK-mTOR-regulated autophagy pathway in diabetic rats. Phytomedicine. (2019) 64:153060. doi: 10.1016/j.phymed.2019.153060

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Zhou W, Zhou X, Zhang Y, Wang Y, Wu W, Zhang H, et al. A network pharmacology-based approach for exploring key active compounds and pharmacological mechanisms of tangshen formula for treatment of diabetic nephropathy. J Diabetes Res. (2021) 2021:1–7. doi: 10.1155/2021/8833688

CrossRef Full Text | Google Scholar

61. Yang JH, Kim SC, Kim KM, Jang CH, Cho SS, Kim SJ, et al. Isorhamnetin attenuates liver fibrosis by inhibiting TGF-β/Smad signaling and relieving oxidative stress. Eur J Pharmacol. (2016) 783:92–102. doi: 10.1016/j.ejphar.2016.04.042

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Lu X, Liu T, Chen K, Xia Y, Dai W, Xu S, et al. Isorhamnetin: A hepatoprotective flavonoid inhibits apoptosis and autophagy via P38/PPAR-α pathway in mice. Biomedicine Pharmacotherapy. (2018) 103:800–11. doi: 10.1016/j.biopha.2018.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Akhlaghipour I, Shad AN, Askari VR, Maharati A, Rahimi VB. How caffeic acid and its derivatives combat diabetes and its complications: A systematic review. J Funct Foods. (2023) 110:105862. doi: 10.1016/j.jff.2023.105862

CrossRef Full Text | Google Scholar

64. Xu W, Luo Q, Wen X, Xiao M, Mei Q. Antioxidant and anti-diabetic effects of caffeic acid in a rat model of diabetes. Trop J Pharm Res. (2020) 19:1227–32. doi: 10.4314/tjpr.v19i6.17

CrossRef Full Text | Google Scholar

65. Salem AM, Ragheb AS, Hegazy MG, Matboli M, Eissa S. Caffeic acid modulates miR-636 expression in diabetic nephropathy rats. Indian J Clin Biochem. (2019) 34:296–303. doi: 10.1007/s12291-018-0743-0

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Bhattacharya S, Oksbjerg N, Young JF, Jeppesen PB. Caffeic acid, naringenin and quercetin enhance glucose-stimulated insulin secretion and glucose sensitivity in INS-1E cells. Diabetes Obes Metab. (2014) 16:602–12. doi: 10.1111/dom.12236

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Lekomtseva Y, Zhukova I, Wacker A. Rhodiola rosea in subjects with prolonged or chronic fatigue symptoms: results of an open-label clinical trial. Complementary Med Res. (2017) 24:46–52. doi: 10.1159/000457918

CrossRef Full Text | Google Scholar

68. El-Sharkawy HM, Anees MM, Van Dyke TE. Propolis improves periodontal status and glycemic control in patients with type 2 diabetes mellitus and chronic periodontitis: a randomized clinical trial. J periodontology. (2016) 87:1418–26. doi: 10.1902/jop.2016.150694

CrossRef Full Text | Google Scholar

69. Chao PC, Hsu CC, Yin MC. Anti-inflammatory and anti-coagulatory activities of caffeic acid and ellagic acid in cardiac tissue of diabetic mice. Nutr Metab. (2009) 6:1–8. doi: 10.1186/1743-7075-6-33

CrossRef Full Text | Google Scholar

70. Rodríguez-Rodríguez C, Torres N, Gutiérrez-Uribe JA, Noriega LG, Torre-Villalvazo I, Leal-Díaz AM, et al. The effect of isorhamnetin glycosides extracted from Opuntia ficus-indica in a mouse model of diet induced obesity. Food Funct. (2015) 6:805–15. doi: 10.1039/C4FO01092B

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Shi R, Niu Z, Wu B, Zhang T, Cai D, Sun H, et al. Nomogram for the risk of diabetic nephropathy or diabetic retinopathy among patients with type 2 diabetes mellitus based on questionnaire and biochemical indicators: a cross-sectional study. Diabetes Metab Syndrome Obes. (2020) 20:1215–29. doi: 10.2147/DMSO.S244061

CrossRef Full Text | Google Scholar

72. Habte ML, Melka DS, Degef M, Menon MK, Yifter H, Feyisa TO. Comparison of lipid profile, liver enzymes, creatine kinase and lactate dehydrogenase among type II diabetes mellitus patients on statin therapy. Diabetes Metab Syndrome Obes. (2020) 18:763–73. doi: 10.2147/DMSO.S234382

CrossRef Full Text | Google Scholar

73. Atiba AS, Oparinde DP, Babatunde OA, Niran-Atiba T, Jimoh AK, Adepeju A. Liver enzymes and lipid profile among type 2 diabetic patients in Osogbo, Nigeria. Greener J Med Sci. (2013) 3:174–8. doi: 10.15580/GJMS

CrossRef Full Text | Google Scholar

74. Al-Jameil N, Khan FA, Arjumand S, Khan MF, Tabassum H. Associated liver enzymes with hyperlipidemic profile in type 2 diabetes patients. Int J Clin Exp pathology. (2014) 7:4345. doi: 10.1039/c4fo01092b

CrossRef Full Text | Google Scholar

75. Zhong JB, Yao YF, Zeng GQ, Zhang Y, Ye BK, Dou XY, et al. A closer association between blood urea nitrogen and the probability of diabetic retinopathy in patients with shorter type 2 diabetes duration. Sci Rep. (2023) 13:9881. doi: 10.1038/s41598-023-35653-z

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Takeuchi M, Imano H, Muraki I, Shimizu Y, Hayama-Terada M, Kitamura A, et al. Serum creatinine levels and risk of incident type 2 diabetes mellitus or dysglycemia in middle-aged Japanese men: a retrospective cohort study. BMJ Open Diabetes Res Care. (2018) 6:e000492. doi: 10.1136/bmjdrc-2017-000492

CrossRef Full Text | Google Scholar

77. Li R, He M, Yang Q, Liang Z, Li Y, Huang L, et al. Association between serum creatinine and type 2 diabetes in the Chinese population: a retrospective cohort study. Sci Rep. (2023) 13:6806. doi: 10.1038/s41598-023-33878-6

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Albadawy R, Hasanin AH, Agwa SH, Hamady S, Aboul-Ela YM, Raafat MH, et al. Rosavin ameliorates hepatic inflammation and fibrosis in the NASH rat model via targeting hepatic cell death. Int J Mol Sci. (2022) 23:10148. doi: 10.3390/ijms231710148

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Pollack RM, Donath MY, LeRoith D, Leibowitz G. Anti-inflammatory agents in the treatment of diabetes and its vascular complications. Diabetes Care. (2016) 39:S244–52. doi: 10.2337/dcS15-3015

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Bellucci PN, González Bagnes MF, Di Girolamo G, González CD. Potential effects of nonsteroidal anti-inflammatory drugs in the prevention and treatment of type 2 diabetes mellitus. J Pharm Practice. (2017) 30:549–56. doi: 10.1177/0897190016649551

CrossRef Full Text | Google Scholar

81. Murthy AM, Robinson N, Kumar S. Crosstalk between cGAS–STING signaling and cell death. Cell Death Differentiation. (2020) 27:2989–3003. doi: 10.1038/s41418-020-00624-8

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Mosallanejad K, Kagan JC. Control of innate immunity by the cGAS-STING pathway. Immunol Cell Biol. (2022) 100:409–23. doi: 10.1111/imcb.12555

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Luo Z, Ji Y, Gao H, Dos Reis FC, Bandyopadhyay G, Jin Z, et al. CRIg+ macrophages prevent gut microbial DNA–containing extracellular vesicle–induced tissue inflammation and insulin resistance. Gastroenterology. (2021) 160:863–74. doi: 10.1053/j.gastro.2020.10.042

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Luo L, An Y, Geng K, Wan S, Zhang F, Tan X, et al. High glucose-induced endothelial STING activation inhibits diabetic wound healing through impairment of angiogenesis. Biochem Biophys Res Commun. (2023) 668:82–9. doi: 10.1016/j.bbrc.2023.05.081

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Ma XM, Geng K, Law BY, Wang P, Pu YL, Chen Q, et al. Lipotoxicity-induced mtDNA release promotes diabetic cardiomyopathy by activating the cGAS-STING pathway in obesity-related diabetes. Cell Biol toxicology. (2023) 39:277–99. doi: 10.1007/s10565-021-09692-z

CrossRef Full Text | Google Scholar

86. Yang X, Chen Z, Luo Z, Yang D, Hao Y, Hu J, et al. STING deletion alleviates podocyte injury through suppressing inflammation by targeting NLRP3 in diabetic kidney disease. Cell Signalling. (2023) 15:110777. doi: 10.1016/j.cellsig.2023.110777

CrossRef Full Text | Google Scholar

87. Qiao J, Zhang Z, Ji S, Liu T, Zhang X, Huang Y, et al. A distinct role of STING in regulating glucose homeostasis through insulin sensitivity and insulin secretion. Proc Natl Acad Sci. (2022) 119:e2101848119. doi: 10.1073/pnas.2101848119

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Rong Y, Yang H, Xu H, Li S, Wang P, Wang Z, et al. Bioinformatic analysis reveals hub immune-related genes of diabetic foot ulcers. Front Surgery. (2022) 9:878965. doi: 10.3389/fsurg.2022.878965

CrossRef Full Text | Google Scholar

89. Yoneyama M, Kikuchi M, Natsukawa T, Shinobu N, Imaizumi T, Miyagishi M, et al. The RNA helicase RIG-I has an essential function in double-stranded RNA-induced innate antiviral responses. Nat Immunol. (2004) 5:730–7. doi: 10.1038/ni1087

PubMed Abstract | CrossRef Full Text | Google Scholar

90. Lei L, Bai Y, Fan Y, Li Y, Jiang H, Wang J. Comprehensive diagnostics of diabetic nephropathy by transcriptome RNA sequencing. Diabetes Metab Syndrome Obesity: Targets Ther. (2022) 1:3069–80. doi: 10.2147/DMSO.S371026

CrossRef Full Text | Google Scholar

91. Frietze KK, Brown AM, Das D, Franks RG, Cunningham JL, Hayward M, et al. Lipotoxicity reduces DDX58/Rig-1 expression and activity leading to impaired autophagy and cell death. Autophagy. (2022) 18:142–60. doi: 10.1080/15548627.2021.1920818

PubMed Abstract | CrossRef Full Text | Google Scholar

92. Song Q, Fan Y, Zhang H, Wang N. Z-DNA binding protein 1 orchestrates innate immunity and inflammatory cell death. Cytokine Growth Factor Rev. (2024) 109. doi: 10.1016/j.cytogfr.2024.03.005

CrossRef Full Text | Google Scholar

93. Shen J, San W, Zheng Y, Zhang S, Cao D, Chen Y, et al. Different types of cell death in diabetic endothelial dysfunction. Biomedicine Pharmacotherapy. (2023) 168:115802. doi: 10.1016/j.biopha.2023.115802

PubMed Abstract | CrossRef Full Text | Google Scholar

94. Song Q, Qi Z, Wang K, Wang N. Z-nucleic acid sensor ZBP1 in sterile inflammation. Clin Immunol. (2024) 10:109938. doi: 10.1016/j.clim.2024.109938

CrossRef Full Text | Google Scholar

95. Lei Y, Vanportfliet JJ, Chen YF, Bryant JD, Ragan KB, Tong CW, et al. ZBP1 sequesters cGAS in the cytoplasm and sustains type I interferon responses to mitochondrial DNA. BioRxiv. (2022) 30:2022–05. doi: 10.1038/ni1087

CrossRef Full Text | Google Scholar

96. Jiao H, Wachsmuth L, Kumari S, Schwarzer R, Lin J, Eren RO, et al. Z-nucleic-acid sensing triggers ZBP1-dependent necroptosis and inflammation. Nature. (2020) 580:391–5. doi: 10.1038/s41586-020-2129-8

PubMed Abstract | CrossRef Full Text | Google Scholar

97. Larabi A, Barnich N, Nguyen HT. New insights into the interplay between autophagy, gut microbiota and inflammatory responses in IBD. Autophagy. (2020) 16:38–51. doi: 10.1080/15548627.2019.1635384

PubMed Abstract | CrossRef Full Text | Google Scholar

98. Lin X, Peng C, Greenbaum J, Li ZF, Wu KH, Ao ZX, et al. Identifying potentially common genes between dyslipidemia and osteoporosis using novel analytical approaches. Mol Genet Genomics. (2018) 293:711–23. doi: 10.1007/s00438-017-1414-1

PubMed Abstract | CrossRef Full Text | Google Scholar

99. Jieru S, San W, Zheng Y, Zhang S, Cao D, Chen Y, et al. High glucose promotes proliferation and invasiveness in pancreatic cancer cells by upregulating RET, a proto-oncogene encoding a receptor tyrosine kinase. doi: 10.1016/j.biopha.2023.115802

CrossRef Full Text | Google Scholar

100. Allison SJ. A lncRNA and miRNA megacluster in diabetic nephropathy. Nat Rev Nephrol. (2016) 12:713. doi: 10.1038/nrneph.2016.151

CrossRef Full Text | Google Scholar

101. Özcan S. Minireview: microRNA function in pancreatic β cells. Mol endocrinology. (2014) 28:1922–33. doi: 10.1101/2022.05.30.493783

CrossRef Full Text | Google Scholar

102. Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. (2012) 22:1775–89. doi: 10.1101/gr.132159.111

CrossRef Full Text | Google Scholar

103. Guo J, Liu Z, Gong R. Long noncoding RNA: an emerging player in diabetes and diabetic kidney disease. Clin Science. (2019) 133:1321–39. doi: 10.1042/CS20190372

CrossRef Full Text | Google Scholar

104. Laurent GS, Wahlestedt C, Kapranov P. The Landscape of long noncoding RNA classification. Trends Genet. (2015) 31:239–51. doi: 10.1016/j.tig.2015.03.007

CrossRef Full Text | Google Scholar

105. Ali HS, Kamel MM, Agwa SH, Hakeem MS, Meteini MS, Matboli M. Analysis of mRNA-miRNA-lncRNA differential expression in prediabetes/type 2 diabetes mellitus patients as potential players in insulin resistance. Front Endocrinology. (2023) 14:1131171. doi: 10.3389/fendo.2023.1131171

CrossRef Full Text | Google Scholar

106. Chen G, Hu J, Huang Z, Yang L, Chen M. MicroRNA-1976 functions as a tumor suppressor and serves as a prognostic indicator in non-small cell lung cancer by directly targeting PLCE1. Biochem Biophys Res Commun. (2016) 473:1144–51. doi: 10.1016/j.bbrc.2016.04.030

CrossRef Full Text | Google Scholar

107. Wang J, Li M, Han X, Wang H, Wang X, Ma G, et al. MiR-1976 knockdown promotes epithelial–mesenchymal transition and cancer stem cell properties inducing triple-negative breast cancer metastasis. Cell Death Disease. (2020) 11:500. doi: 10.1038/s41419-020-2711-x

CrossRef Full Text | Google Scholar

108. Renna M, Bento CF, Fleming A, Menzies FM, Siddiqi FH, Ravikumar B, et al. IGF-1 receptor antagonism inhibits autophagy. Hum Mol Genet. (2013) 22:4528–44. doi: 10.1093/hmg/ddt300

CrossRef Full Text | Google Scholar

109. Bao T, Liu J, Leng J, Cai L. The cGAS–STING pathway: more than fighting against viruses and cancer. Cell Bioscience. (2021) 11:1–20. doi: 10.1186/s13578-021-00724-z

CrossRef Full Text | Google Scholar

110. Yamamoto A, Cremona ML, Rothman JE. Autophagy-mediated clearance of huntingtin aggregates triggered by the insulin-signaling pathway. J Cell Biol. (2006) 172:719–31. doi: 10.1083/jcb.200510065

CrossRef Full Text | Google Scholar

111. Feng Y, Xu W, Zhang W, Wang W, Liu T, Zhou X. LncRNA DCRF regulates cardiomyocyte autophagy by targeting miR-551b-5p in diabetic cardiomyopathy. Theranostics. (2019) 9:4558. doi: 10.7150/thno.31052

CrossRef Full Text | Google Scholar

112. Liu Y, Levine B. Autosis and autophagic cell death: the dark side of autophagy. Cell Death Differentiation. (2015) 22:367–76. doi: 10.1038/cdd.2014.143

CrossRef Full Text | Google Scholar

113. Taheri M, Eghtedarian R, Ghafouri-Fard S, Omrani MD. Non-coding RNAs and type 2 diabetes mellitus. Arch Physiol Biochem. (2023) 129:526–35. doi: 10.1080/13813455.2020.1843498

CrossRef Full Text | Google Scholar

114. Kokkinopoulou I, Maratou E, Mitrou P, Boutati E, Sideris DC, Fragoulis EG, et al. Decreased expression of microRNAs targeting type-2 diabetes susceptibility genes in peripheral blood of patients and predisposed individuals. Endocrine. (2019) 66:226–39. doi: 10.1007/s12020-019-02062-0

CrossRef Full Text | Google Scholar

115. Samidurai A, Ockaili R, Cain C, Roh SK, Filippone SM, Kraskauskas D, et al. Differential regulation of mTOR complexes with miR-302a attenuates myocardial reperfusion injury in diabetes. Iscience. (2020) 23. doi: 10.1016/j.isci.2020.101863

CrossRef Full Text | Google Scholar

116. Samidurai A, Salloum FN, Durrant D, Chernova OB, Kukreja RC, Das A. Chronic treatment with novel nanoformulated micelles of rapamycin, Rapatar, protects diabetic heart against ischaemia/reperfusion injury. Br J Pharmacol. (2017) 174:4771–84. doi: 10.1111/bph.14059

CrossRef Full Text | Google Scholar

117. Yuan Q, Chen Z, Santulli G, Gu L, Yang ZG, Yuan ZQ, et al. Functional role of Calstabin2 in age-related cardiac alterations. Sci Rep. (2014) 4:7425. doi: 10.1038/srep07425

CrossRef Full Text | Google Scholar

118. Gao G, Chen W, Yan M, Liu J, Luo H, Wang C, et al. Rapamycin regulates the balance between cardiomyocyte apoptosis and autophagy in chronic heart failure by inhibiting mTOR signaling. Int J Mol Med. (2020) 45:195–209. doi: 10.1016/j.biopha.2020.111176

CrossRef Full Text | Google Scholar

119. Jia M, Wei Z, Liu P, Zhao X. Silencing of ABCG2 by microRNA-3163 inhibits multidrug resistance in retinoblastoma cancer stem cells. J Korean Med science. (2016) 31:836–42. doi: 10.3346/jkms.2016.31.6.836

CrossRef Full Text | Google Scholar

120. Yang B, Wang C, Xie H, Wang Y, Huang J, Rong Y, et al. MicroRNA-3163 targets ADAM-17 and enhances the sensitivity of hepatocellular carcinoma cells to molecular targeted agents. Cell Death disease. (2019) 10:784. doi: 10.1038/s41419-019-2023-1

CrossRef Full Text | Google Scholar

121. Defronzo RA. Banting Lecture. From the triumvirate to the ominous octet: a new paradigm for the treatment of type 2 diabetes mellitus. . Diabetes. (2009) 58:773–95. doi: 10.2337/db09-9028

CrossRef Full Text | Google Scholar

Keywords: type 2 diabetes, therapeutic targets, machine learning, drug response, rats

Citation: Matboli M, Al-Amodi HS, Khaled A, Khaled R, Roushdy MMS, Ali M, Diab GI, Elnagar MF, Elmansy RA, TAhmed HH, Ahmed EME, Elzoghby DMA, M.Kamel HF, Farag MF, ELsawi HA, Farid LM, Abouelkhair MB, Habib EK, Fikry H, Saleh LA and Aboughaleb IH (2024) Comprehensive machine learning models for predicting therapeutic targets in type 2 diabetes utilizing molecular and biochemical features in rats. Front. Endocrinol. 15:1384984. doi: 10.3389/fendo.2024.1384984

Received: 11 February 2024; Accepted: 03 May 2024;
Published: 24 May 2024.

Edited by:

Yun Shen, Pennington Biomedical Research Center, United States

Reviewed by:

Yumin Ma, Yangzhou University, China
Dan Xue, University of Pittsburgh, United States
Zhe Sang, Icahn School of Medicine at Mount Sinai, United States
Song Yi, First Affiliated Hospital of Zhengzhou University, China

Copyright © 2024 Matboli, Al-Amodi, Khaled, Khaled, Roushdy, Ali, Diab, Elnagar, Elmansy, TAhmed, Ahmed, Elzoghby, M.Kamel, Farag, ELsawi, Farid, Abouelkhair, Habib, Fikry, Saleh and Aboughaleb. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marwa Matboli, RHJNYXJ3YV9NYXRib2x5QG1lZC5hc3UuZWR1LmVn; Ibrahim H. Aboughaleb, ZWhlNDNAaG90bWFpbC5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.