Abstract

Peroxisome proliferator-activated receptorδ (PPARδ) regulates a multiplicity of physiological processes associated with glucose and lipid metabolism, inflammation, and proliferation. One or more of these processes likely create risk factors associated with the ability of PPARδ agonists to promote tumorigenesis in some organs. In the present study, we describe a new gastric tumor mouse model that is dependent on the potent and highly selective PPARδ agonist GW501516 following carcinogen administration. The progression of gastric tumorigenesis was rapid as determined by magnetic resonance imaging and resulted in highly metastatic squamous cell carcinomas of the forestomach within two months. Tumorigenesis was associated with gene expression signatures indicative of cell adhesion, invasion, inflammation, and metabolism. Increased PPARδ expression in tumors correlated with increased PDK1, Akt, β-catenin, and S100A9 expression. The rapid development of metastatic gastric tumors in this model will be useful for evaluating preventive and therapeutic interventions in this disease.

1. Introduction

Gastric cancer is the second leading cause of cancer-related death globally, particularly in developing countries, due to its metastatic nature at the time of diagnosis [1, 2]. Among the many risk factors for gastric cancer are diet, smoking, and inflammation associated with H. pylori infection [35] that are likely exacerbated in patients with proinflammatory gene polymorphisms [6].

PPARs are ligand-dependent nuclear receptors that regulate expression of a multiplicity of genes associated with metabolic disorders, such as type II diabetes and lipodystrophies [7, 8]. PPARs consist of the α, γ and δ isotypes that regulate not only glucose and lipid metabolism, but also proliferation, inflammation, and angiogenesis [913]. PPARδ expression is increased in breast, colon, and head and neck cancers [9, 1417] and is associated with a more aggressive phenotype in breast cancer cells [18]. The selective PPARδ agonist GW501516 acts as a tumor promoter in mammary carcinogenesis [19] and colon tumorigenesis [15, 20, 21], whereas disruption of PPARδ expression blocks mammary [22] and colon tumorigenesis [23, 24]. PPARδ is regulated by risk factors implicated in cancer progression, including K-Ras [25], Wnt [26], and Pges/Cox2 [22, 27], and is associated with activation of angiogenesis [20, 2830]. PPARδ regulates proinflammatory signaling through NFκB and IL-1β [31], and is activated by metabolites of arachidonic acid metabolism that serve as PPAR ligands [14, 29, 32]. Notwithstanding these results, there are several reports showing that azoxymethane-induced colon carcinogenesis is increased in PPARδ nullizygous mice [33, 34] and by the selective PPARδ ligand GW0742 [35] (reviewed in [36]). Differences between the various nullizygous models that may account for some of these disparities have been discussed [9, 21].

Here we demonstrate that activation of PPARδ by a selective agonist, GW501516, rapidly induces highly metastatic gastric tumors following carcinogen administration, which expressed a markedly increased inflammatory gene expression signature. These findings suggest a crucial role for PPARδ in the progression of this disease.

2. Materials and Methods

2.1. Materials

GW501516 was synthesized as previously described [37] and was provided by the Chemoprevention Branch, National Cancer Institute, NIH, Bethesda, MD.

2.2. Carcinogenesis

DMBA (Sigma-Aldrich, St Louis, MO) was dissolved in cottonseed oil at a concentration of 10 mg/ml, and six week-old female FVB mice were administered 4 weekly doses of 1 mg DMBA by gavage. Animals were fed a standard diet or a diet supplemented with 0.005% (w/w) GW501516 beginning one day after the last dose of DMBA as previously described [19]. Mice were euthanized when gastric tumors reached ≥2cm3 as determined by MRI imaging or if mice became moribund. All protocols were approved by the Georgetown University Animal Care and Use Committee.

2.3. Magnetic Resonance Imaging

