ETM Experimental and Therapeutic Medicine 1792-0981 1792-1015 D.A. Spandidos 10.3892/etm.2017.4042 ETM-0-0-4042 Articles Genome-scale transcriptional analysis reveals key genes associated with the development of type II diabetes in mice ZhangYuchi 1 * HanDongwei 1 * YuPengyang 1 HuangQijing 1 GePengling 1 2 Department of Pharmacology, School of Basic Medical Sciences, Heilongjiang University of Chinese Medicine, Harbin, Heilongjiang 150040, P.R. China The Key Laboratory of Myocardial Ischemia, Harbin Medical University of Chinese Ministry of Education, Harbin, Heilongjiang 150086, P.R. China Correspondence to: Professor Pengling Ge, Department of Pharmacology, School of Basic Medical Sciences, Heilongjiang University of Chinese Medicine, 24 Heping Road, Harbin, Heilongjiang 150040, P.R. China, E-mail: penglingge@126.com

Contributed equally

 03 2017 13 01 2017 13 3 1044 1150 22092015 09112016 Copyright © 2017, Spandidos Publications 2017

Diabetes mellitus is one of the primary diseases that pose a threat to human health. The focus of the present study is type II diabetes (T2D), which is caused by obesity and is the most prevalent type of diabetes. However, genome-scale transcriptional analysis of diabetic liver in the development process of T2D is yet to be further elucidated. Microassays were performed on liver tissue samples from three-, six- and nine-week-old db/db mice with diabetes and db/m mice to investigate differentially expressed mRNA. Based on the results of genome-scale transcriptional analysis, five genes were screened in the present study: chromobox 8 (CBX8), de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1 (DDA1), Phosphoinositide-3-kinase regulatory subunit 6 (PIK3R6), WD repeat domain 41 (WDR41) and Glycine Amidinotransferase (GATM). At three weeks of age, no significant differences in levels or ratios of expression were observed. However, at six and nine weeks, expression of CBX8, DDA1, PIK3R6 and WDR41 was significantly upregulated (P<0.05) in the db/db model group compared with the control group, whereas GATM expression was significantly downregulated (P<0.05). These results suggest that T2D-related differential expression of genes becomes more marked with age, which was confirmed via reverse transcription-quantitative polymerase chain reaction. Genome-scale transcriptional analysis in diabetic mice provided a novel insight into the molecular. events associated with the role of mRNAs in T2D development, with specific emphasis upon CBX8, DDA1, PIK3R6, GATM and WDR41. The results of the present study may provide rationale for the investigation of the target genes of these mRNAs in future studies.

 type II diabetes mRNA gene expression microarray Introduction

Modern lifestyle changes, such as an increased prevalence of high-calorie diets and low levels of physical activity have lead to obesity becoming a global pandemic (1). Obesity may lead to insulin resistance (IR), which in turn may lead to diabetes (2). Type II diabetes (T2D), which is caused by chronic IR and loss of functional islet β-cells is now one of the most prevalent types of diabetes (3,4). In humans and animal models, obesity has been shown to be a major causative factor in the development of IR (2), which is typically associated with dysfunctional metabolism and leads to increased intracellular fat content in liver (5). The liver serves a key role in energy metabolism, as does IR. Therefore, analyzing the changes in gene expression profiles of the liver may elucidate the process by which T2D develops.

However, the underlying molecular mechanisms remain to be determined (6). The identification of a suitable model animal is critical for the study of T2D (6). The db/db mouse is one of the best animal models for the study of human T2D as some of its phenotypes, including polyuria, hyperglycemia and glucosuria, are similar to those observed in human patients with T2D. Db/db mice were discovered by Hummel et al in the inbred line of C57BLKS/J (BKS) in 1966, which is caused by mutations in the gene of leptin receptor (Lepr) gene located in mouse chromosome 4 (7). In db/db mice this typically presents as a hypothalamic defect; wherein satiety is not produced due to a lack of the satiety substance (leptin), lack of reaction to anabolic greater than catabolism, and gradually develop severe T2D accompanied by hyperglycemia due to the accumulation of fat (8) Mice and humans with T2D exhibit similar clinical symptoms, such as polydipsia, polyphagia, urorrhagia, obesity, hyperglycemia, hyperglycemia, IR and lipid metabolism abnormality (9). Using the db/db mice model, it may be possible to investigate the underlying molecular mechanisms and elucidate the complex intrinsic metabolic processes associated with T2D.

