- Research
- Open access
- Published:
miRDM-rfGA: Genetic algorithm-based identification of a miRNA set for detecting type 2 diabetes
BMC Medical Genomics volume 16, Article number: 195 (2023)
Abstract
Background
Type 2 diabetes mellitus (T2DM) affects approximately 451 million adults globally. In this study, we identified the optimal combination of marker candidates for detecting T2DM using miRNA-Seq data from 95 samples including T2DM and healthy individuals.
Methods
We utilized the genetic algorithm (GA) in the discovery of an optimal miRNA biomarker set. We discovered miRNA subsets consisting of three miRNAs for detecting T2DM by random forest-based GA (miRDM-rfGA) as a feature selection algorithm and created six GA parameter settings and three settings using traditional feature selection methods (F-test and Lasso). We then evaluated the prediction performance to detect T2DM in the miRNA subsets derived from each setting.
Results
The miRNA subset in setting 5 using miRDM-rfGA performed the best in detecting T2DM (mean AUROC = 0.92). Target mRNA identification and functional enrichment analysis of the best miRNA subset (hsa-miR-125b-5p, hsa-miR-7-5p, and hsa-let-7b-5p) validated that this combination was involved in T2DM. We also confirmed that the targeted genes were negatively correlated with the clinical variables related to T2DM in the BxD mouse genetic reference population database.
Conclusions
Using GA in miRNA-Seq data, we identified the optimal miRNA biomarker set for T2DM detection. GA can be a useful tool for biomarker discovery and drug-target identification.
Background
Type 2 diabetes mellitus (T2DM) is characterized by abnormalities in carbohydrate, lipid, and protein metabolism pathways. Dysregulation of insulin secretion and response in T2DM results in hyperglycemia.
According to the International Diabetes Federation, 451 million adults worldwide have diabetes; this number is expected to reach 693 million by 2045 [1]. Diabetes is among the top 10 causes of death globally, and the risk of all-cause mortality increases by approximately 2–threefold among individuals with diabetes [2].
T2DM is a representative disease that leads to diabetic nephropathy, retinopathy, neuropathy, and other complications, including colon and liver cancers [3,4,5,6,7,8,9].
The main objective of this study was to identify the optimal biomarkers for detecting T2DM and identifying drug targets to treat T2DM. Among the various putative drug targets, epigenetic mechanisms such as microRNAs (miRNAs) may contribute to the development of common diseases, including T2DM [10, 11]. miRNAs, which are short (~ 22 nucleotides) noncoding RNAs, have emerged as key cell type-specific regulators of gene expression, operating primarily to inhibit target genes post-transcription by binding with complementary mRNA [12].
In T2DM, miRNAs target various genes related to glucose and fatty acid metabolism and the insulin signaling pathway in diverse tissues (e.g., skeletal muscle, pancreas, adipocytes, and liver), thereby affecting physiological functions [10, 13, 14].
Because of the importance of the regulation of miRNAs in T2DM, several studies have tried to identify miRNA biomarkers for T2DM by characterizing differentially expressed miRNAs (in blood, pancreas, adipocytes, skeletal muscle, and liver) in T2DM patients [3,4,5,6,7,8,9, 15,16,17,18,19]. In addition, miRNA biomarker discovery has confirmed the negative correlation between the expression of discovered miRNAs and their target mRNAs [3,4,5,6,7,8,9,10,11, 14,15,16,17,18,19].
Biomarker discovery based on feature selection and applying the machine learning method is a promising approach [20,21,22,23,24,25]. However, these methods use genetic data to create a model that classifies T2DM patients and to derive important features that affect the classification performance of the model. In this case, it is possible to calculate the importance of miRNAs that affect disease diagnosis. However, the optimal number of entries in the biomarker combination for disease diagnosis is difficult to define using these methods. In other words, it is difficult to determine an optimal combination of marker candidates. As marker candidates belong to differentially expressed (DE) genes (or miRNAs), and DE genes are usually numerous, identification of the best combination of select markers is challenging [26]. Even though the combination of machine learning and genetic algorithm (GA) techniques can be adopted for addressing these challenges, it still has not been widely explored in the context of T2DM biomarker discovery.
To identify an optimal miRNA set for T2DM, we integrated GA and random forest to develop a novel feature selection algorithm. We first obtained public miRNA-Seq datasets from T2DM patients and healthy controls (HC) from the Gene Expression Omnibus (GEO). We then compared diverse feature selection methods (i.e., F-test in analysis of variance [ANOVA] and least absolute shrinkage and selection operator [Lasso]) under three settings with our GA-based feature selection (named miRDM-rfGA) under six settings. In each setting, we evaluated the performance of each biomarker set in detecting T2DM. We then obtained publicly available mouse phenotype data to show biological associations between T2DM-related phenotypes and the biomarker set [27]. In this study, we demonstrated the utility of GA for the discovery of an optimal miRNA biomarker set. This study not only emphasizes the significance of miRNAs in T2DM but also provides the novel application of GA and machine learning techniques in the discovery of optimal combinations of disease biomarkers.
Materials and methods
Data collection
We searched for and acquired public miRNA-Seq data in the blood tissue for biomarker discovery (Fig. 1). First, we accessed Sequence Read Archive (SRA) [28] and obtained relevant datasets by using the search terms “Diabetes” and “miRNA”. After acquiring the results from the query, two datasets (SRP151126 and SRP093728) of miRNA profiling related to T2DM in human blood were obtained (Supplementary Figure S1 and Supplementary Table S1). Thus, we selected 95 samples (56 HCs and 39 T2DM patients) to download fastq sequences and perform further analyses (Supplementary Table S1).
Pre-processing of miRNA-Seq data and dataset preparation
For the 95 miRNA-Seq samples, we pre-processed the data using FastQC v0.11.7 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and Cutadapt v1.3 [29]. Samples with a Phred quality score of less than 20 were removed. We used Illumina Universal Adapter and Illumina Small RNA 3’ Adapter for adapter trimming. Subsequently, we trimmed the reads with lengths of 18–25 bp. The miRdeep2 mapper (v 2.0.0.8) [30] was used for sequence alignment of the human reference genome version hg19, and miRNAs were annotated using miRbase version 21 [31]. For quantifying the miRNA expression, reads per million (RPM) were used and Z-normalization was used to display the miRNA expression.
To remove sparse miRNAs, miRNAs with a ratio of 90% or more in a sample with an expression value of 0 for each miRNA from a total of 95 samples were removed, yielding a total of 1169 miRNAs.
The training and test sets were created in a ratio of 7:3. The training set was used for biomarker discovery using the GA or traditional feature selection method, and the test set was used to evaluate the performance of the T2DM patient classification model using biomarkers selected by each method.
Biomarker set discovery using GA and traditional feature selection methods
According to previous studies [32,33,34,35,36], the optimal number of entries in biomarker sets range from two to nine. Inspired by a previous study using three miRNAs as biomarker combinations [36], we aimed to discover an optimal biomarker set consisting of three miRNAs for detecting T2DM. Therefore, we used the feature selection method to identify a biomarker set comprising three miRNAs. To this end, GA was adopted, and we used the F-test in ANOVA and Lasso as a traditional feature selection method for comparison with GA.
GA is a heuristic algorithm that uses natural selection to find the best solution [37,38,39]. To simulate natural selection, the GA generates a population and a set of individuals (i.e., chromosomes and solutions). Individuals consist of a set of features (or genes) that represent a set of solutions to a predefined problem.
For each individual, the fitness score was calculated to show how each solution, indicating selected features, was optimal. The best individual (i.e., the best solution) in the population was selected and included in a new population for the next generation, and the other individuals were produced by crossover and mutation. Thus, the GA generates individuals repeatedly, assesses their fitness, and terminates when the given goal is met or when some stopping criteria are met.
In this study, we implemented a random forest (RF)-based GA feature selection algorithm (miRDM-rfGA) with miRNA-Seq data as input, and determined the three most optimal miRNA subsets evaluated by a fitness score using RF.
For miRDM-rfGA, five phases were considered: (1) generation of the initial population, (2) evaluation of each individual by a fitness function that explains whether the solution is good or not, (3) selection of individuals with the highest fitness score, indicating the best RF performance (e.g., area under the receiver operating characteristic [AUROC]), (4) crossover, and (5) mutation for new individuals in the next generation (Supplementary Figure S2).
In the first step of miRDM-rfGA, we generated the initial population to undergo successive evolution through the GA. A large number of individuals with diverse solutions within the initial population is necessary for successfully deriving optimal solutions. Therefore, we randomly generated 1,000 individuals, limited by our computational power available. The randomly generated 1000 individuals have diverse solutions for accessing the miRNA expression dataset. Within the initial population, each individual carried a specific number of miRNAs (described as “genes”) as features (Supplementary Figure S2A). The features of each individual were mapped to 1139 miRNA indices in the pre-processed dataset. In other words, an individual (ind) consists of 1139 features. the features were represented as ind = (x1, x2, …, x1139). xi has a binary variable (i.e., {0,1}) where 1 indicates selection of the i-th miRNA and 0 indicates non-selection of the miRNA in the dataset. Subsequently, miRDM-rfGA randomly selects features from the 1139 indices according to the predefined number of selected features for each individual in the initial population (Supplementary Figure S2A). Then, given the individual, the expression profiles of the selected features (miRNAs) were extracted from the miRNA-Seq input data (i.e., the expression matrix of samples by miRNAs), generating a subset of the input data as a training set for the RF model. With the selected features in each individual, we calculated a fitness score as follows (fitness function) (Supplementary Figure S2B):
where AUCk is the AUROC from the RF model for classifying T2DM and HC in the k-th fold during M fold cross-validation in the training set; \(x\) is the number of selected features (miRNAs) in the individual; W is the penalty weight; and b is the optimal number of selected genes for the optimal biomarker combination.
Each individual was evaluated with a fitness score, and the individual with the highest fitness score was included in the population. The individual was then included in the next generation, and the other individuals in the next generation were produced by crossover and mutation (Supplementary Figure S2C). miRDM-rfGA iterates these phases for G generations and derives the best individual among these generations (Supplementary Figure S2D). Based on miRDM-rfGA, we created six parameter settings (settings 1 to 6), and the parameters N, W, b, G, crossover rate, and mutation rate are listed in Table 1. For the GA process, “DEAP” v.1.3.1 in Python 3.7 were used [40].
For traditional feature selection using the F-test in ANOVA, scikit-learn's f_classif function was used. Feature importance was calculated for each feature using f_classif, the top three miRNAs with the highest feature importance were selected through SelectKBest (k = 3), and the RF method was used to evaluate the discrimination power of T2DM for the selected miRNA biomarker set (setting 7).
For Lasso, we applied logistic regression using L1-regularization and used the SelectedFromModel (k = 3) function. Then, Lasso (setting 8) and RF (setting 9) were used as models for classifying T2DM using the selected miRNA biomarker set. Detailed information on settings 7, 8, and 9 is presented in Table 2.
Performance comparison
For each trained model in each setting, the mean AUROC score was calculated using threefold cross-validation in the test set. The performance of each setting was compared using the mean AUROC.
Principal component analysis
For the three best miRNA biomarker sets, we devised a method to distinguish between T2DM and HC using the most optimum miRNA biomarker set. Therefore, principal component analysis (PCA) was performed with three components using the miRNA-Seq data.
Target mRNA identification and pathway enrichment analysis
MIENTURNET (http://userver.bio.uniroma1.it/apps/mienturnet/, accessed on 12 April, 2021) [41] was used for identifying the putative target genes of the differentially expressed circulating miRNAs and analyzing their pathway enrichment.
Correlation analysis in mouse population data
The purpose of this analysis was to determine the correlation between diabetes-related clinical indicators and gene expression targeted by selected miRNAs using publicly available mouse population data [27]. For this analysis, two gene expression datasets in mouse liver tissue were obtained from the BxD mouse high-fat diet (HFD) [EPFL/LISP BXD HFD Liver Affy Mouse Gene 1.0 ST (Aug 18) RMA] and chow diet (CD) [EPFL/LISP BXD CD Liver Affy Mouse Gene 1.0 ST (Aug 18)] cohort from the GeneNetwork 1 database (http://gn1.genenetwork.org) [27].
The mRNA expression levels of IGF1R, IRS2, and PIK3CD were confirmed from mRNA gene expression array data in the liver tissue of the HFD and CD cohorts. In the gene expression data, we compared the top 25% of each gene with the bottom 25% of the group and clinical parameters.
Spearman’s correlation analysis was performed between gene expression data and the clinical parameters insulin response (IR) and oral glucose tolerance test (OGTT).
Statistical analysis
Student’s t-test was used for analysis of differentially expressed miRNAs and DE genes in BxD mouse data between groups (T2DM vs. HC). For the DE analysis in the GEO gene and miRNA expression array dataset, we used GEO2R [42]. Correlations were evaluated using Spearman’s correlation coefficient. All reported P values were statistically significant when less than 0.05.
Results
Performance of feature selection using miRDM-rfGA and traditional feature selection
In this study, we identified the optimal miRNA biomarker combination. For this, we used mirDM-rfGA and traditional feature selection methods and evaluated the prediction performance of the miRNA subsets derived from each method. The workflow for determining the optimal miRNA features is presented in Fig. 1.
First, we derived the three miRNA subsets using mirDM-rfGA under six settings (settings 1 to 6) (Fig. 2A, B, C, D, E, F and Table 1). Using each miRNA subset data derived from each setting, we constructed an RF model classifying T2DM patients. Each RF model was compared with the mean AUROC (± 1 standard deviation) through threefold cross-validation using the test set.
Among the six settings using mirDM-rfGA, setting 5 showed the best performance. The AUROC of the model was calculated through the test set, and the mean AUROC was 0.92 ± 0.04. The miRNAs subset using setting 5 included 'hsa-let-7b-5p,' 'hsa-miR-125b-5p,' and 'hsa-miR-7-5p,' (Fig. 2E and Table 1).
We also applied traditional feature selection using univariate feature selection methods (F-test in ANOVA and Lasso; settings 7 to 9; Fig. 2G, H, I and Table 2) for comparison with the settings using mirDM-rfGA.
As a result, the mean AUROC value of setting 7 was 0.72 ± 0.08 ('hsa-miR-6820-5p,' 'hsa-miR-29b-2-5p,' and 'hsa-miR-1307-3p'); that of setting 8 was 0.64 ± 0.05 ('hsa-miR-22-3p,' 'hsa-miR-92a-3p,' and 'hsa-miR-181a-5p'), and that of setting 9 was 0.52 ± 0.02 ('hsa-miR-22-3p,' 'hsa-miR-92a-3p,' and 'hsa-miR-181a-5p') (Table 2).
In summary, the setting 5 using mirDM-rfGA showed the best performance in detecting T2DM, and the miRNA set derived by the setting 5 included 'hsa-let-7b-5p,' 'hsa-miR-125b-5p,' and 'hsa-miR-7-5p'.
The log2 fold change of T2DM vs. HC for each miRNA in setting 5 was 0.468, − 0.853, and 0.953, respectively (Fig. 3A and Table 3). The P values were also statistically significant at 4.33e − 5, 4.59e − 4, and 1.75e − 4, respectively. The DE analysis results and statistical significance of miRNAs derived from other settings are described in Fig. 3B–H and Table 3.
The relationship between the best miRNA biomarker set and their target mRNAs
We used PCA to determine how the optimal miRNA biomarker subset (hsa-miR-125b-5p, hsa-miR-7-5p, and hsa-let-7b-5p) distinguished between T2DM patients and HCs. The T2DM and HC groups were differentiated using only the expression levels of the three miRNAs (Fig. 4A). After confirming that these three miRNAs can discriminate against the diabetic group, functional enrichment analysis of the miRNA-target gene was conducted based on the miRTarbase database [43, 44] to determine the signaling pathway to which the target genes of the three miRNAs belong. All three miRNAs were found to be most enriched in insulin receptor substrate (IRS)-related events triggered by IGF1R in the Reactome database (Fig. 4B) [45].
In this pathway, the target genes were PIK3CB, IGF1R targeted by hsa-miR-125b-5p, IRS2, IGF1R, and PIK3CD targeted by hsa-miR-7-5p, and IRS2 and IGF1R targeted by hsa-let-7b-5p.
Interestingly, all three miRNAs target IGF1R in the pathway network (Fig. 4B).
Next, we confirmed the expression patterns of the genes targeted by the miRNA biomarker in a public dataset of T2DM studies.
In the skeletal muscle tissue (GSE22309), we also confirmed that the gene expression levels of IRS2, IGF1R, and PIK3CD were decreased in patients with T2DM after insulin treatment compared with those in healthy participants after insulin treatment (Fig. 4C and Supplementary Table S2). However, the regulation pattern of hsa-miR-125b-5p on its targets (IGF1R and PIK3CB) was not observed (Fig. 4C).
Correlation with diabetes-related clinical variables according to the gene expression targeted by the miRNA biomarkers
We also investigated the correlation between mRNA targeted by miRNA biomarkers in IRS-related events by the IGF1R pathway and clinical indicators related to diabetes in a publicly available BxD mouse database [27].
Correlation analysis was performed using the GeneNetwork BxD mouse database. Among various clinical variables, when the expression levels of IGF1R, PIK3CD, and IRS2 were low in HFD and CD mice, both glucose levels in the OGTT and IR values during OGTT increased, and these results showed a negative correlation (Fig. 5 and Supplementary Table S3).
Discussion
In this study, we identified a miRNA biomarker set consisting of three miRNAs for detecting T2DM using public miRNA-Seq data. To this end, we devised miRDM-rfGA, which is an GA-based feature selection algorithm, and created biomarker discovery settings using GA (six settings) and traditional feature selection methods (3 settings using F-test and Lasso) by constructing RF classifier models to detect T2DM, using miRNAs obtained through each setting. We then compared the classification performance of the biomarker sets from all settings. As a result, setting 5 using miRDM-rfGA (mean AUROC: 0.92) outperformed traditional feature selection approaches in identifying these biomarkers. From the setting 5, we determined that a set of biomarkers consisting of three miRNAs (hsa-miR-125b-5p, hsa-miR-7-5p, and hsa-let-7b-5p) had the most optimal classification performance for T2DM (Fig. 2E). To confirm the association between the discovered miRNA biomarker set and diabetes, we performed a functional enrichment analysis of miRNA-targeted mRNAs. As a result, IRS-related signaling by IGF1R, which is directly related to diabetes, was derived (IRS2, IGF1R, PIK3CD, and PIK3CB) (Fig. 4C). In addition, correlation analysis was performed between the expression levels of the corresponding mRNA and clinical variables related to the detection of diabetes using the public data of BxD mice (Fig. 5) [27]. These results confirmed that the gene expression levels in diabetes were consistent with those from our study.
In addition, we confirmed that there are various studies supporting that three miRNAs (hsa-miR-125b-5p, hsa-miR-7-5p, and hsa-let-7b-5p) are related to diabetes.
Let-7b-5p is a miRNA belonging to the let-7 family, which has a common seed region up to nucleotides 2–8, and is known to perform similar functions [46]. In this study, the relationship between let-7 and diabetes showed that glucose tolerance was inhibited in let-7 overexpressing mice, and let-7 knockdown via let-7 anti-miR restored glucose intolerance caused by obesity [46]. In addition, functional recovery of insulin signaling was observed in the muscle and liver following let-7 anti-miR treatment [46]. In another study, a Lin28 transgenic mouse model showed that let-7 downregulation effectively reduced glucose levels in a high-fat diet and activated the Insulin-PI3K-mTOR signaling pathway through let-7 suppression [47]. In our study (Fig. 3E), let-7b-5p was overexpressed in the group of patients with T2DM as compared to that in the HCs. The results of this study confirmed that the overexpression of let-7 was consistent with the study results of miRNAs adversely affecting glucose tolerance and insulin sensitivity.
miR-7 is known as a miRNA related to insulin secretion and glucose homeostasis in the insulin signaling pathway. Agbu et al. reported that circulating glucose and triglyceride levels were increased in Drosophila insulin-producing cells (IPCs) overexpressing miR-7 as compared to those in wild-type cells [48]. Additionally, miR-7-5p inhibits glucose uptake and insulin-induced AKT phosphorylation by affecting the downregulation of IRS-1 and IRS-2 [16, 49]. Expression patterns of miR-7-5p and IRS2 were also observed in the present study (Figs. 3E and 4C).
According to Gong et al., the expression of hsa-miR-125b-5p is downregulated in the retinas of rats with STZ-induced diabetes [50]. They also confirmed that miR-125b-5p was downregulated with the progression of diabetic retinopathy (DR), suggesting that hsa-miR-125b-5p might be used as an effective T2DM and DR treatment, and the expression pattern of hsa-miR-125b-5p was identical to that observed in our study (Fig. 3E).
The genes targeted by the three miRNAs for the optimal miRNA biomarker set of T2DM were IRS2, IGF1R, and PIK3CD. IRS2 is an insulin receptor substrate 2 involved in insulin sensitivity as an insulin signaling pathway, and whose expression levels are decreased in T2DM [51]. In a mouse model study, when IRS2 was knocked out, obesity and insulin sensitivity decreased. As a result, glucose tolerance was induced and developed into T2DM. This led to hyperinsulinemia and β-cell damage [52, 53].
IGF1R functions as an insulin receptor by forming a dimer with the insulin receptor as an insulin-like growth factor-1 receptor [54]. Dong et al. confirmed that the expression level of IGF1R was decreased in the liver of diabetic rats, and the regulation of this gene by miRNAs could play a role in the improvement of insulin resistance [55]. Razny et al. compared the expression levels of whole blood mRNA between an obese group with insulin resistance and an obese group without insulin resistance [56], and found that the gene expression level of IGF1R had decreased. Therefore, the decreasing pattern in the gene expression level of IGF1R coincided with the increasing pattern in the expression level of hsa-let-7b-5p and hsa-miR-7-5p targeting IGF1R.
PIK3CD refers to Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit delta, which is a catalytic subunit of PI3K. PIK3CD participates in PI3-Kinase signaling and affects the AKT pathway, and the gene expression level of PIK3CD is diminished in the skeletal muscle tissue of diabetic patients [57]. In addition, inhibition of PI3K signaling in skeletal muscle tissue in mouse models results in insulin resistance and systemic glucose intolerance. Further, free fatty acid and triglyceride levels in the blood are elevated [58]. Insulin resistance occurs when PI3K signaling is inhibited [59, 60] and has been observed in adipocytes [61, 62], muscle cells [57], the liver [63], and blood [56]. Du et al. reported that PIK3CD-targeting miRNA can inhibit the insulin signaling pathway and thus become a target gene that can regulate insulin resistance [63].
By devising miRDM-rfGA, we identified a set of putative diagnostic and treatment biomarkers for T2DM using GA. Our study was limited by the number of samples (95). However, to compensate for this limitation, we used correlation analysis with the public BxD mouse database to confirm that the mRNAs targeted by the miRNA biomarker set were also related to T2DM [27]. In addition, our study was limited to discovering diagnostic and treatment biomarkers composed of three miRNAs. Based on this, an extended follow-up study may be conducted.
In addition, while our method could have allowed for an extended study into diabetic complications along with diabetes, our study was also focused on only type 2 diabetes. However, using a publicly available miRNA dataset (GEO accession: GSE51674) covering ‘diabetic nephropathy with type 2 diabetes’ (T2DN) in kidney tissue (Supplementary Table S1), hsa-let-7b, hsa-miR-125b were up-regulated in T2DN (Supplementary Table S4 and Supplementary Figure S3). Among these miRNAs, we may indicate hsa-let-7b as a biomarker candidate not only for T2DM but also for T2DN.
This study suggested that the feature selection method combining genetic algorithms (GA) and machine learning is superior to traditional feature selection methods, and also that this approach can be useful in deriving optimal biomarker combinations that can be applied across various diseases.
Conclusions
We derived an optimal miRNA biomarker set that could detect T2DM using GA to process miRNA-Seq data. Thus, GA can be used as an effective method for biomarker discovery.
Availability of data and materials
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI Sequence Read Archive (SRA), accession no: PRJNA476995 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA476995) and PRJNA354381 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA354381); NCBI Gene Expression Omnibus (GEO), accession no: GSE22309 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE22309) and GSE51674 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE51674). The source code for the miRDM-rfGA in setting 5 is available at GitHub (https://github.com/labnams/miRDM-rfGA).
Abbreviations
- T2DM :
-
Type 2 diabetes mellitus
- HC :
-
Healthy control
- miRNA :
-
MicroRNA
- SRA :
-
NCBI Sequence Read Archive
- GEO :
-
NCBI Gene Expression Omnibus
- R PM :
-
Reads per million
- DE :
-
Differentially expressed
- GA :
-
Genetic algorithm
- AUROC :
-
Area under receiver operating characteristic
- P CA :
-
Principal component analysis
- R F :
-
Random forest
- HFD :
-
High fat diet
- CD :
-
Chow diet
- IR :
-
Insulin response
- OGTT :
-
Oral glucose tolerance test
- IPC :
-
Drosophila insulin production cell
- DR :
-
Diabetic retinopathy
- T2DN :
-
Diabetic nephropathy with type 2 diabetes
References
Lin X, Xu Y, Pan X, Xu J, Ding Y, Sun X, et al. Global, regional, and national burden and trend of diabetes in 195 countries and territories: an analysis from 1990 to 2025. Sci Rep. 2020;10(1):14790.
Yang JJ, Yu D, Wen W, Saito E, Rahman S, Shu X-O, et al. Association of Diabetes With All-Cause and Cause-Specific Mortality in Asia: A Pooled Analysis of More Than 1 Million Participants. JAMA Netw Open. 2019;2(4):e192696.
Chen M, Sun Q, Giovannucci E, Mozaffarian D, Manson JE, Willett WC, et al. Dairy consumption and risk of type 2 diabetes: 3 cohorts of US adults and an updated meta-analysis. BMC Med. 2014;12:215-.
Chien H-Y, Lee T-P, Chen C-Y, Chiu Y-H, Lin Y-C, Lee L-S, et al. Circulating microRNA as a diagnostic marker in populations with type 2 diabetes mellitus and diabetic complications. J Chin Med Assoc. 2015;78(4):204–11.
de Candia P, Spinetti G, Specchia C, Sangalli E, La Sala L, Uccellatore A, et al. A unique plasma microRNA profile defines type 2 diabetes progression. PLoS One. 2017;12(12):e0188980.
Miao C, Zhang G, Xie Z, Chang J. MicroRNAs in the pathogenesis of type 2 diabetes: new research progress and future direction. Can J Physiol Pharmacol. 2017;96(2):103–12.
Vasu S, Kumano K, Darden CM, Rahman I, Lawrence MC, Naziruddin B. MicroRNA Signatures as Future Biomarkers for Diagnosis of Diabetes States. Cells. 2019;8(12):1533.
Yang J-S, Lu C-C, Kuo S-C, Hsu Y-M, Tsai S-C, Chen S-Y, et al. Autophagy and its link to type II diabetes mellitus. Biomedicine (Taipei). 2017;7(2):8.
Yaribeygi H, Katsiki N, Behnam B, Iranpanah H, Sahebkar A. MicroRNAs and type 2 diabetes mellitus: Molecular mechanisms and the effect of antidiabetic drug treatment. Metabolism. 2018;87:48–55.
Eliasson L, Esguerra JLS. MicroRNA Networks in Pancreatic Islet Cells: Normal Function and Type 2 Diabetes. Diabetes. 2020;69(5):804–12.
Rosado AJ, Diez-Bello R, Salido MG, Jardin I. Fine-tuning of microRNAs in Type 2 Diabetes Mellitus. Curr Med Chem. 2019;26(22):4102–18.
Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. 2007;8(2):93–103.
Calderari S, Diawara MR, Garaud A, Gauguier D. Biological roles of microRNAs in the control of insulin secretion and action. Physiol Genomics. 2016;49(1):1–10.
Kostyniuk DJ, Marandel L, Jubouri M, Dias K, de Souza RF, Zhang D, et al. Profiling the rainbow trout hepatic miRNAome under diet-induced hyperglycemia. Physiol Genomics. 2019;51(9):411–31.
Deng J, Guo F. MicroRNAs and type 2 diabetes. ExRNA. 2019;1(1):36.
Feng J, Xing W, Xie L. Regulatory Roles of MicroRNAs in Diabetes. Int J Mol Sci. 2016;17(10):1729.
Jiménez-Lucena R, Rangel-Zúñiga OA, Alcalá-Díaz JF, López-Moreno J, Roncero-Ramos I, Molina-Abril H, et al. Circulating miRNAs as predictive biomarkers of type 2 diabetes mellitus development in coronary heart disease patients from the CORDIOPREV study. Mol Ther Nucleic Acids. 2018;12:146–57.
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(2):226–39.
Pordzik J, Jakubik D, Jarosz-Popek J, Wicik Z, Eyileten C, De Rosa S, et al. Significance of circulating microRNAs in diabetes mellitus type 2 and platelet reactivity: bioinformatic analysis and review. Cardiovasc Diabetol. 2019;18(1):113.
Shahrjooihaghighi A, Frigui H, Zhang X, Wei X, Shi B, Trabelsi A. An ensemble feature selection method for biomarker discovery. Proc IEEE Int Symp Signal Proc Inf Tech. 2017;2017:416–21.
He Z, Yu W. Stable feature selection for biomarker discovery. Comput Biol Chem. 2010;34(4):215–25.
Shi Z, Wen B, Gao Q, Zhang B. Feature selection methods for protein biomarker discovery from proteomics or multiomics data. Mol Cell Proteomics. 2021;20:100083.
Saeys Y, Inza I, Larrañaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23(19):2507–17.
Abeel T, Helleputte T, Van de Peer Y, Dupont P, Saeys Y. Robust biomarker identification for cancer diagnosis with ensemble feature selection methods. Bioinformatics. 2010;26(3):392–8.
Dessì N, Pascariello E, Pes B. A comparative analysis of biomarker selection techniques. Biomed Res Int. 2013;2013:387673.
Vandewater L, Brusic V, Wilson W, Macaulay L, Zhang P. An adaptive genetic algorithm for selection of blood-based biomarkers for prediction of Alzheimer’s disease progression. BMC Bioinformatics. 2015;16(18):S1.
Wu C-C, Huang H-C, Juan H-F, Chen S-T. GeneNetwork: an interactive tool for reconstruction of genetic networks using microarray data. Bioinformatics. 2004;20(18):3691–3.
Leinonen R, Sugawara H, Shumway M. International nucleotide sequence database C. The sequence read archive. Nucleic Acids Res. 2011;39(Database issue):D19–21.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal: Next Generat Sequencing Data Analys. 2011;17:1.
Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.
Griffiths-Jones S. miRBase: The MicroRNA Sequence Database. In: Ying S-Y, editor. MicroRNA Protocols. Totowa, NJ: Humana Press; 2006. p. 129–38.
Huang J, Khademi M, Fugger L, Lindhe Ö, Novakova L, Axelsson M, et al. Inflammation-related plasma and CSF biomarkers for multiple sclerosis. Proc Natl Acad Sci U S A. 2020;117(23):12952–60.
Bi Z, Qiu P-F, Zhang Y, Song X-G, Chen P, Xie L, et al. A Three lncRNA Set. AC009975.1, POTEH-AS1 and AL390243.1 as nodal efficacy biomarker of neoadjuvant therapy for HER-2 positive breast cancer. Front Oncol. 2021;11:779140-.
Fortino V, Scala G, Greco D. Feature set optimization in biomarker discovery from genome-scale data. Bioinformatics. 2020;36(11):3393–400.
Yu H, Liu Y, He B, He T, Chen C, He J, et al. Platelet biomarkers for a descending cognitive function: a proteomic approach. Aging Cell. 2021;20(5):e13358.
Liu C, Yu Z, Huang S, Zhao Q, Sun Z, Fletcher C, et al. Combined identification of three miRNAs in serum as effective diagnostic biomarkers for HNSCC. EBioMedicine. 2019;50:135–43.
Sumida BH, Houston AI, McNamara JM, Hamilton WD. Genetic algorithms and evolution. J Theor Biol. 1990;147(1):59–84.
Trevino V, Falciani F. GALGO: an R package for multivariate variable selection using genetic algorithms. Bioinformatics. 2006;22(9):1154–6.
Katoch S, Chauhan SS, Kumar V. A review on genetic algorithm: past, present, and future. Multimedia Tools Appl. 2021;80(5):8091–126.
Fortin F-A, Rainville F-MD, Gardner M-AG, Parizeau M, Gagné C. DEAP: evolutionary algorithms made easy. J Mach Learn Res. 2012;13(70):2171–5.
Licursi V, Conte F, Fiscon G, Paci P. MIENTURNET: an interactive web tool for microRNA-target enrichment and network-based analysis. BMC bioinformatics. 2019;20(1):545.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41(D1):D991–5.
Hsu S-D, Lin F-M, Wu W-Y, Liang C, Huang W-C, Chan W-L, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39(Database issue):D163–9.
Huang H-Y, Lin Y-C-D, Li J, Huang K-Y, Shrestha S, Hong H-C, et al. miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 2020;48(D1):D148–54.
Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46(D1):D649–55.
Frost RJA, Olson EN. Control of glucose homeostasis and insulin sensitivity by the Let-7 family of microRNAs. Proc Natl Acad Sci U S A. 2011;108(52):21075–80.
Zhu H, Shyh-Chang N, Segrè AV, Shinoda G, Shah SP, Einhorn WS, et al. The Lin28/let-7 axis regulates glucose metabolism. Cell. 2011;147(1):81–94.
Agbu P, Cassidy JJ, Braverman J, Jacobson A, Carthew RW. MicroRNA miR-7 Regulates Secretion of Insulin-Like Peptides. Endocrinology. 2020;161(2):bqz040.
Fernández-de Frutos M, Galán-Chilet I, Goedeke L, Kim B, Pardo-Marqués V, Pérez-García A, et al. MicroRNA 7 impairs insulin signaling and regulates Aβ levels through posttranscriptional regulation of the insulin receptor substrate 2, insulin receptor, insulin-degrading enzyme, liver X receptor pathway. Mol Cell Biol. 2019;39(22):e00170-e219.
Gong Q, Xie Jn, Liu Y, Li Y, Su G. Differentially expressed MicroRNAs in the development of early diabetic retinopathy. J Diabetes Res. 2017:4727942.
Brady MJ. IRS2 takes center stage in the development of type 2 diabetes. J Clin Invest. 2004;114(7):886–8.
Kubota T, Kubota N, Kadowaki T. Imbalanced insulin actions in obesity and type 2 diabetes: key mouse models of insulin signaling pathway. Cell Metab. 2017;25(4):797–810.
Oliveira JM, Rebuffat SA, Gasa R, Gomis R. Targeting type 2 diabetes: lessons from a knockout model of insulin receptor substrate 2. Can J Physiol Pharmacol. 2014;92(8):613–20.
Cohen DH, LeRoith D. Obesity, type 2 diabetes, and cancer: the insulin and IGF connection. Endocr Relat Cancer. 2012;19(5):F27–45.
Dong L, Hou X, Liu F, Tao H, Zhang Y, Zhao H, et al. Regulation of insulin resistance by targeting the insulin-like growth factor 1 receptor with microRNA-122-5p in hepatic cells. Cell Biol Int. 2019;43(5):553–64.
Razny U, Polus A, Goralska J, Zdzienicka A, Gruca A, Kapusta M, et al. Effect of insulin resistance on whole blood mRNA and microRNA expression affecting bone turnover. Eur J Endocrinol. 2019;181(5):525–37.
Palsgaard J, Brøns C, Friedrichsen M, Dominguez H, Jensen M, Storgaard H, et al. Gene expression in skeletal muscle biopsies from people with type 2 diabetes and relatives: differential regulation of insulin signaling pathways. PLoS One. 2009;4(8):e6575.
Luo J, Sobkiw CL, Hirshman MF, Logsdon MN, Li TQ, Goodyear LJ, et al. Loss of class IA PI3K signaling in muscle leads to impaired muscle growth, insulin response, and hyperlipidemia. Cell Metab. 2006;3(5):355–66.
Huang X, Liu G, Guo J, Su Z. The PI3K/AKT pathway in obesity and type 2 diabetes. Int J Biol Sci. 2018;14(11):1483–96.
Molinaro A, Becattini B, Mazzoli A, Bleve A, Radici L, Maxvall I, et al. Insulin-Driven PI3K-AKT signaling in the hepatocyte is mediated by redundant PI3Kα and PI3Kβ activities and is promoted by RAS. Cell Metab. 2019;29(6):1400-9.e5.
Kursawe R, Eszlinger M, Narayan D, Liu T, Bazuine M, Cali AMG, et al. Cellularity and adipogenic profile of the abdominal subcutaneous adipose tissue from obese adolescents: association with insulin resistance and hepatic steatosis. Diabetes. 2010;59(9):2288–96.
Sales V, Patti M-E. The ups and downs of insulin resistance and type 2 diabetes: lessons from Genomic analyses in humans. Curr Cardiovasc Risk Rep. 2013;7(1):46–59.
Du X, Li X, Chen L, Zhang M, Lei L, Gao W, et al. Hepatic miR-125b inhibits insulin signaling pathway by targeting PIK3CD. J Cell Physiol. 2018;233(8):6052–66.
Acknowledgements
None.
Funding
This work was supported by the Gachon University research fund of 2020 (GCU-202008430003 to SN); the Industrial Technology Innovation program (NO. 20016417, AI prediction platform development for lung and gastric cancer with Korean genetic data and its servitization) funded by the Ministry of Trade, Industry & Energy (MOTIE (KATS)/KEIT, Korea); and the Gachon University Gil Medical Center (Grant number: FRD2023-13 to SN).
Author information
Authors and Affiliations
Contributions
Conceptualization, S.N.; methodology, S.N. and A.P.; formal analysis, A.P.; investigation, A.P.; data curation, A.P.; writing—original draft preparation, A.P.; writing—review and editing, S.N.; visualization, A.P.; supervision, S.N.; funding acquisition, S.N. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable as we used the publicly available datasets.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Description of public miRNA-Seq and gene expression array data used in this study. Table S2. GEO2R analysis of mRNAs (PIK3CD, IGF1R, IRS2, NRAS, and PIK3CB) targeted by the best miRNA biomarker set in skeletal muscle tissue (GSE22309). Table S3. Spearman correlation matrix between hepatic mRNA expression of the targeted genes in BxD mouse and clinical parameters related to type 2 diabetes. Table S4. Differential expression analysis of miRNAs (hsa-let-7b, hsa-miR-125b, and hsa-miR-7) in kidney tissues of T2DN vs. HC (GEO accession: GSE51674). Fig. S1. Prisma diagram of miRNA-Seq dataset acqusition in this study. In SRA, we acquired 2 datasets of type 2 diabetes study in blood tissue. Hence, 95 samples were used for biomarker discovery. Fig. S2. Workflow of miRDM-rfGA for biomarker discovery of optimal miRNA to classify T2DM and HC. (A) In the initial population, We generated 1000 individuals, and features in each individual are mapped to 1139 miRNA indices. miRDM-rfGA randomly select features among 1139 miRNA. (B) Fitness score is calculated for each individual based on AUROC score from RF classifier, then (C) miRDM-rfGA chooses the individual with the highest fitness score. The individual is included in the next generation and the other individuals are produced for next generation. Crossover and mutation are applied during the production of new individuals. (D) miRDM-rfGA iterates these phases for the number of generation G and derives the best individual among G generations. Fig. S3. Difference of the selected miRNA biomarkers expression levels between ‘diabetic nephropathy with type 2 diabetes’ (T2DN) and HC group. We obtained a publicly available miRNA expression dataset (GEO accession: GSE51674) to inspect the possibility of prediction of the risk of diabetic complications. We compared the selected miRNA expression levels between T2DN and HC groups.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Park, A., Nam, S. miRDM-rfGA: Genetic algorithm-based identification of a miRNA set for detecting type 2 diabetes. BMC Med Genomics 16, 195 (2023). https://doi.org/10.1186/s12920-023-01636-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-023-01636-2