Skip to main content

2-amino-1-methyl-6-phenylimidazo(4,5-b) pyridine (PhIP) induces gene expression changes in JAK/STAT and MAPK pathways related to inflammation, diabetes and cancer

Abstract

Background

2-amino-1-methyl-6-phenylimidazo(4,5-b)pyridine (PhIP), a heterocyclic aromatic amine (HCA) formed in meat that is cooked at high temperatures and then ingested, can potentially be retained in human adipose tissues.

Methods

To determine if PhIP is bioactive in the adipocyte, we exposed a human adipocyte cell line,HepG2 and Caco-2 cells to low dose PhIP. Uptake and retention of PhIP was determined and cytotoxicity was assessed by the TUNEL assay. Relative expression of PhIP-activating genes (CYP1A1, CYP1A2, SULT1A1 and UGT1A1) was determined by RT-PCR and global expression changes were also examined.

Results

The percent retention of 0.1 μCi [14C]-PhIP over a 24 h period was significantly higher in the adipocyte than the HepG2 (p = 0.0001) and Caco-2 (p = 0.0007) cell lines. Cytotoxicity rates were 14.4 and 2.6 % higher compared to controls in Caco-2 and HepG2 cells (p < 0.001 and 0.054, respectively); no significant differences were detected in adipocyte cells (p = 0.18). Caco-2 and HepG2 cells, respectively, had significantly higher basal expression of CYP1A1 (p = 0.001, p = 0.003), SULT1A1 (p = 0.04, p < 0.001) and UGT1A1 (p < 0.001, p = 0.01) compared to the adipocyte. Exposure to 5nM PhIP did not significantly induce expression of these genes in any of the cell lines. Global gene expression analysis of mature adipocytes exposed to 5nM PhIP for 72 h resulted in statistically significant changes in 8 genes (ANGPTL2, CD14, CIDEA, EGR1, FOS, IGFBP5, PALM and PSAT1). Gene-gene interaction and pathway analysis indicates that PhIP modulates genes controlled by the STAT3 transcriptional factor and initiates leptin signaling via the JAK/STAT and MAPK pathway cascades. Early growth response 1 (EGR1) and prostaglandin synthase 2 (COX-2) were down-regulated via c-Fos, while insulin binding protein 5 (IBP5) was up regulated. Expression of transcription factors (ANGPTL2, HP, LEP, SAA1, SAA2), genes related to inflammation (SAA1, LEP), diabetes (IGFBP5) and cancer risk (SAA2) were also elevated upon exposure to 5 nM PhIP..

Conclusions

PhIP mediates gene expression changes within the adipocyte, and the pathways most affected are related to cancer and other chronic diseases. Further studies are needed on the relationship between dietary carcinogens such as PhIP with cancer, obesity and diabetes.

Background

Heterocyclic aromatic amines (HCAs) are dietary chemicals formed when muscle meat is cooked at high temperatures, such as broiling, frying or grilling. Lipids make up 70 % of white adipose tissue in humans and this endocrine organ represents a reservoir for many lipophilic contaminants [1, 2]. The most abundant lipophilic HCA in cooked meat is 2-amino-1-methyl-6-phenylimidazo[4,5-b] pyridine (PhIP) [3]. PhIP is a procarcinogen that is metabolically activated by cytochromes P450 (CYP) 1A1 and 1A2 to form the metabolite, N-hydroxy-PhIP (N 2-OH-PhIP) [46]. N 2-OH-PhIP can then undergo detoxification reactions (i.e. glucuronidation) or reactions that generate reactive metabolites (i.e. sulfation) [79]. N 2–OH-PhIP metabolites induce chromosomal aberrations [10], form PhIP-DNA adducts in colon mucosa [10], and causes colon and mammary carcinomas in animal models [1116]. The current consensus in the literature is that consumption of well-done meat is associated with both elevated cancer risk and mortality [17]. However, case control studies examining meat consumption, HCA content and cancer risk lack concordance. Conflicting reports in the literature are almost certainly influenced by the difficulty of quantifying exposures, since HCA content can vary 100-fold depending on cooking temperature [18]. Additionally, the interaction between adiposity, HCA exposure and cancer risk is poorly understood.

Obesity is overtaking under-nutrition and infectious diseases as the most significant causes of ill-health [19]. Obesity is also predicted to have the major influence on cancer incidence, since overweight and obesity are recognized risk factors for various cancers [20]. In addition to its role as an endocrine organ, adipose tissue acts as a reservoir for a large fraction of lipophilic contaminants, including persistent organic pollutants, or POP [21]. A study of the distribution of radiolabeled PhIP in mice demonstrated that after administration of a single dietary equivalent dose of PhIP, the highest concentrations were detected in the stomach, intestine and adipose tissue [22]. Others have reported that human mammary lipid can be mutagenic, and although the identity of the mutagens was not confirmed, the lipid was isolated using a solid-phase tandem extraction procedure originally developed for extraction of heterocyclic amines. The presence of substances in extracts of essentially normal human mammary lipid induced point mutation in bacteria in the presence of metabolizing enzymes and chromosomal damage (signaled by micronucleus formation) in metabolically competent human cells [23].

In the current study, we compared the uptake and retention of PhIP into a human adipocyte cell line (differentiated from human mesenchymal stem cells, (hMSC)), in comparison with a hepatocellular cell line (HepG2) and a colorectal cell line (Caco-2) in order to determine if there were differences. We also compared the ability of PhIP to induce DNA damage and subsequent cell death in these three cell lines. We examined the ability of PhIP to induce gene expression changes in PhIP metabolizing enzymes in each of the cell lines. Finally, global changes in gene expression in the PhIP-exposed adipocyte cell line were also explored.

