Phenological synchrony shapes pathology in host-parasite systems

Broadscale spatial synchrony in a West Nile virus mosquito vector across multiple timescales

Insects often exhibit irruptive population dynamics determined by environmental conditions. We examine if populations of the Culex tarsalis mosquito, a West Nile virus (WNV) vector, fluctuate synchronously over broad spatial extents and multiple timescales and whether climate drives synchrony in Cx. tarsalis, especially at annual timescales, due to the synchronous influence of temperature, precipitation, and/or humidity. We leveraged mosquito collections across 9 National Ecological Observatory Network (NEON) sites distributed in the interior West and Great Plains region USA over a 45-month period, and associated gridMET climate data. We utilized wavelet phasor mean fields and wavelet linear models to quantify spatial synchrony for mosquitoes and climate and to calculate the importance of climate in explaining Cx. tarsalis synchrony. We also tested whether the strength of spatial synchrony may vary directionally across years. We found significant annual synchrony in Cx. tarsalis, and short-term synchrony during a single period in 2018. Mean minimum temperature was a significant predictor of annual Cx. tarsalis spatial synchrony, and we found a marginally significant decrease in annual Cx. tarsalis synchrony. Significant Cx. tarsalis synchrony during 2018 coincided with an anomalous increase in precipitation. This work provides a valuable step toward understanding broadscale synchrony in a WNV vector.

Into the danger zone: How the within-host distribution of parasites controls virulence

Despite the importance of virulence in epidemiological theory, the relative contributions of host and parasite to virulence outcomes remain poorly understood. Here, we use reciprocal cross experiments to disentangle the influence of host and parasite on core virulence components-infection and pathology-and understand dramatic differences in parasite-induced malformations in California amphibians. Surveys across 319 populations revealed that amphibians' malformation risk was 2.7× greater in low-elevation ponds, even while controlling for trematode infection load. Factorial experiments revealed that parasites from low-elevation sites induced higher per-parasite pathology (reduced host survival and growth), whereas there were no effects of host source on resistance or tolerance. Parasite populations also exhibited marked differences in within-host distribution: ~90% of low-elevation cysts aggregated around the hind limbs, relative to <60% from high-elevation. This offers a novel, mechanistic basis for regional variation in parasite-induced malformations while promoting a framework for partitioning host and parasite contributions to virulence.

Parasite-mediated selection on host phenology

The timing of seasonal activity, or phenology, is an adaptive trait that maximizes individual fitness by timing key life events to coincide with favorable abiotic factors and biotic interactions. Studies on the biotic interactions that determine optimal phenology have focused on temporal overlaps among positively-interacting species such as mutualisms. Less well understood is the extent that negative interactions such as parasitism impact the evolution of host phenology. Here, we present a mathematical model demonstrating the evolution of host phenological patterns in response to sterilizing parasites. Environments with parasites favor hosts with shortened activity periods or greater distributions in emergence timing, both of which reduce the temporal overlap between hosts and parasites and thus reduce infection risk. Although host populations with these altered phenological patterns are less likely to mature and reproduce, the fitness advantage of parasite avoidance can be greater than the cost of reduced reproduction. These results illustrate the impact of parasitism on the evolution of host phenology and suggest that shifts in host phenology could serve as a strategy to mitigate the risk of infection.

A reconstruction of parasite burden reveals one century of climate-associated parasite decline

Long-term data allow ecologists to assess trajectories of population abundance. Without this context, it is impossible to know whether a taxon is thriving or declining to extinction. For parasites of wildlife, there are few long-term data-a gap that creates an impediment to managing parasite biodiversity and infectious threats in a changing world. We produced a century-scale time series of metazoan parasite abundance and used it to test whether parasitism is changing in Puget Sound, United States, and, if so, why. We performed parasitological dissection of fluid-preserved specimens held in natural history collections for eight fish species collected between 1880 and 2019. We found that parasite taxa using three or more obligately required host species-a group that comprised 52% of the parasite taxa we detected-declined in abundance at a rate of 10.9% per decade, whereas no change in abundance was detected for parasites using one or two obligately required host species. We tested several potential mechanisms for the decline in 3+-host parasites and found that parasite abundance was negatively correlated with sea surface temperature, diminishing at a rate of 38% for every 1 °C increase. Although the temperature effect was strong, it did not explain all variability in parasite burden, suggesting that other factors may also have contributed to the long-term declines we observed. These data document one century of climate-associated parasite decline in Puget Sound-a massive loss of biodiversity, undetected until now.

