Applications of metabolomics for understanding the action of peroxisome proliferator-activated receptors (PPARs) in diabetes, obesity and cancer

The peroxisome proliferator-activated receptors (PPARs) are a set of three nuclear hormone receptors that together play a key role in regulating metabolism, particularly the switch between the fed and fasted state and the metabolic pathways involving fatty-acid oxidation and lipid metabolism. In addition, they have a number of important developmental and regulatory roles outside metabolism. The PPARs are also potent targets for treating type II diabetes, dyslipidemia and obesity, although a number of individual agonists have also been linked to unwanted side effects, and there is a complex relationship between the PPARs and the development of cancer. This review examines the part that metabolomics, including lipidomics, has played in elucidating the roles PPARs have in regulating systemic metabolism, as well as their role in aspects of drug-induced cancer and xenobiotic metabolism. These studies have defined the role PPARδ plays in regulating fatty-acid oxidation in adipose tissue and the interaction between aging and PPARα in the liver. The potential translational benefits of these approaches include widening the role of PPAR agonists and improved monitoring of drug efficacy.

variety of other cancers (Table 1). This is an example of the complications of targeting transcription factors that have such a wide variety of roles in the body. More recently, thiazolidinedione drugs, compounds derived from a ring structure containing a sulfur and two ketone groups that target PPARγ, have been linked to cardiovascular disease, and specifically myocardial ischemia [4], prompting the United States (US) Food and Drug Administration (FDA) to issue a safety alert on one of the most widely prescribed classes of anti-diabetes drugs. However, since the publication of this alert, the American Heart Association and American College of Cardiology Foundation have published an advisory note stating that health professionals should weigh up the benefits from improved insulin sensitivity and glycemic control afforded by thiazolidinedione drugs such as rosiglitazone and pioglitazone against the inconclusive role in ischemic heart disease [5].
Metabolomics, also referred to by some as metabonomics, is the comprehensive analysis of the metabolic complement of a cell, tissue, biofluid or organism. The approach makes use of techniques of analytical chemistry to profile as wide a range of metabolites as possible, in a quantitative or semi-quantitative manner, to follow a variety of physiological and pathophysiological stimuli such as disease progression, gene function or a toxic insult. Metabolomics has been used as part of personalized medicine to identify populations that will respond to a particular drug intervention [6], identify predictive biomarkers of disease presence and future prognosis [7][8][9][10], conduct functional genomics [11][12][13] and identify genes involved in diseases such as diabetes [14,15].
In this review we assess the use of metabolomics to elucidate the part PPARs play in regulating systemic meta bolism. Although physiology-based approaches have been highly successful at defining the broad actions of these receptors, much of the detail in terms of how regulation is exerted has still to be defined. Metabolomics is viewed in its broadest sense to also include lipidomic approaches, which focus on the analysis of intact lipids, as well as fluxomic approaches, which assess a range of pathway activities using stable isotope techniques. We begin by briefly reviewing the important role animal models have had in the understanding of PPAR activities, before describing the advantages of a metabolomics approach and the insights gained. The focus of much of the research into PPARs has been on their role in treating type II diabetes, obesity and dyslipidemia, but it has also been known since the original discovery of these receptors that agonists for all three receptors alter the relative risks for a variety of cancers. In addition we are now beginning to understand how the PPAR signaling systems interact with other molecular signaling pathways to modulate cellular metabolism, proliferation and inflammation, indicating their profound and diverse roles in regulation within the cell (Figure 1).

