A Data-Driven Approach to Estimating Occupational Inhalation Exposure Using Workplace Compliance Data

A growing list of chemicals are approved for production and use in the United States and elsewhere, and new approaches are needed to rapidly assess the potential exposure and health hazard posed by these substances. Here, we present a high-throughput, data-driven approach that will aid in estimating occupational exposure using a database of over 1.5 million observations of chemical concentrations in U.S. workplace air samples. We fit a Bayesian hierarchical model that uses industry type and the physicochemical properties of a substance to predict the distribution of workplace air concentrations. This model substantially outperforms a null model when predicting whether a substance will be detected in an air sample, and if so at what concentration, with 75.9% classification accuracy and a root-mean-square error (RMSE) of 1.00 log10 mg m-3 when applied to a held-out test set of substances. This modeling framework can be used to predict air concentration distributions for new substances, which we demonstrate by making predictions for 5587 new substance-by-workplace-type pairs reported in the US EPA's Toxic Substances Control Act (TSCA) Chemical Data Reporting (CDR) industrial use database. It also allows for improved consideration of occupational exposure within the context of high-throughput, risk-based chemical prioritization efforts.

A multi-scale approach for identification of potential pesticide use sites impacting vernal pool critical habitat in California

Spatially explicit ecological risk assessment (ERA) requires estimating the overlap between chemical and receptor distribution to evaluate the potential impacts of exposure on nontarget organisms. Pesticide use estimation at field level is prone to error due to inconsistencies between ground-reporting and geospatial data coverage; attempts to rectify these inconsistencies have been limited in approach and rarely scaled to multiple crop types. We built upon a previously developed Bayesian approach to combine multiple crop types for a probabilistic determination of field-crop assignments and to examine co-occurrence of critical vernal pool habitats and bifenthrin application within a 5-county area in California (Madera, Merced, Sacramento, San Joaquin, and Stanislaus counties). We incorporated a multi-scale repeated sampling approach with an area constraint to improve the delineation of field boundaries and better capture variability in crop assignments and rotation schemes. After comparing the accuracy of the spatial probabilistic approach to USDA Census of Agriculture crop acreage data, we found our approach allows more specificity in the combination of crop types represented by the potential application area and improves acreage estimates when compared to traditional deterministic approaches. In addition, our multi-scale sampling scheme improved estimates of bifenthrin acreage variability for co-occurrence analysis and allowed for estimates of crop rotations that were previously uncaptured. Our approach could be leveraged for more realistic, spatially resolved exposure and effects models both in and outside of California.

Estimating dermal contact soil exposure for amphibians

Chemical exposure estimation through the dermal route is an underemphasized area of ecological risk assessment for terrestrial animals. Currently, there are efforts to create exposure models to estimate doses from this pathway for use in ecological risk assessment. One significant limitation has been insufficient published data to characterize exposure and to support the selection and parameterization of appropriate models, particularly for amphibians in terrestrial habitats. Recent publications measuring pesticide doses to terrestrial-phase amphibians have begun to rectify this situation. We collated and summarized available measurements of terrestrial amphibian dermal exposure to pesticides from 11 studies in which researchers measured tissue concentrations associated with known pesticide experimental application rates. This data set included tissue concentrations in 11 amphibian species and 14 different pesticides. We then compared the results of two screening exposure models that differed based on surface area scaling approaches as a function of body weight (one based on birds as surrogates for amphibians and another amphibian-specific) to the measured tissue residue concentrations. We define a false-negative rate for each screening model as the proportion of amphibians for which the predicted concentration is less than the observed concentration (i.e., underestimate), contrary to the intent of screening models, which are intended to have a bias for higher exposure concentrations. The screening model that uses birds as surrogates did not have any instances where estimated expected avian doses were less than measured amphibian body burdens. When using the amphibian-specific exposure model that corrected for differences between avian and amphibian surface area, measured concentrations were greater than model estimates for 11.3% of the 1158 comparisons. The database of measured pesticide concentrations in terrestrial amphibians is provided for use in calculating bioconcentration factors and for future amphibian dermal exposure model development. Integr Environ Assess Manag 2023;19:9-16. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Differentiating metabolomic responses of amphibians to multiple stressors

One of the biggest challenges in ecological risk assessment is determining the impact of multiple stressors on individual organisms and populations in real world scenarios. Frequently, data derived from laboratory studies of single stressors are used to estimate risk parameters and do not adequately address scenarios where other stressors exist. Emerging 'omic technologies, notably metabolomics, provide an opportunity to address the uncertainties surrounding ecological risk assessment of multiple stressors. The objective of this study was to use metabolomic profiling to investigate the effect of multiple stressors on amphibian metamorphs. We exposed post-metamorphosis (180 days) southern leopard frogs (Lithobates sphenocephala) to the insecticide carbaryl (480 μg/L), predation stress, and a combined pesticide and predation stress treatment. Corticosterone analysis revealed mild support for an induction in response to predation stress alone but strongly suggests that carbaryl exposure, alone or in combination with predation cues, can significantly elevate this known biomarker in amphibians. Metabolomics analysis accurately classed, based on relative nearness, carbaryl and predation induced changes in the hepatic metabolome and biochemical fluxes appear to be associated with a similar biological response. Support vector machine analysis with recursive feature elimination of the acquired metabolomic spectra demonstrated 85-96% classification accuracy among control and all treatment groups when using the top 75 ranked retention time bins. Biochemical fluxes observed in the groups exposed to carbaryl, predation, and the combined treatment include amino acids, sugar derivatives, and purine nucleotides. Ultimately, this methodology could be used to interpret short-term toxicity assays and the presence of environmental stressors to overall metabolomic effects in non-target organisms.

