Skip to main content

Combined analysis of whole-exome sequencing and RNA sequencing in type 2 diabetes mellitus patients with thirst and fatigue



The principal objective of this study was to gain a better understanding of the mechanisms of type 2 diabetes mellitus (T2DM) patients with fatigue (D-T2DM) through exome and transcriptome sequencing.


After whole-exome sequencing on peripheral blood of 6 D-T2DM patients, the consensus mutations were screen out and analyzed by a series of bioinformatics analyses. Then, we combined whole-exome sequencing and transcriptome sequencing results to find the important genes that changed at both the DNA and RNA levels.


The results showed that a total of 265,393 mutation sites were found in D-T2DM patients compared with normal individuals, 235 of which were consensus mutations shared with D-T2DM patients. These genes significantly enriched in HIF-1 signaling pathway and sphingolipid signaling pathway. At the RNA level, a total of 375 genes were identified to be differentially expressed. After the DNA-RNA joint analysis, eight genes were screened that changed at both DNA and RNA levels. Among these genes, FUS and LMNA were related to carbohydrate metabolism, energy metabolism, and mitochondrial function. Subsequently, we predicted the herbs, including Qin Pi and Hei Zhi Ma, that might play a therapeutic role in D-T2DM through the SymMap database.


These findings have significant implications for understanding the mechanisms of D-T2DM and provide potential targets for D-T2DM diagnosis and treatment.


Diabetes is a severe worldwide problem threatening the health of millions of people. The number of individuals with type 2 diabetes mellitus (T2DM) is rapidly increasing worldwide. Current estimates indicate that by 2040, approximately 642 million people worldwide will be suffered from T2DM [1]. The occurrence of T2DM has been associated with genetic and acquired factors [2]. In Traditional Chinese Medicine (TCM), there are several different types of diabetes. One type of diabetes is called “Dual Deficiency of Qi and Yin Syndrome”, a common type of diabetes [3]. The main symptoms of this type are dry throat and mouth, fatigued spirit, and lack of strength. In modern medicine, some T2DM patients are also accompanied by thirst and frailty [4,5,6,7,8]. Hence, the diagnosis and treatment of this type of T2DM (D-T2DM) will become a problem facing the world.

Nowadays, with the emergence of next-generation sequencing technology, we can overcome the limitations of traditional genetic disease research methods. Whole-exome sequencing technology (WES) is a method that performs high-throughput analysis of all exon regions to sensitively identify rare and low-frequency disease-related mutations [9, 10]. WES involves three main steps to find the pathogenic genes: capture and enrichment of exome, high-throughput sequencing, and bioinformatic analysis. Several previous studies used WES to uncover variation associated with complex human traits. Some trials have explored the association between T2DM and exonic variants through WES, such as PAX4, RREB1, and PPP1R3A [11,12,13]. So far, however, there has been no detailed investigation of the changes in the exome of D-T2DM patients.

Therefore, this present study using the WES to find the specific mutation sites in D-T2DM patients compared with normal people. In addition, we performed functional enrichment and pathway analysis to explore the potential functions of the mutation gene. Finally, we combined exome and RNA sequencing to explore the pathogenesis of D-T2DM. The results may provide references for clinical application and introduce a new strategy for the treatment of D-T2DM.


Sample information description

Ethical approval accorded with the Declaration of Helsinki for this study was obtained from the Beijing University of Chinese Medicine ethics committee and Beijing Hepingli Hospital ethics committee. The primary inclusion and exclusion criteria are listed in Table 1. Finally, a total of 12 individuals, including six D-T2DM patients and six healthy subjects, were included in this study and provided informed consent. These individuals were separated into the D-T2DM group (patient ID: QYD001, QYD002, QYD003, QYD004, QYD005, and QYD006) and the control group (patient ID: WZC001, WZC002, WZC003, WZC004, WZC005, and WZC006). Then, fasting peripheral blood was collected from the participants and stored at – 80 ℃ for subsequent analysis.

Table 1 Inclusion and exclusion criteria

DNA extraction and sequencing

After DNA was extracted, Agarose gel electrophoresis was used to check the quality of DNA, and Qubit 3.0 was used to determine the concentration of DNA. Next, in brief, DNA samples were fragmented by the Covaris instrument, then the sequencing library was constructed by end repair, A-tailing, and adapter ligation. After that, exome DNA was captured by Agilent SureSelect Human All Exon V6 Kit (Agilent). Then, the DNA libraries were attained after purification and PCR amplification. Finally, pair-end 150 bp sequencing was done using the Illumina sequencing platform.

Bioinformatics analysis

