Ozone inhalation modifies the rat liver proteome.

Ozone (O3) is a serious public health concern. Recent findings indicate that the damaging health effects of O3 extend to multiple systemic organ systems. Herein, we hypothesize that O3 inhalation will cause downstream alterations to the liver. To test this, male Sprague-Dawley rats were exposed to 0.5ppm O3 for 8h/day for 5 days. Plasma liver enzyme measurements showed that 5 day O3 exposure did not cause liver cell death. Proteomic and mass spectrometry analysis identified 10 proteins in the liver that were significantly altered in abundance following short-term O3 exposure and these included several stress responsive proteins. Glucose-regulated protein 78 and protein disulfide isomerase increased, whereas glutathione S-transferase M1 was significantly decreased by O3 inhalation. In contrast, no significant changes were detected for the stress response protein heme oxygenase-1 or cytochrome P450 2E1 and 2B in liver of O3 exposed rats compared to controls. In summary, these results show that an environmentally-relevant exposure to inhaled O3 can alter the expression of select proteins in the liver. We propose that O3 inhalation may represent an important unrecognized factor that can modulate hepatic metabolic functions.

Introduction

Approximately 50% of the U.S. population resides in areas where ambient ozone (O3) concentrations exceed the current 0.075 ppm 8 h time average set by the National Ambient Air Quality Standards (NAAQS) [1]. O3 is a primary component of photochemical smog and is a serious public health concern especially when air quality is poor. Persons particularly sensitive to O3 exposure include the elderly, young children, and those with pre-existing pulmonary diseases such as asthma and chronic obstructive pulmonary disease. It is established that exposure to high levels of O3 decreases lung function, increases pulmonary hyper-reactivity, causes airway epithelial cell damage and remodeling, and increases epithelial permeability [2-4].

In addition to adverse pulmonary effects, emerging evidence shows that O3 inhalation can cause tissue injury and altered metabolism in other systemic organ systems. For example, studies in nonhuman primates show that postnatal episodic O3 exposure only during infancy (30 days to 6 months of age) resulted in long-term effects to the systemic innate immune system that were detectable up to 1 yr of age [5]. Recent epidemiologic studies report a positive correlation of air pollution exposure with increased cardiovascular-related morbidity and mortality resulting in increased hospital admissions related to cardiac events [6-9]. Additional studies reported that successive days of high ambient O3 exposure correlated with increased blood pressure, blood lipids, and decreased glucose tolerance in humans [8]. Ballinger, Postlethwait, and colleagues showed that exposure to 0.5 ppm O3 induced vascular dysfunction and increased aortic mitochondrial DNA damage in healthy wild-type mice, and increased progression of atherosclerosis in mice genetically predisposed to cardiovascular disease [10]. Moreover, a single 5-day regimen of O3 also increased mtDNA damage in the abdominal aorta of infant nonhuman primates [10]; a model more closely mimicking human exposures. This work is important because not only does it provide direct evidence that O3 inhalation damages the cardiovascular system, but also shows that O3 has the potential to exacerbate cardiovascular disease including atherosclerosis [11].

O3 inhalation has also been shown to induce metabolic changes in other peripheral organ systems. For example, pentobarbital-induced sleeping times increased after O3 exposure, suggesting alterations in hepatic drug metabolism and clearance mechanisms [12,13]. Since this early work, other groups have addressed the effect of inhaled O3 on the liver. Last and colleagues, using a microarray approach, examined the effects of O3 (1.0 ppm) on the liver transcriptome and showed that O3 inhalation significantly decreased the mRNA levels of several xenobiotic, carbohydrate, and fatty acid metabolism genes [14]. Other groups have shown that acute exposures of O3 (1.0–2.0 ppm) for 3 h increased nitric oxide production in isolated hepatocytes, as well as increased rates of protein synthesis [15]. Herein, our goal was to determine whether inhalation of O3 at a concentration lower than used in previous studies, alters the hepatic proteome and specific xenobiotic metabolizing enzymes in an animal model without pre-existing disease.

Materials and methods

O3 exposure protocol