It is widely accepted that if only a single gene is mutated, the incidence of T2D is low (10). However, environmental factors may lead to the modification of DNA. When mRNA expression reaches a certain level as a result of this, the body may never recover to the normal state (11). DNA microarray technology allows for simple and precise detection of significant changes in gene expression (11); however, searching large quantities of data to identify small changes in the association between biological molecules and disease is challenging. To overcome this problem in the present study, a comprehensive analysis of differences in mRNA expression in the db/db mice liver tissue was performed by integrating genome-scale microarray gene expression analysis with bioinformatics analysis. Several novel potential biomarkers associated with T2D were also investigated, with a view to facilitating the future development of a more in-depth study of the pathogenesis of T2D.

 Materials and methods <sec> <title>Ethics statement

All animal protocols were approved by the Animal Experimental Ethical Committee of Heilongjiang University of Chinese Medicine (Heilongjiang, China).

 Animal models

Male BKS.Cg-m+/+Leprdb/J (db/db) mice (Nanjing Biomedical Research Institute of Nanjing University, Nanjing, China) aged three (body weight, 7.20±0.51 g), six (body weight, 28.97±0.93 g) and nine weeks old (body weight, 40.27±0.63 g) were used. Age-matched male heterozygote mice (db/m), a non-penetrant genotype (Nanjing Biomedical Research Institute of Nanjing University), were used as the control animals (three weeks old, 7.23±0.62 g; six weeks old, 18.06±0.84 g; nine weeks old, 23.05±0.35 g). A total of 30 mice were individually housed at a constant temperature (20±1°C), humidity (50±5%) and light (12 h/day) conditions with a standard pellet diet and water provided ad libitum.

A total of 15 male db/m mice were divided into three groups: Three-, six- and nine-week-old control models (n=5 each). Similarly, a total of 15 male db/db mice were also divided into three groups: Three-, six- and nine-week-old model groups (n=5 each). All mice were weighed and sacrificed via inhaled 99% ether (Tianjin Tianhe Chemical Co., Ltd., Tianjin China). Blood serum and livers were harvested snap frozen in liquid nitrogen and stored at −80°C for further analysis.

 RNA isolation and microarray profiling

Liver tissues were thawed and homogenized on ice in TRIzol (Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA) and total RNA was extracted using the TRIzol kit (cat. no. 15596-026; Invitrogen; Thermo Fisher Scientific, Inc.), according to the manufacturer's protocol. Concentration and purity of RNA were determined using the ASP-3700 spectrophotometer (ACT Gene, Inc., Piscataway, NJ, USA).

Whole-genome expression profiling was performed using Agilent mouse 4×44 K microarrays (Kang Chen Bio-Tech, Inc., Shanghai, China). Briefly, RNA samples were amplified and labeled with the Agilent One-Color RNA Spike-In Kit (cat. no. 5188-5282; Agilent Technologies, Inc., Santa Clara, CA, USA), and cRNA was hybridized to the arrays in the Agilent Hybridization Chamber (Agilent Technologies, Inc.). Hybridization and washing were performed with the Gene Expression Wash Buffer kit (cat. no. 5188-5327; Agilent Technologies, Inc.), and arrays were scanned with the GenePix 4000B microarray scanner (Molecular Devices LLC, Sunnyvale, CA, USA). Image analysis for grid alignment and the expression data was performed with Nimble Scan software (version 2.5; Roche NimbleGen, Inc., Madison, WI, USA). Volcano plots (GeneSpring Software, version 7.2; Aligent Technologies, Inc.) were used to analyze the raw data files.

 Assessment of IR in the development of T2D

Fasting blood-glucose (FBG) was assayed using a blood glucose meter (Roche Diagnostics, Basel, Switzerland). Fasting insulin (FINS) was assayed using a commercially available Quantikine ELISA kit (cat. no. QC107; R&D Systems, Inc., Minneapolis, MN, USA). The fasting insulin sensitivity index (FISI) was calculated using the following formula: FISI=ln [1/(FBG × FINS)].

 Reverse transcription-quantitative polymerase chain reaction (RT-qPCR)

RT-qPCR was performed to verify gene expression microarray data. cDNA was synthesized via reverse transcription of total liver RNA using the AccuPower Rocket Script RT Premix (Bioneer Corporation, Daejeon, Korea), according to the manufacturer's protocol. Primers used for amplification are presented in Table I. Data were normalized against the expression levels of the housekeeping gene β-actin, and relative expression values were calculated using the 2−∆∆Cq method. PCR products were analyzed by agarose gel electrophoresis to ensure specificity.

 Statistical analysis

