It is unclear how suitable human-made wetlands are for supporting wildlife and how they impact wildlife disease risk. Natural wetlands (those that were created without human actions) can support more diverse and resilient communities that are at lower risk of disease outbreaks. We compared frog community composition and infection with the pathogenic fungus Batrachochytrium dendrobatidis (Bd) between human-made and natural wetlands in Tippecanoe County, Indiana, US. We conducted visual encounter surveys of frog communities and quantified Bd infection prevalence at four natural and five human-made wetlands. Water parameters associated with human practices (e.g., pH, salinity) and surrounding land use were also compared across sites. We found higher Bd infection prevalence at human-made sites than at natural sites, with monthly differences showing highest infection in spring and fall, and decreasing infection with increasing water temperature. However, we found no differences between human-made and natural sites regarding amphibian community composition, water quality, or surrounding land use. Further, we found frog density increased with distance to nearest roads among both human-made and natural sites. These findings might suggest that human-made wetlands can support frog communities similar to natural wetlands, but pose a greater risk of Bd infection.

The fungal pathogen Batrachochytrium dendrobatidis (Bd) causes the emerging infectious disease, chytridiomycosis, in amphibians (Longcore et al. 1999). Chytridiomycosis is a leading factor in global population declines and species extinctions of amphibians in recent decades (Scheele et al. 2019). Batrachochytrium dendrobatidis infections can vary spatially within and among species due to differences in habitat quality, temperature, and host community composition (Kriger and Hero 2007; Searle et al. 2011). Exploring amphibian communities and their diseases across varying levels of human activity is important for developing methods to control infectious disease and manage productive wetlands.

Anthropogenic practices can influence Bd infection dynamics through direct effects on the pathogen, often decreasing disease risk for hosts (Kriger and Hero 2007; Stockwell et al. 2015). For example, elevated levels of salinity (Stockwell et al. 2012, 2015), and pH at or below 4 and at or above 8 reduce growth and suggest reduced transmission (Piotrowski et al. 2004). Similarly, pesticides from agricultural fields can inhibit in vitro Bd growth and abundance, and increase survival rates of infected individuals (Reeves et al. 2017). Water temperature, influenced by runoff or canopy shading, can alter Bd growth, which grows best between 17 C and 25 C (Piotrowski et al. 2004; Becker et al. 2012; Heard et al. 2014). Human-caused habitat loss in the tropics was correlated with decreased Bd infection prevalence and intensity, possibly due to increased temperature from reduced canopy cover (Becker and Zamudio 2011). Conversely, increased human population densities and traffic can increase Bd introductions, thereby spreading or maintaining disease risk (Morgan et al. 2007; Murray et al. 2011). The anthropogenic practices that can impact Bd infection dynamics are also known to influence host population structure and susceptibility.

Environmental factors that negatively affect Bd proliferation can also reduce amphibian survival and fitness (Relyea 2004). Human-made wetlands are established for water retention, water regulation, and aesthetics (Brand and Snodgrass 2010), and are more prone to anthropogenic disturbances such as runoff pollution or fragmentation. Natural wetlands are associated with vegetation or forest buffers from anthropogenic disturbances that are experienced in human-made wetlands. Exposure to pesticides and heavy metals in human-made wetlands can have immediate and long-term effects on growth, survival, reproductive fitness, and the population stability of amphibians (Sanzo and Hecnar 2006; Harless et al. 2011). Wood frog (Lithobates sylvaticus) tadpoles exposed to increased concentrations of chloride showed decreased survivorship, weight, and activity, while also experiencing faster time to metamorphosis and increased rates of developmental abnormalities (Sanzo and Hecnar 2006). Additionally, contaminant exposure can cause immunosuppression in frogs, increasing their susceptibility to pathogen infection (Belden and Kiesecker 2005; Peterson et al. 2013; McMahon et al. 2017; Rollins-Smith 2017). Human-induced changes to habitat connectivity and patch size can indirectly affect host-pathogen interactions through altered species density and diversity (Cumming 2002; Molyneaux et al. 2002; Power and Mitchell 2004; Johnson et al. 2008). Habitat connectivity can support frog metacommunities, in part by providing environmental refugia from pathogens (Heard et al. 2015). Investigating how anthropogenic disturbances affect both wildlife and their pathogens will clarify of how humans can impact ecological communities.

Human-made wetlands, both those established unintentionally or for water retention, are typically more prone to anthropogenic influence (e.g., increased runoff pollution) than are wetlands established naturally (Borne 2014; Chumchal and Drenner 2015). Such anthropogenic impacts can negatively affect local amphibian species and influence disease dynamics. It is not known if human-made and natural wetlands present the same risk of Bd infection for amphibian hosts.

We compared Bd infection prevalence between human-made versus natural wetlands to better understand the role of anthropogenic influence on Bd infection dynamics. We recorded frog density, community composition, water parameters, and prevalence of Bd infection at nine sites. We hypothesized that prevalence of Bd infection in human-made sites would be higher than in natural sites because human-related environmental impacts would negatively affect frogs' ability to resist infection. We also hypothesized that frog community composition would differ between natural sites because anthropogenic disturbances in human-made sites would be detrimental to frog population growth and stability. Examining the habitat characteristics that impact Bd infections can shed light on the risk of disease outbreaks in various locations and the environmental factors keeping infections at endemic levels.

Study sites

We sampled nine sites in Tippecanoe County, Indiana, USA from May 2017 to October 2017 (Fig. 1 and Table 1). Only Mulvey Pond was sampled three times (May, June, and July), whereas all other sites except one were sampled five or six times (Ross Reserve was sampled seven times). All sites fall within the Central Till Plain Natural Region of the Wabash Drainage System. Site elevations ranged from 177 m to 234 m (mean=200 m). Average monthly rainfall for this area during sampling was 12.7 cm, with an average monthly temperature of 19 C (range=14–22 C), which overlaps with optimum growth temperature for Bd (17–25 C; Voyles et al. 2012; US Climate Data 2020). The months during which our sampling occurred also overlapped with the breeding phenology of the known frog species inhabiting these sites, and provide a temporal range to examine changes in frog community and infection parameters. All sampling events were conducted at night between 2100 hours and 0300 hours. Sites were selected if they were permanent wetlands and supported at least one frog species. We identified site origin (human-made or natural) by consulting with landowners and managers. Human-made wetlands were established by complete or partial development with human involvement and consisted of dugout drainage ditches, water retention ponds, and other means to regulate river flow (Fig. 2A, B; Brand and Snodgrass 2010). Wetlands were classified as natural if there was no evidence of human involvement in their creation or water flow control, and included ponds, marshes, and prairie pothole wetlands (Fig. 2C, D). Characteristics of water quality were collected during each visit along a randomly chosen transect to document microhabitat characteristics of sampled frogs, and included water and air temperature, salinity, total dissolved solids (TDS), conductivity, and pH with a Multi-Parameter PCSTestr (Oakton Instruments, Vernon Hills, Illinois, USA) at 5–8 cm below the water surface. We disinfected all equipment between sites to minimize transmission of Bd and other pathogens with 10% bleach, or by rinsing off organic matter and allowing equipment to completely dry for 24 hr between sites (Johnson et al. 2003).

