Once a reservoir, always a reservoir? Seasonality affects the pathogen maintenance potential of amphibian hosts.

Host species that can independently maintain a pathogen in a host community and contribute to infection in other species are important targets for disease management. However, the potential of host species to maintain a pathogen is not fixed over time, and an important challenge is understanding how within- and across-season variability in host maintenance potential affects pathogen persistence over longer time scales relevant for disease management (e.g., years). Here, we sought to understand the causes and consequences of seasonal infection dynamics in leopard frogs (Rana sphenocephala and Rana pipiens) infected with the fungal pathogen Batrachochytrium dendrobatidis (Bd). We addressed three questions broadly applicable to seasonal host-parasite systems. First, to what degree are observed seasonal patterns in infection driven by temperature-dependent infection processes compared to seasonal host demographic processes? Second, how does seasonal variation in maintenance potential affect long-term pathogen persistence in multi-host communities? Third, does high deterministic maintenance potential relate to the long-term stochastic persistence of pathogens in host populations with seasonal infection dynamics? To answer these questions, we used field data collected over 3 years on >1400 amphibians across four geographic locations, laboratory and mesocosm experiments, and a novel mathematical model. We found that the mechanisms that drive seasonal prevalence were different from those driving seasonal infection intensity. Seasonal variation in Bd prevalence was driven primarily by changes in host contact rates associated with breeding migrations to and from aquatic habitat. In contrast, seasonal changes in infection intensity were driven by temperature-induced changes in Bd growth rate. Using our model, we found that the maintenance potential of leopard frogs varied significantly throughout the year and that seasonal troughs in infection prevalence made it unlikely that leopard frogs were responsible for long-term Bd persistence in these seasonal amphibian communities, highlighting the importance of alternative pathogen reservoirs for Bd persistence. Our results have broad implications for management in seasonal host-pathogen systems, showing that seasonal changes in host and pathogen vital rates, rather than the depletion of susceptible hosts, can lead to troughs in pathogen prevalence and stochastic pathogen extirpation.

Evaluating environmental DNA as a tool for detecting an amphibian pathogen using an optimized extraction method.

Environmental DNA (eDNA) detection is a valuable conservation tool that can be used to identify and monitor imperiled or invasive species and wildlife pathogens. Batrachochytrium pathogens are of global conservation concern because they are a leading cause of amphibian decline. While eDNA techniques have been used to detect Batrachochytrium DNA in the environment, a systematic comparison of extraction methods across environmental samples is lacking. In this study, we first compared eDNA extraction methods and found that a soil extraction kit (Qiagen PowerSoil) was the most effective for detecting Batrachochytrium dendrobatidis in water samples. The PowerSoil extraction had a minimum detection level of 100 zoospores and had a two- to four-fold higher detection probability than other commonly used extraction methods (e.g., QIAamp extraction, DNeasy+Qiashredder extraction method, respectively). Next, we used this extraction method on field-collected water and sediment samples and were able to detect pathogen DNA in both. While field-collected water filters were equivalent to amphibian skin swab samples in detecting the presence of pathogen DNA, the seasonal patterns in pathogen quantity were different between skin swabs and water samples. Detection rate was lowest in sediment samples. We also found that detection probability increases with the volume of water filtered. Our results indicate that water filter eDNA samples can be accurate in detecting pathogen presence at the habitat scale but their utility for quantifying pathogen loads in the environment appears limited. We suggest that eDNA techniques be used for early warning detection to guide animal sampling efforts.

Optimized Batrachochytrium dendrobatidis DNA extraction of swab samples results in imperfect detection particularly when infection intensities are low.