Understanding PPARs through animal models: what classic physiology tells us about the receptors
Of the three PPARs, the first discovered was PPARα. It was found in rodents following the observation that it was the target of drugs that caused peroxisomes (organelles that primarily break down lipids) to proliferate in the liver [16]. Remarkably, despite the important regu latory roles that PPARα plays in regulating metabolism in key organs, PPARα-knockout mice have been generated [17]. In addition to investigating the role PPARα plays in peroxisome proliferation, the mouse also became a useful research tool for following its role in regulating metabolism. Kersten and co-workers [18] either fed the null mice a high-fat diet or fasted them for 24 hours. Both challenges induced an accumulation of lipid in the liver, but fasting demonstrated that the PPARα-null mouse had reduced capacity for gluconeogenesis and ketogenesis, experienced hypothermia and had increased levels of circulating free fatty acids. Similarly, Leone and colleagues [19] used the knockout mouse to investigate how fat metabolism is regulated by PPARα during the fasting response. They found that fasting induced steatosis (retention of fatty acids) in the liver and heart, with a failure to upregulate the ketogenic and glucogenic pathways in the liver. These and numerous other studies indicate that PPARα regulates a number of catabolic pathways; by contrast, the activity of PPARγ appears to be mainly one of lipid storage and synthesis, as well as playing an important role in adipogenesis [6,7]. PPARγ is also the target of a number of potent insulin-sensitizing agents, including the PPARγ agonists the thiazolidinediones. The expression of PPARγ is also regulated by diet, with marked reductions in its expression detected in adipose tissue during fasting, but induction by high-fat feeding, demonstrating its role in fat storage [20]. There are three major isoforms of PPARγ, with PPARγ1 found in nearly all tissue apart from muscle, PPARγ2 found mostly in adipose tissue [21] and PPARγ3 found in macrophages, the large intestine and adipose tissue. Failure to express PPARγ causes embryonic lethality in mice, but heterozygous mice and humans with a Pro12Ala polymorphism in PPARγ2 appear protected from the development of type II diabetes, suggesting that PPARγ is part of the thrifty genotype thought to predispose the individual to insulin resistance [22]. Although PPARγ is required for life, PPARγ2-null mice are viable and have been used to investigate its role in insulin sensitivity, Figure 1. A schematic diagram illustrating the cross-talk of the three PPAR receptors and the metabolic pathways they interact with. For each PPAR, the initial step involves the receptor binding a ligand to activate it. Then, the retinoid X receptor (RXR), a nuclear receptor activated by 9-cis retinoic acid, heterodimerizes with a PPAR prior to the subsequent binding of the complex with DNA. Key: NSAID, non-steroidal anti-inflammatory drug; TZD, thiazolidinedione, a class of drugs that bind to PPARγ and have insulin-sensitizing properties; VLDL, very low density lipoprotein (used to transport lipids in the blood). where it appears to be an important factor in insulin sensitivity in adipose tissue [23].

Fats in VLDL chylomicrons
PPARδ has been much less studied than the other PPARs. Despite the important role the receptor is thought to play in controlling systemic metabolism, selective agonists have only recently been developed for PPARδ. This receptor has a role in adipocyte differ entiation, lipoprotein metabolism, atherosclerosis and βoxidation in adipocytes, skeletal muscle and the heart [24]. It also, in part, regulates fat metabolism in skeletal muscle, and over-expression of PPARδ in the mouse almost doubles this animal's endurance exercise capacity [25]. PPARδ also regulates glucose homeostasis, with PPARδ-null mice having reduced capabilities for glucose handling; PPARδ agonists also reduce insulin resistance and weight in obese obob mice [26], which, as they lack the leptin receptor, have such a high appetite that they develop obesity and type II diabetes as a result of over-nutrition.
Of course the role of PPARs in regulating metabolism in humans has been the subject of great interest, particularly within the pharmaceutical industry. It would be beyond the scope of this review to do full justice to the pharmacology of PPAR agonists in humans, and instead it seems more profitable to consider the rare human mutations of these receptors. Dominant negative mutations have been reported in PPARγ in humans [27]. Savage and colleagues [27] reported a dominant negative mutation of PPARγ arising from a proline-467-leucine mutation, which resulted in severe peripheral and hepatic insulin resistance, partial lipodystrophy (the loss/ reduction of certain fat depots) and hepatic steatosis. Although relatively rare, individuals with this mutation have proved to be invaluable in investigating the role of PPARγ in humans, and a useful contrast with their rodent counterparts.

