Simulated developmental and reproductive impacts on amphibian populations and implications for assessing long-term effects

Exposure to Low Concentrations of AMPA Influences Morphology and Decreases Survival During Larval Development in a Widespread Amphibian Species

Glyphosate's primary metabolite, AMPA (aminomethylphosphonic acid), is one of the most widely detected anthropogenic substance in surface waters worldwide. However, ecotoxicological studies on the potential effects of this metabolite at environmental concentrations on wildlife are scarce. Yet, due to its chemical properties, AMPA is likely to affect non-target species. In this study, we investigated sublethal effects of environmental concentrations of AMPA on the larval development of a widespread amphibian species, the spined toad Bufo spinosus. We performed a factorial experiment to study the effect of concentration and the timing of exposure (during embryonic development, larval development or both) to AMPA on the morphology, rate of development and survival of tadpoles. AMPA and timing of exposure interactively affected tadpole size (individuals exposed to AMPA after hatching were transitorily smaller, while individuals exposed to AMPA before hatching were longer), but not duration of development. Most of these effects were linked to exposure during embryonic development. Such effects in individuals exposed during embryonic development solely were long-lasting and persisted until the latest larval stages. Finally, we found that exposure to AMPA after hatching (during the larval stage) increased mortality. Exposure to low environmental concentrations of AMPA could have long-lasting consequences on fitness and population persistence. These findings are especially important to take into account at a time when multiple threats can interact to affect wildlife.

Naphthenic acid fraction compounds reduce the reproductive success of wood frogs (Rana sylvatica) by affecting offspring viability

Understanding the toxicity of organic compounds in oil sands process-affected water (OSPW) is necessary to inform the development of environmental guidelines related to wastewater management in Canada's oil sands region. In the present study, we investigated the effects of naphthenic acid fraction compounds (NAFCs), one of the most toxic components of OSPW, on mating behaviour, fertility, and offspring viability in the wood frog (Rana sylvatica). Wild adult wood frogs were exposed separately from the opposite sex to 0, 5, or 10 mg/L of OSPW-derived NAFCs for 24 h and then combined in outdoor lake water mesocosms containing the same NAFC concentrations (n = 2 males and 1 female per mesocosm, n = 3 mesocosms per treatment). Mating events were recorded for 48 h and egg masses were measured to determine adult fertility. NAFC exposure had no significant effect on mating behaviour (probability of amplexus and oviposition, amplexus and oviposition latency, total duration of amplexus and number of amplectic events) or fertility (fertilization success and clutch size). Tadpoles (50 individuals per mesocosm at hatching, and 15 individuals per mesocosm from 42 d post-hatch) were reared in the same mesocosms under chronic NAFC exposure until metamorphic climax (61-85 d after hatching). Offspring exposed to 10 mg/L NAFCs during development were less likely to survive and complete metamorphosis, grew at a reduced rate, and displayed more frequent morphological abnormalities. These abnormalities included limb anomalies at metamorphosis, described for the first time after NAFC exposure. The results of this study suggest that NAFCs reduce wood frog reproductive success through declines in offspring viability and therefore raise the concern that exposure to NAFCs during reproduction and development may affect the recruitment of native amphibian populations in the oil sands region.

Projected Climate and Hydroregime Variability Constrain Ephemeral Wetland-Dependent Amphibian Populations in Simulations of Southern Toads

Amphibian populations are threatened globally by stressors, including diminishing availability of suitable wetland breeding sites, altered hydroregimes driven by changing weather patterns, and exposure to contaminants. Ecological risk assessment should encompass spatial and temporal scales that capture influential ecological processes and demographic responses. Following the PopGUIDE framework of population model development for risk assessment, we used matrix population models, in conjunction with existing hydroregime predictions, under a climate change scenario to evaluate the effects of environmental stochasticity and aquatic pesticide exposure on amphibians that are dependent on ephemeral wetlands. Using southern toads (Anaxyrus terrestris) as an example, we simulated population dynamics with breeding success dependent on hydroregime suitability. Years were defined as optimal, marginal, or insufficient for successful toad recruitment, based on the duration of their potential breeding season and rate of larval development to metamorphosis. We simulated both probabilistic and chronologically specific population projections, including variable annual fecundity, based on hydroregime suitability and reduced larval survival from carbaryl exposure. In our simulations, populations were more negatively impacted by prolonged drought, and consequently multiple sequential years of reproductive failure, than by aquatic pesticide exposure. These results highlight the necessity of reliable climate projections to accurately represent the effects of altered hydroregimes on amphibian populations. Risk assessment approaches could be improved with flexible modifications that allow inclusion of various extrinsic stressors and identification of demographic and ecological vulnerabilities when precise data are lacking.