Accurate detection of the amphibian fungal pathogen Batrachochytrium dendrobatidis (Bd) is critical for wildlife disease research; however, false negatives in detection do occur. Here we compared different DNA extraction methods to determine the threshold for Bd detection and identify an optimal extraction method to improve detection and quantification of the pathogen. We extracted both lab-created cell suspension standards using PrepMan Ultra, Chelex resin, and 3 spin column DNA extraction kits (Qiagen DNeasy Blood and Tissue, Zymo Quick DNA miniprep, and IBI gMAX mini kit), and further compared extraction methods using field-collected samples. We found that when extracting Bd DNA from cells in lab-created culture, the spin column extraction methods and PrepMan Ultra were equivalent, while the resin method detected higher Bd DNA quantities, especially at higher loads. However, when swabs from live animals were analyzed, low Bd quantities were more than twice as likely to be detected using a spin column extraction than with the PrepMan Ultra extraction method. All tested spin column extraction methods performed similarly across both field and lab samples. Samples containing low Bd quantities yielded inconsistent detection and quantification of Bd DNA copies regardless of extraction method. To manage imperfect detection of Bd, we suggest that presence/absence analyses are more informative than attempting to quantify Bd DNA when quantities are low. Overall, we recommend that a cost-benefit analysis of target species susceptibility and epidemiology be taken into consideration when designing an experiment to determine the most appropriate DNA extraction method to be used, because sometimes detecting low Bd quantities is imperative to the study, whereas in other situations, detecting low DNA quantities is less important.

Defectlike structures and localized patterns in the cubic-quintic-septic Swift-Hohenberg equation.

We study numerically the cubic-quintic-septic Swift-Hohenberg (SH357) equation on bounded one-dimensional domains. Under appropriate conditions stripes with wave number k≈1 bifurcate supercritically from the zero state and form S-shaped branches resulting in bistability between small and large amplitude stripes. Within this bistability range we find stationary heteroclinic connections or fronts between small and large amplitude stripes, and demonstrate that the associated spatially localized defectlike structures either snake or fall on isolas. In other parameter regimes we also find heteroclinic connections to spatially homogeneous states and a multitude of dynamically stable steady states consisting of patches of small and large amplitude stripes with different wave numbers or of spatially homogeneous patches. The SH357 equation is thus extremely rich in the types of patterns it exhibits. Some of the features of the bifurcation diagrams obtained by numerical continuation can be understood using a conserved quantity, the spatial Hamiltonian of the system.

Pattern analysis in a benthic bacteria-nutrient system.

We study steady states in a reaction-diffusion system for a benthic bacteria-nutrient model in a marine sediment over 1D and 2D domains by using Landau reductions and numerical path following methods. We point out how the system reacts to changes of the strength of food supply and ingestion. We find that the system has a stable homogeneous steady state for relatively large rates of food supply and ingestion, while this state becomes unstable if one of these rates decreases and Turing patterns such as hexagons and stripes start to exist. One of the main results of the present work is a global bifurcation diagram for solutions over a bounded 2D domain. This bifurcation diagram includes branches of stripes, hexagons, and mixed modes. Furthermore, we find a number of snaking branches of stationary states, which are spatial connections between homogeneous states and hexagons, homogeneous states and stripes as well as stripes and hexagons in parameter ranges, where both corresponding states are stable. The system under consideration originally contains some spatially varying coefficients and with these exhibits layerings of patterns. The existence of spatial connections between different steady states in bistable ranges shows that spatially varying patterns are not necessarily due to spatially varying coefficients. The present work gives another example, where these effects arise and shows how the analytical and numerical observations can be used to detect signs that a marine bacteria population is in danger to die out or on its way to recovery, respectively. We find a type of hexagon patches on a homogeneous background, which seems to be new discovery. We show the first numerically calculated solution-branch, which connects two different types of hexagons in parameter space. We check numerically for bounded domains whether the stability changes for hexagons and stripes, which are extended homogeneously into the third dimension. We find that stripes and one type of hexagons have the same stable range over bounded 2D and 3D domains. This does not hold for the other type of hexagons. Their stable range is shorter for the bounded 3D domain, which we used here. We find a snaking branch, which bifurcates when the hexagonal prisms loose their stability. Solutions on this branch connects spatially between hexagonal prisms and a genuine 3D pattern (balls).

Heterozygosity predicts clutch and egg size but not plasticity in a house sparrow population with no evidence of inbreeding.