MRI was performed on a 7.0 Tesla Bruker horizontal spectrometer/imager with a 20 cm bore equipped with 100 gauss/cm microimaging gradients and run by Paravision 4.0 software in the Preclinical Imaging Research Laboratory, Lombardi Comprehensive Cancer Center. Mice were anesthetized using 1.5% isoflurane and 30% nitrous oxide, positioned in a custom-made stereotaxic animal holder with temperature and respiration control, and imaged in a 72 mm birdcage radiofrequency volume coil as previously described [38]. The imaging protocol used was a T2-weighted RARE (rapid acquisition with refocused echos) two-dimensional sequence with the following parameters: Matrix: 256, TR: 5822 msec, TE: 36 msec, number of averages: 4, RARE factor: 8, FOV: 4.0 cm, number of slices: 4, slice thickness: 0.5 mm, and with respiratory gating to account for breathing movement.

2.4. Histopathology and Immunohistochemistry (IHC)

Stomach and tumors were excised, and formalin-fixed, paraffin-embedded sections were prepared for H&E staining and IHC. Antigen retrieval was carried out by incubation of tissue sections in 10 mM sodium citrate buffer (pH 6.0) for 20 min at a subboiling temperature in an electric steamer as previously described [19, 39]. Endogenous peroxidase activity was quenched with 3% hydrogen peroxide for 10 min and incubated for 30 min with blocking solution (10% goat serum in Tris-buffered saline), followed by incubation overnight at 4°C with the appropriate primary antibody diluted in blocking solution. Biotin-conjugated secondary antibodies were diluted in TBS containing 0.1% Tween-20 and incubated for 30 min at room temperature using the ABC Vectastain (Vector Laboratories) detection system and diaminobenzidine (Pierce), and slides were counterstained with Harris-modified hematoxylin (Fisher Scientific) and dehydrated and mounted in Permount (Fisher Scientific). The following antibodies and their dilutions were used: anti-CK14 (1 : 200, ms-115-P1, Neomarkers), anti-CK18 (1 : 200, sc-7197, Santa Cruz Biotechnology), PDK1 (1 : 50, sc-17765, Santa Cruz Biotechnology), anti-pS473 Akt (1 : 200, 40515, Cell Signaling Biotechnology), anti-β-catenin (1 : 50, sc-7963, Santa-Cruz Biotechnology) and anti-S100a9 (1 : 50, sc-65580, Santa Cruz Biotechnology).

2.5. Gene Microarray Analysis

Three groups of tissues were analyzed: (1) gastric tumors, (2) forestomach (nonglandular) from GW501516-treated mice, and (3) forestomach from DMBA-treated mice. Tissue was excised, washed in phosphate-buffered saline, and stored in RNAlater (Ambion) at 20°C until RNA extraction. Tissue was snap-frozen in liquid nitrogen, pulverized in a mortar and pestle, and RNA extracted using an RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. RNA purity was assessed by an A260/A280 ratio of ≥1.9, and by the integrity of 18S and 28S rRNA using an Agilent microfluidic chip. Array analysis was carried out on cRNA prepared from equal amounts of RNA (1 μg) pooled from 5 mice per group as previously described [19, 40]. Biotin-labeled cRNA was fragmented at 94°C for 35 min and hybridized overnight to an Affymetrix mouse 430A 2.0 GeneChip representing approximately 14,000 annotated mouse genes. GeneChips were scanned with an Agilent Gene Array scanner, and grid alignment and raw data generation with the Affymetrix GeneChip Operating software 1.1. A noise value ( ) based on the variance of low-intensity probe cells was used to calculate a minimum threshold for each GeneChip. Samples were averaged and data refined by eliminating genes with signal intensities <300 in both comparison groups, and heat maps were generated from ≥3-fold changes in gene expression normalized to untreated forestomach using unsupervised hierarchical cluster analysis as previously described [41]. Gene ontology analysis utilized Ariadne Pathway Studio version 7.1

2.6. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

