REVIEW article

Front. Ecol. Evol., 13 March 2024

Sec. Ecophysiology

Volume 12 - 2024 | https://doi.org/10.3389/fevo.2024.1336747

Framework for multi-stressor physiological response evaluation in amphibian risk assessment and conservation

  • JA

    Jill A. Awkerman 1*

  • DA

    Donna A. Glinski 2

  • WM

    W. Matthew Henderson 2

  • RV

    Robin Van Meter 3

  • ST

    S. Thomas Purucker 4

  • 1. Center for Ecosystem Measurement and Modeling, Office of Research and Development, US Environmental Protection Agency, Gulf Breeze, FL, United States

  • 2. Center for Ecosystem Measurement and Modeling, Office of Research and Development, US Environmental Protection Agency, Athens, GA, United States

  • 3. Environmental Science and Studies, Washington College, Chestertown, MD, United States

  • 4. Center for Computational Toxicology and Exposure, Office of Research and Development, US Environmental Protection Agency, Durham, NC, United States

Article metrics

View details

11

Citations

3,6k

Views

1,1k

Downloads

Abstract

Controlled laboratory experiments are often performed on amphibians to establish causality between stressor presence and an adverse outcome. However, in the field, identification of lab-generated biomarkers from single stressors and the interactions of multiple impacts are difficult to discern in an ecological context. The ubiquity of some pesticides and anthropogenic contaminants results in potentially cryptic sublethal effects or synergistic effects among multiple stressors. Although biochemical pathways regulating physiological responses to toxic stressors are often well-conserved among vertebrates, different exposure regimes and life stage vulnerabilities can yield variable ecological risk among species. Here we examine stress-related biomarkers, highlight endpoints commonly linked to apical effects, and discuss differences in ontogeny and ecology that could limit interpretation of biomarkers across species. Further we identify promising field-based physiological measures indicative of potential impacts to health and development of amphibians that could be useful to anuran conservation. We outline the physiological responses to common stressors in the context of altered functional pathways, presenting useful stage-specific endpoints for anuran species, and discussing multi-stressor vulnerability in the larger framework of amphibian life history and ecology. This overview identifies points of physiological, ecological, and demographic vulnerability to provide context in evaluating the multiple stressors impacting amphibian populations worldwide for strategic conservation planning.

1 Stressors

Multiple common sources of physiological stress contribute to the ubiquitous threats to amphibian populations worldwide, including disease, pollution, and habitat loss as well as combinations of these stressors (Stuart et al., 2004; Wake and Vredenburg, 2008; Foden et al., 2013; Grant et al., 2016; Green et al., 2020). Stressor impacts can be detected at the organismal level before long-term population decline is apparent. Habitat constraints are frequently observed as higher density resource competition inhibiting metamorphosis, recruitment, or reproductive success in some species (Harper and Semlitsch, 2007; Rittenhouse and Semlitsch, 2007). Disease transmission often presents as an immunological response prior to mass mortality (Ohmer et al., 2021). Pollution, likewise, can result in reduced reproductive success or growth in addition to mortality, and the chronic effects of these stressors can often be detected as systemic responses within the organism that precede impacts apparent at the population level (Thambirajah et al., 2019; Trudeau et al., 2020). Oxidative stress, compromised immunity, endocrine disruption, and altered metabolic activity are some physiological indications of perturbations in biological function that can lead to phenotypical impacts on individual fitness, with implications for population dynamics.

1.1 Habitat degradation

Habitat conversion, degradation, and fragmentation are the primary global causes of terrestrial biodiversity loss (Haddad et al., 2015; Newbold et al., 2015). Though global amphibian declines are linked to multiple stressors and their interactions, habitat loss typically plays an outsized role due to impacts on survival, gene flow, and dispersal (Sodhi et al., 2008). Spatial range, dispersal rates, or seasonal constraints influence population connectivity, and the abiotic conditions limiting habitat availability are projected to be less favorable in response to climate change (Sodhi et al., 2008; Funk et al., 2021). Warmer and drier conditions produced from changing climatic trends provide a direct thermal stressor and are expected to accelerate habitat loss of ephemeral wetlands (Blaustein et al., 2010; Lertzman-Lepofsky et al., 2020). Thermal stressors can geographically constrain or shift suitable aquatic (Duarte et al., 2012) and terrestrial (Hoffmann et al., 2021) ranges, particularly for cold-adapted species and microclimate-dependent life history stages with limited acclimation capacity (Frishkoff et al., 2015).

Sources of anthropogenic modifications linked to amphibian habitat loss are driven by deforestation and urbanization (Cordier et al., 2021). Continued fragmentation of amphibian populations based on their hydroregime dependency has demonstrated that periods of drought effectively isolate numerous endangered species (Zamberletti et al., 2018; Allen et al., 2020). Further, conservation of breeding wetlands is insufficient to overcome the challenges presented by anthropogenically or climatically modified habitats (Allen et al., 2020), particularly the anticipated reduction in temporary wetland inundation (Brice et al., 2022). Additionally, wetland protection depends on legal decisions that are subject to amendment or revision. Even with ample habitat available, environmental stochasticity increases variance in juvenile recruitment for species dependent on ephemeral wetlands (Greenberg et al., 2017), particularly for species with high dispersal rates and/or an energetically costly metamorphosis (Funk et al., 2021; Brooks and Kindsvater, 2022). Hydroperiod duration could have a greater impact on metapopulation persistence than pathogen or contaminant exposure (Smalling et al., 2019), specifically anomalous deluge events or multiple years of drought (Walls et al., 2013; Awkerman and Greenberg, 2022; cf. Moss et al., 2021).

1.2 Pathogens

In addition to the limitations of habitat availability, amphibian populations are also regulated by disease and predation. Many species require fish-free breeding ponds for sufficient reproductive success and are vulnerable to predation by aquatic insects (Ohba, 2011). Anuran species and life stages vary in inherent susceptibility and ecological likelihood of exposure to waterborne pathogens such as ranaviruses and Batrachochytrium dendrobatidis (Bd) (Haislip et al., 2011; Hoverman et al., 2011). Bd, the fungus responsible for chytridiomycosis, is found in cooler, lentic waterbodies (Spitzen-van der Sluijs et al., 2017), and prevalence is often highest among amphibian larvae, with later life stages more resistant to infection (Li et al., 2021). Amphibian response to chytridiomycosis often involves the complement system, in an immunological response to the pathogen, and is frequently detected through bacteria-killing assays (BKA; Rodriguez and Voyles, 2020). Ranavirus is often detected in amphibian communities with greater species diversity (Bienentreu et al., 2022). Pathogen effects can be exacerbated by transmission via more resilient invasive species that spread disease in addition to competing for the diminishing habitat of native species. For example, the American bullfrog (Lithobates catesbeianus) is a particularly invasive species that is less susceptible to ranavirus and chytridiomycosis-induced lethality, and therefore acts as an influential vector facilitating world-wide transmission of ranavirus (Hossack et al., 2023). The global trade and subsequent farming of this species for human consumption have resulted in the detection of ranavirus in native populations from previously uncontaminated regions such as those of Brazil and Mexico (Ruggeri et al., 2019 and Saucedo et al., 2019, respectively). International trade of the invasive Xenopus laevis has also contributed to the spread of chytridiomycosis (Fisher and Garner, 2020).

Amphibian species differ in their response to the fungal pathogen Bd with some species showing downregulation of cellular and metabolic functions and upregulation of adaptive immune gene response; however, such responses are ultimately insufficient to prevent high microbial infection loads (e.g., Eskew et al., 2018). Other species with more diverse dermal antimicrobial peptide communities showed minimal response to infection (Eskew et al., 2018). Lower temperatures may increase inflammation-related responses as opposed to warmer temperatures increasing adaptive immune responses (Ellison et al., 2020). More bacterial reads, presumably from frog microbiomes, were found in populations with a history of ranavirus (Campbell et al., 2018). Differential impacts of changing climate on host and pathogen further complicate strategies to prevent transmission (Blaustein et al., 2012). It is likely that warming climates will impact viral loads, as observed in juveniles at warmer temperatures with less intense but persistent infections (Brunner et al., 2019), and bacteria-killing ability is reduced at higher temperatures in some species (Rodriguez and Voyles, 2020). Coinfection of ranavirus and chytrid in several endemic tadpoles underscores the importance of understanding the etiology and interactions of these pathogens for effective conservation of amphibians and other aquatic vertebrates (Warne et al., 2016).

1.3 Pesticides

Agricultural and residential pesticide use has also been implicated as a contributing factor in declining amphibian populations (Hayes et al., 2010; Brühl et al., 2011, 2013) with agriculture identified as the most common cause of extinction threats for amphibians and other terrestrial vertebrates (Munstermann et al., 2022). A meta-analysis of pesticide effects revealed moderate impacts on survival and decreased mass and relatively greater impacts from deformities not associated with phylogeny. Although contaminants of emerging concern were underrepresented in pollutant studies, pesticide effects were comparable with those of wastewater, less impactful compared to deicer effects, and relatively greater than those of metals and phosphorus compounds (Egea-Serrano et al., 2012). Additionally, transgenerational impacts, lethal and sublethal, have been demonstrated from exposure to environmentally relevant pesticide concentrations (Karlsson et al., 2021; Usal et al., 2021). The amphibian life cycle allows complex exposure dynamics in both aquatic and terrestrial environments, and recommended application rates of many pesticides result in high mortality from terrestrial exposure (Brühl et al., 2013), although terrestrial effects are less frequently documented. Indirect effects of pesticide use at lower concentrations than those toxic to amphibians potentially impact the full lifecycle of amphibians through reduction of resources, although aquatic food web effects are more frequently reported than terrestrial food web effects (Relyea and Diecks, 2008; Relyea, 2009). Overall, aquatic pesticide exposure can alter various endocrine functions important to development and reproduction and result in a variety of systemic impacts in amphibians (Thambirajah et al., 2019). A recent review of endocrine disruption by agrochemicals summarized changes in lipid and energy metabolism among fungicides; effects on metabolism, metamorphic success, and gonadal development for some herbicides; and reduced metamorphosis from fertilizer and other pesticide exposures (Trudeau et al., 2020). Evaluating the non-lethal effects of pesticides is complicated by timing of exposure and sample collection as well as tissue type, such that measured effects vary depending on species, mechanism of action, route of exposure, and the concentration of the compound (Rohr and McCoy, 2010; Glinski et al., 2021; Seim et al., 2022). Even with the abundance of scientific support correlating pesticide exposure to declining amphibian populations, it is unrealistic that impacts of pesticide exposure will be reversed, given the moderate generation times of most amphibians, the complexity of potential exposure based on their life cycle, and the substantial proportion of croplands in protected areas associated with continuing tradeoffs between food security and conservation (Vijay and Armsworth, 2021).

1.4 Stressor interactions

Uncertainty surrounding individual response, species vulnerability, and exposure regime complicates risk assessment determinations of multiple stressor impacts at the landscape level (Relyea and Hoverman, 2006). For instance, co-stressors such as heat, pesticides, and parasites impact amphibian immune responses and can have synergistic effects on fecundity and post-recruitment survival (Kiesecker, 2002, Thompson et al., 2022). When anthropogenic stressors and abiotic factors synergize, the immune system is challenged (Kiesecker, 2011), and early stress experienced during development can affect resilience in later life stages (Kohli et al., 2019; Le Sage et al., 2022). Disease susceptibility can increase following herbicide exposure (Rohr et al., 2013), and lower microbiota diversity, a common result of pesticide exposure, is associated with reduced parasite resistance (Knutie et al., 2017). Anticipating potential long-term effects in response to various stressors and their interactions, which can promulgate into subsequent life stages, challenges both establishing in situ causality from single stressors needed for tighter regulations and effective conservation management. Ultimately, ecological risk assessment is complicated not only by a deficit of toxicological data, but also a lack of ecological data to document changes in land use, species abundance and distribution, and disease transmission that are necessary for adaptive management approaches (Womack et al., 2022). Extrinsic stressors are presented in Table 1 along with physiological measurements of these effects, endogenous processes affecting the same biochemical pathways, and life stages in which departures from typical functions are detectable and/or problematic (Figure 2).

Table 1

