Sensitivity of the Indo-Pacific coral Acropora millepora to aromatic hydrocarbons.

The risks posed by petroleum spills to coral reefs are poorly understood and quantifying acute toxicity thresholds for aromatic hydrocarbons to reef-building corals is required to assess their sensitivity relative to other taxa. In this study, we exposed Acropora millepora to toluene, naphthalene and 1-methylnaphthalene (1-MN) in a flow-through system and assessed survivorship and sublethal responses including growth, colour and the photosynthetic performance of symbionts. Median 50% lethal concentrations (LC50s) decreased over the 7-d exposure period, reaching asymptotic values of 22,921, 5,268, 1167 µg L-1 for toluene, naphthalene and 1-MN, respectively. Corresponding toxicokinetic parameters (εLC50) defining the time progression of toxicity were 0.830, 0.692, and 0.256 d-1, respectively. Latent effects after an additional 7-d recovery in uncontaminated seawater were not observed. Effect concentrations (EC50s) for 50% growth inhibition were 1.9- to 3.6-fold lower than the LC50s for each aromatic hydrocarbon. There were no observed effects of aromatic hydrocarbon exposure on colour score (a proxy for bleaching) or photosynthetic efficiency. Acute and chronic critical target lipid body burdens (CTLBBs) of 70.3 ± 16.3 and 13.6 ± 18.4 µmol g-1 octanol (± standard error) were calculated for survival and growth inhibition based on 7-d LC50 and EC10 values, respectively. These species-specific constants indicate adult A. millepora is more sensitive than other corals reported so far but is of average sensitivity in comparison with other aquatic taxa in the target lipid model database. These results advance our understanding of acute hazards of petroleum contaminants to key habitat-building tropical coral reef species.

Adopting a toxic unit model paradigm in design, analysis and interpretation of oil toxicity testing.

The lack of a conceptual understanding and unifying quantitative framework to guide conduct and interpretation of laboratory oil toxicity tests, has led investigators to divergent conclusions that can confuse stakeholders and impede sound decision-making. While a plethora of oil toxicity studies are available and continue to be published, due to differences in experimental design, results between studies often cannot be compared. Furthermore, much resulting data fails to advance quantitative effect models that are critically needed for oil spill risk and impact assessments. This paper discusses the challenges posed when evaluating oil toxicity test data based on traditional, total concentration-based exposure metrics and offers solutions for improving the state of practice by adopting a unifying toxic unit (TU) model framework. Key advantages of a TU framework is that differences in test oil composition, sensitivity of the test organism/endpoint, and toxicity test design (i.e., type of test) can be taken into quantitative account in predicting aquatic toxicity. This paradigm shift is intended to bridge the utility of laboratory oil toxicity tests with improved assessment of effects in the field. To illustrate these advantages, results from literature studies are reassessed and contrasted with conclusions obtained based on past practice. Using instructive examples, model results are presented to explain how dissolved oil composition and concentrations and resulting TUs vary in WAFs prepared using variable loading or dilution test designs and the important role that unmeasured oil components contribute to predicted oil toxicity. Model results are used to highlight how the TU framework can serve as a valuable aid in designing and interpreting empirical toxicity tests and provide the data required to validate/refine predictive toxicity models. To further promote consistent exposure and hazard assessment of physically and chemically dispersed oil toxicity tests recommendations for advancing the TU framework are presented.

Bridging the lab to field divide: Advancing oil spill biological effects models requires revisiting aquatic toxicity testing.

Oil fate and exposure modeling addresses the complexities of oil composition, weathering, partitioning in the environment, and the distributions and behaviors of aquatic biota to estimate exposure histories, i.e., oil component concentrations and environmental conditions experienced over time. Several approaches with increasing levels of complexity (i.e., aquatic toxicity model tiers, corresponding to varying purposes and applications) have been and continue to be developed to predict adverse effects resulting from these exposures. At Tiers 1 and 2, toxicity-based screening thresholds for assumed representative oil component compositions are used to inform spill response and risk evaluations, requiring limited toxicity data, analytical oil characterizations, and computer resources. Concentration-response relationships are employed in Tier 3 to quantify effects of assumed oil component mixture compositions. Oil spill modeling capabilities presently allow predictions of spatial and temporal compositional changes during exposure, which support mixture-based modeling frameworks. Such approaches rely on summed effects of components using toxic units to enable more realistic analyses (Tier 4). This review provides guidance for toxicological studies to inform the development of, provide input to, and validate Tier 4 aquatic toxicity models for assessing oil spill effects on aquatic biota. Evaluation of organisms' exposure histories using a toxic unit model reflects the current state-of the-science and provides an improved approach for quantifying effects of oil constituents on aquatic organisms. Since the mixture compositions in toxicity tests are not representative of field exposures, modelers rely on studies using single compounds to build toxicity models accounting for the additive effects of dynamic mixture exposures that occur after spills. Single compound toxicity data are needed to quantify the influence of exposure duration and modifying environmental factors (e.g., temperature, light) on observed effects for advancing use of this framework. Well-characterized whole oil bioassay data should be used to validate and refine these models.