Total RNA was extracted as described above, and equal amounts of RNA (1 μg) were pooled from each group (five samples per group) and reverse transcribed with the Omniscript RT kit (Qiagen) in a total volume of 20 μL as previously described [19, 39]. PCR was performed in triplicate using an ABI-Prism 7700 (Applied Biosystems, Foster City, CA) with SYBRGreen I detection (Qiagen) according to the manufacturer's protocol. Amplification using the appropriate primers (see Table S1 in supplementary material available at http://dx.doi.org/10.1155/2010/571783) was confirmed by ethidium bromide staining of the PCR products on an agarose gel. The expression of each target gene was normalized to GAPDH and is presented as the ratio of the target gene to GADPH expression calculated using the formula, , where .

3. Results

3.1. PPARδ Agonist GW501516 Rapidly Promotes Gastric Tumorigenesis

Mice maintained on a diet supplemented with PPARδ agonist GW501516 following carcinogen administration resulted in the rapid development of gastric tumors in 12/15 animals, whereas treatment with either GW501516 or DMBA alone was not tumorigenic (Table 1). To follow the onset and progression of tumorigenesis more precisely, five mice were monitored by MRI (Figure 1(a)). Tumors were visible as early as 19 days after beginning the GW501516 diet and appeared to initiate in the forestomach (Figure 1(a)). By 50 days, tumor had filled the stomach lumen, and by 56 days it had extravasated through the gastric wall (Figure 1(a)). Gross inspection of the stomach confirmed that the tumor was confined within the stomach at day 20 but had invaded through the stomach wall forming local metastases by day 56 (Figure 1(b)). Mice showed signs of morbidity between days 63 and 70 (mean survival 67 days), where metastases were present throughout the mesentery and adjacent serosal organ surfaces including the abdominal wall (Figure 1(b)).

Primary tumors and metastases were uniformly squamous cell carcinomas with varying degrees of keratinization (Figure 2(A)). Animals fed the GW501516 diet for six months without prior DMBA treatment did not exhibit hyperplasia or dysplasia (Figure 2(B)), and DMBA treatment alone produced squamous cell hyperplasia of the forestomach without signs of dysplasia (Figures 2(A) and 2(B)). No changes occurred in the gastric mucosa resulting from DMBA and GW501516 treatment alone (data not shown), and esophageal squamous epithelium was unaffected by DMBA treatment (Figure 2(B)). Gastric tumors were positive for the squamous basal cell marker CK14, and negative for the columnar epithelial cell marker CK18 [42] (Figure 2(C)).

3.2. Differential Gene Expression in Gastric Tumors and Stomach

Tumors manifested a marked inflammatory phenotype denoted by increased expression of chemokines Cxcl-1,-2,-5,-9,-14 and Ccl-2,-3,-8, S100a8, S100a9, and S100a3 and interleukins IL-1β, IL-6, and IL-24 (Table 2, Figure 3(b), and Table S2). In addition to these changes, increased expression of prostaglandin metabolism genes Ptgs2/Cox2 and Ptges and reduced expression of PPARγ and PPARα were noted. To determine if these changes were tumor specific, gene expression was assessed in stomach tissue after treatment with either GW501516 for seven days (Figure 3(a), Table S3) or DMBA for four weeks (Table S4). GW501516 increased expression of only five genes ≥3-fold, Angptl4, Cyp2b10, Cfd/Adipsin, Adipoq and Chi3l4 and markedly reduced expression of Gast, Ccla3, Glycam1, Spp1, Serpina1a, Cela1, Cldn2, and Fabp2 (Table S3). DMBA increased expression of S100a8, S100a9, and Ccl8 4–10-fold and reduced expression of the same subset of genes as GW501516 (Table S4). Thus, changes in Gast, Ccla3, Glycam1, Spp1, Serpina1a, Cela1, Cldn2, and Fabp2 are a result of GW501516 treatment and are not tumor specific. On the other hand, DMBA produced an inflammatory response denoted by increase S100A8, Ccl8, and S100A9 although it was an order of magnitude less than in tumors. The increase in Krt6a by DMBA is consistent with increased keratinization in the stomach (Figure 2(B)) but was less pronounced than in tumors. Real-time qRT-PCR analysis generally confirmed the changes in expression of several genes, including Cldn8, Cxcl1, Cxcl5, Foxg1, S100a8, Angptl4, Cyp2b10, Vegfα and Spp1, Gast, Dkk1, Bmp3, Bmp4, PPARα, and PPARγ (Figure 3(c)).

