Diabetes 54:1854-1862, 2005 © 2005 by the American Diabetes Association, Inc. Pharmacogenetic Analysis of Rosiglitazone-Induced Hepatosteatosis in New Mouse Models of Type 2 Diabetes
1 The Jackson Laboratory, Bar Harbor, Maine
Although thiazolidinediones suppress hyperglycemia in diabetic (NON x NZO)F1 males, these mice exhibit unusual sensitivity to drug-induced exacerbation of an underlying hepatosteatosis only rarely experienced in human patients. To establish the pharmacogenetic basis for this sensitivity, a panel of recombinant congenic strains (RCSs) with varying degrees of obesity and diabetes was generated by fixing selected NZO HlLt alleles on the diabetes- and hepatosteatosis-resistant NON/Lt background. Four new strains in this panel were exposed to chronic rosiglitazone treatment. Only one, NONcNZO8 (designated RCS8), exhibited an F1-like hepatosteatotic response. In both the F1 and RCS8 males, this adverse effect correlated with rosiglitazone suppression of already impaired hepatic phosphatidylcholine biosynthetic enzymes in both arms of the biosynthetic pathway, the phosphatidylethanolamine methyl- transferase pathway, and the CDP-choline pathway, including choline kinase and CTP-cholinephosphate cytidylyltransferase. This adverse response was not reproduced by CL316,243, a ß3-adrenergic receptor agonist with potent antihyperlipemic effects. Genome comparison showed that RCS8 differed from the other strains in carrying NZO-derived genome on virtually all of chromosome 16 and in smaller segments on chromosomes 6, 14, and 17. Thus, these RCSs present a panel of new mouse models exhibiting differential levels of obesity and diabetes as well as different drug responses. This panel can be used to screen for treatments for type 2 diabetes and its complications.
Address correspondence and reprint requests to Edward H. Leiter, PhD, The Jackson Laboratory, 600 Main St., Bar Harbor, ME 04609. E-mail: ehl{at}jax.org
Abbreviations:
PEMT, phosphatidylethanolamine methyltransferase; PPAR
The most common forms of human type 2 diabetes are polygenic in origin with contributions from both sides of the pedigree. In mice, obesity-associated diabetes (diabesity) also has a complex polygenic basis. Admixing genomes of unrelated strains of inbred mice provides insight as to how a complex disease such as type 2 diabetes can become more common as genetic heterogeneity increases in an outbreeding human population. As strains of laboratory mice become increasingly more inbred, genomes are selected that are capable of sustaining reproductive fitness in the face of increased mutational load. As mutations occur that affect the expression or function of multiple quantitative trait loci (QTLs), continued reproductive fitness requires co-adaptation within the overall genetic architecture for sets of QTLs that must interact to maintain metabolic homeostasis. Variation among inbred strains in physiologic parameters associated with glucose homeostasis (e.g., nonfasting plasma glucose or insulin concentrations) reflects strain-dependent systemic adaptations to strain-unique genetic polymorphisms at multiple QTLs. To model how a genetic outcross can destabilize such metabolic adaptations, we previously crossed two unrelated inbred strains, nonobese nondiabetic (NON) and New Zealand obese (NZO). The NON/Lt strain carries QTLs conferring latent type 2 diabetes susceptibility, whereas NZO/HlLt mice carry numerous QTLs contributing to male diabesity in a threshold fashion (1). Whereas diabesity spontaneously developed in 0% of NON/Lt males and in The hyperglycemia, hyperinsulinemia, and hyperlipidemia that developed in diabetic F1 males were effectively suppressed by chronic treatment with the thiazolidinedione (TZD) compound, rosiglitazone, incorporated into the diet. However, these antidiabetic effects were accompanied by marked exacerbation of an underlying hepatosteatosis (3). Because choline deficiency is known to produce hepatosteatosis in mice (4,5), we tested the lipid composition of milk from untreated F1 lactating dams; this analysis showed a marked deficiency in all classes of phosphatidylcholine (S. Watkins, Lipomics, Sacramento, CA, personal communication). This result suggested to us that the genetic milieu of the F1 mice combined QTLs that adversely affect the two phosphatidylcholine biosynthetic pathways in liver and that this effect is influenced by treatment with TZDs.
In mouse liver, the phosphatidylethanolamine methyltransferase (PEMT) pathway contributes
Our objective is to use mouse models of diabesity to understand how natural allelic variation at multiple genetic loci can affect responsiveness to drug-mediated therapy. The polygenic obesity/diabesity syndromes in mice are particularly relevant because the "metabolic syndrome" associated with most human diabetogenic obesities also has a complex genetic basis. To elucidate the pharmacogenetic basis for the adverse responses of the F1 liver to TZD exposure, we compared the hepatic activities of the phosphatidylcholine biosynthetic enzymes in both parental strains with the F1 and in a panel of genetically characterized recombinant congenic strains (RCSs). The RCSs were produced by backcrossing the F1 for two cycles onto the parental NON/Lt background, with selection for different subsets of diabesity QTLs sorting into each RCS (12). For example, RCS10 contains the greatest number of diabesity contributions from both parental backgrounds (on chromosomes 1, 4, 5, 11, 12, and 18) and develops the highest (F1-like) diabetes frequency (90–100%) of any of the RCSs. In contrast, RCS1, selected to contain the diabetogenic QTLs mainly on chromosomes 1 and 15, showed only 50% diabetes frequency. RCS2, a strain lacking most of the NZO-derived diabesity contribution on chromosome 15 while otherwise carrying the same genetic composition as RCS1, showed an even lower (25%) diabesity frequency (12). In the present study, we show that hybridizing the two parental genomes produced an unusual lowering of hepatic phosphatidylcholine biosynthetic enzyme activities in control F1 males and that rosiglitazone treatment further reduced these enzymatic activities. Screening a selected panel of RCSs showed that only one, NONcNZO8, developing a diabesity frequency of
(NON x NZO)F1, NON/Lt, and four previously reported RCSs (12) (NONcNZO1 [designated RCS1], NONcNZO2 [RCS2], NONcNZO8 [RCS8], and NONcNZO10 [RCS10]) were bred in our research vivarium. They were housed two to four per pen in double-pen Plexiglas boxes on shaved pine bedding and given free access to food and acidified water. All mice shared the same mouse room with controlled temperature and humidity and a 14-/10-h light/dark cycle. Because diabesity is male sex limited in the models studied (1,12), only males were used. Eight to 12 males of each strain at 8 weeks of age were split into two groups. A control group of males continued to receive the maintenance (control) diet (NIH-31 containing 6% fat; Purina Test Diets, Richmond, IN). The remaining one-half received the same diet supplemented with 50 mg rosiglitazone/kg (a gift from Dr. S. Smith, GlaxoSmithKline). Mice were maintained on the irradiation-sterilized diets for 14 weeks. At experiment termination (22 weeks of age), the mice were killed by CO2 inhalation, and blood was collected by heart puncture using a heparinized 25-gauge needle and syringe. Livers were collected, weighed, and homogenized immediately for enzyme assay. A piece of each liver was isolated for histologic analysis. For comparative purposes, a separate study group of five F1 males at 8 weeks of age were placed for 6 weeks on the NIH-31 diet alone or supplemented with 0.001% CL316,243, a ß3-adrenergic receptor agonist (gift of Dr. Kurt Steiner, Wyeth Research, Princeton, NJ) previously shown to exert potent antihyperglycemic/antihyperlipemic effects in the NZO/HlLt strain (13). Also, to compare compounds within the TZD class, troglitazone (a gift of Dr. H. Horikoshi, Sankyo, Tokyo) was prepared in pelleted NIH-31 diet at a concentration of 2 g/kg diet based on a previous report (14). Groups of five NON/Lt, (NZO x NON)F1, and RCS8 males were placed on this diet from 6 to 16 weeks of age.
Clinical markers.
Histology studies.
Cell fractionation and enzyme assays.
Western blot analysis.
Unusual suppression of both arms of the phosphatidylcholine biosynthetic pathway in F1 livers. Although F1 hybrid mice generated by outcross of two inbred progenitor strains are generally considered more "robust," the combination of NON/Lt and NZO/HlLt genomes actually produced a "negative heterosis" in two of three critical enzymes in hepatic phosphatidylcholine biosynthesis. As shown in Fig. 1A and B, both PEMT-specific and choline kinase-specific activities were significantly lower in total homogenates of F1 liver than in total homogenate from either the NON/Lt or NZO/HlLt parental strains. A comparable depression of either of these enzymatic activities could not be modeled by a 1:1 mixing of liver total homogenates from NON and NZO directly. Instead, an intermediate activity for both PEMT and choline kinase was observed (Fig. 1A and B), as would be expected for additive activity contributions from each parental source. As shown in Fig. 1C, this unexpected hybrid inhibitory effect on PEMT and choline kinase activities did not extend to cholinephosphate cytidylyltransferase (the rate-limiting enzyme in the CDP-choline pathway). There were no significant activity differences distinguishing either parental strain from the F1. Despite the significantly reduced PEMT and choline kinase enzymatic activities distinguishing F1 from the two parental strains, Western blots using antibodies against PEMT or against each of the two choline kinase subunits, choline kinase- and choline kinase-ß, showed no differences in enzyme protein concentration among the three strains (Fig. 1D). Thus, the F1 metabolic milieu was selectively producing post-transcriptional changes potentially affecting phosphatidylcholine production via the PEMT pathway and downregulating activity of one enzyme (choline kinase) in the CDP-choline pathway.
Rosiglitazone treatment suppressed all three key phosphatidylcholine-biosynthetic enzymes in F1 livers. The F1 males are distinguished from both parental strain males by a higher frequency of diabesity development (1). Although diabetic hyperglycemia and hyperlipemia were reversed by supplementing the diet with rosiglitazone, this therapy was accompanied by marked hepatosteatotic side effects (3). This effect was not limited to rosiglitazone, because troglitazone fed over the same time frame at 2 g/kg diet produced the same degree of hepatosteatosis (data not shown). The finding that two out of three phosphatidylcholine biosynthetic enzyme activities assayed in F1 liver were unexpectedly reduced below both parental strains, coupled with our previous observation of increased triglyceride accumulation in livers of rosiglitazone-treated F1 males (3), suggested that increased F1 hepatosteatotic sensitivity after TZD treatment entailed a further drug-mediated impairment of these enzymatic functions. Data in Fig. 2 strongly support this. All three enzyme activities measured in untreated ("control") F1 liver homogenates were significantly inhibited by chronic exposure to rosiglitazone at a dosage of 50 mg/kg diet (Fig. 2A). This included not only PEMT and choline kinase activities already significantly reduced in untreated "control" F1 livers (Fig. 1A and B), but also cholinephosphate cytidylyltransferase activity that, in samples of untreated control F1 liver, did not differ significantly from parental activities. Western blot showed that rosiglitazone markedly suppressed concentrations of both choline kinase subunits (choline kinase- and choline kinase-ß) in the F1 liver, whereas PEMT concentration were quite comparable between control and rosiglitazone group (Fig. 2C, left panel). In contrast to rosiglitazone, CL316,243 significantly stimulated the expression of choline kinase- subunit and slightly increased PEMT concentration (Fig. 2C, right panel). The regulation of protein expression correlates well with the actual enzyme activities. (Fig. 2A and B). Even though the rosiglitazone concentration used in the present study was fourfold lower than reported previously (3), histologic assessment of samples of the F1 livers used in the current experiments revealed macrovesicular hepatosteatosis (Fig. 3) as severe as reported previously using the higher drug dosage. A null mutation in the PEMT gene in mice produces increased sensitivity to development of hepatic steatosis (4,5). None of the phosphatidylcholine biosynthetic enzymes measured were completely ablated in rosiglitazone-treated F1 mice. Yet the ability of rosiglitazone to impair already inhibited PEMT and choline kinase activities, as well as to significantly reduce activity of cholinephosphate cytidylyltransferase, reveals an unusual genotype-drug interaction that compromises the livers phosphatidylcholine biosynthetic pathways. This impairment was limited to TZD, because treatment with CL316,243, an even more effective antihyperlipemic and antihyperglycemic agent in the F1 model (13), did not impair activity of any of the three enzymes. Although mean activities of all three were elevated above male mice on control diet, the differences did not achieve statistical significance after the 6-week treatment period (Fig. 2B).
Phenotypic screen of RCS sensitivity to rosiglitazone-enhanced steatosis. RCSs represent unique tools for dissecting drug-genome interactions. Introgression of different combinations of selected NZO-derived QTLs onto the diabesity- and steatosis-resistant NON parental background produced a panel of different gene combinations for testing the genetic basis of the steatosis in the rosiglitazone-treated mice. After 14 weeks of rosiglitazone treatment, NON males on control diet maintained normoglycemia. Rosiglitazone treatment did not produce an increase in body weight, although plasma insulin was significantly reduced (Table 1). However, with the exception of RCS1, most RCS strains, similar to the F1 males, responded to rosiglitazone-stimulated body weight gain (Table 1). As previously reported, RCS2 showed resistance to diabesity development, whereas RCS1, RCS8, and RCS10 all exhibited diabesity syndromes of differential severity as manifested by increases in both plasma glucose and insulin concentrations (Table 1). Also, in all three, the insulin-sensitizing action of rosiglitazone was reflected by significant decreases in plasma insulin. Nevertheless, differences in antihyperglycemic effects were noted. Drug-mediated plasma glucose decreases were significant in both RCS1 and RCS10 mice, but the decline in treated RCS8 mice did not achieve significance (Table 1). Moreover, plasma insulin was much more variable in control RCS8 males. Although rosiglitazone treatment reduced mean plasma insulin concentrations and reduced variability in the measurement, the difference again failed to achieve statistical significance. This suggested that the diabesity syndrome developing in RCS8 was different from that in RCS1 and RCS10.
Differential sensitivity of RCS to rosiglitazone-mediated hepatic steatosis. Results of morphometric profiling of hepatic fat infiltration into livers at the end of rosiglitazone treatment are shown in Fig. 3. Comparisons are with the NON parental males that exhibited minimal lipidosis that was not exacerbated by the 14-week period of rosiglitazone treatment. This absence of a drug effect in NON is in contrast to the extreme sensitivity of the F1 to this adverse side effect. Livers of F1 males on control diet exhibited an elevated mean basal concentration of fat (mostly microvesicular). Basal microvesicular fat in RCS males fed control diet were all comparable with the NON basal level (Fig. 3), despite marked differences in the severity of the diabesity syndromes in the RCS mice (Table 1). Neither the diabesity-sensitive RCS1 nor the diabesity-resistant RCS2 differed from NON in terms of resistance to rosiglitazone-exacerbated lipidosis. The most diabesity-prone strains, RCS8 and RCS10, showing comparable basal microvesicular lipidosis on control diet (14.4 and 18%), both exhibited significant increases in macrovesicular fat accumulation in response to rosiglitazone treatment. However, only RCS8, which did not respond as well as RCS10 to the glucose-lowering action of rosiglitazone, showed the extreme hepatosteatotic response indistinguishable from that observed in drug-hypersensitive F1 males (Fig. 3A). The rosiglitazone-mediated increase in mean fatty area in RCS8 liver (45% in rosiglitazone treated versus 17.3% in control) compared well with that in the F1 liver (51% in rosiglitazone treated versus 26% in control, Fig. 3B). In RCS10, basal lipidosis was comparable with RCS8 (14.4%), whereas the rosiglitazone-mediated increase in mean fatty area (to 27.8%) was 37% lower than the rosiglitazone-mediated increase observed in RCS8 (Fig. 3B). This result correlated with the finding that treated RCS8 males exhibited the maximal elevation in plasma alanine:leucine aminotransferase, a clinical marker diagnostic of hepatotoxic stress (data not shown). Hence, among the panel of RCSs tested, RCS8 has the most potent combination of QTLs, fixed in the homozygous state, that predispose to the adverse drug response originally observed in F1.
Only RCS8 responds to rosiglitazone-mediated inhibition of phosphatidylcholine biosynthetic enzymes.
Genomic differences in RCSs. The genome-wide scans distinguishing the panel of RCSs used in this study are available online (http://www.jax.org/staff/leiter/labsite/type2_genomics.html), and the known diabesity QTL markers have been published previously (12). The PEMT-encoding gene maps to one of these diabesity QTL regions on chromosome 11 (originally designated as Nidd3). Two linkage peaks on this chromosome were detected; this second peak contained another enzyme associated withphosphatidylcholine metabolism, phosphatidylcholinetransfer protein, Pctp, that we have recently found to be mutated in NZO (H.-J.P., E.H.L., unpublished data). The choline kinase-ß subunit encoding gene maps to another diabesity QTLs on chromosome 15 ("Nidd" designations were discontinued because of the large numbers of the complex epistatic interactions discovered among a multiplicity of QTLs identified [2]). Because the RCS8 response to rosiglitazone was unique among the panel tested (RCS1, RCS2, RCS8, and RCS10), comparative analysis of the genome-wide differences distinguishing RCS8 from the other three RCSs could potentially delineate the genetic basis for the increased sensitivity of this strain. The RCS8-unique set of NZO-derived genomic segments included virtually all of chromosome 16 (markers spanning 8.6–55 cM) and considerably smaller segments marked by single polymorphic simple sequence repeats on chromosome 6 (marked only by D6Mit275, 25.5 cM), chromosome 14 (markers D14Mit195-D14Mit170; 44.3–63 cM), and chromosome 17 (markers D17Mit16-D17Mit11; 17.4–22.8 cM). Hence, the adverse drug response in this strain may likely be controlled by a gene or genes on only one of these chromosomes, most likely chromosome 16. Reciprocal outcrosses between RCS2 (hepatosteatosis resistant) and RCS8 (hepatosteatosis susceptible) show that the F1 males inherit the high hepatosteatotic responsiveness to rosiglitazone in a dominant fashion (data not shown). This would be consistent with the high hepatosteatosis sensitivity of F1 males generated by the reciprocal NON/Lt x NZO/HlLt outcrosses.
Highly inbred strains of mice, such as NON/Lt and NZO/HlLt, cannot accurately model for patient-specific drug responses in an outbred human population, whereas a genetically heterozygous F1 hybrid would be more representative in this regard. The present study used RCSs to differentially distribute genetic susceptibility to TZD-exacerbated hepatosteatosis into a genetic background (NON/Lt) that is resistant to these effects. It is conceivable that the heightened sensitivity of the F1 and RCS8 males to TZD-exacerbated hepatosteatosis represents a component of insulin resistance uniquely present in these two stocks. We compared the antihyperglycemic and insulin-sensitizing action of troglitazone at 2 g/kg diet with rosiglitazone at 50 mg/kg diet in both the F1 and the deriviative RCS8 stock. Troglitazone was much less effective either in suppressing hyperglycemia or in promoting pancreatic ß-cell regranulation with insulin (data not shown). Nevertheless, both TZDs produced hepatosteatosis of equal severity. Detailed analysis of insulin resistance in TZD-treated RCS8 versus RCS10 males has not yet been done. Euglycemic-hyperinsulinemic clamp studies have shown that untreated RCS10 males are extremely insulin resistant (18), yet males of this strain are significantly more resistant than RCS8 to hepatosteatosis. Thus, we have no evidence at this point to conclude that intractable insulin resistance alone distinguishes steatosis-prone from steatosis-resistant RCSs.
Troglitazone, the first TZD marketed for treatment of type 2 diabetes, was withdrawn because of extreme hepatotoxicity in a very small percentage of patients taking the compound (19). Drugs of the later generation TZDs, of which rosiglitazone is an example, have been shown to reduce hepatic triglycerides in most patients rather than promote triglyceride accumulation and steatosis (20). Thus, availability of the panel of genetically characterized RCS models of diabesity with differential TZD sensitivities is useful for understanding the genetic basis for susceptibility to adverse drug responses. Such pharmacogenetic knowledge might distinguish that small percentage of patients who should not take this class of compounds. RCS8 and RCS10 exhibit comparable levels of obesity and basal hepatic lipidosis with a control diet, but only RCS8 experiences a drug-exacerbated steatosis as extreme as originally described in the F1. Interestingly, RCS10 rather than RCS8 is more comparable with the F1 in terms of the numbers of shared diabesity QTLs and, consequently, diabetes frequency (90–100%) and disease severity defined by nonfasting plasma glucose levels. RCS8 develops a lower frequency of diabesity with hyperglycemia not establishing consistently before 20 weeks of age (12). Thus, it was surprising that the moderate shifts in the glycemic status of RCS8 males before 20 weeks were not as responsive to rosiglitazone treatment as were the more hyperglycemic RCS10 males. Because the NON/Lt strain is completely resistant to rosiglitazone-exacerbated steatosis, the steatosis-promoting gene or genes must be NZO in origin. We surmise that this NZO-derived locus or loci promoting the phenotype of drug-exacerbated steatosis may be distinct from the set of NZO-contributed diabesity QTLs that synergize with certain NON-derived QTLs to increase frequency of diabetic hyperglycemia development in the F1 male. This inference is reinforced by the fact that the three steatosis-resistant RCSs collectively carry the spectrum of known diabesity QTLs. RCS8 has been fixed for known diabesity QTLs on chromosomes 1 and 11 (from NZO) and chromosomes 4 and 18 (from NON). The most notable genetic difference distinguishing RCS8 from the others is its NZO origin of chromosome 16 in its entirety. The cholinephosphate cytidylyltransferase-
Although the PEMT pathway only accounts for 30% of phosphatidylcholine biosynthesis in liver, it is required for secretion of lipoproteins (8), so that PEMT inhibition significantly impairs the bulk incorporation of triglycerides into lipoprotein particles for export (26). In livers of mice with a genetically disrupted Pemt allele fed a standard diet, a diet deficient in choline, or a choline-enriched diet, there was hepatic steatosis and significantly higher frequency of apoptotic cells compared with wild-type controls (5). Similarly, reductions in multiple classes of plasma lipoproteins were also observed in mice with a disrupted cholinephosphate cytidylyltransferase-
It remains to be established whether this unusual pharmacogenetic effect is mediated via direct effects of peroxisome proliferator–activated receptor- The insulin-resistant NZO parental strain donating obesity/diabesity QTLs to the RCS panel shows a dissociation of insulins effects on genes involved in glucose and lipid metabolism. Insulin-suppressible genes associated with hepatic gluconeogenesis and glycogenolysis are not suppressed, but insulin-inducible genes associated with lipogenesis are markedly upregulated (29). A combined gene transcription and lipid metabolome analysis of rosiglitazone-treated F1 livers indicated that lipid uptake and rosiglitazone-stimulated de novo lipid biosynthesis were not adequately compensated by lipid export (3). In the case of the RCSs tested in this report, the physiologic profiling of plasma glucose and insulin responses to rosiglitazone in Table 1 suggest different degrees of insulin resistance. Unlike RCS10, whose insulin resistance was clearly diminished by rosiglitazone treatment as reflected by a restoration of normoglycemia, RCS8 appeared less responsive despite a considerably milder mean hyperglycemia over the time period studied.
In summary, the RCS analysis reported herein has allowed biochemical dissection of the marked sensitivity of the F1 diabesity model to TZD-exacerbated steatosis. Our results indicate a constitutive impairment of hepatic phosphatidylcholine biosynthetic enzyme functions that is further exacerbated by TZD treatment. The fact that, unlike RCS8, the NON/Lt parental background strain, RCS1, RCS2, and RCS10 do not exhibit these extreme responses to TZD establishes this RCS panel as potentially useful pharmacogenetic screening tools. Clearly, there are well-known differences between humans and mice in terms of rosiglitazone effects on liver fat accumulation (20). Rosiglitazone treatment of humans is associated with decreases, not increases, in liver fat (20). In mice, inability to express PPAR
E.H.L. has received support from National Institutes of Health Grant DK-56853. H.-J.P. was supported by a mentor-based fellowship from the American Diabetes Associ-ation. D.E.V. has received a grant from the Canadian Institutes of Health Research. TJL Institutional shared services were supported by National Cancer Institute Center Support Grant CA-34196. D.E.V. holds the Canada Research Chair in Molecular and Cell Biology of Lipids and is a Medical Scientist of the Alberta Heritage Foundation for Medical Research. We thank Pam Stanley and Sandra Ungarian for technical assistance. We thank Drs. Rick Woychick and Jürgen Naggert at The Jackson Laboratory for critical reading of the manuscript.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication December 17, 2004 and accepted in revised form February 18, 2005
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||