Figure 1

Map of nine sites sampled for fungal pathogen infection of Batrachochytrium dendrobatidis in frog communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017. Circles (•) represent sites of human-made origin (n=5) and triangles (▲) represent sites of natural origin (n=4). Surrounding land use was categorized into one of three groups from the National Land Cover Database: Agriculture, containing cultivated-pasture classes; forest, containing forest classes; and urban, containing developed classes. Inset map of the US shows state of Indiana in black, with white rectangle representing Tippecanoe County where all sampling sites are located.

Figure 1

Map of nine sites sampled for fungal pathogen infection of Batrachochytrium dendrobatidis in frog communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017. Circles (•) represent sites of human-made origin (n=5) and triangles (▲) represent sites of natural origin (n=4). Surrounding land use was categorized into one of three groups from the National Land Cover Database: Agriculture, containing cultivated-pasture classes; forest, containing forest classes; and urban, containing developed classes. Inset map of the US shows state of Indiana in black, with white rectangle representing Tippecanoe County where all sampling sites are located.

Close modal
Table 1

Characteristics of nine sites in Tippecanoe County, Indiana, USA where frog populations were sampled for the prevalence of the fungal pathogen Batrachochytrium dendrobatidis from May 2017 to October 2017. The sites were categorized by wetland origin: human-made (HM) or natural (NAT). Percent forest, urban, and agriculture are percent of surrounding land use within a 500-m buffer around each wetland and road proximity, measured as distance to nearest paved road, were collected to explore impacts of surrounding land-use characteristics on B. dendrobatidis infections and frog community attributes.

Characteristics of nine sites in Tippecanoe County, Indiana, USA where frog populations were sampled for the prevalence of the fungal pathogen Batrachochytrium dendrobatidis from May 2017 to October 2017. The sites were categorized by wetland origin: human-made (HM) or natural (NAT). Percent forest, urban, and agriculture are percent of surrounding land use within a 500-m buffer around each wetland and road proximity, measured as distance to nearest paved road, were collected to explore impacts of surrounding land-use characteristics on B. dendrobatidis infections and frog community attributes.
Characteristics of nine sites in Tippecanoe County, Indiana, USA where frog populations were sampled for the prevalence of the fungal pathogen Batrachochytrium dendrobatidis from May 2017 to October 2017. The sites were categorized by wetland origin: human-made (HM) or natural (NAT). Percent forest, urban, and agriculture are percent of surrounding land use within a 500-m buffer around each wetland and road proximity, measured as distance to nearest paved road, were collected to explore impacts of surrounding land-use characteristics on B. dendrobatidis infections and frog community attributes.
Figure 2

Representative examples of the two site types (categorized by site origin), human-made (A, B) and natural (C, D) wetlands, that we sampled for fungal pathogen infection of Batrachochytrium dendrobatidis in frog communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017. (A) McCormick Woods (Birk Boilermaker Golf Complex), established in 2017; (B) Martell Forest Water Retention Pond, established in 2002; (C) Martell Forest Pond; (D) Purdue Wildlife Area.

Figure 2

Representative examples of the two site types (categorized by site origin), human-made (A, B) and natural (C, D) wetlands, that we sampled for fungal pathogen infection of Batrachochytrium dendrobatidis in frog communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017. (A) McCormick Woods (Birk Boilermaker Golf Complex), established in 2017; (B) Martell Forest Water Retention Pond, established in 2002; (C) Martell Forest Pond; (D) Purdue Wildlife Area.

Close modal

In addition to on-site data, we characterized surrounding land-use type using a buffer analysis within a 500 m radius of each site in R version 3.4.1 (packages: raster, rgdal; Hijmans and van Etten 2012; Bivand et al. 2014; R Core Team 2018). This allowed us to calculate the proportion of surrounding land use based on the 2011 National Land Cover Database (Homer et al. 2015). Based on the National Land Cover Database, we grouped land use into an agriculture category that contained cultivated-pasture classes, an urban category that contained developed classes, and a forest category that contained forest classes (Fig. 1). Nearest road proximity was measured using aerial photos with Google Earth Pro 7.3 (Google Inc., Mountain View, California, USA) as the shortest distance (m) between a paved road and the shoreline of the wetland.

Visual encounter survey

To measure frog communities, three to four visual encounter surveys (VES) were conducted at every site during each survey period, depending on wetland size. Each VES was 30 m in length and evenly distributed across each wetland. We captured up to 20 individuals per species with gloved hands or large landing net at each site. Gloves were changed and discarded between each captured individual during VES and swabbing, and nets were cleaned of organic debris. Upon capture, frogs were immediately placed in separate plastic bags to record individual data (species, age, sex, and body length [snout-vent length]), and to conduct skin swabbing for Bd infections.

Batrachochytrium dendrobatidis collection and quantitative (q)PCR analysis

To collect samples for Bd testing, one MW 113 rayon-tipped culture swab (Medical Wire and Equipment Co. Ltd., Corsham, England) was passed along three areas of each frog for a total of 40 swipes (10 on the ventrum, 10 on each inner thigh, one under each toe) while wearing latex-free nitrile gloves (Hyatt et al. 2007). The swab was immediately placed in a 1.5 mL micro-centrifuge tube and allowed to air dry away from continued swabbing for approximately 15 min to reduce potential inhibitory effects of water in the sample. Swabs were then stored in a –20 C freezer until DNA extraction.

Extraction and quantification of Bd from swabs was performed following standard protocols for a qPCR Taqman assay (Boyle et al. 2004; Hyatt et al. 2007) using gBlocks (Integrated DNA Technologies, Coralville, Iowa, USA) for Bd internal transcribed spacer genes as standards, which included four serial dilutions in duplicate on each plate. Each plate was run on a CFX Connect Real-Time System machine (Bio-Rad Laboratories, Hercules, California, USA). To detect inhibition within the reactions, TaqMan Exogenous Internal Positive Controls (Applied Biosystem, Foster City, California, USA) were included in each sample (Hyatt et al. 2007). When no amplification occurred and if the qPCR reaction was not inhibited, the sample was considered negative for the presence of Bd.

