Animal treatment. Adult (25 weeks old) male mice of MutaMouse strain (transgenic mouse strain 40.6) were exposed to BaP (Sigma Aldrich, Canada) as described previously (Malik
et al., 2012 (
link)). In brief, animals were dosed daily via oral gavage for 28 days with varying doses of BaP (0, 25, 50, and 75mg/kg body weight/day) dissolved in olive oil. Each dose group contained five animals. Mice were sacrificed by cardiac puncture under isoflorane anesthesia 72h following the final exposure. The right lobe of the lung was excised, flash frozen in liquid nitrogen, and stored at −80°C until use. For the duration of the experiment, food (2014 Teklad Global standard rodent diet) and water were provided
ad libitum, and mice were caged individually in plastic film isolators (Harlan Isotec, U.K.) on a 12-h light/12-h dark cycle. Mice were bred, maintained, and treated in accordance with the Canadian Council for Animal Care Guidelines and approved by Health Canada’s Animal Care Committee.
Tissue selection. The major exposures to BaP occur via the oral route (drinking and feed) (Hettemer-Frey and Travis, 1991 (
link)). The organs directly affected by BaP as a consequence of oral exposure include stomach, esophagus, tongue, and larynx (Culp
et al., 1998 (
link)). However, studies conducted by Stoner
et al. (1984) (
link) and Wattenberg and Leong (1970) (
link) showed lung and liver to be equally impacted by BaP (oral gavage) and, moreover, revealed sensitivity of lung for tumor development in comparison with liver. In alignment with these reports, our previous work (Halappanavar
et al., 2011 (
link); Yauk
et al., 2011 (
link)) revealed lung-specific regulation of the biological processes known to be associated with cancer formation. Thus, in the present study, we have investigated transcriptional responses in lung and liver. The results of the study may also help address the importance of considering the multiorgan toxicity in calculating the risk associated with chemicals, such as BaP.
Tissue DNA extraction. The frozen lung and liver tissues were sliced randomly. Genomic DNA was isolated from a random tissue section for measuring the levels of DNA adducts and transgene mutant frequency. In brief, lung tissue was minced and degassed to remove all traces of air in the alveoli. Livers were homogenized in ice cold TMST buffer (50mM Tris, pH 7.6, 3mM magnesium acetate, 250mM sucrose, 0.2% Triton X-100) as described in the study by Douglas
et al. (1994) (
link). The minced tissue was washed twice in cold PBS and lysed in 10mM Tris, pH 7.6, 10µM EDTA, 100µM NaCl, and 1% SDS overnight at 37°C on a rotating platform. The lysate was digested with proteinase K (1mg/ml lysis buffer). DNA was isolated using a serial phenol/chloroform/isoamyl alcohol (25:24:1) and chloroform/isoamyl alcohol (24:1) extraction (Renault
et al., 1997 (
link)). DNA was precipitated in ethanol and dissolved in Tris-EDTA buffer (10mM Tris, pH 7.6, 500mM EDTA). DNA was stored at 4°C until use.
DNA adduct analysis. DNA adduct formation in each sample was determined using the nuclease P1 digestion enrichment version of the
32P-postlabeling assay as described previously (Phillips and Arlt, 2007 (
link)) with minor modifications. Briefly, 4 µg of DNA was digested overnight with micrococcal nuclease (288 mUnits, Sigma, cat. no. N3755) and calf spleen phosphodiesterase (1.2 mUnits, MP Biomedicals, cat. no. 100977), enriched, and labeled as described elsewhere (Phillips and Arlt, 2007 (
link)). Radiolabeled adducted nucleotide biphosphates were separated by thin-layer chromatography on polyethyleneimine-cellulose plates (Macherey-Nagel, Düren, Germany) with the following chromatographic conditions (Arlt
et al., 2008 (
link)): D1, 1.0M sodium phosphate, pH 6; D3, 4.0M lithium formate, 7.0M urea, pH 3.5; D4, 0.8M LiCl, 0.5M Tris, 8.5M urea, pH 8. Chromatographs were scanned using a Packard Instant Imager (Canberra Packard, Downers Grove USA), and DNA adduct levels (relative adduct labeling) were calculated from the adduct counts per minute (cpm), the specific activity of [γ-
32P]ATP, and the amount of DNA (pmol of DNA-P) used. An external BPDE-DNA standard was used for identification of BaP-DNA adducts. Results are expressed as DNA adducts/10
8 nucleotides.
LacZ
mutant frequency (positive selection). The transgene
lacZ mutant frequency in lungs was determined using the P-gal (phenyl-β-D-galactopyranoside) positive selection assay as described in the study by Vijg and Douglas (1996) and Lambert
et al. (2005) (
link). The λgt10
lacZ DNA was rescued from the genomic DNA using the Transpack lambda packaging system (Stratagene). The packaged phage particles were mixed with host bacterium (
Escherichia coli lacZ
−,
galE
−,
recA
−, pAA119 with
galT and
galK), plated on minimal medium containing 0.3% (w/v) P-Gal, and incubated overnight at 37°C. Total plaque-forming units (pfu) were measured on concurrent titers that did not contain P-Gal. Mutant frequency is expressed as the ratio of the number of mutant pfu to total pfu. Mutant frequency data analysis was performed as described previously (Malik
et al., 2012 (
link)).
RNA extraction and purification. Total RNA was isolated for gene expression analysis and qRT-PCR validation as described previously (Halappanavar
et al., 2011 (
link)). Briefly, total RNA was extracted from the lungs using TRIzol reagent (Invitrogen) and purified using RNeasy Mini Kit (Qiagen, Canada). The RNA quantity and purity were checked using a NanoDrop Spectrophotometer (ThermoFisher Scientific, Canada). The RNA integrity was determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Canada). The samples showing A
260/A
280 ratios between 2.1 and 2.2 and having RNA integrity number above 7.5 were used for further analysis.
Microarray hybridization and analysis. Double-stranded cDNA and cyanine-labeled cRNA were synthesized (Agilent Linear Amplification Kits, Agilent Technologies) from 250ng of total RNA from each sample and universal reference total RNA (Stratagene, Canada). Cyanine-labeled cRNA targets were
in vitro transcribed using T7 RNA polymerase and purified by RNeasy Mini Kit (Qiagen). From each (sample and reference) labeled sample, 825ng of cRNA was hybridized to Agilent 4×44K oligonucleotide microarrays (Agilent Technologies) at 60°C overnight (16h) in the Agilent SureHyb hybridization chamber. Arrays were washed and scanned on an Agilent G2505B Scanner according to manufacturer’s recommendations. Feature extraction software version 10.7.3.1 (Agilent Technologies) was used to extract the data.
A reference design (Kerr and Churchill, 2001 (
link)) was used to analyze mRNA expression as described previously (Malik
et al., 2012 (
link)). All analyses were conducted in the R (R-Development-Core-Team, 2010 ) environment using the MAANOVA library (Wu, 2010 ). The background fluorescence was measured using the negative control (−)3×SLv1 probes; probes with median signal intensities less than the trimmed mean (trim = 5%) plus three trimmed standard deviations of (−)3×SLv1 probe were flagged as absent (within the background signal). Probes were considered present if at least four of the five samples within a condition had signal intensities greater than three trimmed standard deviations above the trimmed mean of the (−)3×SLv1 probes (background signal). Data were normalized using the transform.madata() function using the glowess option with a span of 0.1. Ratio intensity plots and heat maps for the raw and normalized data were constructed to identify outliers. One sample (50mg/kg group) was removed from the analysis based on clustering. The statistical model for this analysis included fixed effects of array and treatment condition and was applied to the log
2 of the absolute intensities. Differentially expressed transcripts (upregulated or downregulated relative to olive oil–treated control mouse lung samples) were determined using the Fs statistic option in the matest() function. The
p values for all statistical tests were estimated by the permutation method with residual shuffling, and false discovery rate (FDR) adjusted
p values were estimated using the adjPval() function. The fold change calculations were estimated as described previously (Malik
et al., 2012 (
link)). Significant genes were selected based on a FDR adjusted
p value < 0.05 for any BaP exposed versus control contrast.
qRT-PCR array validation. Mouse pathway-specific PCR array (cancer-PAMM-033, SABiosciences) and custom PCR arrays consisting of 172 genes in total were employed to validate the microarray results. Genes for the custom array were selected based on their implication in biological processes relevant to lung carcinogenesis. These genes included statistically significant differentially expressed genes (FDR adjusted
p < 0.05), differentially expressed genes that exhibited high fold changes (fold rank only) but were not statistically significant, and genes that were differentially regulated (FDR adjusted
p < 0.05) in the livers from the same mice (Malik
et al., 2012 (
link)). In brief, 0.8 µg of total lung tissue RNA (
n = 5/group) from each sample was reverse transcribed using RT
2 First Strand Kit (SABiosciences). Real-time PCR was performed using RT
2 SYBR Green PCR Master Mix on a CFX96 real-time detection system (Bio-Rad, Canada). Threshold cycle values for each well were averaged. Relative gene expression was determined according to the comparative Ct method and normalized to reference RNAs
Hprt1,
Gapdh, and
Actb housekeeping genes for the Cancer Pathway array and to reference RNAs
Gusb and
Hprt1 housekeeping genes for the custom arrays. Transcripts were further normalized by subtracting the median delta Ct value for each sample. Differential expression was determined with a two-sample bootstrap test using R software (R-Development-Core-Team, 2010 ). The fold change was estimated using the ratio of the arithmetic mean of the treated sample to the mean of the control samples. Standard errors for the fold change values were estimated using the bootstrap test (Efron and Tibshirani, 1993 ).
Bioinformatics. All mRNA data are deposited in the NCBI Gene Expression Omnibus database under accession numbers GSE35718 (lung) and GSE24910 (liver). Following normalization, biological functions perturbed in response to BaP were identified using functional annotation clustering in the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (Huang da
et al., 2009 (
link)) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto, 2000 (
link)). The biological and molecular functions of genes that were significantly differentially expressed (FDR adjusted
p ≤ 0.05 or FDR adjusted
p ≤ 0.1 and fold change ≥ 1.5) following exposure to BaP treatment were analyzed and explored in Ingenuity Pathway Analysis (IPA, Ingenuity Systems, Redwood City, CA). Molecular relational networks of genes modulated by BaP in lung tissue enriched for cancer function were generated using IPA. Each molecule was overlaid onto a global molecular network developed from information contained in the Ingenuity Knowledge Base. Networks were generated based on their connectivity. All relationships are supported by at least one reference from the literature.
Labib S., Yauk C., Williams A., Arlt V.M., Phillips D.H., White P.A, & Halappanavar S. (2012). Subchronic Oral Exposure to Benzo(a)pyrene Leads to Distinct Transcriptomic Changes in the Lungs That Are Related to Carcinogenesis. Toxicological Sciences, 129(1), 213-224.