Advantages of the metabolomic approach to understanding the roles of the PPARs
Before examining in detail the information gained from metabolomic studies, it seems prudent to make a few points concerning the relevance of the approaches to studying functional genomics, and in particular to investigating systems that have a strong metabolic component. Firstly, any functional genomic study focusing on a system with a strong metabolic component is well suited to metabolomics, as many of the key outputs and regulatory stages will involve a significant number of metabolites. Secondly, the technologies used in metabolomics have largely been developed as high-throughput analytical chemistry tools, capable of analyzing relatively large numbers (50 to 200) of samples per day, with relatively low consumable costs, relying on relatively cheap solvents and derivatization agents for mass spectrometry and off-the-peg chromatography. Being cheap on a persample basis allows the approach to be used to analyze the range of tissues and biofluids necessary to study the function of transcription factors across the whole organism as part of a systems biology approach.
This system-wide approach is important for the PPARs because these nuclear hormone receptors are expressed across a range of tissues, and each plays a key role in regulating systemic metabolism in underfed or fasted conditions. In addition, metabolomics allows the analysis of large numbers of replicates of the same sample type. Although most animal models of type II diabetes and obesity have centered on rare monogenic causes, the disease is generally considered to be polygenic in the vast majority of humans; hence, to understand the role of the PPARs in predisposing individuals to metabolic syndromes, or in the treatment of this disease, it is necessary to either investigate a range of monogenic models or examine crosses between strains that are relatively predisposed to and relatively protected from developing a particular disorder. Both approaches necessitate a large number of samples, which makes metabolomics, an approach that is cheap on a per-sample basis, ideal for performing 'functional genomic triage' . As part of a polyomic approach, it may be too expensive to perform transcriptomics on multiple tissues or too time consuming to perform proteomics on such a large scale. However, high-throughput metabolomics can be used to identify the tissues that are most interesting to follow up as part of a triage process.