Raw DNA sequencing data were aligned to the GRCh37/hg19, used as the reference genome for sequence alignment and subsequent analysis. The bioinformatics analysis involved quality assessment of sequencing data, mutation detection, mutation screening, and disease-related prediction.

Screening of mutation site

The mutation site screening was performed to analyze the single nucleotide polymorphisms (SNP) or insertion and deletion (InDel) information. Firstly, remove the mutation site with the reported frequency of > 1% by comparison with 1000h_all, esp6500si_all, gnomAD_ALL, and gnomAD_EAS database. Secondly, preserve the mutation of coding regions and spice sites. Thirdly, remove synonymous mutation and non-frameshift InDel mutation. Finally, retain the mutations that were harmful variations or cause splicing defects according to the SIFT, Polyphen, MutationTaster, and CADD software.

Mutation site harmful classification

According to the American Society of Medical Genetics and Genomics (ACMG) variant interpretation guidlines [14], the pathogenicity of mutations are classified into the categories as follows: pathogenic, likely pathogenic, uncertain significance, like benign and benign. Next, the number of mutation sites in each category was counted.

Analysis of the harmfulness of copy number variations (CNV)

CNVs are found in the biological genome, and some malignant CNVs can cause various diseases, including neurological disorders and cancer. CNVs were detected and annotated by CoNIFER and ANNOVAR software, respectively. After that, we further annotated the CNV results database using multiple databases, including DGV and CNVD database, to filter out benign CNVs and reserve malignant CNVs. Then, CNVs were classified into four types, including H (high), P (possibly deleterious), M (medium), and L (low).

Advanced analysis and DNA-RNA conjoint analyses