The expression of PPARδ and factors known to be associated with its signaling were assessed in tumors and forestomach after GW501516 treatment (Figure 4). GW501516 increased nuclear localization of PPARδ in gastric squamous epithelium and tumors, in contrast to its diffuse cytoplasmic staining in untreated gastric tissue. GW501516 also elicited strong pS473Akt and pT308Akt staining in basal cells and in the submucosal layer, as well as in tumor and stromal tissue, which correlated with more intense PDK1 expression. β-Catenin was expressed in the nuclei of basal squamous epithelial cells and was not altered by GW501516 treatment, whereas tumors expressed increased β-catenin at cellular junctions. S100a9 was absent in untreated gastric epithelium but was expressed in endothelial and epithelial cells from GW501516-treated mice. Tumors expressed S100a9 in a diffuse pattern, with strong expression in blood vessels and adjacent epithelial cells.

4. Discussion

The present study describes a new model of metastatic gastric cancer that is dependent on the tumor promoting activity of PPARδ agonist GW510516 following carcinogen administration. In contrast to a previous study reporting a low percentage of squamous cell carcinomas of the forestomach by DMBA [43], our DMBA regimen produced only forestomach hyperplasia without signs of dysplasia up to five months after treatment (Figure 2(B)). This suggests a high sensitivity of mouse forestomach squamous epithelium to dysplasia, and the predilection of GW501516 to promote tumors of this histotype [19]. This model differs from N-methyl-N-nitrosourea-induced gastric tumors in wild-type and APCMin transgenic mice [44, 45] in that it is dependent on both DMBA-induced mutagenesis and the tumor-promoting effects of GW501516. A feature of this model is its short latency of approximately three weeks in comparison to 10 to 20 weeks for NMU-treated wild-type and APCMin mice. An important histopathological distinction, and perhaps disadvantage of the GW501516 tumor model, is that it produces squamous cell carcinomas from the nonglandular forestomach rather than adenocarcinomas from the glandular tissue that comprises the majority of human gastric cancer [46]. Since this model was dependent on the selective PPARδ agonist GW501516 [47], it is important to note that the dose of GW501516 used in the present and previous studies [19] is equivalent to daily oral doses of 3–10 mg/kg that were previously shown to specifically enhance PPARδ-dependent fatty acid oxidation in mice [48]. In addition, PPARδ agonist GW7042, which is almost identical to GW501516 in structure, potency, and specificity, was inactive in inducing gene expression in PPARδ knockout mice [49], suggesting that the tumor promoting effects of GW501516 and GW7042 are not due to off-target effects.