Contrasting temperature responses in seasonal timing of cercariae shedding by Rhipidocotyle trematodes

Global warming is likely to lengthen the seasonal duration of larval release by parasites. We exposed freshwater mussel hosts, Anodonta anatina, from 2 high-latitude populations to high, intermediate and low temperatures throughout the annual cercarial shedding period of the sympatric trematodes Rhipidocotyle fennica and R. campanula, sharing the same transmission pathway. At the individual host level, under warmer conditions, the timing of the cercarial release in both parasite species shifted towards seasonally earlier period while its duration did not change. At the host population level, evidence for the lengthening of larvae shedding period with warming was found for R. fennica. R. campanula started the cercarial release seasonally clearly earlier, and at a lower temperature, than R. fennica. Furthermore, the proportion of mussels shedding cercariae increased, while day-degrees required to start the cercariae shedding decreased in high-temperature treatment in R. fennica. In R. campanula these effects were not found, suggesting that warming can benefit more R. fennica. These results do not completely support the view that climate warming would invariably increase the seasonal duration of larval shedding by parasites, but emphasizes species-specific differences in temperature-dependence and in seasonality of cercarial release.

Strigea robusta (Digenea: Strigeidae) infection effects on the gonadal structure and limb malformation in toad early development

The anomaly P is a mass morphological anomaly reported in some water frog populations across Europe. It was found that polydactyly is only a mild attenuated form of heavy cases of the anomaly P syndrome, which have strong deformations of the hindlimbs and, partly, forelimbs. It was shown that the anomaly P is caused by the trematode Strigea robusta and this syndrome can be considered as a special case of strigeosis in amphibians. The anomaly P for a long time considered to be specific for water frogs of the genus Pelophylax. Herein, we describe polydactyly and heavy forms of the anomaly P syndrome in toads of the genera Bufo and Bufotes, as a result of exposure to S. robusta cercariae. A total of 150 tadpoles of Bufo bufo, 60 tadpoles of Bufotes viridis, and 60 tadpoles of Bufotes baturae were divided into five experimental and four control groups (30 tadpoles in each group). All anomalies in the toads were similar to those observed in water frogs. The survival of tadpoles in the experimental groups was 76%. The anomaly P was observed in 57.9% of toad tadpoles (51.8% of mild forms and 6.1% of heavy forms). The occurrence of the anomaly P varied among groups from 19% to 78%. Heavy forms of the anomaly P were found in all experimental groups. We described rare asymmetrical cases of the anomaly P. According to severe modification of limb morphology, we supposed changes of gonadal morphology (any modifications of the germ and somatic cells). The gonadal development of infected tadpoles was however the same as in uninfected toad tadpoles, and heterochromatin distribution within gonocytes had no differences as well. It seems like the parasite doesn't have any effect on the gonadal development of the toads. The lack of heavy forms in natural populations of toads, as well as a development of gonads were discussed.

Host phenology can drive the evolution of intermediate virulence strategies in some obligate-killer parasites

Traditional mechanistic trade-offs between transmission and virulence are the foundation of nearly all theory on parasite virulence evolution. For obligate-host killer parasites, evolution toward intermediate virulence depends on a trade-off between virulence (time to death) and transmission (the number of progeny released upon death). Although several ecological factors impact optimal virulence strategies constrained by trade-offs, these factors have been insufficient to explain the intermediate virulence levels observed in nature. The timing of seasonal activity, or phenology, is a factor that commonly influences ecological interactions but is difficult to incorporate into virulence evolution studies. We present a mathematical model of a seasonal obligate-killer parasite to study the impact of host phenology on virulence evolution. The model demonstrates that host phenology can select for intermediate parasite virulence even when a traditional mechanistic trade-off between transmission and virulence is omitted. The optimal virulence strategy is impacted by both the host activity period duration and the host emergence timing variation. Parasites with lower virulence strategies are favored in environments with longer host activity periods and when hosts emerge synchronously. The results demonstrate that host phenology can be sufficient to select for intermediate virulence strategies, providing an alternative driver of virulence evolution in some natural systems.