Metabolomic studies of PPARα: targeting gluconeogenesis to treat diabetes and examining the interaction with aging
The nuclear receptor PPARα plays a central role in regulating the oxidation of fatty acids in the liver and skeletal muscle and gluconeogenesis in the liver, particularly during the fasted state or following high-fat feeding. As already mentioned, one of the surprises of the PPARα-null mouse was that this mouse was viable, and indeed phenotypes were only apparent under relatively large physiological stresses. The role of PPARα in terms of regulating gluconeogenesis was elegantly demonstrated by Xu and co-workers [28], using gas chromatographic mass spectrometry (GC-MS) analysis of labeled isotopic isomers (called isotopomers) following subcutaneous infusion, via Alza miniosmotic pumps, of either [2-13 C]glycerol, [U-13 C 3 ]lactate or [U-13 C 6 ]glucose. Following a 17-hour fast, while hepatic glucose production was found to be increased in the PPARα-null mouse, the stable isotope approach demonstrated reduced capacity for converting lactate into glucose, but increased capacity for the conversion of glycerol into glucose via gluco neogenesis. This identified the perturbation in gluconeogenesis as being at the start of the pathway, with the expression of pyruvate kinase found to be reduced 16-fold in the fasted PPARα-null mouse compared with the fasted wildtype control. In a follow-up study, Xu and colleagues [29] used the same stable isotope approach to demonstrate that a failure to express PPARα reduces futile cycling (the cyclic conversion between two substrates, which consumes adenosine triphosphate (ATP) during the process) in the liver but, to compensate for this, peripheral uptake of glucose is increased.
To examine the role of PPARα across the whole animal, Atherton and co-workers [30] used a combination of high-resolution 1 H nuclear magnetic resonance (NMR) spectroscopy, GC-MS and liquid chromatograph-mass spectrometry (LC-MS) to profile liver, skeletal muscle, smooth muscle, cardiac tissue and white adipose tissue in one-month-old PPARα-null mice and compared them to wild-type controls. Despite conducting the study in mice fed on standard chow, the metabolomic approach readily distinguished the null mouse from the wild-type control in all tissues, with the most marked effects detected in the liver associated with alterations in glycolysis, the citric acid cycle, amino acid metabolism and fatty acid oxidation. One of the biggest surprises was that white adipose tissue demonstrated metabolic alterations, despite very low expression of PPARα in control adipose tissue. In support of this, Knauf and colleagues [31] report that glucose clearance by adipose tissue is increased in the PPARα-null mouse, even when expression is restored in the liver, and may even be mediated by the brain.
Expanding their metabolomic approach to investigate the interaction of PPARα expression with age, Atherton and co-workers [32] modeled the metabolic changes between 3 and 13 months of age in the PPARα-null mouse, using multivariate statistics ( Figure 2). Hepatic glycogen and glucose decreased with age for both the mutant and wild-type mice, but the tissue from the null mouse had lower concentrations, with evidence of ageonset hepatic steatosis on a normal chow diet. This approach demonstrated that any phenotype should be placed within the context of the normal aging process, to understand the role of the gene across the lifetime of the organism.
An important area in drug development is the identification of potential markers of drug efficacy, to gauge how the individual responds to a pharmacological stimulus. To address this, Zhen and co-workers [33] used LC-MS to profile urinary metabolites in control and PPARα-null mice and examined the action of Wy-14,643 ([4-chloro-6-(2,3-xylidino)-2-pyrimidinylthio] acetic acid), a potent PPARα agonist. Using principal compo nents analysis (PCA) to interrogate the data, the PPARα-null groups of mice co-clustered (occupied similar positions on the PCA plot, indicating similar metabolic changes across the group compared with the control group), regardless of drug treatment, while the control groups formed separate clusters. This agonist increased urinary concentrations of 11beta-hydroxy-3,20-dioxopregn-4-en-21-oic acid, 11beta,20-dihydroxy-3-oxopregn-4-en-21oic acid, nicotinamide, nicotinamide 1-oxide, 1-methylnico tinamide, hippuric acid, and 2,8-dihydroxyquinolinebeta-d-glucuronide, and decreased concentrations of xanthurenic acid, hexanoylglycine, phenylpropionyl glycine and cinnamoylglycine. This highlighted the action of PPARα in regulating tryptophan, corticosterone and fatty-acid metabolism and on glucuronidation. This approach was extended by Patterson and colleagues [34], who examined the use of fenofibrate (a fibrate drug that targets PPARα; Table 1) treatment in humans, identifying pantothenic acid and short-chain acyl carnitines as markers of increased β-oxidation, and validating their results using the PPARα-null mouse. Following up the production of 11beta-hydroxy-3,20-dioxopregn-4-en-21oic acid and altered steroid production by PPARα, Wang and co-workers [35] demonstrated cross-talk between the receptor, neuropeptide Y and agouti-related protein interacting with the hypothalamic-pituitary-adrenal axis.
In addition to its role in regulating the fasted response, PPARα also plays an important role in toxicology, and the original drive for the production of the PPARα-null mouse was to investigate the role PPARα plays in druginduced hepatocarcinogenicity. Ohta and colleagues [36] investigated the early stages of fenofibrate toxicity in the Fischer rat. High doses of fenofibrate are known to be carcinogenic in the liver, via a non-genotoxic mechanism. In this study, Ohta and colleagues used a combination of LC-MS and GC-MS to highlight that part of the toxicity was caused by the generation of reactive oxygen species in the liver. Manna and co-workers [37] used metabolomics and the PPARα-null mouse to investigate alcoholinduced hepatotoxicity. The PPARα-null mouse is particularly predisposed to alcohol-induced hepatotoxicity, as a failure to express PPARα reduces the conversion of tryptophan to NAD + (nicotine adenine dinucleotide + ), which exacerbates the action of alcohol and its reducing potential on the liver. The combination of these events produced a characteristic increase in indole-3-lactic acid, which could be followed in the urine of these animals [37]. This suggests that PPARα agonists may be a viable treatment for alcohol-induced hepatic steatosis, and urinary indole-3-lactic acid would provide a marker of drug efficacy.