Statistical data was expressed as the mean ± standard error of the mean. Statistical significance of the differences between model groups and age-matched control groups was analyzed using a Student's t-test and model groups at different time points were analyzed by one-way analysis of variance with Tukey's multiple comparison post hoc tests. Data were analyzed by SPSS 2.0 statistical software (IBM Corp., Armonk, NY, USA). P<0.05 was considered to indicate a statistically significant difference.

 Results <sec> <title>Physiological parameters of mice with T2D of different ages

FINS, FBG and FISI were measured in db/db mice and controls at three, six and nine weeks old to ensure that they had developed T2D. Although FBG and weight were similar in three-week-old db/db mice and controls they were significantly increased at six (FBG, P<0.05; weight, P<0.001) and nine weeks (both P<0.001; Fig. 1A and B), and there was a significant reduction in FISI in db/db mice compared with their respective controls at three, six and nine weeks old (P<0.001; Fig. 1C). FINS in the db/db mice group was significantly greater than in the control group at three, six and nine weeks old (P<0.001; Fig. 1D). These physiological parameters indicated that db/db mice had begun to develop T2D at three, six and nine weeks old.

 Microarray analysis of global gene expression changes in liver tissues from mice at different ages

One of the aims of the present study was to identify patterns of gene expression that may be associated with the development of T2D. Insight into the molecular changes may ultimately guide clinical treatment of patients prior to the establishment of advanced T2D. Genes that were significantly differentially expressed in mice of different ages are shown in Table II. Global gene expression profiling was performed on RNA prepared from the liver tissue harvested from five db/db model mice and five control mice for each age group (three, six and nine weeks). Results were selected according to the raw data data of the global genome microarray data of mice at each age as shown in Fig. 2A-C.

 Validation of microarray data for five differentially expressed mRNA in liver tissues by RT-qPCR

To validate the results of microarray gene expression profiling, five mRNAs with different expression patterns in the development of T2D, as identified by microarray analysis, were selected for expression verification by RT-qPCR. All expression levels were consistent with microarray predictions. mRNA expression levels of chromobox 8 (CBX8), de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1 (DDA1), Phosphoinositide-3-kinase regulatory subunit 6 (PIK3R6) and WD repeat domain 41 (WDR41) by RT-qPCR increased significantly (P<0.05) and Glycine Amidinotransferase (GATM) decreased significantly (P<0.05) in db/db mice compared with controls (Fig. 3A-E). The upregulation and downregulation of these key genes in db/db mice was markedly greater with increased age of mice (Fig. 4A-C).

 Discussion

T2D is the most prevalent metabolic disease affecting various physiological functions and leading to diverse health complications (12). Although the base causes of T2D appear simple, global molecular changes must occur in response to the inability to regulate glucose. Notably, db/db mice, which are deficient in the leptin receptor and thus leptin signaling, eventually develop hyperphagia and obesity regardless of strain background (13), which induces T2D. Such mice develop insulin resistance at an early age, as well as high FBG levels, insulin levels and FISI (14,15). The results of the present study demonstrate that there were significant differences in FISI, FINS, FBG and weight between the control and model groups for mice at six and nine weeks old, with differences becoming more marked with age. This demonstrates that the clinical symptoms of T2D in model rats become more apparent with age, which provides rationale for the use of db/db mice as a reliable model for T2D.

In the present study, a number of differentially expressed mRNAs in liver tissue were detected in db/db mice compared with normal controls. Differential gene data from mice at three, six and nine weeks old demonstrated that there were hundreds of differentially expressed genes. The five genes that were most differentially expressed were identified and screened. RT-qPCR was subsequently performed and analyzed to verify the differential expression of the five mRNAs from the microarray results. Differential expression of CBX8, DDA1, PIK3R6, GATM and WDR41 in liver tissues was confirmed, which may contribute to the occurrence and development of T2D. The present study may provide useful information pertaining to the association between mRNAs and T2D development and provide rationale for further study.

To elucidate the molecular mechanisms responsible for T2D development, different ages of db/db mice, which is a frequently used and reliable animal model, were used to identify the differentially expressed mRNAs in model and normal control liver tissues and to predict the target genes that these mRNAs regulate. This analysis demonstrated that the differentially expressed mRNAs were responsible for regulating cell proliferation, cell apoptosis and the genes responsible for hormone stimulus (16). This indicates that differential expression of these mRNAs may induce changes in liver tissue and lead to T2D development. The key genes associated with these processes include CBX8, DDA1, PIK3R6, GATM and WDR41; however, further studies are required to confirm their functions with T2D.

