Skip to main content

miRDM-rfGA: Genetic algorithm-based identification of a miRNA set for detecting type 2 diabetes

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.

Peer Review reports

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).

Fig. 1
figure 1

Overview of this study. T2DM and HC samples in blood were collected from NCBI GEO. Feature selection using GA with RF was used for miRNA biomarker discovery. We also compared the classification performance of T2DM with other traditional feature selection methods (F-test in ANOVA and Lasso)

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):

$$Fitness=100\times \frac{\sum_{k=1}^{M}{AUC}_{k}}{M}-W \times \left|x-b\right|$$
(1)

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].

Table 1 Description of miRDM-rfGA parameter settings and the mean AUROC score in each setting

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.

Table 2 Description of the traditional feature selection methods and the mean AUROC score in each setting

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.

Fig. 2
figure 2

Performance of classifying T2DM and HC using selected miRNA biomarker set from nine feature selection settings. A Settings 1, (B) 2, (C) 3, (D) 4, (E) 5, and (F) 6 were configured based on miRDM-rfGA. G Settings 7, (H) 8, and (I) 9 were configured using traditional feature selection methods (F-test and Lasso). In the test set, threefold cross validation was used, and mean AUROC and standard deviation was calculated

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.

Fig. 3
figure 3

Difference of miRNA biomarker expression levels between T2DM and HC group. We derived miRNA biomarker set using each setting for feature selection. We compared miRNA expression levels between T2DM and HC groups. A Settings 5 (the best setting), (B) 1, (C) 2, (D) 3, (E) 4, and (F) 6 were configured based on miRDM-rfGA. G Setting 7, and (H) settings 8 and 9 were configured using traditional feature selection methods (F-test and Lasso). The same miRNAs were selected in settings 8 and 9. Each miRNA expression was converted into z-score

Table 3 DE analysis of miRNA biomarkers derived by each setting

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].

Fig. 4
figure 4

The best miRNA biomarker set in setting 5. A Principal Component Analysis (B) Functional enrichment analysis of the miRNA biomarker set. C Pathway of IRS-related events triggered by IGF1R signaling. DE analysis result of the best miRNA biomarker set and the targeted mRNAs in skeletal muscle. The circle size of each point is proportional to the |log2 of fold change in T2DM vs. healthy control|

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).

Fig. 5
figure 5

Correlation and DE analysis between targeted mRNAs and clinical parameters related to T2DM, using public BxD mouse database: GeneNetwork. A Spearman’s correlation analysis between hepatic IRS2 expression (x-axis) in BxD high-fat diet (HFD) mice and Insulin (y-axis) during OGTT and (B) the comparison between the highest 25% mice (HFD) of IRS2 gene expression with the bottom 25% mice (HFD) of the gene expression. C Spearman’s correlation analysis between hepatic IRS expression in BxD chow diet (CD) mice and insulin response (IR) during OGTT and (D) the comparison between the top 25% mice (CD) of IRS2 gene expression with the bottom 25% mice (CD) of the gene expression. E Spearman’s correlation analysis between hepatic IGF1R expression level in BxD HFD mice and glycemia level during OGTT and (F) the comparison between the top 25% mice (HFD) of IGF1R gene expression with the bottom 25% mice of the gene. G Spearman’s correlation analysis between hepatic PIK3CD expression level in BxD HFD mice and IR during OGTT and (H) the comparison between the top 25% mice (HFD) of PIK3CD gene expression with the bottom 25% mice (HFD) of the gene

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

  1. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. 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.

    PubMed  PubMed Central  Google Scholar 

  3. 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-.

    PubMed  PubMed Central  Google Scholar 

  4. 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.

  5. 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.

    PubMed  PubMed Central  Google Scholar 

  6. 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.

    PubMed  Google Scholar 

  7. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. 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.

    PubMed  Google Scholar 

  9. 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.

    CAS  PubMed  Google Scholar 

  10. Eliasson L, Esguerra JLS. MicroRNA Networks in Pancreatic Islet Cells: Normal Function and Type 2 Diabetes. Diabetes. 2020;69(5):804–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    CAS  PubMed  Google Scholar 

  12. Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. 2007;8(2):93–103.

    CAS  PubMed  Google Scholar 

  13. 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.

    PubMed  Google Scholar 

  14. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Deng J, Guo F. MicroRNAs and type 2 diabetes. ExRNA. 2019;1(1):36.

    Google Scholar 

  16. Feng J, Xing W, Xie L. Regulatory Roles of MicroRNAs in Diabetes. Int J Mol Sci. 2016;17(10):1729.

    PubMed  PubMed Central  Google Scholar 

  17. 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.

    PubMed  PubMed Central  Google Scholar 

  18. 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.

    CAS  PubMed  Google Scholar 

  19. 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.

    PubMed  PubMed Central  Google Scholar 

  20. 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.

    PubMed  Google Scholar 

  21. He Z, Yu W. Stable feature selection for biomarker discovery. Comput Biol Chem. 2010;34(4):215–25.

    CAS  PubMed  Google Scholar 

  22. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Saeys Y, Inza I, Larrañaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23(19):2507–17.

    CAS  Google Scholar 

  24. 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.

    CAS  PubMed  Google Scholar 

  25. Dessì N, Pascariello E, Pes B. A comparative analysis of biomarker selection techniques. Biomed Res Int. 2013;2013:387673.

    PubMed  PubMed Central  Google Scholar 

  26. 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.

  27. 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.

    CAS  PubMed  Google Scholar 

  28. Leinonen R, Sugawara H, Shumway M. International nucleotide sequence database C. The sequence read archive. Nucleic Acids Res. 2011;39(Database issue):D19–21.

    Google Scholar 

  29. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal: Next Generat Sequencing Data Analys. 2011;17:1.

    Google Scholar 

  30. 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.

    PubMed  Google Scholar 

  31. Griffiths-Jones S. miRBase: The MicroRNA Sequence Database. In: Ying S-Y, editor. MicroRNA Protocols. Totowa, NJ: Humana Press; 2006. p. 129–38.

    Google Scholar 

  32. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 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-.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Fortino V, Scala G, Greco D. Feature set optimization in biomarker discovery from genome-scale data. Bioinformatics. 2020;36(11):3393–400.

    CAS  PubMed  Google Scholar 

  35. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Sumida BH, Houston AI, McNamara JM, Hamilton WD. Genetic algorithms and evolution. J Theor Biol. 1990;147(1):59–84.

    CAS  PubMed  Google Scholar 

  38. Trevino V, Falciani F. GALGO: an R package for multivariate variable selection using genetic algorithms. Bioinformatics. 2006;22(9):1154–6.

    CAS  PubMed  Google Scholar 

  39. Katoch S, Chauhan SS, Kumar V. A review on genetic algorithm: past, present, and future. Multimedia Tools Appl. 2021;80(5):8091–126.

    Google Scholar 

  40. 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.

  41. 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.

    PubMed  PubMed Central  Google Scholar 

  42. 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.

    CAS  PubMed  Google Scholar 

  43. 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.

  44. 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.

    Google Scholar 

  45. 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.

    CAS  PubMed  Google Scholar 

  46. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Agbu P, Cassidy JJ, Braverman J, Jacobson A, Carthew RW. MicroRNA miR-7 Regulates Secretion of Insulin-Like Peptides. Endocrinology. 2020;161(2):bqz040.

    PubMed  Google Scholar 

  49. 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.

    PubMed  PubMed Central  Google Scholar 

  50. 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.

  51. Brady MJ. IRS2 takes center stage in the development of type 2 diabetes. J Clin Invest. 2004;114(7):886–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 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.

    CAS  PubMed  Google Scholar 

  53. 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.

    CAS  PubMed  Google Scholar 

  54. Cohen DH, LeRoith D. Obesity, type 2 diabetes, and cancer: the insulin and IGF connection. Endocr Relat Cancer. 2012;19(5):F27–45.

    CAS  PubMed  Google Scholar 

  55. 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.

    CAS  PubMed  Google Scholar 

  56. 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.

    CAS  PubMed  Google Scholar 

  57. 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.

    PubMed  PubMed Central  Google Scholar 

  58. 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.

    CAS  PubMed  Google Scholar 

  59. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. 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.

    CAS  PubMed  Google Scholar 

  61. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. 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.

    PubMed  Google Scholar 

  63. 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.

    CAS  PubMed  Google Scholar 

Download references

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

Authors

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

Correspondence to Seungyoon Nam.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-023-01636-2

Keywords