StageSystemic responseEndogenous activityExogenous stressorsAssessment methodsDemographic endpoint
Embryolarval developmentMetabolismEnergy production, DNA synthesis, protein synthesis↑, alanine ↓ aspartate ↓, α-ketoglutamine ↑ (Vastag et al., 2011)Aquatic conditions, including chemical pollution, pathogens, predation (Kiesecker, 2011) Whole-organism metabolites (Vastag et al., 2011)Embryo mortality; cessation of development
MetamorphosisMetabolismIncreasing energy needs, anabolic activity, tail apoptosis; greater dehydration from fasting (Rowland et al., 2023); purine, arginine and pyrimidine, urea cycle metabolites, arginine and purine/pyrimidine, cysteine/methionine, sphingolipid, and eicosanoid metabolism (Ichu et al., 2014)↑ galactose metabolism and lactose degradation with xenobiotic exposure (Glinski et al., 2021) or reduced resources; ↓ glutathione (Ichu et al., 2014); galactose predictive of chytrid (Wang et al., 2021); predation and pesticide exposure alter aminoacyl-tRNA biosynthesis, galactose and glutathione metabolism, arginine biosynthesis (Snyder et al., 2022); pesticide exposure impacts serine and threonine, histadine, linoleic acid, and sphingolipid metabolismWhole-organism or tissue metabolomics (Ichu et al., 2014)Delayed development, reduced transition to juvenile stage
Redox signallingLipid peroxidation ↑, glutathione ↓, catalase ↓, SOD, CAT, MDA expression altered during intestinal development and tail resorption, ascorbic acid ↑ for collagen synthesis (Menon and Rozman, 2007; Guo et al., 2022); glutathione peroxidase ↓, GST ↓, sulfhydryl groups ↓ (Petrović et al., 2021)Lower antioxidant activity and increased lipid peroxidation to xenobiotics or environmental conditions (abiotic or density effects; Burraco et al., 2017; Petrović et al., 2021); ↑ thiol and CAT in pesticide and nematode infection (Marcogliese et al., 2021)ROS production in tissues; antioxidant enzymatic responses of SOD, CAT, MDA, GST (Menon and Rozman, 2007; Chen et al., 2017; Guo et al., 2022); decreased expression of GSH; increased TBARSDelayed development; reduced transition to juvenile stage
Endocrine responseCS in response to ↑ TH, regulate development via diodination, glucuronidation, sulfation, affecting HPT, HPA, HPG axes (Denver, 2009; Duarte-Guterman et al., 2014; Thambirajah et al., 2019)GC ↑ to some xenobiotics (Burraco and Gomez-Mestre, 2016; Trudeau et al., 2020), environmental conditions (Sachs and Buchholz, 2019; Thambirajah et al., 2019), predators (Narayan et al., 2013) ; neurogenerative, oxidative, mitochondrial, teratological effects (Di Lorenzo et al., 2020)CS and TH levels in tissue or immersed water (Gabor et al., 2013a); tissue/organism enzyme activity or DGE in AR, TR, tra, trb, dio2, dio3 (Thambirajah et al., 2022); ambient water assay; size at metamorphosis (Rowland et al., 2023); vitellogenin indicative of feminization (Venturino and de D’Angelo, 2005)Time to metamorphosis; cohort sex ratio; carryover to juvenile immunity, survival, fecundity (Kiesecker, 2002, Denver, 2009; Kohli et al., 2019; Ruthsatz et al., 2020; Le Sage et al., 2022)
ImmunityEndocrine-driven development of immunity; immunosuppression at metamorphosis (Rollins-Smith, 2017)Viral loads, resistance, and parasite prevalance affected by pesticides and abiotic factors (Kerby et al., 2011; Kiesecker, 2011; Knutie et al., 2017; Pochini and Hoverman, 2017); Microbiome in tadpoles impacted by xenobiotic exposureGut microbiome diversity; at advanced developmental stages – blood leukocytes, white cell lymphocytes and granulocytes (basophils, neutrophils, eosinophils); DGE (Row et al., 2016)Reduced survival due to pathogens and parasites (Kiesecker, 2011)
Juvenile maturation to adultEndocrineTRH influences TSH (Paul et al., 2022)Food constraints ↓ CORT (Prokić et al., 2021); variance in CORT along latitudinal cline (Le Sage et al., 2022)Dermal swab, fecal content, tissue or ambient water assay of CS (Gabor et al., 2013a)Behavioral responses to stressors, reduced dispersal
ImmunityGut microbiome linked to resistance of parasites (Knutie et al., 2017), skin microbiome linked to resistance of pathogens (Krynak et al., 2017; McCoy and Peralta, 2018; Jiménez et al., 2021); possibly compromised by shortened developmental hydroperiod (Brannelly et al., 2019); Lower juvenile immunity relative to mature adultsDermal microbiome and pathogen vulnerability impacted by xenobiotic exposure (Krynak et al., 2017; McCoy and Peralta, 2018; Jiménez et al., 2021); habitat degradation affects vulnerability to pathogens (Stevens and Baguette, 2008; Costa et al., 2021; Becker et al., 2023)Microbiome diversity in skin mucosa (Neely et al., 2022); antimicrobial peptides (Huang et al., 2016); white cell lymphocytes and granulocytes (basophils, neutrophils, eosinophils); B and T cells in organs, MHC-II; antibodies - IgA/X, IgD, IgF, IgM, IgYSusceptibility to pathogens, reduced juvenile survival or limited dispersal due to disease or deformities (Kiesecker, 2002, Rohr et al., 2006; Krynak et al., 2017; Kohli et al., 2019)
MetabolismRelated to endocrine activity; longer hydroperiod ↑ lipid stores (Scott et al., 2007)Influenced by temperature, water loss, xenobiotics; food deprivation ↓ CORT; pesticides altered sucrose and starch pathway regulation (Zaya et al., 2011, Dornelles and Oliveria, 2016, Van Meter et al., 2018)Body condition; energy metabolism in tissueReduced juvenile survival
Redox signallingIncreased antioxidants during estivation in preparation for oxidative stress; lower oxidative metabolism enzyme activity during estivation (Rowland et al., 2023)Food constraints ↑ lipid peroxide; ↓ SOD, glutathione peroxidase, GST, glutathione and sulfhydryl groups (Prokić et al., 2021); pesticides and pathogens ↑ thiol; nematode infections ↑ thiol, ↑ catalase (Marcogliese et al., 2021)ROS production; anti-oxidant enzymatic responses of SOD, glutathione peroxidase, glutathione, GST, thiol, catalase (Marcogliese et al., 2021; Prokić et al., 2021)Survival to following breeding season, potentially a function of size/condition at end of season
Adult fecundityEndocrine activityGonadotropins released by pituitary; estrogen, androgen, progestogen regulate reproduction (Duarte-Guterman et al., 2014)Endocrine disrupting compounds can disrupt gonadal development, sexual differentiation (Lambert et al., 2015; Marlatt et al., 2022); temperature impact on sex determination (Lambert et al., 2015; Ruiz-García et al., 2021); density-dependent resource availability (Kissel et al., 2020)ER/AR binding; Aromatase inhibition; impairment of steroidogenesis; vitellogenin expression in response to xenoestrogen exposure; zona radiata, zona pellucida, DGE in er, bteb, tra, trb, thbzip,Altered population sex ratio (Roco et al., 2021; Baranek et al., 2022); reduced fecundity

Endogenous activity associated with transitional phases of the anuran lifecycle, exogenous stressors that alter biological processes, methods to assess organismal effects, and potential demographic impacts. ↑ upregulation or increased expression; ↓ downregulation or decreased expression.

Endogenous activity associated with transitional phases of the anuran lifecycle, exogenous stressors that alter biological processes, methods to assess organismal effects, and potential demographic impacts. ↑ upregulation or increased expression; ↓ downregulation or decreased expression.

Figure 2

2 Lifestage-specific physiology

Effective adaptive conservation management strategies target vulnerable life stages and critical threats to wildlife populations. The biphasic life cycle of anuran amphibians makes them particularly vulnerable to extrinsic stressors because of their dependence on variable aquatic habitat resources as well as terrestrial environmental quality (Nolan et al., 2023). Their complex life history strategy and multiple potential drivers of population decline require a more nuanced approach to targeting spatial and temporal variability in stressors relative to life stage (Walls and Gabor, 2019; Awkerman et al., 2020). Distinguishing stage-specific endogenous variation in physiological processes enables anticipation of compromised physical condition in response to common stressors (Brooks and Kindsvater, 2022; Nolan et al., 2023). Here we focus on stressor impacts on transition between life stages (F, Te, Tl, Tj in Figure 1) but present potential effects on survival and development as well.

Figure 1

2.1 Development (embryo transition to tadpole stage; Te)

Survival during the relatively brief stage of embryo development is largely dependent on a suitable environment to avoid predators, pathogens, or pollution, and the costs associated with such defenses differ among amphibian species and developmental mode (Brooks and Kindsvater, 2022). Amphibian clutches can experience high mortality from pathogens or predation, depending on the geographic location and ecological community composition (Kiesecker, 2011), such that habitat characteristics and regional observations are most informative in identifying these threats (Wake and Vredenburg, 2008). Pesticides aggregated as runoff in wetlands provide another potential stressor for developing embryos (Smalling et al., 2015) with lethal or sublethal impacts on individuals. Ultimately, a systematic review revealed that the time to hatching for embryos was influenced more by taxonomy and exposure to pollution, rather than experimental setting (lab vs. field; Egea-Serrano et al., 2012). Singly or in combinations, stressors during embryogenesis can lead to delayed, wide-ranging effects, resulting in a diverse array of phenotypic outcomes associated with aspects of developmental plasticity that are not observed until later life history stages (Jonsson et al., 2022).

Given the relatively brief duration of this stage in most anuran species, and rapidly changing metabolism, identifying potential stressors based on organismal condition or response could be a challenging diagnostic approach, compared to assessment of anomalous response during later life stages. Endogenous variation in embryonic metabolite levels is suggestive of energy production primarily, presumably for DNA synthesis (Vastag et al., 2011). Contaminant levels in egg masses that are linked to deformities and reduced offspring viability can result from maternal transfer of contaminants rather than indicating direct environmental exposure alone (Todd et al., 2011; Metts et al., 2013). Determining physiological response to a variety of stressors (e.g., contaminant mixtures and abiotic factors), is a complex challenge that might be approached by evaluating exposure-based epigenetic changes (e.g., DNA methylation, histone acetylation) in developing embryos (Fogliano et al., 2023) or simply assessing differential responses in later life stages.

2.2 Metamorphosis (transition from larval to juvenile stage; Tl)

The morphological restructuring for transition from larval stage to juvenile stage is dependent on endocrine drivers, specifically surges in thyroid hormones (TH), regulated by thyroid hormone receptors and retinoic acid receptors (TR and RXR, respectively; reviewed in Paul et al., 2022). Endocrine regulation and body morphogenesis during the larval stage are controlled by the hypothalamic-pituitary-thyroid (HPT) and hypothalamic-pituitary-adrenal/interrenal axes (HPA/HPI) as well as the hypothalamus-pituitary-gonadal (HPG) axis (Duarte-Guterman et al., 2014). Development and metamorphosis are regulated largely by the release of the thyroid hormones thyroxine (T4) and tri-iodothyronine (T3) and modulated by the corticosteroids (CS) corticosterone (CORT) and aldosterone (ALDO; Denver, 2009). Regulation of TH signal involves cellular processes of deiodination, glucuronidation, and sulfation (Thambirajah et al., 2019). Metabolism and cardiac functions associated with development and metamorphosis are also regulated by CS. Ruthsatz et al. (2020) showed that certain metamorphic stages were significantly more susceptible to changes in growth and development due to increased TH levels, with high TH levels associated with reduced weight and size in tadpole and froglet stages as well as increased heart rate and reduced energy stores across all stages.

TH inhibition or impairment can delay development, while CS production is often associated with accelerated metamorphosis in response to pond drying or other stressors (Sachs and Buchholz, 2019; Thambirajah et al., 2019), although the role of CORT as a homeostatic response to stress is complex. The corticotropic releasing hormone (CRH) regulates the HPA axis as well as the HPT axis, thereby contributing to additional crosstalk between these pathways and circulating hormone levels (Thambirajah et al., 2022). CORT levels in southern leopard frogs increased with exposure to multiple aquatic stressors, specifically a nitrogenous fertilizer, a pesticide, and salt (Adelizzi et al., 2019). However, relatively elevated CORT levels were associated with populations less tolerant to contaminant exposure, such that differences in stress response could be indicative of exposure history (Shidemantle et al., 2022). Predator detection can also elicit an increased CORT response (Narayan et al., 2013). Signals of agrochemical disruption of endocrine function among interactions of the thyroid, gonadal, and metabolic axes in amphibians was reviewed in detail by Trudeau et al. (2020). Early life stage stressors that elevate CS production can alter endocrine response throughout the lifecycle of the individual (Denver, 2009).