CBX8, part of the polycomb group repressive complex 1, belongs to the CBX protein family, which is homologous to the Drosophila polycomb protein (17) and regulates the proliferation of diploid human and mouse fibroblasts through direct binding to the INK4A-ARF locus (18). Furthermore, a previous study demonstrated that the ectopic expression of CBX8 leads to cellular immortalization (19). The activation of auxin response factor and mouse double limit 2 homolog (MDM2) and subsequent protein binding weakens the E3 ubiquitin ligase function, inhibiting p53 protein degradation (20–23). Additionally, a recent study demonstrated that human abdominal adipose tissue from obese subjects and adipose tissue of mice fed a high-fat diet exhibited a reduced expression of the E3 ubiquitin ligase and exhibited IR (24). DDA1 is a DDB1-binding protein that links to a negative regulator of E3 ubiquitin ligase (25). Retinoblastoma and p53 proteins cross link to protect cells against abnormal proliferation signals (26–29). The protein kinase B (Akt) pathway is considered to increase cell survival as the interaction of Akt with the phosphorylation of MDM2 in the PI3K-AKT pathway is associated with the survival, proliferation, growth and death of cells (30). It is well known that the PI3K-AKT signaling pathway is associated with the development of T2D (31). RT-qPCR analysis demonstrated an upregulation in mRNA in model liver tissues compared with normal liver tissues. This suggests that CBX8 and DDA1 are associated with the development of T2D.

PIK3R6 is a phosphoinositide 3-kinase (PI3K)γ subunit (32). PI3Ks are classified as class I, II or III based on substrate binding and sequence homology. Class I PI3Ks are subdivided into α, β, γ and δ. PI3K signaling is well known for its association with various functions, including vesicle trafficking, cell metabolism, cell growth and cell survival (33). It has previously been demonstrated that PI3Kγ serves an essential role in the formation of sarcomas induced by a viral G-protein-coupled-receptor encoded by Kaposi's sarcoma herpes virus (34). PI3Kγ also serves important roles in the heart (35). These studies have demonstrated that PI3K may be used as a target for the treatment of certain inflammatory and cardiac disorders (36). PIK3R6 is able to regulate the activity of PI3K, and further regulate the signal transduction of the PI3K-AKT signaling pathway (37). Phosphorylation of PI3K-AKT has been shown to improve glucose uptake and promote glucose transporter type 4 translocation from the cytoplasm to the plasma membrane, which may stimulate cellular uptake of glucose (38). Furthermore, unpublished data by the present authors demonstrated that the expression of PIK3R6 was upregulated, which may disturb the conduction of the PI3K-AKT signaling pathway and indirectly lead to the development of T2D.

GATM is an enzyme required in the rate-limiting step in the regulation of creatine biosynthesis (39). In the present study, activation of the compensatory mechanisms to increase serine, glycine and the levels of creatine may be induced by the differential expression of phosphoserine phosphatase, GATM and GLDC in the skeletal muscle of individuals with diabetes (40). Additionally, a recent study demonstrated that reduced GATM expression may diminish the capacity for phosphocreatine storage, and modify cellular energy storage and adenosine monophosphate-activated protein kinase (AMPK) signaling pathway conduction (41). The PI3K-AKT and AMPK pathways may be potential targets for the regulation of glucose metabolism associated with IR in patients with T2D and obesity (42). Glucose metabolism occurs primarily in the liver (43). In the present study, it was demonstrated that expression of GATM was downregulated in model liver tissues compared with normal liver tissues. This downregulation may affect the AMPK pathway in the liver, inducing glucose metabolism disorders, which may further affect the development of T2D.

WDR41 contains a WD40 domain consisting of six WD40 repeats. WD40 domains are associated with various cellular functions, including signal transduction, vesicular trafficking, cell cycle control, cell apoptosis, chromatin dynamics and DNA damage response (44–46). These functions are prominent features in proteins that mediate diverse protein-protein interactions and coordinate downstream events, including ubiquitination and histone methylation (45,46). Histone methylation is a reversible process catalyzed by specific and general histone methyltransferases and demethylases, which in turn rely on metabolic coenzymes and respond to changes in energy supply and metabolic status (47). The results of the present study demonstrated that the genome-wide changes in T2D were directly associated with the specific histone modifications in the process of histone methylation (48). The WDR41 gene was upregulated in db/db mice at three, six and nine weeks of age, suggesting that WDR41 is associated with histone methylation and serves a function in the development of T2D.

