Sean Timothy Francis O’Connor1 Jiaqi Lan2 Matthew North3 Alexandre Loguinov3 Luoping Zhang4 Martyn T. Smith4 April Z. Gu2 Chris Vulpe3*- 1Vulpe Laboratory, Graduate Group in Environmental Health Sciences, University of California, Berkeley, CA, USA
- 2Gu Laboratory, Department of Civil and Environmental Engineering, Northeastern University, Boston, MA, USA
- 3Vulpe Laboratory, Department of Nutritional Science and Toxicology, University of California, Berkeley, CA, USA
- 4Genes and Environmental Laboratory, Division of Environmental Health Sciences, School of Public Health, University of California, Berkeley, CA, USA
Benzo[a]pyrene (BaP) is a ubiquitous, potent, and complete carcinogen resulting from incomplete organic combustion. BaP can form DNA adducts but other mechanisms may play a role in toxicity. We used a functional toxicology approach in S. cerevisiae to assess the genetic requirements for cellular resistance to BaP. In addition, we examined translational activities of key genes involved in various stress response pathways. We identified multiple genes and processes involved in modulating BaP toxicity in yeast which support DNA damage as a primary mechanism of toxicity, but also identify other potential toxicity pathways. Gene ontology enrichment analysis indicated that DNA damage and repair as well as redox homeostasis and oxidative stress are key processes in cellular response to BaP suggesting a similar mode of action of BaP in yeast and mammals. Interestingly, toxicant export is also implicated as a potential novel modulator of cellular susceptibility. In particular, we identified several transporters with human orthologs (solute carrier family 22) which may play a role in mammalian systems.
Introduction
Benzo[a]pyrene (BaP) is a potent, ubiquitous, and persistent animal carcinogen as well as a probable human carcinogen (International Agency for Research on Cancer class I carcinogen) produced by incomplete organic matter combustion. It readily distributes to soil and dust particles (Hattemer-Frey and Travis, 1991; Chung et al., 2011). The major routes of exposure are food, both due to endogenous contamination and cooking chemistry, as well as ambient air (indoor and outdoor) where vehicular exhaust and traditional combustion cooking and heating processes are the main source of BaP exposure aside from smoking (Mitchell, 1982; Wester et al., 1990; Laurent et al., 2002; Roseiro et al., 2011). Inhalation is the most effective route of BaP absorption, followed by ingestion, and then dermal uptake; however, the largest exposure is likely by ingestion (Wester et al., 1990; Hattemer-Frey and Travis, 1991; Edoardo, 1992; Laurent et al., 2002; Fromme et al., 2004; Moody et al., 2007).
Once in the body, BaP preferentially distributes to the kidneys and liver, with lower levels partitioning to the testes, brain, and spleen (Mitchell, 1982). BaP is difficult to excrete owing in part to its affinity for adipose tissues and its interaction with serum proteins, both of which can cause the chemical to remain in the body for days despite metabolic action (Shu and Bymun, 1983). BaP is metabolized by cytochrome P-450 enzymes (CYPs) to more excretable yet highly reactive and toxic metabolites (Baird et al., 2005). One of these, benzo[a]pyren-7,8-dihydrodiol-9,10-epoxide (BPDE) is thought to form DNA adducts and facilitate the induction of cancer via guanine to thymine transversion mutations, even though cancer promotion pathways remain unclear (Kapitulnik et al., 1978; Bartosiewicz et al., 2001; Baird et al., 2005; Hockley et al., 2009). BaP can also form hydroquinones that undergo a redox cycle generating genotoxic reactive oxygen species such as superoxide and hydrogen peroxide (Kim and Lee, 1997; Baird et al., 2005). Other potential mechanisms of toxicity and cancer promotion involve oxidative damage, immune suppression, activation of transposons, estrogen receptor activation, and disruption of cell-signaling (Romero et al., 1997; Vondrácek et al., 2002; Baird et al., 2005; Stribinskis and Ramos, 2006). Currently, the temporal ordering, specificity, and relative importance of these mechanisms in BaP toxicity is not fully elucidated.
In yeast, including S. cerevisiae, several cytochrome P-448 monooxygenase enzymes can perform Phase I metabolism on compounds like BaP and activate them to their toxic form (i.e., BPDE) in a manner analogous to mammalian microsomes (King et al., 1983; Käppeli, 1986). These yeast microsomes are indeed similar in weight and function to those observed in other mammalian species, including humans, but are enzymatically less efficient and the microsomes more rapidly lose activity when isolated in vitro (Azari and Wiseman, 1981; Käppeli, 1986). Moreover, the genetic expression of this rudimentary metabolism is not regulated by chemical exposure, but rather by a complex and unclear interaction between the presence of oxygen and the type of metabolism the yeast is undergoing at the time of growth (King et al., 1983). While it is clear that S. cerevisiae have the enzymatic capability to activate BaP into the toxic metabolites observed in other organisms including humans, the extent of activation has not been examined in vivo and the value of using a microsomal supernatant S-9 fraction solution as a metabolic adjuvant for bioactivation has not been assessed (Käppeli, 1986).
We used a functional genomic approach in yeast to identify the genes important for resistance or sensitivity to BaP exposure. Yeast is an useful model for functional assessment of the genetic requirements in toxicology. The yeast genome is fully characterized, primary biological functions are conserved and the relatively rapid reproductive rate of yeast enables exposures over multiple generations. In addition, deletion strains were generated through systematic insertion of “molecular barcodes” that individually identifies each gene knockout. As a result, functional genomic approaches in which entire sets of deletion strains can be pooled and simultaneously assayed for growth effects are possible (Giaever et al., 2002). In contrast to correlative genomic approaches such as transcriptomics, functional genomics directly identifies the genes tied to a phenotypic outcome (such as growth or toxicity). We have previously used this approach with several toxicants and shown the capability to translate results from yeast to mammalian systems (Jo et al., 2009a,b; Zhang et al., 2010; Ren et al., 2011). We also performed temporal protein expression profiling of yeast stress response. We used a selected set of yeast strains expressing full-length, chromosomally tagged green fluorescent protein fusion (GFP) proteins with genes encoding general stress response, oxidative stress response, chemical stress response, protein stress response, and DNA stress response (Huh et al., 2003). This approach allows measurement of real-time protein expression changes in response to stress. The selected stress response genes are present and highly conserved in most cell types of metazoans and are activated at significantly lower toxicant concentrations than those causing overt cellular injury (Kultz, 2005; Simmons et al., 2009).
Materials and Methods
Functional Genomic Analysis
Yeast strains
Diploid yeast deletion strains based on the BY4743 background (Invitrogen Corporation, Carlsbad, CA, USA) were used both for the parallel analysis pools (n = 4,757) and individual strain growth assays. For both the deletion strain pools and the individual strains, growth assays were performed in liquid media (2% peptone, 2% dextrose, and 1% dextrose by volume).
Yeast growth curves
In order to assay yeast response to BaP, yeast strains were grown to mid-log phase, diluted to an optical density of 0.0165 at 600 nm (OD600), and dispensed to a 48-well plate (non-treated polystyrene, Grenier Bio-One, Monroe, NC, USA). BaP (Sigma-Aldrich, St. Louis, MO, USA) was subsequently added to each well to the desired final concentration. At least two replicates for each dose and three biological replicates for each strain were performed. Each plate was subsequently incubated in a GENios microplate reader (Tecan, Durham, NC, USA) set to 30°C with intermittent shaking and OD595 measurements at 595 nm were recorded every 15 min for 24 h. The raw absorbance data obtained was corrected for background, averaged across replicates, and plotted against time. Excel 2010 (Microsoft Corporation, Redmond, WA, USA) was used to calculate the area under the curve (AUC) as a measure of growth (the outcome of interest) and to perform statistical analyses, namely, comparing AUC values by means of Student’s t-test, assuming constant variance between samples.
Benzo[a]pyrene stock solutions
BaP was solubilized in dimethylsulfoxide (DMSO) and the final concentration of DMSO in the yeast media never exceeded one percent by volume, the standard DMSO concentration at which no observable effect occurred on yeast. Backgrounds and controls included one percent DMSO to account for any potential effect. The BaP stock solutions were stored in the chemical hood shielded from light. As long as it is shielded from light, BaP is stable in DMSO solutions, with less than 20% degradation over 3 months at 25°C (Dabrowska et al., 2008).
IC20 determination
Initially, using the method developed by Pierce et al. (2007), the concentration of BaP which resulted in a 20% decrease in growth in wild type yeast (IC20) was determined. The IC20 point was chosen for its ability to affect yeast deletion strains measurably and differentially without incurring the loss of toxicant specificity that occurs at higher concentrations (Pierce et al., 2007).
Yeast parallel deletion analysis
Targeted inactivation of a gene, “knockout,” is often used to investigate the gene’s function. Similarly, this same technique can be used to identify how a gene or lack thereof affects survival of a model organism in response to a toxicant (Giaever et al., 2002). Using a gene knockout library of S. cerevisiae strains that represent 96% of the yeast genome’s open reading frames (ORF), we can thus determine in parallel the effect of each gene on the growth outcome as it relates to BAP exposure (Giaever et al., 2002). This is accomplished by way of unique DNA “molecular barcodes” that correspond to each knockout strain. After varying time and dose exposures to BAP, the pooled yeast knockout library cells can then be harvested and the extracted DNA barcodes hybridized to oligonucleotide arrays, thereby allowing for a quantitative determination and comparison of growth outcome for each individual strain/gene (Giaever et al., 2002).
To this end, after determining the IC20, pool growth, genomic DNA extraction, and array hybridization were conducted as described by Pierce et al. (2007), with a few minor changes described herein. In summary, viable homozygous diploid deletion mutant strains of yeast (n = 4,757) were continuously subjected to three different BaP concentrations (IC20, 50% of the IC20, and 25% of the IC20) over a growth period of 5 and 15 generations (5G and 15G). In this way, three total biological replicates were obtained for each dose and generation time period. After collecting the cells, the genomic DNA was extracted with the YDER kit (Pierce Biotechnology, Rockford, IL, USA). The unique barcodes in the DNA, specific to each yeast strain, were amplified by PCR using a set of biotinylated primers, and the resulting products were then hybridized to TAG4 arrays (Affymetrix, Santa Clara, CA, USA).
Differential strain sensitivity analysis
In order to investigate genes required for optimal growth of yeast in the presence of BaP, differential strain sensitivity analysis (DSSA) was used to identify differential growth of viable yeast knockout strains relative to wild type yeast in the presence of the toxicant. To this end, the data obtained from the TAG4 array was log2 transformed, adjusted for signal saturation as described in Pierce et al. (2007) as well as adjusted for mean chip background by means of robust location and scale estimators for log2 transformed amplitudes of null features (18,000 total, evenly spread on the array). To control for variability in strain growth, results from each treatment array were matched to data from 12 controls (5G or 15G) for analysis. Treatment-control pairs were normalized with locally weighted scatterplot smoothing (global normalization), and the difference in growth between strains was identified using an alpha-outlier approach (Loguinov et al., 2004). Data from three biological replicates were combined, resulting in 36 treatment-control data pairs per treatment group. Residual variances (with a robust/* scale estimator) of log2 (treatment/control) for each 36 pairs were inspected using box plots. The “effective pairs” were then determined by excluding pairs with abnormally high residual variance, or with suspected serial correlation in variability (regular patterns in the box plots).
Significant yeast deletion strains (i.e., genes) were statistically inferred by an exact binomial test, with the assumption that the outcomes for each gene in all treatment-control pairs were independent binary variables with the same probability of success for all Bernoulli trials. For a particular gene n, outcomes were considered as “successful” if they were significant in the outlier analysis with q-values ≤0.05 in each of all effective pairs with log2 ratios of the same sign, simultaneously (Loguinov et al., 2004). The corresponding raw p values based on the exact binomial test were then corrected for multiplicity of comparisons using a q-value approach and only the genes with q-value ≤0.05 were considered for further analysis. This approach does not apply a scale estimator and, as a result, it did not require between-chip pair normalization for the statistical inference. Subsequently, a fitness score (ave[log2{Y|X = exposed}] – ave[log2{Y|X = control}]), the difference of the means of the (log2) hybridization signal between exposed and control was determined for each yeast deletion strain. Negative scores indicate sensitivity of the strain to BaP exposure; conversely, a positive score indicates resistance to BaP exposure. Consequently, in the case of a negative score, the gene product deleted from that strain is likely required for resistance to the toxicant. Similarly, in the case of a positive score, the absent gene product may serve to sensitize the cell to BaP (North et al., 2011).
Gene ontology
From the significant genes responsible for resistance or sensitivity to BaP in yeast identified in the DSSA method described above, overrepresented Gene Ontology1 (GO) and MIPS (Munich Information Center for Protein Sequences2) categories were identified to determine which specific (if any) gene clusters related to biological characteristics were relevant to S. cerevisiae’s response to BaP exposure (Ashburner et al., 2000; Mewes et al., 2002). To this end, the Functional Specification resource, FunSpec3, was used in an enrichment analysis to determine the significantly represented GO and MIPS biological categories among the sensitive/resistant strains (Robinson et al., 2002). We used the hypergeometric distribution of FunSpec with a p-value significance threshold of 0.01.
Selected confirmation of yeast strain significance
In order to confirm the data previously obtained, and based on the results from the DSSA and GO, select relevant significant yeast gene deletion strains were chosen for individual testing. Yeast growth curves, as described above, were performed at two different BaP concentrations (100% IC20, 25% IC20) for the chosen deletion strains.
S-9 metabolic adjuvant
S-9 (Molecular Toxicology, Inc., Boone, NC, USA) was obtained as 30% by volume S-9 rat liver microsome mix (NADP, D-glucose-6-phosphate, MgCl2, and KCl co-factors; pH = 7.4 sodium phosphate buffer) containing all co-factors for a sustained enzymatic in vitro metabolic activation. The S-9 was reconstituted with nuclease-free water and stored at −20°C. The S-9 was prepared and mixed with the liquid YPD media immediately prior to dispensing into each of the 48-well plates. The final concentration in each well was either 1 or 2% by volume, based on a prior study which suggested 1.85% by volume as the optimal S-9 concentration for the in vitro activation of BaP with S-9 (Hakura et al., 2001). The yeast growth curve protocol was followed as described previously and the results controlled for background absorption and/or toxicity due to S-9. Our statistical analysis of the results included the use of two-way ANOVA to assess BAP by S-9 interaction as well as the relevant Tukey post hoc tests.
Proteomic Analysis of Stress Response Pathways
GFP-fused yeast strains and stress response ensemble
We employed a library of in frame GFP fusion proteins (Invitrogen, no. 95702) of S. cerevisiae (ATCC 201388), constructed by oligonucleotide-directed homologous recombination to tag each ORF with Aequrea victoria GFP (S65T) in its chromosomal location at the 3′ end. A set of stress response pathway genes, with a total of selected 123 ORFs in yeast, was selected to cover five different functional categories of stress response (Table 1).
Table 1. Cellular stress response pathways ensemble of ORFs for toxicity response in GFP-tagged yeast.
BaP solution
BaP was solubilized in dimethylsulfoxide (DMSO) as described in Section “Benzo[a]Pyrene Stock Solutions,” with rat liver extraction S-9 (Invitrogen, Grand Island) added to obtain a final concentration of 1.4% for bioactivation before exposure as previously discussed. Sub-cytotoxic doses from 0.01 to 1000 mg/L (below IC5 = 3198.89 mg/L, with >95% survival percentage for 24 h exposure by AUC method in Section “Yeast Growth Curves”) were applied.
Real-time transcriptional analysis
A total of six dose concentrations across 6-logs below IC5 were tested using the GFP-tagged library and shown in a heat map. To allow for comparison, we focused our discussion mainly with two concentrations, 10 and 100 mg/L, since these two concentrations are close to those used (83.26 mg/L, 20.82 mg/L) in the deletion library investigation. This study aims to reveal the toxicity mechanism of BAP, so dose concentrations below IC5 but sufficiently high enough to lead to observable molecular disturbance were chosen. We note that there is bio-accumulation of BAP in organisms of four orders of magnitude on average, so cells are likely exposed to much higher concentrations than those typically found in the environment (Janikowska and Wardas, 2002; Wang and Wang, 2006).
GFP tagged yeast strains selected were grown in clear bottom black 384-well plates (Costar) with SD medium at 30°C until the cultures reached early exponential growth (OD600 about 0.2–0.4). Ten microliters of BaP-S-9 mixture was added per well to obtain the final concentration from 1000 to 0.01 mg/L with 10× dilution. PBS-S-9 mixture of equal volume was used as vehicle control. Plates were put in a Micro plate Reader (Synergy H1 Hybrid Multi-Mode, Biotek, Winooski, VT, USA) to read continuously for absorbance (OD600 for cell growth) and GFP signal (filters 485 nm, 528 nm for protein expression) every 5 min for 2 h (with orbital shake at high speed for 1 min before each reading). All tests were performed in triplicate.
Data analysis
OD and GFP raw data from measurement were corrected by background OD and GFP signal of medium control (with or without chemicals). The protein expression for each measurement is then normalized by cell number (ODcorrected) as P = (GFPcorrected/ODcorrected). The P level was corrected with internal control for plate normalization (details not included here). Then the alteration in expression for a given protein at each time point due to chemical exposure, also referred as induction factor I, was represented as I = Pexperiment/Pvehicle by corrected values.
To quantify the chemical-induced protein expression level changes of a treatment, the real-time protein expression profile was then integrated into the Protein Expression Level Index (PELI) as following:
For protein i, the accumulative protein expression over their 2 h exposure period, which is directly related to its function, is calculated as
An assay noise cut-off is determined using a signal-to-noise ratio of 1.5. For triplicates PELIORF were evaluated by Mean ± SD.
Gene ontology
Gene ontology analysis was performed with novel Network Ontology Analysis method (NOA) for enrichment analysis to determine the significantly represented GO biological categories, and to analyze functions of gene network, as it allows enrichment analysis with user defined reference set (Wang et al., 2011). In this study the whole stress library was used as the reference set, and activated ORFs were used as test set, which is defined as PELIORF > 1.5 based on signal-to-noise ratio of the test. A p-value significance threshold of 0.01 was used similar to the deletion library.
Results
BaP IC20 Determination
Based on a literature search of available toxicity data and empirical determination, we used the yeast growth curve assay to analyze the growth response of wild type BY4743 yeast in liquid YPD media as a function of BaP concentration by measuring OD595 every 15 min for 24 h. Total cell growth at each chemical dose was then determined by calculating the AUC of each OD versus time growth curve plot (Figure A1 in Appendix; Figure 1). Visual inspection of Figure A1 in Appendix reveals that most of the BaP toxicity dose-response is expressed later during the stationary growth phase, following exponential growth rather than from initial exposure. Moreover, from this assay, the BaP IC20, or exposure concentration resulting in a 20% inhibition of growth relative to the absence of chemical treatment, was determined to be 330 μM (Figure A2 in Appendix). From the IC20 value thus calculated, we defined two other concentrations to be used subsequently in DSSA as exposure concentrations for the pooled homozygous diploid yeast deletion mutants exposure periods of five and 15 yeast generations (5G and 15G; Table 2).
Figure 1. Average total cell growth. Average total cell growth over the 24 h period was obtained by calculating the area under the curve (AUC) for each dose and averaging the result across three biological replicates. Bars are standard errors of three biological replicates; significance by t-test comparison between each two doses sequentially; *represents p < 0.05.
Table 2. Concentrations used in DSSA pooled homozygous diploid yeast deletion mutant exposures to BaP.
Differential Strain Sensitivity Analysis
Differential strain sensitivity analysis, described in the methods, was performed at the concentrations specified in Table 1 for 5G and 15G exposure periods to identify genes, each represented by a yeast deletion strain, required for resistance or sensitivity to BaP exposure. The results from this assay indicate that the number of yeast genes conferring resistance or sensitivity did not vary by dose for a 5G or 15G exposure, which may be indicative of a sharp threshold response (Figure 2). We note approximately 17-fold more strains exhibiting BaP sensitivity relative to BaP resistance in the 5G exposure (Figure 2) while in the 15G exposure there are more resistant strains than sensitive ones (Figure 2). Overall, more strains responded to BaP exposure during a 5G than during a 15G exposure (Figure 2).
Figure 2. Number of sensitive and resistant yeast deletion strains for each exposure at three BaP concentrations and for five (5G) and 15 (15G) generations.
We found many strains commonly sensitive or resistant in each dose exposure group (5G or 15G) but there was a much smaller overlap between the 5G and 15G exposures (Figure 3).
Figure 3. Venn diagrams of the amount and overlap of sensitive/resistant strains observed across overall exposure duration categories (A) and for each dose group (B).
Overall, 310 yeast mutant deletion strains were resistant and 753 were sensitive to BaP under at least one of the six treatment conditions (three doses, two time periods). Of these sensitive and resistant strains, 209 were resistant and 528 were sensitive under at least three treatment conditions (Tables S1–S4 in Supplementary Material). The treatment results for the most significant deletion strains are reported in Tables 3 and 4.
Table 3. Table of the 25 most sensitive yeast strains exhibiting significantly negative log2 values during a 5G exposure treatment.
Table 4. Table of the 25 most sensitive yeast strains exhibiting significantly negative log2 values during a 15G exposure treatment.
Gene Ontology
From the strains identified in Tables S1 and S2 in Supplementary Material, an enrichment analysis was performed to identify significantly overrepresented categories of genes among the sensitive/resistant strains. The overrepresented categories are summarized in Table S5 in Supplementary Material. Based on the results of Table S5 in Supplementary Material, the significance of the categories represented, and a review of BaP toxicity, we identified select strains associated with biological categories of potential relevance to BaP toxicity. The strains of interest are reported in Table 5, a subset of Table S5 in Supplementary Material. Genes related to transmembrane transport, mitigating reactive oxygen species, and DNA repair were identified as of particular relevance to BaP toxicity.
Table 5. Overrepresented (p-value < 0.01) biological categories based on the MIPS and GO databases of the Functional Specification Instrument for 5G and 15G exposures of select significantly resistant and sensitive strains, as well as any observed overlap strains between 5G and 15G (Ashburner et al., 2000; Mewes et al., 2002; Robinson et al., 2002).
Validation of DSSA Data
From the DSSA results previously described, 18 yeast deletion strains that exhibited significant resistance/sensitivity to BaP exposure were chosen to validate results using individual assays. These 18 deletion strains were independently exposed to BaP concentrations of 0% IC20 (0 μM), 25% IC20 (82.5 μM), and 100% IC20 (330 μM) for 24 h, following the same protocol as the yeast growth curve assay used to identify the IC20. Subsequently, growth inhibition at the various concentrations was calculated in relation to the AUC of each growth curve, in an analogous manner as the IC20 determination. Results and statistical significance of the validation analysis are included in Figure A3 in Appendix. All of the strains tested demonstrated statistically significant inhibition using at least one BaP concentration relative to the 0% IC20 control, with the exception of ste20Δ, ctt1Δ, and oye2Δ whose growth inhibition was not statistically significant. Moreover, of the fifteen strains that demonstrated statistically significant inhibition, only glr1Δ, prx1Δ, and yrm1Δ did not show a dose-dependent growth inhibition response with inhibition increasing concomitantly with dose. The strains whose sensitivity to BaP was confirmed in this assay are the same as the bolded strain subset listed in Table 5.
S-9 Metabolic Activation
In order to assess if BaP activation via S-9 influences toxicity in yeast, we performed wild type yeast growth curve assay at BaP concentrations of 0% IC20 (0 μM), 25% IC20 (82.5 μM), 50% IC20 (165 μM), and 100% IC20 (330 μM) with S-9 extract concentrations of 0, 1, and 2% by volume (Figure 4). Only the 165 μM BaP concentration showed statistically significant more growth inhibition with increasing S-9 extract concentrations (Figure 4). At 82.5 μM, we noted a very modest but significant decrease between the 1 and 2% S-9 but not between 0 and 1% nor 0 and 2%. No difference was seen at the highest BaP dose. Despite the ambiguity of the effect of S-9 on BaP toxicity, there was demonstrable statistical evidence for interaction between BAP and S-9 (p-value = 0.000323). In addition, while BAP toxicity is influenced by S-9 in a non-additive way, the significant effect on growth is sensitive to highly specific BAP and S-9 dose combinations, with maximal consequence centered at 165 μM BAP and 2% S-9 and resulting in an approximately 25% growth reduction relative to the control.
Figure 4. Wild type yeast exposure to various combinations of BaP and S-9 concentrations. *Indicates significance at 0.01 < p < 0.05, **indicates significance at p < 0.01. Bars are standard errors of three biological replicates.
Real-Time Protein Expression in GFP-Infused Stress Pathway Ensemble Yeast Library
We examined real-time differential protein expression profiles (in comparison to control) of BaP-induced stress response over 2 h exposure times for six different concentrations (Figure 5).
Figure 5. Real-time protein expression profile of BaP exposure based on natural log of induction factor. ln I as average of triplicate, Red spectrum colors indicate upregulation, green spectrum colors indicate downregulation in yeast. X-axis bottom: testing time, every 5 min for 2 h. Y-axis left: clusters of proteins by stress response pathways, Y-axis right: ORFs tested.
The profiles are consistent with a dose-dependent cellular response, with varying protein activation profiles and dynamic toxic responses at different dose concentrations. To further quantify and reveal the dominant stress response categories involved in exposure to different chemical concentrations, an expression level rank-based protein expression enrichment analysis was conducted and the results are shown in Figure 6. The proteins and associated stress categories that had higher expression change levels (ranking) were more clearly visualized, indicating that the cellular response shifted with concentration change, from DNA stress at low levels to chemical and protein stress at higher levels, and general stress response for the highest level. At relatively lower dose concentration (<0.1 mg/L), DNA damage repair pathways, primarily base excision repair, were activated. As the dose increases, stress response activation shifts to transport and membrane function, including a wide-ranging activation of the ABC and MFS transporters involved in chemical export. In addition, the protein stress response, including a series of heat shock proteins involved in recognition and repair of damaged/misfolded proteins, is activated. Finally, there is evidence of oxidative damage as indicated by activation of oxidant sensors and enzymes (SLN1, TRX2, and GSH1). At higher concentration (>100 mg/L), stress responses shifted to more general stress and were down regulated, while proteins for apoptosis increased, suggesting that cells transitioned from stress response to cell death.
Figure 6. Enrichment analysis of protein expression according to various stress groups for different concentrations of BaP exposure. For each concentration, the ORFs are ranked in decreasing order by their PELI score (PELIORF, average of triplicates). Different colors indicate different stress groups (see Table 1).
According to the signal-to-noise ratio defined, 24 proteins were activated by a 10 mg/L BaP exposure and 14 proteins for 100 mg/L of BaP. An enrichment analysis was performed based on GO categories to identify significantly overrepresented categories among these activated ORFs, with the whole stress library selected as the reference in Table 6. The results suggested that the main biological process induced by a 10 mg/L BaP exposure focused on transmembrane transport and protein folding. Comparison of the activated proteins in the stress response GFP-tagged library with sensitive genes in the deletion library showed that there were few common proteins, possibly related to the difference in testing period (1G versus 5G) and the approach applied (protein expression versus gene deletion). Also, the GFP-tagged library only focused on stress response with limited ORFs tested (123 ORFs in the GFP-tagged library versus 4,757 ORFs in the gene deletion library), so there may a significant amount of sensitive genes not tested in the GFP-tagged library.
Table 6. Overrepresented (p-value < 0.01) biological categories based on GO databases for a two hour exposure at 10 mg/L and 100 mg/L BaP in yeast (test set) with PELIORF > 1.5, with the stress library serving as a reference set.
Discussion
Although, BaP is a potent carcinogen, our understanding of its toxicity pathways is limited. Indeed, the toxicant can form DNA adducts and reactive oxygen species metabolites that are thought to be the main mechanisms of toxicity. However, these pathways do not fully account for the observed damage and it is likely that other pathways are concurrently at play. Moreover, since DNA adducts constitute only part of the parent BaP toxicity, the paucity of reliable biomarkers specific to the toxicant also motivates the identification of genes responsive to exposure. Although BaP toxicity was less in yeast than for human cells, this was consistent with similar prior studies involving benzene. Resistance to BaP seems to also increase with duration of exposure. Many of the pathways elicited were consistent with prior BaP toxicity studies, which identified reactive oxygen species and DNA damage, but our study also identified solute carriers and transmembrane transport as responsive to BaP exposure and as possible biomarkers. Moreover, addition of S-9 metabolic extract in addition to BaP had limited to no effect on the latter’s toxicity in S. cerevisiae.
BaP Less Toxic to Yeast than Human Cells
BaP was approximately 50–100 times less toxic in S. cerevisiae than has been observed in cultured human cell samples from individuals exposed to BaP (Aust et al., 1980; Vayssier-Taussat et al., 2001). Yeast lacks certain peroxidases found in humans, which are required to transform BaP metabolites, such as hydroquinone, to more toxic forms (North et al., 2011). Moreover, BaP causes human cell proliferation at low levels, largely depending on the proportion of metabolites produced, which could introduce variability in identifying the IC20 (Burchiel et al., 2007). In addition, although the BaP parent compound is relatively stable in the environment, the toxic metabolite(s) is (are) more reactive (Dabrowska et al., 2008). Nevertheless, the BaP IC20 is lower than was found for the three benzene metabolites hydroquinone, catechol, and 1,2,4-Benzentriol independently tested using the same method (North et al., 2011).
We found more sensitive strains in the 5G assay and, conversely more resistant strains in the 15G exposure. Interestingly many of the strains that showed sensitivity in the 5G analysis were resistant in the 15G assays (Figure 7). Cellular responses may increase resistance in these yeast deletion strains over time. Indeed, there are relatively few strains in common between the 5G resistant and the 15G sensitive assays.
Figure 7. Venn diagram overlap between sensitive and resistant strains across generations.
Gene Ontology and Confirmations: Deletion and GFP-Tagged Library Exhibit Concordant Mechanisms of Toxicity and Detoxification
Following from the enrichment analysis of the sensitive/resistant strains and activated ORFs, we identified four biological categories relevant to BaP toxicity. Namely, we found genes related to cellular oxidative stress, DNA damage and repair, and drug transmembrane transport and detoxification by export to be implicated in the response to BaP in both deletion and GFP-tagged libraries. While the toxicity pathways are similar to what has been identified in prior studies of BaP, the implication of transmembrane transport is novel and could prove useful in identifying future biomarkers of exposure.
Cellular oxidative stress and protein stress
BaP perturbs cellular redox homeostasis by producing quinones (BaP 6,12-quinone, BaP 6,12-hydroquinone, BaP 3,6-hydroquinone, BaP 3,6-quinone; Miller and Ramos, 2001). Further reduction of the quinones creates reactive oxygen species (ROS) in the form of superoxide and hydrogen peroxide (Miller and Ramos, 2001). From our results, we observe a perturbation of the cell redox homeostasis, as indicated by the sensitive deletion strains related to oxygen and radical detoxification (PRX1, TRX2, and AHP1), cellular redox homeostasis (PRX1, TRX2, POR2, AHP1, and GLR1), and enzymes necessary for glutathione – involved in mitigating oxidative stress – recycling (OYE3). It is important to note that the genes responsible for glutathione recycling were only significant in the 5G assay, indicating that NADPH recycling is more important in reducing oxidative stress for earlier generations, as was also found for the quinone metabolites of benzene in a prior yeast study (North et al., 2011). This is further confirmed by the activated oxidative stress ORFs in the GFP-tagged library, such as SLN1 and TRX2 in the 10 mg/L exposure, which imply that oxidative stress plays an important role in BaP toxicity. In addition, protein stress was observed in the GFP-tagged library as indicated by activation of protein folding and refolding categories for unfolded proteins such as HSP104, SSE2, SSA4, UBC8, etc., in five of the six concentrations. However, genes involved in response to protein misfolding were not observed in the deletion library, suggesting that the misfolding damage may be a temporary subsequent response of oxidative or chemical stress or that sufficient redundancy exists to mitigate the requirement for any particular gene product.
DNA damage and repair
Many of the confirmed deletion strains that exhibit sensitivity to BaP are associated with DNA damage and cell cycle disruptions that are concordant with observations in human in vitro studies. Indeed, in S. cerevisiae, RAD6, a ubiquitin-conjugating enzyme also involved in post-replicative repair, RAD50, implicated in the repair of double-stranded DNA breaks and homologous/non-homologous end-joining, DNL4, a DNA ligase for non-homologous end-joining, and ELC1, important in polyubiquitylating Rpb1 during global genome repair, were required for yeast tolerance to BaP (Jentsch et al., 1987; Chen et al., 2001; Hoege et al., 2002; Zhang and Paull, 2005; LeJeune et al., 2009). Non-homologous and homologous end-joining are known as the most important aspect of repairing BaP-induced DNA adducts in humans, so the significance of DNL4 in yeast suggests some concordance in response to BaP (Natarajan and Palitti, 2008). Moreover, other genes in yeast appear to be important in responding to other types of DNA damage caused by ROS (e.g., hydroquinones) such as DTD1, involved in nonsense mutation suppression and ELC1, involved in global genomic repair (von der Haar and Tuite, 2007; LeJeune et al., 2009). While SHM2 is important in yeast for synthesizing DNA bases and its homozygous deletion actually promotes resistance to most other toxicant surveyed in other studies, this was the reverse for BaP (Huang et al., 2008). Indeed, both direct and indirect BaP interference with DNA replication and repair are sources of cellular growth inhibition in yeast (Xie et al., 2003). In the GFP-tagged library, DNA stress response is stronger at lower concentrations, e.g., activation of double strand break repair (RAD51 and HTA2) for a 10 μg/L exposure. For 10 mg/L, only base excision repair (OGG1 and APN2) were activated for base damage repair related to base oxidation, suggesting the DNA damage induced by BaP at higher concentration focuses on bases, possibly related to oxidative damage.
Drug transmembrane transport and detoxification by export
Two yeast strains lacking drug transmembrane transport and detoxification by export, qdr2Δ, and yrm1Δ demonstrated significant growth inhibition in the presence of BaP. Both of these genes encode solute carriers (from family 22) and orthologs were noted to be upregulated in response to bisphenol A and benzo[a]pyrene in the livers of other organisms such as rats and mice as well as human cell cultures (Staal et al., 2006). Moreover, the expression of such genes has been proposed as a potential biomarker of exposure to chemicals and is an important response to chemical injury in human liver cells (Borlark, 2005). In both yeast and people, these solute carriers are specific to various organic anion molecules including glutathione, sulfate, and glucuronide conjugates, which are the ultimate products of BaP metabolism prior to excretion in humans (Pritchard et al., 2005). Interestingly, other solute carrier families have been implicated in human BaP toxicity, such as SLC38A5 which is thought to indirectly favor cis DNA adduct formation by increasing the chloride ion concentration thus catalyzing SN2 attack of the BPDE carbocation metabolite (Wolfe et al., 1994). However, literature concerning the role of solute family 22 in human BaP toxicity is limited and warrants further investigation given the potential role played by the orthologous gene in yeast BaP exposure and hydroquinone metabolite toxicity (Keum et al., 2006). Similarly, in the GFP-tagged assay, the two hour 10 mg/L exposure also demonstrated strong activation of proteins for detoxification by export, including transmembrane transport (such as ATM1, BPT1, QDR2, etc.) and transport (GLK1, TRX2, ATM1, etc.), which supports chemical export as a significant detoxification mechanism.
Lack of Growth Response to S-9 Metabolic Adjuvant and BaP Dose Combinations
The combination of rat S-9 and BaP did not seem to alter growth inhibition in wild type yeast. Only the two percent by volume addition of rat S-9 had any significant effect, but this was only observed at the two lowest concentrations and only resulted in a maximum 1.3-fold increase in toxicity. Our results are consistent with the study by Hakura et al. (2001) that found a minimum concentration of approximately 1.85% S-9 by volume was needed to metabolize BaP. Indeed, 2% S-9 increased toxicity significantly more than would be expected from the 1% concentration, as illustrated by the statistical interaction between BAP and S-9 in our study. Rat S-9 is approximately 100 times more potent than S-9 reagents derived from human microsomes, so we would not expect human S-9 to increase the bioactivation of BaP and toxicity in yeast (Hakura et al., 2001). Very few studies have assayed the utility of S-9 in yeast; indeed, one study involving chlorinated ethylenes found no difference between treatments with and without S-9 (Koch et al., 1988). One possibility for the observed inconsistent effect of S-9 may be that the metabolites produced are more reactive and may not be able to effectively enter the yeast cells. The milieu (i.e., within yeast cells or in solution) where the BaP metabolism reaction occurs may also influence the amount and ratio of the metabolites produced and absorbed by yeast cells (Käppeli, 1986). Moreover, yeast cytochromes produced as part of the biosynthesis pathway for ergosterol, which is an important component of yeast membranes, have been shown to be able to metabolize exogenous chemicals as a side reaction (Alexander et al., 1974; Aoyama and Yoshida, 1978). Specifically, these enzymes are capable of oxidizing compounds to the toxic quinones and epoxides, which are products of BAP metabolism (Wiseman et al., 1978). In addition, these enzymes have also been shown to metabolize pro-carcinogens to their carcinogenic form (Callen et al., 1980). Consequently, while S. cerevisiae does not have the dedicated bona fide metabolism observed in higher organisms, some of its biosynthetic pathways may serve this purpose. Given our inconclusive results on the effect of S-9 on the toxicity of BAP, the known reactivity of the mammalian BAP metabolites, and the demonstrated metabolic action of endogenous yeast enzymes, S-9 may not be necessary for bioactivation in yeast, at least for this particular toxicant. It is also important to note, that the specific metabolites generated by the yeast system may not be identical to mammalian metabolites. Nevertheless, the statistically demonstrable interaction between BAP and S-9 suggests that using S-9 bioactivation, especially at concentrations of at least 2% by volume, may prove useful for other compounds whose independent toxicity is more subtle or ambiguous. Alternatively, anaerobic, high glucose concentrations (4–20%), inhibitors of mitochondrial protein synthesis, or respiration-impaired mutants have been shown to favor cytochrome production in S. cerevisiae and these conditions could also be used as part of a model for investigating other bio-activated toxicants (Wiseman and King, 1982).
Human Ortholog Genes to Yeast Sensitive/Resistant Genes Suggest Possible Biomarkers
From the sensitive and resistant genes, we identified orthologous human genes using the YOGY database (Table 7; Penkett et al., 2006). VDAC1 has been found to be involved in the cellular apoptotic response to BaP in mice, rats, and human cell lines (Huc et al., 2007; Andreau et al., 2012). Perez et al. (2008) did not find RAD50 expression to be upregulated by BaP exposure in normal human keratinocytes, but Brevik et al. (2012) determined that paternal exposure to BaP did up-regulate expression in mouse embryos. In this instance, RAD50 may serve as a biomarker of developmental toxicity and exposure to the toxicant. While SHMT2 has been implicated in tumor genesis in mice, no prior study has linked it to BaP exposure directly (Nilsson et al., 2012). Similarly, LIG4 upregulation has been associated with tobacco-related cancers and not BaP specifically, though BaP is a component of tobacco smoke (Wu et al., 2004). In addition, rodent studies of rats and mice found ortholog upregulation of UBE2B and PRDX5 as well as downregulation of DTD1 (Chen et al., 2001; van Kesteren et al., 2011). While none of the genes in Table 7, with the exception of select solute carriers previously discussed and VDAC1, have been conclusively linked to adverse outcomes in human cells and BaP exposure, they represent candidate susceptibility loci as they are related to carcinogenesis pathways.
Table 7. Select yeast genes with significant response to BaP and their human orthologs.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by the United States National Institute of Environmental Health Sciences Superfund Research Program through grant P42 ES004705 as well as the National Science Foundation (NSF; EEC-0926284 and CAREER CBET-0953633). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors would also like to acknowledge Mani Tagmount for imparting his experience with yeast cultures, Brandon Gaytán for assistance with learning how to use the Funspec analysis instrument, as well as Nina Holland for her contribution to the revision process.
Supplementary Material
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Toxicogenomics_/10.3389/fgene.2012.00316/abstract
Footnotes
References
Alexander, K. T. W., Mitropoulos, K. A., and Gibbons, G. F. (1974). A possible role for cytochrome P-450 during the biosynthesis of zymosterol from lanosterol by Saccharomyces cerevisiae. Biochem. Biophys. Res. Commun. 60, 460–467.
Andreau, K., Melanie, L., and Aida, B. (2012). Health and cellular impact of air pollutants: from cytoprotection to cytotoxicity. Biochem. Res. Int. 2012, 1–18.
Aoyama, Y., and Yoshida, Y. (1978). Interaction of lanosterol to cytochrome P-450 purified from yeast microsomes: evidence for contribution of cytochrome P-450 to lanosterol metabolism. Biochem. Biophys. Res. Commun. 82, 33–38.
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29.
Aust, A. E., Falahee, K. J., Maher, V. M., and McCormick, J. J. (1980). Human cell-mediated benzo(a)pyrene cytotoxicity and mutagenicity in human diploid fibroblasts. Cancer Res. 40, 4070–4075.
Azari, M. R., and Wiseman, A. (1981). Purification and characterization of the cytochrome P-448 component of a benzo(a)pyrene hydroxylase from Saccharomyces cerevisiae. Anal. Biochem. 122, 129–138.
Baird, W. M., Hooven, L. A., and Mahadevan, B. (2005). Carcinogenic polycyclic aromatic hydrocarbon-DNA adducts and mechanism of action. Environ. Mol. Mutagen. 45, 106–114.
Bartosiewicz, M., Penn, S., and Buckpitt, A. (2001). Applications of gene arrays in environmental toxicology: fingerprints of gene regulation associated with cadmium chloride, benzo(a)pyrene, and trichloroethylene. Environ. Health Perspect. 109, 71–74.
Borlark, J. (2005). Handbook of Toxicogenomics: A Strategic View of Current Research and Applications. Hannover: Wiley-VCH.
Brevik, A., Birgitte, L., Vendula, R., Ann-Karin, O., Gunnar, B., and Nur, D. (2012). Paternal benzo[a]pyrene exposure affects gene expression in the early developing mouse embryo. Toxicol. Sci. 129, 157–165.
Burchiel, S. W., Thompson, T. A., Lauer, F. T., and Oprea, T. I. (2007). Activation of dioxin response element (DRE)-associated genes by benzo(a)pyrene 3,6-quinone and benzo(a)pyrene 1,6-quinone in MCF-10A human mammary epithelial cells. Toxicol. Appl. Pharmacol. 221, 203–214.
Callen, D. F., Wolf, C. R., and Philpot, R. M. (1980). Cytochrome P-450 mediated genetic activity and cytotoxicity of seven halogenated aliphatic hydrocarbons in Saccharomyces cerevisiae. Mutat. Res. 77, 55–63.
Chen, L., Trujillo, K., Ramos, W., Sung, P., and Tomkinson, A. E. (2001). Promotion of Dnl4-catalyzed DNA end-joining by the Rad50/Mre11/Xrs2 and Hdf1/Hdf2 complexes. Mol. Cell 8, 1105–1115.
Chung, S. Y., Yettella, R. R., Kim, J. S., Kwon, K., Kim, M. C., and Min, D. B. (2011). Effects of grilling and roasting on the levels of polycyclic aromatic hydrocarbons in beef and pork. Food Chem. 129, 1420–1426.
Dabrowska, D., Kot-Wasik, A., and Namiesnik, J. (2008). Stability studies of selected polycyclic aromatic hydrocarbons in different organic solvents and identification of their transformation products. Polish J. Environ. Stud. 17, 17–24.
Edoardo, M. (1992). Urban air pollution by polycyclic aromatic hydrocarbons: levels and sources of variability. Sci. Total Environ. 116, 109–135.
Fromme, H., Lahrz, T., Piloty, M., Gebhardt, H., Oddoy, A., and Rüden, H. (2004). Polycyclic aromatic hydrocarbons inside and outside of apartments in an urban area. Sci. Total Environ. 326, 143–149.
Giaever, G., Chu, A. M., Ni, L., Connelly, C., Riles, L., Véronneau, S., et al. (2002). Functional profiling of the Saccharomyces cerevisiae genome. Nature 418, 387–391.
Hakura, A., Suzuki, S., Sawada, S., Motooka, S., and Satoh, T. (2001). An improvement of the Ames test using a modified human liver S9 preparation. J. Pharmacol. Toxicol. Methods 46, 169–172.
Hattemer-Frey, H. A., and Travis, C. C. (1991). Benzo-a-Pyrene: environmental partitioning and human exposure. Toxicol. Ind. Health 7, 141–157.
Hockley, S. L., Mathijs, K., Staal, Y. C., Brewer, D., Giddings, I., van Delft, J. H., and Phillips, D. H. (2009). Interlaboratory and interplatform comparison of microarray gene expression analysis of HepG2 cells exposed to benzo(a)pyrene. OMICS 13, 115–125.
Hoege, C., Pfander, B., Moldovan, G. L., Pyrowolakis, G., and Jentsch, S. (2002). RAD6-dependent DNA repair is linked to modification of PCNA by ubiquitin and SUMO. Nature 419, 135–141.
Huang, B., Lu, J., and Byström, A. S. (2008). A genome-wide screen identifies genes required for formation of the wobble nucleoside 5-methoxycarbonylmethyl-2-thiouridine in Saccharomyces cerevisiae. RNA 14, 2183–2194.
Huc, L., Xavier, T., Jørn, A., Holme, M. R., Anita, S., Claire, G., et al. (2007). c-Jun NH2-terminal kinase-related Na+/H+ exchanger isoform 1 activation controls hexokinase II expression in benzo(a)pyrene-induced apoptosis. Cancer Res. 67, 1696–1705.
Huh, W.-K., James, F. V., Luke, G. C., Adam, S. C., Russell, W. H., Jonathan, S. W., et al. (2003). Global analysis of protein localization in budding yeast. Nature 425, 686–691.
Janikowska, G., and Wardas, W. (2002). Bioconcentration of benzo(a)pyrene in Clorella BB cells. Polish J. Environ. Stud. 11, 345–348.
Jentsch, S., McGrath, J. P., and Varshavsky, A. (1987). The yeast DNA repair gene RAD6 encodes a ubiquitin-conjugating enzyme. Nature 329, 131–134.
Jo, W. J., Alex, L., Henry, W., Michelle, C., Allan, H. S., Dave, K., et al. (2009a). Comparative funtional genomic analysis identifies distinct and overlapping sets of genes required for resistance to monomethylarsonous acid (MMAIII) and arsenite (ASIII) in yeast. Toxicol. Sci. 111, 424–436.
Jo, W. J., Ren, X., Chu, F., Aleshin, M., Wintz, H., Burlingame, A., et al. (2009b). Acetylated H4K16 by MYST1 protects UROtsa cells from arsenic toxicity and is decreased following chronic arsenic exposure. Toxicol. Appl. Pharmacol. 241, 294–302.
Kapitulnik, J., Wislocki, P. G., Levin, W., Yaqi, H., Jerina, D. M., and Conney, A. H. (1978). Tumorigenicity studies with diol-epoxides of benzo(a)pyrene which indicate that (±)-trans-7β,8α-Dihydroxy-9α,10α-epoxy-7,8,9,10-tetrahydrobenzo(a)pyrene is an ultimate carcinogen in newborn mice. Cancer Res. 38, 354–358.
Keum, Y.-S., Han, Y.-H., Liew, C., Kim, J.-H., Xu, C., Yuan, X., et al. (2006). Induction of heme oxygenase-1 (HO-1) and NAD[P]H: quinone oxidoreductase 1 (NQO1) by a phenolic antioxidant, butylated hydroxyanisole (BHA) and its metabolite, tert-butylhydroquinone (tBHQ) in primary-cultured human and rat hepatocytes. Pharm. Res. 23, 2586–2594.
Kim, K. B., and Lee, B. M. (1997). Oxidative stress to DNA, protein, and antioxidant enzymes (superoxide dismutase and catalase) in rats treated with benzo(a)pyrene. Cancer Lett. 113, 205–212.
King, D. J., Wiseman, A., and Wilkie, D. (1983). Studies on the genetic regulation of cytochrome P-450 production in Saccharomyces cerevisiae. Mol. Gen. Genet. 192, 466–470.
Koch, R., Schlegelmilch, R., and Wolf, H. U. (1988). Genetic effects of chlorinated ethylenes in the yeast Saccharomyces cerevisiae. Mutat. Res. 206, 209–216.
Kultz, D. (2005). Molecular and evolutionary basis of the cellular stress response. Annu. Rev. Physiol. 67, 225–257.
Laurent, C., Feidt, C., Grova, N., Mpassi, D., Lichtfouse, E., Laurent, F., et al. (2002). Portal absorption of 14C after ingestion of spiked milk with 14C-phenanthrene, 14C-benzo[a]pyrene or 14C-TCDD in growing pigs. Chemosphere 48, 843–848.
LeJeune, D., Chen, X., Ruggiero, C., Berryhill, S., Ding, B., and Li, S. (2009). Yeast Elc1 plays an important role in global genomic repair but not in transcription coupled repair. DNA Repair (Amst.) 8, 40–50.
Loguinov, A., Mian, I. S., and Vulpe, C. D. (2004). Exploratory differential gene expression analysis in microarray experiments with no or limited replication. Genome Biol. 5, R18.
Mewes, H. W., Frishman, D. U., Güldener, M. G., Mayer, K., Mokrejs, M., Morgenstern, B. M., et al. (2002). MIPS: a database for genomes and protein sequences. Nucleic Acids Res. 30, 31–34.
Miller, K. P., and Ramos, K. S. (2001). Impact of cellular metabolism on the biological effects of benzo[a]pyrene and related hydrocarbons. Drug Metab. Rev. 33, 1–35.
Mitchell, C. E. (1982). Distribution and retention of benzo(a)pyrene in rats after inhalation. Toxicol. Lett. 11, 35–42.
Moody, R. P., Joncas, J., Richardson, M., and Chu, I. (2007). Contaminated soils (I): in vitro dermal absorption of benzo[a]pyrene in human skin. J. Toxicol. Environ. Health A 70, 1858–1865.
Natarajan, A. T., and Palitti, F. (2008). DNA repair and chromosomal alterations. Mutat. Res. 657, 3–7.
Nilsson, L. M., Tacha Zi, P. F., Sara, R., Christiane, K., Walter, P., Georg, W. B., et al. (2012). Mouse genetics suggests cell-context dependency for myc-regulated metabolic enzymes during tumorigenesis. PLoS Genet. 8:1–11. doi:10.1371/journal.pgen.1002573
North, M., Tandon, V. J., Thomas, R., Loguinov, A., Gerlovina, I., Hubbard, A. E., et al. (2011). Genome-wide functional profiling reveals genes required for tolerance to benzene metabolites in yeast. PLoS ONE 6:e24205. doi:10.1371/journal.pone.0024205
Penkett, C. J., Morris, J. A., Wood, V., and Bähler, J. (2006). YOGY: a web-based, integrated database to retrieve protein orthologs and associated Gene Ontology terms. Nucleic Acids Res. 34, W330–W334.
Perez, D. S., Robert, J. H., Raymond, Y. S. H., and Julie, A. C. (2008). Gene expression changes associated with altered growth and differentiation in benzo[a]pyrene or arsenic exposed normal human epidermal keratinocytes. J. Appl. Toxicol. 28, 491–508.
Pierce, S. E., Davis, R. W., Nislow, C., and Giaever, G. (2007). Genome-wide analysis of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures. Nat. Protoc. 2, 2958–2974.
Pritchard, J. B., Lowes, S., and Miller, D. S. (2005). “The molecular basis of xenobiotic transport and metabolism in choroid plexus,” in The Blood-Cerebrospinal Fluid Barrier, eds A. Chodobski and W. Zheng (Boca Raton, FL: CRC Press), 147–173.
Ren, X., Aleshin, M., Jo, W. J., Dills, R., Kalman, D. A., Vulpe, C. D., et al. (2011). Involvement of N-6 adenine-specific DNA methyltransferase 1 (N6AMT1) in arsenic biomethylation and its role in arsenic-induced toxicity. Environ. Health Perspect. 119, 771–777.
Robinson, M., Grigull, J., Mohammad, N., and Hughes, T. (2002). FunSpec: a web-based cluster interpreter for yeast. BMC Bioinformatics 3:35. doi:10.1186/1471-2105-3-35
Romero, D. L., Mounho, B. J., Lauer, F. T., Born, J. L., and Burchiel, S. W. (1997). Depletion of glutathione by benzo(a)pyrene metabolites, ionomycin, thapsigargin, and phorbol myristate in human peripheral blood mononuclear cells. Toxicol. Appl. Pharmacol. 144, 62–69.
Roseiro, L. C., Gomes, A., and Santos, C. (2011). Influence of processing in the prevalence of polycyclic aromatic hydrocarbons in a Portuguese traditional meat product. Food Chem. Toxicol. 49, 1340–1345.
Shu, H. P., and Bymun, E. N. (1983). Systemic excretion of benzo(a)pyrene in the control and microsomally induced rat: the influence of plasma lipoproteins and albumin as carrier molecules. Cancer Res. 43, 485–490.
Simmons, S. O., Fan, C., and Ramabhadran, R. (2009). Cellular stress response pathway system as a sentinel ensemble in toxicological screening. Toxicol. Sci. 111, 202–225.
Staal, Y. C. M., van Herwijnen, M. H. M., van Schooten, F. J., and van Delft, J. H. (2006). Modulation of gene expression and DNA adduct formation in HepG2 cells by polycyclic aromatic hydrocarbons with different carcinogenic potencies. Carcinogenesis 27, 646–655.
Stribinskis, V., and Ramos, K. S. (2006). Activation of human long interspersed nuclear element 1 retrotransposition by benzo(a)pyrene, an ubiquitous environmental carcinogen. Cancer Res. 66, 2616–2620.
van Kesteren, P. C., Zwart, P. E., Pennings, J. L., Gottschalk, W. H., Kleinjans, J. C., van Delft, J. H., et al. (2011). Deregulation of cancer-related pathways in primary hepatocytes derived from DNA repair-deficient Xpa−/−p53+/− mice upon exposure to benzo[a]pyrene. Toxicol. Sci. 123, 123–132.
Vayssier-Taussat, M., Camilli, T., Aron, Y., Meplan, C., Hainaut, P., Polla, B. S., et al. (2001). Effects of tobacco smoke and benzo[a]pyrene on human endothelial cell and monocyte stress responses. Am. J. Physiol. Heart Circ. Physiol. 280, H1293–H1300.
von der Haar, T., and Tuite, M. F. (2007). Regulated translational bypass of stop codons in yeast. Trends Microbiol. 15, 78–86.
Vondrácek, J., Kozubík, A., and Machala, M. (2002). Modulation of estrogen receptor-dependent reporter construct activation and G0/G1–S-phase transition by polycyclic aromatic hydrocarbons in human breast carcinoma MCF-7 cells. Toxicol. Sci. 70, 193–201.
Wang, J., Qiang, H., Zhi-Ping, L., Yong, W., Ling-Yun, W., Luonan, C., et al. (2011). NOA: a novel network ontology analysis method. Nucleic Acids Res. 39, e87.
Wang, X., and Wang, W.-X. (2006). Bioaccumulation and transfer of benzo(a)pyrene in a simplified marine food chain. Mar. Ecol. Prog. Ser. 312, 101–111.
Wester, R. C., Maibach, H. I., Bucks, D. A., Sedik, L., Melendres, J., Liao, C., et al. (1990). Percutaneous absorption of [14C]DDT and [14C]benzo[a]pyrene from soil. Fundam. Appl. Toxicol. 15, 510–516.
Wiseman, A., and King, D. J. (1982). “Microbial oxygenases and their potential application,” in Topics in Enzyme and Fermentation Biotechnology, ed. A. Wiseman (Chichester: Ellis Horwood Ltd.), 151–206.
Wiseman, A., Lim, T. K., and Woods, L. F. J. (1978). Regulation of the biosynthesis of cytochrome P-450 in brewers [sic] yeast. Role of cyclic AMP. Biochim. Biophys. Acta 544, 615–623.
Wolfe, A. R., Yamamoto, J., and Meehan, T. (1994). Chloride ions catalyze the formation of cis adducts in the binding of anti-benzo[a]pyrene diol epoxide to nucleic acids. Proc. Natl. Acad. Sci. U.S.A. 91, 1371–1375.
Wu, X., Hua, Z., Rebecca, S., and David, C. C. (2004). Genetic susceptibility to tobacco-related cancer. Oncogene 23, 6500–6523.
Xie, Z., Braithwaite, E., Guo, D., Zhao, B., Geacintov, N. E., and Wang, Z. (2003). Mutagenesis of benzo[a]pyrene diol epoxide in yeast: requirement for DNA polymerase ζ and involvement of DNA polymerase η. Biochemistry 42, 11253–11262.
Zhang, L., McHale, C. M., Rothman, N., Li, G., Ji, Z., Vermeulen, R., et al. (2010). Systems biology of human exposure. Chem. Biol. Interact. 184, 86–93.
Zhang, X., and Paull, T. T. (2005). The Mre11/Rad50/Xrs2 complex and non-homologous end-joining of incompatible ends in S. cerevisiae. DNA Repair 4, 1281–1294.
Appendix
Figure A1. Growth curve assay for BY4743 wild type S. cerevisiae exposed to increasing concentrations of BAP in YPD liquid media. Optical density at 595 nm was measured every 15 min for a period of 24 h; each point in the curve represents a mean of three replicates. Standard error is omitted for clarity.
Figure A2. BAP Dose-response curve. Each “bp” green x mark is a point estimate of three biological replicate exposures at five concentrations: 100, 200, 300, 400, and 500 μm. The dotted line represents 95% confidence interval to the fitted line. See Figure A1 for growth curves and total growth data for each of the five doses.
Figure A3. Confirmation of inhibition of selected yeast strains exposed individually to BAP. *Represents significance at 0.01<p<0.05; **represents significance at p<0.01. Bars are standard errors.
Keywords: benzo[a]pyrene, toxicity, yeast, stress, resistance, sensitivity, S-9, ontology
Citation: O’Connor STF, Lan J, North M, Loguinov A, Zhang L, Smith MT, Gu AZ and Vulpe C (2013) Genome-wide functional and stress response profiling reveals toxic mechanism and genes required for tolerance to benzo[a]pyrene in S. cerevisiae. Front. Gene. 3:316. doi: 10.3389/fgene.2012.00316
Received: 01 November 2012; Paper pending published: 17 November 2012;
Accepted: 20 December 2012; Published online: 08 February 2013.
Edited by:
Michael Aschner, Vanderbilt University Medical Center, USAReviewed by:
Garry Wong, University of Eastern Finland, FinlandMarcelo Farina, Universidade Federal de Santa Catarina, Brazil
Copyright: © 2013 O’Connor, Lan, North, Loguinov, Zhang, Smith, Gu and Vulpe. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Chris Vulpe, Vulpe Laboratory, Department of Nutritional Science and Toxicology, 119 Morgan Hall, University of California, Berkeley, CA 94720, USA. e-mail: vulpe@berkeley.edu