Methods

Cell culture and differentiation

The identity of all cell lines was confirmed by submission to American Type Culture Collection (ATCC) (Manassas, VA) for analysis by short tandem repeat methodology on December 9, 2015. Human Mesenchymal Stem Cells (hMSC) were obtained from AllCells (Emeryville, CA). The hMSCs were expanded in growth medium (GM) consisting of alpha Minimum Essential Medium (αMEM), supplemented with 16.5 % fetal bovine serum (FBS) (Tissue Cultures Biologicals, Long Beach, CA), 2 mM L-glutamine (Sigma-Aldrich, St. Louis, MO), and penicillin/streptomycin at a concentration of 100 μg/ml each (Thermo Fisher Scientific, Watham, MA), in a humidified incubator at 37 °C, with 5 % CO2. Cells were seeded 2.0 × 105 cells/well in 6 well plates (for RNA isolation) or 2.0 × 104 /well in 24 well plates (for uptake/retention experiments). Fat differentiation medium (FDM), composed of GM supplemented with 0.5 μM dexamethasone, 0.5 μM isobutylmethylxanthine and 50 μM indomethacin (Sigma-Aldrich, St. Louis, MO), was used to initiate adipocyte differentiation. Fresh FDM was applied every 3-4 days. Differentiation reached a lipid maximum around day 21 at which time PhIP treatments began. Correlation of differentiated vs undifferentiated hMSC was determined by AdipoRed™ Assay Reagent (Lonza Walkersville, Inc. Walkersville, MD) and Trypan Blue (Sigma-Aldrich, St. Louis, MO) readings using a SpectraMax M5 Molecular Devices (VWR, Atlanta, GA) (R 2 = 0.9826).

HepG2 and Caco-2 cell lines were purchased from ATCC and maintained following the protocol described by ATCC at 37 °C in an incubator containing 5 % CO2. HepG2 and Caco-2 cells were cultured in EMEM + L-glutamine (ATCC) supplemented with 10 % FBS and EMEM (ATCC) 20 % FBS, respectively. Penicillin/streptomycin at a concentration of 100 μg/ml (Thermo Fisher, Watham, MA) was added to both media.

Uptake of [14C]-PhIP into human cell lines

All cell lines were treated with 5 μM unlabeled PhIP plus 0.1 μCi [14C]-PhIP (Toronto Research Chemicals, Toronto, Ontario, Canada) in 0.5 ml FDM using dimethyl sulfoxide (DMSO) (Sigma-Aldrich, St. Louis, MO) as vehicle for all PhIP exposures. The PhIP solution was diluted in culture medium and cells were exposed for 0, 5, 7, 15, 60 and 300 s (s). The cells were then lysed with 100 μl RNAqueous lysis buffer solution (Ambion by Life Technologies Corp., Austin, TX) and the lysate was transferred to a Polyethylene Terephthalate (PET) 24 well sample counting plate containing 300 μl OptiPhase SuperMix scintillation cocktail (Perkins Elmer, Waltham, MA). Sample scintillation counting in corrected cycles per minute (CCPM) was read using the MicroBeta Trilux 1450 LSC and Luminescence Counter (Perkin Elmer, Waltham, MA). Samples were run in triplicate and comparisons between cell types were assessed by t-test.

Release of [14C]-PhIP from human cell lines

Cell lines were plated as described for the uptake studies. Each well was treated with 5 μM unlabeled PhIP plus 0.1 μCi 14C-PhIP in FDM using DMSO as vehicle and allowed to reach equilibrium for 30 min. The medium (100 μl) was then transferred to PET 24 well sample counting plates containing 300 μl OptiPhase SuperMix scintillation cocktail for counting. Fresh medium (0.5 ml) was applied to the cells, and then removed after 0, 0.5, 1, 2, 5, 15, 60 and 1440 min of exposure. 100 μl aliquots of the fresh medium were transferred to scintillation plates. Cells were then lysed with 100 μl lysis buffer and the lysate was transferred to scintillation plates and CCPM were recorded.

TUNEL assay

Cells were plated at 1.0 × 104 per well onto Lab-Tek 8 well Chamber Slides (Nalge Nunc International, Rochester, NY) and allowed to grow overnight in a humidified incubator at 37 °C, with 5 % CO2. Adipocytes were allowed to differentiate using medium described above for 21 days. Quadruplicate cell samples were exposed to vehicle or PhIP (5 μM) up to 72 h. TUNEL assays were performed using the In Situ Cell Death Detection Kit (Roche Diagnostics, Indianapolis, IN) following the manufacturer’s protocol. A set of short band filters with respective emission and excitation wavelengths were used for the green spectrum (FITC) and blue spectrum (DAPI) for the TUNEL-positive objects and nuclei, respectively.

Effect of PhIP exposure on expression of PhIP-metabolizing enzymes in human cell lines

Changes in the expression of PhIP- metabolizing enzymes were assessed in HepG2, Caco-2 and hMSC-derived adipocytes seeded at 2 × 105 cells/well in 6- well plates. Medium containing vehicle (DMSO) or PhIP at a final concentration of 5 nM was applied, and the cells were returned to the incubator. All assays were performed in triplicate. Medium with vehicle or PhIP was replaced daily. After 72 h, the cells were washed in ice-cold Dulbecco’s phosphate-buffered saline (DPBS) (Cellgro, Thermo Fisher Scientific, Watham, MA), lysed and total RNA was isolated using the RNAqueous Kit (Ambion). RNA quality and quantity was measured using a NanoDropTM 8000 Spectrophotometer (Thermo Fisher Scientific, Watham, MA). The resulting RNA absorption 260/280 nm ratio was 2.0, indicating that the RNA prep was pure. cDNA was generated using 1 μg RNA with High Capacity RNA to cDNA Reverse Transcription Reagent (Life Technologies Corp., Austin, TX) and thermo cycler conditions of 37 °C for 60 min, followed by 95 °C for 5 mins using the GeneAMP 7900 PCR System (Applied Biosystems by Life Technologies Corp., Austin, TX). The resulting cDNA was diluted 1:20 prior to use in real-time RT-PCR.