Stressor perturbations in endocrine functions are particularly impactful in metamorphosing amphibians and can influence immunity, survival, and fecundity in subsequent terrestrial life stages (Kiesecker, 2002; Kohli et al., 2019). A meta-analysis determining effects on time to metamorphosis found taxonomy, pollution, and timing of exposure to be more influential than the experimental setting (Egea-Serrano et al., 2012). Pond drying constraints influencing larval development are expected to be impacted in various ways by climate change, depending on regional location (Walls et al., 2013). The duration of larval stage and developmental mode, combined with community dynamics between larval competitors and predators, can distinguish species resilience and response to such unpredictable environmental stressors (Belden et al., 2010; Moss et al., 2021; Brooks and Kindsvater, 2022). Interannual variability in hydrologic regime at temporary wetlands determines the length of the developmental period and the density of developing anuran larvae (e.g., Pechmann et al., 1989). Developmental plasticity in metamorphic climax allows variable phenotypic response to interannual conditions and is driven by the neuroendocrine processes responsible for the development of the immune system, highlighting a potential tradeoff between accelerated development and resistance to disease and parasites (Kohli et al., 2019). Likewise, tradeoffs between development and microbiota diversity or immunology are demonstrated later in life with increased susceptibility to pathogens (Le Sage et al., 2022). As northern leopard frog tadpoles approach metamorphic climax, tail tissue decreases expression of mitochondrial energy genes and upregulates expression of immunity genes (Row et al., 2016). Post-metamorphic immune function may be compromised in amphibians experiencing shorter hydrologic regimes (Brannelly et al., 2019), which may further exacerbate disease susceptibility. Therefore, interannual variance in aquatic habitat suitability can have lasting impacts on cohort fitness.

During metamorphic climax, metabolic activity changes, reflecting increasing energy needs, anabolic requirements, and tail apoptosis; these energetic requirements and fasting effects create a vulnerable transition from larva to juvenile in which contaminant body burdens can amplify (Rowland et al., 2023). The aquatic phase of the amphibian life cycle is also susceptible to reduced growth in response to pathogens that have been introduced to waterbodies, although survival is rarely impacted at this stage (Nolan et al., 2023). Although phylogeny and abiotic environmental variables determine the initial likelihood of Bd or ranavirus occurrence in areas of viral compatibility, other stressors such as pesticides can further influence viral loads and the resistance of the host population as well as the prevalence of parasite-induced deformities (e.g., Kerby et al., 2011; Kiesecker, 2011; Pochini and Hoverman, 2017).

2.3 Maturity (juvenile transition to reproductive adult; Tj)

The literature on this critical amphibian life stage is scarce, due in part to the complexity of rearing and maintaining juvenile amphibian populations in a laboratory setting through maturation, as well as the challenge of monitoring individual juvenile amphibians from metamorphosis through reproduction in a field setting (see Petrovan and Schmidt, 2019). Furthermore, a systematic analysis of factors affecting survival found pollutants and the experimental setting (lab vs. field) to be more influential than taxonomic group or developmental stage in the study (Egea-Serrano et al., 2012). Although pre-metamorphic environmental conditions directly influence post-metamorphic life stages, there may also be distinctly different age or stage-specific stress responses in amphibians, as evidenced by variations in sucrose and starch pathway regulation following pesticide exposure in larval and juvenile amphibians (Zaya et al., 2011; Dornelles and Oliveria, 2016; Van Meter et al., 2018). Survival to reproduction was positively correlated with lipid stores at metamorphosis among two Ambystoma salamander species, and lipid stores were greater among individuals emerging from longer hydroperiod wetlands. Furthermore, total rainfall during years of juvenile development was also positively associated with survival to reproduction (Scott et al., 2007).

Potential carry-over effects from compromised development can exist (Kiesecker, 2002; Rohr et al., 2006; Kohli et al., 2019) with additional risk from stressors in the terrestrial environment. Terrestrial habitat degradation and habitat fragmentation influences the connectivity of both amphibian populations and their potential pathogens (Stevens and Baguette, 2008; Costa et al., 2021; Becker et al., 2023). Reduced skin and gut microbiota in the larval stage can also reduce parasite resistance in adults (Knutie et al., 2017). Amphibian skin contains antimicrobial peptides linked to immunity and defense functions as well as to biosynthesis and metabolism (Huang et al., 2016). Skin secretions have demonstrated antimicrobial antioxidant properties and can be beneficial to healing (Wang et al., 2020). Bacterial and fungal taxonomy in skin secretions is associated with disease resistance (Bates et al., 2022), with the skin microbiome affected by the same abiotic factors that influence Bd occurrence (Ruthsatz et al., 2020). Some species’ secretions contain sufficient toxins to be lethal to predators, thereby reducing mortality via predation (Liu et al., 2022). Compromised skin microbiome diversity is implicated in important physiological functions such as electrolyte and hydration loss, disease susceptibility, and increased pesticide effects. Pesticide exposure has been linked to disruption of the skin microbiome and antimicrobial peptides of amphibians, thereby affecting vulnerability to pathogens (Krynak et al., 2017; McCoy and Peralta, 2018; Jiménez et al., 2021). Habitat disturbance has also been associated with skin microbiome diversity, primarily via pathogen dispersal (Neely et al., 2022), underscoring the ecological complexity of proximate mechanisms of multiple stressors and their potential interactions. Enhanced data collection efforts on juvenile amphibians are essential to improve risk assessment and inform management decisions at the local scale.

2.4 Fecundity (adult production of embryos; F)

Reproductive failure associated with insufficient hydroperiod is a determinant of lifetime reproductive success in species dependent on ephemeral wetlands for breeding (Taylor et al., 2006; Stevens and Baguette, 2008). In years with suitable hydroregime, terrestrial density dependence and sex ratio can affect fecundity within a population (Kissel et al., 2020; Baranek et al., 2022). Effects of endocrine disruption in developmental phases as well as during gamete production could also reduce fecundity via altered gonadal development, or a high incidence of intersex individuals in the population (Lambert et al., 2015; Marlatt et al., 2022). Sex-specific age at maturation could further restrict operational sex ratio in amphibian populations (Baranek et al., 2022). In addition to the endocrine disruption associated with xenobiotic exposure, temperature can affect sex determination, with potential impact on operational sex ratio following extended periods of anomalous temperatures (Lambert et al., 2015; Ruiz-García et al., 2021). The lasting impact of such shifts will vary depending on the species life history and genetic sex determination (Bókony et al., 2017).

3 Field-based measures of stressor response

Assessing the status of a wildlife population or relative condition of an individual within its habitat is a challenge complicated by the heterogeneity of both organismal response and stressor distribution. Acquiring an adequate sample size for a meaningful detection of environmental or stressor effects could limit the practical scope of most field-based efforts, while standardization of conditions can bias the interpretation of stressor response in most laboratory or mesocosm settings. Stage-specific physiology, along with ecological or life history vulnerabilities, provides additional context for interpretation of potential stressor effects (Venturino et al., 2003). For example, intestinal development and tail resorption in larvae are coincident with signs of oxidative stress (Menon and Rozman, 2007). Establishing baseline physiology with common biomarkers provides context of endogenous variability during the amphibian life cycle. Identifying these informative endpoints and sensitive stages can preclude the need for extensive accounting of stressor-specific effects and interactions.

3.1 Physiological processes

Systematic responses to stress include endocrine disruption, oxidative stress, metabolic perturbation, and compromised immunity (Venturino et al., 2003). Specifically, elevated CORT and standard metabolic rate as well as decreased antioxidant enzymes are common signals of abiotic and xenobiotic physiological stress (Burraco and Gomez-Mestre, 2016). Endocrine disruption in the interconnected hormonal axes can also trigger responses in other systems, such as neurodegenerative effects, oxidative damage, impairment of mitochondrial function, and teratological effects (Di Lorenzo et al., 2020). Taxonomic family and pollutant exposure were significant determinants of developmental abnormalities in a systematic review of ecotoxicological studies (Egea-Serrano et al., 2012), and specific phenotypical abnormalities can be ascribed to different classes of chemicals (Venturino et al., 2003).

Endocrine-driven developmental processes are highly conserved in vertebrates (Paul et al., 2022), as are many physiological endpoints associated with both homeostatic and lethal responses to pesticides and contaminant exposure. Larval amphibians are especially susceptible to endocrine disruption due to their reliance on hormonal cues for initiation and timing of metamorphosis and sex determination (Duarte-Guterman et al., 2014). Crosstalk between hormonal axes includes an evolutionary history of HPT and HPG signaling (Thambirajah et al., 2022). Genetic sex determination during developmental stages varies between and within amphibian species due to rapid turnover of genes such that either males or females can be heterozygotic, with some species having three sex chromosomes (Miura, 2017). Sex reversal in response to external temperature or steroid hormones can also affect the sex ratio of a cohort (Roco et al., 2021). Although estrogenic and androgenic effects have been studied much more extensively in mammals, intersex amphibian larvae resulting from reproductive steroid hormone exposure have been associated with effects on the androgen receptor (AR) and thyroid receptor (TR) and altered expression of dio1, dio2, dio3, and thrb (Thambirajah et al., 2022). Increased vitellogenin production is a common indication of feminization (Venturino and de D’Angelo, 2005), and increased formic acid has been suggested as an indicator of androgen receptor binding and anti-androgenic effects in larvae (Melvin et al., 2018).

Endogenous changes in metabolism are also associated with lifecycle-dependent physiological processes (Ichu et al., 2014). During metamorphosis, metabolic pathways are dramatically altered in the liver and the tail as a result of lipid and carbohydrate metabolism (Zhu et al., 2020). Metabolomic changes during metamorphosis demonstrate physiological processes associated with morphological restructuring in metabolic pathways, including the urea cycle as well as arginine and purine/pyrimidine, cysteine/methionine, sphingolipid, and eicosanoid metabolism; however, similar metabolite expression in humans is associated with disease (Ichu et al., 2014). As the tadpole tail regresses and intestines restructure, lipid peroxidation is increased; depleted catalase (CAT) and glutathione contribute to oxidative stress, as demonstrated by CAT, superoxide dismutase (SOD), and malondialdehyde (MDA) expression; and the antioxidant ascorbic acid increases as organs develop (Menon and Rozman, 2007; Guo et al., 2022). Epidermal galactose levels, and specifically 25 uniquely expressed genes, are predictive of chytrid outbreaks and are life stage dependent, with higher expression at metamorphosis (Wang et al., 2021). Food constraints in the juvenile stage were associated with higher lipid peroxidase and lower SOD, glutathione S-transferase (GST), glutathione peroxidase, glutathione and sulfhydryl groups (Prokić et al., 2021).

Antioxidant system response (AOS) and oxidative stress is highest at metamorphic peak, and associated with lower glutathione, CAT, glutathione peroxidase, GST, and sulfhydryl groups, and oxidative stress is exacerbated by decreasing water levels (Petrović et al., 2021). Hepatic GST activity has been proposed as a biomarker indicative of TH signaling imbalance and developmental effects (Chen et al., 2017). Upregulated pathways include transamination and the urea cycle because of hepatic catabolism, TCA cycle and oxidative phosphorylation resulting from energy metabolism (although these are downregulated in the tail), and hepatic glycogen phosphorylation and gluconeogenesis (Zhu et al., 2020). Decreased activity occurred in β-oxidation and the pentose phosphate pathway, and downregulation of glycolysis, β-oxidation, and transamination in the tail accompanied reduced protein synthesis and lower energy production and consumption, although glycogenesis, fatty acid elongation and desaturation, and lipid synthesis were maintained (Zhu et al., 2020).

Indication of oxidative stress is a common detoxification response to many chemical classes and is characterized by altered expression of mixed function oxidases (MFO; e.g., CYP1A, EROD, demethylase), GSH, lipid peroxides, and antioxidant enzymes (CAT, SOD; Venturino and de D’Angelo, 2005). Oxidative stress and lipid peroxidation, as demonstrated by increased SOD and CAT activity were also associated with hepatotoxicity resulting from increasing organophosphate exposure, although GST activity was unchanged, and MDA decreased (Li et al., 2017). Common indications of oxidative stress as a detoxification response include glutathione deficits and production of GST (Venturino and de D’Angelo, 2005). Interactive oxidative stress effects of pesticide concentration and parasite abundance were observed in thiol levels of recent metamorphs, with nematode infection related to elevated thiol and catalase expression (Marcogliese et al., 2021).