Integrative exposomic, transcriptomic, epigenomic analyses of human placental samples links understudied chemicals to preeclampsia

BACKGROUND

Environmental health research has recently undergone a dramatic shift, with ongoing technological advancements allowing for broader coverage of exposure and molecular biology signatures. Approaches to integrate such measures are still needed to increase understanding between systems-level exposure and biology.

OBJECTIVES

We address this gap by evaluating placental tissues to identify novel chemical-biological interactions associated with preeclampsia. This study tests the hypothesis that understudied chemicals are present in the human placenta and associated with preeclampsia-relevant disruptions, including overall case status (preeclamptic vs. normotensive patients) and underlying transcriptomic/epigenomic signatures.

METHODS

A non-targeted analysis based on high-resolution mass spectrometry was used to analyze placental tissues from a cohort of 35 patients with preeclampsia (n = 18) and normotensive (n = 17) pregnancies. Molecular feature data were prioritized for confirmation based on association with preeclampsia case status and confidence of chemical identification. All molecular features were evaluated for relationships to mRNA, microRNA, and CpG methylation (i.e., multi-omic) signature alterations involved in preeclampsia.

RESULTS

A total of 183 molecular features were identified with significantly differentiated abundance in placental extracts of preeclamptic patients; these features clustered into distinct chemical groupings using unsupervised methods. Of these features, 53 were identified (mapping to 40 distinct chemicals) using chemical standards, fragmentation spectra, and chemical metadata. In general, human metabolites had the largest feature intensities and strongest associations with preeclampsia-relevant multi-omic changes. Exogenous drugs were second most abundant and had fewer associations with multi-omic changes. Other exogenous chemicals (non-drugs) were least abundant and had the fewest associations with multi-omic changes.

CONCLUSIONS

These global data trends suggest that human metabolites are heavily intertwined with biological processes involved in preeclampsia etiology, while exogenous chemicals may still impact select transcriptomic/epigenomic processes. This study serves as a demonstration of merging systems exposures with systems biology to better understand chemical-disease relationships.

Investigating Vernal Pool Fairy Shrimp Exposure to Organophosphate Pesticides: Implications for Population-Level Risk Assessment

Vernal pool fairy shrimp, Branchinecta lynchi, is a freshwater crustacean endemic to California and Oregon, including California's Central Valley. B. lynchi is listed as a Federally Threatened species under the US Endangered Species Act, and as a vulnerable species on the IUCN Red List. Threats that may negatively impact vernal pool fairy shrimp populations include pesticide applications to agricultural land use (e.g., agrochemicals such as organophosphate pesticides) and climate changes that impact vernal pool hydrology. Pop-GUIDE (Population model Guidance, Use, Interpretation, and Development for Ecological risk assessment) is a comprehensive tool that facilitates development and implementation of population models for ecological risk assessment and can be used to document the model derivation process. We employed Pop-GUIDE to document and facilitate the development of a population model for investigating impacts of organophosphate pesticides on vernal pool fairy shrimp populations in California's Central Valley. The resulting model could be applied in combination with field assessment and laboratory-based chemical analysis to link effects from pesticide exposure to adverse outcomes in populations across their range. B. lynchi has a unique intra-annual life cycle that is largely dependent upon environmental conditions. Future deployment of this population model should include complex scenarios consisting of multiple stressors, whereby the model is used to examine scenarios that combine chemical stress resulting from exposure to pesticides and climate changes.

Satellite-derived cyanobacteria frequency and magnitude in headwaters & near-dam reservoir surface waters of the Southern U.S

Reservoirs are dominant features of the modern hydrologic landscape and provide vital services. However, the unique morphology of reservoirs can create suitable conditions for excessive algae growth and associated cyanobacteria blooms in shallow in-flow reservoir locations by providing warm water environments with relatively high nutrient inputs, deposition, and nutrient storage. Cyanobacteria harmful algal blooms (cyanoHAB) are costly water management issues and bloom recurrence is associated with economic costs and negative impacts to human, animal, and environmental health. As cyanoHAB occurrence varies substantially within different regions of a water body, understanding in-lake cyanoHAB spatial dynamics is essential to guide reservoir monitoring and mitigate potential public exposure to cyanotoxins. Cloud-based computational processing power and high temporal frequency of satellites enables advanced pixel-based spatial analysis of cyanoHAB frequency and quantitative assessment of reservoir headwater in-flows compared to near-dam surface waters of individual reservoirs. Additionally, extensive spatial coverage of satellite imagery allows for evaluation of spatial trends across many dozens of reservoir sites. Surface water cyanobacteria concentrations for sixty reservoirs in the southern U.S. were estimated using 300 m resolution European Space Agency (ESA) Ocean and Land Colour Instrument (OLCI) satellite sensor for a five year period (May 2016-April 2021). Of the reservoirs studied, spatial analysis of OLCI data revealed 98% had more frequent cyanoHAB occurrence above the concentration of >100,000 cells/mL in headwaters compared to near-dam surface waters (P < 0.001). Headwaters exhibited greater seasonal variability with more frequent and higher magnitude cyanoHABs occurring mid-summer to fall. Examination of reservoirs identified extremely high concentration cyanobacteria events (>1,000,000 cells/mL) occurring in 70% of headwater locations while only 30% of near-dam locations exceeded this threshold. Wilcoxon signed-rank tests of cyanoHAB magnitudes using paired-observations (dates with observations in both a reservoir's headwater and near-dam locations) confirmed significantly higher concentrations in headwater versus near-dam locations (p < 0.001).

Using metabolomic profiling to inform use of surrogate species in ecological risk assessment practices

The U.S. EPA frequently uses avian or fish toxicity data to set protective standards for amphibians in ecological risk assessments. However, this approach does not always adequately represent aquatic-dwelling and terrestrial-phase amphibian exposure data. For instance, it is accepted that early life stage tests for fish are typically sensitive enough to protect larval amphibians, however, metamorphosis from tadpole to a terrestrial-phase adult relies on endocrine cues that are less prevalent in fish but essential for amphibian life stage transitions. These differences suggest that more robust approaches are needed to adequately elucidate the impacts of pesticide exposure in amphibians across critical life stages. Therefore, in the current study, methodology is presented that can be applied to link the perturbations in the metabolomic response of larval zebrafish (Danio rerio), a surrogate species frequently used in ecotoxicological studies, to those of African clawed frog (Xenopus laevis) tadpoles following exposure to three high-use pesticides, bifenthrin, chlorothalonil, or trifluralin. Generally, D. rerio exhibited greater metabolic perturbations in both number and magnitude across the pesticide exposures as opposed to X. laevis. This suggests that screening ecological risk assessment surrogate toxicity data would sufficiently protect amphibians at the single life stage studied but care needs to be taken to understand the suite of metabolic requirements of each developing species. Ultimately, methodology presented, and data gathered herein will help inform the applicability of metabolomic profiling in establishing the risk pesticide exposure poses to amphibians and potentially other non-target species.

Induced Hepatic Glutathione and Metabolomic Alterations Following Mixed Pesticide and Fertilizer Exposures in Juvenile Leopard Frogs (Lithobates sphenocephala)

The increasing use of agrochemicals, alone and in combination, has been implicated as a potential causative factor in the decline of amphibians worldwide. Fertilizers and pesticides are frequently combined into single-use tank mixtures for agricultural applications to decrease costs while meeting the food demands of a growing human population. Limited data are available on the effects of increased nitrogen levels in nontarget species, such as amphibians, and therefore investigating alterations in the nitrogen cycle and its impacts on amphibians needs to be considered in best management practices going forward. The objective of the present study was to elucidate the impact of fertilizer (urea) and herbicide (atrazine and/or alachlor) tank mixtures on the hepatic metabolome of juvenile leopard frogs as well as to investigate alterations in oxidative stress by relating these changes to glutathione (GSH) levels. Herbicide exposure only moderately increased this parameter in amphibians, however, urea alone and in combination with either atrazine or alachlor statistically elevated GSH levels. Interestingly, urea also inhibited pesticide uptake: calculated bioconcentration factors were greatly decreased for atrazine and alachlor when urea was present in the exposure mixture. Metabolomic profiling identified fluxes in hepatic metabolites that are involved in GSH and carbohydrate metabolic processes as well as altered intermediates in the urea cycle. Ultimately, understanding the biological impacts of nitrogenous fertilizers alone and in combination with pesticide exposure will inform best management practices to conserve declining amphibian populations worldwide. Environ Toxicol Chem 2022;41:122-133. © 2021 SETAC.

Inferring pesticide toxicity to honey bees from a field-based feeding study using a colony model and Bayesian inference

Honey bees are crucial pollinators for agricultural crops but are threatened by a multitude of stressors including exposure to pesticides. Linking our understanding of how pesticides affect individual bees to colony-level responses is challenging because colonies show emergent properties based on complex internal processes and interactions among individual bees. Agent-based models that simulate honey bee colony dynamics may be a tool for scaling between individual and colony effects of a pesticide. The U.S. Environmental Protection Agency (USEPA) and U.S. Department of Agriculture (USDA) are developing the VarroaPop + Pesticide model, which simulates the dynamics of honey bee colonies and how they respond to multiple stressors, including weather, Varroa mites, and pesticides. To evaluate this model, we used Approximate Bayesian Computation to fit field data from an empirical study where honey bee colonies were fed the insecticide clothianidin. This allowed us to reproduce colony feeding study data by simulating colony demography and mortality from ingestion of contaminated food. We found that VarroaPop + Pesticide was able to fit general trends in colony population size and structure and reproduce colony declines from increasing clothianidin exposure. The model underestimated adverse effects at low exposure (36 µg/kg), however, and overestimated recovery at the highest exposure level (140 µg/kg), for the adult and pupa endpoints, suggesting that mechanisms besides oral toxicity-induced mortality may have played a role in colony declines. The VarroaPop + Pesticide model estimates an adult oral LD₅₀ of 18.9 ng/bee (95% CI 10.1-32.6) based on the simulated feeding study data, which falls just above the 95% confidence intervals of values observed in laboratory toxicology studies on individual bees. Overall, our results demonstrate a novel method for analyzing colony-level data on pesticide effects on bees and making inferences on pesticide toxicity to individual bees.

Route of exposure influences pesticide body burden and the hepatic metabolome in post-metamorphic leopard frogs

Pesticides are being applied at a greater extent than in the past. Once pesticides enter the ecosystem, many environmental factors can influence their residence time. These interactions can result in processes such as translocation, environmental degradation, and metabolic activation facilitating exposure to target and non-target species. Most anurans start off their life cycle in aquatic environments and then transition into terrestrial habitats. Their time in the aquatic environment is generally short; however, many important developmental stages occur during this tenure. Post-metamorphosis, most species spend many years on land but migrate back to the aquatic environment for breeding. Due to the importance of both the aquatic and terrestrial environments to the life stages of amphibians, we investigated how the route of exposure (i.e., uptake from contaminated soils vs. uptake from contaminated surface water) influences pesticide bioavailability and body burden for four pesticides (bifenthrin (BIF), chlorpyrifos (CPF), glyphosate (GLY), and trifloxystrobin (TFS)) as well as the impact on the hepatic metabolome of adult leopard frogs (Gosner stage 46 with 60-90 days post-metamorphosis). Body burden concentrations for amphibians exposed in water were significantly higher (ANOVA p < 0.0001) compared to amphibians exposed to contaminated soil across all pesticides studied. Out of 80 metabolites that were putatively identified, the majority expressed a higher abundance in amphibians that were exposed in pesticide contaminated water compared to soil. Ultimately, this research will help fill regulatory data gaps, aid in the creation of more accurate amphibian dermal uptake models and inform continued ecological risk assessment efforts.

Honey Bees and Neonicotinoid-Treated Corn Seed: Contamination, Exposure, and Effects

Most corn (Zea mays) seeds planted in the United States in recent years are coated with a seed treatment containing neonicotinoid insecticides. Abrasion of the seed coating generates insecticide-laden planter dust that disperses through the landscape during corn planting and has resulted in many "bee-kill" incidents in North America and Europe. We investigated the linkage between corn planting and honey bee colony success in a region dominated by corn agriculture. Over 3 yr we consistently observed an increased presence of corn seed treatment insecticides in bee-collected pollen and elevated worker bee mortality during corn planting. Residues of seed treatment neonicotinoids, clothianidin and thiamethoxam, detected in pollen positively correlated with cornfield area surrounding the apiaries. Elevated worker mortality was also observed in experimental colonies fed field-collected pollen containing known concentrations of corn seed treatment insecticides. We monitored colony growth throughout the subsequent year in 2015 and found that colonies exposed to higher insecticide concentrations exhibited slower population growth during the month of corn planting but demonstrated more rapid growth in the month following, though this difference may be related to forage availability. Exposure to seed treatment neonicotinoids during corn planting has clear short-term detrimental effects on honey bee colonies and may affect the viability of beekeeping operations that are dependent on maximizing colony size in the springtime. Environ Toxicol Chem 2021;40:1212-1221. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Comparison and Evaluation of Gridded Precipitation Datasets in a Kansas Agricultural Watershed Using SWAT

Gridded precipitation datasets are becoming a convenient substitute for gauge measurements in hydrological modeling; however, these data have not been fully evaluated across a range of conditions. We compared four gridded datasets (Daily Surface Weather and Climatological Summaries [DAYMET], North American Land Data Assimilation System [NLDAS], Global Land Data Assimilation System [GLDAS], and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) as precipitation data sources and evaluated how they affected hydrologic model performance when compared with a gauged dataset, Global Historical Climatology Network-Daily (GHCN-D). Analyses were performed for the Delaware Watershed at Perry Lake in eastern Kansas. Precipitation indices for DAYMET and PRISM precipitation closely matched GHCN-D, whereas NLDAS and GLDAS showed weaker correlations. We also used these precipitation data as input to the Soil and Water Assessment Tool (SWAT) model that confirmed similar trends in streamflow simulation. For stations with complete data, GHCN-D based SWAT-simulated streamflow variability better than gridded precipitation data. During low flow periods we found PRISM performed better, whereas both DAYMET and NLDAS performed better in high flow years. Our results demonstrate that combining gridded precipitation sources with gauge-based measurements can improve hydrologic model performance, especially for extreme events.

A sensitivity analysis of pesticide concentrations in California Central Valley vernal pools

Vernal pools are ephemeral wetlands that provide critical habitat to many listed species. Pesticide fate in vernal pools is poorly understood because of uncertainties in the amount of pesticide entering these ecosystems and their bioavailability throughout cycles of wet and dry periods. The Pesticide Water Calculator (PWC), a model used for the regulation of pesticides in the US, was used to predict surface water and sediment pore water pesticide concentrations in vernal pool habitats. The PWC model (version 1.59) was implemented with deterministic and probabilistic approaches and parameterized for three agricultural vernal pool watersheds located in the San Joaquin River basin in the Central Valley of California. Exposure concentrations for chlorpyrifos, diazinon and malathion were simulated. The deterministic approach used default values and professional judgment to calculate point values of estimated concentrations. In the probabilistic approach, Monte Carlo (MC) simulations were conducted across the full input parameter space with a sensitivity analysis that quantified the parameter contribution to model prediction uncertainty. Partial correlation coefficients were used as the primary sensitivity metric for analyzing model outputs. Conditioned daily sensitivity analysis indicates curve number (CN) and the universal soil loss equation (USLE) parameters as the most important environmental parameters. Therefore, exposure estimation can be improved efficiently by focusing parameterization efforts on these driving processes, and agricultural pesticide inputs in these critical habitats can be reduced by best management practices focused on runoff and sediment reductions.

A unified approach for protecting listed species and ecosystem services in isolated wetlands using community-level protection goals

The protection of listed species through the Ecological Risk Assessment (ERA) process is encumbered by the number and diversity of species that need protection and the limited data available to inform assessments. Ecological communities within isolated ecosystems often contain a number of biologically diverse endemic, endangered, and threatened species, as well as providing numerous ecosystem services (ES). We propose an approach that develops community-level protection goals using isolated wetlands that includes both listed species and Service Providing Units (SPUs) that drive ES for ecological risk assessments (ERAs). Community-level protection goals are achieved by developing a protection community and weighing lines of evidence to determine a set of focal species within that community upon which to base the assessment. Lines of evidence include chemical mechanism of action, likely routes of exposure, and taxa susceptibility, as well as relationships among species, and other ecological factors. We demonstrate the process using case studies of chlorpyrifos in California vernal pools and coal ash effluent in Carolina bays. In the California vernal pool case study, listed species were the primary SPUs for the ES provided by the critical habitat. The weight of evidence demonstrated the honey bee as the focal species for the terrestrial environment and the vernal pool fairy shrimp as the focal species for the aquatic environment. The protection community within the Carolina bay case study was more taxonomically diverse than vernal pools for both listed species and SPUs, with amphibians identified as the focal species for which to target mitigation goals and hazard levels. The approach presented here will reduce the time and resource investment required for assessment of risk to listed species and adds an ES perspective to demonstrate value of assessments beyond listed species concerns.

Endogenous and exogenous biomarker analysis in terrestrial phase amphibians (Lithobates sphenocephala) following dermal exposure to pesticide mixtures

Pesticide mixtures are frequently co-applied throughout an agricultural growing season to maximize crop yield. Therefore, non-target ecological species (e.g., amphibians) may be exposed to several pesticides at any given time on these agricultural landscapes. The objectives of this study were to quantify body burdens in terrestrial phase amphibians and translate perturbed metabolites to their corresponding biochemical pathways affected by exposure to pesticides as both singlets and in combination. Southern leopard frogs (Lithobates sphenocephala) were exposed either at maximum or 1/10^th maximum application rate to single, double, or triple pesticide mixtures of bifenthrin (insecticide), metolachlor (herbicide), and triadimefon (fungicide). Tissue concentrations demonstrate both facilitated and competitive uptake of pesticides when in mixtures. Metabolomic profiling of amphibian livers identified metabolites of interest for both application rates, however; magnitude of changes varied for the two exposure rates. Exposure to lower concentrations demonstrated down regulation in amino acids, potentially due to their being utilized for glutathione metabolism and/or increased energy demands. Amphibians exposed to the maximum application rate resulted in up regulation of amino acids and other key metabolites likely due to depleted energy resources. Coupling endogenous and exogenous biomarkers of pesticide exposure can be utilized to form vital links in an ecological risk assessment by relating internal dose to pathophysiological outcomes in non-target species.

Analysis of pesticides in surface water, stemflow, and throughfall in an agricultural area in South Georgia, USA

To study spray drift contributions to non-targeted habitats, pesticide concentrations in stemflow (water flowing down the trunk of a tree during a rain event), throughfall (water from tree canopy only), and surface water in an agriculturally impacted wetland area near Tifton, Georgia, USA were measured (2015-2016). Agricultural fields and sampling locations were on the University of Georgia's Gibbs Research Farm, Tifton, GA. Samples were screened for more than 160 pesticides, and cumulatively, 32 different pesticides were detected across matrices. Data indicate that herbicides and fungicides were present in all types of environmental samples analyzed while insecticides were only detected in surface water samples. The highest pesticide concentration observed was 10.50 μg/L of metolachlor in an August 2015 surface water sample. Metolachlor, tebuconazole, and fipronil were the most frequently detected herbicide, fungicide, and insecticide, respectively, regardless of sample origin. The most frequently detected pesticide in surface water and stemflow samples was metolachlor (0.09-10.5 μg/L), however, the most commonly detected pesticide in throughfall samples was biphenyl (0.02-0.07 μg/L). These data help determine the importance of indirect chemical exposures to non-targeted habitats by assessing inputs from stemflow and throughfall into surface waters.

Effect of hydration status on pesticide uptake in anurans following exposure to contaminated soils

In this study, the impact of hydration status on dermal uptake of pesticides in two species of amphibians is examined. Absorption of pesticides in anurans occurs primarily through a highly vascularized dermal seat patch; however, pesticides can also enter through the superficial dermis following exposure. Despite the growing body of literature on dermal exposure in amphibians, little is known on how hydration status influences uptake. Thus, the objective of this study was to investigate the influence of hydration status on absorption of pesticides (atrazine, triadimefon, metolachlor, chlorothalonil, and imidacloprid) in southern leopard frogs (Lithobates sphenocephala) and Fowler's toads (Anaxyrus fowleri). Amphibian treatments included dehydration periods of 0, 2, 4, 6, 8, or 10 h prior to exposure to pesticide-contaminated soils for 8 h. Following exposure, soil and whole-body homogenates were extracted and analyzed by LC-MS/MS. Dehydration time was then regressed against post-exposure concentrations to infer the impact of dehydration on dermal pesticide uptake. Increased dehydration time resulted in significantly lowered pesticide concentrations in both species (F_{6, 293} = 67.66, p = 0.007) for the five pesticides studied. This phenomenon could be due to an energy and/or dilution effect.

Influence of exposure to pesticide mixtures on the metabolomic profile in post-metamorphic green frogs (Lithobates clamitans)

Pesticide use in agricultural areas requires the application of numerous chemicals to control target organisms, leaving non-target organisms at risk. The present study evaluates the hepatic metabolomic profile of one group of non-target organisms, amphibians, after exposure to a single pesticide and pesticide mixtures. Five common-use pesticide active ingredients were used in this study, three herbicides (atrazine, metolachlor and 2,4-d), one insecticide (malathion) and one fungicide (propiconazole). Juvenile green frogs (Lithobates clamitans) were reared for 60-90days post-metamorphosis then exposed to a single pesticide or a combination of pesticides at the labeled application rate on soil. Amphibian livers were excised for metabolomic analysis and pesticides were quantified for whole body homogenates. Based on the current study, metabolomic profiling of livers support both individual and interactive effects where pesticide exposures altered biochemical processes, potentially indicating a different response between active ingredients in pesticide mixtures, among these non-target species. Amphibian metabolomic response is likely dependent on the pesticides present in each mixture and their ability to perturb biochemical networks, thereby confounding efforts with risk assessment.

Sensitivity analyses for simulating pesticide impacts on honey bee colonies

We employ Monte Carlo simulation and sensitivity analysis techniques to describe the population dynamics of pesticide exposure to a honey bee colony using the VarroaPop+Pesticide model. Simulations are performed of hive population trajectories with and without pesticide exposure to determine the effects of weather, queen strength, foraging activity, colony resources, and Varroa populations on colony growth and survival. The daily resolution of the model allows us to conditionally identify sensitivity metrics. Simulations indicate queen strength and forager lifespan are consistent, critical inputs for colony dynamics in both the control and exposed conditions. Adult contact toxicity, application rate and nectar load become critical parameters for colony dynamics within exposed simulations. Daily sensitivity analysis also reveals that the relative importance of these parameters fluctuates throughout the simulation period according to the status of other inputs.

Using in vitro derived enzymatic reaction rates of metabolism to inform pesticide body burdens in amphibians

Understanding how pesticide exposure to non-target species influences toxicity is necessary to accurately assess the ecological risks these compounds pose. To assess the potential metabolic activation of broad use pesticides in amphibians, in vitro and in vivo metabolic rate constants were derived from toad (Anaxyrus terrestris) livers in experiments measuring the depletion of atrazine (ATZ), triadimefon (TDN), and fipronil (FIP) as well as formation of their metabolites. To determine the predictability of these in vitro derived rate constants, Fowler's toads (Anaxyrus fowleri) were exposed to soil contaminated with each of the pesticides at maximum application rate. Desethyl atrazine (DEA) and deisopropyl atrazine (DIA), both metabolites of ATZ, exhibited similar velocities (V_max) while the K_M constant for DIA was two times higher than DEA. TDN was metabolized into two diastereomers of triadimenol (TDL A and TDL B), where TDL B had a V_max around two times higher than TDL A. The metabolite fipronil sulfone's V_max and K_M were 150 pmol min^-1 mg^-1 and 29 μM, respectively. While intrinsic clearance rates for the pesticides ranged from 0.54 to 38.31 mL min^-1 kg^-1. Thus, gaining knowledge on differences in metabolism of pesticides within amphibians is important in estimating risk to these non-target species since the inherent toxicity of metabolites can differ from the parent compound.

Spatiotemporal sensitivity analysis of vertical transport of pesticides in soil

Environmental fate and transport processes are influenced by many factors. Simulation models that mimic these processes often have complex implementations, which can lead to over-parameterization. Sensitivity analyses are subsequently used to identify critical parameters whose uncertainties can be further reduced or better described and prediction variability minimized. In this study, a variance decomposition based global sensitivity analysis technique (Sobol' method) is conducted based on estimated concentrations in vertical soil compartments using the Pesticide Root Zone Model (PRZM). Daily simulations are performed that explore the input parameter space. Estimated concentrations are compared to data collected over the course of a growing season from an experimental site in Georgia. Our results suggest that model sensitivity is conditional and should be examined at appropriate spatial and temporal resolution to avoid omitting important parameters. This approach can yield a better understanding about the interplay between sensitivity/uncertainty and model dynamics in non-monotonic, non-linear systems.

Biomarker analysis of American toad (Anaxyrus americanus) and grey tree frog (Hyla versicolor) tadpoles following exposure to atrazine

The objective of the current study was to use a biomarker-based approach to investigate the influence of atrazine exposure on American toad (Anaxyrus americanus) and grey tree frog (Hyla versicolor) tadpoles. Atrazine is one of the most frequently detected herbicides in environmental matrices throughout the United States. In surface waters, it has been found at concentrations from 0.04-2859μg/L and thus presents a likely exposure scenario for non-target species such as amphibians. Studies have examined the effect of atrazine on the metamorphic parameters of amphibians, however, the data are often contradictory. Gosner stage 22-24 tadpoles were exposed to 0 (control), 10, 50, 250 or 1250μg/L of atrazine for 48h. Endogenous polar metabolites were extracted and analyzed using gas chromatography coupled with mass spectrometry. Statistical analyses of the acquired spectra with machine learning classification models demonstrated identifiable changes in the metabolomic profiles between exposed and control tadpoles. Support vector machine models with recursive feature elimination created a more efficient, non-parametric data analysis and increased interpretability of metabolomic profiles. Biochemical fluxes observed in the exposed groups of both A. americanus and H. versicolor displayed perturbations in a number of classes of biological macromolecules including fatty acids, amino acids, purine nucleosides, pyrimidines, and mono- and di-saccharides. Metabolomic pathway analyses are consistent with findings of other studies demonstrating disruption of amino acid and energy metabolism from atrazine exposure to non-target species.

Soil organic matter content effects on dermal pesticide bioconcentration in American toads (Bufo americanus)

Pesticides have been implicated as a major factor in global amphibian declines and may pose great risk to terrestrial phase amphibians moving to and from breeding ponds on agricultural landscapes. Dermal uptake from soil is known to occur in amphibians, but predicting pesticide availability and bioconcentration across soil types is not well understood. The present study was designed to compare uptake of 5 current-use pesticides (imidacloprid, atrazine, triadimefon, fipronil, and pendimethalin) in American toads (Bufo americanus) from exposure on soils with significant organic matter content differences (14.1% = high organic matter and 3.1% = low organic matter). We placed toads on high- or low-organic matter soil after applying individual current-use pesticides on the soil surface for an 8-h exposure duration. Whole body tissue homogenates and soils were extracted and analyzed using liquid chromatography-mass spectrometry to determine pesticide tissue and soil concentration, as well as bioconcentration factor in toads. Tissue concentrations were greater on the low-organic matter soil than the high-organic matter soil across all pesticides (average ± standard error; 1.23 ± 0.35 ppm and 0.78 ± 0.23 ppm, respectively), and bioconcentration was significantly higher for toads on the low-organic matter soil (analysis of covariance p = 0.002). Soil organic matter is known to play a significant role in the mobility of pesticides and bioavailability to living organisms. Agricultural soils typically have relatively lower organic matter content and serve as a functional habitat for amphibians. The potential for pesticide accumulation in amphibians moving throughout agricultural landscapes may be greater and should be considered in conservation and policy efforts. Environ Toxicol Chem 2016;35:2734-2741. © 2016 SETAC.

Rainfall-induced release of microbes from manure: model development, parameter estimation, and uncertainty evaluation on small plots

A series of simulated rainfall-runoff experiments with applications of different manure types (cattle solid pats, poultry dry litter, swine slurry) was conducted across four seasons on a field containing 36 plots (0.75 × 2 m each), resulting in 144 rainfall-runoff events. Simulating time-varying release of Escherichia coli, enterococci, and fecal coliforms from manures applied at typical agronomic rates evaluated the efficacy of the Bradford-Schijven model modified by adding terms for release efficiency and transportation loss. Two complementary, parallel approaches were used to calibrate the model and estimate microbial release parameters. The first was a four-step sequential procedure using the inverse model PEST, which provides appropriate initial parameter values. The second utilized a PEST/bootstrap procedure to estimate average parameters across plots, manure age, and microbe, and to provide parameter distributions. The experiment determined that manure age, microbe, and season had no clear relationship to the release curve. Cattle solid pats released microbes at a different, slower rate than did poultry dry litter or swine slurry, which had very similar release patterns. These findings were consistent with other published results for both bench- and field-scale, suggesting the modified Bradford-Schijven model can be applied to microbial release from manure.

Pesticide Uptake Across the Amphibian Dermis Through Soil and Overspray Exposures

For terrestrial amphibians, accumulation of pesticides through dermal contact is a primary route of exposure in agricultural landscapes and may be contributing to widespread amphibian declines. To show pesticide transfer across the amphibian dermis at permitted label application rates, our study was designed to measure pesticide body burdens after two simulated exposure scenarios. We compared direct exposures, where amphibians were present when spraying occurred, to indirect exposures, where amphibians were exposed to soils after pesticide application. During summer 2012, we reared barking (Hyla gratiosa) and green treefrogs (H. cinerea) through 60-90 days post-metamorphosis at a United States Environmental Protection Agency research laboratory. We tested exposure for 8 h to five pesticide active ingredients (imidacloprid, atrazine, triadimefon, fipronil, or pendimethalin) in glass aquaria lined with soil in the laboratory. We quantified total pesticide body burden and soil concentrations using liquid chromatography-mass spectrometry. All individuals in both treatments had measurable body burdens at the end of the study. A randomized block design analysis of variance (n = 18) showed that body burdens (p = 0.03) and bioconcentration factors (BCFs) (p = 0.01) were significantly greater in the direct overspray treatment relative to the indirect soil spray treatment for both species and tested pesticides. BCFs ranged from 0.1 to 1.16 and from 0.013 to 0.78 in the direct and indirect treatments, respectively. Our study shows dermal uptake for multiple pesticides from both direct spray and indirect soil exposures and provides empirical support for the degree to which terrestrial phase amphibians have higher body burdens after overspray pesticide exposure.

Estimating terrestrial amphibian pesticide body burden through dermal exposure

Dermal exposure presents a potentially significant but understudied route for pesticide uptake in terrestrial amphibians. Our study measured dermal uptake of pesticides of varying hydrophobicity (logKow) in frogs. Amphibians were indirectly exposed to one of five pesticide active ingredients through contact with contaminated soil: imidacloprid (logKow = 0.57), atrazine (logKow = 2.5), triadimefon (logKow = 3.0), fipronil (logKow = 4.11) or pendimethalin (logKow = 5.18). All amphibians had measurable body burdens at the end of the exposure in concentrations ranging from 0.019 to 14.562 μg/g across the pesticides tested. Atrazine produced the greatest body burdens and bioconcentration factors, but fipronil was more permeable to amphibian skin when application rate was considered. Soil partition coefficient and water solubility were much better predictors of pesticide body burden, bioconcentration factor, and skin permeability than logKow. Dermal uptake data can be used to improve risk estimates of pesticide exposure among amphibians as non-target organisms.

Brazilian free-tailed bats as insect pest regulators in transgenic and conventional cotton crops

During the past 12000 years agricultural systems have transitioned from natural habitats to conventional agricultural regions and recently to large areas of genetically engineered (GE) croplands. This GE revolution occurred for cotton in a span of slightly more than a decade during which a switch occurred in major cotton production areas from growing 100% conventional cotton to an environment in which 95% transgenics are grown. Ecological interactions between GE targeted insects and other insectivorous insects have been investigated. However, the relationships between ecological functions (such as herbivory and ecosystem transport) and agronomic benefits of avian or mammalian insectivores in the transgenic environment generally remain unclear, although the importance of some agricultural pest management services provided by insectivorous species such as the Brazilian free-tailed bat, Tadarida brasiliensis, have been recognized. We developed a dynamic model to predict regional-scale ecological functions in agricultural food webs by using the indicators of insect pest herbivory measured by cotton boll damage and insect emigration from cotton. In the south-central Texas Winter Garden agricultural region we find that the process of insectivory by bats has a considerable impact on both the ecology and valuation of harvest in Bacillus thuringiensis (Bt) transgenic and nontransgenic cotton crops. Predation on agricultural pests by insectivorous bats may enhance the economic value of agricultural systems by reducing the frequency of required spraying and delaying the ultimate need for new pesticides. In the Winter Garden region, the presence of large numbers of insectivorous bats yields a regional summer dispersion of adult pest insects from Bt cotton that is considerably reduced from the moth emigration when bats are absent in either transgenic or non-transgenic crops. This regional decrease of pest numbers impacts insect herbivory on a transcontinental scale. With a few exceptions, we find that the agronomics of both Bt and conventional cotton production is more profitable when large numbers of insectivorous bats are present.