Combined analysis of whole-exome sequencing and RNA sequencing in type 2 diabetes mellitus patients with thirst and fatigue
Diabetology & Metabolic Syndrome volume 14, Article number: 111 (2022)
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 . The occurrence of T2DM has been associated with genetic and acquired factors . 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 . 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.
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.
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 , 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 (http://genemania.org/). 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 (https://www.symmap.org/). Predicted candidate herbs for each gene were presented and visualized with Cytoscape 3.8.2.
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).
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).
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.
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).
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.
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 . 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).
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.
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.
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).
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 . 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 . 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 . Thus, mitochondrial dysfunction will result in fatigue. In Traditional Chinese Medicine, the lack of “Qi” causes fatigue, a hallmark symptom of mitochondrial disease . In addition, impaired mitochondrial function plays an important role in the development of diabetes and insulin resistance . 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 . Previous studies established that HIF-1 regulates several critical pathways in the adaptive responses of cells to hypoxia . In addition, HIF-1 could influence mitochondrial function by suppressing both the TCA cycle and respiration and controls mitochondrial biogenesis and autophacy . 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 . The dysregulated sphingolipid metabolism leads to mitochondrial dysfunction and disease . 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 . In addition, NF-kappa-B is a physiological regulator of mitochondrial respiration and plays a role in metabolic adaptation in cells . 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 . Furthermore, mitochondria are one of the major targets that Aβ oligomers negatively impact, including the impairment of fast transport and fragmentation . 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 . 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 . The previous study has demonstrated that mitochondrial ATP production is impaired in FUS-expressing cells . 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 . 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 . 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 . 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 . 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
Copy number variations
Type 2 diabetes mellitus patients with fatigue
Insertion and deletion
Kyoto encyclopedia of genes and genomes
Single nucleotide polymorphisms
Type 2 diabetes mellitus
Whole-exome sequencing technology
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. https://doi.org/10.1038/s41598-017-14114-4.
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. https://doi.org/10.1016/j.metabol.2012.03.020.
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. https://doi.org/10.1155/2014/492193.
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. https://doi.org/10.1016/s0168-8227(05)80033-6.
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. https://doi.org/10.1042/cs0740599.
Morley JE, Malmstrom TK, Rodriguez-Manas L, Sinclair AJ. Frailty, sarcopenia and diabetes. J Am Med Dir Assoc. 2014;15(12):853–9. https://doi.org/10.1016/j.jamda.2014.10.001.
Yanase T, Yanagita I, Muta K, Nawata H. Frailty in elderly diabetes patients. Endocr J. 2018;65(1):1–11. https://doi.org/10.1507/endocrj.EJ17-0390.
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. https://doi.org/10.1016/j.jdiacomp.2017.05.003.
Veltman JA, Brunner HG. De novo mutations in human genetic disease. Nat Rev Genet. 2012;13(8):565–75. https://doi.org/10.1038/nrg3241.
Ng SB, Buckingham KJ, Lee C, et al. Exome sequencing identifies the cause of a mendelian disorder. Nat Genet. 2010;42(1):30–5. https://doi.org/10.1038/ng.499.
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. https://doi.org/10.2337/db18-0361.
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. https://doi.org/10.3390/genes9080415.
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. https://doi.org/10.1007/s00438-018-1453-2.
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. https://doi.org/10.1038/gim.2015.30.
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. https://doi.org/10.1093/nar/gkq537.
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. https://doi.org/10.1074/jbc.M414526200.
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. https://doi.org/10.1016/j.bbadis.2017.04.024.
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. https://doi.org/10.1016/j.freeradbiomed.2017.04.018.
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. https://doi.org/10.1073/pnas.0707060104.
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. https://doi.org/10.1016/j.bbacli.2014.04.001.
Maiese K, Chong ZZ, Shang YC. Mechanistic insights into diabetes mellitus and oxidative stress. Curr Med Chem. 2007;14(16):1729–38. https://doi.org/10.2174/092986707781058968.
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. https://doi.org/10.1038/nrm1837.
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. https://doi.org/10.1073/pnas.92.12.5510.
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. https://doi.org/10.1016/j.cmet.2006.01.012.
Patwardhan GA, Liu YY. Sphingolipids and expression regulation of genes in cancer. Prog Lipid Res. 2011;50(1):104–14. https://doi.org/10.1016/j.plipres.2010.10.003.
Russo SB, Ross JS, Cowart LA. Sphingolipids in obesity, type 2 diabetes, and metabolic disease. Handb Exp Pharmacol. 2013;216:373–401. https://doi.org/10.1007/978-3-7091-1511-4_19.
Patwardhan GA, Beverly LJ, Siskind LJ. Sphingolipids and mitochondrial apoptosis. J Bioenerg Biomembr. 2016;48(2):153–68. https://doi.org/10.1007/s10863-015-9602-3.
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. https://doi.org/10.1016/s1734-1140(09)70111-2.
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. https://doi.org/10.1038/ncb2324.
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. https://doi.org/10.3233/jad-161192.
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. https://doi.org/10.1186/s13041-016-0261-z.
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. https://doi.org/10.15252/embj.201796797.
Deng J, Yang M, Chen Y, et al. FUS interacts with HSP60 to promote mitochondrial damage. PLoS Genet. 2015;11(9): e1005357. https://doi.org/10.1371/journal.pgen.1005357.
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. https://doi.org/10.15252/embr.201541726.
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. https://doi.org/10.1080/19491034.2015.1050568.
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. https://doi.org/10.1161/circresaha.117.311401.
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.
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.
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. https://doi.org/10.1016/j.biochi.2012.10.008.
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.
Pharmacopoeia Commission of the PRC. Pharmacopoeia of The People’s Republic of China, vol. 1. Beijing: People’s Medical Publishing House; 2005.
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. https://doi.org/10.1016/j.bbrc.2018.10.195.
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. https://doi.org/10.1021/acs.jafc.8b04210.
This work was supported by Grants from the National Natural Science Foundation of China (NSFC 81774171).
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
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). https://doi.org/10.1186/s13098-022-00884-z