Gene profiling data from the present study may illustrate the contributions of differentially expressed mRNAs to the underlying mechanisms of T2D in the mouse model; however, compared with human models, there are distinct physiological differences. In conclusion, the present study documents the potential molecular mechanisms associated with T2D in model mice. Novel genes associated with T2D have been identified, and future studies should investigate two aspects; to validate expression and functions of these key mRNAs in T2D human and mice, and to investigate how the knockout or knockin effect of these key mRNA expression levels may affect T2D development.

 Acknowledgements

The authors would like to thank Dr. Wei Liu of the Key Laboratory of Myocardial Ischemia (Harbin Medical University) of the Chinese Ministry of Education for guidance and assistance for technical support. The present study was supported by grants from the National Natural Science Foundation of China (grant no. 81273650); the Chinese Ministry of Science and Technology (grant no. 2012ZX09103201-018); the Natural Science Foundation of Heilongjiang province (grant no. LC2011C03); the Harbin Science and Technology Bureau of Heilongjiang Province (grant no. 2011RFLXS024); Heilongjiang University of Chinese Medicine ‘Excellent Creative Talents Support Program’ (grant no. 2012RCD19) and the Key Laboratory of Myocardial Ischemia, Harbin Medical University, Chinese Ministry of Education (grant no. KF201319).

 References HaslamDWJamesWPObesityLancet36611971209200510.1016/S0140-6736(05)67483-116198769 LangeveldMAertsJMGlycosphingolipids and insulin resistanceProg Lipid Res48196205200910.1016/j.plipres.2009.03.00219303901 AiJWangNYangMDuZMZhangYCYangBFDevelopment of Wistar rat model of insulin resistanceWorld J Gastroenterol1136753679200510.3748/wjg.v11.i24.367515968719 StoffersDAThe development of beta-cell mass: Recent progress and potential role of GLP-1Horm Metab Res36811821200410.1055/s-2004-82616815655713 FlückCESlotboomJNuofferJMKreisRBoeschCMullisPENormal hepatic glycogen storage after fasting and feeding in children and adolescents with type 1 diabetesPediatr Diabetes47076200310.1034/j.1399-5448.2003.00015.x14655262 KuangHHanDXieJYanYLiJGePProfiling of differentially expressed microRNAs in premature ovarian failure in an animal modelGynecol Endocrinol305761201410.3109/09513590.2013.85065924188450 HummelKPDickieMMColemanDLDiabetes, a new mutation in the mouseScience15311271128196610.1126/science.153.3740.11275918576 YunKURyuCSLeeJYNohJRLeeCHLeeHSKangJSParkSKKimBHKimSKHepatic metabolism of sulfur amino acids in db/db miceFood Chem Toxicol53180186201310.1016/j.fct.2012.11.04623220616 DavisRCCastellaniLWHosseiniMBen-ZeevOMaoHZWeinsteinMMJungDYJunJYKimJKLusisAJPéterfyMEarly hepatic insulin resistance precedes the onset of diabetes in obese C57BLKS-db/db miceDiabetes5916161625201010.2337/db09-087820393148 PermuttMAWassonJCoxNGenetic epidemiology of diabetesJ Clin Invest11514311439200510.1172/JCI2475815931378 BonnefondAFroguelPVaxillaireMThe emerging genetics of type 2 diabetesTrends Mol Med16407416201010.1016/j.molmed.2010.06.00420728409 PanzerCLauerMSBriekeABlackstoneEHoogwerfBAssociation of fasting plasma glucose with heart rate recovery in healthy adults: A population-based studyDiabetes51803807200210.2337/diabetes.51.3.80311872683 LutzTAWoodsSCOverview of animal models of obesityCurr Protoc Pharmacol Chapter 5: Unit5.61201210.1002/0471141755.ph0561s58 LeeWHamJKwonHCKimYKKimSNAnti-diabetic effect of amorphastilbol through PPARα/γ dual activation in db/db miceBiochem Biophys Res Commun4327379201310.1016/j.bbrc.2013.01.08323376064 PuffRDamesPWeiseMGökeBSeisslerJParhoferKGLechnerAReduced proliferation and a high apoptotic frequency of pancreatic beta cells contribute to genetically-determined diabetes susceptibility of db/db BKS miceHorm Metab Res43306311201110.1055/s-0031-127181721412687 DweepHStichtCKharkarAPandeyPGretzNParallel analysis of mRNA and microRNA microarray profiles to explore functional regulatory patterns in polycystic kidney disease: Using PKD/Mhm rat modelPLoS One8e53780201310.1371/journal.pone.005378023326503 MaertensGNEl Messaoudi-AubertSRacekTStockJKNichollsJRodriguez-NiedenführMGilJPetersGSeveral distinct polycomb complexes regulate and co-localize on the INK4a tumor suppressor locusPLoS One4e6380200910.1371/journal.pone.000638019636380 BrackenAPKleine-KohlbrecherDDietrichNPasiniDGargiuloGBeekmanCTheilgaard-MönchKMinucciSPorseBTMarineJCThe Polycomb group proteins bind throughout the INK4A-ARF locus and are disassociated in senescent cellsGenes Dev21525530200710.1101/gad.41550717344414 DietrichNBrackenAPTrinhESchjerlingCKKosekiHRappsilberJHelinKHansenKHBypass of senescence by the polycomb group protein CBX8 through direct binding to the INK4A-ARF locusEMBO J2616371648200710.1038/sj.emboj.760163217332741 StottFJBatesSJamesMCMcConnellBBStarborgMBrookesSPalmeroIRyanKHaraEVousdenKHPetersGThe alternative product from the human CDKN2A locus, p14(ARF), participates in a regulatory feedback loop with p53 and MDM2EMBO J1750015014199810.1093/emboj/17.17.50019724636 KamijoTWeberJDZambettiGZindyFRousselMFSherrCJFunctional and physical interactions of the ARF tumor suppressor with p53 and Mdm2Proc Natl Acad Sci USA9582928297199810.1073/pnas.95.14.82929653180 ZhangYXiongYYarbroughWGARF promotes MDM2 degradation and stabilizes p53: ARF-INK4a locus deletion impairs both the Rb and p53 tumor suppression pathwaysCell92725734199810.1016/S0092-8674(00)81401-49529249 PomerantzJSchreiber-AgusNLiégeoisNJSilvermanAAllandLChinLPotesJChenKOrlowILeeHWThe Ink4a tumor suppressor gene product, p19Arf, interacts with MDM2 and neutralizes MDM2's inhibition of p53Cell92713723199810.1016/S0092-8674(00)81400-29529248 YangSWangBHumphriesFHoganAEO'SheaDMoynaghPNThe E3 ubiquitin ligase Pellino3 protects against obesity-induced inflammation and insulin resistanceImmunity41973987201410.1016/j.immuni.2014.11.01325526310 PickELauOSTsugeTMenonSTongYDohmaeNPlafkerSMDengXWWeiNMammalian DET1 regulates Cul4A activity and forms stable complexes with E2 ubiquitin-conjugating enzymesMol Cell Biol2747084719200710.1128/MCB.02432-0617452440 BatesSPhillipsACClarkPAStottFPetersGLudwigRLVousdenKHp14ARF links the tumour suppressors RB and p53Nature395124125199810.1038/258679744267 de StanchinaEMcCurrachMEZindyFShiehSYFerbeyreGSamuelsonAVPrivesCRousselMFSherrCJLoweSWE1A signaling to p53 involves the p19(ARF) tumor suppressorGenes Dev1224342442199810.1101/gad.12.15.24349694807 ZindyFEischenCMRandleDHKamijoTClevelandJLSherrCJRousselMFMyc signaling via the ARF tumor suppressor regulates p53-dependent apoptosis and immortalizationGenes Dev1224242433199810.1101/gad.12.15.24249694806 RadfarAUnnikrishnanILeeHWDePinhoRARosenbergNp19(Arf) induces p53-dependent apoptosis during abelson virus-mediated pre-B cell transformationProc Natl Acad Sci USA951319413199199810.1073/pnas.95.22.131949789064 NacerddineKBeaudryJBGinjalaVWestermanBMattiroliFSongJYvan der PoelHPonzOBPritchardCCornelissen-SteijgerPAkt-mediated phosphorylation of Bmi1 modulates its oncogenic potential, E3 ligase activity and DNA damage repair activity in mouse prostate cancerJ Clin Invest12219201932201210.1172/JCI5747722505453 KitagishiYNakanishiAMinamiAAsaiYYasuiMIwaizakoASuzukiMOnoYOguraYMatsudaSCertain diet and lifestyle may contribute to islet β-cells protection in type-2 diabetes via the modulation of cellular PI3K/AKT PathwayOpen Biochem J87482201410.2174/1874091X0140801007425400709 BiethahnKOrinskaZVigoritoEGoyeneche-PatinoDAMirghomizadehFFögerNBulfone-PausSmiRNA-155 controls mast cell activation by regulating the PI3Kγ pathway and anaphylaxis in a mouse modelAllergy69752762201410.1111/all.1240724734904 EngelmanJALuoJCantleyLCThe evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolismNat Rev Genet7606619200610.1038/nrg187916847462 MartinDGalisteoRMolinoloAAWetzkerRHirschEGutkindJSPI3Kγ mediates kaposi's sarcoma-associated herpesvirus vGPCR-induced sarcomagenesisCancer Cell19805813201110.1016/j.ccr.2011.05.00521665152 PrasadSV NagaLaporteSAChamberlainDCaronMGBarakLRockmanHAPhosphoinositide 3-kinase regulates beta2-adrenergic receptor endocytosis by AP-2 recruitment to the receptor/beta-arrestin complexJ Cell Biol158563575200210.1083/jcb.20020211312163475 VadasODboukHAShymanetsAPerisicOBurkeJEAbi SaabWFKhalilBDHarteneckCBresnickARNürnbergBMolecular determinants of PI3Kγ-mediated activation downstream of G-protein-coupled receptors (GPCRs)Proc Natl Acad Sci USA1101886218867201310.1073/pnas.130480111024190998 GeeringBCutillasPRNockGGharbiSIVanhaesebroeckBClass IA phosphoinositide 3-kinases are obligate p85-p110 heterodimersProc Natl Acad Sci USA10478097814200710.1073/pnas.070037310417470792 ZhuSSunFLiWCaoYWangCWangYLiangDZhangRZhangSWangHCaoFApelin stimulates glucose uptake through the PI3K/Akt pathway and improves insulin resistance in 3T3-L1 adipocytesMol Cell Biochem353305313201110.1007/s11010-011-0799-021461612 WalkerJBCreatine: Biosynthesis, regulation, and functionAdv Enzymol Relat Areas Mol Biol501772421979386719 HansenJSZhaoXIrmlerMLiuXHoeneMSchelerMLiYBeckersJHrabĕ de AngelisMHäringHUType 2 diabetes alters metabolic and transcriptional signatures of glucose and amino acid metabolism during exercise and recoveryDiabetologia5818451854201510.1007/s00125-015-3584-x26067360 ChoeCUNabuursCStockebrandMCNeuANunesPMorelliniFSauterKSchillemeitSHermans-BorgmeyerIMarescauBL-arginine: Glycine amidinotransferase deficiency protects from metabolic syndromeHum Mol Genet22110123201310.1093/hmg/ddt32023026748 SharmaBRKimHJRhyuDYCaulerpa lentillifera extract ameliorates insulin resistance and regulates glucose metabolism in C57BL/KsJ-db/db mice via PI3K/AKT signaling pathway in myocytesJ Transl Med1362201510.1186/s12967-015-0412-525889508 MithieuxGGautier-SteinARajasFZitounCContribution of intestine and kidney to glucose fluxes in different nutritional states in ratComp Biochem Physiol B Biochem Mol Biol143195200200610.1016/j.cbpb.2005.11.00716412674 SmithTFDiversity of WD-repeat proteinsSubcell Biochem482030200810.1007/978-0-387-09595-0_318925368 XuCMinJStructure and function of WD40 domain proteinsProtein Cell2202214201110.1007/s13238-011-1018-121468892 StirnimannCUPetsalakiERussellRBMüllerCWWD40 proteins propel cellular networksTrends Biochem Sci35565574201010.1016/j.tibs.2010.04.00320451393 TeperinoRSchoonjansKAuwerxJHistone methyl transferases and demethylases; can they link metabolism and transcription?Cell Metab12321327201010.1016/j.cmet.2010.09.00420889125 JufvasASjödinSLundqvistKAminRVenerAVStrålforsPGlobal differences in specific histone H3 methylation are associated with overweight and type 2 diabetesClin Epigenetics515201310.1186/1868-7083-5-1524004477