Real-time PCR of target genes known to metabolize PhIP was performed on the 7900HT cycler (Applied Biosystems, by Life Technologies Corp., Austin, TX) under the following conditions: 95 °C for 10 min, 40 cycles of 95 °C for 1 s, followed by 60 °C for 20 s. A dissociation cycle of 95 °C for 15 s, 60 °C for 15 s and 95 °C for 15 s was performed at the end of each PCR reaction to ensure the amplicon specificity. Fast SYBR Green PCR Master Mix (Applied Biosystems, by Life Technologies Corp., Austin, TX) was used along with 5 μM of primer purchased from Integrated DNA Technologies (Coralville, IA). Primer sequences are available in Table 1. Expression was normalized to the endogenous control reference gene HRPT1. HRPT1 was identified as “Best Gene” by Normfinder software [24] using gene stability values compared with β-Actin and GAPDH in all three cell types over treatments. Expression was quantified using the relative gene expression 2 -∆∆CT method [25].

Table 1 Primer sequences

Global gene expression analysis

hMSC cells were differentiated and exposed to varying concentrations of PhIP for 72 h. RNA was isolated by the method described above. RNA integrity was analyzed on an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Total RNA (500 ng) was used to generate cRNA for hybridization with the Illumina TotalPrep™ RNA Amplification Kit (Life Technologies) according to the manufacturer’s instructions. Microarray analysis was performed at the University of Arkansas for Medical Sciences Genomics Core. 750 ng cRNA was used for hybridization to Illumina HumanHT-12_v4_BeadChip for Gene Expression according to the manufacturer’s protocol. BeadChips were scanned with the Illumina iScan system, and data collected using GenomeStudio software v2011.1. Data was generated using the gene expression module of Illumina Genome Studio software after quality control analysis. The expression intensities were median-normalized after log transformation. Differentially expressed genes were determined using an empirical Bayes test implemented in Linear Models for Microarray Data (Limma) package. A fold-change > 1.5 and p-value < 0.05 were chosen as the cutoff criteria for screening differentially expressed genes. The p-values were adjusted for multiple testing by controlling the false discovery rate according to the Benjamini Hochberg method. All analyses were conducted using the R statistical environment (R Foundation for Statistical Computing, Vienna, Austria).

Pathway analysis

Functional profiling of differentially affected biological processes, functions, pathways, and networks for the differentially expressed genes (DEGs) were evaluated using open source tools and the commercial pathway packages Metacore (MetacoreTM, www.genego.com), DAVID and STRING databases. Differentially expressed transcripts were mapped to the top most differentially expressed pathways and process networks. For network analysis of differentially expressed genes, Metacore was employed to identify any direct interactions amongst the DEGs. STRING database was also used to explore the gene-gene interactions based on experimental evidence and predictions based on homology, text mining, and co-occurrence.

Statistical analysis

Uptake and retention data was reported as the percent of total [14C]-PhIP that remained in cell lysate at various time points. Percent means ± SD of triplicate measurements and unpaired Student t-tests were employed to identify significant differences between the cell lines. Differential gene expression was determined by the 2-∆∆CT method [24] and reported as mean ± SD of triplicate measurements. The comparisons of significance were carried out by unpaired Student’s t-test. Cytotoxicity was assessed by unpaired t-tests

Results

Uptake and retention of [14C] Up-PhIP in cell lines

Uptake of [14C]-PhIP into Caco-2, HepG2 and hMSC derived adipocyte (at 70 % differentiation) cells was rapid in all cell lines. Maximal uptake of [14C]-PhIP was apparent after 5 s of exposure in HepG2 and Caco-2 cell lines and remained unchanged over the time course of the experiment. Overall, Caco-2 and HepG2 cells did not differ significantly. Uptake of [14C]-PhIP into adipocytes was significantly higher than into Caco-2 and HepG2 cells at 60 s and 300 s, respectively, and did not appear to plateau throughout the time course of the experiment (Fig. 1a). When levels of [14C]-PhIP retained in the cell lines were examined, adipocyte cell [14C]-PhIP content was significantly higher than in Caco-2 and HepG2 cells at all of the time points (Fig. 1b).

Fig. 1
figure 1

a Uptake of 5 mM PhIP + 0.1 mM [14C]-PhIP in Caco-2, HepG2 and adipocytes. (*p < 0.05) (b) Retention of 5 mM PhIP + 0.1 mM [14C]-PhIP in Caco-2, HepG2 and adipocytes. Data presented as Mean ± SD. (*p < 0.05 across time course)

Cytotoxicity of PhIP to human cell lines

PhIP and other HCAs are genotoxic and cytotoxic at high doses and induction of cell death (apoptosis) is cell-type specific [2628]. In order to determine if PhIP had adverse effects at low concentrations and in adipocytes, we exposed Caco-2, HepG2 and adipocyte cells to 5 μM PhIP for 72 h. Irreversible cell death associated with genotoxic DNA fragmentation was assessed by the TUNEL assay. Figure 2 is a representative area of the slides used by TUNEL assay for analysis (Fig. 2a = Caco-2 cells, Fig. 2b = HepG2 cells and Fig. 2c = Adipocyte cells). Unpaired t-test revealed statistical differences from controls (p < 0.05) in the Caco-2 cells at 14.4 % (p = 0.002) and borderline significance in the HepG2 cells at 2.6 % (p = 0.054). There were no significant differences between control and PhIP-treated adipocyte cells (p = 0.18).