Male Sprague-Dawley rats were purchased from Harlan Laboratories (Barrier 217 VAF, Indianapolis, IN) and were provided standard rat chow and water ad libitum. Animals were housed two per cage under barrier conditions and maintained using a standard 12 h light-dark cycle. Rats were exposed to either filtered air (FA) or ozone (O3) at 0.5 ppm for 8 h/day for 5 days between 9:00 AM to 5:00 PM in the University of Alabama at Birmingham (UAB) Environmental Exposure Facility. Standard laboratory rat chow was removed from cages before the start of exposures to prevent rats from ingesting food that may contain oxidized nutrients (e.g., lipids, thiols, antioxidants) as a consequence of O3 reaction. Rat chow was also removed from cages of FA exposed rats. Fresh food was returned to cages after exposures and present throughout the entire dark period when rats are active and consume food. The exposure protocol had no effect on body weight (FA: 348±4 g and O3: 341±15 g, p=0.69). Prior to exposures, animals were acclimated to the chambers for at least 72 h before beginning FA and O3 exposures to minimize stress associated with a novel environment. Importantly, rats remained in their home cage with bedding and water, with cages (stainless steel wire mesh coverings) placed into the exposure chamber. O3 was generated from 100% O2 using a model OZ1PCS-V/SW O3 generator (Ozotech Inc, Yreka, CA), mixed with FA, and flowed into 0.8 m3 stainless steel chambers (~22 °C, 50% relative humidity; 30 volume changes/hr). O3 concentrations were continuously monitored using a Thermo-Environmental Model UV Photometric analyzer (Thermo-Environmental, Franklin, MA). Measurements obtained inside the cages demonstrated that O3 concentrations in the animals breathing zones were equivalent to the chamber bulk phase concentrations. FA controls were housed in separate chambers receiving the same clean air supply at equal air exchange rates throughout the duration of the exposure regimen. Animals were anesthetized with a 50 mg/kg body weight (i.p.) injection of sodium pentobarbital and euthanized via exsanguination. Tissues were harvested within 1 h of cessation of the exposure. All procedures were approved by the University of Alabama at Birmingham Institutional Animal Care and Use Committee.

Bronchoalveolar lavage and cell differential analysis

A cannula was inserted into the trachea via a midline tracheotomy, the chest cavity opened via a midline thoracotomy, and 9 mL of warmed (37 °C) phosphate buffered saline (pH 7.0; 310 mOsm) was gently instilled and withdrawn 3 times to yield bronchoalveolar lavage (BAL) [16]. The BAL was subjected to a standard cytospin procedure (StatSpin Cytofuge, Norwood, MA) to determine cell counts and types. Slides were stained with a modified Wright-Giemsa stain and 300 cells (100 cells per lane) were counted under a light microscope using characteristics unique to each cell type [17]. Total cell counts were averaged for each cell type and percent of total cells calculated. The remaining BAL was centrifuged to generate BAL fluid (BALF) and protein concentrations in the BALF were measured to determine whether altered epithelial permeability [18] was detectable after the 5th period of O3.

Plasma chemistries for liver enzymes and histology

Blood was collected via the abdominal aorta using a heparin-coated syringe and centrifuged at 4 °C at 2000g for 10 min to obtain plasma. Alanine and aspartate aminotransferase (ALT and AST, respectively) activities were measured in plasma using a spectrophotometric assay per manufacturer's directions (Pointe Scientific, Inc, Canton MI). Briefly, NAD+ oxidation was recorded under constant temperature conditions (37 °C) over a time period of 10 min at 340 nm. The rate was measured and ALT and AST levels reported as international units/liter (IU/L). Liver tissue was fixed in formalin and paraffin embedded. Sections were mounted on slides and stained with hematoxylin and eosin (H&E). Slides were scored for injury, steatosis, and inflammation by a pathologist blinded to the experimental groups.

Two dimension isoelectric focusing/SDS-PAGE (2D IEF/SDS-PAGE)

Livers were excised and homogenized in ice-cold 0.25 M sucrose, 5 mM Tris–HCl, and 1 mM EDTA, pH 7.4, containing protease inhibitors [phenylmethylsulfonyl fluoride (40 µg/mL), leupeptin (5 µg/mL), and pepstatin A (7 µg/mL)] [19]. Protease inhibitors were included to prevent sample degradation prior to proteomic analyses. For proteomic studies there were six FA (control) and six O3 exposed rats per group. Post-nuclear supernatant fraction was prepared by centrifugation of liver homogenates at 568g for 10 min at 4°C. Protein concentrations were determined using the Bradford protein assay and bovine serum albumin (BSA) as a standard [20]. Proteomic analyses were performed by methods as previously described [21]. Post-nuclear supernatant (100 µg) from liver homogenates was added to IEF gel strip rehydration buffer containing 7 M urea, 2 M thiourea, 2% (w/v) CHAPS, 0.5% (w/v) n-dodecyl-ß-D-maltoside, 0.002% (w/v) bromophenol blue, ampholine electrophoresis reagent (Sigma, St. Louis, MO, range pH 3–10), 0.04 M DTT and 2 mM tributylphosphine. Following protein extraction, samples were applied to IEF gel strips (Invitrogen ZOOM Strips, pH 3–10, Carlsbad, CA) and rehydration of IEF strips was done overnight. For SDS-PAGE, IEF gel strips were placed horizontally on top of a 10% resolving gel with 4% stacking gel, and sealed into place using warm agarose (1%, w/v) and gels were run at 100 V for . After electrophoresis, gels were stained with Sypro Ruby™ (Invitrogen, Carlsbad, CA) for total protein. Protein stained gels were imaged using a Bio-Rad ChemiDoc XRS imaging system (Bio-Rad Laboratories, Inc, Hercules, CA).

2D gel image analyses

Methods for gel analyses are described in our previous work [21]. Differences in protein density were performed using PDQuest Image Analysis software (Bio-Rad, Hercules, CA). Individual protein spots on all twelve 2D gels (n=6 for FA and n=6 for O3) were identified using the program software and visually matched for accuracy. A master gel was created to serve as the reference gel for FA and O3 groups, which the PDQuest software program uses automatically to match the spots across all gels. Manual verification was then used to correct for incorrect spot matching to the reference gel or spots not initially detected by software. In order to correct for any differences in protein loading, the density for all spots in each gel was normalized to total density in verified protein spots for that particular gel. The protein densities were then entered into a spreadsheet and statistically significant differences were determined using an unpaired t-test (two-tail).

Protein identification with matrix assisted laser desorption ionization time of flight (MALDI-TOF) mass spectrometry

Protein identification was performed as previously described [21]. Protein spots were cut from gels and processed using standard methods in the UAB mass spectrometry shared facility (http://www.uab.edu/proteomics). Briefly, protein samples were de-stained with three successive 30 min washes in 1:1 solution of 50 mM NH4HCO3 and acetonitrile. Samples were treated with 10 mM dithiothreitol in 50 mM NH4HCO3 for 60 min at 60 °C to reduce cysteine residues, which was followed by alkylation of free cysteines with 55 mM iodoacetamide in 50 mM NH4HCO3 for 60 min. This was followed by 16 h incubation at 37 °C with trypsin to digest proteins. The resulting solution was extracted by two successive 30 min washes in a 1:1 solution of 5% formic acid and acetonitrile, supernatants were collected, and dried using a Savant SpeedVac. Lyophilized samples were resuspended in 0.1% formic acid, desalted (C18 ZipTips, Millipore), and diluted 1:10 with a saturated solution of α-cyano-4-hydroxycinnamic acid matrix before being applied to MALDI-TOF plates. Samples were analyzed with a Voyager De-Pro mass spectrometer in the positive mode and resulting spectra were analyzed using Voyager Explorer software. Peptide masses were submitted to the Mascot database for protein identification (www.matrixscience.com). Mascot uses a probability based system of the Mowse algorithm to determine significant results [22]. Mowse is an algorithmic calculation used to assign statistical weight to each peptide match; thus, a higher Mowse score implies higher statistical likelihood of the match being correct. Herein, the Mascot results for proteins IDs ranged from 132–467. All scores were>40 indicating extensive homology and significant matches for this study (p≤0.05).

Western blotting

Based on gene expression results from other studies [12-14] we decided to examine the effect of inhaled O3 on select cytochrome P450 (CYP450) enzymes and the stress-inducible enzyme heme oxygenase 1 (HO-1) in the liver. Equal amounts of protein from post-nuclear supernatant fractions (40 or 50 µg) were loaded and separated on 12% SDS-PAGE gels then transferred for 1 h to nitrocellulose membranes. Membranes were blocked overnight in 5% (w/v) non-fat milk in Tris base saline solution with Tween-20 (TBS-T) while shaking at 5 °C. Levels of HO-1 protein were detected using a 1:5,000 dilution of primary antibody (Stressgen, Ann Arbor, MI). Levels of CYP2E1 and CYP2B were detected using 1:1,000 primary antibody dilutions (Chemicon International, Temecula, CA). After successive washes in TBS-T, blots were incubated in horseradish peroxidase (HRP) conjugated anti-rabbit secondary antibody (1:10,000) in 5% (w/v) non-fat milk in TBS-T for 1 h (Sigma, St. Louis, MO). Bands were visualized using luminol-based chemiluminescence reagents with a ChemiDoc XRS System (Bio-Rad, Hercules, CA) and quantification of band densities were achieved using Quantity One software (Bio-Rad, Hercules, CA) [23]. To ensure equal protein loading of lanes, GAPDH was used as a loading control at a dilution of 1:1000 anti-GAPDH primary antibody. A representative example of a CYP2B blot reprobed for GAPDH is shown in Fig. 2D. Further, equal protein loading of gels was confirmed using total protein staining (e.g., SYPRO® Ruby protein gel stain: FA=34,875±4241 and O3=36,079±3337 arbitrary density units; n=6, p=0.59) for blots not reprobed with GAPDH.

Fig. 2

Protein levels of HO-1, CYP2E1, and CYP2B in liver of O3 and FA rats. Rats were exposed to filtered air (FA) or O3 (0.5 ppm) for 8 h/day for 5 days. After exposures, livers were removed and homogenates were analyzed via western blot technique for (A) HO-1, (B) CYP2E1, (C) CYP2B, and (D) GAPDH as a loading control. O3 exposure did not significantly alter protein expression for all 4 proteins measured (p>0.05). Full-size images of western blots were cropped to highlight protein bands of interest. Continuous gel images are shown with approximate molecular weights provided next to blot image. Data represent mean±SD for n=6 animals per group. HO-1, heme oxygenase-1; CYP2E1, cytochrome P450 2E1; CYP2B, cytochrome P450 2B; and GAPDH, glyceraldehyde 3-phosphate dehydrogenase.

Statistical analyses

Statistical significance for all measured parameters was assessed using an unpaired t-test and significance was determined with a p value≤0.05.

Results

Lung cell differential and BALF protein

BALF from FA and O3 exposed rats was assessed for epithelial permeability and inflammatory cell infiltration as described in Methods. No significant change in BALF protein concentration were observed between FA and O3 exposed groups; 63.0±27.9 and 72.9±46.2 μg/mL, respectively (p>0.05). Cell differential analysis also revealed that macrophage, lymphocyte, neutrophil, or basophil airspace populations in O3 exposed rats did not differ from FA controls at the end of the 5 days of O3 exposure protocol (Table 1).

Table 1

Bronchoalveolar lavage cell differential analysis.

	Filtered air	Ozone	p value
Volume recovered (mL)	6.54±0.52	6.23±0.48	0.32
% Yield	73±7.00	69±7.00	0.32
Cell counta
Macrophages	99.2±0.76	97.5±1.95	0.10
Neutrophils	1.30±1.37	2.20±0.68	0.75
Lymphocytes	0.33±0.82	0.42±1.90	0.23
Basophils	0.00±0.00	0.17±0.40	0.34

Lavage fluid was recovered as described in Methods.

Cell counts were generated using the average of 3 independent counts of 100 total cells/count and are represented as a percentage of total counted cells.

Effect of O3 exposure on plasma liver enzymes and histology

Plasma ALT and AST enzyme activities were measured as indicators of hepatocellular damage due to O3 exposure. After 5 days of 0.5 ppm O3 exposure, AST enzyme activity levels for FA and O3 groups were 47.2±2.3 and 55.2±13.2 IU/L, respectively. Similarly, ALT enzymes activity levels for FA and O3 groups were 48.9±3.7 and 51.0±8.9 IU/L, respectively. No statistically significant changes were observed in liver enzyme activities in the O3 treated group compared to the FA group (p>0.05). Additionally H&E-stained liver sections scored for changes to liver morphology, inflammation, and injury showed no pathology and no differences between O3 and FA groups (data not shown). These results demonstrate that this inhaled O3 exposure regimen did not cause detectable toxicity to the liver.

Effect of ozone exposure on liver proteome

Representative 2D IEF/SDS-PAGE gel images from both FA and O3 treatment groups are shown in Fig. 1. A global 2D proteomic approach was performed to determine whether O3 inhalation altered liver protein abundance. Following separation of proteins in two dimensions (isoelectric point and molecular weight), proteins were analyzed for differences in abundance using PDQuest software tools as described in Methods. Using this mini-gel 2D system, we were able to detect approximately 275 individual proteins spots on each gel. O3 exposure resulted in a significant change in the expression of 23 proteins as compared to FA control. Of these 23 proteins, mass spectrometry analyses identified 10 proteins. A master map showing the location of the identified proteins and their assigned spot numbers is provided in Fig. 1A. A comparison of the protein densities of these 10 proteins is shown graphically in Fig. 1B. Unique protein identifications with statistical analysis (% change and p values) and mass spectrometry details are provided in Table 2. Each protein spot shown in the master map matched the expected molecular weight and isoelectric point of each identified protein (Table 2). Total protein density from 2D gels is presented in Fig. 1C and shows that changes in individual protein density were not due to differences in total protein loading of gels. The 10 identified proteins were categorized into 4 broad groups: cytoskeletal, energy metabolism, drug metabolism, and protein folding/ER stress (Table 3).

Fig. 1

Master map of liver proteins differentially altered by inhaled O3. Rats were exposed to filtered air (FA) or O3 (0.5 ppm) for 8 h/day for 5 days. After exposures, livers were removed and the post-nuclear supernatant was analyzed for global protein expression using 2D gel proteomics and mass spectrometry techniques. (A) Shows the representative 2D gels of proteins from liver of FA and O3 exposed rats. Note that the circled protein “spots” were identified and are listed in Table 2 with respective p values. (B) Shows the change in abundance (increase or decrease) of liver proteins found to be altered by inhaled O3. The number below each pair of bar graphs corresponds to the protein spot shown in the master map (panel A) and match numbers used in Table 2. (C) Total protein density calculated from 2D gels generated from FA and O3 groups. This result shows equal protein loading across all gels (p>0.05). Data represent mean±SD of n=6 animals per group.

Table 2

Hepatic proteins significantly altered in abundance as a result of 0.5 ppm O3 exposure: results from 2D IEF/SDS-PAGE and mass spectrometry.

Spot #	Protein identification	Mass (kDa)	Isoelectric point (pI)	% Change	p value	MOWSE score	Peptides matched
1	Cytochrome b5 type A (microsomal)	11.4	5.26	41.6	0.038	166	6
2	Catechol-O-methyltransferase	29.8	5.41	50.9	0.025	132	5
3	Galactokinase 1	42.8	5.24	29.7	0.036	286	8
4	Actin, beta like 2	42.2	5.30	20.6	0.014	277	13
5	Protein disulfide-isomerase	57.3	4.82	32.3	0.002	784	34
6	Tubulin, alpha 1C	50.6	4.96	41.1	0.003	336	16
7	Glucose-regulated protein 78	72.5	5.07	52.0	0.002	467	15
8	Glycerol kinase	58.2	5.49	40.6	0.01	220	7
9	Glutathione S-transferase mu 1	25.8	8.27	−27.5	0.021	350	15
10	D-β-hydroxybutyrate dehydrogenase, type I	38.7	8.93	−31.6	0.033	187	8

Proteins identified as significantly altered were matched to all Filtered Air (FA) and Ozone (O3) 2D IEF/SDS-PAGE gels. Spot # is the same number used to identify circled proteins in Fig. 1A and densities shown in Fig. 1B. Data represent p values determined using a Student's unpaired t-test (two-tail).

Table 3

Description of cellular pathways and function of identified hepatic proteins altered O3 inhalation.

Protein Identification	Accession #	Pathway	Function
Cytochrome b5	GenBank:AAB67609	Drug metabolism	Cypb5 mediates the rate of P450 dependent mono-oxygenation reactions through the transfer of the second electron from NADPH.
Catechol-O-methyltransferase	GenBank:NP_036663	Drug metabolism/L-dopa metabolism	COMT enzymatically O-methylates catechol containing compounds. Both over expression and loss of activity has been implicated in vascular diseases as well as neuronal disorders.
Galactokinase 1	GenBank:NP_001008283	Cell cycle, energy metabolism	First enzyme required for the conversion of galactose to glucose forming galactose-1-phosphate.
Actin, beta like 2	GenBank:NP_001099879	Cytoskeletal	Important for cell motility, structure, and integrity.
Protein disulfide-isomerase	GenBank:NP_037130	Protein Folding/ER stress	PDI is responsible for inter and intra molecular disulfide bond formation within the endoplasmic reticulum.
Tubulin, alpha 1C	GenBank:NP_001011995	Cytoskeletal	Alpha tubulin is one of two components of tubulin which is assembled to form microtubules. These are important for cell movement and maintenance of shape.
Glucose-regulated protein 78	GenBank:NP_037215	Protein Folding/ER stress	Located within the ER; GRP78/BiP binds and inactivates the proteins responsible for UPR activation. Also directly binds misfolded proteins to assist in proper folding or clearance.
Glycerol kinase	GenBank:NP_077357	Energy metabolism	Catalyzes formation of glycerol 3-phosphate. Overexpression in H4IIE cells resulted in increased fat storage and alters activity of PPAR-alpha and other transcription factors.
Glutathione S-transferase mu 1	GenBank:AAA41286	Phase II drug metabolism	Facilitates conjugation of reduced glutathione to electrophilic compounds for increased excretion.
D-β-hydroxybutyrate dehydrogenase	GenBank:AAB59684	Energy/lipid metabolism	Oxidoreductase, mitochondrial. Involved in fatty acid catabolism. BDH1 facilitates formation of the ketone body 3-hydroxybutyrate.

Protein folding/ER stress proteins

Two proteins involved in ER stress were increased by O3 exposure. Protein disulfide isomerase (PDI) increased by 32% and glucose-regulated protein 78 (GRP78) increased by 52% in liver of O3 exposed rats compared to protein levels measured in FA controls. PDI assists in the proper folding and disulfide bond formation of proteins within the ER [24]. GRP78, a chaperone protein located in the ER, plays an important role in the regulation of the unfolded protein response activated during ER stress [25].

Drug metabolism proteins

Proteins that demonstrated significant expression changes due to O3 exposure compared to FA control samples included microsomal cytochrome b5, catechol-O-methyltransferase (COMT), and glutathione-S-transferase mu 1 (GSTM 1). Cytochrome b5 increased in expression by 43% in liver of O3 exposed rats compared to FA controls. This protein is important for CYP450 enzymatic reactions involving fatty acid desaturation, metabolism of xenobiotics, and cholesterol synthesis and breakdown [26]. COMT protein, a highly abundant liver protein, increased by 50% in liver of O3 exposed rats. COMT catalyzes O-methylation of catechol containing molecules [27] and has been implicated in pathologies such as vascular disease and Parkinson's disease due to the generation of reactive catecholamine semiquinone species [28]. GSTM 1 decreased in expression by 27% in liver of O3 exposed rats compared to FA controls. GSTM 1 is a phase II xenobiotic metabolism enzyme responsible for the conjugation of glutathione groups to endogenous and exogenous compounds for excretion [29]. This enzyme is largely considered a hepatoprotective enzyme due to its ability to conjugate reactive metabolites produced as byproducts from oxidative stress mechanisms [29].

Energy metabolism and cytoskeletal proteins

Three proteins involved in energy metabolism, galactokinase, glycerol kinase, and β-hydroxybutyrate dehydrogenase 1 (BDH1) were increased by O3 exposure. Both galactokinase and glycerol kinase increased in expression by 41% and 30%, respectively, in response to O3, whereas BDH significantly decreased in expression by 32% in liver of O3 exposed rats compared to FA controls. BDH1 is a mitochondrial protein that catalyzes the first step of ketone body metabolism, converting β-hydroxybutyrate to acetoacetate [30]. Cytoskeletal proteins that significantly increased in response to O3 exposure were β-actin and α-tubulin, which are involved in cell cycle, division, and motility [31].

Effect of O3 exposure on CYP450 and HO-1 protein levels

Previous studies showed that high O3 exposure levels (1–3 ppm) prolonged pentobarbital-induced sleeping times in rodents, suggesting a decrease in xenobiotic metabolism [13,32]. Herein, protein levels of CYP2E1 and 2B were measured to determine if lower, more environmentally relevant exposures of O3 (0.5 ppm) altered key drug metabolizing enzymes. After 5 days of 0.5 ppm O3 exposure, protein expression levels of CYP2E1 or CYP2B were not significantly changed as compared to FA controls (Fig. 2B and C). Additionally, the stress inducible form of heme oxygenase, HO-1, was measured in this study. Analysis showed no significant difference in HO-1 protein levels between FA and O3 groups at the end of 5 exposure periods (Fig. 2A).

Discussion

A new disease paradigm has emerged in the field of environmental health sciences suggesting that inhalation of environmental pollutants, like O3, can exacerbate diseases that extend beyond the pulmonary system; the initial site of exposure. Importantly, O3 inhalation may initiate a cascade of events leading to tissue perturbations in other organ systems, including cardiovascular, hematopoietic, and hepatic systems. In support of this new concept, particulate matter (PM) and other gaseous pollutants generated from industry and diesel exhaust exacerbate not only pulmonary injury, but also cardiovascular disease [6]. Epidemiological studies also report increased cardiovascular risk from high ambient O3 days with significant changes in blood lipids, blood pressure, and factors important for glucose metabolism [8]. While human exposures to inhaled pollutants typically involve heterogeneous mixtures of PM and O3, animal studies show that O3 alone can damage cardiovascular tissues and cause hepatic changes [10,14]. Importantly, in nonhuman primates exposure to O3 also induces cardiovascular injury [10] and induces long-term alterations in immune function [5].

The molecular events responsible for O3 induced systemic effects are poorly defined. The reaction of O3 within the lung is complex and results in the production of a diverse array of biologically active products. O3 reacts with mono-unsaturated fatty acids to produce multiple aldehyde and hydroxyhydroperoxide species, as well as small amounts of the Criegee ozonide [33,34]. Studies by Murphy, White, and colleagues demonstrated that the reaction of O3 with calf lung surfactant resulted in the formation of biologically active oxidized phospholipid species that reduced cell viability by both necrotic and apoptotic mechanisms [35]. Moreover, Friedman and colleagues have demonstrated that multiple lipid ozonation products activate cellular signaling pathways and induce inflammatory mediators [36].

As most air pollution episodes last several days, single event exposures in experimental laboratory studies are not environmentally relevant. Most exposures are also not continuous but occur for several hours per day for several consecutive days. For these reasons, we designed an O3 exposure model that encompassed these features; i.e., 8 h per day for 5 consecutive days. Further, the exposure regimen used in our current study better replicates a typical day where ambient O3 concentrations have increased beyond the current NAAQS standard of 0.075 ppm. The concentration of O3 used in the current study is 0.5 ppm, which is higher than the NAAQS standard. Notably, most experimental animal studies use O3 concentrations in the 0.8–3.0 ppm range to examine health impacts. Higher doses of O3 are used in experimental studies for a variety of reasons [11]. First, it is generally accepted that because rats are obligate nose breathers and that they have a more complex nasopharyngeal architecture that facilitates removal of inhaled reactive gases; concentrations above those encountered by humans during typical ambient exposures should be employed. Second, O3 exposures with animals are typically done during the day when rats are less active/asleep (i.e., the light period of a 12 h light:12 h dark cycle). Consequently, the combination of nose breathing and low physical activity function to reduce the dose delivered to the lower respiratory track of rats. These factors will ultimately result in a lower composite dose of O3 reaching the lungs in experimental animals as compared to typical human O3 exposures.

Following a 5 day O3 exposure protocol, we did not observe any significant changes to overall protein levels or inflammatory cells within BALF extracted from lungs (Table 1). Increased protein in BALF would indicate significant epithelial damage and loss of barrier integrity. It has been shown that repeated O3 exposures increase neutrophils, lymphocytes, and macrophages in the lungs. In rodents, neutrophilic inflammation generally begins to resolve by day 3 following a 1.0 ppm O3 exposure [37]. We also did not detect changes in plasma levels of ALT or AST enzyme activities; common measurements for hepatocyte necrosis. Taken together, these findings indicate that by the 5th day of O3 exposure, injury may have resolved sufficiently to reduce injury metrics below detectable levels. Similarly, we observed no change in CYP2E1 and 2B protein levels following O3 exposure, which supports recent results showing no change in gene expression for these two CYP450 enzymes following inhaled O3 exposure [14]. In addition to looking at CYP450s, we measured expression of another key stress responsive protein HO-1; the inducible member of the heme-oxygenase enzyme family. Herein, we saw no change in HO-1 protein after a 5 day O3 exposure regimen. The absence of any change in HO-1 expression could be anticipated considering HO-1 up-regulation typically occurs very early in response to stress. Further, it is important to note that we measured HO-1 levels only at one time point. If O3 exposure elicited an early hepatic stress response (e.g., HO-1 induction) at any time before 5 days it would have been missed in this study.

Despite the absence of overt lung and liver toxicity at Day 5, O3 exposure did alter the levels of other key proteins in the liver as evidenced by our proteomics analyses. Considering that only a few studies have investigated the effects of inhaled O3 on the liver [13-15] we took a global proteomic approach to determine whether O3 exposure alters expression of liver proteins. Moreover, we propose that investigating biological events and/or alterations that occur after pulmonary inflammation has resolved is a more relevant question for understanding human health impacts from inhaled pollutants. Specifically, it is more important to determine whether systemic alterations and/or pathological events induced by O3 are still obvious after pulmonary inflammation has resolved versus simply looking at downstream systemic sequelae that occur at the same time as the toxic insult.

Using 2D gel electrophoresis techniques and mass spectrometry we detected over 20, and successfully identified 10, proteins that were significantly altered by O3 inhalation. These proteins were grouped into 4 different categories: cytoskeletal, drug metabolism, energy metabolism, and ER stress-related proteins; with 8 proteins significantly increased and 2 decreased in response to inhaled O3 exposure. Two of the most intriguing changes observed from our mass spectrometry data are the increases in PDI and GRP78. Both proteins are located in the endoplasmic reticulum and are involved in the proper folding of proteins and maintenance of ER homeostasis. The ER is a complex intracellular membranous network acting as the site for proper folding and maturation of newly synthesized proteins. It is also sensitive to a variety of signals and tightly regulated in order to maintain proper protein and lipid homeostasis. One such signal is excess accumulation of misfolded proteins, which activates the unfolded protein response (UPR). GRP78 binds all three proteins responsible for the activation of UPR (IRE1α, ATF6, and PERK), rendering the signaling pathway inactive [24]. GRP78 is believed to be released from the previous three proteins responsible for UPR activation and is sequestered in the ER lumen where it binds misfolded proteins and subsequently the UPR is activated [38]. An upregulation of GRP78 suggests that a mild form of ER stress occurred by the Day 5 time point in order to bring the ER lumen back to homeostasis. Indeed, overexpression of GRP78 protects cells from known stressors such as a high fat diet; alleviating injury and promoting restoration of ER homeostasis [39,40]. Additionally, studies investigating the effect of PM exposure show induction of ER stress pathways in the liver of animal models with pre-existing liver disease [41].

Despite no detectable effects on CYP2E1 and 2B protein levels we did observe changes in other drug metabolism enzymes. We saw a 50% increase in COMT protein levels. COMT, a phase II drug metabolism enzyme, O-methylates catechol-containing compounds including dietary phytochemicals and metabolites of toxic aryl hydrocarbons [28]. Increases in COMT may be a protective adaptive response because accumulation of endogenous catecholamines can spontaneously or enzymatically be converted to semiquinone and quinone intermediates that are cytotoxic [28,42]. In contrast, we saw a decrease in liver GSTM 1 in response to O3 exposure. The GST super-family consists of phase II drug metabolism enzymes and GSTM 1 is one of seven GST M isoforms [29]. Substrates for GSTs include environmental chemicals, like acrolein and DDT, as well as, numerous drugs and endogenous reactive lipid species [43]. Conjugation of GSH with these compounds typically produces less reactive products that are subsequently excreted; thus, ameliorating cellular oxidative stress and toxicity. GST enzymes are typically increased in response to stress and considered a conserved adaptive response as it occurs throughout the animal kingdom [44]. Interestingly, we observed that hepatic GSTM 1 protein was significantly decreased by O3 inhalation. While several studies have reported that the GSTM 1 deficiency is associated with decreased lung function and airway inflammation following O3 exposure [45]; the impact of O3 and GSTM 1 deficiency on the liver is not known. Animal models null for GSTM 1 display altered hepatic drug clearance [46]. In the context of the current study, decreased GSTM 1 protein in the liver after O3 inhalation could impair adaptive stress responses and/or impair hepatic drug metabolism and clearance mechanisms. Importantly, this could manifest as a serious problem in individuals with GSTM 1 polymorphisms exposed to O3. Indeed, individuals with the GSTM 1 null polymorphism are at higher risk for hepatocellular carcinoma, alcoholic liver disease, and obesity associated fatty liver disease [47-50]. A proteomic analysis of alcohol-induced liver injury found that chronic, binge ethanol administration decreased GSTM 1 levels in liver [51]. Therefore, the failure of key adaptive stress responses may be critical in how liver responds and adapts to O3 [52]. Finally, β-actin and α-tubulin, were found to be altered by O3. Adduction of these proteins by reactive lipid species and acetaldehyde disrupts the cytoskeletal network in hepatocytes [31] with increases in transcript and protein observed in an animal model of pancreatic injury [53]. It was postulated that these increases play a protective role by increasing the amount of cytoskeletal protein needed to replace damaged or modified β-actin and α-tubulin in the cell. Taken together, these results suggest that inhaled O3 alters protein expression patterns in an organ (liver) distal from initial site of exposure.

Conclusion

In summary, we report that an environmentally-relevant inhalation exposure to O3 alters the liver proteome. We observed significant changes in the expression of key stress response and drug metabolism proteins in the liver following O3 inhalation. While it is well established that O3 can exacerbate pulmonary diseases, an emerging literature indicates that inhaled O3 may have significant effects on other organ systems, including the liver. Thus, it will be important to better understand how these systemic effects could impact health. Importantly, the systemic effects observed with O3 alone may be exacerbated when coupled to pre-existing diseases and pathologies. For example, animals predisposed to heart disease display greater vascular injury when exposed to O3 [10]. In addition, rodents administered acetaminophen, a known hepatotoxicant, have increased liver injury following O3 exposure [54]. Also, genetic predisposition to obesity in multiple animal models modulates airway hypersensitivity and inflammatory cytokines as compared to unexposed obese controls [55-57]. Together, the results from the current study provides new insight into the growing field of environmental-related disease and health impacts from inhaled pollutants and suggest that the liver is a systemic site impacted by O3.

Funding

This work was supported by the National Institute of Health grants to SMB (R01 AA15172 and R01 AA18841) and EMP (P01 ES111617 and R01 HL54696).