Amphibian physiological responses to environmental stressors have been well documented (Park and Do, 2023), and the various threats that impact amphibian populations can elicit similar physical effects. For example, xenobiotic exposure or threatening environmental conditions (e.g., pond drying or predator presence) is commonly associated with oxidative stress and production of reactive oxygen species (ROS; Burraco et al., 2017). Habitat fragmentation and degradation, coincident with anthropogenic infringement and climate change, contribute to invasive species introduction, disease outbreak, and increased pollution, multiplying threats to immunocompetency (Kiesecker, 2011). Immune functions impacted by common amphibian stressors and their interactions are indicated in various matrices and measurements. For example, glucocorticoids (GC) are CS hormones influencing the immune system, tissue inflammation, and cardiovascular response (Rollins-Smith, 2017), and frequently indicate physiological stress. However, some stressors, e.g., food deprivation, can yield differential endocrine responses, with reduced CORT in juveniles conserving energy resources as opposed to increased CORT levels in food-deprived tadpoles (Crespi and Denver, 2005).

3.2 Omics technologies

Stressor-specific measurements of organismal response introduce complexity to both laboratory and field-based assessment approaches, as well as to the interpretation of multi-stressor scenarios. Evaluating biomarker expression can help identify biochemical perturbations indicative of systemic stress to environmental conditions. The suite of ‘omics technologies, including genomics, transcriptomics, proteomics, and metabolomics, can shed light on the underlying biological processes and provide a means to identify specific genes, metabolites, and pathways that are affected in an amphibian ecological risk assessment.

Comparative genomics is increasingly recognized as a valuable tool for conservation purposes (Kosch et al., 2023). This includes the use of reference genomes in eDNA approaches for monitoring populations (Breton et al., 2022; Saeed et al., 2022), informing genetic rescue efforts for threatened amphibians (Kosch et al., 2023), and using sequence information to predict protein similarity and infer ecotoxicological implications across species (LaLone et al., 2023). However, compared to other vertebrate classes, genomic coverage for amphibians is currently recognized as lacking (Hotaling et al., 2021). Kosch et al. (2023) provided a review of the 32 available amphibian reference genomes and found variable annotation quality for the available genomes and uneven coverage across amphibian families, with genomic comparison further complicated by the presence of large, repetitive genomes. This limited availability of amphibian reference genomes presents challenges for generalizing ecotoxicological results to higher taxonomic levels within the class Amphibia. Despite these challenges, genetic approaches can provide conservation insights. This is particularly true for amphibian species with cryptic habits and biphasic life cycles, which complicate traditional field-based measurements of survival, fecundity, and migration (Mazerolle et al., 2007; Funk et al., 2021). Landscape genetics, for instance, can reveal connectivity within a population as well as isolated subpopulations (Watts et al., 2015). This information can then be used to prioritize conservation targets for threatened amphibians (e.g., Forester et al., 2022). The pressing need for increasing knowledge of amphibian genomes to assist in conservation efforts was highlighted by Calboli et al. (2011). It is hypothesized that only a relative few, simple mechanisms of gene alterations are indicated in amphibians’ response to numerous environmental stressors. Functional genomics has been used to probe the molecular underpinnings of field observations concerning the sexual differentiation in amphibians (Bögi et al., 2002), fragmentation of populations (McCartney-Melstad et al., 2018), and pathogen-host interactions (Zamudio et al., 2020).

In the laboratory, transcriptomics approaches leverage differential gene expression (DGE) approaches by contrasting the expression level of transcripts in stressed individuals versus control individuals. Changes in gene expression can reveal which genes are upregulated or downregulated, thereby identifying perturbations in specific biochemical pathways regardless of the origin of the stressor. The magnitude of the response could indicate functional points of departure (e.g., Ewald et al., 2022; Mittal et al., 2022). Transcriptomics data, generated from controlled laboratory exposures, provide a comprehensive view of gene expression changes comparable to traditional apical endpoints. The large volume of data, coupled with the fact that the expression responses are specifically associated with the mechanism of the stressor, suggests the possibility of developing expression-based “fingerprints” or signatures resulting from single and multi-stressor exposures. These can be used to determine if the magnitude of an exposure to a toxicant or stressor of a particular mode of action is likely to elicit biological perturbations that can be linked to or predictive of apical effects. Furthermore, high-throughput transcriptomics (HTTr) methods have been developed to observe changes in gene expression in cells, rather than in test species, after exposure to chemicals (Schirmer et al., 2010). These methods are less resource-intensive than traditional toxicity testing and can be used to determine at what concentration chemicals impact cellular biology and to develop adverse outcome pathways (AOPs). For vertebrates, such regulatory testing programs aim to evaluate the potential endocrine-disrupting effects of chemicals, utilizing the conservation of certain endocrine pathways among vertebrate classes to evaluate the feasibility of extrapolating data across taxa. These approaches integrate functional genomics with transcriptomics to establish the confidence levels in pathway conservation while identifying the specific needs for additional data to advance read-across methods for estrogen, androgen, thyroid, and steroidogenesis pathways in vertebrate ecological receptors (McArdle et al., 2020).

Metabolomics technology may also provide a means to address the uncertainties surrounding chemical risk assessment of single and multiple stressors. Available technology measures the changes in hundreds (if not thousands) of metabolites simultaneously, effectively capturing a metabolomic fingerprint as a snapshot of an organism’s altered physiology. This metabolomic fingerprint of subcellular biological responses often represents immediate or early response within the organism to stresses and can be associated with signaling networks that are linked to adverse outcomes at higher levels of biological organization. Successful application of metabolomics to differentiate multi-stressor response was achieved by Snyder et al. (2022). Similarly, the use of transcriptomics and proteomics for advancing amphibian toxicogenomic studies was reviewed in Helbing (2012). Relying on ‘omics technologies to identify meaningful suites of stressor response and target demographic vulnerabilities for sample collection could offer a comparative physiology approach for detecting impacted individuals and populations.

Exogenous impacts of xenobiotic exposure can exacerbate stressors that accompany particular life stages. Reduction in body size during metamorphosis and fasting during hibernation result in increased metabolic demands and greater body burdens of contaminants due to biomagnification (Rowland et al., 2023). Aquatic exposures of various pesticides were associated with increased galactose metabolism and lactose degradation, indicating effects on energy metabolism (Glinski et al., 2021). Pathways associated with glucogenesis and glycolysis were also indicators of energy metabolism impacts in terrestrial juvenile frog exposures (Van Meter et al., 2022). The urea cycle was frequently impacted by various pesticides, and the purine metabolism pathway was also enriched, indicating increased energy production as a response to toxicity. Reduction in glutathione levels is another common result of pesticide exposure indicative of oxidative stress in both larval and juvenile amphibians (Ichu et al., 2014). Reduced citrate, α-ketaglutarate, and fumarate were also proposed as oxidative stress biomarkers, as intermediates of the tricarboxylic acid cycle (Melvin et al., 2018).

The magnitude of altered metabolite regulation during later life stage terrestrial exposures was not indicative of bioaccumulation, and exposure to combinations of pesticide did not always have a synergistic effect on juvenile toads (Van Meter et al., 2018). In fact, extrinsic sources of urea as fertilizer at low concentrations can counteract combined pesticide effects, presumably by facilitating excretion and detoxification, although excessive doses can be detrimental (Van Meter et al., 2022). Hepatic metabolome analyses revealed altered pathways indicating stress caused by both predation and pesticide exposure; these include aminoacyl-tRNA biosynthesis, galactose metabolism, glutathione metabolism, and arginine biosynthesis (Snyder et al., 2022). Transgenerational fertility effects of endocrine disruption due to pesticide exposure were associated with greater mass, increased palmitoleic:palmitic acid ratio, and decreased glucose (Karlsson et al., 2021). As studies of multistressor deviations from normal metabolite activity continue, identification of meaningful pathway perturbations could provide a systematic method of identifying locations of environmental impacts without prerequisite knowledge of specific land use changes or fate and exposure of particular pollutants.

3.3 Sampling strategies

Traditional measures of contaminant impacts on amphibians focus on body burden concentrations and somatic indices of body condition (e.g., Băncilă et al., 2010) as well as general indicators of genotoxic and mutagenic impacts (i.e., comet and micronucleus assays) and more targeted analyses of cellular-level endpoints. For larger taxa (e.g., birds and mammals), marking individuals and collecting tissue samples can be a routine, noninvasive procedure conducted in the field to inform physiological condition. However, the small body size of amphibians hinders repeated sampling of blood or plasma for various analyses. For smaller amphibians, sampling techniques might involve toe-clipping for both individual identification as well as tissue collection, or sample collection might require sacrifice of individuals, particularly at earlier life stages.

Many common assays measure hematological enzymatic responses typical of exposure to specific xenobiotic contaminants (Ohmer et al., 2021). For example, esterase activity (acetylcholinesterase, butyryl cholinesterase, and carboxyl esterase) can indicate potential developmental effects, but response varies greatly within and between species (Venturino et al., 2003; Venturino and de D’Angelo, 2005). Hematological measures representative of immune response include leukocytes, neutrophil/lymphocyte ratio, bacterial killing assays, and delayed hypersensitivity assays (Ohmer et al., 2021). Changes in neutrophils and lymphocytes are often proportional with increased glucocorticoid levels, indicating physiological stress; however, neonicotinoid exposure altered leukocyte profiles relative to neutrophils or eosinophils but did not affect CORT levels in northern leopard frogs (Gavel et al., 2021). Neutrophil to lymphocyte ratios were also a good indication of environmental stressors and were associated with total dissolved solid levels in aquatic habitats that impacted growth and development (Ruso et al., 2021). Combinations of physiological indices are also informative to link endpoints with individual condition (e.g., Park et al., 2021).

Several minimally or non-invasive techniques that may be more effective for amphibians are now routinely used including urinalysis, fecal sampling, waterborne sampling, salivary swabs, and dermal swabs (Narayan et al., 2019). Among these techniques, saliva has only been validated in adults and not juveniles. Janin et al. (2012) concluded that CORT concentrations in saliva were highly correlated with urine measurements in toads. Urinalysis has been previously used to track endocrine response such as the reproductive hormones estradiol, progesterone, and testosterone within both captive and wild caught amphibians (Narayan, 2013). Additionally, CORT levels can be quantified from urine to identify differences in stress response due to predation risk or pathogen prevalence (Kindermann et al., 2012; Narayan et al., 2019). While urine samples have the distinct advantage of being highly concentrated for measuring endocrine functions along with physiological stress, the method is not always ideal for smaller amphibians producing lower volumes (Narayan, 2013; Baugh et al., 2018).

Another minimally invasive technique for endocrine analysis for amphibians of any size is immersing the individual in a clean container of water for a designated length of time, after which the released hormones can be quantified from the water. Most studies have used this technique to measure CORT release rates, which have been validated with circulating plasma levels (Gabor et al., 2013a). Waterborne sampling has evaluated CORT release rates in the presence or absence of Bd, predation, or pesticide exposure (Gabor et al., 2013b; Van Meter et al., 2019; Snyder et al., 2022) and can also be indicative of nitrate stress in amphibians (Ruthsatz et al., 2023). CORT is a potentially useful biomarker for amphibians to indicate stress, and non-invasive sampling methods offer the potential of serially sampling the same individual (Narayan et al., 2019; Tornabene et al., 2021). Therefore, within a short timeframe baseline values and acute elevation in CORT are quantifiable (Hammond et al., 2018). Szymanski et al. (2006) collected feces of adult anurans to quantify reproductive hormones, enabling sex identification. In addition to CORT, two other reproductive hormones, progesterone and estradiol, have also been quantified in water from amphibians (Baugh et al., 2018).

While dermal swabbing is most notable for detecting the presence of pathogens in amphibians (e.g., Standish et al., 2018), more recent studies have expanded on what can be tested in amphibian mucus, such as DNA collection and glucocorticoids (Poorten et al., 2017; Narayan et al., 2019; Santymire et al., 2021). A lab-based salamander study examined the difference between liver metabolomics and dermal swab metabolomics, to determine if similar pathways are impacted in the presence of pesticides, measuring for presence/absence of pathogens and glutathione as well (Van Meter et al., in press). Dermal swabs enable in situ field sampling with minimal handing time, which is particularly advantageous for threatened and endangered species and allows serial sampling of the same individual or environment (Santymire et al., 2018, 2021; Scheun et al., 2019; Tornabene et al., 2021).