FBG, weight, FISI and FINS in the control and diabetic model groups for mice at three, six and nine weeks old (n=5 each). (A) Significant increases in FBG were observed in the model group at six and nine weeks old compared with their respective control groups. (B) Weights of mice in the model group were significantly greater than those in the control group for mice at six and nine weeks old, whereas weight was similar in mice at three weeks old. (C) Significant increases in FISI were observed in model groups for mice at three, six and nine weeks old compared with the control group. (D) Significant increases in FINS were observed in the model group in three, six and nine weeks old compared with the control group. All data are presented as the mean ± standard error of the mean. *P<0.05 vs. control; **P<0.01 vs. control; ***P<0.001 vs. control. FBG, fasting blood-glucose; FISI, fasting insulin sensitivity index; FINS, fasting insulin; model, BKS.Cg-m+/+Leprdb/J mice.

Globally significant genes as filtered by volcano plot between control and model groups at three, six and nine weeks old. Genome-scale differential expression levels of mRNA between the BKS.Cg-m+/+Leprdb/J model (group B) and control group (group A) at (A) three, (B) six and (C) nine weeks old were screened using a volcano plot. Each dot represents one gene. Red dots indicate genes that were significantly differentially expressed in groups A and B, P<0.05.

RT-qPCR confirmation of microarray results. The expression levels of five genes selected from the microarray data were validated via RT-qPCR. (A) Expression of CBX8 increased in model mice compared with control mice at three, six and nine weeks old. Expression of (B) DDA1, (C) PIK3R6 and (D) WDR41 increased in model mice compared with their respective controls at six and nine weeks old. (E) GATM expression decreased in model mice compared with control mice at six and nine weeks old. All data are presented as the mean ± standard error of the mean. *P<0.05, **P<0.01 and ***P<0.001, vs. control groups. RT-qPCR, reverse transcription-quantitative polymerase chain reaction; CBX8, chromobox 8; DDA1, de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1; PIK3R6, phosphoinositide-3-kinase regulatory subunit 6; WDR41, WD repeat domain 41; GATM, glycine amidinotransferase; model, BKS.Cg-m+/+Leprdb/J mice.

