Next Article in Journal
Adapting Multiple Distributions for Bridging Emotions from Different Speech Corpora
Next Article in Special Issue
An Information-Theoretic Bound on p-Values for Detecting Communities Shared between Weighted Labeled Graphs
Previous Article in Journal
Dynamic Multiscale Information Spillover among Crude Oil Time Series
Previous Article in Special Issue
A Comparison of Partial Information Decompositions Using Data from Real and Simulated Layer 5b Pyramidal Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis

1
School of Mathematics and Statistics, Henan University of Science and Technology, Luoyang 471023, China
2
College of Agriculture, Henan University of Science and Technology, Luoyang 471023, China
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(9), 1249; https://doi.org/10.3390/e24091249
Submission received: 19 May 2022 / Revised: 18 August 2022 / Accepted: 2 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Information Theory in Computational Biology)

Abstract

:
Type 2 diabetes mellitus (T2DM) is a metabolic disease caused by multiple etiologies, the development of which can be divided into three states: normal state, critical state/pre-disease state, and disease state. To avoid irreversible development, it is important to detect the early warning signals before the onset of T2DM. However, detecting critical states of complex diseases based on high-throughput and strongly noisy data remains a challenging task. In this study, we developed a new method, i.e., degree matrix network entropy (DMNE), to detect the critical states of T2DM based on a sample-specific network (SSN). By applying the method to the datasets of three different tissues for experiments involving T2DM in rats, the critical states were detected, and the dynamic network biomarkers (DNBs) were successfully identified. Specifically, for liver and muscle, the critical transitions occur at 4 and 16 weeks. For adipose, the critical transition is at 8 weeks. In addition, we found some “dark genes” that did not exhibit differential expression but displayed sensitivity in terms of their DMNE score, which is closely related to the progression of T2DM. The information uncovered in our study not only provides further evidence regarding the molecular mechanisms of T2DM but may also assist in the development of strategies to prevent this disease.

1. Introduction

The development of many diseases, including T2DM, can be regarded as a nonlinear dynamic process which is generally divided into three states: normal state, pre-disease state/critical state, and disease state. The normal state is stable, characterized by stability and robustness, whereby any changes occur slowly. The disease state is a new stable state that represents a phase of deterioration in which more obvious signs and symptoms of the disease appear and many patients begin to receive treatment, but it is difficult to return to the normal state [1]. The pre-disease state/critical state corresponds to a critical point before the system transitions to an irreversible disease state and is accompanied by a drastic change in system dynamics [2]. In this state, the system is usually reversible, and reversion to the normal state is possible if there is appropriate treatment. Many diseases, including T2DM, are diagnosed upon reaching the disease stage. Therefore, it is important to identify the pre-disease state/critical state. However, this is a difficult task because there may be little difference between the normal and the pre-disease state/critical state to permit their clear distinction [3].
Traditional biomarkers can be used to distinguish disease states from normal states mainly based on the differential expression of individual molecules or a group of molecules, but such biomarkers cannot be used to detect the pre-disease state/critical state owing to their static nature [2]. In order to solve the challenge of detecting critical points, a new concept, called dynamic network biomarker (DNB), was proposed [3]. The expression of DNBs reflects the presence or severity of the disease state, and they are required to have constant values that are different in the respective disease and normal states such that they can be used as a signal for detection in the early stages of development of a complex disease. When a biological system in a normal state approaches the critical state, there are three statistical conditions that allow detection of a critical point based on DNBs: correlations between any two members of a DNB group rapidly increase, correlations between one member of a DNB group and any other non-DNB molecule rapidly decrease, and standard deviations of the variables in the DNB group drastically increase [4]. Based on these three necessary conditions, the DNB method can be used to identify a critical state. At present, research on combining the DNB method and information theory to detect the critical state of complex diseases is attracting attention. Based on the hidden Markov model (HMM), an inconsistent index algorithm was proposed to identify the critical state before disease deterioration [5]. Single-cell graph entropy (SGE) can explore the gene–gene associations among cell populations based on single-cell RNA sequencing data [6]. The single-sample-based Jensen–Shannon divergence (sJSD) method is used to detect the early-warning signals of complex diseases before critical transitions based on individual single-sample data [7]. The temporal network flow entropy (TNFE) method, which is based on network fluctuation of molecules, can detect the critical states of complex diseases on the basis of each individual [8].
In this study, we propose a new method, i.e., the degree matrix network entropy (DMNE) algorithm, based on sample-specific network (SSN) and network entropy. Dai et al. presented a new method to construct a cell-specific network (CSN) for each single cell, which transforms the data from “unstable” gene expression form to “stable” gene association form on a single-cell basis. This can also be directly applied to construct an individual network of each single sample [9]. Therefore, we applied this method to construct an SSN for the sample data and propose a method combining DMNE with network entropy that does not use single-cell data. Firstly, SSNs of reference samples and perturbed samples are constructed, respectively, and the network degree matrix (NDM) and local network are obtained. Then, the difference between the DMNE of the reference sample and the perturbed sample is calculated to quantify the differential network information flow. Finally, the DMNE score is calculated to characterize the molecular collective fluctuation or network fluctuation caused by case samples. By transforming floating gene expression datasets into stable network entropy, the DMNE method provides a new method to detect the critical state based on the SSN, which has the following advantages: (i) By calculating the DMNE score based on the difference of the network entropy at each stage, the DMNE method exhibits the differential changes in diseases at the network level. (ii) Based on the DMNE method, the critical states can be successfully detected, and the DNBs of the critical state can be effectively identified. (iii) Based on the DMNE method, “dark genes” are found, which play important roles during the development of T2DM, and drug targets against these are obtained.

2. Materials and Methods

2.1. Data Progression and Functional Analysis