Temporal chemical composition changes in water below a crude oil slick irradiated with natural sunlight.

Photooxidation can alter the environmental fate and effects of spilled oil. To better understand this process, oil slicks were generated on seawater mesocosms and exposed to sunlight for 8 days. The molecular composition of seawater under irradiated and non-irradiated oil slicks was characterized using ion mobility spectrometry-mass spectrometry and polyaromatic hydrocarbons analyses. Biomimetic extraction was performed to quantify neutral and ionized constituents. Results show that seawater underneath irradiated oil showed significantly higher amounts of hydrocarbons with oxygen- and sulfur-containing by-products peaking by day 4-6; however, concentrations of dissolved organic carbon were similar. Biomimetic extraction indicated toxic units in irradiated mesocosms increased, mainly due to ionized components, but remained <1, suggesting limited potential for ecotoxicity. Because the experimental design mimicked important aspects of natural conditions (freshly collected seawater, natural sunlight, and relevant oil thickness and concentrations), this study improves our understanding of the effects of photooxidation during a marine oil spill.

Oil Irradiation Experiments Document Changes in Oil Properties, Molecular Composition, and Dispersant Effectiveness Associated with Oil Photo-Oxidation.

While chemical dispersants are a powerful tool for treating spilled oil, their effectiveness can be limited by oil weathering processes such as evaporation and emulsification. It has been suggested that oil photo-oxidation could exacerbate these challenges. To address the role of oil photo-oxidation in dispersant effectiveness, outdoor mesocosm experiments with crude oil on seawater were performed. Changes in bulk oil properties and molecular composition were quantified to characterize oil photo-oxidation over 11 days. To test relative dispersant effectiveness, oil residues were evaluated using the Baffled Flask Test. The results show that oil irradiation led to oxygen incorporation, formation of oxygenated hydrocarbons, and higher oil viscosities. Oil irradiation was associated with decreased dispersant efficacy, with effectiveness falling from 80 to <50% in the Baffled Flask Test after more than 3 days of irradiation. Increasing photo-oxidation-induced viscosity seems to drive the decreasing dispersant effectiveness. Comparing the Baffled Flask Test results with field data from the Deepwater Horizon oil spill showed that laboratory dispersant tests underestimate the dispersion of photo-oxidized oil in the field. Overall, the results suggest that prompt dispersant application (within 2-4 days), as recommended by current oil spill response guidelines, is necessary for effective dispersion of spilled oil.

Toxicity of two representative petroleum hydrocarbons, toluene and phenanthrene, to five Atlantic coral species.

Coral reefs are keystone coastal ecosystems that can be exposed to petroleum hydrocarbons from multiple sources, and when selecting spill response methods to limit environmental damages, corals represent one of the highest valued resources for protection. Because previous research to characterize the sensitivity of coral species to petroleum hydrocarbon exposures is limited, a continuous-flow passive dosing system and toxicity testing protocol was designed to evaluate the acute effects of two representative petroleum compounds, toluene and phenanthrene, on five coral species: Acropora cervicornis, Porites astreoides, Siderastera siderea, Stephanocoenia intersepta, and Solenastrea bournoni. Using analytically confirmed exposures, sublethal and lethal endpoints were calculated for each species, and used as model inputs to determine critical target lipid body burdens (CTLBBs) for characterizing species sensitivity. Further, quantification of the time-dependent toxicity of single hydrocarbon exposures is described to provide model inputs for improved simulation of spill impacts to corals in coastal tropical environments.

Miniaturised marine tests as indicators of aromatic hydrocarbon toxicity: Potential applicability to oil spill assessment.