Tumors induced by GW501516 exhibited a distinct inflammatory gene expression signature comprised predominantly of chemokine, MMP, and S100 genes (Table S2). This was unexpected in view of the lack of a similar signature after treatment with GW501516 (Table S3), and the fact that GW501516 induces an anti-inflammatory response in macrophages [50] and protects the heart against oxidative stress [51, 52]. Gene ontology analysis of gene expression in the gastric tumors indicates that PPARδ, MMP12, MMP13, Cxcl1, Cxcl5, S100A8, and S100A9 share both common and disparate pathway interactions that likely contributed to the tumorigenic phenotype (Figure 5). PPARδ is associated with activation of genes related to proliferation (EGFR, Akt1) and adhesion (Itgb2), whereas S100A9 is associated with angiogenesis (Fgf2) and inflammation (Ager). Cxcl1 activates proliferation (Mapk3, Mapk14, and Akt1), angiogenesis (Fgf2), and invasion (MMP2, MMP9), and Cxcl5 activates other chemokines (Cxcl1, Cxcl2, and Cxcl3) and Ptgs2/Cox2. This scheme reiterates the ability of S100A8 and S100A9 to act as ligands for Ager (advanced glycation end-product receptor), which mediates acute and chronic inflammation, tumor development, and metastasis [53, 54]. This paradigm is also consistent with GW501516-induced activation of Ptges and Ptgs2/Cox-2 expression [55, 56], which initiate the production of prostacyclins [14] and arachidonic acid metabolites [29] that serve as PPARδ ligands [57]. Since Cox-2 inhibitors reduce inflammation-related gastrointestinal carcinogenesis [58], and overexpression or deletion of Ptgs2 increases or suppresses tumorigenesis, respectively [59, 60]; this suggests cooperativity between PPARδ and inflammatory signaling pathways in gastric tumorigenesis. The ability of PPARδ to have an anti-inflammatory effect in normal cells [51, 52] and a proinflammatory effect in tumors is reminiscent of the dual roles of TGF-β in tumor cells [61, 62]. TGFβ can function as a proinflammatory cytokine to activate S100A8 and S100A9 expression in the presence of activated Ras [63] but acts as a repressor of inflammation-induced PPARδ expression in normal cells [64, 65]. The increased expression of the PPARδ target gene, Agptl4, the TGFβ-activated genes Runx1 and Runx2, and S100A8 and S100A9 in the gastric tumors indeed suggests a duality of function of both PPARδ and TGFβ signaling in gastric tumorigenesis.

PPARδ is ubiquitously expressed in gastrointestinal tissue and gastric tumors [27], and GW501516 elicited increased PPARδ nuclear staining and elevated pAkt in gastric epithelium and tumors. PPARδ-dependant activation of Akt is required for the growth-promoting and antiapoptotic effects of PPARδ [6668], as shown by the delayed wound-healing response of PPARδ-deficient keratinocytes [57, 69]. Enhanced Krt6a and Krt16 expression in tumors further suggests that PPARδ plays an important role in gastric squamous cell differentiation and tissue renewal.

Tumors also exhibited reduced PPARγ and PPARα expression that may have resulted, in part, from the negative regulation of PPARγ by PPARδ [70]. PPARγ suppresses the growth and invasion of human colon [71] and gastric [72, 73] and esophageal carcinoma cells [74], and both PPARγ [75] and PPARα have anti-inflammatory actions [76]. Thus, reduction of PPARα and PPARγ expression may be an additional mechanism for facilitating the proinflammatory and tumor-promoting effects of GW501516.

In summary, we describe a rapidly developing metastatic gastric cancer model dependent on the tumor-promoting effects of GW501516 following carcinogen treatment, which suggests a proinflammatory switch in PPARδ function. This animal model will therefore be useful to delineate the role of PPARδ in tumor initiation and progression and as a possible target for early intervention.

Acknowledgments

This work was supported by Grant no. R01 CA111482 and Contract no. N01 CN43309 from the National Institutes of Health, Bethesda, MD. This investigation was conducted using the Animal Research, Histopathology and Tissue, Genomics and Epigenomics, and Preclinical Imaging Shared Resources supported by Research Facilities Improvement Grant no. C06 RR14567 from the National Center for Research Facilities, and by Cancer Center Support Grant no. 1P30-CA-51008 from the National Cancer Institute. The authors thank Yi Chien Lee for his help with the MRI.

Supplementary Materials

Supplementary material contains qRT-PCR (Table S1) primers and differential gene expression data in gastric tumors (Table S2), and in the stomach of GW501516-treated (Table S3) and DMBA-treated (Table S4) animals.

  1. Supplementary Table S1
  2. Supplementary Table S2
  3. Supplementary Table S3
  4. Supplementary Table S4