The DMNE method was applied to detect the critical transitions during T2DM development and progression accompanying insulin resistance in adipose tissue, gastrocnemius muscle, and liver for rats (diabetes rats: Goto–Kakizaki (GK) rats; control rats: Wistar–Kyoto (WKY) rats). The high-throughput experimental datasets were downloaded from the NCBI GEO database (access ID: GSE13268, GSE13269, and GSE13270) (www.ncbi.nlm.nih.gov/geo (accessed on 1 October 2021)), comprising time-course gene expression data obtained from age-specific rats corresponding to well-designed time series (five samples each from 4, 8, 12, 16, and 20 weeks). These experimental datasets were provided by Almon et al. [10]. All of the GEO samples we used are provided in Table A1. In addition, pancreatic adenocarcinoma (PAAD) survival data and the gene expression values in different samples were downloaded from the cancer genome atlas (TCGA) database (https://cancergenome.nih.gov/ (accessed on 1 March 2022)).
The enrichment analysis of DNBs is based on Gene Ontology Consortium [11] (GOC, http://geneontology.org (accessed on 1 February 2022)), DAVID Bioinformatics Resources [12] (https://david.ncifcrf.gov/ (accessed on 1 February 2022)), and Circos [13] (http://www.circos.ca/ (accessed on 20 February 2022)). Protein–protein interaction (PPI) networks were drawn by STRING (https://string-db.org/ (accessed on 20 February 2022)) and the client software Cytoscape (https://cytoscape.org/ (accessed on 25 February 2022)). The drug targets were downloaded from DrugBank [14] (https://go.drugbank.com/ (accessed on 20 January 2022)), Therapeutic Target Database [15] (TTD, http://db.idrblab.net/ttd/ (accessed on 20 January 2022)), and Pharmacogenomics Knowledgebase [16] (PharmGKB, https://www.pharmgkb.org/ (accessed on 20 January 2022)).

2.2. Theoretical Background

DNBs form an observable subnetwork or molecular group in a disease system or disease network. They can be used to identify the critical state/pre-disease state before the sudden deterioration of complex diseases [17]. When the system approaches the critical state from the normal state, DNB molecules (i.e., genes or proteins) have the following three statistical properties [8]:
  • The correlation between any two members of a DNB group rapidly increases;
  • The correlation between one member in a DNB group and any other non-DNB molecule sharply decreases;
  • The standard deviation for any member in the DNB group drastically increases.
These theoretical results have allowed many achievements in detecting critical points and in discussing related biological processes. According to these three properties, DNBs can be reliably identified, and significant early warning signals can be extracted.
The CSN method [9] provides a promising strategy for analyzing genes and gene associations at the single-cell level. This can also be applied to non-single-cell datasets for constructing an individual network of each single sample in a similar way [9]. We thus applied CSN to each sample to construct an SSN of reference samples (the healthy control samples) and perturbed samples (the diabetic samples), respectively. Therefore, network fluctuations can be quantified based on the differences. We propose a network entropy algorithm based on SSN to predict the critical state during the development of T2DM.

2.3. Data Preprocessing

We analyzed three preprocessed gene expression datasets measured in adipose, gastrocnemius muscle, and liver from diabetic rats (GK) and control rats (WKY) provided the same normal diet. The datasets were downloaded from the NCBI GEO database (access ID: GSE13268, GSE13269, and GSE13270). The three datasets were not raw data, and they were provided by Almon et al. [10]. The datasets comprised 31,099 probes, measured using the Affymetrix Microarray Suite 5.0 (Affymetrix), for 25 WKY controls and 25 spontaneously diabetic GK rats at 4, 8, 12, 16, and 20 weeks of age [10].
Then, we transformed the downloaded matrixes into required gene expression matrixes (GEMs) through ID conversion and the deletion of duplicate genes and null values. In the case of multiple probes corresponding to the same gene, the values were individually averaged for each to obtain three GEMs containing 15,246, 15,343, and 15,246 genes. Finally, the logarithm log(1 + x) was applied to normalize the three GEMs.

2.4. Algorithm to Detect the Tipping Point and Identify DNBs of T2DM Based on DMNE

After extracting some reference samples (samples from a normal cohort that are regarded as the background that represents the healthy or control individuals), we then developed and applied the following algorithm to detect the tipping point using the data of the diabetic case samples and healthy control samples. The algorithm flow is shown in Figure 1.
The steps of the degree matrix network entropy (DMNE) algorithm are specified as follows:
[Step 1] Construct the SSN of reference samples and perturbed samples, respectively.
Based on n reference samples { S 1 , S 2 , , S n } , the R S S N (the SSN of reference sample) is constructed as follows.
Firstly, we normalize the initial gene expression matrix with m rows/genes and n columns/samples and average the expression for the same gene in n samples. Secondly, we plot scatter diagrams for every two genes in a rectangular coordinate system, where the horizontal axis and the vertical axis are the expression value E i of gene g i and the expression value E j of gene g j , respectively, and these two genes constitute a gene pair ( g i , g j ) . m genes lead to m(m − 1)/2 scatter diagrams. Then, we draw two boxes near E i and E j . According to the research of Dai et al. [9], the influence of box size and p value in different datasets on the clustering effect is tested from the viewpoint of clustering. It is indicated that the optimum box size is about 0.1, and the optimum p value is about 0.01 on average, which are set as the default parameters of the CSN method. The specific test results are given in Supplementary Materials Section A. We choose 0.1 n, which is proportional to the number of samples n, and we can adjust the parameters according to specific needs. The boxes are drawn according to a predetermined integer, and the two boxes overlap to obtain the third box. The value is obtained by calculating the number of points (samples) in the third box. Whether there is an edge between the two genes in the scatter diagram of genes g i and g j is determined by the following statistical dependence index:
r i , j = n ( E i , E j ) n n ( E i ) n n ( E j ) n
where n ( E i ) , n ( E j ) and n ( E i , E j ) represent the number of points (samples) in the vertical box, horizontal box, and overlapping box, respectively. If the statistical dependency index is greater than zero, there is an edge between g i and g j ; otherwise, there is no edge. That is, the R S S N is constructed based on the specific network of the reference samples.
Similarly, the case samples at different time points t = T are added to the reference samples to obtain the disturbed samples { S 1 , S 2 , , S n , S c a s e T } , and the P S S N T (the SSN of the perturbed samples at time point T) can be constructed.
[Step 2] Extract degree matrix and local network.
The NDM of reference samples and disturbed samples can be obtained from R S S N and P S S N T , respectively.
Extract each reference local network R L N d ( d = 1 , , l ) from RSSN. R L N d is composed of genes g i d ( i = 1 , 2 , , M d ) (the central gene is g 1 d , and the rest are neighborhood genes).
As the genes in the reference sample and the perturbed sample are the same, similarly, the perturbed local networks P L N T d ( d = 1 , , l ) are extracted from P S S N T . The P L N T d is composed of genes g i d ( i = 1 , 2 , , M d ) (the central gene is g 1 d , and the rest are neighborhood genes).
[Step 3] Calculate the node probability for every gene by fitting a Gaussian distribution.
Based on the NDM of reference samples, the Gaussian distributions of each gene g i d ( i = 1 , 2 , , M d ) in local network R L N d are fitted. Node probabilities p r e f e r e n c e ( x i d ) of gene g i d are calculated as follows:
p r e f e r e n c e ( x i d ) = 1 σ 2 π x i d e ( u μ i r d ) 2 2 σ i r d 2 d u i = 1 M d 1 σ 2 π x i d e ( u μ i r d ) 2 2 σ i r d 2 d u ,
where x i d is the degree value from NDM of gene g i d in the reference sample, and μ i r d and σ i r d are the mean value and standard deviation of degree value for gene g i d in the reference samples, respectively.
Similarly, based on the NDM of perturbed samples, node probabilities p p e r t u r b e d ( x i T d ) of gene g i d are, respectively, calculated as follows:
p p e r t u r b e d ( x i T d ) = 1 σ 2 π x i T d e ( u μ i p d ) 2 2 σ i p d 2 d u i = 1 M d 1 σ 2 π x i T d e ( u μ i p d ) 2 2 σ i p d 2 d u ,
where x i T d is the degree value from NDM of gene g i d in the perturbed sample at time point t = T, and μ i p d and σ i p d are the mean value and standard deviation of the degree value for gene g i d in perturbed samples at time t = T, respectively.
[Step 4] Calculate the network entropy of R S S N and P S S N T .
Firstly, the network entropy of R L N d ( d = 1 , , l ) is calculated as follows:
H n ( x i d ) = i = 1 M d x i d p r e f e r e n c e ( x i d ) log ( x i d p r e f e r e n c e ( x i d ) ) M d .
Then, the network entropy of the global network of R S S N can be calculated as
H n = d = 1 l H n ( x i d ) l .
The network entropy of P L N T d ( d = 1 , , l ) is
H T n + 1 ( x i T d ) = i = 1 M d x i T d p p e r t u r b e d ( x i T d ) log ( x i T d p p e r t u r b e d ( x i T d ) ) M d .
Then, the network entropy of the global network of P S S N T is
H T n + 1 = d = 1 l H n ( x i T d ) l .
[Step 5] Calculate the DMNE score to quantify the network differences caused by the perturbed sample.
The DMNE of reference samples and perturbed samples is calculated as follows:
Δ H T = H T n + 1 H n
where the Δ H T score reflects the global perturbation caused by the case samples at each time point t = T. The higher the score, the greater the difference between the case samples and the reference samples. The sudden increase in Δ H T can be considered as an early warning signal of important changes in the process of disease progression.

3. Results

3.1. Detecting the Critical State of T2DM

Here, we apply the DMNE method to detect the critical transitions during T2DM development and progression accompanying insulin resistance in rat adipose, gastrocnemius muscle, and liver corresponding to the GSE13268, GSE13269, and GSE13270 datasets, respectively, which each contain 25 reference and 25 case samples. In the case samples, there are five samples at each stage. At each stage, the DMNE score is calculated. The top 200 genes (1.5% of all expressed genes) with the highest DMNE score at the critical state are regarded as DNBs.

3.1.1. The Critical State of GSE13268

As shown in Figure 2a, the dramatic increase in the DMNE score for GK rat adipose appeared at 8 weeks, which indicates the upcoming critical transition. To show the DMNE scores in a local view, the landscape of the DMNE scores of DNBs is illustrated in Figure 2d. It can be seen that around 8 weeks, there is a group of genes whose DMNE scores abruptly increase. This critical phenomenon results from the drastic increase in the correlations between molecules in this group when the system approaches the tipping point. In Figure 2g, we illustrate the evolution of the top DMNE gene group/module, i.e., the protein–protein interaction (PPI) network of DNBs. This figure shows that a significant change in the network structure occurs at 8 weeks, signaling the critical transition into disease state from the molecular network level. The DNBs are given in Supplementary Materials Table S1.

3.1.2. The Critical State of GSE13269

As shown in Figure 2b, two critical transitions are detected for muscle tissue during T2DM. The dramatic increases in the DMNE score for GK rats appeared at 4 and 16 weeks, which indicate the upcoming critical transitions. To show the DMNE scores in a local view, the DMNE scores of DNBs are illustrated in Figure 2e. This figure shows that around 4 and 16 weeks, there is a group of genes whose DMNE scores abruptly increase. This critical phenomenon results from the drastic increase in the correlations between molecules in this group when the system approaches the tipping point. In Figure 2h, we illustrate the evolution of the top DMNE gene group/module, i.e., protein–protein interaction (PPI) network of DNBs. This figure shows that the significant changes in the network structure occur at 4 and 16 weeks, signaling the critical transitions into disease state from the molecular network level. The DNBs are given in Supplementary Materials Table S1.

3.1.3. The Critical State of GSE13270

As shown in Figure 2c, there are also two critical transitions during T2DM consistent with GSE13269. The dramatic increases in the DMNE score for GK rat livers appeared at 4 and 16 weeks, which indicate the upcoming critical transitions. To show the DMNE scores in a local view, the DMNE scores of DNBs are illustrated in Figure 2f. It can be seen that around 4 and 16 weeks, there is a group of genes whose DMNE scores abruptly increase. This critical phenomenon results from the drastic increase in the correlations between molecules in this group when the system approaches the tipping point. In Figure 2i, we illustrate the evolution of the top DMNE gene group/module, i.e., the protein–protein interaction (PPI) network of DNBs. This figure shows that the significant changes in the network structure occur at 4 and 16 weeks, signaling the critical transitions into disease state from the molecular network level. The DNBs are given in Supplementary Materials Table S1.

3.2. Tissue-Specific Analysis

To confirm the significant relation between tissue-specific DNBs and T2DM progression, the GO analysis and KEGG enrichments of DNBs are conducted for each DNB to categorize the genes participating in different biological functions or pathways, as shown in Table 1, Table 2 and Table 3 and Figure 3. In addition, “housekeeping” genes are continuously highly expressed, so their role in any disease is negligible. We found some “housekeeping” genes contained in the DNBs we obtained, such as GADPH, Ldha, and Arhgdia [18]. These genes are not considered in the subsequent analysis [18].
For every dataset, we use two kinds of GO backgrounds for analysis, one of which includes all Rattus norvegicus genes, and the other is only the genes that pass the previous filtering steps (the genes in the NDM). We show the two results of GSE13268 in Table 1 and Table 2, which are analyzed separately. For the GSE13269 and GSE13270 datasets, we only show the results for genes that passed the previous filtering steps as background in Table 3. The results for all genes in GSE13269 and GSE13270 are in Supplementary Materials Table S5 and Section S6.

3.2.1. Analysis of GSE13268

For GK rat adipose, the enrichment results of DNBs based on all genes of the species are provided in Table 1. The biological processes at 8 weeks mainly include response to lipid, stimulation, and insulin secretion regulation. Insulin resistance and inflammation are linked with T2DM and related diseases. Abnormal insulin secretion causes hyperglycemia, which will stimulate the development of hyperinsulinemia, reduce the number of insulin receptors, and aggravate insulin resistance [19]. Moreover, insulin resistance in T2DM is accompanied by the dysfunction of glucose and lipid metabolism as well as protein biosynthesis. Therefore, the related biological processes of DNB genes on T2DM also include the regulation of protein phosphorylation and the response to peptide hormones. For cellular calcium ion homeostasis, calcium ion participates as an activator of lipid metabolism [2]. For regulation of the erk1 and erk2 cascade, the increase in DPP-4 content in diabetic mice specifically activates the ERK1/2 and NF-κB signaling pathways and promotes the calcification of their aortic vessels [20]. Renal artery calcification can enhance progressive renal damage in rats with T2DM nephropathy, further leading to hyperlipidemia [21,22].
For some DNB genes, NOS3 polymorphisms are associated with the progression of kidney and cardiovascular disease in Type 2 diabetic patients. FABP1 regulates the absorption and transportation of fatty acids in the liver by promoting the transportation, storage, and utilization of fatty acids and their acyl-CoA derivatives. The overexpression of FABP1 can disrupt the clearance function in autophagy by inhibiting lysosome functions (including lysosomal protein decomposition and lysosomal acidification maintenance), which promotes liver steatosis [23]. The enzyme subunit expressed by the mitochondrial ND3 gene is an important part of respiratory chain complex I. Deficiency of the ND3 gene in β cells will lead to decreased respiratory chain complex I activity, which will result in islet β cell dysfunction and allow for the possibility of developing diabetes [24].
The enrichment results for DNBs based on the genes that pass the previous filtering steps as background are presented in Table 2. The enrichment pathways represented by DNBs mainly include response to lipopolysaccharide, protein kinase binding, endopeptidase activity, oxidation–reduction process, and response to insulin. Statistically significant enrichment is observed in all cases, with p values less than 0.05. The relationship between abnormal insulin secretion and T2DM has been explained in previous analyses. Here, the results for response to lipopolysaccharide and protein kinase binding are the same as in the above analysis [19]. In addition, we also obtained new pathways, such as inflammatory response. Either in tissue-specific analysis or in cross-tissue analysis, lipid metabolism that appears to be abnormal or inflammatory is a common functional cascade associated with T2DM [2].
Comparing the analysis results against two different backgrounds, it can be seen that although there are some differences in the pathways identified through enrichment, the types of pathways are basically the same. The DNBs of adipose tissue are sensitive to metabolic processes and lipid responses.
By setting two different backgrounds for analysis, both demonstrate that the DNBs obtained by DMNE are closely related to the development of T2DM.

3.2.2. Analysis of GSE13269

For GK rat muscle, the related biological processes at 4 weeks mainly include muscle contraction, glycolytic process, peptidase activator activity, response to glucose, T cell co-stimulation, cellular response to glucose stimulus, glucose homeostasis, and the PI3K/Akt signaling pathway, which are provided in Table 3. The biological processes at 16 weeks mainly include regulation of peptidase activity, negative regulation of phosphorylation, carbohydrate metabolic process, glycolytic process, T cell co-stimulation, regulation of protein catabolic process, the HIF-1 signaling pathway, and the PI3K/Akt signaling pathway, which are also provided in Table 3. B/T cells actively participate in the inflammatory response in the first critical period for muscle tissues [25,26]. For the PI3K/Akt signaling pathway, the gene Sirt1 participates in diabetic myocardial injury and the occurrence and development of early diabetic cardiomyopathy by negatively regulating the PI3K/Akt/MTOR signaling pathway [27]. Through cell transfection and other experiments, previous studies have shown that the HIF-1α/KIM1 signaling pathway can participate in the process of renal fibrosis in diabetic nephropathy by regulating the KIM1 expression in renal tubular epithelial cells under a high-glucose environment. HIF-1α is the key transcription factor of oxygen homeostasis regulation, which can make the body adapt to the external environment oxygen by regulating the expression of target genes. The expression levels of HIF-1α, KIM1, and COL-1mRNA and protein controlled by them in the kidneys of diabetic rats are significantly increased, suggesting that the diabetic nephropathy model is accompanied by changes in the HIF-1α/KIM1 pathway in kidneys [28]. For some DNB genes, PRDX2 has been found to participate in the oxidative stress process through a variety of signaling pathways. The environment of hypoxia that it induces leads to an increase in reactive oxygen species. In order to survive in a high level of reactive oxygen species, cells need to increase antioxidant levels, so the expression of related factors, including PRDX2 protein, is altered. Previous studies suggest that PRDX2 overexpression can prevent pancreatic cell apoptosis induced by oxidative stress and reduce the risk of developing diabetes [29]. Some studies [30,31] have suggested that SPARC is an autocrine and/or paracrine factor of adipose tissue that can inhibit fat formation. Other studies [32] have shown that insulin can increase the expression of SPARC, which can increase the phosphorylation of AKT and PI3-K, indicating that SPARC may be involved in insulin signaling pathways [33].

3.2.3. Analysis of GSE13270

For GK rat liver, the biological processes at 4 weeks mainly include inflammatory response, response to hormone, oxidoreductase activity, catalytic activity, response to interleukin-1, response to lipid, and metabolism of xenobiotics by cytochrome P450 and the AGE/RAGE signaling pathway in diabetic complications, which are provided in Table 3. The biological processes at 16 weeks mainly include liver development, lipid binding, negative regulation of immune response, the NOS2-CD74 complex, vitamin A metabolic process, the PPAR signaling pathway, and the AGE/RAGE signaling pathway in diabetic complications, which are also provided in Table 3. In the first critical state, the metabolism of xenobiotics by the cytochrome P450 pathway plays an important role in the oxidation of organic substances, whose involved enzymes in liver could generally act as metabolic intermediates, e.g., lipids and steroid hormones [34,35]. Vitamin A and its metabolites can inhibit the expression of resting and activated pancreatic stellate cells α-SMA [36]. The peroxisome proliferator-activated receptors (PPARs) modulate several biological processes that are perturbed in obesity, including inflammation, lipid and glucose metabolism, and overall energy homeostasis. PPARs regulate the functions of adipose tissues, such as adipogenesis, lipid storage, and adaptive thermogenesis [37]. For the AGE/RAGE signaling pathway, moderate-intensity aerobic exercise inhibits the AGE/RAGE axis and the NF-κB pathway, which may decrease oxidative stress and inflammation and thus reduce tissue injury for the prevention and treatment of T2DM complications [38]. For some DNB genes, CD74 receptor and other mechanisms are stimulated by multi-effect cytokine MIF to promote inflammatory response in glomerular podocytes. MIF can activate CD74 on the surface of local glomerular podocytes, resulting in the phosphorylation of extracellular signal-regulated kinase 1/2 (ERK1/2) and p38MAPK [39]. SPARC is an autocrine and paracrine factor of adipose tissue that can inhibit adipogenesis. Insulin can increase SPARC expression, while SPARC can increase Akt and PI3-K phosphorylation, indicating that SPARC may be involved in insulin signaling pathways [40]. PRDX2 overexpression can prevent pancreatic β cell apoptosis induced by oxidative stress and reduce the risk of developing diabetes. The changes in oxidative stress in diabetic kidney disease and other pathological processes are related to the expression of PRDX2 and related pathways. For example, PRDX2 can participate in the regulation of tumor development and oxidative stress in treatment through the PI3/AKT-resistant pathway. PRDX2 can regulate the oxidative stress of colon cancer cells through the Wnt/B-catenin signaling pathway and influence the oxidative stress of tumor by regulating some microRNAs [29].

3.3. Cross-Tissue Analysis

From the above enrichment analysis, we observe the critical point of tissue specificity. There are different degrees of differences in the ways DNBs participate in T2DM, but there is also relative consistency. We can infer that the phenotypic changes in GK rats share a cross-organizational functional relationship during the development of T2DM.
In the pre-disease period, the high blood glucose is observed in GK rats fed with the same food in the healthy control WKY rats, although the insulin secretion level in GK rats is similar to the level in WKY rats. Hence, it can be inferred that the insulin receptors in GK rats could not correctly respond to the insulin regulation, such that the glucose cannot be transported successfully in cells. In the preceding analysis, we find that there are abnormal responses to lipids in the early warning stage of lipid tissue and the first early warning stage of liver tissue, and abnormal responses of the PPAR signaling pathway are found in the fourth early warning stage of liver tissue. This pathway is activated by upstream FABP family proteins and is then involved in fatty acid transport. This may explain the remarkable abnormal lipid metabolism in the early stage of the disease due to a lack of glucose intake.
In the pre-transition period, hyperglycemia remains comparatively stable, but the concentration of plasma insulin decreases. The ability to absorb and utilize glucose is reduced. Consequently, other metabolic processes in tissues, such as lipid metabolism and protein metabolism, are activated. Many genes related to lipid transport are dysfunctional in the liver, where disorders in protein digestion and absorption can first be detected. The occurrence of β cell failure at 20 weeks may result from the increasingly serious lack of glucose uptake. The abundance of enrichment pathways in various tissues indicates that with the further development of the disease, complications will increase, and functional proteins will encounter obvious obstacles in their operation, such as in the examples of negative regulation of proteoglycan and regulation of phosphorylation in cancer. Drug metabolism and vitamin metabolism are also abnormal. In the above analysis, it can also be found that the response to stimulation appears throughout the entire development process of T2DM.

3.4. DMNE Reveals Non-Differential “Dark Genes”

In clinical practice and scientific research, differentially expressed genes draw much attention in the early diagnosis of disease, screening drug targets, treating diseases, and developing new drugs. However, some non-differentially expressed genes in the coding region of DNA are called “dark genes” [41]. Based on the DMNE method, we found some “dark genes” that did not exhibit differential expression but whose DMNE scores were especially sensitive. Such genes are usually ignored in traditional analyses.
Many epidemiological studies have found that T2DM is positively associated with an increased risk of PAAD. Long-term Type 2 diabetes increases the risk of PAAD by 1.5 to 2 times. Thirty to forty percent of pancreatic cancer patients have diabetes, and eighty percent have abnormal glucose tolerance [42]. Therefore, some “dark genes” related to the prognosis of PAAD may play important roles in the development of T2DM, such as COL1A1, which is the most significant gene in the extracellular matrix receptor interaction pathway and is linked to hypoglycemic activity [43], which may affect blood glucose fluctuations, T2DM prognosis, and the occurrence and development of chronic complications [44].
To further confirm the effectiveness of the “dark genes”, we subjected them to GO and KEGG enrichments analyses. BCKDHA exists in liver tissue and encodes the BCKDC component. In the study of [45], it was found that knocking out the BCKDK gene, which regulates BCKDC kinase, can significantly reduce the concentration of branched-chain amino acid (BCAA) in blood. BCAA can stimulate mammalian targets of rapamycin and S6 kinase. It can also phosphorylate insulin receptor substrate on serine residues, thus damaging insulin signaling pathways. The decrease in BCAA catabolism will lead to insulin resistance and glucose intolerance, eventually leading to T2DM. Therefore, the increase in BCKDHA expression will cause a decrease in plasma BCAA concentration, subsequently affecting the progression of T2DM, which can explain why the low BCKDHA expression is associated with favorable prognosis [43]. COL1A1 contributes to T2DM via the ECM–receptor interaction pathway. ECM–receptor interactions are microenvironmental regulators of the structure and function of cells and tissues. ECM synthesis is closely linked to the PI3K/Akt pathway and the regulation of T2DM [46,47]. The enrichment analysis of these “dark genes” also showed that their associated pathways provided in Figure 4 are related to T2DM. Figure 4 shows the enrichment pathways of “dark genes” and the number of involved genes.
As pancreatic adenocarcinoma (PAAD) and T2DM are interrelated, and this association could inform our understanding of the mechanisms of both diseases, this relationship has become a topic of interest receiving considerable attention in the literature [42,48,49,50]. At present, due to the lack of corresponding prognostic data in the study of T2DM, the prognostic analysis of T2DM cannot be directly performed. Therefore, we used the PAAD data to analyze the prognosis of T2DM, as detailed in Supplementary Materials Section S2, representing an attempt to reveal the mechanism of T2DM development from another approach.

3.5. Drug Targets

Identifying drug targets is a key task in drug discovery and chemical genomics. In this study, we collected information on the targets and anticancer drugs for T2DM from the online service web pages Drug Bank, Therapeutic Target Database, and Pharmacogenomics Knowledgebase. The number of drug targets in the first, second, third, fourth, and fifth stages of the data of GSE13268 is 20, 29, 21, 23, and 19, respectively. The number of drug targets in the first, second, third, fourth, and fifth stages of the data of GSE13269 is 17, 20, 19, 20, and 15, respectively. The number of drug targets in the first, second, third, fourth, and fifth stages of the data of GSE13270 is 20, 35, 33, 37, and 29, respectively. We integrated the drug target information for critical states (phase 2 of GSE13268, phases 1 and 4 of GSE13269, and phases 1 and 4 of GSE13270).
The critical state of GSE13268 is the second stage, the critical states of GSE13269 are the first and fourth stages, and the critical states of GSE13270 are the first and fourth stages. Figure 5 shows the drug target analysis for the critical states of GSE13270. The drug target analyses [51,52,53,54,55] of GSE13268 and GSE13269 are provided in Supplementary Materials Sections S3 and S4.
In Figure 5, the red circle is the abbreviation of Type 2 diabetes, the blue circle is the drug target of the corresponding stage in the three sets of data of Type 2 diabetes, the green circle is the corresponding target protein, and the yellow circle is the corresponding anticancer drug of the protein that has been approved by the FDA or used in experiments. Infographics of drug targets for other stages of GSE13270 are provided in Supplementary Materials Section S5.
In the first stage of GSE13270, according to the analysis of MaoA and MaoB protein expression in the islets of Type 2 diabetes mellitus in mice and humans, it can be found that the amount of MaoB in β cells of Type 2 diabetes mellitus is significantly reduced, which shows that the loss of Mao causes dysfunction of β cells, resulting in Type 2 diabetes [56]. In other studies, a relationship between LGALS2 and glucose and fasting insulin has also been found. It is this relationship that supports the idea that LTA plays a role in determining the insulin–glucose spectrum. LGALS2 genotypes and LTA pathways may play other roles in insulin metabolism which require further research and exploration [57]. In the fourth stage of GSE13270 (liver tissue), FABP2, the intestinal fatty-acid-binding protein, is considered a candidate gene for diabetes. The reason is that, in the process of the absorption and metabolism of fatty acids (FA), FABP2 is involved in encoding proteins. Therefore, FABP2 may affect insulin sensitivity and glucose metabolism [58]. The gene-encoding alpha2 Heremans Schmid glycoprotein (AHSG) is a candidate for Type 2 diabetes and metabolic syndrome, and variants have been successfully demonstrated in previous attempts to be associated with Type 2 diabetes and obesity when examining Swedish patients, as well as Caucasian patients in France [59].

4. Discussion

In this study, we propose a new method, the DMNE algorithm, to detect the critical state and identify DNBs of T2DM. By applying this method to three rat gene expression datasets involving adipose tissue, gastrocnemius muscle, and liver tissue in healthy control vs. Type 2 diabetes conditions, the early warning signals of T2DM were successfully detected and DNBs were effectively identified. Specifically, for the three tissues, the DMNE method can not only detect the critical state via the DNBs but also identify drug targets. We used two kinds of GO backgrounds for analysis, and both demonstrated that the DNBs obtained by DMNE are closely linked to the development of T2DM. In addition, we mined “dark genes” that are associated with T2DM. In contrast to the information of differential expression used in traditional biomarkers to diagnose disease, DMNE is based on molecular network fluctuation and is able to predict disease.
We detected the critical state by exploiting the high-dimensional information of time series data. Firstly, by constructing the SSN, we converted floating gene expression into a stable SSN. Then, we calculated the local degree matrix network entropy and degree matrix network entropy according to the information of the differential network. The DMNE method reliably quantifies network fluctuation, i.e., collective fluctuation of molecules caused by a perturbed sample against a group of given reference samples, so as to reduce the noise and thus enhance the robustness and effectiveness by exploring dynamical and high-dimensional information of omics data. Based on the DMNE method, we can detect the critical states before the disease occurs and identify the DNBs. Furthermore, based on the DMNE method, “dark genes” related to disease progression can be identified. Although those “dark genes” are not differentially expressed genes, they are associated with disease progression in T2DM patients and were validated by the functional analysis.
Furthermore, it is noteworthy that the degree matrix network entropy (DMNE) is based on the dynamical network biomarker (DNB), and it can detect the critical states by calculating the score of each stage and using the difference in scores to show the changes in the development of diseases. Therefore, it can be directly applied to stage-course disease datasets in theory, such as stage-course tumor datasets and diabetes datasets, only requiring the modification of corresponding parameters based on sample data and the number of genes when constructing the network. However, it may not be possible to detect the critical states for the disease with fewer than three stages, such as the datasets divided into reference and control only but not staged, because DNB theory for early-warning is not applicable to diseases with fewer than three stages. If a dataset includes fewer than three stages, no valid score comparisons can be made.
In summary, we proposed a robust and effective algorithm that can be used not only to detect critical states of T2DM development and identify DNBs by using stage-course data but also to identify the important “dark genes” and pinpoint drug targets.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/e24091249/s1, Figure S1: Comparison of different parameters of CSN; Figure S2: The identification of “dark genes” for T2DM with PAAD data; Figure S3: Drug target analysis for GSE13268: rat adipose tissue; Figure S4: Drug target analysis for GSE13269: rat gastrocnemius muscle tissue; Figure S5: Drug target analysis for GSE13270: rat liver tissue; Figure S6: The DNBs of GSE13269, involving rat muscle tissue, and GSE13270, involving rat liver tissue, are involved in important biological processes of T2DM. (a) GSE13269. (b) GSE13270; Table S1: The DNBs of critical states of T2DM; Table S2: The degree matrix and local network for GSE13268 (rat adipose tissue); Table S3: The degree matrix and local network for GSE13269 (rat gastrocnemius muscle tissue); Table S4: The degree matrix and local network for GSE13270 (rat liver tissue); Table S5: The function enrichment of DNBs from GSE13268 (rat adipose tissue), GSE13269 (rat gastrocnemius muscle tissue), and GSE13270 (rat liver tissue); Table S6: The “dark genes”, which involve Type 2 diabetes, in rat adipose, gastrocnemius muscle, and liver tissue; Table S7: Enrichment results for DNBs of GSE13269 (rat gastrocnemius muscle tissue) and GSE13270 (rat liver tissue) datasets.

Author Contributions

Conceptualization, Y.Y. and Z.T.; methodology, Y.Y.; software, Y.Y., Z.T., M.S., and Z.G.; validation, Z.T., M.S., and C.M.; formal analysis, Z.T. and Z.G.; investigation, M.S.; resources, P.L.; data curation, M.S. and C.M.; writing—original draft preparation, Y.Y. and Z.T.; writing—review and editing, P.L.; visualization, Y.Y., Z.T., and C.M.; supervision, P.L.; project administration, P.L., Y.Y., and Z.T.; funding acquisition, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (Nos. 61673008, 12071112, and 11471102), the Young Backbone Teacher Funding Scheme of Henan (No. 2019GGJS079), Key R&D and Promotion Special Program of Henan Province (No. 212102310988), basic research projects of key scientific research projects in Henan province (No. 20ZX001), and the Key Science and Technology Research Project of Henan Province of China (No. 222102210053), and the Key Scientific Research Project in Colleges and Universities of Henan Province of China (Grant No. 21A510003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The high-throughput experimental datasets were downloaded from the NCBI GEO database (access ID: GSE13268, GSE13269 and GSE13270) (www.ncbi.nlm.nih. gov/geo/ (accessed on 1 October 2021)). Pancreatic adenocarcinoma (PAAD) survival data and the gene expression values in different samples were downloaded from The Cancer Genome Atlas (TCGA) database (https://cancerge nome.nih.gov/ (accessed on 1 March 2022)). The drug targets were downloaded from DrugBank (https://go.drugbank.com/ (accessed on 20 January 2022)), Therapeutic Target Database (TTD, http://db.idrblab.net/ttd/ (accessed on 20 January 2022)), and Pharmacogenomics Knowledgebase (PharmGKB, https://www.pharmgkb.org/ (accessed on 20 January 2022)). The original codes are available at https://github.com/yyk124/DMNE-main.git (accessed on 26 July 2022).

Acknowledgments

We thank Luochen Chen for his kind guidance, and we thank Lin Li for his valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The GEO datasets used are as follows.
Table A1. The GEO datasets used are as follows.
GEO Datasets List
Access IDTissueRat SpeciesConditionGEO Samples
GSE13268AdiposeWistar–Kyoto (WKY) ratsControlGSM334860-GSM334864
GSM334880-GSM334884
GSM334900-GSM334904
GSM334920-GSM334924
GSM334940-GSM334944
Goto–Kakizaki (GK) ratsDiabeticGSM334850-GSM334854
GSM334870-GSM334874
GSM334890-GSM334894
GSM334910-GSM334914
GSM334930-GSM334934
GSE13269Gastrocnemius muscleWistar–Kyoto (WKY) ratsControlGSM334961-GSM334965
GSM334981-GSM334985
GSM335001-GSM335005
GSM335021-GSM335025
GSM335041-GSM335045
Goto–Kakizaki (GK) ratsDiabeticGSM334951-GSM334955
GSM334971-GSM334975
GSM334991-GSM334995
GSM335011-GSM335015
GSM335031-GSM335035
GSE13270LiverWistar–Kyoto (WKY) ratsControlGSM335062-GSM335066
GSM335082-GSM335086
GSM335102-GSM335106
GSM335122-GSM335126
GSM335142-GSM335146
Goto–Kakizaki (GK) ratsDiabeticGSM335052-GSM335056
GSM335072-GSM335076
GSM335092-GSM335096
GSM335112-GSM335116
GSM335132-GSM335136

References

  1. Liu, R.; Wang, X.; Aihara, K.; Chen, L. Early Diagnosis of Complex Diseases by Molecular Biomarkers, Network Biomarkers, and Dynamical Network Biomarkers. Med. Res. Rev. 2014, 34, 455–478. [Google Scholar] [CrossRef]
  2. Li, M.; Zeng, T.; Liu, R.; Chen, L. Detecting tissue-specific early warning signals for complex diseases based on dynamical network biomarkers: Study of type 2 diabetes by cross-tissue analysis. Brief. Bioinform. 2013, 15, 229–243. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, L.; Liu, R.; Liu, Z.P.; Li, M.; Aihara, K. Detecting early-warning signals for sudden deterioration of complex diseases by dynamical network biomarkers. Sci. Rep. 2012, 2, 342. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, R.; Chen, P.; Chen, L. Single-sample landscape entropy reveals the imminent phase transition during disease progression. Bioinformatics 2020, 36, 1522–1532. [Google Scholar] [CrossRef] [PubMed]
  5. Chen, P.; Li, Y.; Liu, X.; Liu, R.; Chen, L. Detecting the tipping points in a three-state model of complex diseases by temporal differential networks. J. Transl. Med. 2017, 15, 217. [Google Scholar] [CrossRef] [PubMed]
  6. Zhong, J.; Han, C.; Zhang, X.; Chen, P.; Liu, R. Predicting cell fate commitment of embryonic differentiation by single-cell graph entropy. Genomics. Proteom. Bioinform. 2020. [Google Scholar] [CrossRef]
  7. Yan, J.; Li, P.; Gao, R.; Li, Y.; Chen, L. Identifying Critical States of Complex Diseases by Single-Sample Jensen-Shannon Divergence. Front. Oncol. 2021, 11, 371–374. [Google Scholar] [CrossRef]
  8. Gao, R.; Yan, J.; Li, P.; Chen, L. Detecting the critical states during disease development based on temporal network flow entropy. Brief. Bioinform. 2022, bbac164. [Google Scholar] [CrossRef]
  9. Dai, H.; Li, L.; Zeng, T.; Chen, L. Cell-specific network constructed by single-cell RNA sequencing data. Nucleic Acids Res. 2019, 47, e62. [Google Scholar] [CrossRef]
  10. Almon, R.R.; DuBois, D.C.; Lai, W.; Xue, B.; Nie, J.; Jusko, W.J. Gene expression analysis of hepatic roles in cause and development of diabetes in Goto-Kakizaki rats. J. Mater. 2009, 200, 331–346. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, J.; Huang, Q.; Liu, Z.P.; Wang, Y.; Wu, L.Y.; Chen, L.; Zhang, X.S. NOA: A novel network ontology analysis method. Nucleic Acids Res. 2011, 39, e87. [Google Scholar] [CrossRef]
  12. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 10, gkac194. [Google Scholar] [CrossRef] [PubMed]
  13. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed]
  14. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef] [PubMed]
  15. Zhou, Y.; Zhang, Y.; Lian, X.; Li, F.; Wang, C.; Zhu, F.; Qiu, Y.; Chen, Y. Therapeutic target database update 2022: Facil-itating drug discovery with enriched comparative data of targeted agents. Nucleic Acids Res. 2022, 50, D1398–D1407. [Google Scholar] [CrossRef] [PubMed]
  16. Thorn, C.F.; Klein, T.E.; Altman, R.B. PharmGKB: The pharmacogenomics knowledge base. In Phar-Macogenomics; Humana Press: Totowa, NJ, USA, 2013; pp. 311–320. [Google Scholar]
  17. Chen, P. The Algorithm of Detecting the Critical Transition Point of Complex Diseases. Ph.D. Thesis, South China University of Technology, Guangzhou, China, 2018. [Google Scholar]
  18. Eisenberg, E.; Levanon, E.Y. Human housekeeping genes are compact. Trends Genet. 2003, 19, 362–365. [Google Scholar] [CrossRef]
  19. Ma, F.H.; Li, X.J. Islet β Cell self-insulin resistance and type 2 diabetes-a new field of insulin resistance research. Int. J. Endocrinol. Metab. 2004, 24, 154–155. [Google Scholar]
  20. Wang, Z.; Wu, C.W.; Chen, H.; Shi, S.; Zeng, H.; Liu, Y. Dipeptidyl peptidase—IV inhibitors activate ERK1/2 and NF-κB signaling pathway promotes aortic vascular calcification in type 2 diabetic mice. J. Transl. Med. 2019, 8, 5. [Google Scholar]
  21. Qin, C.M.; Gong, C.P.; Wei, X.; Wu, W.H.; Liu, J.; Ou, S.T.; Nephrology, D. Renal artery calcification increases progressive renal damage in type 2 diabetic nephropathy rats. Basic Res. Clin. Med. 2017, 2, 217–223. [Google Scholar]
  22. Jiang, J.; Qin, F.; Liu, J. Analysis of the influence factors for lipid metabolic disturbance of the patients with type 2 diabetes mellitus. Acta Univ. Med. Anhui 2004, 39, 284–286. [Google Scholar]
  23. Wang, L.; Guan, X.; Yu, M.; Qin, W.; Li, L. The relationship between serum fatty acid binding protein 1, protein 3 and type 2 diabetes mellitus complicated with nonalcoholic fatty liver disease. Int. J. Lab. Med. 2021, 42, 3009–3017. [Google Scholar]
  24. He, X. The Association of ND3 Gene Variation of Mitochondrial DNA with T2DM. Master’s Thesis, Wenzhou Medical University, Zhejiang, China, 2007. [Google Scholar]
  25. Ferreira, D.M.; Castro, R.E.; Machado, M.V.; Evangelista, T.; Silvestre, A.; Costa, A.; Coutinho, J.; Carepa, F.; Cortez-Pinto, H.; Rodrigues, C.M.P. Apoptosis and insulin resistance in liver and peripheral tissues of morbidly obese patients is as-so-ciated with different stages of non-alcoholic fatty liver disease. Diabetologia 2011, 54, 1788–1798. [Google Scholar] [CrossRef]
  26. Kelley, D.E.; He, J.; Menshikova, E.V.; Ritov, V.B. Dysfunction of mitochondria in human skeletal muscle in type 2 diabetes. Diabetes 2002, 51, 2944–2950. [Google Scholar] [CrossRef]
  27. Sun, X.H.; Wang, Y.; Mou, Y. Alteration of PI3K/Akt/mTOR signaling during development of diabetic cardiomyopathy and regulation of SIRT1. Chin. Pharmacol. Bull. 2017, 33, 793–798. [Google Scholar]
  28. Liu, J.; Hu, C.; Li, H.; Wang, W.; Zhang, W.; Niu, C.; Chen, M. Based on Kim1/HIF—1 α Signaling pathways explore the mechanism of kidney injury in diabetic rats. Chin. J. Integr. Tradit. West. Nephrol. 2020, 21, 786–788. [Google Scholar]
  29. Yao, J.; Jiang, A.; Wang, L.; Wei, F.; Chen, H.; Yu, H.; Meng, J. Triptolide upregulated PRDX2 inhibits oxidative stress in type 2 diabetic nephropathy. J. Tianjin Med. Univ. 2019, 25, 466–470. [Google Scholar]
  30. Nie, J.; Sage, E.H. SPARC functions as all inhibitor of adipogenesis. J. Cell Commun. Signal. 2009, 3, 247–254. [Google Scholar] [CrossRef]
  31. Kos, K.; Wong, S.; Tan, B.; Gummesson, A.; Jernas, M.; Franck, N.; Kerrigan, D.; Nystrom, F.H.; Carlsson, L.M.S.; Randeva, H.S.; et al. Regulation of the fibrosis and angiogenesis promoter SPARC/osteonectin in human adipose tissue by weight change, leptin, insulin, and glucose. Diabetes 2009, 58, 1780–1788. [Google Scholar] [CrossRef] [Green Version]
  32. Shi, Q.; Bao, S.; Maxwell, J.A.; Reese, E.D.; Friedman, H.S.; Bigner, D.D.; Wang, X.-F.; Rich, J.N. Secreted protein acidic, rich in eysteine (SPARC), mediates cellular survival of gliomas through AKT activation. J. Biol. Chem. 2004, 279, 52200–52209. [Google Scholar] [CrossRef]
  33. Jia, G.; Di, F.; Wang, L.; Yu, H.; Li, X. Association between the serum SPARC and FGF21 levels in newly diagnosed T2DM patients. China Diabetes J. 2013, 21, 223–225. [Google Scholar]
  34. Baron, V.; Muriel, P. Role of glutathione, lipid peroxidation and antioxidants on acute bile-duct obstruction in the rat. Bi-ochim. Biophys. Acta 1999, 1472, 173–180. [Google Scholar] [CrossRef]
  35. Hafner, M.; Rezen, T.; Rozman, D. Regulation of hepatic cytochromes p450 by lipids and cholesterol. Curr. Drug Metab. 2011, 12, 173–185. [Google Scholar] [CrossRef] [PubMed]
  36. She, Y.; Gong, W.; Yu, J. Effect of Vitamin A on Islet Fibrois in Type 2 Diabetes Mellitus Animal Model. Lab. Anim. Sci. Adm. 2015, 32, 17–20. [Google Scholar]
  37. Gross, B.; Pawlak, M.; Lefebvre, P.; Staels, B. PPARs in obesity-induced T2DM, dyslipidaemia and NAFLD. Nature reviews. Endocrinology 2017, 13, 36–49. [Google Scholar] [PubMed]
  38. Shen, Z.; Zhu, L.; Xu, W. Effects of aerobic exercise on AGE-RAGE axis and NF-κB pathway in type 2 diabetic rats. Kunming Med. Univ. 2018, 39, 16–19. [Google Scholar]
  39. Liu, J.; Ye, S.; Bi, S.; Zhou, W. Effect of metformin on the expression of macrophage migration inhibitory factor and its receptor CD74 in renal tissue of type 2 diabetic rats. Acta Univ. Med. Anhui 2021, 56, 5. [Google Scholar]
  40. Jia, G.; Di, F.; Wang, L.; Yu, H.; Li, X. Correlation between serum SPARC and FGF21 levels in patients with new type 2 diabetes mellitus. Chin. J. Diabetes 2013, 21, 5. [Google Scholar]
  41. Johnson, J.M.; Edwards, S.; Shoemaker, D.; Schadt, E.E. Dark Matter in the Genome: Evidence of Widespread Transcription Detected by Microarray Tiling Experiments. Trends Genet. 2005, 21, 93–102. [Google Scholar] [CrossRef]
  42. Ben, Q.; Xu, M.; Ning, X.; Liu, J.; Hong, S.; Huang, W.; Zhang, H.; Li, Z. Diabetes mellitus and risk of pancreatic cancer: A meta-analysis of cohort studies. Eur. J. Cancer 2011, 47, 1928–1937. [Google Scholar] [CrossRef]
  43. Lin, G.; Wan, X.; Liu, D.; Wen, Y.; Yang, C.; Zhao, C. COL1A1 as a potential new biomarker and therapeutic target for type 2 diabetes. Pharmacol. Res. 2021, 165, 105436. [Google Scholar] [CrossRef]
  44. Wang, X.; Lu, J. Effect of blood glucose fluctuation on the prognosis of diabetes and the occurrence and development of chronic complications. Sect. Endocrinol. Foreign M Ed. Sci. 2005, 25, 169–171. [Google Scholar] [CrossRef]
  45. Zhang, W.; Wu, Y.; Fan, W.; Chen, H.; Du, H.; Rao, J. The pattern of plasma BCAA concentration and liver Bckdha gene expression in GK rats during T2D progression. Anim. Models Exp. Med. 2018, 1, 305–313. [Google Scholar] [CrossRef]
  46. Zeng, R.; Xiong, Y.; Zhu, F.; Ma, Z.; Liao, W.; He, Y.; He, J.; Li, W.; Yang, J.; Lu, Q. Fenofibrate attenuated glucose-induced mesangial cells proliferation and extracellular matrix synthesis via PI3K/AKT and ERK1/2. PLoS ONE 2013, 8, e76836. [Google Scholar] [CrossRef] [Green Version]
  47. Choi, J.; Kim, K.-J.; Koh, E.-J.; Lee, B.-Y. Gelidium elegans extract ameliorates type 2 diabetes via regulation of MAPK and PI3K/Akt Signaling. Nutrients 2018, 10, 51. [Google Scholar] [CrossRef]
  48. Hwang, A.; Narayan, V.; Yang, Y.X. Type 2 diabetes mellitus and survival in pancreatic adenocarcinoma: A retrospective cohort study. Cancer 2013, 119, 404–410. [Google Scholar] [CrossRef]
  49. Biadgo, B.; Abebe, M. Type 2 Diabetes Mellitus and Its Association with the Risk of Pancreatic Carcinogenesis: A Review. Korean J Gastroenterol. 2016, 67, 168–177. [Google Scholar] [CrossRef]
  50. Rahn, S.; Zimmermann, V.; Viol, F.; Knaack, H.; Stemmer, K.; Peters, L.; Lenk, L.; Ungefroren, H.; Saur, D.; Schäfer, H.; et al. Diabetes as risk factor for pancreatic cancer: Hyperglycemia promotes epithe-li-al-mesenchymal-transition and stem cell properties in pancreatic ductal epithelial cells. Cancer Lett. 2018, 415, 129–150. [Google Scholar] [CrossRef]
  51. Hoshino, F.; Murakami, C.; Sakai, H.; Satoh, M.; Sakane, F. Creatine kinase muscle type specifically interacts with sa-tu-rated fatty acid-and/or monounsaturated fatty acid-containing phosphatidic acids. Biochem. Biophys. Res. Commun. 2019, 513, 1035–1040. [Google Scholar] [CrossRef]
  52. Fu, M.; Damcott, C.M.; Sabra, M.; Pollin, T.I.; Ott, S.H.; Wang, J.; Garant, M.J.; O’connell, J.R.; Mitchell, B.D.; Shuldiner, A.R. Polymorphism in the Calsequestrin 1 (CASQ1) Gene on Chromosome 1q21 Is Associated with Type 2 Diabetes in the Old Order Amish. Diabetes 2004, 53, 3292–3299. [Google Scholar] [CrossRef]
  53. Briata, P.; Bordo, D.; Puppo, M.; Gorlero, F.; Rossi, M.; Perrone-Bizzozero, N.; Gherzi, R. Diverse roles of the nucleic ac-id-binding protein KHSRP in cell differentiation and disease. Wiley Interdiscip. Rev. RNA 2016, 7, 227–240. [Google Scholar] [CrossRef]
  54. Ruby, M.A.; Riedl, I.; Massart, J.; Åhlin, M.; Zierath, J.R. Protein kinase N2 regulates AMP kinase signaling and insulin responsiveness of glucose metabolism in skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 2017, 313, 483–491. [Google Scholar] [CrossRef] [PubMed]
  55. Anja, B.; Anna-Maria, O.; Jürgen, M.; Mratin, H.; Caroline, K.; Fausto, M.; Fritz, S.; Norbert, S.; Andreas, F.; Hans-Ulrich, H. Common genetic variation in the SERPINF1 locus determines overall adiposity, obesity-related insulin resistance, and circulating leptin levels. PLoS ONE 2012, 7, e34035. [Google Scholar]
  56. Ganic, E.; Johansson, J.K.; Bennet, H.; Fex, M.; Artner, I. Islet-specific monoamine oxidase A and B expression depends on MafA transcriptional activity and is compromised in type 2 diabetes. Biochem. Biophys. Res. Commun. 2015, 468, 629–635. [Google Scholar] [CrossRef] [PubMed]
  57. Tian, J.; Hu, S.; Wang, F.; Yang, X.; Li, Y.; Huang, C. PPARG, AGTR1, CXCL16 and LGALS2 polymor-phisms are cor-re-lated with the risk for coronary heart disease. Int. J. Clin. Exp. Pathol. 2015, 8, 3138–3143. [Google Scholar]
  58. Weiss, E.P.; Brown, M.D.; Shuldiner, A.R.; Hagberg, J.M. Fatty acid binding protein-2 gene variants and insulin resistance: Gene and gene-environment interaction effects. Physiol. Genom. 2002, 10, 145–157. [Google Scholar] [CrossRef]
  59. Andersen, G.; Burgdorf, K.S.; Sparso, T.; Borch-Johnsen, K.; Jorgensen, T.; Hansen, T.; Pedersen, O. AHSG tag single nu-cleotide polymorphisms associate with type 2 diabetes and dyslipidemia: Studies of metabolic traits in 7683 white Danish subjects. Diabetes 2008, 57, 1427–1432. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic illustration of the degree matrix network entropy (DMNE) algorithm. (i) Construct the SSN of reference samples (the healthy control samples) and perturbed samples (the diabetic samples), respectively. (ii) Extract degree matrix and local network. (iii) Calculate the node probability for every gene by fitting a Gaussian distribution. (iv) Calculate the network entropy of reference samples and perturbed samples. (v) Calculate the differential network entropy based on the degree matrix to obtain DMNE score.
Figure 1. Schematic illustration of the degree matrix network entropy (DMNE) algorithm. (i) Construct the SSN of reference samples (the healthy control samples) and perturbed samples (the diabetic samples), respectively. (ii) Extract degree matrix and local network. (iii) Calculate the node probability for every gene by fitting a Gaussian distribution. (iv) Calculate the network entropy of reference samples and perturbed samples. (v) Calculate the differential network entropy based on the degree matrix to obtain DMNE score.
Entropy 24 01249 g001
Figure 2. Identification of the critical states of T2DM. (a) DMNE score curve of GSE13268, for rat adipose tissue, during T2DM, which shows the critical state around 8 weeks. (b) DMNE score curve of GSE13269, for rat gastrocnemius muscle tissue, during T2DM, which shows the critical states around 4 and 16 weeks. (c) DMNE score curve of GSE13270, for rat liver tissue, during T2DM, which shows the critical states around 4 and 16 weeks. (df) The dynamic changes in degree matrix network entropy (DMNE) scores demonstrate the landscape of network entropy. (gi) Dynamic evolution of dynamic network biomarkers (DNBs) of three tissues for T2DM, where the red nodes represent DNBs and the blue nodes represent non-DNBs. In (h), it can be seen that the DNBs are more active and the network structure undergoes significant changes at 4 and 16 weeks. By using the DMNE approach, the early warning signals before the critical transition into disease state can be detected by taking a network perspective.
Figure 2. Identification of the critical states of T2DM. (a) DMNE score curve of GSE13268, for rat adipose tissue, during T2DM, which shows the critical state around 8 weeks. (b) DMNE score curve of GSE13269, for rat gastrocnemius muscle tissue, during T2DM, which shows the critical states around 4 and 16 weeks. (c) DMNE score curve of GSE13270, for rat liver tissue, during T2DM, which shows the critical states around 4 and 16 weeks. (df) The dynamic changes in degree matrix network entropy (DMNE) scores demonstrate the landscape of network entropy. (gi) Dynamic evolution of dynamic network biomarkers (DNBs) of three tissues for T2DM, where the red nodes represent DNBs and the blue nodes represent non-DNBs. In (h), it can be seen that the DNBs are more active and the network structure undergoes significant changes at 4 and 16 weeks. By using the DMNE approach, the early warning signals before the critical transition into disease state can be detected by taking a network perspective.
Entropy 24 01249 g002aEntropy 24 01249 g002b
Figure 3. The dynamic network biomarkers (DNBs) of GSE13268, involving rat adipose tissue, are involved in important biological processes of T2DM. (a) All Rattus norvegicus genes as background. (b) The genes that passed the previous filtering steps as background. The left side of the outer ring represents the detected DNB members, and the right side represents detailed biological processes in which these genes are involved. In the inner ring, the color and width of links indicate diverse enrichment pathways and significant levels of gene functions, respectively.
Figure 3. The dynamic network biomarkers (DNBs) of GSE13268, involving rat adipose tissue, are involved in important biological processes of T2DM. (a) All Rattus norvegicus genes as background. (b) The genes that passed the previous filtering steps as background. The left side of the outer ring represents the detected DNB members, and the right side represents detailed biological processes in which these genes are involved. In the inner ring, the color and width of links indicate diverse enrichment pathways and significant levels of gene functions, respectively.
Entropy 24 01249 g003aEntropy 24 01249 g003b
Figure 4. The enrichment pathways of “dark genes” and the number of involved genes.
Figure 4. The enrichment pathways of “dark genes” and the number of involved genes.
Entropy 24 01249 g004
Figure 5. The drug target analysis in GEO study GSE13270 (rat liver tissue) showing the drug targets and corresponding protein and targeted drugs at (a) stage I and (b) stage IV.
Figure 5. The drug target analysis in GEO study GSE13270 (rat liver tissue) showing the drug targets and corresponding protein and targeted drugs at (a) stage I and (b) stage IV.
Entropy 24 01249 g005aEntropy 24 01249 g005b
Table 1. Enrichment results for DNBs based on all genes of the species in GSE13268 (adipose tissue).
Table 1. Enrichment results for DNBs based on all genes of the species in GSE13268 (adipose tissue).
TissueCaseTermp-ValueTerm Name
AdiposeAdipose
8 weeks
GO:00339933.37 × 10−9Response to lipid
GO:00508961.33 × 10−8Response to stimulus
GO:00485832.51 × 10−5Regulation of response to stimulus
GO:00066290.0238Lipid metabolic process
GO:00507960.0238Regulation of insulin secretion
GO:00019320.0049Regulation of protein phosphorylation
GO:00434340.0124Response to peptide hormone
GO:00703720.0159Regulation of erk1 and erk2 cascade
GO:00068740.0235Cellular calcium ion homeostasis
Table 2. Enrichment results for DNBs based on the genes that passed the previous filtering steps in GSE13268 (adipose tissue).
Table 2. Enrichment results for DNBs based on the genes that passed the previous filtering steps in GSE13268 (adipose tissue).
TissueCaseTermp-ValueTerm Name
AdiposeAdipose
8 weeks
GO:00324960.004939298Response to lipopolysaccharide
GO:00199010.008557741Protein kinase binding
GO:00041750.012946747Endopeptidase activity
GO:00551140.013323312Oxidation–reduction process
GO:00328680.016048825Response to insulin
GO:00337000.01966288Phospholipid efflux
GO:00513840.020031683Response to glucocorticoid
GO:00485450.031308927Response to steroid hormone
GO:00069540.031444423Inflammatory response
GO:00082030.035288946Positive regulation of B cell receptor signaling pathway
GO:00550880.03504665Lipid homeostasis
GO:00507290.038857846Positive regulation of inflammatory response
GO:00082030.047115437Cholesterol metabolic process
Table 3. Enrichment results for DNBs of GSE13269 (gastrocnemius muscle tissue) and GSE13270 (liver tissue) datasets.
Table 3. Enrichment results for DNBs of GSE13269 (gastrocnemius muscle tissue) and GSE13270 (liver tissue) datasets.
TissueCaseTermp-ValueTerm Name
MuscleMuscle
4 weeks
GO:00069360.001436732Muscle contraction
GO:00060960.005358956Glycolytic process
GO:00165040.011881481Peptidase activator activity
GO:00097490.013277291Response to glucose
GO:00312950.019842374T cell co-stimulation
GO:00713330.03156758Cellular response to glucose stimulus
GO:00425930.048961141Glucose homeostasis
rno001900.019507611PI3K-Akt signaling pathway
Muscle
16 weeks
GO:00525470.019261847Regulation of peptidase activity
GO:00423260.021460151Negative regulation of phosphorylation
GO:00059750.017537912Carbohydrate metabolic process
GO:00060960.000947641Glycolytic process
GO:00312950.02116544T cell co-stimulation
GO:00421760.037110591Regulation of protein catabolic process
rno040660.007510911HIF-1 signaling pathway
rno041510.006275621PI3K-Akt signaling pathway
LiverLiver
4 weeks
GO:00069540.001396509Inflammatory response
GO:00097250.002627081Response to hormone
GO:00164910.002701852Oxidoreductase activity
GO:00038240.013755671Catalytic activity
GO:00705550.014842206Response to interleukin-1
GO:00339930.027356508Response to lipid
rno009800.017258367Metabolism of xenobiotics by cytochrome P450
rno049330.027963125AGE-RAGE signaling pathway in diabetic complications
Liver
16 weeks
GO:00018890.00703339Liver development
GO:00082890.00291932Lipid binding
GO:00507770.005409903Negative regulation of immune response
GO:00356930.014337447NOS2-CD74 complex
GO:00067760.028470377Vitamin A metabolic process
rno033200.003125368PPAR signaling pathway
rno049330.006912547AGE-RAGE signaling pathway in diabetic complications
rno009820.008912531Drug metabolism—cytochrome P450
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Y.; Tian, Z.; Song, M.; Ma, C.; Ge, Z.; Li, P. Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis. Entropy 2022, 24, 1249. https://doi.org/10.3390/e24091249

AMA Style

Yang Y, Tian Z, Song M, Ma C, Ge Z, Li P. Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis. Entropy. 2022; 24(9):1249. https://doi.org/10.3390/e24091249

Chicago/Turabian Style

Yang, Yingke, Zhuanghe Tian, Mengyao Song, Chenxin Ma, Zhenyang Ge, and Peiluan Li. 2022. "Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis" Entropy 24, no. 9: 1249. https://doi.org/10.3390/e24091249

APA Style

Yang, Y., Tian, Z., Song, M., Ma, C., Ge, Z., & Li, P. (2022). Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis. Entropy, 24(9), 1249. https://doi.org/10.3390/e24091249

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop