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

One of the challenges in exploiting high throughput measurement techniques such as microarrays is the conversion of the vast amounts of data obtained into relevant knowledge. Of particular importance is the identification of the intrinsic response of a transcriptional experiment and the characterization of the underlying dynamics.

Methodology and Findings

The proposed algorithm seeks to provide the researcher a summary as to various aspects relating to the dynamic progression of a biological system, rather than that of individual genes. The approach is based on the identification of smaller number of expression motifs that define the transcriptional state of the system which quantifies the deviation of the cellular response from a control state in the presence of an external perturbation. The approach is demonstrated with a number of data sets including a synthetic base case and four animal studies. The synthetic dataset will be used to establish the response of the algorithm on a “null” dataset, whereas the four different experimental datasets represent a spectrum of possible time course experiments in terms of the degree of perturbation associated with the experiment as well as representing a wide range of temporal sampling strategies. This wide range of experimental datasets will thus allow us to explore the performance of the proposed algorithm and determine its ability identify relevant information.

Conclusions and Significance

In this work, we present a computational approach which operates on high throughput temporal gene expression data to assess the information content of the experiment, identify dynamic markers of important processes associated with the experimental perturbation, and summarize in a concise manner the evolution of the system over time with respect to the experimental perturbation.