Assessing oil spill toxicity in real time is challenging due to dynamic field exposures and lack of simple, rapid, and sensitive tests. We investigated the relative sensitivity of two commercially available marine toxicity tests to aromatic hydrocarbons using the target lipid model (TLM). State of the art passive dosing in sealed vials was used to assess the sensitivity of brine shrimp (Artemia franciscana) and rotifer (Brachionus plicatilis). Organisms were exposed to toluene, 1-methylnaphthalene and phenanthrene for 24 h. Toxicity results were analysed using the TLM to estimate the critical target lipid body burden and support comparison to empirical data for 79 other aquatic organisms. Our findings demonstrate the applicability of passive dosing to test small volumes and indicate that the two rapid cyst-based assays are insensitive in detecting hydrocarbon exposures compared to other aquatic species. Our results highlight the limitations of applying these tests for oil pollution monitoring and decision-making.

Assessing toxicity of hydrophobic aliphatic and monoaromatic hydrocarbons at the solubility limit using novel dosing methods.

Reliable delineation of aquatic toxicity cut-offs for poorly soluble hydrocarbons is lacking. In this study, vapor and passive dosing methods were applied in limit tests with algae and daphnids to evaluate the presence or absence of chronic effects at exposures corresponding to the water solubility for representative hydrocarbons from five structural classes: branched alkanes, mono, di, and polynaphthenic (cyclic) alkanes and monoaromatic naphthenic hydrocarbons (MANHs). Algal growth rate and daphnid immobilization, growth and reproduction served as the chronic endpoints investigated. Results indicated that the dosing methods applied were effective for maintaining mean measured exposure concentrations within a factor of two or higher of the measured water solubility of the substances investigated. Chronic effects were not observed for hydrocarbons with an aqueous solubility below approximately 5 µg/L. This solubility cut-off corresponds to structures consisting of 13-14 carbons for branched and cyclic alkanes and 16-18 carbons for MANHs. These data support reliable hazard and risk evaluation of hydrocarbon classes that comprise petroleum substances and the methods described have broad applicability for establishing empirical solubility cut-offs for other classes of hydrophobic substances. Future work is needed to understand the role of biotransformation on the observed presence or absence of toxicity in chronic tests.

The sensitivity of the deepsea species northern shrimp (Pandalus borealis) and the cold-water coral (Lophelia pertusa) to oil-associated aromatic compounds, dispersant, and Alaskan North Slope crude oil.

This study investigated the sensitivity of two deepsea species using mortality of northern shrimp (Pandalus borealis) and polyp activity of stony coral (Lophelia pertusa) to dispersant, Corexit 9500 and aromatic hydrocarbons (toluene, 2-methylnaphthalene, phenanthrene) in 96-h tests. Resulting hydrocarbon toxicity data were fit to the Target Lipid Model to generate predictive models and determine species sensitivity. Toxicity of chemically enhanced water accommodated fractions of Alaskan North Slope crude oil (ANS-oil) was also investigated with shrimp using nominal loading, total petroleum hydrocarbons and biomimetic extraction (BE) as oil exposure metrics. Coral were more sensitive to dispersant than shrimp while similar sensitivity was observed for hydrocarbons. Study and literature findings indicate deepsea species exhibit acute sensitivities to dispersant, hydrocarbons and oil that are comparable to pelagic species. Results support use of passive sampling methods to quantify dissolved oil for interpreting oil toxicity tests and suggest models for predicting time-dependence of toxicity warrant re-evaluation.

The treatment of biodegradation in models of sub-surface oil spills: A review and sensitivity study.

Biodegradation is important for the fate of oil spilled in marine environments, yet parameterization of biodegradation varies across oil spill models, which usually apply constant first-order decay rates to multiple pseudo-components describing an oil. To understand the influence of model parameterization on the fate of subsurface oil droplets, we reviewed existing algorithms and rates and conducted a model sensitivity study. Droplets were simulated from a blowout at 2000 m depth and were either treated with sub-surface dispersant injection (2% dispersant to oil ratio) or untreated. The most important factor affecting oil fate was the size of the droplets, with biodegradation contributing substantially to the fate of droplets ≤0.5 mm. Oil types, which were similar, had limited influence on simulated oil fate. Model results suggest that knowledge of droplet sizes and improved estimation of pseudo-component biodegradation rates and lag times would enhance prediction of the fate and transport of subsurface oil.