All field sampling was conducted under permission of the Purdue Animal Care and Use Committee (no. 1607001450) and the Indiana Department of Natural Resources Scientific Purposes license (no. 17-020).

Statistical analyses

First, we tested the influence of biotic and abiotic variables on infection status of frogs. We used a generalized linear mixed-effects model (GLMM) with binomial error distribution to identify factors that best predicted the response variable of infection status (infected or not infected) of each individual sampled for Bd. We fit a GLMM model using the glmer function in the lme4 package in R (Bates et al. 2015). All explanatory variables were tested for multicollinearity using variance inflation values (VIF), and a cutoff value of five was used to begin removing collinear variables (James et al. 2013). Sex and age had VIF values greater than five. We then removed sex and recalculated VIF values; all remaining values were below 5. Included as fixed effects were site origin, age, sampling month, species, snout-vent-length, and road proximity. We included sampling site nested within sampling date as a random effect to control for the fact that infection status of individuals from the same site was not independent, and sites sampled on the same date might have infection status levels that were more similar to one another than from other dates (Harrison et al. 2018). Starting with this global model, a set of GLMMs was created, followed by model selection (Burnham and Anderson 2002). Model assessment was ranked by Akaike information criterion (AIC), and multiple comparisons of differences within single significant factors were determined using the drop1 function (Bates et al. 2015). A threshold ΔAIC of two was used to distinguish differences among models (Burnham and Anderson 2002).

Because all habitat characteristic parameters were not measured at every pathogen sampling event (all parameters measured at 23/34 events) due to lack of or malfunctioning equipment, we used data from only those sampling events that included all water parameter measurements in a separate binomial GLMM to examine effects on Bd infection status (all other variables collected on same day habitat parameters were still included). We constructed a GLMM with binomial error distribution similar to above. All explanatory variables were scaled (mean=0, SD=1) to bring the size and variability of the different data to a common scale. We tested for multicollinearity of explanatory variables using VIF, and a cutoff value of five was used to remove collinear variables (James et al. 2013). Conductivity, salinity, and TDS had VIF values greater than five. We then removed conductivity and salinity and all remaining values were below five. Included as fixed effects were air temperature, water temperature, and pH. We included sampling site nested within date collected as a random effect. We constructed a set of GLMMs starting with a global model that included all variables, followed by model selection similar to above (Burnham and Anderson 2002).

Next, we examined how frog density was influenced by multiple environmental variables we collected. We used stepwise model selection of constructed linear regression models to determine the most impactful explanatory variables for frog density (mean frogs/m2) and interpreted the best fit model determined by AIC with correction for small sample size. A threshold ΔAIC of two was used to distinguish differences among models (Burnham and Anderson 2002). Explanatory variables analyzed were TDS, salinity, conductivity, pH, surrounding land use, site origin, and distance to nearest road. All explanatory variables were scaled (mean=0, SD=1) to bring the size and variability of the different data to a common scale. We found pH, conductivity, salinity, and TDS were collinear, so we removed all but pH, which had the lowest VIF value. We also found all land-use proportion estimates were collinear and removed all except proportion urban land use. We recalculated VIF values and all remaining values were below five.

To determine if there were differences in water quality in wetlands that were natural versus human-made, we conducted Welch's two sample t-tests for average water parameters (pH, salinity, conductivity, TDS, and water temperature) and amphibian density. Data were log-transformed to meet assumptions of parametric statistics.

Differences in species composition between human-made and natural sites were examined via nonmetric multidimensional scaling (NMDS) with the Bray-Curtis dissimilarity measurement (Venables and Ripley 2002) using the R package vegan (Oksanen et al. 2011). After performing the NMDS, a permutational multivariate analysis of variance to test for statistical differences in species composition between site types (human-made versus natural). To directly compare alpha diversity, we tested Shannon-Wiener diversity between site types (human-made versus natural) via one-way analysis of variance. All statistical analyses were performed in R version 3.5.1 (R Core Team 2018).

Mean Bd infection prevalence was 33% (118/361) across all samples, with five of seven species having at least one infected individual (Table 2 and Supplementary Material Table S1). The best supported GLMM model predicting infection included site origin, species, and month (Table 3). The mean prevalence of Bd infection was higher at human-made sites (38% [87/230]) than at natural sites (24% [31/131]; Fig. 3). Species and month were also significantly related to Bd infection prevalence. Multiple comparison analyses of species' probability to test positive showed American bullfrogs (Lithobates catesbeianus) was most likely to be infected, whereas gray treefrogs (Hyla versicolor), pickerel frogs (Lithobates palustris), and spring peepers (Pseudacris crucifer) were the least likely to be infected (Table 2). Infection prevalence among months showed a bimodal distribution with lowest percent infection in July (Fig. 4 and Table 2). The most robust model predicting infection based on habitat characteristics contained only water temperature (Table 4), in which there was a decreasing probability of infection with increasing temperature.

Table 2

The prevalences and frequencies of the fungal pathogen Batrachochytrium dendrobatidis in samples of frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Infection was compared across species, age, month, and sex. Infection prevalence was calculated as a percent (n positive/n tested), with 95% lower and upper binomial confidence intervals (CIs).

The prevalences and frequencies of the fungal pathogen Batrachochytrium dendrobatidis in samples of frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Infection was compared across species, age, month, and sex. Infection prevalence was calculated as a percent (n positive/n tested), with 95% lower and upper binomial confidence intervals (CIs).
The prevalences and frequencies of the fungal pathogen Batrachochytrium dendrobatidis in samples of frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Infection was compared across species, age, month, and sex. Infection prevalence was calculated as a percent (n positive/n tested), with 95% lower and upper binomial confidence intervals (CIs).
Table 3

Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Best model included site origin (human-made versus natural), species, and month as the predictor variables based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P values derived from drop1 function to test single fixed effects.

Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Best model included site origin (human-made versus natural), species, and month as the predictor variables based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P values derived from drop1 function to test single fixed effects.
Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs from nine human-made (n=5) and natural (n=4) sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Best model included site origin (human-made versus natural), species, and month as the predictor variables based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P values derived from drop1 function to test single fixed effects.
Figure 3

Percent of frogs infected with the fungal pathogen Batrachochytrium dendrobatidis in communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017 between sites of human-made or natural origin. Percent infection was pooled across all species, sites, and months. Asterisks (*) denote mean percent infection. Human-made sites had significantly higher infection prevalence than sites of natural origin.

Figure 3

Percent of frogs infected with the fungal pathogen Batrachochytrium dendrobatidis in communities in Tippecanoe County, Indiana, USA from May 2017 to October 2017 between sites of human-made or natural origin. Percent infection was pooled across all species, sites, and months. Asterisks (*) denote mean percent infection. Human-made sites had significantly higher infection prevalence than sites of natural origin.

Close modal
Figure 4

Percent of frogs infected with Batrachochytrium dendrobatidis in Tippecanoe County, Indiana, USA from May 2017 to October 2017 across sampling months. Percent infection was pooled across all sites and site visits within a month. Asterisks (*) denote mean percent infection. Significant differences are indicated by nonmatching letters above each boxplot.

Figure 4

Percent of frogs infected with Batrachochytrium dendrobatidis in Tippecanoe County, Indiana, USA from May 2017 to October 2017 across sampling months. Percent infection was pooled across all sites and site visits within a month. Asterisks (*) denote mean percent infection. Significant differences are indicated by nonmatching letters above each boxplot.

Close modal
Table 4

Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs. Data are from nine sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. This model only included data from sampling events that included all water parameter measurements (23/34 sampling events). Best model included water temperature as the only predictor variable based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P value derived from drop1 function to test single fixed effects.

Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs. Data are from nine sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. This model only included data from sampling events that included all water parameter measurements (23/34 sampling events). Best model included water temperature as the only predictor variable based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P value derived from drop1 function to test single fixed effects.
Final predictor variables from best supported generalized linear mixed model (GLMM) with binomial error distribution predicting infection of the fungal pathogen Batrachochytrium dendrobatidis in frogs. Data are from nine sites in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. This model only included data from sampling events that included all water parameter measurements (23/34 sampling events). Best model included water temperature as the only predictor variable based on Akaike information criterion values. Top three GLMMs and interaction model of best GLMM listed. P value derived from drop1 function to test single fixed effects.

The best-supported model predicting frog density included only road proximity (Table 5), with denser frog communities at further distances to the nearest road (F(1,7)=11.66, P=0.011, r2=0.571; Fig. 5). Mean amphibian density did not differ between site origins (t(3,5)=–0.57, P=0.601), nor was site origin included in the best-supported model predicting frog density. Welch's t-tests revealed no significant differences between site origin for salinity (t(6.5)=0.95, P=0.37), conductivity (t(6.5)=0.97, P=0.37), pH (t(6.4)=1.01, P=0.35), or TDS (t(10.2)=1.26, P=0.24; Fig. 6), or water temperature (t(4.2)=–0.22, P=0.83).

Table 5

Linear model selection for mean frog density (frogs/m2) and explanatory site variables among all nine sites (human-made [n=5] and natural [n=4]) sampled in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Final global model (after removing collinear variables) included road proximity, pH, site origin, and proportion urban land use (proportion urban only). Best model included road proximity as the only predictor variable based on corrected Akaike information criterion values.

Linear model selection for mean frog density (frogs/m2) and explanatory site variables among all nine sites (human-made [n=5] and natural [n=4]) sampled in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Final global model (after removing collinear variables) included road proximity, pH, site origin, and proportion urban land use (proportion urban only). Best model included road proximity as the only predictor variable based on corrected Akaike information criterion values.
Linear model selection for mean frog density (frogs/m2) and explanatory site variables among all nine sites (human-made [n=5] and natural [n=4]) sampled in Tippecanoe County, Indiana, USA sampled from May 2017 to October 2017. Final global model (after removing collinear variables) included road proximity, pH, site origin, and proportion urban land use (proportion urban only). Best model included road proximity as the only predictor variable based on corrected Akaike information criterion values.
Figure 5

Density of frog communities sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 showing a significant positive relationship with distance to nearest paved road among all sites (P=0.011, r2=0.571).

Figure 5

Density of frog communities sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 showing a significant positive relationship with distance to nearest paved road among all sites (P=0.011, r2=0.571).

Close modal
Figure 6

Estimates of water parameters sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 compared between human-made and natural sites pooled across all sampling events, with mean values indicated by asterisks (*). No significant differences between human-made and natural sites were determined for all water parameters based on Welch's t-tests.

Figure 6

Estimates of water parameters sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 compared between human-made and natural sites pooled across all sampling events, with mean values indicated by asterisks (*). No significant differences between human-made and natural sites were determined for all water parameters based on Welch's t-tests.

Close modal

Based on the NMDS plot of amphibian assemblages (Fig. 7), community composition did not differ between habitat types. The overlap in datapoints and lack of separation suggest the compositions are quite similar between site types (F(1,8)=0.223, R=0.031, P=0.943). Further, we found no difference in Shannon-Wiener diversity between human-made and natural sites (F(1,7)=0.004, P=0.951).

Figure 7

Nonmetric multidimensional scaling plot of nonsignificant variation in frog community assemblage sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 between human-made (n=5) and natural sites (n=4). Ellipses represent 0.75 confidence intervals around the centroid of human-made (gray) and natural (white) points.

Figure 7

Nonmetric multidimensional scaling plot of nonsignificant variation in frog community assemblage sampled in Tippecanoe County, Indiana, USA from May 2017 to October 2017 between human-made (n=5) and natural sites (n=4). Ellipses represent 0.75 confidence intervals around the centroid of human-made (gray) and natural (white) points.

Close modal

We found that Bd infection prevalence in frogs was higher in human-made sites than in natural sites (Fig. 3), consistent with our hypothesis. This suggests that frogs occupying habitats with greater potential of human disturbance are more likely to acquire or maintain infection. Frog density and community composition did not differ between wetlands with different site origins, indicating that human-made wetlands can support frog populations that are similar to natural wetlands. Therefore, human-made wetlands can play a role in sustaining frog communities, but could pose a greater risk for Bd infection than natural wetlands.

The mechanism by which human-made wetlands have increased Bd infection compared to natural wetlands is unclear, but factors associated with human activity can impact wildlife disease dynamics (Bradley and Altizer 2006). In our study, human-made sites had, on average, higher levels of salinity, conductivity, TDS, and pH, although differences were nonsignificant. In contrast, Bd growth and host infection load has been shown to decrease at salinity concentrations over 1 part per thousand (1 g/L) in field and laboratory settings (Stockwell et al. 2012, 2015). Because we did not record salinity above 0.3 g/L (mean=0.15 g/L, SD=0.09), the salinity concentrations (and collinear with conductivity) we observed at our sites likely had little to no direct impact on Bd even if they were having sublethal effects on the amphibian hosts. It is also possible that other immunosuppressive abiotic factors that we did not measure, such as nitrate concentration, dissolved oxygen, pesticide, or heavy metal concentrations, are higher in human-made wetlands and are increasing host susceptibility to Bd infections (Christin et al. 2003; Parris and Baud 2004; Calderon et al. 2019). Examining immunocompromising stressors that are common in human-made wetlands could reduce negative effects associated with pathogen infection.

Differences in infection prevalence among months revealed patterns that are useful for monitoring Bd outbreaks and amphibian population dynamics (Fig. 4). The prevalence of Bd infection typically declines in many regions of North America in the summer due to high temperatures that limit Bd growth or cause Bd mortality, similar to our findings (Piotrowski et al. 2001; Longcore et al. 2007). This pattern was further supported by the negative relationship that we found between water temperature and infection probability, as temperatures increased above the optimal Bd growth range of 17–25 C (Piotrowski et al. 2004). Higher infection prevalence in spring could be due to increased frog densities when breeding, and frogs spending more time in Bd-inoculated water, providing more opportunities for transmission (Carey et al. 2006). Alternatively, stressful conditions caused by springtime reproduction, development, or contaminant runoff with winter snowmelt, can limit frog resistance or tolerance to Bd (Rollins-Smith 2001; Sanzo and Hecnar 2006; Karraker et al. 2008; Denver 2009). High infection prevalence in the fall suggests that environmental conditions become more suitable for Bd, likely driven by lower temperatures compared to summer months (Berger et al. 2004). Alternatively, amphibians experience increased susceptibility due to a depressed immune system following temperature decreases or cold exposure (Maniero and Carey 1997; Greenspan et al. 2017), as well as fall acclimation, potentially creating a period of increased susceptibility to pathogens (Raffel et al. 2006). Studying seasonal differences in infection prevalence is important to identify the time of year when species are most at risk of disease-related effects, and when to quantify infection in amphibian communities.

Differences in biotic conditions or human activities that we did not measure could explain the higher infection prevalence that we found in human-made sites. For example, zooplankton, which can consume Bd zoo-spores, potentially reducing infection risk in amphibians (Buck et al. 2011; Searle et al. 2013; Schmeller et al. 2014), might be negatively affected by undetected aspects of human-made sites. Further, wetland structure and vegetative diversity can influence frog diversity (Wassens et al. 2010), which can itself impact disease dynamics (Searle et al. 2011). Therefore, differences in diversity or structure of aquatic vegetation might have impacted disease dynamics, and should be considered in future studies. Additionally, human-made sites might also have higher rates of anthropogenic activity such that increased traffic at sites with easy public access results in more human-caused introductions of Bd (Morgan et al. 2007). More human activity and impacts can also lead to increased stress and corticosterone levels that depress immune functions (Dhabhar 2002; Padgett and Glaser 2003; Rollins-Smith et al. 2011). Erecting buffers, creating habitat with refugia harmful to pathogens, increasing habitat connectivity, or implementing shoe and equipment cleaning for the public and researchers can help to reduce the human-caused spread of pathogens and stress to wildlife (Phillot et al. 2010; Heard et al. 2018).

Infection prevalence was highest in American bullfrogs and green frogs (Table 2), two species that are suspected reservoirs for Bd with typically high infection loads (Daszak et al. 2004; Longcore et al. 2007; Klemish et al. 2012). Their year-round occupancy in wetlands suggests they might harbor infection through summer and winter, playing an important role in the increased infection prevalence we observed in spring and fall. Percent infection in gray treefrogs was high (19% [10/52]) relative to other studies, some of which find 0% infection (Longcore et al. 2007). Terrestrial species such as gray treefrogs lose infection outside of the breeding season (Longcore et al. 2007), suggesting that spring breeding events increases their risk of Bd infection. Measuring Bd infection in gray treefrogs and other predominately terrestrial species outside of the breeding season can better illustrate patterns of infection in these species.

We found that both pH and road proximity significantly predicted frog density. Although water quality parameters did not differ between site types (Fig. 6), pH (which was collinear with TDS, salinity, and conductivity) was positively correlated with frog density in our final model (Table 5). Low water pH is associated with lower larval amphibian survival, delayed metamorphosis, and increased adult mortality, and thus adults prefer more neutral environments (Sadinski and Dunson 1992; Horne and Dunson 1995; Vatnick et al. 1999). Future studies should explore how these and other habitat quality characteristics (e.g., pesticide types and concentrations), impact disease dynamics and amphibian communities of larvae and postmetamorphic individuals. Our results show distances to the nearest paved road are positively correlated with frog density across all sites (Fig. 5). This could be an important driver of frog community composition. Closer proximity to anthropogenic features, such as roads or residential landscaping, can result in amphibians being less common due to poor migration success from habitat fragmentation between wetlands and upland forests, and impacts of surface water pollution (Lehtinen et al. 1999; Knutson et al. 2000; Gagne and Fahrig 2007; Dodson 2008). Amphibians are at particularly high risk of road mortality in the region we sampled (Tippecanoe County) where a 2–yr-long study found 93% (9,809/10,515) vertebrate road mortalities were amphibians (Glista et al. 2008). Thus, human disturbances pose great risk to frog communities, regardless of whether the sites are natural or human-made. The similarity of these risks across sites could be why we did not observe differences in density or community composition between site types.

Our findings suggest that human-made and natural wetlands can act as comparably suitable habitats for frogs, but the conditions in human-made wetlands might be making frogs more susceptible to pathogen infection. The combined effect of anthropogenic disturbances and frog diseases can ultimately influence frog population numbers, and subsequently cause breakdowns of local food webs and the community structure of wetland habitats (Petchey et al. 1999). The creation of wetlands by humans, whether intended or not, can act as viable habitat for species being displaced by human encroachment. Therefore, ditches, water retention ponds, and other human-made wetlands should receive greater attention by conservation managers to increase and maintain the quality of these wetlands for the preservation of native species while also considering the effects of these activities on infectious disease.

We thank the Purdue University Department of Biological Sciences for funding and resources. We thank the Department of Forestry and Natural Resources at Purdue University and the Ross Biological Reserve for access to field sites. M.C.B. was funded through the Summer Stays program at Purdue University. We thank A. Catenazzi, T. DeBlieux, J. Hoverman, A. Shepack, and H. Zumbado-Ulate for training and processing of swab samples. We thank R. Lim, L. Rosario, and J. Toth for help with sample collection and processing.