Metabolomic studies of PPARγ
Given its central role in regulating lipid storage and its importance as a potent target for restoring insulin sensitivity, PPARγ has been the subject of a number of influential lipidomic studies. One of the first studies, and an elegant demonstration that one does not necessarily need expensive equipment to perform lipidomics, was a study by Watkins and co-workers [38] to investigate the action of rosiglitazone in obese mice. Chronic treatment with rosiglitazone restored insulin sensitivity but increased hepatic steatosis. To assess these lipid changes, a combination of preparative thin-layer chromatography and gas-chromatograph flame-ionization detector (GC-FID) analysis was used to separate different lipid classes and analyze their fatty acid composition in plasma and key tissues. This provided enough detail to follow the increased synthesis of fatty acids within the liver and the remodeling of cardiolipin in the heart. Many of these changes were reflected in the plasma, suggesting that drug efficacy can be monitored effectively using metabolomics.
LC-MS-based lipidomics has been instrumental in determining the role of PPARγ in regulating adipose tissue morphology and expandability. There are three isoforms of PPARγ, with PPARγ2 being the most potent isoform for adipogenesis in vitro. Medina-Gomez and coworkers [39] used LC-MS, magnetic resonance imaging (MRI), transcriptomics and comprehensive physiology to investigate the role of PPARγ2 in regulating lipid storage in adipose tissue. Intriguingly, while ablation of PPARγ2 induced increased insulin resistance in chow-fed mice, this effect was not apparent in animals fed a high-fat diet. The increased basal insulin resistance was not associated with lipodystrophy, although adipocytes were larger in the transgenic mice following high-fat feeding. However, there is a high degree of compensation following the  ablation of PPARγ2. To address this, Medina-Gomez and colleagues [39] crossed the PPARγ2-null mouse with the obese Lep(ob)/Lep(ob) mouse to generate a PPARγ2 (-/-) Lep(ob)/Lep(ob) mouse, referred to as the POKO mouse. This animal has severe insulin resistance, β-cell failure, reduced fat mass and dyslipidemia. This demonstrates the importance of the PPARγ2 isoforms in regulating normal adipose tissue expandability. The LC-MS analysis of tissues and blood plasma also demonstrated that the isoform plays an important role in reducing potentially lipotoxic intermediates. Medina-Gomez and colleagues [40] followed up this initial study with a further characterization of the β-cell failure in the pancreas of the POKO mouse. They followed changes in the pancreas between 4 and 16 weeks of age, highlighting that the lipidome changes of the pancreas were relatively mild in comparison to the large changes detected in the plasma, liver, muscle and adipose tissue of the mice. Lipotoxic species were only detected at the later time point, suggesting the POKO mouse may be a useful tool to follow β-cell dysfunction secondary to peripheral insulin resistance or lipid infiltration.