Fig. 2
figure 2

TUNEL assay of Vehicle or 5 μM PhIP exposure for 72 h cytotoxic damage in (a) Caco-2 (14.4 % higher than controls, p = 0.0019) (b) HepG2 (2.6 % higher than controls p = 0.0536) and (c) adipocyte (no differences from controls, p = 0.1790). Green spectrum are TUNEL positive objects, blue spectrum are nuclei

Relative expression of PhIP-metabolizing enzymes in cell lines

Since low level PhIP exposure did not induce cell death in the adipocyte cell line, but did in Caco-2 and HepG2 cells, we then examined the expression levels of genes involved in the metabolic activation of PhIP, and determined if exposure of the cell lines to PhIP could induce the expression of these genes. Exposure of each cell line to 5 nM PhIP did not significantly increase the expression of PhIP-metabolizing genes in any cell line tested (Fig. 3a). We then compared baseline levels of CYP1A1, CYP1A2, SULT1A1 and UGT1A1 in adipocytes, HepG2 and Caco-2 cell lines. Expression of CYP1A1 was 10- and 100-fold higher in HepG2 and Caco-2 cell lines compared to adipocytes, where expression of CYP1A1 was negligible. CYP1A2 expression was negligible in all cell lines, while SULT1A1 and UGT1A1 were most highly expressed in HepG2 cells (Fig. 3b). Exposure of the cell lines to 5 nM PhIP did not result in any changes in this observed pattern (Fig. 3c).

Fig. 3
figure 3

a Relative gene expression of PhIP metabolizing enzymes compared to controls after exposure to 5 nM PhIP for 72 h. Data is presented as Mean ± SD *p < 0.05. Expression of CYP1A1 (p = 0.0511, p = 0.5459, p = 0.1615), CYP1A2 (p = 0.1194, p = 0.6966, p = 0.4045), SULT1A1 (p = 0.0866, p = 0.4488, p = 0.6952) and UGT1A1 (p = 0.1719, p = 0.5599, p = 0.3261) in Adipocyte, HepG2 and Caco-2 cells respectively. There were no significant differences from controls. b Relative gene expression of baseline PhIP metabolizing enzymes normalized to adipocyte controls. Data is presented as Mean ± SD *p < 0.05. Expression of CYP1A1 (p = 0.0030, p = 0.0055), CYP1A2 (p = 0.0044, p = 0.0047), SULT1A1 (p = 0.0001, p = 0.0358) and UGT1A1 (p = 0.0104, p = 0.0001) in HepG2 and Caco-2 cells respectively. HepG2 and Caco-2 cells had significantly higher baseline PhIP metabolizing enzymes than adipocytes with the exception of CYP1A2. c Relative gene expression of PhIP metabolizing enzymes after 5nM PhIP exposure, normalized to adipocyte treatment. Data is presented as Mean ± SD *p < 0.05. Expression of CYP1A1 (p = 0.0030, p = 0.0037), CYP1A2 (p = 0.0140, p = 0.3050), SULT1A1 (p = 0.0005, p = 0.0516) and UGT1A1 (p = 0.0001, p = 0.0754) in HepG2 and Caco-2 cells respectively. HepG2 cell expression was significantly higher than the adipocyte in all four PhIP metabolizing enzymes while Caco-2 was significantly higher in the CYP1A1 enzyme only

Global gene expression changes in h-MSC-derived adipocytes exposed to PhIP

PhIP was taken up and retained to a greater degree by adipocytes compared to HepG2 and Caco2; thus, we sought to determine whether PhIP within the adipocyte cell line remained inert, or exhibited biological effects. Others have reported interactions of PhIP with both the androgen and estrogen receptor at physiologically achievable concentrations to induce the expression of steroid-responsive genes [2931]. Therefore, we then examined global gene expression in mature adipocytes. We exposed mature adipocytes to varying concentrations of PhIP (5 nM – 2.5 μM) for 72 h. Exposure to 5nM PhIP resulted in significant differential expression of 64 genes in an unadjusted model (Additional file 1: Table S1, see end of document). After adjustment for multiple comparisons, 8 genes (ANGPTL2, CD14, CIDEA, EGR1, FOS, IGFBP5, PALM and PSAT1) remained significant at the p = 0.05 level (Table 2). Of the eight genes, four of them (ANGPTL2, EGR1, FOS and IBP5) persisted to be important in the gene-gene interaction analysis.

Table 2 Differential expression of genes in h-MSC-derived adipocytes exposed to PhIP

Gene-gene interaction and pathway analysis

To functionally validate the roles of these genes in biological pathways, the top significantly induced genes compared to vehicle-treated controls (at p = 0.01) in the unadjusted model were then subjected to gene-gene interaction and pathway analysis. From this analysis, we found that PhIP modulated genes controlled by the STAT3 transcriptional factor (Fig. 4). Among genes controlled by STAT3, exposure to PhIP resulted in the down-regulation of early growth response 1 (EGR1) and prostaglandin-endoperoxide synthase 2 (PTGS2) via c-Fos, while insulin binding protein 5 (IBP5) was up regulated. The expression of other transcription factors (ANGPTL2, HP, LEP, SAA1, and SAA2) were up-regulated by PhIP exposure. When direct interactions were assessed, two genes were found to be down-regulated (i.e. cFos, EGR1; p = 0.61E-05 and 1.25 E-05, respectively), and IBP5 was upregulated (p = 0.014). Analysis also revealed two genes that are part of the gene-gene interaction network, which are predicted to be down regulated (i.e. PTGS2 and SGK1). Through gene-gene interaction analysis we were able to reveal that the down-regulated transcriptional factors cFOS and EGR1 that are downstream of MAPK signaling display direct physical DNA-DNA interaction as predicted by Metacore software (Fig. 5) and elsewhere.