SUPPLEMENTARY MATERIAL

Supplementary material for this article is online at http://dx.doi.org/10.7589/2019-09-220.

Bates
D
,
Mächler
M
,
Bolker
BM
,
Walker
SC
.
2015
.
Fitting linear mixed-effects models using lme4.
J Stat Softw
67
:
1
48
.
Becker
CG
,
Rodriguez
D
,
Longo
A
,
Talaba
A
,
Zamudio
K
.
2012
.
Disease risk in temperate amphibian populations is higher at closed-canopy sites.
PLoS One
7
:
e48205
.
Becker
CG
,
Zamudio
KR
.
2011
.
Tropical amphibian populations experience higher disease risk in natural habitats.
Proc Natl Acad Sci U S A
108
:
9893
9898
.
Belden
LK
,
Kiesecker
JM
.
2005
.
Glucocorticosteroid hormone treatment of larval treefrogs increases infection by Alaria sp. trematode cercariae.
J Parasitol
91
:
686
688
.
Berger
L
,
Speare
R
,
Hines
HB
,
Marantelli
G
,
Hyatt
AD
,
McDonald
KR
,
Skerratt
LF
,
Olsen
V
,
Clarke
JM
,
Gillespie
G
, et al.
2004
.
Effect of season and temperature on mortality in amphibians due to chytridiomycosis.
Aust Vet J
82
:
31
36
.
Bivand
R
,
Keitt
T
,
Rowlingson
B
.
2014
.
rgdal: Bindings for the “geospatial” data abstraction library. R package version 1.
.
Borne
KE.
2014
.
Floating treatment wetland influences on the fate and removal performance of phosphorus in stormwater retention ponds.
Ecol Eng
69
:
76
82
.
Boyle
DG
,
Boyle
DB
,
Olsen
V
,
Morgan
JAT
,
Hyatt
AD
.
2004
.
Rapid quantitative detection of chytridiomycosis (Batrachochytrium dendrobatidis) in amphibian samples using real-time Taqman PCR assay.
Dis Aquat Org
60
:
141
148
.
Bradley
CA
,
Altizer
S
.
2006
.
Urbanization and the ecology of wildlife diseases.
Trends Ecol Evol
22
:
95
102
.
Brand
AB
,
Snodgrass
JW
.
2010
.
Value of artificial habitats for amphibian reproduction in altered landscapes.
Conserv Biol
24
:
295
301
.
Buck
JC
,
Truong
L
,
Blaustein
AR
.
2011
.
Predation by zooplankton on Batrachochytrium dendrobatidis: Biological control of the deadly amphibian chytrid fungus?
Biodivers Conserv
20
:
3549
3553
.
Burnham
KP
,
Anderson
DR
.
2002
.
Model selection and multimodel inference: A practical information-theoretic approach.
Springer Science and Business Media
,
New York, New York
,
490
pp.
Calderon
MR
,
Almeida
CA
,
González
P
,
Jofré
MB
.
2019
.
Influence of water quality and habitat conditions on amphibian community metrics in rivers affected by urban activity.
Urban Ecosyst
22
:
743
755
.
Carey
C
,
Bruzgul
JE
,
Livo
LJ
,
Walling
ML
,
Kuehl
KA
,
Dixon
BF
,
Pessier
AP
,
Alford
RA
,
Rogers
KB
.
2006
.
Experimental exposures of boreal toads (Bufo boreas) to a pathogenic chytrid fungus (Batrachochytrium dendrobatidis).
EcoHealth
3
:
5
21
.
Christin
M
,
Gendron
AD
,
Brousseau
P
,
Ménard
L
,
Marcogliese
DJ
,
Cyr
D
,
Ruby
S
,
Fournier
M
.
2003
.
Effects of agricultural pesticides on the immune system of Rana pipiens and on its resistance to parasitic infection.
Environ Toxicol Chem
22
:
1127
1133
.
Chumchal
MM
,
Drenner
RW
.
2015
.
An environmental problem hidden in plain sight? Small human-made ponds, emergent insects, and mercury contamination of biota in the great plains.
Environ Toxicol Chem
34
:
1197
1205
.
Cumming
GS.
2002
.
Comparing climate and vegetation as limiting factors for species ranges of African ticks.
Ecology
83
:
255
268
.
Daszak
P
,
Strieby
A
,
Cunningham
AA
,
Longcore
JE
,
Brown
CC
,
Porter
D
.
2004
.
Experimental evidence that the bullfrog (Rana catesbeiana) is a potential carrier of chytridiomycosis, an emerging fungal disease of amphibians.
Herpetol J
14
:
201
207
.
Denver
RJ.
2009
.
Stress hormones mediate environment-genotype interactions during amphibian development.
Gen Comp Endocrinol
164
:
20
31
.
Dhabhar
FS.
2002
.
A hassle a day may keep the doctor away: Stress and the augmentation of immune function.
Integr Comp Biol
42
:
556
564
.
Dodson
SI.
2008
.
Biodiversity in southern Wisconsin storm-water retention ponds: Correlations with watershed cover and productivity.
Lake Reservoir Manag
24
:
370
380
.
Gagné
SA
,
Fahrig
L
.
2007
.
Effect of landscape context on anuran communities in breeding ponds in the National Capital Region, Canada.
Landsc Ecol
22
:
205
215
.
Glista
DJ
,
Devault
TL
,
DeWoody
JA
.
2008
.
Vertebrate road mortality predominantly impacts amphibians.
Herpetol Conserv Biol
3
:
77
87
.
Greenspan
SE
,
Bower
DS
,
Webb
RJ
,
Berger
L
,
Rudd
D
,
Schwarzkopf
L
,
Alford
RA
.
2017
.
White blood cell profiles in amphibians help to explain disease susceptibility following temperature shifts.
Dev Comp Immunol
77
:
280
286
.
Harless
ML
,
Huckins
CJ
,
Grant
JB
,
Pypker
TG
.
2011
.
Effects of six chemical deicers on larval wood frogs (Rana sylvatica).
Environ Toxicol Chem
30
:
1637
1641
.
Harrison
XA
,
Donaldson
L
,
Correa-Cano
ME
,
Evans
J
,
Fisher
DN
,
Goodwin
CED
,
Robinson
BS
,
Hodgson
DJ
,
Inger
R
.
2018
.
A brief introduction to mixed effects modelling and multi-model inference in ecology.
PeerJ
6
:
e4794
.
Heard
GW
,
Scroggie
MP
,
Clemann
N
,
Ramsey
DSL
.
2014
.
Wetland characteristics influence disease risk for a threatened amphibian.
Ecol Appl
24
:
650
662
.
Heard
GW
,
Scroggie
MP
,
Ramsey
DS
,
Clemann
N
,
Hodgson
JA
,
Thomas
CD
.
2018
.
Can habitat management mitigate disease impacts on threatened amphibians?
Conserv Lett
11
:
e12375
.
Heard
GW
,
Thomas
CD
,
Hodgson
JA
,
Scroggie
MP
,
Ramsey
DSL
,
Clemann
N
.
2015
.
Refugia and connectivity sustain amphibian metapopulations afflicted by disease.
Ecol Lett
18
:
853
863
.
Hijmans
RJ
,
van Etten
J
.
2012
.
raster: Geographic analysis and modeling with raster data. R package version 2.0-05.
.
Homer
C
,
Dewitz
J
,
Yang
L
,
Jin
S
,
Danielson
P
,
Xian
G
,
Coulston
J
,
Herold
N
,
Wickham
J
,
Megown
K
.
2015
.
Completion of the 2011 National Land Cover Database for the conterminous United States—Representing a decade of land cover change information.
Photogramm Eng Remote Sens
81
:
345
354
.
Horne
MT
,
Dunson
WA
.
1995
.
Effects of low pH, metals, and water hardness on larval amphibians.
Arch Environ Contam Toxicol
29
:
500
505
.
Hyatt
AD
,
Boyle
DG
,
Olsen
V
,
Boyle
DB
,
Berger
L
,
Obendorf
D
,
Dalton
A
,
Kriger
K
,
Hero
M
,
Hines
H
, et al.
2007
.
Diagnostic assays and sampling protocols for the detection of Batrachochytrium dendrobatidis.
Dis Aquat Org
73
:
175
192
.
James
G
,
Witten
D
,
Hastie
T
,
Tibshirani
R
.
2013
.
An introduction to statistical learning with applications in R.
Springer Publishing Company
,
New York, New York
,
431
pp.
Johnson
ML
,
Berger
L
,
Philips
L
,
Speare
R
.
2003
.
Fungicidal effects of chemical disinfectants, UV light, desiccation and heat on the amphibian chytrid Batrachochytrium dendrobatidis.
Dis Aquat Org
57
:
255
260
.
Johnson
PT
,
Carpenter
SR
,
Ostfeld
R
,
Keesing
F
,
Eviner
V
.
2008
.
Influence of eutrophication on disease in aquatic ecosystems: Patterns, processes, and predictions.
In
:
Infectious disease ecology: Effects of ecosystems on disease and of disease on ecosystems
,
Ostfield
RS
,
Keesing
F
,
Eviner
VT
, editors.
Princeton University Press
,
Princeton, New Jersey
, pp.
71
79
.
Karraker
NE
,
Gibbs
JP
,
Vonesh
JR
.
2008
.
Impacts of road deicing salt on the demography of vernal pool-breeding amphibians.
Ecol Appl
18
:
724
734
.
Klemish
JL
,
Johnson
BL
,
Siddons
SR
,
Wild
ER
.
2012
.
Occurrence of Batrachochytrium dendrobatidis among populations of Lithobates clamitans and L. pipiens in Wisconsin, USA.
Herpetol Rev
43
:
282
288
.
Knutson
M
,
Sauer
JR
,
Olsen
DA
,
Mossman
MJ
,
Hemesath
LM
,
Lannoo
MJ
.
2000
.
Landscape associations of frog and toad species in Iowa and Wisconsin, USA.
J Iowa Acad Sci
107
:
134
145
.
Kriger
KM
,
Hero
JM
.
2007
.
The chytrid fungus Batrachochytrium dendrobatidis is non-randomly distributed across amphibian breeding habitats.
Divers Distrib
13
:
781
788
.
Lehtinen
RM
,
Galatowitsch
SM
,
Tester
JR
.
1999
.
Consequences of habitat loss and fragmentation for wetland amphibian assemblages.
Wetlands
19
:
1
12
.
Longcore
JR
,
Longcore
JE
,
Pessier
AP
,
Halteman
WA
.
2007
.
Chytridiomycosis widespread in anurans of northeastern United States.
J Wildl Manag
71
:
435
444
.
Longcore
JE
,
Pessier
AP
,
Nichols
DK
.
1999
.
Batrachochytrium dendrobatidis gen. et sp. nov., a chytrid pathogenic to amphibians.
Mycologia
91
:
219
227
.
Maniero
GD
,
Carey
C
.
1997
.
Changes in selected aspects of immune function in the leopard frog, Rana pipiens, associated with exposure to cold.
J Comp Physiol B
167
:
256
263
.
McMahon
TA
,
Boughton
RK
,
Martin
LB
,
Rohr
JR
.
2017
.
Exposure to the herbicide atrazine nonlinearly affects tadpole corticosterone levels.
J Herpetol
51
:
270
273
.
Molyneaux
DH.
2002
.
Vector-borne infections and health related to landscape changes.
In
:
Conservation medicine: Ecological health in practice
,
Aguirre
AA
,
Ostfeld
RS
,
Tabor
GM
,
House
C
,
Pearl
MC
, editors.
Oxford Univerity Press
,
Oxford, UK
, pp
194
206
.
Morgan
JAT
,
Vredenburg
VT
,
Rachowicz
LJ
,
Knapp
RA
,
Stice
MJ
,
Tunstall
T
,
Bingham
RE
,
Parker
JM
,
Longcore
JE
,
Moritz
C
, et al.
2007
.
Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis.
Proc Natl Acad Sci U S A
104
:
13845
13850
.
Murray
KA
,
Retallick
RWR
,
Puschendorf
R
,
Skerratt
LF
,
Rosauer
D
,
McCallum
HI
,
Berger
L
,
Speare
R
,
VanDerWal
J
.
2011
.
Assessing spatial patterns of disease risk to biodiversity: Implications for the management of the amphibian pathogen, Batrachochytrium dendrobatidis.
J Appl Ecol
48
:
163
173
.
Oksanen
J
,
Blanchet
FG
,
Kindt
R
,
Legendre
P
,
Minchin
PR
,
O'Hara
RB
,
Wagner
H
.
2011
.
Vegan: Community ecology package. R package version 1.7-4.
.
Padgett
DA
,
Glaser
R
.
2003
.
How stress influences the immune response.
Trends Immunol
24
:
444
448
.
Parris
MJ
,
Baud
DR
.
2004
.
Interactive effects of a heavy metal and chytridiomycosis on gray treefrog larvae (Hyla chrysoscelis).
Copeia
2004
:
344
350
.
Petchey
O
,
McPhearson
PT
,
Casey
TM
,
Morin
PJ
.
1999
.
Environmental warming alters food-web structure and ecosystem function.
Nature
402
:
69
72
.
Peterson
JD
,
Steffen
JE
,
Reinert
LK
,
Cobine
PA
,
Appel
A
,
Rollins-Smith
L
,
Mendonça
MT
.
2013
.
Host stress response is important for the pathogenesis of the deadly amphibian disease, chytridiomycosis, in Litoria caerulea.
PLoS One
8
:
e62146
.
Phillot
AD
,
Speare
R
,
Hines
HB
,
Skerratt
LF
,
Meyer
E
,
McDonald
KR
,
Cashins
SD
,
Mendez
D
,
Berger
L
.
2010
.
Minimising exposure of amphibians to pathogens during field studies.
Dis Aquat Org
92
:
175
185
.
Piotrowski
JS
,
Annis
SL
,
Longcore
JE
.
2001
.
Physiology, zoospore behavior, and enzyme production of Batrachochytrium dendrobatidis, a chytrid pathogenic to amphibians.
Phytopathology
91
:
S121
.
Piotrowski
JS
,
Annis
SL
,
Longcore
JE
.
2004
.
Physiology of Batrachochytrium dendrobatidis, a chytrid pathogen of amphibians.
Mycologia
96
:
9
15
.
Power
AG
,
Mitchell
CE
.
2004
.
Pathogen spillover in disease epidemics.
Am Nat
164
:
S79
S89
.
R Core Team
.
2018
.
R: A language and environment for statistical computing.
R Foundation for Statistical Computing
,
Vienna, Austria
.
http://www.R-project.org. Accessed October 2018
.
Raffel
TR
,
Rohr
JR
,
Kiesecker
JM
,
Hudson
PJ
.
2006
.
Negative effects of changing temperature on amphibian immunity under field conditions.
Funct Ecol
20
:
819
828
.
Reeves
RA
,
Pierce
CL
,
Vandever
MW
,
Muths
E
,
Smalling
KL
.
2017
.
Amphibians, pesticides, and the amphibian chytrid fungus in restored wetlands in agricultural landscapes.
Herpetol Conserv Biol
12
:
68
77
.
Relyea
RA.
2004
.
Growth and survival of five amphibian species exposed to combinations of pesticides.
Environ Toxicol Chem
23
:
1737
1742
.
Rollins-Smith
L.
2001
.
Neuroendocrine-immune system interactions in amphibians.
Immunol Res
23
:
273
280
.
Rollins-Smith
LA.
2017
.
Amphibian immunity—stress, disease, and climate change.
Dev Comp Immunol
66
:
111
119
.
Rollins-Smith
LA
,
Ramsey
JP
,
Pask
JD
,
Reinert
LK
,
Woodhams
DC
.
2011
.
Amphibian immune defenses against chytridiomycosis: Impacts of changing environments.
Integr Comp Biol
51
:
552
562
.
Sadinski
WJ
,
Dunson
WA
.
1992
.
A multilevel study of effects of low pH on amphibians of temporary ponds.
J Herpetol
26
:
413
422
.
Sanzo
D
,
Hecnar
SJ
.
2006
.
Effects of road de-icing salt (NaCl) on larval wood frogs (Rana sylvatica).
Environ Pollut
140
:
247
256
.
Scheele
BC
,
Pasmans
F
,
Skerratt
LF
,
Berger
L
,
Martel
A
,
Beukema
W
,
Acevedo
AA
,
Burrowes
PA
,
Carvalho
T
,
Catenazzi
A
, et al.
2019
.
Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity.
Science
363
:
1459
1463
.
Schmeller
DS
,
Blooi
M
,
Martel
A
,
Garner
TWJ
,
Fisher
MC
,
Azemar
F
,
Clare
FC
,
Leclerc
C
,
Jäger
L
,
Guevara-Nieto
M
, et al.
2014
.
Microscopic aquatic predators strongly affect infection dynamics of a globally emerged pathogen.
Curr Biol
24
:
176
180
.
Searle
CL
,
Biga
LM
,
Spatafora
JW
,
Blaustein
AR
.
2011
.
A dilution effect in the emerging amphibian pathogen Batrachochytrium dendrobatidis.
Proc Natl Acad Sci U S A
108
:
16322
16326
.
Searle
CL
,
Mendelson
JR
III
,
Green
LE
,
Duffy
MA
.
2013
.
Daphnia predation on the amphibian chytrid fungus and its impacts on disease risk in tadpoles.
Ecol Evol
3
:
4129
4138
.
Stockwell
MP
,
Clulow
J
,
Mahony
MJ
.
2012
.
Sodium chloride inhibits the growth and infective capacity of the amphibian chytrid fungus and increases host survival rates.
PLoS One
7
:
e36942
.
Stockwell
MP
,
Clulow
J
,
Mahony
MJ
.
2015
.
Evidence of a salt refuge: Chytrid infection loads are suppressed in hosts exposed to salt.
Oecologia
177
:
901
910
.
US Climate Data
.
2020
.
Weather history West Lafayette 2017.
www.usclimatedata.com. Accessed January 2020
.
Vatnick
I
,
Brodkin
MA
,
Simon
MP
,
Grant
BW
,
Conte
CR
,
Gleave
M
,
Myers
R
,
Sadoff
MM
.
1999
.
The effects of exposure to mild acidic conditions on adult frogs (Rana pipiens and Rana clamitans): Mortality rates and pH preferences.
J Herpetol
33
:
370
374
.
Venables
WN
,
Ripley
BD
.
2002
.
Modern applied statistics with S.
4th Ed.
Springer
,
New York, New York
,
496
pp.
Voyles
J
,
Johnson
LR
,
Briggs
CJ
,
Cashins
SD
,
Alford
RA
,
Berger
L
,
Skerratt
LF
,
Speare
R
,
Rosenblum
EB
.
2012
.
Temperature alters reproductive life history patterns in Batrachochytrium dendrobatidis, a lethal pathogen associated with the global loss of amphibians.
Ecol Evol
2
:
2241
2249
.
Wassens
S
,
Hall
A
,
Osborne
W
,
Watts
R
.
2010
.
Habitat characteristics predict occupancy patterns of the endangered amphibian Litoria raniformis in flow-regulated flood plain wetlands.
Aust Ecol
35
:
944
955
.

Supplementary data