Expression levels of (A) DDA1, (B) CBX8 and GATM, and (C) PIK3R6, and WDR41 in BKS.Cg-m+/+Leprdb/J mice model mice at three, six and nine weeks old. Model groups (n=5 each) were analyzed by one-way analysis of variance post-hoc tests. *P<0.05 vs. 3-week old model group; **P<0.01 vs. 3-week old model group; ***P<0.001 vs. 3-week old model group. DDA1, de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1; CBX8, chromobox 8; GATM, glycine amidinotransferase; PIK3R6, phosphoinositide-3-kinase regulatory subunit 6; WDR41, WD repeat domain 41.

Primers used in reverse transcription-quantitative polymerase chain reaction.

Gene Forward primer sequence Reverse primer sequence
CBX8 5′-ATTCGCAAAGGACGCATGGAA-3′ 5′-CCTCGCTTTTTGGGGCCATA-3′
DDA1 5′-GGCTTGCCCGTCTACAACAA-3′ 5′-CCGCAGAAGTATGTTCGTCTTTT-3′
PIK3R6 5′-GACGCCGTGTCATCATTCC-3′ 5′-TGGTGCAGAAGGTATAAGCTCTA-3′
GATM 5′GCTTCCTCCCGAAATTCCTGT-3′ 5′-CCTCTAAAGGGTCCCATTCGT-3′
WDR41 5′-CAGGGCCTAGCCGAGAAAG-3′ 5′-CCAGCAGTTCAGTATAGGGGTT-3′

CBX8, chromobox 8; DDA1, de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1; PIK3R6, phosphoinositide-3-kinase regulatory subunit 6; WDR41, WD repeat domain 41; GATM, glycine amidinotransferase.

Significantly differentially expressed genes in mice at different ages.

3 weeks old 6 weeks old 9 weeks old
Gene Fold-change P-value Fold-change P-value Fold-change P-value
CBX8 2.21 0.0350 2.48 0.0297 2.71 0.0009
DDA1 0.74 0.0477 3.37 0.0085 3.96 0.0009
PIK3R6 0.56 0.5040 2.1 0.0126 2.58 0.0317
GATM −1.02 0.2521 −3.89 0.0125 −4.44 0.0077
WDR41 0.58 0.0713 4.25 0.0024 4.36 0.0011

List of significant genes filtered where fold-change is >2.0 or ≤-2.0 and P<0.05. CBX8, chromobox 8; DDA1, de-etiolated homolog 1 and damage specific DNA binding protein 1 associated 1; PIK3R6, phosphoinositide-3-kinase regulatory subunit 6; WDR41, WD repeat domain 41; GATM, glycine amidinotransferase.