Fig. 4
figure 4

5 nM PhIP activates STAT3-mediated transcriptional regulation. Genes within the STAT3 pathway were modulated by exposure to 5 nM PhIP for 72 h. One star denotes significance at <0.05 in the unadjusted model and two stars denote significance at <0.05 after adjustment for multiple comparisons

Fig. 5
figure 5

Direct interactions between genes induced or suppressed by 5nM PhIP exposure. One star denotes significance at <0.05 in the unadjusted model and two stars denote significance at <0.05 after adjustment for multiple comparisons

Discussion

This is the first report of PhIP mediating changes in gene expression within the adipocyte. Along with other heterocyclic amines, PhIP is lipophilic [22]. We hypothesized that PhIP would be preferentially taken up by the adipocyte and associate with the lipid droplet and not undergo metabolic activation that leads to DNA damage. In agreement with this hypothesis, PhIP in fact is not activated in the adipocyte and does not cause cell death measurable by TUNEL. Additionally, the adipocyte cell line employed in these studies had negligible expression of the enzymes necessary for metabolic activation, and exposure to PhIP did not induce their expression. However, we did find that PhIP induces changes in gene expression in various genes that are related to inflammation (SAA1, LEP), diabetes (IBP5) and cancer risk (SAA2).

Reports in the literature support the role of adipose tissue as sites of interaction with environmental chemicals. It has long been recognized that persistent organic pollutants (POPs) are sequestered in adipose tissue, and this sequestration offers initial protection of other organs and tissues from the harmful effects of these chemicals. Over time, however, adiposity can increase the body burden of POPs; weight loss and normal lipolysis can then release these chemicals, resulting in tissue and organ damage. Adipose is also a target of POPs and is thought to participate in the mediation of the metabolic effects of POPs (reviewed in [32]). In addition to POPs, there are recent reports that Bisphenol-A (BPA), a known endocrine disrupting chemical, can regulate the expression of genes within adipocytes that participate in metabolic homeostasis and obesity [33]. Efforts to limit exposures to these chemicals have received much attention. There have also been efforts to reduce the consumption of red and processed meats based on their association with cancer, cardiovascular disease and diabetes [34]. And while obesity is a risk factor for these diseases as well, the potential effect of food-borne HCAs such as PhIP on adipocyte function has been unrecognized. The observation that PhIP induces gene expression changes in important pathways could have profound implications in human health, particularly in Western countries where the consumption of well-done cooked meat is common. And analysis of world meat consumption patterns documents a surge in meat consumption by emerging countries such as China and Brazil [35].

Analysis of global gene expression revealed that low levels of PhIP perturbed pathways related to obesity, inflammation, diabetes and risk of certain cancers. For example, cFOS positively binds and activates STAT3, which positively activates EGR1, cFOS, and PTGS2 (COX-2). PhIP exposure suppressed this pathway, while it induced the expression of IBP5, a biomarker of obesity [36, 37], suggesting that the effects of PhIP exposure are complex. However, STAT3 can physically bind to and positively activate serum amyloid A1 binding proteins A1 and A2 (SAA1, SAA2) and these proteins were up-regulated by PhIP exposure. SAA1 encodes a member of the serum amyloid A family of apolipoproteins, which is a major acute phase protein that is highly expressed in response to inflammation and tissue injury [38, 39]. SAA1 also plays an important role in high density lipoprotein (HDL) metabolism and cholesterol homeostasis [40]. High levels of this protein are associated with chronic inflammatory diseases [38]. SAA2 may also be a potential biomarker for certain tumors [41, 42]. In the present study, these genes were significantly modulated in the unadjusted models, and some retained significance when adjusted for multiple comparisons. Future studies in our laboratory will explore the impact of PhIP exposure under varied experimental conditions to determine the exact contributions of pathway genes to the effects of PhIP on adipocytes.

The magnitude of changes in gene expression was modest in this report, but a strength of the current study is that the PhIP concentrations examined approach dietary relevance in humans. PhIP daily intake estimates range from 0 to 630 ng/day according to meat-specific food frequency questionnaires (FFQs) and diet diaries coupled with the NCI Computerized Heterocyclic Amines Resource for Research in Epidemiology of Disease (CHARRED), [43, 44]. By examining low doses of PhIP, we more closely approach levels that could be expected in human populations. While we only reported on the effects of 5 nM PhIP, concentrations ranging up to 100 μM were also examined. Interestingly, as PhIP concentrations increased, the numbers of genes and pathways that were significantly changed decreased. This suggests that in vitro studies employing high, non-physiological concentrations may not identify important interactions. We also measured changes that occurred after relatively short time intervals; effects of long-term, low-dose exposures are unclear and are the subject of further investigation in our laboratory.

Nonetheless, the observation that PhIP is mediating gene expression changes within the adipocyte, and that the pathways most affected are related to cancer and other chronic disease is intriguing. The established paradigm of PhIP-induced carcinogenesis involves metabolic activation [45], DNA adduct formation resulting in deleterious mutations that initiate carcinogenesis [28, 4648]. However, there are reports in the literature of PhIP interacting with both the estrogen and androgen receptor at physiologically achievable concentrations [2931], indicating that DNA adduct formation is only one pathway of PhIP effects, and perturbation of the microenvironment produced by adiposity could generate changes conducive to cancer and other chronic diseases.