Metabolomic studies of PPARδ/PPARβ
Metabolomic and lipidomic studies of PPARδ have been relatively rare, as much less work has focused on PPARδ, in part because of difficulties in developing selective ligands that do not target the other PPARs as well, and also because of the importance of PPARδ in normal development, complicating its genetic manipulation in rodent models. However, Roberts and co-workers [41] have used a comprehensive metabolomic and lipidomic approach to study the role of PPARδ in regulating systemic metabolism, by focusing on the obese and diabetic obob mouse and a highly selective PPARδ agonist. Given the high expression of PPARδ in both liver and skeletal muscle, both tissues were examined, alongside plasma, to understand how tissue changes influence the composition of plasma. This was contrasted with the metabolic changes induced in these tissues by a highly selective PPARγ ligand. While both ligands restored insulin sensitivity, only the agonist for PPARδ reduced the body weight of the obese obob mice. Both ligands induced an increase in glycolysis in both tissues and a decrease in gluconeogensis in the liver, partly explaining how both compounds reduce blood glucose levels. However, PPARδ induced increased citric acid cycle activity, and mobilized free fatty acids, ketone bodies and triglycerides and activity in the linoleic and α-linolenic acid essentialfatty-acid pathways. Thus, PPARδ has a profound role in hepatic and skeletal muscle fatty-acid oxidation.
One problem with metabolomic approaches that rely on taking tissue samples at set time points is that it is difficult to talk about fluxes from the total concentration changes that are measured using this approach. To address this, Roberts and colleagues [42] used 3T3-L1 adipocytes in conjunction with stable isotope labeling, as well as an in vivo study in the obob mouse, to investigate the role of PPARδ in adipocytes (Figure 3). Stimulation of the receptor using a highly selective agonist induced increased oxygen consumption in 3T3-L1 adipocytes, and in both adipose tissue and cultured cells there were increases in fatty acid β-oxidation, tricarboxylic acid cycle rate and oxidation of extracellular branch chain amino acids. This highlights the potential of PPARδ agonists to act as both anti-diabetes and anti-obesity treatments, addressing both these consequences of dyslipidemia.
PPARδ agonists have also been studied in moderately overweight humans, using a combination of stable isotope techniques to monitor changes in flux and classic physiology measures [43]. Treatment with GW501516, a selective PPARδ agonist, resulted in a 30% reduction in plasma triglycerides, 20% reduction in liver fat and 30% reduction in pro-inflammatory isoprostanes in the urine. The stable isotope measurements were able to show an increase in fatty-acid oxidation, thus demonstrating that PPARδ was exerting its beneficial effects by upregulating β-oxidation, with skeletal muscle being a major target for the drug action.