We investigated the link between heterozygosity and the reaction norm attributes of reproductive performance in female house sparrows (Passer domesticus). We collected data on clutch size, egg size, hatching success and nestling survival in 2816 nesting attempts made by 791 marked individuals over a 16-year period. Pedigree analysis revealed no evidence of inbreeding. Neither parent-offspring regression nor an animal model revealed significant heritability in clutch or egg size. We selected 42 females that laid at least seven clutches at our study site and used a survey of 21 autosomal microsatellite loci to estimate heterozygosity for each female. We controlled for phenotypic plasticity and found that both clutch and egg size showed significant positive correlations with heterozygosity. We found no evidence that heterozygosity influenced the slope of individual reaction norms. Further analysis suggested that clutch size was affected by heterozygosity across the genome, but egg size had more complex relationships, with evidence favouring the influence of multiple loci. Given the apparent lack of inbreeding and large population size, our results suggest associative overdominance as the likely mechanism for the impact of heterozygosity, but also created a puzzle about the process producing associations between neutral markers and the genes affecting clutch size or egg size. One possible explanation is a long-term residual effect of the historical bottleneck that occurred when house sparrows were introduced into North America. The existence of heterozygosity-fitness correlations in a population with considerable phenotypic plasticity and little inbreeding implies that the effects of heterozygosity may be more significant than previously thought.

Individual variation in parental care reaction norms: integration of personality and plasticity.

Personality (consistent differences between individuals in behavior) and plasticity (changes within individuals in behavior) are often viewed as separate and opposing phenomena. We tested this idea by analyzing parental care reaction norms in a bird that exhibits biparental care. Personality in provisioning behavior existed (r(ic) = 0.11) and persisted despite being reduced after accounting for individual differences in environment. Plasticity was also evident and differed between the sexes. Male visit rate was associated with changes in brood size and time of day, but female visit rate was associated with changes in nestling age and date. In both sexes changes in visit rate were positively correlated with changes in their partner's visit rate. Both sexes also exhibited multidimensional reaction norms; interaction terms revealed that within-individual visit rates increased more steeply with brood size when nestlings were older, and the effect of the partner's visit rate was sensitive to variation in date, precipitation, and the focal bird's age. Individuals also varied in how they responded (reaction norm slope) to changes in nestling age and partner visits. Moreover, parental personality was interdependent with individual plasticity in several ways. Individuals of both sexes with a high visit rate also responded more positively to changes in nestling age, and males also showed this pattern with changes in partner visit rate. Explicit use of the behavioral reaction norm integrated personality and plasticity, revealed that these are not opposing concepts, and stimulated new hypotheses about sexual conflict over care and provisioning as a life-history trait.

Heterozygosity and extra-pair paternity: biased tests result from the use of shared markers.

Recent studies of extra-pair paternity have found support for the idea that heterozygous males have an advantage in siring offspring. Most studies use DNA microsatellite loci to determine paternity and then use the same loci to estimate individual heterozygosity. However, because the likelihood of detecting extra-pair offspring depends on the combinations of parental alleles, it is possible that biases arise from particular allele combinations. This might produce false support for the influence of heterozygosity on mating behaviour. We used a simulation model to assess how large this bias might be. We found two sources of bias. First, we found a bias in the null hypothesis of a simple statistical test commonly used to test several predictions of the heterozygosity hypothesis. The use of randomization tests could eliminate this bias. Second, we found that using the same loci for both paternity and heterozygosity can cause an increase in results supporting the heterozygosity hypothesis when no effect of heterozygosity actually exists. This bias is reduced through the use of more markers with higher levels of polymorphism and heterozygosity, but can be eliminated entirely by using a separate set of markers to determine paternity and assess heterozygosity. The two sources of bias reduce evidence favouring the heterozygosity hypothesis, but do not negate all of the studies that support it. We suggest that further studies of heterozygosity and extra-pair paternity are important and likely to be informative, but our recommendations should be incorporated by researchers to improve the reliability of their conclusions.