Sampling techniques that are well-established at the individual level can also provide a comparison, through DGE or hepatic metabolites, of localized stressor response indicative of differential population-level effects. Such response-based metrics could preclude the necessity to anticipate synergies, compensation, and interactions between ubiquitous stressors that might be heterogeneously distributed. Rather than spatially explicit analysis of stressors within the species distribution, identifying variance in relative response within the population could target conservation concerns more rapidly within a diverse landscape while accounting for baseline fluxes in physiology. For example, the complex physiological changes during metamorphosis comprise endocrine interactions and changes in energy allocation as tail resorption and leg growth occurs, such that tissues sampled could vary in the cellular-level activity. Stage-specific fluxes in endocrine activity also affect response observed in individuals, such that standardizing measurements within consistent developmental stages is advisable when sampling larvae. Field-based observations could also be influenced by the environment of the individuals, making observations about water quality, larval density, and community composition relevant to evaluating stress response. As measures of individual response are considered within appropriate life stage time points, comparable evaluation of location-specific perturbations in baseline physiological functions can guide more targeted conservation actions.

4 Discussion

As amphibians are impacted by multiple stressors and their interactions, the capability to assess cumulative impacts on biochemical pathways within an organism’s native habitat facilitates quantification of exposure risk and possible additive or synergistic effects. However, even at the organismal level, amphibians often lack sufficient toxicology data for evaluation of cellular-level responses (Marlatt and Martyniuk, 2017), and variance in individual response and chronic exposure obscure definitive metrics of detrimental effects on the population. Additional research is needed to identify reliable biomarkers that are consistently indicative of points of departure from normal cellular function in response to environmental stressors. Standardized indices of biomarker perturbations in response to stressors enables identification of reliable predictors of long-term impacts (Pham and Sokolova, 2023); however, caution must be taken to verify cellular responses are linked to demographic effects, rather than a homeostatic response to stressors (Forbes et al., 2006). Evaluating potential threats linked to synergistic exposure effects (e.g., reduced dermal microbiota) or multiple exposure routes (i.e., aquatic and terrestrial) requires situation-specific ecological context. Implementation of weight of evidence effects could further classify cumulative threat levels of variable biomarker responses (Cecchetto et al., 2023). An initial step towards multi-stressor risk assessment is outlined here, namely by exploring stage-specific variance in biochemical pathways and identifying points of physiological vulnerability in the life cycle as a screening-level conservation approach.

Statements

Author contributions

JA: Writing – review & editing, Writing – original draft, Conceptualization. DG: Writing – review & editing. WH: Writing – review & editing. RM: Writing – review & editing. SP: Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Acknowledgments