Further advanced analysis was conducted to identify the true mutations responsible for the disease from all mutation results. The biological function and involving pathway of mutation genes were enriched by Gene Ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. In addition, we collected disease-associated genes from the DisGeNet database and constructed the gene-disease correlation map. Then, protein interaction networks of the mutant genes were constructed using the gene network prediction tool GeneMania ( Finally, we screened significantly differentially expressed genes among the shared mutation genes and performed functional enrichment analysis with GO and KEGG pathway categories.

Herb-gene interaction network

To predict herbs associated with the DE shared mutations, an herb-gene interaction network was constructed base on the SymMap database ( Predicted candidate herbs for each gene were presented and visualized with Cytoscape 3.8.2.

Statistical analysis

All statistical analyses were accomplished by Student's independent-samples t-test using the SPSS software (Version 20.0). Results were expressed as mean ± SEM, and with P < 0.05 considered statistically significant.


Participants information description

At the beginning of the study, we recruited six D-T2DM patients (three males and three females) and six normal individuals (six females). All the patients fulfilled diagnostic criteria for D-T2DM. Figure 1 shows the characteristics of all individuals. There were no apparent differences in age and BMI between the two groups (P > 0.05). In addition, FPG levels were higher in the D-T2DM group when compared to the control group (P < 0.01).

Fig. 1
figure 1

Characteristics of study subjects

Whole-exome sequencing data summary

Whole-exome sequencing was performed on fasting peripheral blood of six D-T2DM patients and six normal individuals. Approximately 153.4 Gb of raw data with a 0.1% average error rate were obtained, and about 91.63% of raw reads had Qphred quality scores of > 30 (Additional file 1: Table S1). The clean reads were aligned and sorted to reference genome GRCh37/hg19 using BWA tools and SAMtools. 99.85% of the mapped content and 97.5% of the target fraction were covered with at least 10 × in depth. The average depth of sequencing was 124.96 (Additional file 1: Table S2).

The results of variation detection

After obtaining mapped data, SAMtools was used to detect and filter the single nucleotide variant (SNV) sites. We detect a total of 265,393 exonic SNVs with 123,095 missense SNV, 886 stop-gain, and 94 stop-less. These SNVs lead to changes in amino acids during translation. Next, ANNOVAR software was used to annotate the SNPs. The results of annotation included the location, type, and conservative prediction of SNPs. Insertion and deletion (InDel) in the coding region and splice site can also cause acid changes during translation. Thus, we determined the count of InDel with different types on the coding region and the genome (Additional file 1: Tables S3 and S4).

Screening and classification of the mutation sites

After fundamental analysis, the SNP/InDel information was screened for mutation sites, and finally, 3720 mutation sites were obtained (Table 2). According to ACMG's evidence, the mutation sites were classified into pathogenic, likely pathogenic, VUS, likely benign, and benign categories (Table 3). Lastly, a harmfulness analysis of CNV was performed and found a total of 160 malignant CNVs (Additional file 1: Table S5).

Table 2 The results of mutation sites screening (partial)
Table 3 Harmful classification screening results of each patient

Screening of shared mutated sites between samples

After filtering the harmful mutation sites, the mutated genes shared by more than two patients were screened. At the same time, follow the principle that 90% of normal individuals did not share. We think the more patients have the same mutated gene, the more likely this mutation is to be associated with D-T2DM. The results indicated that there were 235 shared mutated sites in D-T2DM patients, and the top 15 were listed in Table 4. Among the results, the mutation of NT5DC4 was present in all six D-T2DM patients, and in addition, the PARP1 gene, which is closely related to T2DM, was mutated in three D-T2DM patients.

Table 4 The top 15 shared mutated genes and their annotation results

GO and KEGG pathway enrichment analysis of the mutated sites

Since different genes usually cooperate to perform their biological functions, especially in a complex disease like T2DM. To reveal the biological functions of the mutated genes, gene ontology (GO) analysis was performed to classify genes according to their biological functions. GO analysis enriched the mutated genes in three different categories, including biological process (BP), cellular component (CC), and molecular function (MF). In our results, the total number of significantly enriched BP, CC, and MF were 422, 89, and 81 (P < 0.05), respectively. As well, the most enriched term in BP, CC, and MF was cellular process (GO: 0098656), cytoplasm (GO: 0005737), and molecular function (GO: 0003674), respectively (Fig. 2A, B, C).

Fig. 2
figure 2

Scatter plot of the GO and KEGG pathway enrichment analysis of mutated sites. A GO category of BP. B GO category of CC. C GO category of MF. D KEGG pathway

In addition, to further investigate the biological functions of the mutated genes in D-T2DM patients, we performed pathway enrichment analysis based on the KEGG database. The results showed that a total of 6 pathways were enriched significantly (P < 0.05). The scatterplot in Fig. 2D demonstrates the top 10 enriched pathways. Among them, the HIF-1 signaling pathway (hsa04066) and sphingolipid signaling pathway (hsa04071) are related to T2DM and mitochondrial function.

Gene-disease phenotype correlation analysis

Based on the DisGeNte database V5.0 (a database of gene-disease associations with 561,119 gene-disease and 135,588 mutation-disease association records), we studied the associations between T2DM and our sequencing results. The results showed that in our sequencing results, there were 23,230 genes associated with the development of T2DM. Next, we constructed the gene-phenotype-D-T2DM interaction network based on the first 49 genes of the candidate genes (Fig. 3). Then, we subjected the candidate genes to Phenolyzer Software to generate a ranked gene list. As illustrated in Fig. 4, the higher the rank, the higher the correlation between candidate genes and T2DM, and the top two genes with the highest rank were RELA and APP.

Fig. 3
figure 3

Gene-phenotype-disease association network. The circle size for each of the indicated genes represents the strength of the correlation with the disease. Nodes with a green color represent genes related to diabetes according to existing reports or databases; nodes with an orange color represent genes related to green color genes according to various associations

Fig. 4
figure 4

The top 20 genes in the network, which were ranked by association degree

Protein-protein interaction analysis

To further explore potential interactions among the shared mutated genes, we performed protein functional interaction network analysis using the online software GeneMania [15]. It includes protein–protein, protein-DNA-genetic interactions, pathways, reactions, gene-protein expression data, protein domains-phenotypic screening profiles. Afterwards, the protein–protein interaction network was constructed by Cytoscape software (Fig. 5).

Fig. 5
figure 5

Protein–protein interaction network of the shared mutated genes

DNA and RNA joint analysis results

Screening of the differentially expressed consensus mutated genes in D-T2DM patients

We thought that if a gene is a shared mutated gene in D-T2DM patients and significantly differentially express in D-T2DM patients compared with normal people, this gene may be critical for the development of D-T2DM. Therefore, we analyzed the relative expression of the shared mutated genes at the RNA level. Differential gene expression was assessed using the Cuffdiff software. Then, based on the sequencing depth and gene length for the reads count, FPKM was chosen to express the expression values of genes. The results were shown in Fig. 6 and Table 5: eight shared mutate genes were screened with significant differential expression. Among them, the CDH23 gene mutated in three patients, the rest of the mutated genes were shared with two patients. Moreover, CHD23, LMNA, LRTOMT, VSIG2, and XRCC3 were up-regulated, well FUS, GRAMD1C, and SFTPB were down-regulated in D-T2DM patients compared with normal individuals.

Fig. 6
figure 6

Clustering heatmap of differentially expressed consensus mutated genes

Table 5 Mutant genes and differentially expressed genes in D-T2DM patients

GO and KEGG pathway analysis of the differentially expressed consensus mutant genes in D-T2DM patients

In order to identify a potential biological function for the key genes which were shared mutated genes in DNA level and differentially expressed in RNA level, GO and KEGG pathway analysis were proceed on these key genes. As shown in Fig. 7, the most significantly enriched BP, CC, and MF entry was the biological process, cellular component, and molecular function, respectively. The pathways in cancer, endocytosis, and ECM-receptor interaction were significantly enriched in the KEGG analysis.

Fig. 7
figure 7

The histogram of GO and the scatterplot of the KEGG pathway analysis. A Biological process. B Cellular Component. C Molecular Function. D KEGG pathway

Drug prediction based on the interaction network of herbs and the differentially expressed consensus mutant genes

To further predict herbs that might act on the differentially expressed consensus mutant genes, we constructed an interaction network of these genes with herbs through the SymMap database. As shown in Fig. 8, four genes could be predicted with their associated herbs through the database. A total of 26, 39, 52, and 34 herbs were associated with CDH23, FUS, LMNA, and SFTPB, respectively. Among these herbs, 23 were related to two of the four genes, and four were related to three of the four genes. These four herbs were Qin Pi (Cortex Fraxini, bark of Largeleaf Chinese Ash), Ling Ling Xiang (Lysimachiae Foenigraeci Herba), Hei Zhi Ma (Semen Sesami Nigrum, Black Sesame), and Di Er Cao (Herba Hyperici Japonici, all-grass of Japanese St. Johnswort).

Fig. 8
figure 8

Herb-gene interaction network. Green squares represent the differentially expressed consensus mutant genes; Circles represent the predicted herbs; Different colors represent the different number of herbs associated with the node, yellow nodes associated with one herb, orange nodes associated with two herbs, red nodes associated with three herbs


In our work, we subjected 12 individuals (six D-T2DM patients and six normal individuals) to whole-exome sequencing. After comparison with the existing database, the mutation site was screened out. Then we classified the pathogenicity of these mutations, screened the mutant genes shared in patients, and performed a series of bioinformatics analyses like GO, KEGG, and gene-disease phenotypic correlation analysis. Next, we combined alteration of genetic information of DNA level with RNA level to understand the pathogenesis of D-T2DM further and provide new targets for D-T2DM diagnosis and treatment.

After sample preparation, sequencing procedure and analysis, consensus mutations were screened out. There was a total of 235 consensus mutations. Among them, NT5DC4 has mutations in all 6 D-T2DM patients. Although there have been no studies showing the relationship between NT5DC4 and T2DM, we speculated that NT5DC4 might play a role in the pathogenesis of D-T2DM. In addition, PARP1, poly(ADP-ribose) polymerase 1, is related to immune response, inflammation, and infection responses in diabetes [16]. Studies have shown that PARP1 regulates MMP-9 expression to maintain mitochondrial homeostasis through manipulating the binding of NF-kB/AP-1 at the MMP-9 promoter [17]. Nuclear PARP1 rapidly triggers mitochondrial dysfunction, and the inhibition of PARP can protect mitochondria and reduce ROS production via PARP-1-ATF4-MKP-1-MAPK retrograde pathway [16, 18]. Mitochondria are essential organelles that provide energy in the form of adenosine triphosphate (ATP), including energy needed for muscle contraction, nerve impulses and the synthesis of all complex molecules in the body [19]. Thus, mitochondrial dysfunction will result in fatigue. In Traditional Chinese Medicine, the lack of “Qi” causes fatigue, a hallmark symptom of mitochondrial disease [20]. In addition, impaired mitochondrial function plays an important role in the development of diabetes and insulin resistance [21]. Therefore, the mechanism of D-T2DM may be related to mitochondrial dysfunction. In this study, we identified three missense mutations of PARP1, including G1537A in exon 10, C659T in exon 5, and C14G in exon 1. Hence, it could conceivably be hypothesized that the mutation of PARP1 is possibly associated with the mitochondrial dysfunction in D-T2DM patients.

To further reveal the metabolic pathways associated with the consensus mutations of D-T2DM patients, we performed a KEGG pathway enrichment analysis. The results showed that there was significant enrichment of 6 pathways. Based on previous studies, we found that some of these pathways may be associated with the pathogenesis of D-T2DM, including the metabolic pathway, HIF-1 signaling pathway, and sphingolipid signaling pathway. HIF-1 signaling pathway has well-established roles in insulin secretion and glucose homeostasis [22]. Previous studies established that HIF-1 regulates several critical pathways in the adaptive responses of cells to hypoxia [23]. In addition, HIF-1 could influence mitochondrial function by suppressing both the TCA cycle and respiration and controls mitochondrial biogenesis and autophacy [24]. Sphingolipids maintain the structural integrity of cell membrane and regulate multiple critical cellular processes through signal transduction and gene regulation, which are related to the occurrence and development of various diseases such as diabetes, inflammatory bowel disease, asthma, etc [25, 26]. Furthermore, sphingolipids metabolize nearby mitochondria and regulate mitochondrial structure and function [27]. The dysregulated sphingolipid metabolism leads to mitochondrial dysfunction and disease [27]. In our study, the consensus mutation gene in D-T2DM patients is significantly enriched in HIF-1 and sphingolipid signaling pathways. Consequently, we speculated that the consensus mutant gene might influence the development of D-T2DM through functioning in pathways associated with glycan metabolism, energy metabolism, and mitochondrial disease.

In disease-related research, it is indispensable to identify the correlation between candidate genes and human disease. Therefore, by comparing with the database, we constructed a gene-phenotype-D-T2DM correlation network. In the results, a total of 23 230 genes were found. In addition, we ranked these genes base on their relevance with T2DM. The results showed that the top two relevant genes were RELA (RELA proto-oncogene, NF-kB subunit) and APP (amyloid-beta precursor protein). The combination of NFKB1and RELA constitutes the most abundant form of NF-kappa-B, a ubiquitous transcription factor involved in several biological processes. The activation of NF-kappa-B is a crucial event early in the pathogenesis of diabetes [28]. In addition, NF-kappa-B is a physiological regulator of mitochondrial respiration and plays a role in metabolic adaptation in cells [29]. APP, a precursor protein of amyloid-beta (Aβ), is first cleaved by either α-secretase or β-secretase to produce CTFs 83aa (C83) or 99aa (C99) long, respectively. Then, C83 and C99 are cleaved by PS1 and PS2 to produce p3 or Aβ, respectively. A recent study has reviewed that Aβ could affect insulin sensitivity, reduce glucose-dependent insulin secretion and contribute to the onset of diabetes [30]. Furthermore, mitochondria are one of the major targets that Aβ oligomers negatively impact, including the impairment of fast transport and fragmentation [31]. In addition, C99, the intermediate of Aβ, was a driver of mitochondrial dysfunction in Alzheimer's disease and mediated by the loss of sphingolipid homeostasis [32]. Therefore, since RELA and APP were closely related to T2DM and mitochondrial dysfunction, we believe that RELA and APP may serve as a clinical predictor of the diagnosis in patients with D-T2DM.

D-T2DM is a multifactor disease caused by a combination of environmental and genetic factors. In recent years, with the development of high-throughput sequencing technology, significant breakthroughs in the study on the occurrence and development of D-T2DM have been achieved. But a single-omic approach can only reveal the changes at an individual level. Accordingly, to deepen our understanding of the development mechanisms of D-T2DM, our efforts should focus on integrating multi-omics datasets. The sequencing of the genome contains important information regarding the underlying genetic variations in D-T2DM. The transcriptome sequencing using high-throughput methods enables studies on gene expression of D-T2DM patients and links genome to proteome. Consequently, we combined the sequencing of DNA and RNA level to establish the relationship between the information of genomic mutation and the changes of transcriptome expression.

After the combined analysis, we screened out eight genes altered on both DNA and RNA levels. Among these genes, FUS (FUS RNA binding protein) could interact with HSP60, a mitochondrial chaperonin, to promote mitochondrial damage [33]. The previous study has demonstrated that mitochondrial ATP production is impaired in FUS-expressing cells [34]. In our research, the FUS gene showed missense mutation at rs201533156 sites in the exonic region and significantly down-regulated (log2foldchange = − 1.24573) in D-T2DM patients compared with normal individuals. Lamin A/C (LMNA) plays a role in intracellular redox homeostasis. A study by Tom Sieprath et al. demonstrated that persistent LMNA depletion elevates reactive oxygen species (ROS) levels [35]. ROS, as prime modulators of cellular dysfunction contribution to disease pathophysiology, can act in the mitochondrial energy metabolism and the regulation of metabolic/inflammatory diseases such as diabetes [36]. In this study, we found that the LMNA gene mutated at rs200917748 and significantly up-regulated in mRNA expression levels (log2foldchange = 1.04216) in D-T2DM patients compared with normal individuals. Thus, we hypothesized that the FUS and LMNA genes might play a role in the pathogenesis of D-T2DM. Furthermore, the other six screened genes are also differentially expressed genes, and all have common mutations in D-T2DM patients. Although no previous studies have provided information on the relationship between these genes and T2DM or fatigue, we believe that these genes may have a role in the pathogenesis of D-T2DM because they changed at both the DNA and RNA levels.

To guide subsequent drug studies, we further predicted Chinese herbs associated with the key genes through the database. In total, four herbs were more correlated with differentially expressed shared mutations in the interaction network. Among them, Qin Pi is known to possess anti-inflammatory and anti-oxidative stress effects [37, 38]. In the previous researches by Prabakaran et al., esculetin, one of the main active ingredients in Qin Pi, can exert antioxidant and anti-hyperglycemic capacity in diabetic rats [39, 40]. Hei Zhi Ma is rich in nutrients and has the effect of tonifying the liver and kidney and replenishing vital and blood in TCM theory [41]. Many recent studies have shown that Hei Zhi Ma extract has the ability to improve insulin resistance and reduce blood glucose [42, 43]. In addition, according to the Chinese Pharmacopoeia, Ling Ling Xiang and Di Er Cao can be used to treat fatigue [41]. Accordingly, these herbs were likely to be able to exert the therapeutic effect on D-T2DM. Further experiments need to be done to investigate the detailed mechanisms.


Taken together, this study set out to gain a greater understanding of the mechanisms of occurrence and development of D-T2DM. After whole-exome sequencing on peripheral blood of D-T2DM patients, the consensus mutations were screened out and analyzed by a series of bioinformatics analyses. Then, we combined the results of WES and transcriptome sequencing, and eight genes, including FUS and LMNA, were found that changed at both the DNA and RNA levels. Subsequently, we predicted the herbs that might play a therapeutic role in D-T2DM through the database. However, the pathological mechanism of D-T2DM is complex and multifactorial. The peripheral blood selected as the subject in this study can only partially explain the mechanism of D-T2DM. Therefore, in the future, we intend to investigate the muscle, islet and adipose tissues of D-T2DM patients to further elucidate the developmental mechanism of D-T2DM. In summary, our findings have important implications for understanding the mechanisms of D-T2DM and provide potential targets for D-T2DM diagnosis and treatment.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.





The American Society of Medical Genetics and Genomics


Adenosine triphosphate


Biological process


CTFs 83aa


CTFs 99aa


Cellular component


Copy number variations


Type 2 diabetes mellitus patients with fatigue


Gene Ontology


Insertion and deletion


Kyoto encyclopedia of genes and genomes


Molecular function


Single nucleotide polymorphisms


Type 2 diabetes mellitus


Whole-exome sequencing technology


  1. Premilovac D, Gasperini RJ, Sawyer S, et al. A new method for targeted and sustained induction of type 2 diabetes in rodents. Sci Rep. 2017;7(1):14158.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. Cortez-Espinosa N, Garcia-Hernandez MH, Reynaga-Hernandez E, et al. Abnormal expression and function of Dectin-1 receptor in type 2 diabetes mellitus patients with poor glycemic control (HbA1c>8%). Metabolism. 2012;61(11):1538–46.

    CAS  Article  PubMed  Google Scholar 

  3. Guo J, Chen H, Song J, Wang J, Zhao L, Tong X. Syndrome differentiation of diabetes by the traditional chinese medicine according to evidence-based medicine and expert consensus opinion. Evid Based Complement Alternat Med. 2014;2014: 492193.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Kamoi K, Ishibashi M, Yamaji T. Thirst and plasma levels of vasopressin, angiotensin II and atrial natriuretic peptide in patients with non-insulin-dependent diabetes mellitus. Diabetes Res Clin Pract. 1991;11(3):195–202.

    CAS  Article  PubMed  Google Scholar 

  5. Thompson CJ, Davis SN, Butler PC, Charlton JA, Baylis PH. Osmoregulation of thirst and vasopressin secretion in insulin-dependent diabetes mellitus. Clin Sci (Lond). 1988;74(6):599–606.

    CAS  Article  Google Scholar 

  6. Morley JE, Malmstrom TK, Rodriguez-Manas L, Sinclair AJ. Frailty, sarcopenia and diabetes. J Am Med Dir Assoc. 2014;15(12):853–9.

    Article  PubMed  Google Scholar 

  7. Yanase T, Yanagita I, Muta K, Nawata H. Frailty in elderly diabetes patients. Endocr J. 2018;65(1):1–11.

    Article  PubMed  Google Scholar 

  8. Sinclair AJ, Abdelhafiz AH, Rodriguez-Manas L. Frailty and sarcopenia—newly emerging and high impact complications of diabetes. J Diabetes Complications. 2017;31(9):1465–73.

    Article  PubMed  Google Scholar 

  9. Veltman JA, Brunner HG. De novo mutations in human genetic disease. Nat Rev Genet. 2012;13(8):565–75.

    CAS  Article  PubMed  Google Scholar 

  10. Ng SB, Buckingham KJ, Lee C, et al. Exome sequencing identifies the cause of a mendelian disorder. Nat Genet. 2010;42(1):30–5.

    CAS  Article  PubMed  Google Scholar 

  11. Kwak SH, Chae J, Lee S, et al. Nonsynonymous variants in PAX4 and GLP1R are associated with type 2 diabetes in an East Asian population. Diabetes. 2018;67(9):1892–902.

    CAS  Article  PubMed  Google Scholar 

  12. Barbitoff YA, Serebryakova EA, Nasykhova YA, et al. Identification of novel candidate markers of type 2 diabetes and obesity in russia by exome sequencing with a limited sample size. Genes (Basel). 2018.

    Article  Google Scholar 

  13. Sanchez-Pozos K, Ortiz-Lopez MG, Pena-Espinoza BI, et al. Whole-exome sequencing in maya indigenous families: variant in PPP1R3A is associated with type 2 diabetes. Mol Genet Genomics. 2018;293(5):1205–16.

    CAS  Article  PubMed  Google Scholar 

  14. Richards S, Aziz N, Bale S, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Warde-Farley D, Donaldson SL, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(2):214–20.

    CAS  Article  Google Scholar 

  16. Cipriani G, Rapizzi E, Vannacci A, Rizzuto R, Moroni F, Chiarugi A. Nuclear poly(ADP-ribose) polymerase-1 rapidly triggers mitochondrial dysfunction. J Biol Chem. 2005;280(17):17227–34.

    CAS  Article  PubMed  Google Scholar 

  17. Mishra M, Kowluru R. Role of PARP-1 as a novel transcriptional regulator of MMP-9 in diabetic retinopathy. Biochim Biophys Acta. 2017;1863(7):1761–9.

    CAS  Article  PubMed Central  Google Scholar 

  18. Hocsak E, Szabo V, Kalman N, et al. PARP inhibition protects mitochondria and reduces ROS production via PARP-1-ATF4-MKP-1-MAPK retrograde pathway. Free Radic Biol Med. 2017;108:770–84.

    CAS  Article  PubMed  Google Scholar 

  19. Zhang D, Liu Z, Choi C, et al. Mitochondrial dysfunction due to long-chain Acyl-CoA dehydrogenase deficiency causes hepatic steatosis and hepatic insulin resistance. Proc Natl Acad Sci USA. 2007;104(43):17075–80.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Filler K, Lyon D, Bennett J, et al. Association of mitochondrial dysfunction and fatigue: a review of the literature. BBA Clin. 2014;1:12–23.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Maiese K, Chong ZZ, Shang YC. Mechanistic insights into diabetes mellitus and oxidative stress. Curr Med Chem. 2007;14(16):1729–38.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Taniguchi C, Emanuelli B, Kahn C. Critical nodes in signalling pathways: insights into insulin action. Nat Rev Mol Cell Biol. 2006;7(2):85–96.

    CAS  Article  PubMed  Google Scholar 

  23. Wang G, Jiang B, Rue E, Semenza G. Hypoxia-inducible factor 1 is a basic-helix-loop-helix-PAS heterodimer regulated by cellular O2 tension. Proc Natl Acad Sci USA. 1995;92(12):5510–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. Papandreou I, Cairns RA, Fontana L, Lim AL, Denko NC. HIF-1 mediates adaptation to hypoxia by actively downregulating mitochondrial oxygen consumption. Cell Metab. 2006;3(3):187–97.

    CAS  Article  PubMed  Google Scholar 

  25. Patwardhan GA, Liu YY. Sphingolipids and expression regulation of genes in cancer. Prog Lipid Res. 2011;50(1):104–14.

    CAS  Article  PubMed  Google Scholar 

  26. Russo SB, Ross JS, Cowart LA. Sphingolipids in obesity, type 2 diabetes, and metabolic disease. Handb Exp Pharmacol. 2013;216:373–401.

    CAS  Article  Google Scholar 

  27. Patwardhan GA, Beverly LJ, Siskind LJ. Sphingolipids and mitochondrial apoptosis. J Bioenerg Biomembr. 2016;48(2):153–68.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Patel S, Santani D. Role of NF-kappa B in the pathogenesis of diabetes and its associated complications. Pharmacol Rep. 2009;61(4):595–603.

    CAS  Article  PubMed  Google Scholar 

  29. Mauro C, Leow SC, Anso E, et al. NF-kappaB controls energy homeostasis and metabolic adaptation by up-regulating mitochondrial respiration. Nat Cell Biol. 2011;13(10):1272–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Bharadwaj P, Wijesekara N, Liyanapathirana M, et al. The link between type 2 diabetes and neurodegeneration: roles for amyloid-β, amylin, and tau proteins. J Alzheimers Dis. 2017;59(2):421–32.

    CAS  Article  PubMed  Google Scholar 

  31. Rui Y, Zheng JQ. Amyloid beta oligomers elicit mitochondrial transport defects and fragmentation in a time-dependent and pathway-specific manner. Mol Brain. 2016;9(1):79.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. Pera M, Larrea D, Guardia-Laguarta C, et al. Increased localization of APP-C99 in mitochondria-associated ER membranes causes mitochondrial dysfunction in Alzheimer disease. EMBO J. 2017;36(22):3356–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. Deng J, Yang M, Chen Y, et al. FUS interacts with HSP60 to promote mitochondrial damage. PLoS Genet. 2015;11(9): e1005357.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Stoica R, Paillusson S, Gomez-Suaga P, et al. ALS/FTD-associated FUS activates GSK-3beta to disrupt the VAPB-PTPIP51 interaction and ER-mitochondria associations. EMBO Rep. 2016;17(9):1326–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. Sieprath T, Corne T, Nooteboom M, et al. Sustained accumulation of prelamin A and depletion of lamin A/C both cause oxidative stress and mitochondrial dysfunction but induce different cell fates. Nucleus (Austin, Tex). 2015;6(3):236–46.

    CAS  Article  Google Scholar 

  36. Forrester S, Kikuchi D, Hernandes M, Xu Q, Griendling K. Reactive oxygen species in metabolic and inflammatory signaling. Circ Res. 2018;122(6):877–902.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Liu SQ, He L, Peng H, Liu J. Effect of ash bark on matrix metalloproteinase 1, nitric oxide and prostaglandin E2 in rabbits with experimental osteorthritis. Chin J Clin Rehab. 2005;9(6):150–2.

    Google Scholar 

  38. Marinova EM, Yanishlieva NV, Kostova IN. Antioxidative action of the ethanolic extract and some hydroxycoumarins of Fraxinus ornus bark. Food Chem. 1994;51(2):125–32.

    CAS  Article  Google Scholar 

  39. Prabakaran D, Ashokkumar N. Protective effect of esculetin on hyperglycemia-mediated oxidative damage in the hepatic and renal tissues of experimental diabetic rats. Biochimie. 2013;95(2):366–73.

    CAS  Article  PubMed  Google Scholar 

  40. Prabakaran D, Ashokkumar N. Antihyperglycemic effect of esculetin modulated carbohydrate metabolic enzymes activities in streptozotocin induced diabetic rats. J Funct Foods. 2012;4(4):776–83.

    CAS  Article  Google Scholar 

  41. Pharmacopoeia Commission of the PRC. Pharmacopoeia of The People’s Republic of China, vol. 1. Beijing: People’s Medical Publishing House; 2005.

    Google Scholar 

  42. Jahagirdar A, Usharani D, Srinivasan M, Rajasekharan R. Sesaminol diglucoside, a water-soluble lignan from sesame seeds induces brown fat thermogenesis in mice. Biochem Biophys Res Commun. 2018;507:155–60.

    CAS  Article  PubMed  Google Scholar 

  43. Yang Y, Wang J, Zhang Y, Li J, Sun W. Black sesame seeds ethanol extract ameliorates hepatic lipid accumulation, oxidative stress, and insulin resistance in fructose-induced nonalcoholic fatty liver disease. J Agric Food Chem. 2018;66(40):10458–69.

    CAS  Article  PubMed  Google Scholar 

Download references


Not applicable.


This work was supported by Grants from the National Natural Science Foundation of China (NSFC 81774171).

Author information

Authors and Affiliations



GJ and XL designed the experiments; BL and XY wrote the manuscript; TA, YW, ZH, and BL performed the experiments; YW, FT, JZ, and YH analyzed the data. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Xiaokun Liu or Guangjian Jiang.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of Beijing University of Chinese Medicine (BUCM) (2017BZHYLL0105) and conducted in accordance with the principles of the Declaration of Helsinki. Consent has been obtained from each patient or subject after full explanation of the purpose and nature of all procedures used.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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.

summary of the quality of data. Table S2. Comparison rate and coverage. Table S3. number of InDel in different regions of the genome. Table S4. number of different types of InDel on the coding area. Table S5. Annotation of mutation cite.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Lv, B., Yang, X., An, T. et al. Combined analysis of whole-exome sequencing and RNA sequencing in type 2 diabetes mellitus patients with thirst and fatigue. Diabetol Metab Syndr 14, 111 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Diabetes
  • Whole-exome sequencing
  • RNA sequencing
  • Fatigue
  • Mitochondria