The sensitivity of a deep-sea fish species (Anoplopoma fimbria) to oil-associated aromatic compounds, dispersant, and Alaskan North Slope crude oil.

A predominant concern following oil spills is toxicity to aquatic organisms. However, few data are available on effects in deep-sea cold water fishes. The present study had 3 major objectives. The first was to investigate the relative sensitivity of the deep-sea species Anoplopoma fimbria (sablefish) to acute effects of 3 aromatic compounds (toluene, 2-methylnaphthalene, and phenanthrene), dispersant alone, and chemically enhanced water accommodated fractions (CEWAFs) of Alaskan North Slope crude oil. The second was to determine the critical target lipid body burden (CTLBB) for sablefish by fitting aromatic hydrocarbon toxicity data to the target lipid model (TLM), which then allowed expression of CEWAF exposures in terms of dissolved oil toxic units. The final aim was to apply a passive sampling method that targets bioavailable, dissolved hydrocarbons as an alternative analytical technique for improved CEWAF exposure assessment. The results indicate that sablefish exhibit sensitivity to Corexit 9500 (96-h median lethal concentration [LC50] = 72.2 mg/L) within the range reported for other fish species. However, the acute CTLBB of 39.4 ± 2.1 µmol/goctanol lies at the lower end of the sensitivity range established for aquatic species. The utility of both toxic units and passive sampling measurements for describing observed toxicity of dispersed oil is discussed. The present study is novel in that a new test species is investigated to address the uncertainty regarding the sensitivity of deep-sea fishes, while also employing modeling and measurements to improve exposure characterization in oil toxicity tests. Environ Toxicol Chem 2018;37:2210-2221. © 2018 SETAC.

A re-evaluation of PETROTOX for predicting acute and chronic toxicity of petroleum substances.

The PETROTOX model was developed to perform aquatic hazard assessment of petroleum substances based on substance composition. The model relies on the hydrocarbon block method, which is widely used for conducting petroleum substance risk assessments providing further justification for evaluating model performance. Previous work described this model and provided a preliminary calibration and validation using acute toxicity data for limited petroleum substance. The objective of the present study was to re-evaluate PETROTOX using expanded data covering both acute and chronic toxicity endpoints on invertebrates, algae, and fish for a wider range of petroleum substances. The results indicated that recalibration of 2 model parameters was required, namely, the algal critical target lipid body burden and the log octanol-water partition coefficient (KOW ) limit, used to account for reduced bioavailability of hydrophobic constituents. Acute predictions from the updated model were compared with observed toxicity data and found to generally be within a factor of 3 for algae and invertebrates but overestimated fish toxicity. Chronic predictions were generally within a factor of 5 of empirical data. Furthermore, PETROTOX predicted acute and chronic hazard classifications that were consistent or conservative in 93 and 84% of comparisons, respectively. The PETROTOX model is considered suitable for the purpose of characterizing petroleum substance hazard in substance classification and risk assessments. Environ Toxicol Chem 2017;36:2245-2252. © 2017 SETAC.

Investigating the role of dissolved and droplet oil in aquatic toxicity using dispersed and passive dosing systems.

Characterization of the aquatic toxicity of oil is needed to support hazard assessment and inform spill response. Natural processes and mitigation strategies involving dispersant use can result in exposures to both dissolved and droplet oil that are not typically differentiated when oil exposures are characterized in toxicity tests. Thus, the impact of droplets on aquatic toxicity is largely uncharacterized. To improve the understanding of the role of droplets, acute toxicity tests with Daphnia magna and Americamysis bahia were performed with Endicott crude oil in low-energy mixing systems with and without Corexit 9500 dispersant. Exposures were also prepared by placing crude oil in silicone tubing and passively dosing test media to provide dissolved oil exposures without droplets. A framework is described for characterizing dissolved phase exposures using both mechanistic modeling and passive sampling measurements. The approach is then illustrated by application to data from the present study. Expression of toxicity in terms of toxic units calculated from modeled dissolved oil concentrations or passive sampling measurements showed similar dose responses between exposure systems and organisms, despite the gradient in droplet oil. These results indicate that droplets do not appreciably contribute to toxicity for the 2 species investigated and further support hazard evaluation of dispersed oil on the basis of dissolved exposure metrics. Environ Toxicol Chem 2017;36:1020-1028. © 2016 SETAC.

Guidance for improving comparability and relevance of oil toxicity tests.

The complex nature and limited aqueous solubility of petroleum substances pose challenges for consistently characterizing exposures in aquatic life hazard assessments. This paper reviews important considerations for the design, conduct and interpretation of laboratory toxicity tests with physically and chemically dispersed oils based on an understanding of the behavior and toxicity of the hydrocarbons that comprise these substances. Guiding principles are provided that emphasize the critical need to understand and, when possible, characterize dissolved hydrocarbon exposures that dictate observed toxicity in these tests. These principles provide a consistent framework for interpreting toxicity studies performed using different substances and test methods by allowing varying dissolved exposures to be expressed in terms of a common metric based on toxic units (TUs). The use of passive sampling methods is also advocated since such analyses provide an analytical surrogate for TUs. The proposed guidance is translated into a series of questions that can be used in evaluating existing data and in guiding design of future studies. Application of these questions to a number of recent publications indicates such considerations are often ignored, thus perpetuating the difficulty of interpreting and comparing results between studies and limiting data use in objective hazard assessment. Greater attention to these principles will increase the comparability and utility of oil toxicity data in decision-making.

Evaluating toxicity of heavy fuel oil fractions using complementary modeling and biomimetic extraction methods.

The toxicity of chemically dispersed heavy fuel oil (HFO) and 3 distillate fractions to rainbow trout (Oncorhynchus mykiss) embryos was evaluated using the PETROTOX model and a biomimetic extraction technique that involved passive sampling of oil-contaminated test media with solid-phase microextraction (SPME) fibers. Test solutions for toxicity testing were generated using a combination of dispersant and high-energy mixing. The resulting water accommodated fractions (WAF) provided complex exposure regimens that included both dissolved hydrocarbons and oil droplets. The toxicity of the various fractions differed by approximately 3 orders of magnitude when expressed on the basis of WAF dilution. Using detailed compositional data, the PETROTOX model predicted the speciation of hydrocarbons between dissolved and oil droplet phases and explained observed toxicity based on computed dissolved phase toxic units (TUs). A key finding from model calculations was that dissolved hydrocarbon exposures and associated TUs were a nonlinear function of WAF dilution, because dissolved hydrocarbons were largely controlled by the dissolution of oil droplets that were transferred in WAF dilutions. Hence, oil droplets served to "buffer" dissolved concentrations in WAF dilutions at loadings greater than 1 mg/L, resulting in higher dissolved concentrations and TUs than expected based on dilution. The TUs computed at each WAF dilution explained the observed toxicity among the HFO and fractions to within a factor of 3. Dissolved material measured by SPME showed a consistent relationship with model-predicted TUs, confirming the utility of this approach for providing an integrated measure of exposure to bioavailable hydrocarbons. These 2 approaches provide complementary tools for better defining bioavailability of complex petroleum substance.

PETRORISK: a risk assessment framework for petroleum substances.

PETRORISK is a modeling framework used to evaluate environmental risk of petroleum substances and human exposure through these routes due to emissions under typical use conditions as required by the European regulation for the Registration, Evaluation, Authorization and Restriction of Chemicals (REACH). Petroleum substances are often complex substances comprised of hundreds to thousands of individual hydrocarbons. The physicochemical, fate, and effects properties of the individual constituents within a petroleum substance can vary over several orders of magnitude, complicating risk assessment. PETRORISK combines the risk assessment strategies used on single chemicals with the hydrocarbon block approach to model complex substances. Blocks are usually defined by available analytical characterization data on substances that are expressed in terms of mass fractions for different structural chemical classes that are specified as a function of C number or boiling point range. The physicochemical and degradation properties of the blocks are determined by the properties of representative constituents in that block. Emissions and predicted exposure concentrations (PEC) are then modeled using mass-weighted individual representative constituents. Overall risk for various environmental compartments at the regional and local level is evaluated by comparing the PECs for individual representative constituents to corresponding predicted no-effect concentrations (PNEC) derived using the Target Lipid Model. Risks to human health are evaluated using the overall predicted human dose resulting from multimedia environmental exposure to a substance-specific derived no-effect level (DNEL). A case study is provided to illustrate how this modeling approach has been applied to assess the risks of kerosene manufacture and use as a fuel.

PETROTOX: an aquatic toxicity model for petroleum substances.

A spreadsheet model (PETROTOX) is described that predicts the aquatic toxicity of complex petroleum substances from petroleum substance composition. Substance composition is characterized by specifying mass fractions in constituent hydrocarbon blocks (HBs) based on available analytical information. The HBs are defined by their mass fractions within a defined carbon number range or boiling point interval. Physicochemical properties of the HBs are approximated by assigning representative hydrocarbons from a database of individual hydrocarbons with associated physicochemical properties. A three-phase fate model is used to simulate the distribution of each structure among the water-, air-, and oil-phase liquid in the laboratory test system. Toxicity is then computed based on the predicted aqueous concentrations and aquatic toxicity of each structure and the target lipid model. The toxicity of the complex substance is computed assuming additivity of the contribution of the individual assigned hydrocarbons. Model performance was evaluated by using direct comparisons with measured toxicity data for petroleum substances with sufficient analytical characterization to run the model. Indirect evaluations were made by comparing predicted toxicity distributions using analytical data on petroleum substances from different product categories with independent, empirical distributions of toxicity data available for the same categories. Predictions compared favorably with measured aquatic toxicity data across different petroleum substance categories. These findings demonstrate the utility of PETROTOX for assessing environmental hazards of petroleum substances given knowledge of substance composition.

Comparing laboratory and field measured bioaccumulation endpoints.

An approach for comparing laboratory and field measures of bioaccumulation is presented to facilitate the interpretation of different sources of bioaccumulation data. Differences in numerical scales and units are eliminated by converting the data to dimensionless fugacity (or concentration-normalized) ratios. The approach expresses bioaccumulation metrics in terms of the equilibrium status of the chemical, with respect to a reference phase. When the fugacity ratios of the bioaccumulation metrics are plotted, the degree of variability within and across metrics is easily visualized for a given chemical because their numerical scales are the same for all endpoints. Fugacity ratios greater than 1 indicate an increase in chemical thermodynamic activity in organisms with respect to a reference phase (e.g., biomagnification). Fugacity ratios less than 1 indicate a decrease in chemical thermodynamic activity in organisms with respect to a reference phase (e.g., biodilution). This method provides a holistic, weight-of-evidence approach for assessing the biomagnification potential of individual chemicals because bioconcentration factors, bioaccumulation factors, biota-sediment accumulation factors, biomagnification factors, biota-suspended solids accumulation factors, and trophic magnification factors can be included in the evaluation. The approach is illustrated using a total 2393 measured data points from 171 reports, for 15 nonionic organic chemicals that were selected based on data availability, a range of physicochemical partitioning properties, and biotransformation rates. Laboratory and field fugacity ratios derived from the various bioaccumulation metrics were generally consistent in categorizing substances with respect to either an increased or decreased thermodynamic status in biota, i.e., biomagnification or biodilution, respectively. The proposed comparative bioaccumulation endpoint assessment method could therefore be considered for decision making in a chemicals management context.

Hepatic gene expression in rainbow trout (Oncorhynchus mykiss) exposed to different hydrocarbon mixtures.

Traditional biomarkers for hydrocarbon exposure are not induced by all petroleum substances. The objective of this study was to determine if exposure to a crude oil and different refined oils would generate a common hydrocarbon-specific response in gene expression profiles that could be used as generic biomarkers of hydrocarbon exposure. Juvenile rainbow trout (Oncorhynchus mykiss) were exposed to the water accommodated fraction (WAF) of either kerosene, gas oil, heavy fuel oil, or crude oil for 96 h. Tissue was collected for RNA extraction and microarray analysis. Exposure to each WAF resulted in a different list of differentially regulated genes, with few genes in common across treatments. Exposure to crude oil WAF changed the expression of genes including cytochrome P4501A (CYP1A) and glutathione-S-transferase (GST) with known roles in detoxification pathways. These gene expression profiles were compared to others from previous experiments that used a diverse suite of toxicants. Clustering algorithms successfully identified gene expression profiles resulting from hydrocarbon exposure. These preliminary analyses highlight the difficulties of using single genes as diagnostic of petroleum hydrocarbon exposures. Further work is needed to determine if multivariate transcriptomic-based biomarkers may be a more effective tool than single gene studies for exposure monitoring of different oils.