Authors thank Jon Haselman and Sandy Raimondo as well as the two reviewers for suggestions on previous drafts.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AdelizziR.PortmannJ.Van MeterR. (2019). Effect of individual and combined treatments of pesticide, fertilizer, and salt on growth and corticosterone levels of larval southern leopard frogs (Lithobates sphenocephala). Arch. Environ. Contam Toxicol.77, 2939. doi: 10.1007/s00244-019-00629-6

  • 2

    AllenC.GonzalesR.ParrottL. (2020). Modelling the contribution of ephemeral wetlands to landscape connectivity. Ecol. Model.419, 108944.

  • 3

    AwkermanJ.RaimondoS.SchmolkeA.GalicN.Rueda-CedielP.KapoK.et al. (2020). Guidance for developing amphibian population models for ecological risk assessment. Integr. Environ. Assess. Manage.16, 223233. doi: 10.1002/ieam.4215

  • 4

    AwkermanJ. A.GreenbergC. H. (2022). Projected climate and hydroregime variability constrain ephemeral wetland-dependent amphibian populations in simulations of southern toads. Ecologies3, 235248. doi: 10.3390/ecologies3020018

  • 5

    BăncilăR. I.HartelT.PlăiaşuR.SmetsJ.CogălniceanuD. (2010). Comparing three body condition indices in amphibians: a case study of yellow-bellied toad Bombina variegata. Amphibia-Reptilia31, 558562.

  • 6

    BaranekC. T.MaynardC.McHenryC.ClulowJ.MahonyM. (2022). Identifying a limiting factor in the population dynamics of a threatened amphibian: The influence of extended female maturation on operational sex ratio. Austral Ecol.47, 239250.

  • 7

    BatesK. A.SommerU.HopkinsK. P.SheltonJ. M. G.WierzbickiC.SergeantC.et al. (2022). Microbiome function predicts amphibian chytridiomycosis disease dynamics. Microbiome10, 44. doi: 10.1186/s40168-021-01215-6

  • 8

    BaughA. T.BastienB.StillM. B.StowellN. (2018). Validation of water-borne steroid hormones in a tropical frog (Physalaemus pustulosus). Gen. Comp. Endocrinol.261, 6780.

  • 9

    BeckerC. G.GreenspanS. E.MartinsR. A.LyraM. L.PristP.MetzgerJ. P.et al. (2023). Habitat split as a driver of disease in amphibians. Biol. Rev.98, 727746. doi: 10.1111/brv.12927

  • 10

    BeldenL. K.WingfieldJ. C.KieseckerJ. M. (2010). Variation in the hormonal stress response among larvae of three amphibian species. J. Exp. Zoology Part A: Ecol. Genet. Physiol.313, 524531.

  • 11

    BienentreuJ. F.SchockD. M.GreerA. L.LesbarrèresD. (2022). Ranavirus amplification in low-diversity amphibian communities. Front. Vet. Sci.9. doi: 10.3389/fvets.2022.755426

  • 12

    BlausteinA. R.GervasiS. S.JohnsonP. T.HovermanJ. T.BeldenL. K.BradleyP. W.et al. (2012). Ecophysiology meets conservation: understanding the role of disease in amphibian population declines. Philos. Trans. R Soc. Lond B Biol. Sci.367, 16881707. doi: 10.1098/rstb.2012.0011

  • 13

    BlausteinA. R.WallsS. C.BancroftB. A.LawlerJ. J.SearleC. L.GervasiS. S. (2010). Direct and indirect effects of climate change on amphibian populations. Diversity2, 281313.

  • 14

    BögiC.LevyG.LutzI.KloasW. (2002). Functional genomics and sexual differentiation in amphibians. Comp. Biochem. Physiol. B Biochem. Mol. Biol.133, 559570. doi: 10.1016/s1096-4959(02)00162-8

  • 15

    BókonyV.KövérS.NemesháziE.LikerA.SzékelyT. (2017). Climate-driven shifts in adult sex ratios via sex reversals: the type of sex determination matters. Philos. Trans. R. Soc. B: Biol. Sci.372, 20160325.

  • 16

    BrannellyL. A.OhmerM. E. B.SaenzV.Richards-ZawackiC. L. (2019). Effects of hydroperiod on growth, development, survival and immune defences in a temperate amphibian. Funct. Ecol.33, 19521961. doi: 10.1111/1365-2435.13419

  • 17

    BretonB. A.BeatyL.BennettA. M.KyleC. J.LesbarrèresD.VilaçaS. T.et al. (2022). Testing the effectiveness of environmental DNA (eDNA) to quantify larval amphibian abundance. Environ. DNA4, 12291240.

  • 18

    BriceE. M.HalabiskyM.RayA. M. (2022). Making the leap from ponds to landscapes: Integrating field-based monitoring of amphibians and wetlands with satellite observations. Ecol. Indic.135, 108559.

  • 19

    BrooksG. C.KindsvaterH. K. (2022). Early development drives variation in amphibian vulnerability to global change. Front. Ecol. Evolution.10, 813414. doi: 10.3389/fevo.2022.813414

  • 20

    BrühlC. A.PieperS.WeberB. (2011). Amphibians at risk? Susceptibility of terrestrial amphibian life stages to pesticides. Environ. Toxicol. Chem.30, 24652472. doi: 10.1002/etc.650

  • 21

    BrühlC. A.SchmidtT.PieperS.AlscherA. (2013). Terrestrial pesticide exposure of amphibians: an underestimated cause of global decline? Sci. Rep.3, 1135. doi: 10.1038/srep01135

  • 22

    BrunnerJ. L.OlsonA. D.RiceJ. G.MeinersS. E.Le SageM. J.CundiffJ. A.et al. (2019). Ranavirus infection dynamics and shedding in American bullfrogs: consequences for spread and detection in trade. Dis. Aquat Organ135, 135150. doi: 10.3354/dao03387

  • 23

    BurracoP.Díaz-PaniaguaC.Gomez-MestreI. (2017). Different effects of accelerated development and enhanced growth on oxidative stress and telomere shortening in amphibian larvae. Sci. Rep.7, 7494. doi: 10.1038/s41598-017-07201-z

  • 24

    BurracoP.Gomez-MestreI. (2016). Physiological stress responses in amphibian larvae to multiple stressors reveal marked anthropogenic effects even below lethal levels. Physiol. Biochem. Zool.89, 462472. doi: 10.1086/688737

  • 25

    CalboliF. C.FisherM. C.GarnerT. W.JehleR. (2011). The need for jumpstarting amphibian genome projects. Trends Ecol. Evol.26, 378379.

  • 26

    CampbellL. J.HammondS. A.PriceS. J.SharmaM. D.GarnerT. W. J.BirolI.et al. (2018). A novel approach to wildlife transcriptomics provides evidence of disease-mediated differential expression and changes to the microbiome of amphibian populations. Mol. Ecol.27, 14131427. doi: 10.1111/mec.14528

  • 27

    CecchettoM.PeruzzaL.GiubilatoE.BernardiniI.RovereG. D.MarcominiA.et al. (2023). An innovative index to incorporate transcriptomic data into weight of evidence approaches for environmental risk assessment. Environ. Res.227, 115745. doi: 10.1016/j.envres.2023.115745

  • 28

    ChenQ.LiC.GongZ.ChanE. C. Y.SnyderS. A.LamS. H. (2017). Common deregulated gene expression profiles and morphological changes in developing zebrafish larvae exposed to environmental-relevant high to low concentrations of glucocorticoids. Chemosphere.172, 429439. doi: 10.1016/j.chemosphere.2017.01.036

  • 29

    CordierJ. M.AguilarR.LescanoJ. N.LeynaudG. C.BoninoA.MilochD.et al. (2021). A global assessment of amphibian and reptile responses to land-use changes. Biol. Conserv.253, 108863.

  • 30

    CostaA.DonderoL.AllariaG.Morales SanchezB. N.RosaG.SalvidioS.et al. (2021). Modelling the amphibian chytrid fungus spread by connectivity analysis: towards a national monitoring network in Italy. Biodivers Conserv.30, 28072825. doi: 10.1007/s10531-021-02224-5

  • 31

    CrespiE. J.DenverR. J. (2005). Roles of stress hormones in food intake regulation in anuran amphibians throughout the life cycle. Comp. Biochem. Physiol. A Mol. Integr. Physiol.141, 381390. doi: 10.1016/j.cbpb.2004.12.007

  • 32

    DenverR. J. (2009). Stress hormones mediate environment-genotype interactions during amphibian development. Gen. Comp. Endocrinol.164, 2031. doi: 10.1016/j.ygcen.2009.04.016

  • 33

    Di LorenzoM.BarraT.RosatiL.ValianteS.CapaldoA.De FalcoM.et al. (2020). Adrenal gland response to endocrine disrupting chemicals in fishes, amphibians and reptiles: A comparative overview. Gen. Comp. Endocrinology.297, 113550.

  • 34

    DornellesM. F.OliveriaG. T. (2016). Toxicity of atrazine, glyphosate, and quinclorac in bullfrog tadpoles exposed to concentrations below legal limits. Environ. Sci. pollut. Res.23, 16101620.

  • 35

    DuarteH.TejedoM.KatzenbergerM.MarangoniF.BaldoD.BeltránJ. F.et al. (2012). Can amphibians take the heat? Vulnerability to climate warming in subtropical and temperate larval amphibian communities. Global Change Biol.18, 412421.

  • 36

    Duarte-GutermanP.Navarro-MartínL.TrudeauV. L. (2014). Mechanisms of crosstalk between endocrine systems: regulation of sex steroid hormone synthesis and action by thyroid hormones. Gen. Comp. Endocrinol.203, 6985. doi: 10.1016/j.ygcen.2014.03.015

  • 37

    Egea-SerranoA.RelyeaR.TejedoM.TorralvaM. (2012). Understanding of the impact of chemicals on amphibians: a meta-analytic review. Ecol. Evol.2, 1382-1397. doi: 10.1002/ece3.249

  • 38

    EllisonA.ZamudioK.LipsK.Muletz-WolzC. (2020). Temperature-mediated shifts in salamander transcriptomic responses to the amphibian-killing fungus. Mol. Ecol.29, 325343.

  • 39

    EskewE. A.ShockB. C.LaDouceurE. E. B.KeelK.MillerM. R.FoleyJ. E.et al. (2018). Gene expression differs in susceptible and resistant amphibians exposed to Batrachochytrium dendrobatidis. R Soc. Open Sci.5, 170910. doi: 10.1098/rsos.170910

  • 40

    EwaldJ. D.BasuN.CrumpD.BoulangerE.HeadJ. (2022). Characterizing variability and uncertainty associated with transcriptomic dose-response modeling. Environ. Sci. Technol.56, 1596015968. doi: 10.1021/acs.est.2c04665

  • 41

    FisherM. C.GarnerT. W. J. (2020). Chytrid fungi and global amphibian declines. Nat. Rev. Microbiol.18, 332343. doi: 10.1038/s41579-020-0335-x

  • 42

    FodenW. B.ButchartS. H.StuartS. N.ViéJ. C.AkçakayaH. R.AnguloA.et al. (2013). Identifying the world’s most climate change vulnerable species: a systematic trait-based assessment of all birds, amphibians and corals. PloS One8, e65427.

  • 43

    FoglianoC.MottaC. M.AcloqueH.AvalloneB.CarotenutoR. (2023). Water contamination by delorazepam induces epigenetic defects in the embryos of the clawed frog Xenopus laevis. Sci. Total Environ.896, 165300.

  • 44

    ForbesV. E.PalmqvistA.BachL. (2006). The use and misuse of biomarkers in ecotoxicology. Environ. Toxicol. Chem.25, 272280. doi: 10.1897/05-257r.1

  • 45

    ForesterB. R.MurphyM.MellisonC.PetersenJ.PilliodD. S.Van HorneR.et al. (2022). Genomics-informed delineation of conservation units in a desert amphibian. Mol. Ecol.31, 52495269. doi: 10.1111/mec.16660

  • 46

    FrishkoffL. O.HadlyE. A.DailyG. C. (2015). Thermal niche predicts tolerance to habitat conversion in tropical amphibians and reptiles. Global Change Biol.21, 39013916.

  • 47

    FunkW. C.ZamudioK. R.CrawfordA. J. (2021). Advancing understanding of amphibian evolution, ecology, behavior, and conservation with massively parallel sequencing. Population genomics: wildlife.2021, 211254.

  • 48

    GaborC. R.BoschJ.FriesJ. N.DavisD. R. (2013a). A non-invasive water-borne hormone assay for amphibians. Amphibia-Reptilia34, 151162.

  • 49

    GaborC. R.FisherM. C.BoschJ. (2013b). A non-invasive stress assay shows that tadpole populations infected with Batrachochytrium dendrobatidis have elevated corticosterone levels. PloS One8, e56054.

  • 50

    GavelM. J.YoungS. D.DaltonR. L.SoosC.McPheeL.ForbesM. R.et al. (2021). Effects of two pesticides on northern leopard frog (Lithobates pipiens) stress metrics: Blood cell profiles and corticosterone concentrations. Aquat. Toxicol.235, 105820.

  • 51

    GlinskiD. A.Van MeterR. J.PuruckerS. T.HendersonW. M. (2021). Route of exposure influences pesticide body burden and the hepatic metabolome in post-metamorphic leopard frogs. Sci. Total Environ.779, 146358. doi: 10.1016/j.scitotenv.2021.146358

  • 52

    GrantE. H. C.MillerD. A. W.SchmidtB. R.AdamsM. J.AmburgeyS. M.ChambertT.et al. (2016). Quantitative evidence for the effects of multiple drivers on continental-scale amphibian declines. Sci. Rep.6, 25625. doi: 10.1038/srep25625

  • 53

    GreenD. M.LannooM. J.LesbarrèresD.MuthsE. (2020). Amphibian population declines: 30 years of progress in confronting a complex problem. Herpetologica.76, 97100.

  • 54

    GreenbergC. H.ZarnochS. J.AustinJ. D. (2017). Weather, hydroregime, and breeding effort influence juvenile recruitment of anurans: implications for climate change. Ecosphere8, e01789. doi: 10.1002/ecs2.1789

  • 55

    GuoW.YangY.ZhouX.MingR.HuD.LuP. (2022). Insight into the toxic effects, bioconcentration and oxidative stress of acetamiprid on Rana nigromaculata tadpoles. Chemosphere305, 135380.

  • 56

    HaddadN. M.BrudvigL. A.ClobertJ.DaviesK. F.GonzalezA.HoltR. D.et al. (2015). Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci. Adv.1, e1500052.

  • 57

    HaislipN. A.GrayM. J.HovermanJ. T.MillerD. L. (2011). Development and disease: how susceptibility to an emerging pathogen changes through anuran development. PloS One6, e22307.

  • 58

    HammondT. T.AuZ. A.HartmanA. C.Richards-ZawackiC. L. (2018). Assay validation and interspecific comparison of salivary glucocorticoids in three amphibian species. Conserv. Physiol.6, coy055.

  • 59

    HarperE. B.SemlitschR. D. (2007). Density dependence in the terrestrial life history stage of two anurans. Oecologia153, 879889. doi: 10.1007/s00442-007-0796-x

  • 60

    HayesT. B.KhouryV.NarayanA.NazirM.ParkA.BrownT.et al. (2010). Atrazine induces complete feminization and chemical castration in male African clawed frogs (Xenopus laevis). Proc. Natl. Acad. Sci. U S A.107, 46124617. doi: 10.1073/pnas.0909519107

  • 61

    HelbingC. C. (2012). The metamorphosis of amphibian toxicogenomics. Front. Genet.14. doi: 10.3389/fgene.2012.00037

  • 62

    HoffmannE. P.CavanoughK. L.MitchellN. J. (2021). Low desiccation and thermal tolerance constrains a terrestrial amphibian to a rare and disappearing microclimate niche. Conserv. Physiol.9, coab027.

  • 63

    HossackB. R.OjaE. B.OwensA. K.HallD.CobosC.CrawfordC. L.et al. (2023). Empirical evidence for effects of invasive American Bullfrogs on occurrence of native amphibians and emerging pathogens. Ecol. Appl.33(2), e2785. doi: 10.1002/eap.2785

  • 64

    HotalingS.KelleyJ. L.FrandsenP. B. (2021). Toward a genome sequence for every animal: Where are we now? Proc. Natl. Acad. Sci.118, e2109019118.

  • 65

    HovermanJ. T.GrayM. J.HaislipN. A.MillerD. L. (2011). Phylogeny, life history, and ecology contribute to differences in amphibian susceptibility to ranaviruses. EcoHealth8, 301319.

  • 66

    HuangL.LiJ.AnboukariaH.LuoZ.ZhaoM.WuH. (2016). Comparative transcriptome analyses of seven anurans reveal functions and adaptations of amphibian skin. Sci. Rep.6, 24069.

  • 67

    IchuT. A.HanJ.BorchersC. H.LesperanceM.HelbingC. C. (2014). Metabolomic insights into system-wide coordination of vertebrate metamorphosis. BMC Dev. Biol.14, 5. doi: 10.1186/1471-213X-14-5

  • 68

    JaninA.LÉNAJ. P.DebloisS.JolyP. (2012). Use of stress-hormone levels and habitat selection to assess functional connectivity of a landscape for an Amphibian. Conserv. Biol.26, 923931.

  • 69

    JiménezR. R.AlvaradoG.RuepertC.BallesteroE.SommerS. (2021). The fungicide chlorothalonil changes the amphibian skin microbiome: A potential factor disrupting a host disease-protective trait. Appl. Microbiol.1, 2637. doi: 10.3390/applmicrobiol1010004

  • 70

    JonssonB.JonssonN.HansenM. M. (2022). Knock-on effects of environmental influences during embryonic development of ectothermic vertebrates. Q. Rev. Biol.97, 95139.

  • 71

    KarlssonO.SvanholmS.ErikssonA.ChidiacJ.ErikssonJ.JernerénF.et al. (2021). Pesticide-induced multigenerational effects on amphibian reproduction and metabolism. Sci. Total Environ.775, 145771. doi: 10.1016/j.scitotenv.2021.145771

  • 72

    KerbyJ. L.HartA. J.StorferA. (2011). Combined effects of virus, pesticide and predator cue on the larval tiger salamander (Ambystoma tigrinum). EcoHealth8, 46–54. doi: 10.1007/s10393-011-0682-1

  • 73

    KieseckerJ. M. (2002). Synergism between trematode infection and pesticide exposure: a link to amphibian limb deformities in nature? Proc. Natl. Acad. Sci.99, 99009904.

  • 74

    KieseckerJ. M. (2011). Global stressors and the global decline of amphibians: tipping the stress immunocompetency axis. Ecol. Res.26, 897908. doi: 10.1007/s11284-010-0702-6

  • 75

    KindermannC.NarayanE. J.HeroJ. M. (2012). Urinary corticosterone metabolites and chytridiomycosis disease prevalence in a free-living population of male Stony Creek frogs (Litoria wilcoxii). Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol.162, 171176.

  • 76

    KisselA. M.TenanS.MuthsE. (2020). Density dependence and adult survival drive dynamics in two high elevation amphibian populations. Diversity12, 478.

  • 77

    KnutieS. A.WilkinsonC. L.KohlK. D.RohrJ. R. (2017). Early-life disruption of amphibian microbiota decreases later-life resistance to parasites. Nat. Commun.8, 86. doi: 10.1038/s41467-017-00119-0

  • 78

    KohliA. K.LindauerA. L.BrannellyL. A.OhmerM. E. B.Richards-ZawackiC.Rollins-SmithL.et al. (2019). Disease and the drying pond: Examining possible links among drought, immune function, and disease development in amphibians. Physiol. Biochem. Zool.92, 339348. doi: 10.1086/703137

  • 79

    KoschT. A.WaddleA. W.CooperC. A.ZengerK. R.GarrickD. J.BergerL.et al. (2023). Genetic approaches for increasing fitness in endangered species. Trends Ecol. Evolution.37, 332345.

  • 80

    KrynakK. L.BurkeD. J.BenardM. F. (2017). Rodeo™ herbicide negatively affects Blancha’d’s cricket frogs (Acris blanchardi) survival and alters the skin-associated bacterial community. J. Herpetology51, 402410.

  • 81

    LaLoneC. A.BlatzD. J.JensenM. A.VlietS. M.MayasichS.MattinglyK. Z.et al. (2023). From protein sequence to structure: The next frontier in cross-species extrapolation for chemical safety evaluations. Environ. Toxicol. Chem.42, 463474.

  • 82

    LambertM. R.GillerG. S.BarberL. B.FitzgeraldK. C.SkellyD. K. (2015). Suburbanization, estrogen contamination, and sex ratio in wild amphibian populations. Proc. Natl. Acad. Sci.112, 1188111886.

  • 83

    Lertzman-LepofskyG. F.KisselA. M.SinervoB.PalenW. J. (2020). Water loss and temperature interact to compound amphibian vulnerability to climate change. Global Change Biol.26, 48684879.

  • 84

    Le SageE. H.OhmerM. E. B.LaBumbardB. C.AltmanK. A.ReinertL. K.BednarkJ. G.et al. (2022). Localized carry-Over effects of pond drying on survival, growth, and pathogen defenses in amphibians. Ecosphere13, e4224. doi: 10.1002/ecs2.4224

  • 85

    LiB.MaY.ZhangY. H. (2017). Oxidative stress and hepatotoxicity in the frog, Rana chensinensis, when exposed to low doses of trichlorfon. J. Environ. Sci. Health Part B.52, 476482.

  • 86

    LiZ.WangQ.SunK.FengJ. (2021). Prevalence of Batrachochytrium dendrobatidis in Amphibians From 2000 to 2021: A global systematic review and meta-analysis. Front. Vet. Sci.8. doi: 10.3389/fvets.2021.791237

  • 87

    LiuN.MengB.BianW.YangM.ShuL.LiuY.et al. (2022). The beneficial roles of poisonous skin secretions in survival strategies of the odorous frog Odorrana andersonii. Sci. Nat.109, 4. doi: 10.1007/s00114-021-01776-4

  • 88

    MarcoglieseD.KingK.BatesK. (2021). Effects of multiple stressors on northern leopard frogs in agricultural wetlands. Parasitology148, 827834. doi: 10.1017/S003118202100038X

  • 89

    MarlattV. L.BayenS.Castaneda-CortèsD.DelbèsG.GrigorovaP.LangloisV. S.et al. (2022). Impacts of endocrine disrupting chemicals on reproduction in wildlife and humans. Environ. Res.208, 112584.

  • 90

    MarlattV. L.MartyniukC. J. (2017). Biological responses to phenylurea herbicides in fish and amphibians: New directions for characterizing mechanisms of toxicity. Comp. Biochem. Physiol. Part C: Toxicol. Pharmacol.194, 921. doi: 10.1016/j.cbpc.2017.01.002

  • 91

    MazerolleM. J.BaileyL. L.KendallW. L.RoyleJ. A.ConverseS. J.NicholsJ. D. (2007). Making great leaps forward: accounting for detectability in herpetological field studies. J. Herpetology.41, 672689.

  • 92

    McArdleM. E.FreemanE. L.StaveleyJ. P.OrtegoL. S.CoadyK. K.WeltjeL.et al. (2020). Critical review of read-across potential in testing for endocrine-related effects in vertebrate ecological receptors. Environ. Toxicol. Chem.39, 739753.

  • 93

    McCartney-MelstadE.VuJ. K.ShafferH. B. (2018). Genomic data from an endangered amphibian reveal unforeseen consequences of fragmentation by roads. bioRxiv27(22):4430–4443. doi: 10.1101/306340

  • 94

    McCoyK. A.PeraltaA. L. (2018). Pesticides could alter amphibian skin microbiomes and the effects of Batrachochytrium dendrobatidis. Front. Microbiol.9, 748.

  • 95

    MelvinS. D.LeuschF. D. L.CarrollA. R. (2018). Metabolite profiles of striped marsh frog (Limnodynastes peronii) larvae exposed to the anti-androgenic fungicides vinclozolin and propiconazole are consistent with altered steroidogenesis and oxidative stress. Aquat Toxicol.199, 232239. doi: 10.1016/j.aquatox.2018.04.004

  • 96

    MenonJ.RozmanR. (2007). Oxidative stress, tissue remodeling and regression during amphibian metamorphosis. Comp. Biochem. Physiol. Part C: Toxicol. Pharmacol.145, 625631.

  • 97

    MettsB. S.BuhlmannK. A.TubervilleT. D.ScottD. E.HopkinsW. A. (2013). Maternal transfer of contaminants and reduced reproductive success of southern toads (Bufo [Anaxyrus] terrestris) exposed to coal combustion waste. Environ. Sci. Technol.47, 28462853.

  • 98

    MittalK.EwaldJ.BasuN. (2022). Transcriptomic points of departure calculated from rainbow trout gill, liver, and gut cell lines exposed to methylmercury and fluoxetine. Environ. Toxicol. Chem.41, 19821992. doi: 10.1002/etc.5395

  • 99

    MiuraI. (2017). Sex determination and sex chromosomes in Amphibia. Sexual Dev.11, 298306.

  • 100

    MossW. E.Mc-Devitt-GallesT.MuthsE.BobzienS.PurificatoJ.JohnsonP. T. J. (2021). Resilience of native amphibian communities following catastrophic drought: Evidence from a decade of regional-scale monitoring. Biol. Cons263, 109352. doi: 10.1016/j.biocon.2021.109352

  • 101

    MunstermannM. J.HeimN. A.McCauleyD. J.PayneJ. L.UphamN. S.WangS. C.et al. (2022). A global ecological signal of extinction risk in terrestrial vertebrates. Conserv. Biol.36, e13852.

  • 102

    NarayanE. J. (2013). Non-invasive reproductive and stress endocrinology in amphibian conservation physiology. Conserv. Physiol.1, cot011.

  • 103

    NarayanE. J.CockremJ. F.HeroJ.-M. (2013). Sight of a predator induces a corticosterone stress response and generates fear in an amphibian. PloS One8, e73564. doi: 10.1371/journal.pone.0073564

  • 104

    NarayanE. J.ForsburgZ. R.DavisD. R.GaborC. R. (2019). Non-invasive methods for measuring and monitoring stress physiology in imperiled amphibians. Front. Ecol. Evol.7, 431.

  • 105

    NeelyW. J.GreenspanS. E.StahlL. M.HeraghtyS. D.MarshallV. M.AtkinsonC. L.et al. (2022). Habitat disturbance linked with host microbiome dispersion and Bd dynamics in temperate amphibians. Microbial Ecol.84, 901910.

  • 106

    NewboldT.HudsonL. N.HillS. L.ContuS.LysenkoI.SeniorR. A.et al. (2015). Global effects of land use on local terrestrial biodiversity. Nature520, 4550.

  • 107

    NolanN.HaywardM. W.Klop-TokerK.MahonyM.LemckertF.CallenA. (2023). Complex organisms must deal with complex threats: How does amphibian conservation deal with biphasic life cycles? Anim. (Basel)13, 1634. doi: 10.3390/ani13101634

  • 108

    OhbaS. (2011). Density-dependent effects of amphibian prey on the growth and survival of an endangered giant water bug. Insects2, 435446.

  • 109

    OhmerM. E.AltonL.CrampR. L. (2021). “Physiology provides a window into how the multi-stressor environment contributes to amphibian declines,” in Conservation physiology: applications for wildlife conservation and management (Oxford, United Kingdom: Oxford University Press), 165182. doi: 10.1093/oso/9780198843610.003.0010

  • 110

    ParkJ.-K.DoY. (2023). Current state of conservation physiology for amphibians: Major research topics and physiological parameters. Animals13, 3162. doi: 10.3390/ani13203162

  • 111

    ParkJ. K.KimJ. B.DoY. (2021). Reference intervals in combined veterinary clinical examinations of male black-spotted pond frogs (Pelophylax nigromaculatus). Anim. (Basel)11, 1407. doi: 10.3390/ani11051407

  • 112

    PaulB.SternerZ. R.BuchholzD. R.ShiY.-B.SachsL. M. (2022). Thyroid and corticosteroid signaling in amphibian metamorphosis. Cells11, 1595. doi: 10.3390/cells11101595

  • 113

    PechmannJ. H. K.ScottD. E.WhitfieldG. J.SemlitschR. D. (1989). Influence of wetland hydroperiod on diversity and abundance of metamorphosing juvenile amphibians. Wetl Ecol. Manag1, 311.

  • 114

    PetrovanS. O.SchmidtB. R. (2019). Neglected juveniles; a call for integrating all amphibian life stages in assessments of mitigation success (and how to do it). Biol. Conserv.236, 252260.

  • 115

    PetrovićT. G.KijanovićA.Kolarov TomaševićN.GavrićJ. P.DespotovićS. G.GavrilovićB. R.et al. (2021). Effects of desiccation on metamorphic climax in Bombina variegata: Changes in levels and patterns of oxidative stress parameters. Anim. (Basel)11, 953. doi: 10.3390/ani11040953

  • 116

    PhamD. N.SokolovaI. M. (2023). Dissecting integrated indices of multiple biomarker responses: Think before use. Integr. Environ. Assess. Manage.19, 302311. doi: 10.1002/ieam.4676

  • 117

    PochiniK. M.HovermanJ. T. (2017). Reciprocal effects of pesticides and pathogens on amphibian hosts: The importance of exposure order and timing. Environ. Poll221, 359366.

  • 118

    PoortenT. J.KnappR. A.RosenblumE. B. (2017). Population genetic structure of the endangered Sierra Nevada yellow-legged frog (Rana sierrae) in Yosemite National Park based on multi-locus nuclear data from swab samples. Conserv. Genet.18, 731744.

  • 119

    ProkićM. D.PetrovićT. G.GavrilovićB. R.DespotovićS. G.GavrićJ. P.KijanovićA.et al. (2021). Carry-over effects of desiccation stress on the oxidative status of fasting anuran juveniles. Front. Physiol.12, 2206.

  • 120

    RelyeaR. A. (2009). A cocktail of contaminants: how mixtures of pesticides at low concentrations affect aquatic communities. Oecologia159, 363376. doi: 10.1007/s00442-008-1213-9

  • 121

    RelyeaR. A.DiecksN. (2008). An unforeseen chain of events: lethal effects of pesticides on frogs at sublethal concentrations. Ecol. Appl.18, 17281742.

  • 122

    RelyeaR.HovermanJ. (2006). Assessing the ecology in ecotoxicology: a review and synthesis in freshwater systems. Ecol. Lett.9, 11571171.

  • 123

    RittenhouseT. A.SemlitschR. D. (2007). Distribution of amphibians in terrestrial habitat surrounding wetlands. Wetlands27, 153161.

  • 124

    RocoÁ.S.Ruiz-GarcíaA.BullejosM. (2021). Testis development and differentiation in amphibians. Genes12, 578. doi: 10.3390/genes12040578

  • 125

    RodriguezK. M.VoylesJ. (2020). The amphibian complement system and chytridiomycosis. J. Exp. Zool.333, 706719. doi: 10.1002/jez.2419

  • 126

    RohrJ. R.McCoyK. A. (2010). A qualitative meta-analysis reveals consistent effects of atrazine on freshwater fish and amphibians. Environ. Health Perspect.118, 2032.

  • 127

    RohrJ. R.RaffelT. R.HalsteadN. T.McMahonT. A.JohnsonS. A.BoughtonR. K.et al. (2013). Early-life exposure to a herbicide has enduring effects on pathogen-induced mortality. Proc. Biol. Sci.280, 20131502. doi: 10.1098/rspb.2013.1502

  • 128

    RohrJ. R.SagerT.SesterhennT. M.PalmerB. D. (2006). Exposure, postexposure, and density-mediated effects of atrazine on amphibians: breaking down net effects into their parts. Environ. Health Perspect.114, 4650. doi: 10.1289/ehp.8405

  • 129

    Rollins-SmithL. A. (2017). Amphibian immunity–stress, disease, and climate change. Developmental & Comparative Immunology, 66, 111-119.

  • 130

    RowJ. R.DonaldsonM. E.LonghiJ. N.SavilleB. J.MurrayD. L. (2016). Tissue-specific transcriptome characterization for developing tadpoles of the northern leopard frog (Lithobates pipiens). Genomics.108, 232240. doi: 10.1016/j.ygeno.2016.10.002

  • 131

    RowlandF. E.MuthsE.Eagles-SmithC. A.StrickerC. A.KrausJ. M.HarringtonR. A.et al. (2023). Complex life histories alter patterns of mercury exposure and accumulation in a pond-breeding amphibian. Environ. Sci. Technol57(10):4133–4142. doi: 10.1021/acs.est.2c04896

  • 132

    RuggeriJ.RibeiroL. P.PontesM. R.ToffoloC.CandidoM.CarrieroM. M.et al. (2019). Discovery of wild amphibians infected with ranavirus in Brazil. J. Wildl Dis.55, 897902.

  • 133

    Ruiz-GarcíaA.RocoÁ.S.BullejosM. (2021). Sex differentiation in amphibians: Effect of temperature and its influence on sex reversal. Sexual Dev.15, 157167.

  • 134

    RusoG. E.HoganN. S.SheedyC.GallantM. J.JardineT. D. (2021). Effects of agricultural stressors on growth and an immune status indicator in wood frog (Lithobates sylvaticus) tadpoles and metamorphs. Environ. Toxicol. Chem.40, 22692281. doi: 10.1002/etc.5107

  • 135

    RuthsatzK.EterovickP. C.BartelsF.MausbachJ. (2023). Contributions of water-borne corticosterone as one non-invasive biomarker in assessing nitrate pollution stress in tadpoles of Rana temporaria. Gen. Comp. Endocrinol.331, 114164. doi: 10.1016/j.ygcen.2022.114164

  • 136

    RuthsatzK.LyraM. L.LambertiniC.BelasenA. M.JenkinsonT. S.da Silva LeiteD.et al. (2020). Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil’s Atlantic Forest treefrogs. Sci. Rep.10, 22311. doi: 10.1038/s41598-020-79130-3

  • 137

    SachsL. M.BuchholzD. R. (2019). Insufficiency of thyroid hormone in frog metamorphosis and the role of glucocorticoids. Front. Endocrinol. (Lausanne)10. doi: 10.3389/fendo.2019.00287

  • 138

    SaeedM.RaisM.AkramA.WilliamsM. R.KellnerK. F.HashshamS. A.et al. (2022). Development and validation of an eDNA protocol for monitoring endemic Asian spiny frogs in the Himalayan region of Pakistan. Sci. Rep.12, 5624.

  • 139

    SantymireR. M.ManjerovicM. B.Sacerdote-VelatA. (2018). A novel method for the measurement of glucocorticoids in dermal secretions of amphibians. Conserv. Physiol.6, coy008. doi: 10.1093/conphys/coy008

  • 140

    SantymireR. M.Sacerdote-VelatA. B.GygliA.KeinathD. A.PooS.HinksonK. M.et al. (2021). Using dermal glucocorticoids to determine the effects of disease and environment on the critically endangered Wyoming toad. Conserv. Physiol.9, coab093. doi: 10.1093/conphys/coab093

  • 141

    SaucedoB.SerranoJ. M.Jacinto-MaldonadoM.LeuvenR. S. E. W.Rocha GarcíaA. A.Méndez BernalA.et al. (2019). Pathogen risk analysis for wild amphibian populations following the first report of a ranavirus outbreak in farmed American Bullfrogs (Lithobates catesbeianus) from Northern Mexico. Viruses.11, 26. doi: 10.3390/v11010026

  • 142

    ScheunJ.GreeffD.MedgerK.GanswindtA. (2019). Validating the use of dermal secretion as a matrix for monitoring glucocorticoid concentrations in African amphibian species. Conserv. Physiol.7, coz022. doi: 10.1093/conphys/coz022

  • 143

    SchirmerK.FischerB. B.MadureiraD. J.PillaiS. (2010). Transcriptomics in ecotoxicology. Analytical bioanalytical Chem.397, 917923.

  • 144

    ScottD. E.CaseyE. D.DonovanM. F.LynchT. K. (2007). Amphibian lipid levels at metamorphosis corelate to post-metamorphic terrestrial survival. Oecologia153, 521532.

  • 145

    SeimR. F.GlinskiD. A.LavelleC. M.AwkermanJ. A.HemmerB. L.HarrisP.et al. (2022). Using metabolomic profiling to inform use of surrogate species in ecological risk assessment practices. Comp. Biochem. Physiol. Part D Genomics Proteomics.41, 100947. doi: 10.1016/j.cbd.2021.100947

  • 146

    ShidemantleG.BussN.HuaJ. (2022). Are glucocorticoids good indicators of disturbance across populations that exhibit cryptic variation in contaminant tolerance? Anim. Conserv.25, 273284. doi: 10.1111/acv.12737

  • 147

    SmallingK. L.Eagles-SmithC. A.KatzR. A.GrantE. H. (2019). Managing the trifecta of disease, climate, and contaminants: Searching for robust choices under multiple sources of uncertainty. Biol. Conserv.236, 153161.

  • 148

    SmallingK. L.ReevesR.MuthsE.VandeverM.BattaglinW. A.HladikM. L.et al. (2015). Pesticide concentrations in frog tissue and wetland habitats in a landscape dominated by agriculture. Sci. Total Environ.502, 8090. doi: 10.1016/j.scitotenv.2014.08.114

  • 149

    SnyderM. N.HendersonW. M.GlinskiD. A.PuruckerS. T. (2022). Differentiating metabolomic responses amphibians to multiple stressors. Sci. Total Environment.838, 155666. doi: 10.1016/j.scitotenv.2022.155666

  • 150

    SodhiN. S.BickfordD.DiesmosA. C.LeeT. M.KohL. P.BrookB. W.et al. (2008). Measuring the meltdown: drivers of global amphibian extinction and decline. PloS One3, e1636. doi: 10.1371/journal.pone.0001636

  • 151

    Spitzen-van der SluijsA.CanessaS.MartelA.PasmansF. (2017). Fragile coexistence of a global chytrid pathogen with amphibian populations is mediated by environment and demography. Philos. Trans. R Soc. B.284, 20171444. doi: 10.1098/rspb.2017.1444

  • 152

    StandishI.LeisE.SchmitzN.CredicoJ.EricksonS.BaileyJ.et al. (2018). Optimizing, validating, and field testing a multiplex qPCR for the detection of amphibian pathogens. Dis. Aquat Organ129, 113. doi: 10.3354/dao03230

  • 153

    StevensV. M.BaguetteM. (2008). Importance of habitat quality and landscape connectivity for the persistence of endangered natterjack toads. Conserv. Biol.22, 11941204.

  • 154

    StuartS. N.ChansonJ. S.CoxN. A.YoungB. E.RodriguesA. S. L.FischmanD. L.et al. (2004). Status and trends of amphibian declines and extinctions worldwide. Science306, 17831786.

  • 155

    SzymanskiD. C.GistD. H.RothT. L. (2006). Anuran gender identification by fecal steroid analysis. Zoo Biology: Published affiliation Am. Zoo Aquarium Assoc.25, 3546.

  • 156

    TaylorB. E.ScottD. E.GibbonsJ. W. (2006). Catastrophic reproductive failure, terrestrial survival, and persistence of the marbled salamander. Conserv. Biol.20, 792801.

  • 157

    ThambirajahA. A.KoideE. M.ImberyJ. J.HelbingC. C. (2019). Contaminant and environmental influences on thyroid hormone action in amphibian metamorphosis. Front. Endocrinol. (Lausanne).10, 405.

  • 158

    ThambirajahA. A.WadeM. G.VerreaultJ.BuisineN.AlvesV. A.LangloisV. S.et al. (2022). Disruption by stealth – Interference of endocrine disrupting chemicals on hormonal crosstalk with thyroid axis function in humans and other animals. Environ. Res.203, 111906. doi: 10.1016/j.envres.2021.111906

  • 159

    ThompsonC. M.SweeneyM. R.PopescuV. D. (2022). Carryover effects of pesticide exposure and pond drying on performance, behavior, and sex ratios in a pool breeding amphibian. J. Zoology.317, 229240.

  • 160

    ToddB. D.BergeronC. M.HepnerM. J.HopkinsW. A. (2011). Aquatic and terrestrial stressors in amphibians: a test of the double jeopardy hypothesis based on maternally and trophically derived contaminants. Environ. Toxicol. Chem.30, 22772284.

  • 161

    TornabeneB. J.HossackB. R.CrespiE. J.BreunerC. W. (2021). Evaluating corticosterone as a biomarker for amphibians exposed to increased salinity and ambient corticosterone. Conserv. Physiol.9, coab049. doi: 10.1093/conphys/coab049

  • 162

    TrudeauV. L.ThomsonP.ZhangW. S.ReynaudS.Navarro-MartinL.LangloisV. S. (2020). Agrochemicals disrupt multiple endocrine axes in amphibians. Mol. Cell Endocrinol.513, 110861. doi: 10.1016/j.mce.2020.110861

  • 163

    UsalM.VeyrencS.DarracqM.RegnaultC.SrodaS.FiniJ. B. (2021). Transgenerational metabolic disorders and reproduction defects induced by benzo [a] pyrene in Xenopus tropicalis. Environmental Pollution, 269, 116109.

  • 164

    Van MeterR. J.AdelizziR.GlinskiD. A.HendersonW. M. (2019). Agrochemical mixtures and amphibians: the combined effects of pesticides and fertilizer on stress, acetylcholinesterase activity, and bioaccumulation in a terrestrial environment. Environ. Toxicol. Chem.38, 10521061.

  • 165

    Van MeterR. J.GlinskiD. A.PuruckerS. T.HendersonW. M. (2018). Influence of exposure to pesticide mixtures on the metabolomic profile in post-metamorphic green frogs (Lithobates clamitans). Sci. Total Environ.624, 13481359. doi: 10.1016/j.scitotenv.2017.12.175

  • 166

    Van MeterR. J.GlinskiD. A.PuruckerS. T.HendersonW. M. (2022). Induced hepatic glutathione and metabolomic alterations following mixed pesticide and fertilizer exposures in juvenile leopard frogs (Lithobates sphenocephala). Environ. Toxicol. Chem.41, 122133. doi: 10.1002/etc.5245

  • 167

    Van MeterR. J.GlinskiD. A.WanatJ. J.PuruckerS. T.HendersonW. M.Validation of salamander dermal mucus swabs as a novel, non-lethal approach for amphibian metabolomics, glutathione and disease analysis following pesticide exposure. Environ. Toxicol. Chem. (In press).

  • 168

    VastagL.JorgensenP.PeshkinL.WeiR.RabinowitzJ. D.KirschnerM. W. (2011). Remodeling of the metabolome during early frog development. PloS One6, e16881. doi: 10.1371/journal.pone.0016881

  • 169

    VenturinoA.de D’AngeloA. M. P. (2005). Biochemical targets of xenobiotics: Biomarkers in amphibian ecotoxicology. Appl. Herpetology2, 335353.

  • 170

    VenturinoA.RosenbaumE.Caballero de CastroA.AnguianoO. L.GaunaL.Fonovich de SchroederT.et al. (2003). Biomarkers of effect in toads and frogs. Biomarkers8, 167186. doi: 10.1080/1354700031000120116

  • 171

    VijayV.ArmsworthP. R. (2021). Pervasive cropland in protected areas highlight trade-offs between conservation and food security. Proc. Natl. Acad. Sci.118, e2010121118.

  • 172

    WakeD. B.VredenburgV. T. (2008). Are we in the midst of the sixth mass extinction? A view from the world of amphibians. Proc. Natl. Acad. Sci.105, 1146611473.

  • 173

    WallsS. C.BarichivichW. J.BrownM. E. (2013). Drought, deluge and declines: the impact of precipitation extremes on amphibians in a changing climate. Biol. (Basel)2, 399418. doi: 10.3390/biology2010399

  • 174

    WallsS. C.GaborC. R. (2019). Integrating behavior and physiology into strategies for amphibian conservation. Front. Ecol. Evol.7, 234.

  • 175

    WangS.FengC.YinS.FengZ.TangJ.LiuN.et al. (2020). A novel peptide from the skin of amphibian Rana limnocharis with potency to promote skin wound repair (Natural product research), 15. doi: 10.1080/14786419.2019.1710702

  • 176

    WangY.VerbruggheE.MeurisL.ChiersK.KellyM.StrubbeD.et al. (2021). Epidermal galactose spurs chytrid virulence and predicts amphibian colonization. Nat. Commun.12, 5788. doi: 10.1038/s41467-021-26127-9

  • 177

    WarneR. W.LaBumbardB.LaGrangeS.VredenburgV. T.CatenazziA. (2016). Co-infection by chytrid fungus and ranaviruses in wild and harvested frogs in the Tropical Andes. PloS One11, e0145864. doi: 10.1371/journal.pone.0145864

  • 178

    WattsA. G.SchlichtingP. E.BillermanS. M.JesmerB. R.MichelettiS.FortinM. J.et al. (2015). How spatio-temporal habitat connectivity affects amphibian genetic structure. Front. Genet.6, 275.

  • 179

    WomackM. C.SteigerwaldE.BlackburnD. C.CannatellaD. C.CatenazziA.CheJ.et al. (2022). State of the Amphibia 2020: A review of five years of amphibian research and existing resources. Ichthyology Herpetology110, 638661.

  • 180

    ZamberlettiP.ZaffaroniM.AccatinoF.CreedI. F.De MicheleC. (2018). Connectivity among wetlands matters for vulnerable amphibian populations in wetlandscapes. Ecol. Model.384, 119127.

  • 181

    ZamudioK. R.McDonaldC. A.BelasenA. M. (2020). High variability in infection mechanisms and host responses: A review of functional genomic studies of amphibian chytridiomycosis. Herpetologica76, 189200.

  • 182

    ZayaR. M.AminiZ.WhitakeA. S.IdeC. F. (2011). Exposure to atrazine affects the expression of key genes in metabolic pathways integral to energy homeostatis in Xenopus laevis tadpoles. Aquat Toxicol.104, 254262.

  • 183

    ZhuW.ChangL.ZhaoT.WangB.JiangJ. (2020). Remarkable metabolic reorganization and altered metabolic requirements in frog metamorphic climax. Front. Zool17, 30. doi: 10.1186/s12983-020-00378-6

Summary

Keywords

amphibian, conservation, risk assessment, stressor, physiology

Citation

Awkerman JA, Glinski DA, Henderson WM, Meter RV and Purucker ST (2024) Framework for multi-stressor physiological response evaluation in amphibian risk assessment and conservation. Front. Ecol. Evol. 12:1336747. doi: 10.3389/fevo.2024.1336747

Received

11 November 2023

Accepted

26 February 2024

Published

13 March 2024

Volume

12 - 2024

Edited by

Mary Ann Ottinger, University of Houston, United States

Reviewed by

Paula Henry, Patuxent Wildlife Research Center (USGS), United States

Caitlin R. Gabor, Texas State University, United States

Updates

Copyright

*Correspondence: Jill A. Awkerman,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics