animals were maintained at the university of michigan in a pathogen - free environment and cared for following the university of michigan committee on the care and use of animals guidelines . mice were given continuous access to food ( purina 5001 , purina mills llc , st . after a 6 h fast , tail blood was analyzed using a standard glucometer ( one touch profile , lifescan , inc . , , ghb was measured using the helena laboratories test kit , glyco - tek affinity column method ( helena laboratories corp . , plasma samples from six mice per group ( db / db and db/+ ) were pooled , and each pool was analyzed by fast - protein liquid chromatography . hind paw withdrawal time and tail flick time were recorded as measures of thermal latency ( 14 ) . sural and sciatic ncvs were recorded and ienf density in footpads of the mice was measured as previously described ( 8,15 ) . at 24 weeks of age , mice were killed by sodium pentobarbital overdose and tissues were harvested as previously described ( 8) . the left sciatic nerve was dissected and stored in 30 l rnalater ( ambion , an applied biosystems business , austin , tx ) . based isolation protocol ( rneasy mini kit , qiagen inc . , valencia , ca ) , including an on - column dnase digestion following the manufacturer s protocol . rna quantity and quality were measured by microfluid electrophoresis using the rna 6000 pico labchip on a 2100 bioanalyzer ( agilent technologies , santa clara , ca ) . rna isolated from 10 db / db and 9 db/+ mice was used for hybridization . of the total rna , 75 ng from each sample was amplified and biotin labeled using the ovation biotin - rna amplification and labeling system ( nugen technologies inc . rna amplification and hybridization was performed in two batches ( batch 1 : db / db [ n = 6 ] , db/+ [ n = 5 ] ; batch 2 : db / db [ n = 4 ] , db/+ [ n = 4 ] ) by the university of michigan comprehensive cancer center affymetrix and microarray core facility using the genechip mouse genome 430 2.0 array ( affymetrix , santa clara , ca ) . image files were quantified using affymetrix genechip software ( mas5 ) to generate cel files . body weight , blood glucose and ghb levels , plasma cholesterol and triglyceride levels , thermal latency measures , ienf density , and ncvs of the mice from the two groups were compared using a t test . data for the mice corresponding to outlier microarrays ( see below ) were excluded from these analyses , leaving six db / db and seven db/+ animals . the affymetrix cel files were processed using a local copy of genepattern , a bioinformatics platform from the broad institute ( 16 ) , and standard libraries from the bioconductor project ( http://www.bioconductor.org ) and the r project ( http://www.r-project.org ) . six arrays ( two db/+ and four db / db ) identified as outliers were removed from further analysis based on the array quality assessment . array intensities were normalized with the robust multichip average method , and probe set to gene mapping was done using the brainarray custom cdf version 12 ( 17 ) . log2-transformed data from the two batches of microarrays were adjusted for batch effect bias using distance weighted discrimination batch effect correction method ( 18 ) . out of 16,394 genes on the chip , 915 genes with intensity lower than background were removed from further analysis . raw and processed microarray data are available in the national center for biotechnology information gene expression omnibus repository ( http://www.ncbi.nlm.nih.gov/geo , accession number gse27382 ) . intensity - based moderated t test ( 19 ) was used to determine significant differential expression , and false discovery rate < 0.05 was used as a threshold to obtain a list of differentially expressed genes ( degs ) . significantly overrepresented biological categories in the degs were identified using the database for annotation , visualization , and integrated discovery ( david ) ( 12 ) version 6.7 , and relationships among the categories were visualized using the conceptgen gene set relationship mapping tool ( 20 ) . degs were clustered with chinese restaurant clustering ( crc ) ( 21 ) , and functional enrichment analysis of degs from each of the clusters was performed using david . the genomatix ( http://www.genomatix.de ) literature cocitation network analysis tool ( bibliosphere ) was used to identify pathways and regulatory networks among the degs . protein protein interactions ( ppis ) among the degs were identified using the michigan molecular interaction ( mimi ) database ( http://mimi.ncibi.org ) and visualized using cytoscape ( http://www.cytoscape.org ) . genomatix frameworker database was used to search the promoter sequences for overrepresented transcription factor binding site ( tfbs ) motifs . a region of 1,000 base pair ( bp ) upstream and 100 bp downstream of a transcription start site was designated as a promoter region . a distance of 1500 bp and a distance variance of 100 bp were permitted between tfbss . the expression levels of five degs ( highly upregulated in the db / db samples ) were measured by real time rt - pcr ( qrt - pcr ) in five db / db and six db/+ samples for technical validation of microarray data . the amount of rna from each sample was first normalized to an endogenous reference ( tata box binding protein [ tbp ] ; ct ) and then relative to the control group ( ct ) and was expressed as 2 . bks.cg-mlepr/j mice ( stock number 000642 ) were purchased from jackson laboratories ( bar harbor , me ) . animals were maintained at the university of michigan in a pathogen - free environment and cared for following the university of michigan committee on the care and use of animals guidelines . mice were given continuous access to food ( purina 5001 , purina mills llc , st . fasting blood glucose levels were measured every 4 weeks . after a 6 h fast , tail blood was analyzed using a standard glucometer ( one touch profile , lifescan , inc . , , ghb was measured using the helena laboratories test kit , glyco - tek affinity column method ( helena laboratories corp . , plasma samples from six mice per group ( db / db and db/+ ) were pooled , and each pool was analyzed by fast - protein liquid chromatography . hind paw withdrawal time and tail flick time were recorded as measures of thermal latency ( 14 ) . sural and sciatic ncvs were recorded and ienf density in footpads of the mice was measured as previously described ( 8,15 ) . at 24 weeks of age , mice were killed by sodium pentobarbital overdose and tissues were harvested as previously described ( 8) . the left sciatic nerve was dissected and stored in 30 l rnalater ( ambion , an applied biosystems business , austin , tx ) . based isolation protocol ( rneasy mini kit , qiagen inc . , valencia , ca ) , including an on - column dnase digestion following the manufacturer s protocol . rna quantity and quality were measured by microfluid electrophoresis using the rna 6000 pico labchip on a 2100 bioanalyzer ( agilent technologies , santa clara , ca ) . rna isolated from 10 db / db and 9 db/+ mice was used for hybridization . of the total rna , 75 ng from each sample was amplified and biotin labeled using the ovation biotin - rna amplification and labeling system ( nugen technologies inc . rna amplification and hybridization was performed in two batches ( batch 1 : db / db [ n = 6 ] , db/+ [ n = 5 ] ; batch 2 : db / db [ n = 4 ] , db/+ [ n = 4 ] ) by the university of michigan comprehensive cancer center affymetrix and microarray core facility using the genechip mouse genome 430 2.0 array ( affymetrix , santa clara , ca ) . image files were quantified using affymetrix genechip software ( mas5 ) to generate cel files . body weight , blood glucose and ghb levels , plasma cholesterol and triglyceride levels , thermal latency measures , ienf density , and ncvs of the mice from the two groups were compared using a t test . data for the mice corresponding to outlier microarrays ( see below ) were excluded from these analyses , leaving six db / db and seven db/+ animals . the affymetrix cel files were processed using a local copy of genepattern , a bioinformatics platform from the broad institute ( 16 ) , and standard libraries from the bioconductor project ( http://www.bioconductor.org ) and the r project ( http://www.r-project.org ) . six arrays ( two db/+ and four db / db ) identified as outliers were removed from further analysis based on the array quality assessment . array intensities were normalized with the robust multichip average method , and probe set to gene mapping was done using the brainarray custom cdf version 12 ( 17 ) . log2-transformed data from the two batches of microarrays were adjusted for batch effect bias using distance weighted discrimination batch effect correction method ( 18 ) . out of 16,394 genes on the chip , 915 genes with intensity lower than background were removed from further analysis . raw and processed microarray data are available in the national center for biotechnology information gene expression omnibus repository ( http://www.ncbi.nlm.nih.gov/geo , accession number gse27382 ) . intensity - based moderated t test ( 19 ) was used to determine significant differential expression , and false discovery rate < 0.05 was used as a threshold to obtain a list of differentially expressed genes ( degs ) . significantly overrepresented biological categories in the degs were identified using the database for annotation , visualization , and integrated discovery ( david ) ( 12 ) version 6.7 , and relationships among the categories were visualized using the conceptgen gene set relationship mapping tool ( 20 ) . degs were clustered with chinese restaurant clustering ( crc ) ( 21 ) , and functional enrichment analysis of degs from each of the clusters was performed using david . the genomatix ( http://www.genomatix.de ) literature cocitation network analysis tool ( bibliosphere ) was used to identify pathways and regulatory networks among the degs . protein interactions ( ppis ) among the degs were identified using the michigan molecular interaction ( mimi ) database ( http://mimi.ncibi.org ) and visualized using cytoscape ( http://www.cytoscape.org ) . twenty - two degs from one crc cluster , all annotated with  lipid metabolism ,  were selected for promoter sequence analysis . genomatix frameworker database was used to search the promoter sequences for overrepresented transcription factor binding site ( tfbs ) motifs . a region of 1,000 base pair ( bp ) upstream and 100 bp downstream of a transcription start site was designated as a promoter region . a distance of 1500 bp and a distance variance of 100 bp were permitted between tfbss . the expression levels of five degs ( highly upregulated in the db / db samples ) were measured by real time rt - pcr ( qrt - pcr ) in five db / db and six db/+ samples for technical validation of microarray data . the amount of rna from each sample was first normalized to an endogenous reference ( tata box binding protein [ tbp ] ; ct ) and then relative to the control group ( ct ) and was expressed as 2 . table 1 lists phenotyping parameters for the db / db and db/+ mice measured at 24 weeks of age . the db / db mice are significantly heavier and exhibit higher levels of fasting blood glucose and ghb than the db/+ mice at 24 weeks of age . plasma cholesterol and triglyceride levels are significantly elevated in the db / db mice , confirming dyslipidemia in these mice . thermal latency in the hind paw and tail is significantly increased in db / db mice , indicative of sensory loss . both ienf density and sciatic motor nerve conduction velocity ( smncv ) are significantly lower in db / db mice . increased thermal latency , loss of ienf density , and decreased smncv confirm the development of dn in the db / db mice at 24 weeks of age . phenotyping parameters for the db / db and db/+ mice diabetes , dyslipidemia , and neuropathy phenotyping parameters measured at 24 weeks of age ( mean  sd ) . figure 1 outlines the steps in the transcriptional analysis of the microarray data . using intensity - based moderated t test false discovery rate < 0.05 as a threshold , we identified 4,017 degs ; 2,122 genes were upregulated and 1,895 genes were downregulated in the db / db relative to the db/+ samples . microarray data were validated using qrt - pcr for 5 selected degs , which demonstrated high correlation between microarray and qrt - pcr expression levels ( supplementary table 1 ) . rna samples obtained from the sciatic nerve ( scn ) of db / db and db/+ mice were hybridized to affymetrix genechip microarrays . enrichment analysis identifies functional categories that have a higher representation in the regulated genes and , thus , are likely to be active in the experimental condition . table 2 lists gene ontology ( go ) terms overrepresented in degs identified using david . overrepresented biological processes with a higher number of upregulated genes are cell cycle , lipid metabolic process , lipid transport , carbohydrate metabolic process , response to stress , and apoptosis . overrepresented cellular component terms are mitochondrion and axon , with a majority of mitochondrial genes upregulated and a majority of axonal genes downregulated . bh p value , benjamini - hochberg corrected p value of enrichment in david ; go_bp , go biological process ; go_cc , go cellular component . gene network analysis is useful for exploring both published and novel gene associations that may be involved in the process being studied . conceptgen is a gene set relation mapping tool that identifies relationships ( significant overlap ) between an experimental gene set and curated gene sets in sources such as go terms , kyoto encyclopedia of genes and genomes ( kegg ) pathways , and medical subject heading ( mesh ) terms ( 20 ) . conceptgen was used to analyze 717 degs annotated with the go biological processes shown in table 1 . the mouse gene identifiers were converted by conceptgen homology conversion tool to 695 human homologs . relationship among selected concepts ( go biological processes , kegg pathways , and mesh terms ) associated with the gene set were visualized using conceptgen network graph tool ( fig . kegg pathways with significant overlap with the deg set include fatty acid metabolism , peroxisome proliferator  activated receptor ( ppar ) signaling , and extracellular matrix  receptor interaction . concepts related to the deg set clustered in three groups : lipid metabolism , apoptosis , and axonogenesis ( fig . network depicting relationships and overlap among significantly overrepresented concepts ( go biological processes , kegg pathways , and mesh terms ) in the degs . node color represents concept type , node size represents number of genes in the concept , and thickness of the edge represents gene overlap between concepts . bibliosphere generates a network of genes based on their cocitation in pubmed abstracts and identifies functionally related subnetworks based on known regulatory relationships and pathway annotation . figure 3a illustrates a bibliosphere network of degs centered on tumor necrosis factor- ( tnf- ) with gene tnf- was chosen as the central node because of its importance in type 2 diabetes and lipid metabolism ( 22 ) . central nodes of subnetworks and associated pathways include hifna ( apoptosis ) , serpine1 ( p53 signaling ) , app ( axonogenesis ) , lyn ( fc  ri signaling ) , and pparg ( ppar signaling and fatty acid oxidation ) ( fig . nodes represent genes , node color represents fold change ( red / yellow : upregulation ; blue : downregulation ) . the network is centered at tnf-. only shortest paths to the central node are shown for clarity . b : enlarged view of the part of the network identified in a illustrates connections between the central node and other highly connected nodes in the network . highly connected nodes in the network include pparg , app , serpine1 , hif1a , irf1 , vdr , lyn , and ripk1 . ppi network of degs was constructed using cytoscape mimi plug - in ( 23 ) . subnetworks were identified based on pathway annotation of genes in mimi ( fig . regulated pathways identified by ppi network are carbohydrate metabolism ( glycolysis and tricarboxylic acid cycle [ tca ] ) , energy metabolism ( oxidative phosphorylation ) , lipid metabolism ( fatty acid metabolism and glycerolipid metabolism ) , glutathione metabolism , cell development ( cell cycle , apoptosis , axon guidance , and neurotrophin signaling ) , cell communication ( cell adhesion and extracellular matrix  receptor interaction ) , and signal transduction pathways ( ppar signaling , transforming growth factor- [ tgf- ] signaling , adipocytokine signaling , wnt signaling , mitogen - activated protein kinase [ mapk ] signaling , erbb signaling , and mammalian target of rapamycin [ mtor ] signaling ) b : enlarged view of part of the network identified in a. node color represents pathway annotation of genes ( for example , blue : oxidative phoshorylation ; orange : tca ; red : apoptosis ) . supplementary table 2 lists degs in the pathways regulated in the db / db nerve . the majority of degs in the glycolysis , tca , oxidative phosphorylation , fatty acid metabolism , glycerolipid metabolism , mitochondrial fatty acid elongation , lipid transport , adipocytokine signaling , ppar signaling , and apoptosis pathways are upregulated , whereas most of the axonogenesis - related genes are downregulated in db / db nerve . crc , a model - based bayesian clustering tool , grouped the degs into 18 clusters based on positive as well as negative correlation among their expression profiles . we performed enrichment analysis for gene sets from individual clusters to identify overrepresented functional categories ( supplementary table 3 ) . mitochondrion ,   lipid metabolic process ,  and  axonogenesis .  gene regulatory sequence elements are organized into motifs of tfbss known as promoter modules , possibly suggesting coordinated regulation by multiple transcription factors ( tfs ) ( 24 ) . we used genomatix frameworker tool to analyze 41 promoter sequences of 22 genes annotated with table 3 lists the top five most significantly overrepresented promoter modules consisting of three vertebrate tf matrix families ( groups of functionally similar tfs ) ( 25 ) . the hoxf family tfs , the hbox family tf meox1 , the lhxf family tfs , and the homf family tf msx1 are downregulated , whereas the etsf family tfs are upregulated in the db / db nerve ( supplementary table 4 ) . b : expression profiles of 22  lipid metabolism  genes selected for promoter sequence analysis . top five most significantly overrepresented tfbs modules five promoter modules consisting of three elements were significantly overrepresented in the promoter sequences of lipid metabolism  related genes . module matrix families , strand orientation , percent of promoter sequences containing the module , and p values of enrichment are shown . the p value is the probability of finding the module in a set of randomly selected promoters and is determined by scanning a background promoter sequence set of 5,000 human promoters . table 1 lists phenotyping parameters for the db / db and db/+ mice measured at 24 weeks of age . the db / db mice are significantly heavier and exhibit higher levels of fasting blood glucose and ghb than the db/+ mice at 24 weeks of age . plasma cholesterol and triglyceride levels are significantly elevated in the db / db mice , confirming dyslipidemia in these mice . thermal latency in the hind paw and tail is significantly increased in db / db mice , indicative of sensory loss . both ienf density and sciatic motor nerve conduction velocity ( smncv ) are significantly lower in db / db mice . increased thermal latency , loss of ienf density , and decreased smncv confirm the development of dn in the db / db mice at 24 weeks of age . phenotyping parameters for the db / db and db/+ mice diabetes , dyslipidemia , and neuropathy phenotyping parameters measured at 24 weeks of age ( mean  sd ) . figure 1 outlines the steps in the transcriptional analysis of the microarray data . using intensity - based moderated t test false discovery rate < 0.05 as a threshold , we identified 4,017 degs ; 2,122 genes were upregulated and 1,895 genes were downregulated in the db / db relative to the db/+ samples . microarray data were validated using qrt - pcr for 5 selected degs , which demonstrated high correlation between microarray and qrt - pcr expression levels ( supplementary table 1 ) . rna samples obtained from the sciatic nerve ( scn ) of db / db and db/+ mice were hybridized to affymetrix genechip microarrays . enrichment analysis identifies functional categories that have a higher representation in the regulated genes and , thus , are likely to be active in the experimental condition . table 2 lists gene ontology ( go ) terms overrepresented in degs identified using david . overrepresented biological processes with a higher number of upregulated genes are cell cycle , lipid metabolic process , lipid transport , carbohydrate metabolic process , response to stress , and apoptosis . overrepresented cellular component terms are mitochondrion and axon , with a majority of mitochondrial genes upregulated and a majority of axonal genes downregulated . bh p value , benjamini - hochberg corrected p value of enrichment in david ; go_bp , go biological process ; go_cc , go cellular component . gene network analysis is useful for exploring both published and novel gene associations that may be involved in the process being studied . conceptgen is a gene set relation mapping tool that identifies relationships ( significant overlap ) between an experimental gene set and curated gene sets in sources such as go terms , kyoto encyclopedia of genes and genomes ( kegg ) pathways , and medical subject heading ( mesh ) terms ( 20 ) . conceptgen was used to analyze 717 degs annotated with the go biological processes shown in table 1 . the mouse gene identifiers were converted by conceptgen homology conversion tool to 695 human homologs . relationship among selected concepts ( go biological processes , kegg pathways , and mesh terms ) associated with the gene set were visualized using conceptgen network graph tool ( fig . kegg pathways with significant overlap with the deg set include fatty acid metabolism , peroxisome proliferator  activated receptor ( ppar ) signaling , and extracellular matrix  receptor interaction . concepts related to the deg set clustered in three groups : lipid metabolism , apoptosis , and axonogenesis ( fig . network depicting relationships and overlap among significantly overrepresented concepts ( go biological processes , kegg pathways , and mesh terms ) in the degs . node color represents concept type , node size represents number of genes in the concept , and thickness of the edge represents gene overlap between concepts . bibliosphere generates a network of genes based on their cocitation in pubmed abstracts and identifies functionally related subnetworks based on known regulatory relationships and pathway annotation . figure 3a illustrates a bibliosphere network of degs centered on tumor necrosis factor- ( tnf- ) with gene tnf- was chosen as the central node because of its importance in type 2 diabetes and lipid metabolism ( 22 ) . central nodes of subnetworks and associated pathways include hifna ( apoptosis ) , serpine1 ( p53 signaling ) , app ( axonogenesis ) , lyn ( fc  ri signaling ) , and pparg ( ppar signaling and fatty acid oxidation ) ( fig . nodes represent genes , node color represents fold change ( red / yellow : upregulation ; blue : downregulation ) . the network is centered at tnf-. only shortest paths to the central node are shown for clarity . b : enlarged view of the part of the network identified in a illustrates connections between the central node and other highly connected nodes in the network . highly connected nodes in the network include pparg , app , serpine1 , hif1a , irf1 , vdr , lyn , and ripk1 . ppi network of degs was constructed using cytoscape mimi plug - in ( 23 ) . subnetworks were identified based on pathway annotation of genes in mimi ( fig . regulated pathways identified by ppi network are carbohydrate metabolism ( glycolysis and tricarboxylic acid cycle [ tca ] ) , energy metabolism ( oxidative phosphorylation ) , lipid metabolism ( fatty acid metabolism and glycerolipid metabolism ) , glutathione metabolism , cell development ( cell cycle , apoptosis , axon guidance , and neurotrophin signaling ) , cell communication ( cell adhesion and extracellular matrix  receptor interaction ) , and signal transduction pathways ( ppar signaling , transforming growth factor- [ tgf- ] signaling , adipocytokine signaling , wnt signaling , mitogen - activated protein kinase [ mapk ] signaling , erbb signaling , and mammalian target of rapamycin [ mtor ] signaling ) . b : enlarged view of part of the network identified in a. node color represents pathway annotation of genes ( for example , blue : oxidative phoshorylation ; orange : tca ; red : apoptosis ) . supplementary table 2 lists degs in the pathways regulated in the db / db nerve . the majority of degs in the glycolysis , tca , oxidative phosphorylation , fatty acid metabolism , glycerolipid metabolism , mitochondrial fatty acid elongation , lipid transport , adipocytokine signaling , ppar signaling , and apoptosis pathways are upregulated , whereas most of the axonogenesis - related genes are downregulated in db / db nerve . crc , a model - based bayesian clustering tool , grouped the degs into 18 clusters based on positive as well as negative correlation among their expression profiles . we performed enrichment analysis for gene sets from individual clusters to identify overrepresented functional categories ( supplementary table 3 ) . gene regulatory sequence elements are organized into motifs of tfbss known as promoter modules , possibly suggesting coordinated regulation by multiple transcription factors ( tfs ) ( 24 ) . we used genomatix frameworker tool to analyze 41 promoter sequences of 22 genes annotated with table 3 lists the top five most significantly overrepresented promoter modules consisting of three vertebrate tf matrix families ( groups of functionally similar tfs ) ( 25 ) . the hoxf family tfs , the hbox family tf meox1 , the lhxf family tfs , and the homf family tf msx1 are downregulated , whereas the etsf family tfs are upregulated in the db / db nerve ( supplementary table 4 ) . b : expression profiles of 22  lipid metabolism  genes selected for promoter sequence analysis . top five most significantly overrepresented tfbs modules five promoter modules consisting of three elements were significantly overrepresented in the promoter sequences of lipid metabolism  related genes . module matrix families , strand orientation , percent of promoter sequences containing the module , and p values of enrichment are shown . the p value is the probability of finding the module in a set of randomly selected promoters and is determined by scanning a background promoter sequence set of 5,000 human promoters . how the metabolic dysregulation caused by insulin resistance leads to local nerve fiber damage in dn is unclear . in this study , we performed global gene expression profiling of sciatic nerve from 24-week - old db / db mice , which exhibit hyperglycemia , dyslipidemia , and dn typical of type 2 diabetes . we detected alterations in the expression of genes involved in carbohydrate and lipid metabolism , lipid transport , stress responses , and apoptosis . furthermore , promoter sequences of coexpressed lipid metabolism  related degs exhibited significantly overrepresented tfbs motifs and are likely to share a common regulatory mechanism relevant to the development of dn . morphological and electrophysiological changes in the bks db / db mice closely mimic the changes observed in dn patients ( 8,9,2628 ) . however , the db / db mouse model does not fully duplicate human dn pathology ; not all structural changes observed in human dn are seen in the db / db mouse ( 27,29 ) . duration of diabetes in mice is likely not long enough for the development of severe long - term degeneration associated with human dn . the development of hyperglycemia and dn in the bks db / db mice are considered to be consequences of hyperphagia and insulin resistance resulting from impaired leptin signaling in the hypothalamus ( 8) . however , the role of leptin signaling in peripheral nerve is not understood , and whether leptin receptor mutation and obesity have a compounding effect on the development of dn in these mice is not known . schwann cells ( sc ) are major contributors to the mrna in the sciatic nerve biopsies , with a small contribution coming from axons , epineural fibroblasts , adipocytes , and vascular endothelial cells ( 30 ) . our analyses revealed that genes encoding myelin structural proteins ( p0 , pmp22 , connexin-32 , mag , e - cadherin , and periaxin ) , myelin structural lipid synthesis genes ( plp , elovl6 , and apolipoprotein d ) , as well as sox10 , a tf required for myelination , are downregulated in db / db nerve . segmental demyelination has been observed in human dn ( 31,32 ) ; mouse studies , however , do not report evidence of demyelination ( 27,29 ) . downregulation of myelin protein - encoding genes and myelin lipid synthesis genes may represent sc abnormalities preceding the structural changes . the tfs c - jun and krox-24 , and other genes expressed by denervated sc to promote axonal regeneration ( gap-43 , ncam , and l1 ) ( 33 ) , are downregulated , as are positive regulators of axonogenesis ( cdh4 , nefl , and mapt ) and axon guidance pathway genes ( netrin g1 and fyn ) , suggesting a lack of axonal regeneration . these changes are consistent with the axonal degeneration observed in morphological studies in the db / db mice ( 27,29 ) . loss of neurotrophic signal in dn may impair nerve generation and cause dying back of axons ( 3,34 ) . neurotrophin ngf ; nt-3 receptor trkc ; gdnf receptors gfra1 , gfra3 , and gfra4 ; neuregulin receptor erbb ; and neurotrophin signaling pathway genes ( akt3 and mapk10 ) are downregulated , indicating impaired neurotrophic support in the db / db nerve . downregulation of ngf , trkc , and gfr1 is consistent with the results of studies in diabetic rats ( 3537 ) . some studies in rats with streptozotocin - induced diabetes report reduced expression of bdnf , nt-3 , nt-4/5 , igf - i , igf - ii , cntf , gdnf , and neurotrophin receptors p75 and trkb approximately 12 weeks after the induction of diabetes ( 3540 ) , whereas others report no change in nt-3 expression in rat sciatic nerve ( 41 ) . the difference in the animal models used in the studies ( rat model of type 1 diabetes vs. mouse model of type 2 diabetes ) may explain this discrepancy . the sc marker s100 as well as antiapoptotic genes ( mapk8ip1 and tgf-2 ) are downregulated , whereas proapoptotic genes ( fas , tnf- , casp8 , and brca1 ) are upregulated in the db / db samples ( supplementary table 2 ) . the downregulation of sc marker and antiapoptotic genes along with upregulation of proapoptotic genes suggests sc apoptosis in db / db nerve ; however , no studies to date have demonstrated sc death in dn . the causes of axonal and sc degradation in dn are not well understood , but several hypotheses have been developed in regard to the pathogenesis of the nerve injury ( 3,34 ) . our analyses indicate increased glucose , energy , and lipid metabolism in the db / db nerve and support the roles of hyperglycemia- and hyperlipidemia - induced oxidative stress , inflammatory response , and vascular ischemia in dn . in addition , our data indicate impaired axonal transport , neurotrophic signal , cell adhesion , and cell communication . in this study , we focused on hyperglycemia- and dyslipidemia - related gene expression changes in the nerve . hyperglycemia is a major factor in the development and progression of dn ; current hypotheses suggest the effect of hyperglycemia is likely to be vascular , metabolic , or a combination of both ( 34,42 ) . the metabolic hypothesis of axonal and sc damage in diabetes suggests that activation of glucose metabolism pathways in hyperglycemia results in oxidative stress . in our data , glycolysis , tca , and oxidative phosphorylation genes are upregulated ( supplementary table 2 ) , suggesting activation of glucose and energy metabolism in the db / db nerve . upregulation of oxidative phosphorylation genes and downregulation of mitochondrial h transporting atp synthases in the db / db nerve suggest increased superoxide production ( 43 ) . the mitochondrial proton carrier ucp2 , known to be upregulated in response to elevated reactive oxygen species ( ross ) ( 44 ) , is highly upregulated ; the oxidative stress - induced growth inhibitor osgin1 is also upregulated in the db / db nerve . upregulation of antioxidant genes ( sod2 , peroxiredoxins , and catalase ) suggests cellular response to increased ros production . ( 3 ) hypothesized that increased cellular oxidative stress results in activation of the poly ( adp - ribose ) polymerase ( parp ) pathway , which in turn induces inflammatory responses in the nerve . upregulation of parp , inflammatory response , and mapk signaling pathway genes in the db / db nerve support the parp pathway  mediated inflammatory response hypothesis . nuclear factor-b ( nf-b ) induced  inflammatory response is also implicated in demyelination and axon degeneration ( 3 ) ; however , nf-b and genes regulated by nf-b , such as inos , are not differentially regulated in our data . the vascular hypothesis of nerve damage in diabetes suggests that activation of glucose metabolism pathways causes functional and structural changes in the neuronal vasculature leading to endothelial hypoxia and axonal ischemia ( 45 ) . ( 46 ) noted thickening of perineurial cell basal lamina and increased endoneurial collagenization around sc in the sural nerve of dn patients . activation of protein kinase c by hyperglycemia induces expression of vascular endothelial growth factor ( vegf)an angiogenic gene and permeability factor  inhibits production of nitric oxide and alters na - k - atpase activity . activation of glucose metabolism pathways induces tgf-1 and serpine1 expression , which results in endothelial fibrosis and thickening of vascular membranes ( 3 ) . vegf - c , tgf-1 , and serpine1 are upregulated ; na - k - atpases ( atp1b2 and atp1b3 ) are downregulated in the db / db nerve . upregulation of hypoxia - inducible factor hif1a may indicate cellular response to ischemia ( 47 ) . imbalance in expression of nitric oxide synthase ( nos ) isoforms has been implicated in diabetic endothelial dysfunction ( 48 ) ; however , these are not differentially expressed in our data . dysregulation of lipid metabolism is implicated in peripheral neuropathy ( 7,30,33,49 ) , but the mechanism of its effect is not well understood . we observed upregulation of genes expressed by adipocytes ( adiponectin , perilipin , lipin1 , and lpl ) , possibly due to increased amounts of adipose tissue in the epineurium of the db / db nerve . ( 30 ) demonstrated that epineurium of adult mouse sciatic nerve contains adipose - like tissue , and both endoneurium and epineurium express lipid metabolism genes ; they suggested that adipocytokine signaling from endoneurial adipocytes may regulate lipid metabolism in sc and axons . our analysis indicates upregulation of fatty acid metabolism , glycerolipid metabolism , lipid transport , ppar signaling , and adipocytokine signaling genes in the db / db nerve ( supplementary table 2 ) . vincent et al . ( 7,50 ) hypothesized that increased production of ros in hyperglycemia results in lipid peroxidation and that the formation of oxidized low density lipoproteins ( oxldl ) contributes to the nerve injury in diabetic dyslipidemia . the apoptotic effect of oxldl on vascular endothelial cells and neurons is mediated via oxldl receptor ( olr1/lox-1 ) , which is induced by hyperglycemia , ros , oxldl , and tnf-. in vitro studies demonstrate that binding of oxldl to olr1 leads to nadph oxidase ( nox ) activation , mitochondrial superoxide production , and neuronal oxidative stress in dorsal root ganglia of mice ( 50 ) . oxldl receptor olr1 is highly overexpressed and nox2 is upregulated in our data , supporting the hypothesis of oxldl - induced nerve injury in dyslipidemia , possibly via activation of nox . oxldl also acts on macrophages via scavenger receptor ( scarb1/cd36 ) and triggers immune responses via mapk and nf-b signaling pathways ( 7 ) . in our data , cd36 and mapk signaling pathway genes ( tnf- , il1a , and tgf-1 ) are upregulated , supporting the hypothesis of cd36-mediated inflammatory response in the db / db nerve . lipid metabolism ,   mitochondrion ,  and  axonogenesis  clustered together based on their correlated expression pattern ( fig . most of the lipid metabolism and mitochondria - associated genes in the cluster are upregulated , whereas axonogenesis genes are downregulated . to investigate the potential coregulation mechanism , we analyzed promoter sequences of 22 degs from this cluster annotated with the term  lipid metabolism  ; 5 of these genes are also annotated with it is interesting that tfs binding to the overrepresented tfbss in promoter sequences of lipid metabolism genes are associated with neuron differentiation , axon guidance , and immune response , suggesting a possible common regulatory mechanism in lipid metabolism , stress response , and axonal degeneration . the hoxf family , the hbox family , and the lhxf family tfs are annotated with  motor axon guidance  and  regulation of neuron differentiation . negative regulation of vascular endothelial growth factor receptor signaling pathway  and  neuron differentiation . immune responseassociated etsf family tfs are upregulated in the db / db nerve ( supplementary table 4 ) . oct1 family tf oct-6 , an sc tf important in myelination ( 30 ) , is not among the degs . in conclusion , the changes are consistent with pathological characteristics of dn , such as axonal degeneration and potential loss of neurotrophic signal . our findings support the role of hyperglycemia - induced oxidative stress and ischemia in nerve injury . our results also support the hypothesis of oxidized lipid  mediated nerve injury and increased mitochondrial oxidative stress in dyslipidemia . in addition , our analyses revealed possible coregulation of lipid metabolism , stress response , and axonal degeneration genes and identified tfs that may modulate this coregulation . further investigation into the signal mediated by these tfs is likely to provide more insight into dyslipidemia - induced nerve injury .
objectivea better understanding of the molecular mechanisms underlying the development and progression of diabetic neuropathy ( dn ) is essential for the design of mechanism - based therapies . 
 we examined changes in global gene expression to define pathways regulated by diabetes in peripheral nerve.research design and methodsmicroarray data for 24-week - old bks db / db and db/+ mouse sciatic nerve were analyzed to define significantly differentially expressed genes ( degs ) ; degs were further analyzed to identify regulated biological processes and pathways . 
 expression profile clustering was performed to identify coexpressed degs . 
 a set of coexpressed lipid metabolism genes was used for promoter sequence analysis.resultsgene expression changes are consistent with structural changes of axonal degeneration . 
 pathways regulated in the db / db nerve include lipid metabolism , carbohydrate metabolism , energy metabolism , peroxisome proliferator  activated receptor signaling , apoptosis , and axon guidance . 
 promoter sequences of lipid metabolism  related genes exhibit evidence of coregulation of lipid metabolism and nervous system development genes.conclusionsour data support existing hypotheses regarding hyperglycemia - mediated nerve damage in dn . 
 moreover , our analyses revealed a possible coregulation mechanism connecting hyperlipidemia and axonal degeneration .