Developing discriminate model and comparative analysis of differentially expressed genes and pathways for bloodstream samples of diabetes mellitus type 2
© Liu et al.; licensee BioMed Central Ltd. 2014
Published: 16 December 2014
Diabetes mellitus of type 2 (T2D), also known as noninsulin-dependent diabetes mellitus (NIDDM) or adult-onset diabetes, is a common disease. It is estimated that more than 300 million people worldwide suffer from T2D. In this study, we investigated the T2D, pre-diabetic and healthy human (no diabetes) bloodstream samples using genomic, genealogical, and phonemic information. We identified differentially expressed genes and pathways. The study has provided deeper insights into the development of T2D, and provided useful information for further effective prevention and treatment of the disease.
A total of 142 bloodstream samples were collected, including 47 healthy humans, 22 pre-diabetic and 73 T2D patients. Whole genome scale gene expression profiles were obtained using the Agilent Oligo chips that contain over 20,000 human genes. We identified 79 significantly differentially expressed genes that have fold change ≥ 2. We mapped those genes and pinpointed locations of those genes on human chromosomes. Amongst them, 3 genes were not mapped well on the human genome, but the rest of 76 differentially expressed genes were well mapped on the human genome. We found that most abundant differentially expressed genes are on chromosome one, which contains 9 of those genes, followed by chromosome two that contains 7 of the 76 differentially expressed genes. We performed gene ontology (GO) functional analysis of those 79 differentially expressed genes and found that genes involve in the regulation of cell proliferation were among most common pathways related to T2D. The expression of the 79 genes was combined with clinical information that includes age, sex, and race to construct an optimal discriminant model. The overall performance of the model reached 95.1% accuracy, with 91.5% accuracy on identifying healthy humans, 100% accuracy on pre-diabetic patients and 95.9% accuract on T2D patients. The higher performance on identifying pre-diabetic patients was resulted from more significant changes of gene expressions among this particular group of humans, which implicated that patients were having profound genetic changes towards disease development.
Differentially expressed genes were distributed across chromosomes, and are more abundant on chromosomes 1 and 2 than the rest of the human genome. We found that regulation of cell proliferation actually plays an important role in the T2D disease development. The predictive model developed in this study has utilized the 79 significant genes in combination with age, sex, and racial information to distinguish pre-diabetic, T2D, and healthy humans. The study not only has provided deeper understanding of the disease molecular mechanisms but also useful information for pathway analysis and effective drug target identification.
Diabetes mellitus of type 2 (T2D) is among 10 most common diseases. According to World Health Organization (WHO), it is estimated that 347 million people suffer from type 2 diabetes. T2D is usually considered not reversible, but if not controlled well, it will eventually lead to fatal complications. However, earlier diagnosis of diabetes requires preventative screening and regular healthcare monitoring, which are not always provided in many countries. Therefore, lower income countries have higher death rates from diabetes. WHO estimated that diabetes deaths will be doubled by the year of 2030. Unfortunately about half of diabetes patients do not realize that they have the disease . In United States of America, diabetes commonly occurs in all ethnic population, and is the seventh most common cause of death. According to American Diabetes Association, the prevalence of the disease is now climbing towards 10% and 30 million people. T2D is also an age related disease. The prevalence in seniors is more than a quarter of the entire population. Because T2D is related to lifestyles, especially foods and exercises, the disease has been and will be consistently increasing and has been estimated to reach one third of the population by 2050 . While type 1 diabetes mellitus is known as juvenile diabetes with strong genetic dispositions or viral involvement, type 2 diabetes is known genetically diverse and contribute to more than 90% of all diabetes . All diabetes can lead to not only long-term but also lethal complications, including cardiovascular, retina, nerve system complication, chronic renal failure, and greater susceptibility of infection. Because T2D is so popular, it has led to the increasing death rate as well as social and economic burdens. While T2D has been known as a genetically complex and multifactorial disease with impaired glucose regulation, such as impaired fasting glucose (IFG) and impaired glucose tolerance (IGT), it is a gradually developing disease and commonly considered as irreversible at our current treatment capabilities. T2D often proceed gradually with aging and can be worsened by many factors, such as hypertension, high cholesterol, lacking of exercise, genetic disposition and family history of diabetes [4–7]. Obesity and immune/inflammatory issues can contribute to the disease [8, 9]. The changes of modern lifestyles in humans have been considered an important factor in the ever-increasing T2D occurrences.
While so far no effective treatment can virtually cure diabetes, significant research progresses have been made in understanding the genetic changes in the development of T2D, especially the mechanisms of gene regulation. New research has led to better prevention of the development of T2D and effective identification of drug targets for blocking or even reversing the disease development . Frayling from the UK's Peninsula College of Medicine and Dentistry found that Single Nucleotide Polymorphism in fat mass and obesity associated gene (FTO) has a strong association with the risk of T2D . Zhao's research showed that analysis of combined gene expression and lipid profiles helped to identify the pathogenesis of T2D . As our research has been unfolded to the study of differential expression of genes and regulatory mechanisms of the genes in diabetes mellitus, we have made significant progresses in finding specific genes and the specific pathways to be targeted by drugs for the purposes of preventing and inversing the disease development.
Pre-diabetic or early stage T2D patients usually do not have any noticeable symptoms and are not diagnosable without blood analysis. Clinical diagnosis of diabetes includes fasting plasma glucose test , hemoglobin A1c test, glucose tolerance test, and clinical screening during physical checkup, but all rely on laboratory blood analysis . Although blood analysis does not provide any genetic, genomic or pathological information about the disease, such information can be useful in assessing the stage, subtype, prognosis, damage and impact of the disease. New research efforts include differential gene expression profiling and genome-wide association studies (GWAS) [15–19] have been made. Yet, the molecular mechanisms of disease are potentially heterogeneous while the limited availability of samples for genome-wide association studies almost prevented effective population genetics analysis and sub-type identification of the disease. In this paper, we assessed and analyzed a number of patients and healthy humans using genomic information including peripheral white blood cell gene expression profile (GEP), and phonemic information including age, gender, and race. Our past studies found that phonemic factors have little influence on GEP. Hierarchical clustering and principal component analysis (PCA) showed that GEP were not directly related to the phonemic factors including gender, blood sugar level, age, and race. However, race and gender are not randomly distributed in the clustering analysis, which implicated that they had potential relevancy with GEP . Therefore, we use age, gender, race, and 79 significant genes as parameters to derive the discriminant model. The model was used successfully to classify the samples of different disease stages respectively with high performance and accuracy. This has led to find molecular mechanisms and genetic diversity for identifying sub-types and pathogenesis of type 2 diabetes mellitus.
Materials and methods
Research objectives, and laboratory and clinical data
Based on our previous research, we now aim to identify characteristic genes and pathways in diabetes. The criteria for the diagnosis of T2D were based on the American Diabetes Association (ADA)  guidelines in accordance with the symptoms in diabetes. The diagnostic criteria are positively correlated with body mass index (BMI) and fasting blood glucose level > 126 mg/dl, or 2H blood glucose level > 200 mg/dl in the oral glucose tolerance test. A total of 142 bloodstream samples were collected, including 47 people from a healthy control group, 22 pre-diabetic, and 73 T2D patients. The experiments were carried out and analyzed using comprehensive information that includes age, gender, race, and GEP. Tougaloo College in Mississippi provided data, and the research was approved by the Institutional Review Board of Tougaloo.
Firstly, Total RNA from 8 - 10 ml of peripheral blood white cells was extracted according to the manufacturer's instructions with LeukoLock ™ general RNA system (Anbion Inc, Austen, Texas, USA). Then the content of the RNA was detected and separated by using Nanodrop spectrophotometer and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, USA). All protocols have been carefully assured and all RNAs have been carefully inspected to ensure that RNAs were not degenerated.
All standard protocols and instructions on handling RNA have been carefully followed. A total of 500 ng of RNA was amplified and labeled by Agilent Low RNA Input Fluorescent Linear Amplification. 850 ng of Cy5- (universal control) and Cy3-labeled (sample) cRNA were mixed and dispersed by the Agilent In Situ Hybridization Kit for every two color array. According to the Agilent 60-mer oligo microarray, hybridizations were put forward for 17 hours in rotating hybridization oven first, and then an Agilent Scanner (G2565AA, Agilent Technologies, Wilmington, Delaware, USA) was used to wash and scan them. Finally, quality control analysis was performed and the Agilent Feature Extraction software (v.22.214.171.124) was used to handle the information from the chips and correction of background noises from default parameters.
Microarray data analysis
Gene chip data was analyzed using GeneSpring 10.0 and quality control was conducted with Pearson correlation coefficients between each sample of the experiment and others pair-wisely. Samples showing less than 80% of correlation with other samples were excluded for further analysis. Scanning probe intensity of less than 5 was directly converted into 5. All probe values using the chip in the 5000 percentile were standardized as per-chip (inner) data. Each gene was standardized using the median value of one gene in all of the samples. The probe characteristics were screened by markers while "Occupying" or "Absence" of the symbol could be used to define the Aglient properties. The Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/) online analysis tool was used to analyze chromosomal localization and function of differentially expressed genes.
Analysis of differentially expressed genes
GeneSpring GX 10.0 (Agilent Technologies, Santa Clara, California, USA) software was used for gene expression analysis. To standardize the data using the method of Lowess , each chip used the 50% as the base, and each gene was standardized by a median reference. The standardized sample data was entered into the GeneSpring GX software. Firstly, quality controls were performed to select the data and samples of the inserted microarray data. Secondly, all the samples were divided into 3 groups: the normal group (ND), pre-diabetic group (PD), and T2D group (D). Thirdly, matched pairwise analysis was conducted and groups were compared pair-wisely. Finally, less significant genes were filtered out. The selection threshold was set as: False discovery rate (FDR) = 5%, p < 0. 05, |fold change| ≥ 2.
Fisher's exact test was used as a discriminant method and a discriminant model was built by all parameters in the experimental group. Statistical software package SPSS 16.0 was used. All samples were finally classified by the discriminant model.
Quality control on samples and entities
We identified 79 significantly differentially expressed genes with fold change ≥ 2.0. In the differentially expressed gene analysis between the pre-diabetes, and T2D groups, there are 24 genes with fold change ≥ 2.0, among which 7 genes were up-regulated and 14 genes were down-regulated. In the comparison between pre-diabetes and normal human groups, there are 74 genes with fold change ≥ 2.0, among which 20 genes were up-regulated and 14 genes were down-regulated. A number of genes were expressed differentially among all groups. For example TFEB gene was not only significant in the T2D and pre-diabetes comparison but also in the healthy and pre-diabetes comparison, with 2.76 and 3.88 times in difference in each pairs of comparisons respectively.
Chromosome location of differentially expressed genes
We performed bioinformatics analysis, and found 79 differentially expressed genes with more than 2 fold changes. Locations of 3 of the79 genes are not known and the other 76 differentially expressed genes were distributed across different chromosomes. Chromosome 1 contains 9 of them (11.3%); and chromosome 2 contains 7 of them (8.86%). These two chromosomes are most abundant in significant genes. No differentially expressed genes were found in chromosomes 4, 9, 20, and the sex Y chromosome. It is reasonable that T2D does not have obvious gender difference. 57 genes were located in the long arm of a chromosome, accounting for 72.2% in total; 19 genes were located in the short arm of a chromosome, accounting for 24.1%. Although chromosome 1 is the longest chromosome in a human genome, differentially expressed genes are still more abundant relatively on chromosomes 1 and 2.
Functional classifications of differentially expressed genes
Classification of Biological Process Categories Based on Gene Ontology (GO) analysis.
Regulation of cell proliferation
LAMB1, ALDH1A2, GPNMB, PTN, IL8RB, CXADR, CTTNBP2, PRAME, IL4, IGFBP5, LAMA1, CXADR
CCR3, IL8RB, CXCL14, IL4, IL8RA
CCR3, IL8RB, CXCL14, IL4, IL8RA
Postive regulation of cell proliferation
LAMB1, ALDH1A2, PTN, IL8RB, PRAME, IL4, LAMA1
LAMA1, APOB, CTTNBP2
Localization of cell
ROPN1B, APOB, NR2F1, IL8RB, CTTNBP2, LAMA1
ROPN1B, APOB, NR2F1, IL8RB, CTTNBP2, LAMA1
ROPN1B, APOB, NR2F1, IL8RB, SEMA6A, CTTNBP2, LAMA1
Rho protein signal transduction
ROPN1B, COL1A2, ARHGAP29
AHNAK, ALDH1A2, ZIC2, NEUROG3, LAMA1
LAMB1, ROPN1B, CCR3, FLRT2, GPNMB, BCAN, CXADR, LAMA1, CXADR
LAMB1, ROPN1B, CCR3, FLRT2, GPNMB, BCAN, CXADR, LAMA1, CXADR
CCR3, IL8RB, CXCL14, IL4, IL8RA
Cellular defense response
CCR3, IL8RB, IL4
CCR3, PTN, IL8RB, CXCL14, IL4, IL8RA
Cell surface receptor linked signal transduction
P2RY14, CCR3, COL1A2, SEMA6A, IL8RB, PTN, IL8RA, LAMA1, FFAR2, ROR1, GPR161, HRH4, STC2
ALDH1A2, COL1A2, ZIC2, LAMA1
Negative regulation of cell proliferation
ALDH1A2, GPNMB, CXADR, CTTNBP2, IGFBP5, CXADR
ALDH1A2, SALL1, CTTNBP2, ZIC2
Morphogenesis of an epithelium
ALDH1A2, ZIC2, LAMA1
A discriminant model was built by 82 parameters (V2-V83) from the 142 samples, and V60 was eliminated in the analysis process. Two typical discrimination functions (Function1, Function2) were extracted, among which Function1 explained 74.4% of all variations and Function2 explained the rest 25.6%. Testing results showed a p-value = 0.000 from Function1 through Function2, which meant the discriminate function has the greatest statistical significance. The functional expressions of discriminate functions were analyzed.
Predicted Group Membership
While Grayson et al.  published an article about the study of human peripheral blood gene chips, their research only showed that the difference of gene expression in T2D played an important role in signal transduction of T cell activation, but the number of samples used in their study (only six cases) was limited, and their samples did not include racial information. When Lei Kong et. al.  discussed the significance of seven microRNA in the serum (miR-9, miR-29a, miR-30d, miR34a, miR-124a, miR146a, and miR375) that are related to diabetes, their samples were 56 healthy controls, 18 newly diagnosed T2D patients (n-T2D), and 19 pre-diabetic patients with known susceptibilities (s-NGT). Canonical discriminant analysis results confirmed 70.6% of n-T2D samples (12/17), while the samples of the s-NGT and pre-diabetic could not be distinguished from each other. Rui Wang-Sattler et al  quantified 140 metabolites of fasting serum samples of 4297 and confirmed candidate biomarkers for pre-diabetes using metabolomic methods to identify three metabolites [glycine, lysophosphatidylcholine (LPC) and acetyl] for prediction of IGT and T2D. Wang et al.  studied the samples of 189 T2D and showed that the increasing content of a small group of essential amino acids [leucine (Leu), valine (Val), isoleucine (Ile)] and aromatic amino acids [phenylalanine (Phe), tyrosine (Tyr)] in serum are associated with risk of T2D by five-fold increment. Our study provides complementary insights into the mechanism of T2D and useful information for better prevention and treatment of the disease and effective identification of drug targets.
This study identified 79 significant genes with more than 2-fold changes in differentially expressed genes using bioinformatics approaches. Differentially expressed genes were mainly distributed in chromosomes 1, 2, 3, 5, and 7, with more abundance in chromosomes 1 and 2. According to Gene Ontology and gene functional analysis, genes which belong to the regulation of cell proliferation were very significant and played important roles in the pathogenesis of T2D. Many genes have multiple functions. For instance, insulin receptor is involved in diabetes and also plays a role in cell proliferation and cancer. Diabetes is a disorder of metabolic syndrome, which will also induce cell proliferation changes on some tissues. T2D patients may have compromised cellular function in absorbing bloodstream sugar, it is reasonable to have elevated gene expression relating to cell proliferation pathways. This study discussed feasibility of combined molecular and bioinformatics methods to distinguish normal humans, pre-diabetic, and T2D effectively. We have analyzed 142 blood samples, including the healthy control group of 47 people, 22 pre-diabetic, and 73 T2D patients. By comparing the gene chip spectrum of these samples, T2D biomarkers can be implicated from the 79 significant genes. Discriminant analysis model showed that combination of 79 genes with three phonemic factors could effectively distinguish healthy human, pre-diabetic, and T2D patients. The results showed that 95.1% of the samples were correctly classified, amongst which 100% was acheived in predicting pre-diabetic samples, 95.9% accuracy was achieved in T2D group, and 91.5% accuracy in healthy human group. The research provided a combined molecular and pedigree analytic method that could potentially lead to an effective screening tool for identifying overall health or illness of humans and prediction of the prognosis of the disease development. The results also showed that 79 genes are significant in diabetes, and these 79 differentially expressed genes have revealed deeper molecular mechanisms of the disease. The research has also led to effective pathway and drug target identification, treatment planning and future therapeutic strategies. In addition, since the discriminant analysis method can separate the pre-diabetes group well from the other two groups, it can lead to the development of new diagnostic tool for the earlier detection of the disease.
The research and publication of the research were supported by the Natural Science Foundation Hubei Province of China (2011CDB236 and 2012FFB04903), Open Research Fund Program of the State Key Laboratory of Virology of China (2012002) and United States National Institutes of Health (NIH). MQY was supported by NIH 5P20GM10342913 and ASTA Award #15-B-23, and YD was supported by NIH R21CA164764.
This article has been published as part of BMC Bioinformatics Volume 15 Supplement 17, 2014: Selected articles from the 2014 International Conference on Bioinformatics and Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/15/S17.
- World Health Organization. [http://www.who.int/diabetes/facts/world_figures/en/index.html]
- Centers for Disease Control and Prevention. [http://www.cdc.gov/diabetes/pubs/pdf/ndfs_20]
- From the Centers for Disease Control and Prevention. [http://www.cdc.gov/mmwr/preview/M]
- Letchuman GR, Wan NW, Wan MW, Chandran LR, Tee GH, Jamaiyah H, Isa MR, Zanariah H, Fatanah I, Ahmad FY: Prevalence of diabetes in the Malaysian National Health Morbidity Survey III 2006. Med J Malaysia. 2010, 65: 180-186.PubMedGoogle Scholar
- Thorpe LE, Upadhyay UD, Chamany S, Garg R, Mandel-Ricci J, Kellerman S, Berger DK, Frieden TR, Gwynn C: Prevalence and control of diabetes and impaired fasting glucose in New York City. Diabetes Care. 2009, 32: 57-62. 10.2337/dc08-0727.PubMed CentralView ArticlePubMedGoogle Scholar
- Escobedo J, Buitron LV, Velasco MF, Ramirez JC, Hernandez R, Macchia A, Pellegrini F, Schargrodsky H, Boissonnet C, Champagne BM: High prevalence of diabetes and impaired fasting glucose in urban Latin America: the CARMELA Study. Diabet Med. 2009, 26: 864-871. 10.1111/j.1464-5491.2009.02795.x.View ArticlePubMedGoogle Scholar
- Qian Y, Lin Y, Zhang T, Bai J, Chen F, Zhang Y, Luo S, Shen H: The characteristics of impaired fasting glucose associated with obesity and dyslipidaemia in a Chinese population. BMC Public Health. 2010, 10: 139-10.1186/1471-2458-10-139.PubMed CentralView ArticlePubMedGoogle Scholar
- Herder C, Brunner EJ, Rathmann W, Strassburger K, Tabak AG, Schloot NC, Witte DR: Elevated levels of the anti-inflammatory interleukin-1 receptor antagonist precede the onset of type 2 diabetes: the Whitehall II study. Diabetes Care. 2009, 32: 421-423. 10.2337/dc08-1161.PubMed CentralView ArticlePubMedGoogle Scholar
- Herder C, Illig T, Rathmann W, Martin S, Haastert B, Muller-Scholze S, Holle R, Thorand B, Koenig W, Wichmann HE, Kolb H: Inflammation and type 2 diabetes: results from KORA Augsburg. Gesundheitswesen. 2005, 67 (Suppl 1): S115-S121.View ArticlePubMedGoogle Scholar
- Kong L, Zhu J, Han W, Jiang X, Xu M, Zhao Y, Dong Q, Pang Z, Guan Q, Gao L: Significance of serum microRNAs in pre-diabetes and newly diagnosed type 2 diabetes: a clinical study. Acta Diabetol. 2011, 48: 61-69. 10.1007/s00592-010-0226-0.View ArticlePubMedGoogle Scholar
- Frayling TM, Hattersley AT, McCarthy A, Holly J, Mitchell SM, Gloyn AL, Owen K, Davies D, Smith GD, Ben-Shlomo Y: A putative functional polymorphism in the IGF-I gene: association studies with type 2 diabetes, adult height, glucose tolerance, and fetal growth in U.K. populations. Diabetes. 2002, 51: 2313-2316. 10.2337/diabetes.51.7.2313.View ArticlePubMedGoogle Scholar
- Singh K, Agrawal NK, Gupta SK, Singh K: Association of Variant rs7903146 (C/T) Single Nucleotide Polymorphism of TCF7L2 Gene With Impairment in Wound Healing Among North Indian Type 2 Diabetes Population: A Case-Control Study. Int J Low Extrem Wounds. 2013Google Scholar
- Reichelt AJ, Spichler ER, Branchtein L, Nucci LB, Franco LJ, Schmidt MI: Fasting plasma glucose is a useful test for the detection of gestational diabetes. Brazilian Study of Gestational Diabetes (EBDG) Working Group. Diabetes Care. 1998, 21: 1246-1249. 10.2337/diacare.21.8.1246.View ArticlePubMedGoogle Scholar
- Ealovega MW, Tabaei BP, Brandle M, Burke R, Herman WH: Opportunistic screening for diabetes in routine clinical practice. Diabetes Care. 2004, 27: 9-12. 10.2337/diacare.27.1.9.View ArticlePubMedGoogle Scholar
- McCarthy MI, Zeggini E: Genome-wide association studies in type 2 diabetes. Curr Diab Rep. 2009, 9: 164-171. 10.1007/s11892-009-0027-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer CN: Novel insights into the etiology of diabetes from genome-wide association studies. Diabetes. 2009, 58: 2444-2447. 10.2337/db09-1153.PubMed CentralView ArticlePubMedGoogle Scholar
- Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007, 445: 881-885. 10.1038/nature05616.View ArticlePubMedGoogle Scholar
- Frayling TM: Genome-wide association studies provide new insights into type 2 diabetes aetiology. Nat Rev Genet. 2007, 8: 657-662. 10.1038/nrg2178.View ArticlePubMedGoogle Scholar
- Zeggini E: A new era for Type 2 diabetes genetics. Diabet Med. 2007, 24: 1181-1186. 10.1111/j.1464-5491.2007.02274.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao C, Mao J, Ai J, Shenwu M, Shi T, Zhang D, Wang X, Wang Y, Deng Y: Integrated lipidomics and transcriptomic analysis of peripheral blood reveals significantly enriched pathways in type 2 diabetes mellitus. BMC Med Genomics. 2013, 6 (Suppl 1): S12-PubMed CentralPubMedGoogle Scholar
- Sakoda H, Ogihara T, Anai M, Funaki M, Inukai K, Katagiri H, Fukushima Y, Onishi Y, Ono H, Yazaki Y: No correlation of plasma cell 1 overexpression with insulin resistance in diabetic rats and 3T3-L1 adipocytes. Diabetes. 1999, 48: 1365-1371. 10.2337/diabetes.48.7.1365.View ArticlePubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.PubMed CentralView ArticlePubMedGoogle Scholar
- Grayson BL, Wang L, Aune TM: Peripheral blood gene expression profiles in metabolic syndrome, coronary artery disease and type 2 diabetes. Genes Immun. 2011, 12: 341-351. 10.1038/gene.2011.13.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang-Sattler R, Yu Z, Herder C, Messias AC, Floegel A, He Y, Heim K, Campillos M, Holzapfel C, Thorand B: Novel biomarkers for pre-diabetes identified by metabolomics. Mol Syst Biol. 2012, 8: 615-PubMed CentralView ArticlePubMedGoogle Scholar
- Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, Lewis GD, Fox CS, Jacques PF, Fernandez C: Metabolite profiles and the risk of developing diabetes. Nat Med. 2011, 17: 448-453. 10.1038/nm.2307.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.