Quantifying phenological diversity: a framework based on Hill numbers theory

Background

Despite the great concern triggered by the environmental crisis worldwide, the loss of temporal key functions and processes involved in biodiversity maintenance has received little attention. Species are restricted in their life cycles by environmental variables because of their physiological and behavioral properties; thus, the timing and duration of species' presence and their activities vary greatly between species within a community. Despite the ecological relevance of such variation, there is currently no measure that summarizes the key temporal aspects of biological diversity and allows comparisons of community phenological patterns. Here, we propose a measure that synthesizes variability of phenological patterns using the Hill numbers-based attribute diversity framework.

Methods

We constructed a new phenological diversity measure based on the aforementioned framework through pairwise overlapping distances, which was supplemented with wavelet analysis. The Hill numbers approach was chosen as an adequate way to define a set of diversity values of different order q, a parameter that determines the sensitivity of the diversity measure to abundance. Wavelet transform analysis was used to model continuous variables from incomplete data sets for different phenophases. The new measure, which we call Phenological Hill numbers (PD), considers the decouplings of phenophases through an overlapping area value between pairs of species within the community. PD was first tested through simulations with varying overlap in phenophase magnitude and intensity and varying number of species, and then by using one real data set.

Results

PD maintains the diversity patterns of order q as in any other diversity measure encompassed by the Hill numbers framework. Minimum PD values in the simulated data sets reflect a lack of differentiation in the phenological curves of the community over time; by contrast, the maximum PD values reflected the most diverse simulations in which phenological curves were equally distributed over time. PD values were consistent with the homogeneous distribution of the intensity and concurrence of phenophases over time, both in the simulated and the real data set.

Discussion

PD provides an efficient, readily interpretable and comparable measure that summarizes the variety of phenological patterns observed in ecological communities. PD retains the diversity patterns of order q characteristic of all diversity measures encompassed by the distance-based Hill numbers framework. In addition, wavelet transform analysis proved useful for constructing a continuous phenological curve. This methodological approach to quantify phenological diversity produces simple and intuitive values for the examination of phenological diversity and can be widely applied to any taxon or community's phenological traits.

Acute mortality in California tiger salamander (Ambystoma californiense) and Santa Cruz long-toed salamander (Ambystoma macrodactylum croceum) caused by Ribeiroia ondatrae (Class: Trematoda)

In early September 2019, a morbidity and mortality event affecting California tiger salamanders (Ambystoma californiense) and Santa Cruz long-toed salamanders (Ambystoma macrodactylum croceum) in late stages of metamorphosis was reported at a National Wildlife Refuge in Santa Cruz County, California, U.S.A. During the postmortem disease investigation, severe integumentary metacercarial (Class: Trematoda) infection, associated with widespread skin lesions, was observed. Planorbid snails collected from the ponds of the refuge within seven days of the mortality event were infected with Ribeiroia ondatrae, a digenetic trematode that can cause malformation and death in some amphibians. We suggest sustained seasonal high-water levels due to active habitat management along with several years of increased rainfall led to increased bird visitation, increased over-wintering of snails, and prolonged salamander metamorphosis, resulting in a confluence of conditions and cascading of host-parasite dynamics to create a hyper-parasitized state.

•

Mortality event in wild endangered salamanders in California.

•

Infection with Ribeiroia ondatrae caused severe fatal skin lesions.

•

Sustaining water levels may increase parasite transmission.

Phenological and intrinsic predictors of mite and haemacoccidian infection dynamics in a Mediterranean community of lizards