Moving beyond Risk Quotients: Advancing Ecological Risk Assessment to Reflect Better, More Robust and Relevant Methods

Under standard guidance for conducting Ecological Risk Assessments (ERAs), the risks of chemical exposure to diverse organisms are most often based on deterministic point estimates evaluated against safety-factor-based levels of concern (LOCs). While the science and guidance for mechanistic effect models (e.g., demographic, population, and agent-based) have long been demonstrated to provide more ecologically relevant effect endpoints upon which risk can be evaluated, their application in ERAs has been limited, particularly in the US. This special issue highlights the state of the science in effect modeling for ERAs through demonstrated application of the recently published Population modeling Guidance, Use, Interpretation, and Development for ERA (Pop-GUIDE). We introduce this issue with a perspective on why it is critical to move past the current application of deterministic endpoints and LOCs. We demonstrate how the current, widely used approaches contain extensive uncertainty that could be reduced considerably by applying models that account for species life histories and other important endogenous and exogenous factors critical to species sustainability. We emphasize that it is long past time to incorporate better, more robust, and ecologically relevant effect models into ERAs, particularly for chronic risk determination. The papers in this special issue demonstrate how mechanistic models that follow Pop-GUIDE better inform ERAs compared to the current standard practice.

Exposure to the synthetic phenolic antioxidant 4,4'-thiobis(6-t-butyl-m-cresol) disrupts early development in the frog Silurana tropicalis

Many chemicals in commonly used household and industrial products are being released into the environment, yet their toxicity is poorly understood. The synthetic phenolic antioxidant, 4,4'-thiobis(6-t-butyl-m-cresol) (CAS 96-69-5; TBBC) is present in many common products made of rubber and plastic. Yet, this phenolic antioxidant has not been tested for potential toxicity and developmental disruption in amphibians, a sensitive and susceptible class. We investigated whether acute and chronic exposure to TBBC would interfere with thyroid hormone-dependent developmental processes in the frog Silurana tropicalis and thus affect its early life-stage development. We exposed S. tropicalis embryos at the Nieuwkoop-Faber (NF) 9-10 stage to TBBC at nominal concentrations (0, 25, 50, 75, 100, 200 and 400 μg/L) to determine the 96h lethal concentrations and sublethal effects. We conducted a chronic exposure starting at stage NF47-48 to three sublethal TBBC nominal concentrations (0, 0.002, 0.1 and 5 μg/L) for 48-52 days to evaluate effects on growth and metamorphosis. The 96h lethal and effective (malformations) TBBC concentrations (LC₅₀ and EC₅₀) were 70.5 and 76.5 μg/L, respectively. Acute exposure to all TBBC concentrations affected S. tropicalis growth and was lethal at 200 and 400 μg/L. Chronic exposure to sublethal TBBC concentrations reduced body size by 8% at 5 μg/L and body mass by 17% at 0.002 μg/L when metamorphosis was completed. This study demonstrates that TBBC is toxic, induces malformations and inhibits tadpole growth after acute and chronic exposures. These findings call for further investigations on the mode of actions of TBBC and related antioxidants for developmental disruption in amphibians.

Aminomethylphosphonic acid (AMPA) alters oxidative status during embryonic development in an amphibian species

Glyphosate's primary metabolite (aminomethylphosphonic acid, AMPA) is known to alter embryonic development at environmentally relevant concentrations in amphibians. However, we have limited understanding of the physiological mechanisms through which AMPA affects organisms. In this study, we tested whether alteration of the oxidative status is one mechanism through which AMPA affects organism performance. To this end, we analysed several oxidative status markers in hatchling tadpoles that were exposed to sublethal concentrations of AMPA during embryonic development (~16 days). We compared the influence of environmentally relevant concentrations of AMPA (from 0.07 to 3.57 μg l^-1) on the relation between developmental traits (i.e, embryonic development duration, embryonic mortality and hatchling size) and oxidative status markers known to alter homeostasis when unbalanced (superoxide dismutase (SOD), glutathione peroxidase (GPx), catalase (CAT), thiols and ratios thereof). We included measures of telomere length as an indicator of physiological state. We found that AMPA concentrations induce non-monotonic effects on some oxidative status markers with hatchlings displaying elevated antioxidant responses (elevated thiols and unbalanced SOD/(GPx + CAT) ratio). The lack of effect of AMPA on the relation between developmental traits, oxidative status and telomere length suggests that selective mortality of embryos susceptible to oxidative stress may have occurred prior to hatching in individuals less resistant to AMPA which display lower hatching success. Future studies are required to disentangle whether oxidative unbalance is a cause or a consequence of AMPA exposition. This study highlights the need to investigate effects of the metabolites of contaminants at environmental concentrations to comprehensively assess impacts of anthropogenic contamination on wildlife.

Divergence in heritable life history traits suggests potential for local adaptation and trade-offs associated with a coal ash disposal site

Globally, human activities have resulted in rapid environmental changes that present unique challenges for wildlife. However, investigations of local adaptation in response to simultaneous exposure to multiple anthropogenic selection pressures are rare and often generate conflicting results. We used an in situ reciprocal transplant design within a quantitative genetic framework to examine how adaptive evolution and phenotypic plasticity contribute to the persistence of an amphibian population inhabiting an environment characterized by high levels of multiple toxic trace elements. We found evidence of phenotypic divergence that is largely consistent with local adaptation to an environment contaminated with multiple chemical stressors, tied to potential trade-offs in the absence of contaminants. Specifically, the population derived from the contaminated environment had a reduced risk of mortality and greater larval growth and in the contaminated environment, relative to offspring from the naïve population. Further, while survival in the uncontaminated environment was not compromised in offspring from the contaminant-exposed population, they did show delayed development and reduced growth rates over larval development, relative to the naïve population. We found no evidence of reduced additive genetic variation in the contaminant-exposed population, suggesting long-term selection in a novel environment has not reduced the evolutionary potential of that population. We also saw little evidence that past selection in the ASH environment had reduced trait plasticity in the resident population. Maternal effects were prominent in early development, but we did not detect any trends suggesting these effects were associated with the maternal transfer of toxic trace elements. Our results demonstrate the potential for adaptation to multiple contaminants in a wild amphibian population, which may have facilitated long-term persistence in a heavily impacted environment.

Cross-Taxa Distinctions in Mechanisms of Developmental Effects for Aquatic Species Exposed to Trifluralin

Standard ecological risk assessment practices often rely on larval and juvenile fish toxicity data as representative of the amphibian aquatic phase. Empirical evidence suggests that endpoints measured in fish early life stage tests are often sufficient to protect larval amphibians. However, the process of amphibian metamorphosis relies on endocrine cues that affect development and morphological restructuring and are not represented by these test endpoints. The present study compares developmental endpoints for zebrafish (Danio rerio) and the African clawed frog (Xenopus laevis), 2 standard test species, exposed to the herbicide trifluralin throughout the larval period. Danio rerio were more sensitive and demonstrated a reduction in growth measurements with increasing trifluralin exposure. Size of X. laevis at metamorphosis was not correlated with exposure concentration; however, time to metamorphosis was delayed relative to trifluralin concentration. Gene expression patterns indicate discrepancies in response by D. rerio and X. laevis, and dose-dependent metabolic activity suggests that trifluralin exposure perturbed biological pathways differently within the 2 species. Although many metabolites were correlated with exposure concentration in D. rerio, nontargeted hepatic metabolomics identified a subset of metabolites that exhibited a nonmonotonic response to trifluralin exposure in X. laevis. Linking taxonomic distinctions in cellular-level response with ecologically relevant endpoints will refine assumptions used in interspecies extrapolation of standard test effects and improve assessment of sublethal impacts on amphibian populations. Environ Toxicol Chem 2020;39:1797-1812. Published 2020. This article is a US government work and is in the public domain in the USA.

Guidance for Developing Amphibian Population Models for Ecological Risk Assessment

Despite widespread acceptance of the utility of population modeling and advocacy of this approach for a more ecologically relevant perspective, it is not routinely incorporated in ecological risk assessments (ERA). A systematic framework for situation-specific model development is one of the major challenges to broadly adopting population models in ERA. As risk assessors confront the multitude of species and chemicals requiring evaluation, an adaptable stepwise guide for model parameterization would facilitate this process. Additional guidance on interpretation of model output and evaluating uncertainty would further contribute to establishing consensus on good modeling practices. We build on previous work that created a framework and decision guide for developing population models for ERA by focusing on data types, model structure, and extrinsic stressors relevant to anuran amphibians. Anurans have a unique life cycle with varying habitat requirements and high phenotypic plasticity. These species belong to the amphibian class, which is facing global population decline in large part due to anthropogenic stressors, including chemicals. We synthesize information from databases and literature relevant to amphibian risks to identify traits that influence exposure likelihood, inherent sensitivity, population vulnerability, and environmental constraints. We link these concerns with relevant population modeling methods and structure in order to evaluate pesticide effects with appropriate scale and parameterization. A standardized population modeling approach, with additional guidance for anuran ERA, offers an example method for quantifying population risks and evaluating long-term impacts of chemical stressors to populations. Integr Environ Assess Manag 2020;16:223-233. © 2019 SETAC.