Conclusions

The observation that PhIP is mediating gene expression changes within the adipocyte, and that the pathways most affected are related to cancer and other chronic disease suggests studies that address mechanisms other than DNA adduct formation are needed to clarify the relationship between dietary carcinogens such as PhIP with cancer, obesity and diabetes.

Abbreviations

ATCC, American Type Culture Collection; Caco-2, colorectal cell line; CCPM, corrected cycles per minute; CHARRED, computerized heterocyclic amines resource for research in epidemiology of disease; CYPs, cytochrome P450s; DEGs, differentially expressed genes; DMSO, dimethyl sulfoxide; EGR1, early growth response 1; FBS, fetal bovine serum; FDM, fat differentiation medium; FFQs, food frequency questionnaires; GM, growth medium; HCAs, heterocyclic aromatic amines; HDL, high density lipoprotein; HepG2, hepatocellular cell line; hMSC, human mesenchymal stem cells; HP, haptoglobin; IBP5, insulin binding protein 5; IGFBP5, insulin-like growth factor binding protein 5; LEP, leptin; Limma, linear models for microarray data; N2-OH-PhIP, N-hydroxy-PhIP; NAT, N-acetyltransferase; NMT, N-methyltransferase; PET, polyethylene terephthalate; PhIP, 2-amino-1-methyl-6-phenylimidazol[4,5-b]pyridine; POP, persistent organic pollutants; PTGS2, prostaglandin-endoperoxide synthase 2; SAA1, SAA2, serum amyloid A1 binding proteins A1 and A2; SULTs, sulfotransferase; UGTs, glucuronidation; αMEM, alpha minimum essential medium