Metabolomic studies of the systems that interact with the PPAR signaling pathways
PPARγ co-activator (PGC)1α and PGC1β are transcriptional co-activators that interact with PPARγ and, in turn, allow interaction with the multiple proteins involved in the regulation of cellular metabolism, including cAMPresponse-element-binding protein (CREB) and nuclear respiratory factors (NRFs). PGC1α has been implicated in mitochondrial biogenesis, muscle fiber-type determination and adaptation to endurance exercise. While the PGC1 co-activators have been studied in depth in terms of physiology, they have been the focus of much less attention in terms of metabolomics and lipidomics. However, a recent study by Gurung and colleagues [44] provided intriguing evidence that PGC1β is involved in cardiac arrhythmia. PGC1β(-/-) mice had increased concentrations of pro-arrhythmic and pro-inflammatory lipids, including lysophosphatidylcholine. The pro-arrhythmogenic phenotype was confirmed by Langendorff heart perfusion, with the hearts showing ventricular tachycardia. In turn this was related to altered calcium handling in the heart tissue of PGC1β(-/-) mice.
The PPARs play a critical physiological role as lipid sensors and regulators of lipid metabolism. Poly unsaturated fatty acids (PUFAs) such as arachidonic acid and docosahexaenoic acid (DHA) are natural ligands with medium affinity for the PPARs and subsequently induce the heterodimerization of the PPAR receptors with the RXR [45,46], linking eicosanoid and related lipid mediator signaling to the PPAR systems. Once the PPAR receptor binds to its ligands, it heterodimerizes and activates the subsequent complex to induce genes for various functions such as adipocyte proliferation [47], lipid homeostasis in cardiomyocytes and skeletal muscle, stem cell proliferation [48] and the regulation of inflammation [49]. There are a number of proposed potential molecular mechanisms for the anti-inflammatory effects of PPARs. For example, PPARα activation regulates inflammation by inhibiting nuclear factor-κB (NF-κB) activation and consequently the inhibition of cyclo oxygenase (COX), which modulates the production of inflammatory lipid mediators. The RXR/PPAR complex can also inhibit expression of other inflammatory mediators, including inducible nitric oxide synthase (iNOS), tumor necrosis factor-α (TNF-α), interferon-γ (IFN-γ) and interleukin-12 (IL-12) [50].
Although the identities of high-affinity natural ligands for the PPARs are currently unknown, there is evidence that fatty acids and their oxidized metabolites can serve as activating ligands for PPARs. In particular, COX-related metabolites have demonstrated their physiological function via binding to PPARs. For instance, prosta glandin D2 (PGD 2 ) metabolites, such as 15-deoxy Δ12,14 prosta glandin J2 (PGJ 2 ), promote adipocyte differentiation via activation of PPARγ [51], whereas prostacyclin (PGI 2 ) activates PPARδ in pregnant mouse uterus during natural implantation of the placenta [52]. In addition, a number of studies have indicated that PPAR signaling is activated by lipoxygenase (LOX)-derived metabolites. For instance, 8(S)-hydroxyeicosatetraenoic acid (8(S)-HETE) has been reported to be an effective activator of PPARα [53]. Linoleicacid-oxidized metabolites, 9-HODE and 13-HODE, have been reported to activate PPARγ-dependent gene expression in monocytes. In another study, PPARγ and 15-LOX were coordinately induced by IL-4 in macrophages, and it has been suggested that the 15-LOXderived metabolites act as endogenous ligands of PPARγ to mediate transcriptional induction of the CD36 gene in these cells [54,55]. Lipidomic approaches are beginning to piece together how these different signaling pathways interact with one another. For example, Bazan and co-workers [56] have recently reported DHA-derived neuroprotectin (NPD1) as a potential PPARγ ligand and demonstrated that NPD1 anti-amyloidogenic bioactivity is mediated via activation of the PPARγ. Kim and co-workers [57] have described a novel metabolomic approach aimed at profiling metaboliteprotein interactions on columns, to look at lipids that bind to PPARγ and PPARα. Following a pull-down assay approach, they ran crude lipid extracts down the columns containing immobilized PPARγ and PPARα, identifying arachidonic acid (C20:4), oleate (C18:1), linoleic acid (C18:2), palmitoleic acid (C16:2) and C18:1 lysophos phatidic acid binding to PPARγ and arachidonic acid, linoleic acid and oleic acid binding to PPARα. Although naturally occurring lipid mediators that activate PPARs have been identified in vitro, establishing the physiological relevance of many of these substances as regulators of PPARs in vivo is challenging, due to their low concentrations and relatively low affinities to PPARs. However, advancements in analytic techniques and profiling of these metabolites could overcome the difficulties of detecting these ligands in vivo and could advance our understanding of their mechanisms of action. Consequently, due to its specificity and sensitivity, lipidomics analysis using LC-MS/MS is currently the most powerful tool for the analysis of lipid mediators (ms/ms: tandem mass spectrometry where mass spectrometry is performed on the fragments of the parent ion) [58].
PPARδ signaling is also intricately linked with lipid mediator signaling, and this has been investigated in the pancreas using high-resolution LC-MS. Cohen and colleagues [59] observed that high glucose levels in the media of INS-1E β-cells activated cPLA(2), releasing arachidonic and linoleic acid from phospholipids in the cell membrane. The high glucose and PPARδ activation both stimulated the generation of 4-hyroxy-2E-nonenal, an endogenous ligand for PPARδ, which in turn amplifies insulin secretion in β-pancreatic cells.

