Effects of high fat diet (HFD) on obesity and, subsequently, on diabetes are highly variable and modulated by genetics in both humans and rodents. In this report, we characterized the response of Goto-Kakizaki (GK) rats, a spontaneous polygenic model for lean diabetes and healthy Wistar-Kyoto (WKY) controls, to high fat feeding from weaning to 20 weeks of age. Animals fed either normal diet or HFD were sacrificed at 4, 8, 12, 16 and 20 weeks of age and a wide array of physiological measurements were made along with gene expression profiling using Affymetrix gene array chips. Mining of the microarray data identified differentially regulated genes (involved in inflammation, metabolism, transcription regulation, and signaling) in diabetic animals, as well as the response of both strains to HFD. Functional annotation suggested that HFD increased inflammatory differences between the two strains. Chronic inflammation driven by heightened innate immune response was identified to be present in GK animals regardless of diet. In addition, compensatory mechanisms by which WKY animals on HFD resisted the development of diabetes were identified, thus illustrating the complexity of diabetes disease progression.

The Goto-Kakizaki (GK) rat, a polygenic non-obese model of type 2 diabetes, is a useful surrogate for study of diabetes-related changes independent of obesity. GK rats and appropriate controls were killed at 4, 8, 12, 16 and 20 weeks post-weaning and differential muscle gene expression along with body and muscle weights, plasma hormones and lipids, and blood cell measurements were carried out. Gene expression analysis identified 204 genes showing 2-fold or greater differences between GK and controls in at least 3 ages. Array results suggested increased oxidative capacity in GK muscles, as well as differential gene expression related to insulin resistance, which was also indicated by HOMA-IR measurements. In addition, potential new biomarkers in muscle gene expression were identified that could be either a cause or consequence of T2DM. Furthermore, we demonstrate here the presence of chronic inflammation evident both systemically and in the musculature, despite the absence of obesity.

Circadian rhythms (24 h cycles) are observed in virtually all aspects of mammalian function from expression of genes to complex physiological processes. The master clock is present in the suprachiasmatic nucleus (SCN) in the anterior part of the hypothalamus and controls peripheral clocks present in other parts of the body. Components of this core clock mechanism regulate the circadian rhythms in genome-wide mRNA expression, which in turn regulate various biological processes. Disruption of circadian rhythms can be either the cause or the effect of various disorders including metabolic syndrome, inflammatory diseases and cancer. Furthermore, circadian rhythms in gene expression regulate both the action and disposition of various drugs and affect therapeutic efficacy and toxicity based on dosing time. Understanding the regulation of circadian rhythms in gene expression plays an important role in both optimizing the dosing time for existing drugs and in development of new therapeutics targeting the molecular clock.

Type 2 diabetes (T2DM) is a heterogeneous group of diseases that is progressive and involves multiple tissues. Goto-Kakizaki (GK) rats are a polygenic model with elevated blood glucose, peripheral insulin resistance, a non-obese phenotype, and exhibit many degenerative changes observed in human T2DM. As part of a systems analysis of disease progression in this animal model, this study characterized the contribution of adipose tissue to pathophysiology of the disease. We sacrificed subgroups of GK rats and appropriate controls at 4, 8, 12, 16 and 20 weeks of age and carried out a gene array analysis of white adipose tissue. We expanded our physiological analysis of the animals that accompanied our initial gene array study on the livers from these animals. The expanded analysis included adipose tissue weights, HbA1c, additional hormonal profiles, lipid profiles, differential blood cell counts, and food consumption. HbA1c progressively increased in the GK animals. Altered corticosterone, leptin, and adiponectin profiles were also documented in GK animals. Gene array analysis identified 412 genes that were differentially expressed in adipose tissue of GKs relative to controls. The GK animals exhibited an age-specific failure to accumulate body fat despite their relatively higher calorie consumption which was well supported by the altered expression of genes involved in adipogenesis and lipogenesis in the white adipose tissue of these animals, including Fasn, Acly, Kklf9, and Stat3. Systemic inflammation was reflected by chronically elevated white blood cell counts. Furthermore, chronic inflammation in adipose tissue was evident from the differential expression of genes involved in inflammatory responses and activation of natural immunity, including two interferon regulated genes, Ifit and Iipg, as well as MHC class II genes. This study demonstrates an age specific failure to accumulate adipose tissue in the GK rat and the presence of chronic inflammation in adipose tissue from these animals.

Background

Circadian rhythms are 24 hour oscillations in many behavioural, physiological, cellular and molecular processes that are controlled by an endogenous clock which is entrained to environmental factors including light, food and stress. Transcriptional analyses of circadian patterns demonstrate that genes showing circadian rhythms are part of a wide variety of biological pathways.

Pathway activity method can identify the significant pattern of the gene expression levels within a pathway. In this method, the overall gene expression levels are translated to a reduced form, pathway activity levels, via singular value decomposition (SVD). A given pathway represented by pathway activity levels can then be as analyzed using the same approaches used for analyzing gene expression levels. We propose to use pathway activity method across time to identify underlying circadian pattern of pathways.

Results

We used synthetic data to demonstrate that pathway activity analysis can evaluate the underlying circadian pattern within a pathway even when circadian patterns cannot be captured by the individual gene expression levels. In addition, we illustrated that pathway activity formulation should be coupled with a significance analysis to distinguish biologically significant information from random deviations. Next, we performed pathway activity level analysis on a rich time series of transcriptional profiling in rat liver. The over-represented five specific patterns of pathway activity levels, which cannot be explained by random event, exhibited circadian rhythms. The identification of the circadian signatures at the pathway level identified 78 pathways related to energy metabolism, amino acid metabolism, lipid metabolism and DNA replication and protein synthesis, which are biologically relevant in rat liver. Further, we observed tight coordination between cholesterol biosynthesis and bile acid biosynthesis as well as between folate biosynthesis, one carbon pool by folate and purine-pyrimidine metabolism. These coupled pathways are parts of a sequential reaction series where the product of one pathway is the substrate of another pathway.

Conclusions

Rather than assessing the importance of a single gene beforehand and map these genes onto pathways, we instead examined the orchestrated change within a pathway. Pathway activity level analysis could reveal the underlying circadian dynamics in the microarray data with an unsupervised approach and biologically relevant results were obtained.