References

  1. Covaci A, de Boer J, Ryan JJ, Voorspoels S, Schepens P. Distribution of organobrominated and organochlorinated contaminants in Belgian human adipose tissue. Environ Res. 2002;88:210–8.

    Article  CAS  PubMed  Google Scholar 

  2. Bluher M. Transgenic animal models for the study of adipose tissue biology. Best Pract Res Clin Endocrinol Metab. 2005;19:605–23.

    Article  PubMed  Google Scholar 

  3. Felton JS, Knize MG, Shen NH, Andresen BD, Bjeldanes LF, Hatch FT. Identification of the mutagens in cooked beef. Environ Health Perspect. 1986;67:17–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Edwards RJ, Murray BP, Murray S, Schulz T, Neubert D, Gant TW, Thorgeirsson SS, Boobis AR, Davies DS. Contribution of CYP1A1 and CYP1A2 to the activation of heterocyclic amines in monkeys and human. Carcinogenesis. 1994;15:829–36.

    Article  CAS  PubMed  Google Scholar 

  5. He X, Feng S. Role of metabolic enzymes P450 (CYP) on activating procarcinogen and their polymorphisms on the risk of cancers. Curr Drug Metab. 2015;16:850–63.

    Article  CAS  PubMed  Google Scholar 

  6. Turesky RJ, Lang NP, Butler MA, Teitel CH, Kadlubar FF. Metabolic activation of carcinogenic heterocyclic aromatic amines by human liver and colon. Carcinogenesis. 1991;12:1839–45.

    Article  CAS  PubMed  Google Scholar 

  7. Dellinger RW, Chen G, Blevins-Primeau AS, Krzeminski J, Amin S, Lazarus P. Glucuronidation of PhIP and N-OH-PhIP by UDP-glucuronosyltransferase 1A10. Carcinogenesis. 2007;28:2412–8.

    Article  CAS  PubMed  Google Scholar 

  8. Malfatti MA, Ubick EA, Felton JS. The impact of glucuronidation on the bioactivation and DNA adduction of the cooked-food carcinogen 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine in vivo. Carcinogenesis. 2005;26:2019–28.

    Article  CAS  PubMed  Google Scholar 

  9. Williams JA, Stone EM, Fakis G, Johnson N, Cordell JA, Meinl W, Glatt H, Sim E, Phillips DH. N-Acetyltransferases, sulfotransferases and heterocyclic amine activation in the breast. Pharmacogenetics. 2001;11:373–88.

    Article  CAS  PubMed  Google Scholar 

  10. Ochiai M, Watanabe M, Kushida H, Wakabayashi K, Sugimura T, Nagao M. DNA adduct formation, cell proliferation and aberrant crypt focus formation induced by PhIP in male and female rat colon with relevance to carcinogenesis. Carcinogenesis. 1996;17:95–8.

    Article  CAS  PubMed  Google Scholar 

  11. Ito N, Hasegawa R, Sano M, Tamano S, Esumi H, Takayama S, Sugimura T. A new colon and mammary carcinogen in cooked food, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP). Carcinogenesis. 1991;12:1503–6.

    Article  CAS  PubMed  Google Scholar 

  12. Hasegawa R, Sano M, Tamano S, Imaida K, Shirai T, Nagao M, Sugimura T, Ito N. Dose-dependence of 2-amino-1-methyl-6-phenylimidazo[4,5-b]-pyridine (PhIP) carcinogenicity in rats. Carcinogenesis. 1993;14:2553–7.

    Article  CAS  PubMed  Google Scholar 

  13. Kashiwabara S, Kashimoto N, Uesaka T, Wakabayashi K, Kamiya K, Watanabe H. Tumor induction by azoxymethane (AOM) and 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) in F344 rat gastric mucosa featuring intestinal metaplasia caused by X-irradiation. J Exp Clin Cancer Res. 2005;24:305–12.

    CAS  PubMed  Google Scholar 

  14. Li G, Wang H, Liu AB, Cheung C, Reuhl KR, Bosland MC, Yang CS. Dietary carcinogen 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine-induced prostate carcinogenesis in CYP1A-humanized mice. Cancer Prev Res (Phila). 2012;5:963–72.

    Article  CAS  Google Scholar 

  15. Qiu C, Shan L, Yu M, Snyderwine EG. Steroid hormone receptor expression and proliferation in rat mammary gland carcinomas induced by 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine. Carcinogenesis. 2005;26:763–9.

    Article  CAS  PubMed  Google Scholar 

  16. Shan L, Yu M, Snyderwine EG. Global gene expression profiling of chemically induced rat mammary gland carcinomas and adenomas. Toxicol Pathol. 2005;33:768–75.

    Article  CAS  PubMed  Google Scholar 

  17. Pan A, Sun Q, Bernstein AM, Schulze MB, Manson JAE, Stampfer MJ, Willett WC, Hu FB. Red meat consumption and mortality: results from two prospective cohort studies. Arch Intern Med. 2012;172:555–63.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Keating GA, Bogen KT. Estimates of heterocyclic amine intake in the US population. J Chromatogr B Analyt Technol Biomed Life Sci. 2004;802:127–33.

    Article  CAS  PubMed  Google Scholar 

  19. Obesity: preventing and managing the global epidemic. Report of a WHO consultation. World Health Organ Tech Rep Ser 2000, 894:i-xii, 1-253.

  20. Calle EE, Rodriguez C, Walker-Thurmond K, Thun MJ. Overweight, obesity, and mortality from cancer in a prospectively studied cohort of U.S. adults. N Engl J Med. 2003;348:1625–38.

    Article  PubMed  Google Scholar 

  21. Mullerova D, Kopecky J. White adipose tissue: storage and effector site for environmental pollutants. Physiol Res. 2007;56:375–81.

    CAS  PubMed  Google Scholar 

  22. Turteltaub KW, Vogel JS, Frantz CE, Shen N. Fate and distribution of 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine in mice at a human dietary equivalent dose. Cancer Res. 1992;52:4682–7.

    CAS  PubMed  Google Scholar 

  23. Martin FL, Carmichael PL, Crofton-Sleigh C, Venitt S, Phillips DH, Grover PL. Genotoxicity of human mammary lipid. Cancer Res. 1996;56:5342–6.

    CAS  PubMed  Google Scholar 

  24. Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245–50.

    Article  CAS  PubMed  Google Scholar 

  25. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

  26. Pfau W, Martin FL, Cole KJ, Venitt S, Phillips DH, Grover PL, Marquardt H. Heterocyclic aromatic amines induce DNA strand breaks and cell transformation. Carcinogenesis. 1999;20:545–51.

    Article  CAS  PubMed  Google Scholar 

  27. Knasmuller S, Schwab CE, Land SJ, Wang CY, Sanyal R, Kundi M, Parzefall W, Darroudi F. Genotoxic effects of heterocyclic aromatic amines in human derived hepatoma (HepG2) cells. Mutagenesis. 1999;14:533–40.

    Article  CAS  PubMed  Google Scholar 

  28. Pezdirc M, Zegura B, Filipic M. Genotoxicity and induction of DNA damage responsive genes by food-borne heterocyclic aromatic amines in human hepatoma HepG2 cells. Food Chem Toxicol. 2013;59:386–94.

    Article  CAS  PubMed  Google Scholar 

  29. Lauber SN, Ali S, Gooderham NJ. The cooked food derived carcinogen 2-amino-1-methyl-6-phenylimidazo[4,5-b] pyridine is a potent oestrogen: a mechanistic basis for its tissue-specific carcinogenicity. Carcinogenesis. 2004;25:2509–17.

    Article  CAS  PubMed  Google Scholar 

  30. Papaioannou MD, Koufaris C, Gooderham NJ. The cooked meat-derived mammary carcinogen 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) elicits estrogenic-like microRNA responses in breast cancer cells. Toxicol Lett. 2014;229:9–16.

    Article  CAS  PubMed  Google Scholar 

  31. Glass-Holmes M, Aguilar BJ, Gragg 3rd RD, Darling-Reed S, Goodman CB. Characterization of 2-amino-1-methyl-6-phenylimidazo[4,5b]pyridine at androgen receptor: mechanistic support for its role in prostate cancer. Am J Cancer Res. 2015;5:191–200.

    CAS  PubMed  Google Scholar 

  32. La Merrill M, Emond C, Kim MJ, Antignac JP, Le Bizec B, Clement K, Birnbaum LS, Barouki R. Toxicological function of adipose tissue: focus on persistent organic pollutants. Environ Health Perspect. 2013;121:162–9.

    PubMed  Google Scholar 

  33. Menale C, Piccolo MT, Cirillo G, Calogero RA, Papparella A, Mita L, Del Giudice EM, Diano N, Crispi S, Mita DG. Bisphenol A effects on gene expression in adipocytes from children: association with metabolic disorders. J Mol Endocrinol. 2015;54:289–303.

    Article  CAS  PubMed  Google Scholar 

  34. Battaglia Richi E, Baumer B, Conrad B, Darioli R, Schmid A, Keller U. Health risks associated with meat consumption: a review of epidemiological studies. Int J Vitam Nutr Res. 2015;85:70–8.

    Article  PubMed  Google Scholar 

  35. Sans P, Combris P. World meat consumption patterns: an overview of the last fifty years (1961-2011). Meat Sci. 2015;109:106–11.

    Article  CAS  PubMed  Google Scholar 

  36. Leyh B, Dittmer A, Lange T, Martens JW, Dittmer J. Stromal cells promote anti-estrogen resistance of breast cancer cells through an insulin-like growth factor binding protein 5 (IGFBP5)/B-cell leukemia/lymphoma 3 (Bcl-3) axis. Oncotarget. 2015;6:39307–28.

    PubMed  PubMed Central  Google Scholar 

  37. Kallio P, Tolppanen AM, Kolehmainen M, Poutanen K, Lindstrom J, Tuomilehto J, Kuulasmaa T, Kuusisto J, Pulkkinen L, Uusitupa M. Association of sequence variations in the gene encoding insulin-like growth factor binding protein 5 with adiponectin. Int J Obes (Lond). 2009;33:80–8.

    Article  CAS  Google Scholar 

  38. Rossmann C, Windpassinger C, Brunner D, Kovacevic A, Schweighofer N, Malli R, Schuligoi R, Prokesch A, Kluve-Beckerman B, Graier WF, et al. Characterization of rat serum amyloid A4 (SAA4): a novel member of the SAA superfamily. Biochem Biophys Res Commun. 2014;450:1643–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Anderberg RJ, Meek RL, Hudkins KL, Cooney SK, Alpers CE, Leboeuf RC, Tuttle KR. Serum amyloid A and inflammation in diabetic kidney disease and podocytes. Lab Invest. 2015;95:250–62.

    Article  CAS  PubMed  Google Scholar 

  40. Poulsen SS, Saber AT, Mortensen A, Szarek J, Wu D, Williams A, Andersen O, Jacobsen NR, Yauk CL, Wallin H, et al. Changes in cholesterol homeostasis and acute phase response link pulmonary exposure to multi-walled carbon nanotubes to risk of cardiovascular disease. Toxicol Appl Pharmacol. 2015;283:210–22.

    Article  CAS  PubMed  Google Scholar 

  41. Biran H, Friedman N, Neumann L, Pras M, Shainkin-Kestenbaum R. Serum amyloid A (SAA) variations in patients with cancer: correlation with disease activity, stage, primary site, and prognosis. J Clin Pathol. 1986;39:794–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Sung HJ, Ahn JM, Yoon YH, Rhim TY, Park CS, Park JY, Lee SY, Kim JW, Cho JY. Identification and validation of SAA as a potential lung cancer biomarker and its involvement in metastatic pathogenesis of lung cancer. J Proteome Res. 2011;10:1383–95.

    Article  CAS  PubMed  Google Scholar 

  43. Deziel NC, Buckley TJ, Sinha R, Abubaker S, Platz EA, Strickland PT. Comparability and repeatability of methods for estimating the dietary intake of the heterocyclic amine contaminant 2-amino-1-methyl-6-phenylimidazo[4,5b]pyridine (PhIP). Food Addit Contam Part A Chem Anal Control Expo Risk Assess. 2012;29:1202–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Sinha R. An epidemiologic approach to studying heterocyclic amines. Mutat Res. 2002;506–507:197–204.

    Article  PubMed  Google Scholar 

  45. Zhao K, Murray S, Davies DS, Boobis AR, Gooderham NJ. Metabolism of the food derived mutagen and carcinogen 2-amino-1-methyl-6-phenylimidazo(4,5-b)pyridine (PhIP) by human liver microsomes. Carcinogenesis. 1994;15:1285–8.

    Article  CAS  PubMed  Google Scholar 

  46. Muckel E, Frandsen H, Glatt HR. Heterologous expression of human N-acetyltransferases 1 and 2 and sulfotransferase 1A1 in Salmonella typhimurium for mutagenicity testing of heterocyclic amines. Food Chem Toxicol. 2002;40:1063–8.

    Article  CAS  PubMed  Google Scholar 

  47. Hein DW, Rustan TD, Ferguson RJ, Doll MA, Gray K. Metabolic activation of aromatic and heterocyclic N-hydroxyarylamines by wild-type and mutant recombinant human NAT1 and NAT2 acetyltransferases. Arch Toxicol. 1994;68:129–33.

    Article  CAS  PubMed  Google Scholar 

  48. Buonarati MH, Turteltaub KW, Shen NH, Felton JS. Role of sulfation and acetylation in the activation of 2-hydroxyamino-1-methyl-6-phenylimidazo[4,5-b]pyridine to intermediates which bind DNA. Mutat Res. 1990;245:185–90.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We want to thank Alena Savenka (University of Arkansas for Medical Sciences) for TUNEL analysis.

Funding

This work was supported by the Applied Biosciences Institute at the University of Arkansas for Medical Sciences, Little Rock, AR.

Availability of data and materials

Not applicable.

Authors’ contributions

LJR performed cell culture, PhIP treatment, RT-PCR experiments, analysis, data interpretation and wrote the manuscript. AGB, MSO, VR and SAK contributed to analysis, interpretation of data and discussion of the results. BN contributed hMSC protocols and discussion of the results. AYB and SAK formulated the experimental design. All authors critically revised the article and approved the final version.

Authors’ information

None.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Susan Kadlubar.

Additional file

Additional file 1: Table S1.

Global gene expression changes in h-MSC-derived adipocytes exposed to PhIP. (DOCX 24 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rogers, L.J., Basnakian, A.G., Orloff, M.S. et al. 2-amino-1-methyl-6-phenylimidazo(4,5-b) pyridine (PhIP) induces gene expression changes in JAK/STAT and MAPK pathways related to inflammation, diabetes and cancer. Nutr Metab (Lond) 13, 54 (2016). https://doi.org/10.1186/s12986-016-0111-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12986-016-0111-0

Keywords