Ectotherms are vulnerable to environmental changes and their parasites are biological health indicators. Thus, parasite load in ectotherms is expected to show a marked phenology. This study investigates temporal host–parasite dynamics in a lizard community in Eastern Spain during an entire annual activity period. The hosts investigated were Acanthodactylus erythrurus, Psammodromus algirus and Psammodromus edwardsianus, three lizard species coexisting in a mixed habitat of forests and dunes, providing a range of body sizes, ecological requirements and life history traits. Habitat and climate were considered as potential environmental predictors of parasite abundance, while size, body condition and sex as intrinsic predictors. Linear models based on robust estimates were fitted to analyse parasite abundance and prevalence. Ectoparasitic mites and blood parasites from two haemococcidian genera were found: Lankesterella spp. and Schellackia spp. Habitat type was the only predictor explaining the abundance of all parasites, being mostly higher in the forest than in the dunes. The results suggest that particularities in each host–parasite relationship should be accounted even when parasites infect close-related hosts under the same environmental pressures. They also support that lizard parasites can be biomarkers of environmental perturbation, but the relationships need to be carefully interpreted for each host–parasite assemblage.

Disease's hidden death toll: Using parasite aggregation patterns to quantify landscape-level host mortality in a wildlife system

World-wide, infectious diseases represent a major source of mortality in humans and livestock. For wildlife populations, disease-induced mortality is likely even greater, but remains notoriously difficult to estimate-especially for endemic infections. Approaches for quantifying wildlife mortality due to endemic infections have historically been limited by an inability to directly observe wildlife mortality in nature. Here we address a question that can rarely be answered for endemic pathogens of wildlife: what are the population- and landscape-level effects of infection on host mortality? We combined laboratory experiments, extensive field data and novel mathematical models to indirectly estimate the magnitude of mortality induced by an endemic, virulent trematode parasite (Ribeiroia ondatrae) on hundreds of amphibian populations spanning four native species. We developed a flexible statistical model that uses patterns of aggregation in parasite abundance to infer host mortality. Our model improves on previous approaches for inferring host mortality from parasite abundance data by (i) relaxing restrictive assumptions on the timing of host mortality and sampling, (ii) placing all mortality inference within a Bayesian framework to better quantify uncertainty and (iii) accommodating data from laboratory experiments and field sampling to allow for estimates and comparisons of mortality within and among host populations. Applying our approach to 301 amphibian populations, we found that trematode infection was associated with an average of between 13% and 40% population-level mortality. For three of the four amphibian species, our models predicted that some populations experienced >90% mortality due to infection, leading to mortality of thousands of amphibian larvae within a pond. At the landscape scale, the total number of amphibians predicted to succumb to infection was driven by a few high mortality sites, with fewer than 20% of sites contributing to greater than 80% of amphibian mortality on the landscape. The mortality estimates in this study provide a rare glimpse into the magnitude of effects that endemic parasites can have on wildlife populations and our theoretical framework for indirectly inferring parasite-induced mortality can be applied to other host-parasite systems to help reveal the hidden death toll of pathogens on wildlife hosts.

Phenological synchrony shapes pathology in host-parasite systems

A key challenge surrounding ongoing climate shifts is to identify how they alter species interactions, including those between hosts and parasites. Because transmission often occurs during critical time windows, shifts in the phenology of either taxa can alter the likelihood of interaction or the resulting pathology. We quantified how phenological synchrony between vulnerable stages of an amphibian host (Pseudacris regilla) and infection by a pathogenic trematode (Ribeiroia ondatrae) determined infection prevalence, parasite load and host pathology. By tracking hosts and parasite infection throughout development between low- and high-elevation regions (San Francisco Bay Area and the Southern Cascades (Mt Lassen)), we found that when phenological synchrony was high (Bay Area), each established parasite incurred a 33% higher probability of causing severe limb malformations relative to areas with less synchrony (Mt Lassen). As a result, hosts in the Bay Area had up to a 50% higher risk of pathology even while controlling for the mean infection load. Our results indicate that host-parasite interactions and the resulting pathology were the joint product of infection load and phenological synchrony, highlighting the sensitivity of disease outcomes to forecasted shifts in climate.