<title>Introduction</title> <p>Cardiovascular disease is the leading cause of morbidity and mortality in patients with type 2 diabetes ##REF##12556246##[1]##. While atherosclerotic coronary artery disease is highly prevalent in many diabetic patients, the occurrence of nonischemic cardiomyopathy (“diabetic cardiomyopathy”) suggests that other processes such as microangiopathy, metabolic factors, or myocardial fibrosis ##REF##12556246##[1]## might be involved. Thiazolidinediones (TZDs) such as rosiglitazone are insulin sensitizers that may also have beneficial properties for diabetic cardiomyopathy. These agonists bind to a subfamily of nuclear hormone receptors, the peroxisome proliferator-activated receptors (PPARs, including α, β/δ and γ isoforms, of which the γ isoform is the most common target), and activate transcription factors that modulate gene expression, ultimately leading to increased insulin sensitivity in peripheral tissues through poorly defined pathways ##UREF##0##[2]##, ##REF##10953021##[3]##. However, a recent meta-analysis of clinical trial outcomes for the PPARγ agonist rosiglitazone found a significant increase in the risk of myocardial infarction, and a borderline significant risk of death from cardiovascular causes ##REF##17517853##[4]##. The reasons for this increased risk are unclear, and a better understanding of the effects of rosiglitazone on the diabetic heart is urgently needed.</p> <p>Despite the large number of studies on rosiglitazone and diabetes, the transcriptional network changes by which PPARγ agonists and other TZDs promote insulin sensitivity are not well characterized. A number of groups have studied gene expression in the <italic>db/db</italic> mouse liver ##REF##16752183##[5]##, as well as transcriptional changes induced by TZDs in adipocytes, kidney, aorta, and pancreatic islet cells in various mouse models ##REF##17003339##[6]##, ##REF##11916924##[7]##, ##REF##15054141##[8]##, ##REF##15126236##[9]##, but none has studied the heart directly. Moreover, determining PPAR-agonist target tissues is complicated by the fact that the three PPAR isoforms are not distributed equally across all tissues. In adipose tissue, PPARγ predominates and promotes differentiation and lipid storage ##REF##11459852##[10]##, resulting in suppression of lipolysis by insulin and reduction of plasma free fatty acid concentrations. In the cardiovascular system, PPARγ is expressed in vascular smooth muscle and vascular endothelium ##REF##9831704##[11]##, but is believed to be in low abundance in cardiomyocytes ##REF##12649261##[12]##.</p> <p>Animal models of type 2 diabetes present an optimal system for studying the effects of rosiglitazone. A well-established murine model is the C57BL/KLS-lepr<sup>db</sup>/lepr<sup>db</sup> (<italic>db/db</italic>) mouse that has a mutation in the leptin receptor. Homozygous <italic>db/db</italic> mice become obese by 3–4 wk of age and develop hyperglycemia at 4–8 wk. Serum insulin levels increase as early as 10–14 days, peak at 6–8 wk, then decrease afterward. These mice continue to be hyperinsulinemic throughout life and ultimately develop cardiomyopathy (contractile dysfunction), as evidenced from metabolic experiments using cultured <italic>db/db</italic> cardiomyocytes ##REF##12540618##[13]##, echocardiographic studies ##REF##12181126##[14]##, isolated working heart preparations ##REF##11052966##[15]##, and direct <italic>in situ</italic> ventricular pressure measurements ##REF##16141388##[16]##. Recent studies by our lab and others have delineated the longitudinal structural and metabolic cardiomyopathic changes in the <italic>db/db</italic> mouse using cardiac magnetic resonance (CMR) and Fluorine-18–2-fluoro-2-deoxy-d-glucose ([<sup>18</sup>F]FDG) PET scanning ##REF##17122193##[17]##.</p> <p>To address the many questions concerning the molecular and functional effects of rosiglitazone on the diabetic heart, we used CMR imaging and echocardiagraphic studies to assess the development of cardiomyopathy in untreated and rosiglitazone-treated <italic>db/db</italic> mice, as well as normal <italic>db/+</italic> mice. These studies were followed by cardiac gene expression profiling of the diabetic and control phenotypes to understand the underlying molecular changes in heart tissue. A comparison of the cardiac transcriptomes of rosiglitazone-treated <italic>db/db</italic> mice with untreated <italic>db/db</italic> and normal <italic>db/+</italic> mice thus provides insights into the effects of rosiglitazone on the heart, either through direct or indirect mechanisms, as well as clues to its potential toxicity.</p>
<title>Methods</title> <title>Animal and sample tissue preparation</title> <p>Homozygous <italic>db/db</italic> mice (Jackson Laboratories, Bar Harbor, ME) were maintained on a normal chow diet and housed in a room with a 12∶12-h light-dark cycle and an ambient temperature of 22°C. For the treatment group, homozygous <italic>db/db</italic> mice were maintained on 5 mg/kg/day (approximately 0.225 mg/mouse/day) rosiglitazone (GlaxoSmithKline, London, UK) that was mixed with normal chow. Unless otherwise stated, heterozygous <italic>db/+</italic> littermates were used as control animals. All protocols were approved by the Administrative Panel on Laboratory Animal Care at Stanford University and were carried out in accordance with the guidelines of the American Association for Accreditation of Laboratory Animal Care. At the end of the study, animals were euthanized with a lethal dose of isoflurane. Immediately after death, wet heart weight (HW) and body weight (BW) were measured. Whole hearts, pancreas, and liver were harvested and preserved in TRIzol reagent (Invitrogen, San Diego, CA) for subsequent mRNA isolation, and a subset of hearts were fixed in 10% formalin for histological evaluation.</p> <title>Insulin tolerance testing</title> <p>Insulin tolerance testing was performed on mice after a 6-h fast. At the time of testing, a bolus of human regular insulin (Eli Lilly, Indianapolis, IN; 0.75 IU/kg) was injected intraperitoneally. In blood derived from a tail nick, glucose levels were then determined with a FreeStyle blood glucose monitoring system (Abbott Laboratories, Abbott Park, IL) at baseline and 30, 60, and 120 min after injection.</p> <title>Insulin Enzyme-Linked ImmunoSorbent Assay (ELISA)</title> <p>After a 6-h fast, roughly 200 µl of blood was obtained from each animal via retroorbital bleeding. Samples were placed in 500- µl tubes containing EDTA (Becton-Dickinson, Franklin Lakes, NJ) and centrifuged at 4°C and 13,200 rpm for 10 min. Approximately 50–75 µl of supernatant (serum) was then collected for further processing. Serum insulin concentrations were measured with a Mercodia mouse insulin enzyme immunoassay kit (Alpco Diagnostics, Salem, NH).</p> <title>Nonesterified fatty acids (NEFA)</title> <p>Plasma nonesterified fatty acid concentrations were measured using a HR series NEFA-HR(2) colorimetric assay kit (Wako Chemicals USA, Richmond, VA).</p> <title>Left ventricular functional analysis with echocardiogram</title> <p>Echocardiography was performed by a blinded investigator (ZL) using the Siemens-Acuson Sequioa C512 system equipped with a multi-frequency (8–14 MHz) 15L8 transducer. Analysis of M-mode images was performed using the Siemens built-in software. Left ventricular end diastolic diameter (EDD) and end-systolic diameter (ESD) were measured and used to calculate fractional shortening (FS) by the following formula: FS = [EDD-ESD]/EDD. LV volume at end diastolic (EDV) and end-systole (ESV) were calculated by the bullet method as follows: EDV = 0.85×CSA(d)×L(d), ESV = 0.85×CSA(s)×L(s), where CSA(d) and (s) are endocardial area in end-diastole and end-systole, respectively, obtained from short-axis view at the level of the papillary muscles. L(d) and L(s) are the LV length (apex to mid-mitral annulus plane) in end-diastole and end-systole, respectively, obtained from the parasternal long-axis view. LV ejection fraction (EF%) was calculated as: EF% = (EDV−ESV)×100/EDV.</p> <title>Cardiac magnetic resonance imaging</title> <p>To prepare for scanning, induction of anesthesia was accomplished with 2% isoflurane and 1 l/min oxygen. Respiratory rate was monitored and used to manually calibrate the maintenance dose of isoflurane at 1.25–1.5%. Platinum needle ECG leads were inserted subcutaneously in the right and left anterior chest wall. Respiration was monitored with a pneumatic pillow sensor positioned along the abdomen. Body temperature was maintained at 36–37°C by a flow of heated air thermostatically controlled by a rectal temperature probe. Heart rate (HR), respiratory rate, and body temperature were recorded every 4 min during image acquisition. Magnetic resonance images were acquired by a blinded investigator (PY) with a 4.7-T magnet (Bruker BioSpin, Fremont, CA) controlled by a Varian Inova Console (Varian, Palo Alto, CA), using a transmit-receive quadrature volume coil with an inner diameter of 3.5 cm. For particularly obese animals, a larger coil with an inner diameter of 6 cm was utilized. Image acquisition was gated to the ECG R-wave (Small Animal Instruments, Stony Brook, NY). Coronal and axial scout images were used to position a two-dimensional imaging plane along the short axis of the left ventricular (LV) cavity. Gated gradient echo sequences were then used to acquire sequential short-axis slices spaced 1 mm apart from apex to base. For each sequence, 12 cine frames encompassing one cardiac cycle were obtained at each slice level with the following sequence parameters: acquisition time (TR) = 100–140 ms, echo time (TE) = 2.8–3.5 ms, number of repeats (NEX) = 8, field of view (FOV) = 30×30 mm, matrix = 128×128, flip angle = 60°. For each short-axis slice, planimetry measurements of LV myocardial area were conducted off-line by tracing the epicardial and endocardial borders at end systole and end diastole with MRVision software (MRVision, Winchester, MA). For these purposes the papillary muscles were considered part of the LV cavity. Anteroseptal and posterior wall thickness measurements were performed on a midventricular slice at the level of the papillary muscles at end diastole. LV mass (LVM) was derived from the sum of the differences between the end-diastolic epicardial and endocardial areas from apex to base, adjusted for the specific gravity of myocardial tissue (1.055 g/ml). LV end-diastolic (LVEDV) and end-systolic (LVESV) volumes were calculated as the sum of the endocardial areas of each slice from the apex to the LV outflow tract at end diastole and end systole, respectively. LV ejection fraction (LVEF) was calculated as (LVEDV−LVESV)/LVEDV. Cardiac output (CO) was calculated as (LVEDV−LVESV)/(HR).</p> <title>RNA isolation and quality control</title> <p>Total RNA was isolated separately from five heart samples for each condition using TRIzol reagent (Invitrogen, San Diego, CA) according to the manufacturer's instructions for a total of 15 distinct RNA samples. Total RNA was purified using RNeasy columns (Qiagen, Chatsworth, CA). RNA concentration was measured by spectrophotometry, and RNA integrity assessed with an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA) with 6000 Nano Chips according to the manufacturer's instructions. RNA was judged as suitable for array hybridization only if samples exhibited intact bands corresponding to 18S and 28S ribosomal RNA subunits and had a RNA Integrity Number (RIN) greater than six. Universal Reference RNAs for mouse (Stratagene, La Jolla, CA) were purchased as microarray reference controls.</p> <title>Labeling reaction and hybridization</title> <p>Using Low RNA Input Fluorescent Linear Amplification Kits (Agilent Technologies, Santa Clara, CA, USA), cDNA was reverse transcribed from each RNA sample, and cRNA was then transcribed and fluorescently labeled from each cDNA sample. The fifteen mouse heart tissue cRNA samples were labeled with Cy5 and the Universal Mouse reference cRNA was labeled with Cy3, yielding a total of 5 biological replicates per condition. The resulting cRNA was purified using an RNeasy kit (Qiagen, Valencia, CA, USA) followed by quantification of the cRNA by spectroscopy using an ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). 825 ng Cy3- and Cy5- labeled and amplified cRNA was mixed and fragmented according to the Agilent technology protocol. cRNA was hybridized to 4×44K whole human genome microarray slides from Agilent (Part G4112F) according to the manufacturer's instructions. The hybridization was carried in a rotating hybridization chamber in the dark at 65°C for 17 h.</p> <title>Washing, scanning, and feature extraction</title> <p>Slides were washed with Gene Expression Wash Buffer 1 and 2 (Agilent Technologies, Santa Clara, CA) followed by Acetonitrile. A final wash in Stabilization and Drying Solution was performed to prevent Cy-5 degradation by ozone. The array was scanned using an Agilent G2505B DNA microarray scanner under extended Dynamic range. The image files were extracted using the Agilent Feature Extraction software version 9.5.1 applying LOWESS background subtraction and dye-normalization. Further analyses were performed with BRB ArrayTools Version 3.4 (National Cancer Institute) and TIGR MeV software (<ext-link ext-link-type="uri" xlink:href="http://www.tm4.org/">http://www.tm4.org/</ext-link>).</p> <title>Real-time quantitative PCR for confirmation of microarray results</title> <p>RNA to cDNA conversion was performed with a SuperScript First-Strand cDNA synthesis kit (Invitrogen). The cDNA was then used as a template in a TaqMan real-time PCR assay with the ABI Prism 7700 Sequence Detection System (Applied Biosystems, Foster City, CA). All samples were run in triplicates. Specialized premade gene expression assay reagents for potassium channel, subfamily K, member 1 (Kcnk1; catalog no. Mm00434624_m1), matrix metalloproteinase 3 (Mmp3; Mm00440295_m1), beta-1,4-N-acetyl-galactosaminyl transferase 1 (B4galnt1; Mm00484653_m1), carboxyl ester lipase (Cel; Mm00486975_m1), and N-acylsphingosine amidohydrolase 2 (Asah2; Mm00479659_m1), purchased from Applied Biosystems, were used for these experiments. Threshold cycles were placed in the logarithmic portion of the amplification curve, and each sample was referenced to 18S RNA (part no. 4319413E) amplification to control for the total amount of RNA. Fold differences between samples (relative quantification) were calculated with the delta-delta method [S1/S2 = 2<sup>−(T1–T2)</sup>], where S1 and S2 represent samples 1 and 2 and T1 and T2 denote the threshold cycles for S1 and S2.</p> <title>Histologic analysis of cardiac tissue</title> <p>At the time of death, a select group of hearts were fixed in 10% formalin, cut at the midventricular level, and embedded in paraffin blocks. The blocks were then sectioned into short-axis slices, which were subsequently stained with hematoxylin and eoxin according to standard protocols. High-power fields magnified to 40× from the midportion of the LV free wall were photographed.</p> <title>Statistical Analysis</title> <p>For all tests, one-way ANOVA was performed. All <italic>P</italic> values &lt;0.05 were considered significant. For microarray analysis, the Significance Analysis of Microarrays (SAM) algorithm was used to identify genes with statistically significant differences in expression among the conditions. SAM is a statistically rigorous test that incorporates a FDR calculation to correct for multiple testing errors. Gene Ontology group overrepresentation analysis was performed using Fisher's exact test with FDR correction through High Throughput GoMiner software.</p>
<title>Results</title> <title>Insulin resistance in rosiglitazone-treated, untreated and control mice</title> <p>We analyzed a number of metabolic parameters to confirm insulin resistance in our <italic>db/db</italic> murine model of diabetes before proceeding with microarrays. Insulin resistance is strongly associated with obesity, and one mechanism may be the generation of metabolic messengers such as free fatty acids by adipose tissue that inhibit insulin action on muscle ##REF##8971073##[18]##, ##REF##11916933##[19]##. We measured plasma non-esterified fatty acid (NEFA) and insulin levels to follow insulin resistance in our three groups of mice. ##FIG##0##
Figure 1a
## shows plasma NEFA levels in the three groups of mice. Compared to untreated <italic>db/db</italic> mice, NEFA levels in the rosiglitazone-treated <italic>db/db</italic> mice decreased dramatically. Fasting insulin levels (##FIG##0##
Figure 1b
##) and insulin tolerance testing (##FIG##0##
Figure 1c
##) further confirmed the improved insulin sensitivity with rosiglitazone treatment.</p> <title>Body and heart morphometric analysis</title> <p>In addition to the obesity associated with type 2 diabetes, rosiglitazone and other TZDs are widely known to cause further weight gain due to fluid retention and altered adipogenesis ##REF##14693998##[20]##. Mean body weights of both treated and untreated <italic>db/db</italic> mice were markedly higher than <italic>db/+</italic> controls at baseline (##FIG##0##
Figure 1d
##). Furthermore, rosiglitazone-treated mice had generally higher mean body weights when compared to untreated mice, as expected. The heart/body weight ratio was significantly higher in the <italic>db/+</italic> group when compared to the treated and untreated <italic>db/db</italic> groups (##FIG##0##
Figure 1e
##). However, histological analysis did not reveal obvious differences at the microscopic level (##FIG##1##
Figure 2
##).</p> <title>Left ventricular functional analysis with echocardiogram</title> <p>Previous groups have observed a significant decrease in left ventricular ejection fraction (LVEF) and fractional shortening (FS) in <italic>db/db</italic> mice when compared with <italic>db/+</italic> controls ##REF##12181126##[14]##, ##REF##14600074##[21]##. In one study, the FS was reduced by as much as 16% at 12 weeks of age ##REF##12181126##[14]##. However, more accurate CMR studies of cardiac contractility have not confirmed the dramatic differences seen in the echocardiographic studies, observing only an approximately 2% reduction in LVEF at 22 weeks of age ##REF##17122193##[17]##. Furthermore, clinical echocardiographic studies of type 2 diabetic patients treated with rosiglitazone found no significant changes in LVEF after 52 weeks of treatment ##REF##17448371##[22]##. To investigate cardiac contractility and verify these previous findings, we assessed cardiac contractility in the three mouse groups using echocardiography (##FIG##2##
Figure 3
##). The mean LVEF and FS levels of the untreated <italic>db/db</italic> group were not significantly different from those in the <italic>db/+</italic> group, though there does appear to be a downward trend in contractility in the <italic>db/db</italic> group. Further, we found no significant LVEF or FS differences between the rosiglitazone-treated <italic>db/db</italic> group and untreated <italic>db/db</italic> group, confirming the results of the clinical study that found no evidence of rosiglitazone-induced changes in cardiac contractility.</p> <title>Cardiac magnetic resonance imaging</title> <p>Using cardiac magnetic resonance (CMR) scanning, our group has previously shown progressive cardiomyopathic changes in <italic>db/db</italic> mice when compared to <italic>db/+</italic> controls ##REF##17122193##[17]##. Left ventricular mass (LVM), interventricular septal thickness (IVST) and posterior wall thickness (PWT) were significantly increased in <italic>db/db</italic> mice, while LVEF was only marginally decreased. For the CMR imaging in this new study (##FIG##3##
Figure 4
##), we included rosiglitazone-treated <italic>db/db</italic> mice in addition to the untreated <italic>db/db</italic> and <italic>db/+</italic> groups, and scanned them at 8 weeks after treatment initiation to look for cardiomyopathic changes. CMR scans revealed significant, albeit subtle, cardiomyopathic changes in untreated <italic>db/db</italic> mice compared to <italic>db/+</italic> controls (##FIG##3##
Figure 4
##). However, rosiglitazone treatment resulted in no significant improvements in cardiac contractility, which confirmed our echocardiographic studies.</p> <title>Gene expression analysis</title> <p>Whole-heart RNA from five mice from each of the three groups after four months with or without treatment was used for microarray analysis. Overall, there were significant transcriptional differences between the <italic>db/+</italic> and both <italic>db/db</italic> groups, while differences between rosiglitazone-treated and untreated <italic>db/db</italic> groups were much less dramatic. ##FIG##4##
Figure 5a
## shows a 3-dimensional scaling plot based on principal component analysis of the expression data, which graphically demonstrates that rosiglitazone did not restore <italic>db/db</italic> transcriptomes to the normal <italic>db/+</italic> state.</p> <p>Genes that exhibited significant differential expression across the three groups were identified using the SAM statistical algorithm (selected genes are listed in ##TAB##0##
Tables 1
##
<bold>, </bold>
##TAB##1##
2
##; full gene txt files are included as <bold>Supplemental ##SUPPL##0##Gene Lists S1##, ##SUPPL##1##S2##</bold>). Interesting genes that were significantly upregulated in the untreated <italic>db/db</italic> group compared to control include PPARγ (Pparg, the receptor for rosiglitazone), FK506 binding proteins (Fkbp5 and Fkbp10), several potassium channel proteins (Kcnk1 and Kcnd3), two matrix metalloproteinases (Mmp3 and Mmp8), a cyclin-dependent kinase inhibitor (Cdkn1), and myostatin (Gdf8) (##TAB##0##
Table 1
##). Specific downregulated genes of interest include the anti-apoptotic gene B-cell lymphoma–leukemia 2 (Bcl2), several tumor necrosis factor-alpha (TNFα) genes and transforming growth factor (TGFβ) genes, vascular endothelial growth factor (Vegfc), adiponectin (Adipoq), and apelin (Apln). We confirmed the microarray results of 4 selected genes with real-time PCR (##FIG##4##
Figures 5b, 5c
##). When we compared the <italic>db/db</italic> rosiglitazone treated group to control, we found that upregulated genes include the previously-mentioned Mmp3 and Kcnk1 transcripts, as well as a potassium channel-interacting gene (Kcnip4), ryanodine receptor (Ryr1), Myosin IIIB (Myo3b), Patched homolog 1, and apelin (##TAB##1##
Table 2
##).</p> <p>While identifying individual genes that are differentially regulated in each of the conditions is informative, it is also useful to ascertain which cellular <italic>processes</italic> are up- or down- regulated; therefore, we performed the Gene Ontology pathway overrepresentation analysis using Fisher's exact test. Lipid and protein metabolism, fatty acid beta-oxidation, cell death, apoptosis, peroxisome organization, and biogenesis were significantly upregulated in untreated <italic>db/db</italic> mice when compared to control <italic>db/+</italic> mice (<bold>Supplemental ##SUPPL##0##Tables S1##, ##SUPPL##1##S2##, ##SUPPL##2##S3##, ##SUPPL##3##S4##</bold>). Among major pathways that were significantly downregulated are cell proliferation and cycle, immune response, blood vessel and vasculature development, and anti-apoptosis. Interestingly, comparing the rosiglitazone-treated <italic>db/db</italic> mice to control mice, we did not see a reversal of the pathways that were upregulated in the untreated <italic>db/db</italic> group. The additional pathways that were significantly upregulated in the rosiglitazone-treated <italic>db/db</italic> hearts, however, included secretion, exocytosis, protein and intracellular transport, cellular protein metabolism, and sphingolipid metabolism.</p>
<title>Discussion</title> <p>Few studies have studied the global cardiac gene expression changes resulting from rosiglitazone treatment. Understanding these molecular changes will help elucidate the potential mechanism(s) for the increased cardiac disease observed in diabetic patients treated with this drug. Our microarray results demonstrate that rosiglitazone does not reverse many of the significant gene expression and pathway changes that occur in untreated diabetic hearts, such as apoptosis and lipid metabolism. Furthermore, our microarray analysis indicates that rosiglitazone did not upregulate important insulin-regulated glucose transporters such as Glut4 in the heart, a finding that has been reported previously in isolated rat cardiomyocytes ##REF##12600885##[23]##. However, we observed a number of unique genes which appear to be differentially regulated between rosiglitazone-treated and untreated animals that may give mechanistic insights into rosiglitazone's effects on the diabetic heart.</p> <p>We were initially puzzled to find that PPARγ expression was upregulated in the hearts of untreated diabetic mice, as it has been reported to be low in cardiomyocytes ##REF##12649261##[12]##, ##REF##9794462##[24]##, ##REF##10696077##[25]##. However, several studies have shown PPARγ expression to be important to cardiac metabolism. One study found that rosiglitazone-treated Zucker fatty rat hearts had improved glucose uptake as well as improved contractility compared to lean rat hearts ##REF##11916933##[19]##. Another study found that mice lacking PPARγ in cardiomyocytes exhibited mild cardiac hypertrophy ##REF##16051889##[26]##. These findings may reflect an indirect effect of rosiligitazone on the heart, as demonstrated by a study in which metabolic changes in PPARγ-treated <italic>db/db</italic> hearts were thought to be secondary to changes in the supply of exogenous substrates ##REF##14600074##[21]##. It appears that cardiac PPARγ expression may play an important though poorly understood role in cardiac physiology, and it is possible that modulation of PPARγ signaling may affect the response to cardiac ischemia.</p> <p>Matrix metalloproteinases are known to play a critical role in atherosclerosis and cardiovascular tissue remodeling, mediating the balance between matrix accumulation and degradation ##REF##15341739##[27]##. We found that rosiglitazone caused significant upregulation of Mmp3, which encodes the matrix metalloproteinase stromelysin. Mmp3 is typically expressed by macrophages within the atherosclerotic plaques found in coronary artery disease ##REF##1896464##[28]##, and has been shown to promote plaque rupture, myocardial infarction and aneurysm ##REF##12821236##[29]##. It is therefore possible that the <italic>in vivo</italic> upregulation of Mmp3 by rosiglitazone will promote plaque rupture, leading to increased rates of myocardial infarction, as well as increasing tissue remodeling after ischemia.</p> <p>In both untreated and rosiglitazone-treated <italic>db/db</italic> mice, we found several potassium channel-related genes that were highly upregulated, including Kcnip4 and Kcnk1. These findings support the hypothesis that cardiomyocytes from <italic>db/db</italic> hearts exhibit electrophysiological alterations such as attenuated outward K<sup>+</sup> currents ##REF##14600074##[21]##, ##REF##11711563##[30]##. In general, glycolytic ATP production is important for the normal function of cardiac membrane ion channels and pumps, and inhibition of glycolysis causes altered intracellular ion concentrations that may affect cardiac action potentials. Specifically, Kcnip4 modulates A-type potassium channels and has been shown to increase expression of cardiac A-type Kv-4 channels ##REF##8831489##[31]##. Kcnk1 encodes a two-pore potassium channel, TWIK-1, which plays a major role in setting the resting membrane potential in many cell types ##UREF##1##[32]##. Although it is unclear to what degree cardiac arrhythmias might have played a role in the mortality shown in the meta-analysis of rosiglitazone trials ##REF##17517853##[4]##, the changes we observed in cardiomyocyte potassium channel expression raise some concerns.</p> <p>Rosiglitazone treatment also upregulated four genes related to sphingolipid metabolism: B4galnt1, Cel, Asah2, and St3gal5. Found in eukaryotic plasma membranes, sphingolipids have emerged as a class of lipid mediators believed to be important for endogenous modulation of cardiovascular function ##REF##16555682##[33]##, ##REF##11717151##[34]##. Generation of sphingolipid metabolites such as ceramide, sphingosine, and sphingosine-1-phosphate can activate smooth muscle cell proliferation, endothelial cell differentiation, apoptosis and cell death, migration, vasoconstriction, and dilation. Specifically, ceramide is an important mediator of lipotoxicity in the heart, and accumulation of this metabolite has been associated with cardiomyocyte apoptosis ##REF##12531555##[35]##. Recent evidence indicates that the TZD pioglitazone causes increased <italic>de novo</italic> ceramide synthesis in rat hearts ##REF##17259076##[36]##, although another group has reported decreased ceramide accumulation and subsequent reduction in apoptosis ##REF##10677535##[37]##. Regardless of the specific mechanism involved, our microarray data suggest a significant role for sphingolipid metabolism in the rosiglitazone-treated diabetic heart which we are actively investigating further.</p> <p>The downregulation of the hormone adiponectin in the diabetic state has been well established in adipocytes and serum ##REF##15780820##[38]##, but it was not until 2005 that Pineiro <italic>et al.</italic> showed evidence that cardiomyocytes also synthesize adiponectin ##REF##16140297##[39]##. Confirming a recent report by Natarajan <italic>et al.</italic>
##REF##18287886##[40]##, our results demonstrate a reduction of adiponectin in the diabetic heart, though we did not observe restoration of its expression after rosiglitazone treatment. In contrast, the signaling peptide apelin, a recently discovered regulator of cardiac and vascular function ##REF##12963638##[41]##, was downregulated in the untreated <italic>db/db</italic> group but was restored to normal levels after rosiglitazone treatment. Apelin is expressed in the endothelium of heart, kidney, and lung, and its receptor (APJ) is expressed by myocardial cells and some vascular smooth muscle cells. Apelin has inotropic effects on cardiac contractility and has been shown to decrease systemic vascular resistance, and is therefore likely beneficial during heart failure. In fact, apelin signaling may augment inhibition of the renin-angiotensin system ##REF##16919293##[42]##, which is known to exacerbate heart failure when left unchecked. Though not well understood, the altered expression of apelin we observed in the diabetic heart in response to rosiglitazone may present a novel pathway for the study of diabetic cardiomyopathy.</p> <p>In conclusion, this study has mapped the genomic expression patterns in the hearts of the <italic>db/db</italic> murine model of diabetes and the effects of rosiglitazone on these patterns. Overall, as we expected, the <italic>db/db</italic> gene expression signature was markedly different from control, but to our initial surprise was not significantly reversed with rosiglitazone. In fact, many of the transcriptional changes induced by diabetic disease on the heart appeared to be further exacerbated by rosiglitazone. In particular, we have highlighted a number of rosiglitazone modulated genes and pathways that may play a role in the pathophysiology of the increase in cardiac mortality seen in the meta-analysis by Nissen <italic>et al.</italic> The cumulative upregulation of (1) a matrix metalloproteinase gene that has previously been implicated in plaque rupture, (2) potassium channel genes involved in membrane potential maintenance and action potential generation, and (3) sphingolipid and ceramide metabolism related genes, all give cause for concern over rosiglitazone's safety. Interestingly, a second meta-analysis has suggested that pioglitazone, also a TZD, may not have the same negative cardiovascular effects as rosiglitazone ##REF##17848652##[43]##. We believe the overall expression changes caused by rosiglitazone in the heart will likely also be found with pioglitazone treatment, though with some differences. This is based on previous microarray studies of adipocytes ##REF##17963725##[44]## and hepatocytes ##REF##17031537##[45]## that found similar, but not identical, expression profiles after treatment with pioglitazone, rosiglitazone, and other TZDs. Future investigations of how these two drugs differ at the transcriptional level in the heart will further elucidate the networks that govern diabetic cardiovascular function and disease.</p>
<p>Conceived and designed the experiments: PT JW. Performed the experiments: KW ZL PY GN MH. Analyzed the data: TQ RW KW JW ZL PY PGY. Contributed reagents/materials/analysis tools: TQ PT PY GN DH. Wrote the paper: RW KW JW. Other: The Human Immune Monitoring Center at Stanford University facility played a role in the design of the experiments and discussion of the results. The Human Immune Monitoring Center at Stanford University reviewed the final draft of the paper. Director of the Human Immune Monitoring Center at Stanford University that supervised a technician that conducted all the Agilent microarray experiments and data analysis: DH.</p> <title>Background</title> <p>Recently, the type 2 diabetes medication, rosiglitazone, has come under scrutiny for possibly increasing the risk of cardiac disease and death. To investigate the effects of rosiglitazone on the diabetic heart, we performed cardiac transcriptional profiling and imaging studies of a murine model of type 2 diabetes, the C57BL/KLS-lepr<sup>db</sup>/lepr<sup>db</sup> (<italic>db/db</italic>) mouse.</p> <title>Methods and Findings</title> <p>We compared cardiac gene expression profiles from three groups: untreated <italic>db/db</italic> mice, <italic>db/db</italic> mice after rosiglitazone treatment, and non-diabetic <italic>db/+</italic> mice. Prior to sacrifice, we also performed cardiac magnetic resonance (CMR) and echocardiography. As expected, overall the <italic>db/db</italic> gene expression signature was markedly different from control, but to our surprise was not significantly reversed with rosiglitazone. In particular, we have uncovered a number of rosiglitazone modulated genes and pathways that may play a role in the pathophysiology of the increase in cardiac mortality as seen in several recent meta-analyses. Specifically, the cumulative upregulation of (1) a matrix metalloproteinase gene that has previously been implicated in plaque rupture, (2) potassium channel genes involved in membrane potential maintenance and action potential generation, and (3) sphingolipid and ceramide metabolism-related genes, together give cause for concern over rosiglitazone's safety. Lastly, <italic>in vivo</italic> imaging studies revealed minimal differences between rosiglitazone-treated and untreated <italic>db/db</italic> mouse hearts, indicating that rosiglitazone's effects on gene expression in the heart do not immediately turn into detectable gross functional changes.</p> <title>Conclusions</title> <p>This study maps the genomic expression patterns in the hearts of the <italic>db/db</italic> murine model of diabetes and illustrates the impact of rosiglitazone on these patterns. The <italic>db/db</italic> gene expression signature was markedly different from control, and was not reversed with rosiglitazone. A smaller number of unique and interesting changes in gene expression were noted with rosiglitazone treatment. Further study of these genes and molecular pathways will provide important insights into the cardiac decompensation associated with both diabetes and rosiglitazone treatment.</p>
<title>Supporting Information</title>
<fig id="pone-0002609-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.g001</object-id><label>Figure 1</label><caption><title>Insulin resistance and mean body/heart weights of rosiglitazone-treated, untreated, and control mice.</title><p>(A) Plasma NEFA levels at 0, 1, and 3 months in untreated <italic>db/db</italic> mice (n = 10), rosiglitazone-treated <italic>db/db</italic> mice (n = 10), and <italic>db/+</italic> mice (n = 10). NEFA levels were higher in both the treated and untreated <italic>db/db</italic> groups at baseline when compared to <italic>db/+</italic> control mice. Compared to untreated <italic>db/db</italic> mice, NEFA levels in the rosiglitazone-treated <italic>db/db</italic> group decreased significantly during the first month (960.9 vs. 1501.0 mmol/L, <italic>P</italic>&lt;0.05), and decreased further by the third month to near <italic>db/+</italic> control levels. (B) Fasting insulin levels in the three groups at 0, 1, and 3 months. At baseline, both treated and untreated <italic>db/db</italic> groups had increased levels of insulin compared to <italic>db/+</italic> control. Insulin levels dramatically increased in untreated <italic>db/db</italic> mice at 1 month, but progressively decreased in rosiglitazone-treated mice. (C) Insulin tolerance testing at 3 months. Serum glucose levels after insulin injection were significantly higher in untreated <italic>db/db</italic> mice when compared to <italic>db/+</italic> mice; the rosiglitazone-treated <italic>db/db</italic> group had moderately elevated glucose levels after insulin administration. These results corroborate the NEFA and insulin studies that demonstrate improved insulin sensitivity with rosiglitazone treatment. (D) Mean body weights of the three groups of mice at 0, 1, and 3 months. Weights of both treated and untreated <italic>db/db</italic> mice were markedly higher than <italic>db/+</italic> controls at baseline. Furthermore, rosiglitazone-treated mice had generally higher mean body weights when compared to untreated mice. (E) Mean heart/body weight ratios at sacrifice (4 months after treatment initiation). No significant difference between treated and untreated groups (3.27 vs. 2.58, <italic>P</italic> = 0.13), though there does appear to be a downward trend in the data. Heart/body weight ratio was significantly higher in the <italic>db/+</italic> group compared to both <italic>db/db</italic> groups. Values are mean±SEM. *<italic>P</italic>&lt;0.05 vs. age-matched <italic>db/+</italic>; #<italic>P</italic>&lt;0.05 vs. age-matched untreated <italic>db/db</italic>.</p></caption></fig> <fig id="pone-0002609-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.g002</object-id><label>Figure 2</label><caption><title>High-power light microscopy slides of myocardial tissue from rosiglitazone-treated and untreated <italic>db/db</italic> mice, and <italic>db/+</italic> controls.</title><p>Individual cardiomyocytes from all three groups appear qualitatively similar based on morphology and size, though the untreated <italic>db/db</italic> cardiomyocytes exhibit very mild thickening. No increases in inflammatory cells, cell death, fibrosis, or other processes were observed in the two <italic>db/db</italic> groups compared to <italic>db/+</italic> controls.</p></caption></fig> <fig id="pone-0002609-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.g003</object-id><label>Figure 3</label><caption><title>Left ventricular functional analysis with echocardiogram.</title><p>The mean left ventricular ejection fraction (LVEF) and fractional shortening (FS) of the untreated <italic>db/db</italic> group (72.4±2.3% and 36.2±1.7%, respectively) were not significantly different from the <italic>db/+</italic> group (76.9±4.3%, <italic>P</italic> = 0.11 and 39.7±3.6%, <italic>P</italic> = 0.11, respectively), though there appears to be a downward trend in contractility between these two groups. Further, there were no significant LVEF or FS differences between the rosiglitazone-treated <italic>db/db</italic> group (74.3±3.9% and 37.8±3.1%, respectively) and untreated <italic>db/db</italic> group (<italic>P</italic> = 0.43 and 0.40, respectively). Values are mean±SEM.</p></caption></fig> <fig id="pone-0002609-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.g004</object-id><label>Figure 4</label><caption><title>Cardiac magnetic resonance imaging.</title><p>CMR imaging at 8 weeks after treatment initiation of rosiglitazone-treated <italic>db/db</italic> mice (n = 4), untreated <italic>db/db</italic> mice (n = 4) and <italic>db/+</italic> groups (n = 4). Left ventricular mass (LVM), interventricular septal thickness (IVST) and posterior wall thickness (PWT) were significantly increased in both treated and untreated <italic>db/db</italic> mice relative to <italic>db/+</italic> controls, while LVEF was decreased (though not significantly in the rosiglitazone-treated group). No significant differences in LVM, IVST, PWT, or LVEF between rosiglitazone-treated and untreated <italic>db/db</italic> mice were detected (<italic>P</italic> = 0.82, 0.66, 0.14, and 0.88, respectively). Values are mean±SEM. *<italic>P</italic>&lt;0.05 vs. age-matched <italic>db/+</italic>.</p></caption></fig> <fig id="pone-0002609-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.g005</object-id><label>Figure 5</label><caption><title>Multidimensional scaling plot of the expression data and real-time PCR analysis.</title><p>(A) 3-dimensional scaling plot provides a graphical representation of high-dimensional expression data in low dimensions. Each point within a “cloud” represents a single microarray, and the similarity within a set of microarrays is indicated by their physical proximity to one other. As evidenced from the figure, each group (<italic>db/+</italic> control, untreated <italic>db/db</italic>, and rosiglitazone-treated <italic>db/db</italic>) clusters into a distinct grouping, indicating that each group has a similar transcriptome that can be distinguished from the two other groups. Taqman real-time PCR of four selected genes normalized to 18s shows (B) significant upregulation of Kcnk1 in untreated <italic>db/db</italic> vs. <italic>db/+</italic> mice and (C) significant upregulation of all four genes in rosiglitazone-treated vs. untreated <italic>db/db</italic> mice. These data confirm the altered regulatory patterns of these four genes in the microarray data. Values are mean±SEM. *<italic>P</italic>&lt;0.05 vs. control.</p></caption></fig>
<table-wrap id="pone-0002609-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.t001</object-id><label>Table 1</label><caption><title>Untreated <italic>db/db</italic> vs. <italic>db/+</italic> control mice.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene Symbol</td><td align="left" rowspan="1" colspan="1">RefSeq ID</td><td align="left" rowspan="1" colspan="1">UGCluster</td><td align="left" rowspan="1" colspan="1">Name</td><td align="left" rowspan="1" colspan="1">Score(d)</td><td align="left" rowspan="1" colspan="1">Fold Change</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Kcnk1</td><td align="left" rowspan="1" colspan="1">NM_008430</td><td align="left" rowspan="1" colspan="1">Mm.10800</td><td align="left" rowspan="1" colspan="1">Potassium channel, subfamily K, member 1</td><td align="left" rowspan="1" colspan="1">7.82</td><td align="left" rowspan="1" colspan="1">26.03</td></tr><tr><td align="left" rowspan="1" colspan="1">Gdf8</td><td align="left" rowspan="1" colspan="1">NM_010834</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">growth differentiation factor 8 (Gdf8)</td><td align="left" rowspan="1" colspan="1">5.12</td><td align="left" rowspan="1" colspan="1">15.93</td></tr><tr><td align="left" rowspan="1" colspan="1">Egfbp2</td><td align="left" rowspan="1" colspan="1">NM_010115</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">epidermal growth factor binding protein type B (Egfbp2)</td><td align="left" rowspan="1" colspan="1">4.77</td><td align="left" rowspan="1" colspan="1">9.47</td></tr><tr><td align="left" rowspan="1" colspan="1">Fkbp5</td><td align="left" rowspan="1" colspan="1">NM_010220</td><td align="left" rowspan="1" colspan="1">Mm.276405</td><td align="left" rowspan="1" colspan="1">FK506 binding protein 5</td><td align="left" rowspan="1" colspan="1">9.76</td><td align="left" rowspan="1" colspan="1">6.11</td></tr><tr><td align="left" rowspan="1" colspan="1">Mmp8</td><td align="left" rowspan="1" colspan="1">NM_008611</td><td align="left" rowspan="1" colspan="1">Mm.16415</td><td align="left" rowspan="1" colspan="1">Matrix metallopeptidase 8</td><td align="left" rowspan="1" colspan="1">5.22</td><td align="left" rowspan="1" colspan="1">5.93</td></tr><tr><td align="left" rowspan="1" colspan="1">Edn3</td><td align="left" rowspan="1" colspan="1">NM_007903</td><td align="left" rowspan="1" colspan="1">Mm.9478</td><td align="left" rowspan="1" colspan="1">Endothelin 3</td><td align="left" rowspan="1" colspan="1">6.59</td><td align="left" rowspan="1" colspan="1">5.72</td></tr><tr><td align="left" rowspan="1" colspan="1">Cdkn1a</td><td align="left" rowspan="1" colspan="1">NM_007669</td><td align="left" rowspan="1" colspan="1">Mm.195663</td><td align="left" rowspan="1" colspan="1">Cyclin-dependent kinase inhibitor 1A (P21)</td><td align="left" rowspan="1" colspan="1">5.13</td><td align="left" rowspan="1" colspan="1">3.74</td></tr><tr><td align="left" rowspan="1" colspan="1">Kcnj4</td><td align="left" rowspan="1" colspan="1">NM_008427</td><td align="left" rowspan="1" colspan="1">Mm.140760</td><td align="left" rowspan="1" colspan="1">Potassium inwardly-rectifying channel, subfamily J, member 4</td><td align="left" rowspan="1" colspan="1">5.06</td><td align="left" rowspan="1" colspan="1">2.68</td></tr><tr><td align="left" rowspan="1" colspan="1">Map3k6</td><td align="left" rowspan="1" colspan="1">NM_016693</td><td align="left" rowspan="1" colspan="1">Mm.36640</td><td align="left" rowspan="1" colspan="1">Mitogen-activated protein kinase kinase kinase 6</td><td align="left" rowspan="1" colspan="1">9.39</td><td align="left" rowspan="1" colspan="1">2.54</td></tr><tr><td align="left" rowspan="1" colspan="1">Map3k15</td><td align="left" rowspan="1" colspan="1">BC031147</td><td align="left" rowspan="1" colspan="1">Mm.386889</td><td align="left" rowspan="1" colspan="1">Mitogen-activated protein kinase kinase kinase 15</td><td align="left" rowspan="1" colspan="1">5.88</td><td align="left" rowspan="1" colspan="1">2.10</td></tr><tr><td align="left" rowspan="1" colspan="1">Fkbp10</td><td align="left" rowspan="1" colspan="1">NM_010221</td><td align="left" rowspan="1" colspan="1">Mm.3894</td><td align="left" rowspan="1" colspan="1">FK506 binding protein 10</td><td align="left" rowspan="1" colspan="1">5.21</td><td align="left" rowspan="1" colspan="1">2.09</td></tr><tr><td align="left" rowspan="1" colspan="1">Mmp3</td><td align="left" rowspan="1" colspan="1">NM_010809</td><td align="left" rowspan="1" colspan="1">Mm.4993</td><td align="left" rowspan="1" colspan="1">Matrix metallopeptidase 3</td><td align="left" rowspan="1" colspan="1">4.30</td><td align="left" rowspan="1" colspan="1">1.90</td></tr><tr><td align="left" rowspan="1" colspan="1">Kcnd3</td><td align="left" rowspan="1" colspan="1">NM_019931</td><td align="left" rowspan="1" colspan="1">Mm.44530</td><td align="left" rowspan="1" colspan="1">Potassium voltage-gated channel, Shal-related family, member 3</td><td align="left" rowspan="1" colspan="1">5.93</td><td align="left" rowspan="1" colspan="1">1.77</td></tr><tr><td align="left" rowspan="1" colspan="1">Ryr3</td><td align="left" rowspan="1" colspan="1">XM_619795</td><td align="left" rowspan="1" colspan="1">Mm.436657</td><td align="left" rowspan="1" colspan="1">Ryanodine receptor 3</td><td align="left" rowspan="1" colspan="1">5.45</td><td align="left" rowspan="1" colspan="1">1.74</td></tr><tr><td align="left" rowspan="1" colspan="1">Pparg</td><td align="left" rowspan="1" colspan="1">NM_011146</td><td align="left" rowspan="1" colspan="1">Mm.3020</td><td align="left" rowspan="1" colspan="1">Peroxisome proliferator activated receptor gamma</td><td align="left" rowspan="1" colspan="1">5.51</td><td align="left" rowspan="1" colspan="1">1.42</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfrsf22</td><td align="left" rowspan="1" colspan="1">NM_023680</td><td align="left" rowspan="1" colspan="1">Mm.261384</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor receptor superfamily, member 22</td><td align="left" rowspan="1" colspan="1">4.06</td><td align="left" rowspan="1" colspan="1">1.40</td></tr><tr><td align="left" rowspan="1" colspan="1">Aifm2</td><td align="left" rowspan="1" colspan="1">NM_153779</td><td align="left" rowspan="1" colspan="1">Mm.286309</td><td align="left" rowspan="1" colspan="1">Apoptosis-inducing factor, mitochondrion-associated 2</td><td align="left" rowspan="1" colspan="1">5.62</td><td align="left" rowspan="1" colspan="1">1.39</td></tr><tr><td align="left" rowspan="1" colspan="1">Il4ra</td><td align="left" rowspan="1" colspan="1">NM_001008700</td><td align="left" rowspan="1" colspan="1">Mm.233802|Mm.441865</td><td align="left" rowspan="1" colspan="1">interleukin 4 receptor, alpha (Il4ra)</td><td align="left" rowspan="1" colspan="1">4.96</td><td align="left" rowspan="1" colspan="1">1.37</td></tr><tr><td align="left" rowspan="1" colspan="1">Sod1</td><td align="left" rowspan="1" colspan="1">BC057074</td><td align="left" rowspan="1" colspan="1">Mm.431677</td><td align="left" rowspan="1" colspan="1">superoxide dismutase 1, soluble (cDNA clone IMAGE:5697222), partial cds.</td><td align="left" rowspan="1" colspan="1">4.02</td><td align="left" rowspan="1" colspan="1">1.34</td></tr><tr><td align="left" rowspan="1" colspan="1">Adipor2</td><td align="left" rowspan="1" colspan="1">NM_197985</td><td align="left" rowspan="1" colspan="1">Mm.291826</td><td align="left" rowspan="1" colspan="1">Adiponectin receptor 2</td><td align="left" rowspan="1" colspan="1">4.42</td><td align="left" rowspan="1" colspan="1">1.22</td></tr><tr><td align="left" rowspan="1" colspan="1">Myo9a</td><td align="left" rowspan="1" colspan="1">AK029836</td><td align="left" rowspan="1" colspan="1">Mm.249545</td><td align="left" rowspan="1" colspan="1">Myosin IXa</td><td align="left" rowspan="1" colspan="1">−4.70</td><td align="left" rowspan="1" colspan="1">−1.26</td></tr><tr><td align="left" rowspan="1" colspan="1">Casp7</td><td align="left" rowspan="1" colspan="1">NM_007611</td><td align="left" rowspan="1" colspan="1">Mm.35687</td><td align="left" rowspan="1" colspan="1">Caspase 7</td><td align="left" rowspan="1" colspan="1">−4.46</td><td align="left" rowspan="1" colspan="1">−1.27</td></tr><tr><td align="left" rowspan="1" colspan="1">Casp2</td><td align="left" rowspan="1" colspan="1">NM_007610</td><td align="left" rowspan="1" colspan="1">Mm.3921|Mm.433648</td><td align="left" rowspan="1" colspan="1">caspase 2 (Casp2)</td><td align="left" rowspan="1" colspan="1">−5.10</td><td align="left" rowspan="1" colspan="1">−1.34</td></tr><tr><td align="left" rowspan="1" colspan="1">Vegfc</td><td align="left" rowspan="1" colspan="1">NM_009506</td><td align="left" rowspan="1" colspan="1">Mm.1402</td><td align="left" rowspan="1" colspan="1">Vascular endothelial growth factor C</td><td align="left" rowspan="1" colspan="1">−6.20</td><td align="left" rowspan="1" colspan="1">−1.39</td></tr><tr><td align="left" rowspan="1" colspan="1">Il13ra1</td><td align="left" rowspan="1" colspan="1">NM_133990</td><td align="left" rowspan="1" colspan="1">Mm.24208</td><td align="left" rowspan="1" colspan="1">Interleukin 13 receptor, alpha 1</td><td align="left" rowspan="1" colspan="1">−4.24</td><td align="left" rowspan="1" colspan="1">−1.40</td></tr><tr><td align="left" rowspan="1" colspan="1">Gja1</td><td align="left" rowspan="1" colspan="1">NM_010288</td><td align="left" rowspan="1" colspan="1">Mm.378921</td><td align="left" rowspan="1" colspan="1">Gap junction membrane channel protein alpha 1</td><td align="left" rowspan="1" colspan="1">−4.85</td><td align="left" rowspan="1" colspan="1">−1.49</td></tr><tr><td align="left" rowspan="1" colspan="1">Il13ra1</td><td align="left" rowspan="1" colspan="1">NM_133990</td><td align="left" rowspan="1" colspan="1">Mm.24208</td><td align="left" rowspan="1" colspan="1">Interleukin 13 receptor, alpha 1</td><td align="left" rowspan="1" colspan="1">−4.23</td><td align="left" rowspan="1" colspan="1">−1.50</td></tr><tr><td align="left" rowspan="1" colspan="1">Myo1b</td><td align="left" rowspan="1" colspan="1">NM_010863</td><td align="left" rowspan="1" colspan="1">Mm.3390</td><td align="left" rowspan="1" colspan="1">Myosin IB</td><td align="left" rowspan="1" colspan="1">−6.11</td><td align="left" rowspan="1" colspan="1">−1.57</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfaip8</td><td align="left" rowspan="1" colspan="1">NM_134131</td><td align="left" rowspan="1" colspan="1">Mm.27740</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor, alpha-induced protein 8</td><td align="left" rowspan="1" colspan="1">−4.65</td><td align="left" rowspan="1" colspan="1">−1.71</td></tr><tr><td align="left" rowspan="1" colspan="1">Adrb2</td><td align="left" rowspan="1" colspan="1">NM_007420</td><td align="left" rowspan="1" colspan="1">Mm.5598</td><td align="left" rowspan="1" colspan="1">Adrenergic receptor, beta 2</td><td align="left" rowspan="1" colspan="1">−4.91</td><td align="left" rowspan="1" colspan="1">−1.81</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfaip8l1</td><td align="left" rowspan="1" colspan="1">NM_025566</td><td align="left" rowspan="1" colspan="1">Mm.2312</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor, alpha-induced protein 8-like 1</td><td align="left" rowspan="1" colspan="1">−5.52</td><td align="left" rowspan="1" colspan="1">−1.88</td></tr><tr><td align="left" rowspan="1" colspan="1">Il16</td><td align="left" rowspan="1" colspan="1">NM_010551</td><td align="left" rowspan="1" colspan="1">Mm.10137</td><td align="left" rowspan="1" colspan="1">Interleukin 16</td><td align="left" rowspan="1" colspan="1">−24.20</td><td align="left" rowspan="1" colspan="1">−2.42</td></tr><tr><td align="left" rowspan="1" colspan="1">Myo1d</td><td align="left" rowspan="1" colspan="1">AK037051</td><td align="left" rowspan="1" colspan="1">Mm.151948</td><td align="left" rowspan="1" colspan="1">Myosin ID</td><td align="left" rowspan="1" colspan="1">−4.19</td><td align="left" rowspan="1" colspan="1">−2.45</td></tr><tr><td align="left" rowspan="1" colspan="1">Apln</td><td align="left" rowspan="1" colspan="1">NM_013912</td><td align="left" rowspan="1" colspan="1">Mm.29262</td><td align="left" rowspan="1" colspan="1">Apelin</td><td align="left" rowspan="1" colspan="1">−4.51</td><td align="left" rowspan="1" colspan="1">−2.54</td></tr><tr><td align="left" rowspan="1" colspan="1">Bcl2</td><td align="left" rowspan="1" colspan="1">NM_009741</td><td align="left" rowspan="1" colspan="1">257460|Mm.462969</td><td align="left" rowspan="1" colspan="1">B-cell leukemia/lymphoma 2 (Bcl2), transcript variant 1</td><td align="left" rowspan="1" colspan="1">−7.36</td><td align="left" rowspan="1" colspan="1">−2.56</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfrsf13c</td><td align="left" rowspan="1" colspan="1">NM_028075</td><td align="left" rowspan="1" colspan="1">Mm.240047</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor receptor superfamily, member 13c</td><td align="left" rowspan="1" colspan="1">−6.33</td><td align="left" rowspan="1" colspan="1">−4.54</td></tr><tr><td align="left" rowspan="1" colspan="1">Adipoq</td><td align="left" rowspan="1" colspan="1">NM_009605</td><td align="left" rowspan="1" colspan="1">Mm.3969</td><td align="left" rowspan="1" colspan="1">Adiponectin, C1Q and collagen domain containing</td><td align="left" rowspan="1" colspan="1">−5.31</td><td align="left" rowspan="1" colspan="1">−4.59</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002609-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002609.t002</object-id><label>Table 2</label><caption><title>Rosiglitazone-treated vs. untreated <italic>db/db</italic> mice.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene Symbol</td><td align="left" rowspan="1" colspan="1">RefSeq ID</td><td align="left" rowspan="1" colspan="1">UGCluster</td><td align="left" rowspan="1" colspan="1">Name</td><td align="left" rowspan="1" colspan="1">Score(d)</td><td align="left" rowspan="1" colspan="1">Fold Change</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Kcnip4</td><td align="left" rowspan="1" colspan="1">NM_030265</td><td align="left" rowspan="1" colspan="1">Mm.160172</td><td align="left" rowspan="1" colspan="1">Kv channel interacting protein 4</td><td align="left" rowspan="1" colspan="1">4.46</td><td align="left" rowspan="1" colspan="1">4.76</td></tr><tr><td align="left" rowspan="1" colspan="1">Ptch1</td><td align="left" rowspan="1" colspan="1">NM_008957</td><td align="left" rowspan="1" colspan="1">Mm.228798</td><td align="left" rowspan="1" colspan="1">Patched homolog 1</td><td align="left" rowspan="1" colspan="1">6.67</td><td align="left" rowspan="1" colspan="1">3.32</td></tr><tr><td align="left" rowspan="1" colspan="1">Ryr1</td><td align="left" rowspan="1" colspan="1">NM_009109</td><td align="left" rowspan="1" colspan="1">Mm.439745</td><td align="left" rowspan="1" colspan="1">Ryanodine receptor 1, skeletal muscle</td><td align="left" rowspan="1" colspan="1">5.27</td><td align="left" rowspan="1" colspan="1">3.17</td></tr><tr><td align="left" rowspan="1" colspan="1">Kcnk1</td><td align="left" rowspan="1" colspan="1">NM_008430</td><td align="left" rowspan="1" colspan="1">Mm.10800</td><td align="left" rowspan="1" colspan="1">Potassium channel, subfamily K, member 1</td><td align="left" rowspan="1" colspan="1">6.38</td><td align="left" rowspan="1" colspan="1">2.84</td></tr><tr><td align="left" rowspan="1" colspan="1">Myo3b</td><td align="left" rowspan="1" colspan="1">AK033795</td><td align="left" rowspan="1" colspan="1">Mm.99648|Mm.458853</td><td align="left" rowspan="1" colspan="1">adult male epididymis cDNA, RIKEN full-length enriched library, clone:9230110G05</td><td align="left" rowspan="1" colspan="1">4.46</td><td align="left" rowspan="1" colspan="1">2.58</td></tr><tr><td align="left" rowspan="1" colspan="1">Mmp3</td><td align="left" rowspan="1" colspan="1">NM_010809</td><td align="left" rowspan="1" colspan="1">Mm.4993</td><td align="left" rowspan="1" colspan="1">Matrix metallopeptidase 3</td><td align="left" rowspan="1" colspan="1">4.70</td><td align="left" rowspan="1" colspan="1">2.03</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfsf18</td><td align="left" rowspan="1" colspan="1">NM_183391</td><td align="left" rowspan="1" colspan="1">Mm.276823</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor (ligand) superfamily, member 18</td><td align="left" rowspan="1" colspan="1">3.82</td><td align="left" rowspan="1" colspan="1">1.86</td></tr><tr><td align="left" rowspan="1" colspan="1">Apln</td><td align="left" rowspan="1" colspan="1">NM_013912</td><td align="left" rowspan="1" colspan="1">Mm.29262</td><td align="left" rowspan="1" colspan="1">Apelin</td><td align="left" rowspan="1" colspan="1">3.52</td><td align="left" rowspan="1" colspan="1">1.80</td></tr><tr><td align="left" rowspan="1" colspan="1">Anxa13</td><td align="left" rowspan="1" colspan="1">NM_027211</td><td align="left" rowspan="1" colspan="1">Mm.237985</td><td align="left" rowspan="1" colspan="1">Annexin A13</td><td align="left" rowspan="1" colspan="1">4.04</td><td align="left" rowspan="1" colspan="1">1.64</td></tr><tr><td align="left" rowspan="1" colspan="1">St3gal5</td><td align="left" rowspan="1" colspan="1">NM_011375</td><td align="left" rowspan="1" colspan="1">Mm.38248</td><td align="left" rowspan="1" colspan="1">ST3 beta-galactoside alpha-2,3-sialyltransferase 5</td><td align="left" rowspan="1" colspan="1">5.56</td><td align="left" rowspan="1" colspan="1">1.62</td></tr><tr><td align="left" rowspan="1" colspan="1">Col22a1</td><td align="left" rowspan="1" colspan="1">XM_193814</td><td align="left" rowspan="1" colspan="1">Mm.322500</td><td align="left" rowspan="1" colspan="1">Collagen, type XXII, alpha 1</td><td align="left" rowspan="1" colspan="1">3.65</td><td align="left" rowspan="1" colspan="1">1.45</td></tr><tr><td align="left" rowspan="1" colspan="1">Col4a1</td><td align="left" rowspan="1" colspan="1">NM_009931</td><td align="left" rowspan="1" colspan="1">Mm.738</td><td align="left" rowspan="1" colspan="1">Procollagen, type IV, alpha 1</td><td align="left" rowspan="1" colspan="1">4.98</td><td align="left" rowspan="1" colspan="1">1.26</td></tr><tr><td align="left" rowspan="1" colspan="1">Kras</td><td align="left" rowspan="1" colspan="1">NM_021284</td><td align="left" rowspan="1" colspan="1">Mm.383182</td><td align="left" rowspan="1" colspan="1">V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog</td><td align="left" rowspan="1" colspan="1">4.40</td><td align="left" rowspan="1" colspan="1">1.23</td></tr><tr><td align="left" rowspan="1" colspan="1">B4galnt1</td><td align="left" rowspan="1" colspan="1">BC022180</td><td align="left" rowspan="1" colspan="1">Mm.386762</td><td align="left" rowspan="1" colspan="1">Beta-1,4-N-acetyl-galactosaminyl transferase 1</td><td align="left" rowspan="1" colspan="1">7.83</td><td align="left" rowspan="1" colspan="1">1.44</td></tr><tr><td align="left" rowspan="1" colspan="1">Cel</td><td align="left" rowspan="1" colspan="1">NM_009885</td><td align="left" rowspan="1" colspan="1">Mm.236017</td><td align="left" rowspan="1" colspan="1">Carboxyl ester lipase</td><td align="left" rowspan="1" colspan="1">6.31</td><td align="left" rowspan="1" colspan="1">1.38</td></tr><tr><td align="left" rowspan="1" colspan="1">Pparbp</td><td align="left" rowspan="1" colspan="1">NM_013634</td><td align="left" rowspan="1" colspan="1">Mm.12926</td><td align="left" rowspan="1" colspan="1">Peroxisome proliferator activated receptor binding protein</td><td align="left" rowspan="1" colspan="1">3.68</td><td align="left" rowspan="1" colspan="1">1.29</td></tr><tr><td align="left" rowspan="1" colspan="1">Tnfrsf10b</td><td align="left" rowspan="1" colspan="1">NM_020275</td><td align="left" rowspan="1" colspan="1">Mm.193430</td><td align="left" rowspan="1" colspan="1">Tumor necrosis factor receptor superfamily, member 10b</td><td align="left" rowspan="1" colspan="1">3.54</td><td align="left" rowspan="1" colspan="1">1.27</td></tr><tr><td align="left" rowspan="1" colspan="1">Asah2</td><td align="left" rowspan="1" colspan="1">NM_018830</td><td align="left" rowspan="1" colspan="1">Mm.104900</td><td align="left" rowspan="1" colspan="1">N-acylsphingosine amidohydrolase 2</td><td align="left" rowspan="1" colspan="1">4.38</td><td align="left" rowspan="1" colspan="1">1.21</td></tr><tr><td align="left" rowspan="1" colspan="1">Myh9</td><td align="left" rowspan="1" colspan="1">NM_022410</td><td align="left" rowspan="1" colspan="1">Mm.29677</td><td align="left" rowspan="1" colspan="1">Myosin, heavy polypeptide 9, non-muscle</td><td align="left" rowspan="1" colspan="1">4.03</td><td align="left" rowspan="1" colspan="1">1.19</td></tr><tr><td align="left" rowspan="1" colspan="1">Bdkrb1</td><td align="left" rowspan="1" colspan="1">NM_007539</td><td align="left" rowspan="1" colspan="1">Mm.377078</td><td align="left" rowspan="1" colspan="1">Bradykinin receptor, beta 1</td><td align="left" rowspan="1" colspan="1">−6.72</td><td align="left" rowspan="1" colspan="1">−1.96</td></tr><tr><td align="left" rowspan="1" colspan="1">Hsf1</td><td align="left" rowspan="1" colspan="1">NM_008296</td><td align="left" rowspan="1" colspan="1">Mm.347444</td><td align="left" rowspan="1" colspan="1">Heat shock factor 1</td><td align="left" rowspan="1" colspan="1">−5.70</td><td align="left" rowspan="1" colspan="1">−1.17</td></tr><tr><td align="left" rowspan="1" colspan="1">Col23a1</td><td align="left" rowspan="1" colspan="1">AK162470</td><td align="left" rowspan="1" colspan="1">Mm.154093|Mm.392367</td><td align="left" rowspan="1" colspan="1">12 days embryo embryonic body between diaphragm region and neck cDNA, RIKEN full-length enriched library, clone:9430076L03</td><td align="left" rowspan="1" colspan="1">−5.70</td><td align="left" rowspan="1" colspan="1">−1.13</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002609.s001"><label>Table S1</label><caption><p>Upregulated Gene Ontology pathways with rosiglitazone-treated db/db mice.</p><p>(0.03 MB RTF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002609.s002"><label>Table S2</label><caption><p>Downregulated Gene Ontology pathways with rosiglitazone-treated db/db mice.</p><p>(0.02 MB RTF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002609.s003"><label>Table S3</label><caption><p>Upregulated Gene Ontology pathways in untreated db/db mice.</p><p>(0.04 MB RTF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002609.s004"><label>Table S4</label><caption><p>Downregulated Gene Ontology pathways in untreated db/db mice.</p><p>(0.05 MB RTF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002609.s005"><label>Gene List S1</label><caption><p>txt file of microarray results from untreated db/db vs. db/+ analysis</p><p>(0.15 MB TXT)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002609.s006"><label>Gene List S2</label><caption><p>txt file of microarray results from db/db+rosiglitazone vs db/db analysis.</p><p>(0.04 MB TXT)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>Selected genes that were differentially expressed across untreated <italic>db/db</italic> mice and <italic>db/+</italic> controls. All of the genes presented in the tables have statistically significant expression changes. Listed are the most important parameters, including the standard identifications of the gene, a statistical measure of significance (“Score(d)”), and magnitude of the change (“Fold Change”).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><label/><p>Selected genes that were differentially expressed across rosiglitazone-treated and untreated <italic>db/db</italic> mice.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>Funding for this work was supplied in part by research grants from Burroughs Wellcome Foundation Career Award in Medical Scientist (JCW). KDW is a recipient of a Stanford University Bio-X Graduate Student Fellowship and Society of Nuclear Medicine Student Fellowship. This study was designed by JCW in collaboration with his academic colleagues. This investigation did not have separate industry funding. The funding sources had no role in the study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.</p></fn></fn-group>
<graphic xlink:href="pone.0002609.g001"/> <graphic xlink:href="pone.0002609.g002"/> <graphic xlink:href="pone.0002609.g003"/> <graphic xlink:href="pone.0002609.g004"/> <graphic xlink:href="pone.0002609.g005"/> <graphic id="pone-0002609-t001-1" xlink:href="pone.0002609.t001"/> <graphic id="pone-0002609-t002-2" xlink:href="pone.0002609.t002"/>
<media xlink:href="pone.0002609.s001.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002609.s002.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002609.s003.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002609.s004.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002609.s005.txt"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002609.s006.txt"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Vamecq", "Latruffe"], "given-names": ["J", "N"], "year": ["1999"], "article-title": ["Medical significance of peroxisome proliferator-activated receptors."], "source": ["The Lancet"], "volume": ["354"], "fpage": ["141"], "lpage": ["148"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Lesage", "Lazdunski"], "given-names": ["F", "M"], "year": ["2000"], "article-title": ["Molecular and functional properties of two-pore-domain potassium channels."], "source": ["American Journal of Physiology Renal, Fluid and Electrolyte Physiology"], "volume": ["279"], "fpage": ["F793"], "lpage": ["F801"]}]
PLoS One. 2008 Jul 9; 3(7):e2609