Metabolomic studies of PPAR-induced cancer
The PPARs have a long and complex history in terms of cancer biology. Agonists for all three receptors have been linked to increased relative risks of certain cancers in carcinogenicity studies in rodents. However, for many of these agonists, the increased risks for certain types of cancer are offset by a reduction in the relative risk for other types of cancers. However, one area that has received much attention from those using metabolomics is the role peroxisome proliferation plays in PPARαinduced cancer.
Peroxisome proliferators (PPs), many of them substrates of PPARα, are non-DNA reactive, non-genotoxic (epigenetic) carcinogens. They cause peroxisome proliferation, hyperplasia and ultimately hepatocarcino ma in rats and mice. As described earlier, the property of PPs to cause hepatocellular carcinoma eventually led to the discovery of PPARα [16]. The mechanisms underlying the carcinogenesis are not fully understood, although the hepatocarcinogenicity of PPs provides one of the beststudied [60][61][62][63][64] and perhaps the most informative system for understanding non-genotoxic mechanisms, as well as key aspects of general metabolism.
Many efforts in the pharmaceutical industry have focused on the early detection of carcinogenic potential, in an attempt to replace the costly, long-term rodent bioassays and histopathology for the detection of PPinduced non-genotoxic carcinogenicity. Non-invasive tech niques, such as urinary metabolite profiling, are particularly attractive in terms of low costs for both time and animals. Studies by Delaney and colleagues describe urinary measurements of N-methylnicotinamide (NMN) and N-methyl-4-pyridone-3-carboxamide (4PY) (endproducts of the tryptophan-NAD + pathway) as being corre lated with PP as measured by electron microscopy [65]. Urinary metabolites were determined by highresolution 1 H NMR spectroscopy and NMN, 2PY and 4PY were determined using high-pressure liquid chro mato graphy (HPLC)-MS/MS for accurate quantification. A follow-up report by Connor and co-workers [66] further characterized the NMR spectra obtained in the same animal study.
The discovery that peroxisome proliferation leads to epigenetic events gave grounds for the investigation of Sadenosylmethionine (SAMe) and S-adenosyl homo cysteine (SAH) [67], as well as other metabolites in the methionine cycle, and even phospholipid metabolism following peroxisomal proliferation [68]. The increased production of H 2 O 2 following peroxisome proliferation evidently links oxidative stress to non-genotoxic carcinogenesis and to metabolites of oxidative stress, including glutathione (GSH)/oxidized glutathione (GSSG) [69], epoxides, uric acid [70], peroxisomal and mitochondrial ω-oxidation (carnitine and acyl-carnitine species) and other products of lipid oxidation [71]. It has been well documented that PPs mediate liver P450s, which carry out ω-hydroxylation of lauric acid (and other lipids) [72]. Metabolism of lauric acid has since been of interest as a useful tool for the detection of possible toxicity associated ω-hydroxylation [73,74].
Metabolomics has also been used to investigate the role PPARδ plays in breast cancer [70]. In this study, to induce tumors 3-phosphoinositide-dependent kinase-1 (PDK1) was expressed in mammary epithelium. The expression of the gene induced increased expression of PPARδ. In addition, treatment of both the transgenic mouse and its wild-type control with GW501516 (a selective PPARδ agonist) increased rates of tumor formation, particularly in the transgenic animal. LC-MS-based metabolomics of the methanol extract of breast tumors demonstrated that the agonist induced a distinctive metabolomic profile in breast tumors [70]. Intriguingly, this was associated with increased phospholipid and fatty acid accumulation, in contrast with the action of PPARδ agonists in healthy liver, skeletal muscle and adipose tissue, where increased β-oxidation appears central.

Conclusions
Because of the diverse roles of the PPAR transcription factors, and in particular their key part in regulating systemic metabolism, this area of biology has been a useful 'sandbox' for researchers to develop metabolomic techniques. In the opinion of the authors, there is adequate evidence that metabolomics has an integral part to play alongside more established physiology and transcriptomic techniques. This has included monitoring the role individual PPARs play in regulating systemic metabolism and how the normal aging process interacts with these receptors. The use of metabolomics has also opened up novel areas of research, such as the role PPARα plays in alcohol-induced steatohepatitis or the induction of liver cancer by non-genotoxic carcinogenicity. However, the interest in PPARs as a direct target for drugs has somewhat diminished recently, due to concerns of unwanted side effects. The capability of metabolomics to monitor some of the changes in lipid mediators may provide alternative targets either to bypass the PPAR receptor itself, or to provide a more controlled manner of stimulation for these key regulators of cellular metabolism. Metabolomics-derived biomarkers of drug efficacy may also allow the deployment of drugs as part of a personalized medicine approach and ensure that some of the side effects are avoided in their future use.