Experimental investigation of alternative transmission functions: Quantitative evidence for the importance of nonlinear transmission dynamics in host-parasite systems

Diverging effects of host density and richness across biological scales drive diversity-disease outcomes

Understanding how biodiversity affects pathogen transmission remains an unresolved question due to the challenges in testing potential mechanisms in natural systems and how these mechanisms vary across biological scales. By quantifying transmission of an entire guild of parasites (larval trematodes) within 902 amphibian host communities, we show that the community-level drivers of infection depend critically on biological scale. At the individual host scale, increases in host richness led to fewer parasites per host for all parasite taxa, with no effect of host or predator densities. At the host community scale, however, the inhibitory effects of richness were counteracted by associated increases in total host density, leading to no overall change in parasite densities. Mechanistically, we find that while average host competence declined with increasing host richness, total community competence remained stable due to additive assembly patterns. These results help reconcile disease-diversity debates by empirically disentangling the roles of alternative ecological drivers of parasite transmission and how such effects depend on biological scale.

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.

Frequency-dependent transmission of Batrachochytrium salamandrivorans in eastern newts

Transmission is the fundamental process whereby pathogens infect their hosts and spread through populations, and can be characterized using mathematical functions. The functional form of transmission for emerging pathogens can determine pathogen impacts on host populations and can inform the efficacy of disease management strategies. By directly measuring transmission between infected and susceptible adult eastern newts (Notophthalmus viridescens) in aquatic mesocosms, we identified the most plausible transmission function for the emerging amphibian fungal pathogen Batrachochytrium salamandrivorans (Bsal). Although we considered a range of possible transmission functions, we found that Bsal transmission was best explained by pure frequency dependence. We observed that >90% of susceptible newts became infected within 17 days post‐exposure to an infected newt across a range of host densities and initial infection prevalence treatments. Under these conditions, we estimated R₀ = 4.9 for Bsal in an eastern newt population. Our results suggest that Bsal has the capability of driving eastern newt populations to extinction and that managing host density may not be an effective management strategy. Intervention strategies that prevent Bsal introduction or increase host resistance or tolerance to infection may be more effective. Our results add to the growing empirical evidence that transmission of wildlife pathogens can saturate and be functionally frequency‐dependent.

How parasite exposure and time interact to determine Australapatemon burti (Trematoda: Digenea) infections in second intermediate hosts (Erpobdella microstoma) (Hirudinea: Erpodellidae)

Australapatemon spp. are cosmopolitan trematodes that infect freshwater snails, aquatic leeches, and birds. Despite their broad geographic distribution, relatively little is known about interactions between Australapatemon spp. and their leech hosts, particularly under experimental conditions and in natural settings. We used experimental exposures to determine how Australapatemon burti cercariae dosage (number administered to leech hosts, Erpobdella microstoma) affected infection success (fraction to encyst as metacercariae), infection abundance, host survival, and host size over the 100 days following exposure. Interestingly, infection success was strongly density-dependent, such that there were no differences in metacercariae load even among hosts exposed to a 30-fold difference in cercariae. This relationship suggests that local processes (e.g., resource availability, interference competition, or host defenses) may play a strong role in parasite transmission. Our results also indicated that metacercariae did not become evident until ~4 weeks post exposure, with average load climbing until approximately 13 weeks. There was no evidence of metacercariae death or clearance over the census period. Parasite exposure had no detectable effects on leech size or survival, even with nearly 1,000 cercariae. Complementary surveys of leeches in California revealed that 11 of 14 ponds supported infection by A. burti (based on morphology and molecular sequencing), with an average prevalence of 32% and similar metacercariae intensity as in our experimental exposures. The extended development time and extreme density dependence of A. burti has implications for studying naturally occurring host populations, for which detected infections may represent only a fraction of cercariae to which animals have been exposed. Future investigation of these underlying mechanisms would be benefical in understanding host-parasite relationships.

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.

Patterns of Clinostomum marginatum infection in fishes and amphibians: integration of field, genetic, and experimental approaches

Digenetic trematodes of the genus Clinostomum are cosmopolitan parasites infecting fishes, amphibians, reptiles, and snails as intermediate hosts. Despite the broad geographical distribution of this genus, debate about the number of species and how they vary in host use has persisted. To better understand patterns of infection among host species and across life stages, we used large-scale field surveys and molecular tools to examine five species of amphibians and seven species of fishes from 125 California ponds. Among the 12,360 examined hosts, infection was rare, with an overall prevalence of 1.7% in amphibians and 9.2% in fishes. Molecular evidence indicated that both groups were infected with Clinostomum marginatum. Using generalized linear mixed effects models, host species identity and host life stage had a strong influence on infection status, such that Lepomis cyanellus (green sunfish) (49.3%) and Taricha granulosa (rough skinned newt) (9.2%) supported the highest overall prevalence values, whereas adult amphibians tended to have a higher prevalence of infection relative to juveniles (13.3% and 2.5%, respectively). Experimentally, we tested the susceptibility of two amphibian hosts (Pseudacris regilla [Pacific chorus frog] and Anaxyrus boreas [western toad]) to varying levels of cercariae exposure and measured metacercariae growth over time. Pseudacris regilla was 1.3× more susceptible to infection, while infection success increased with cercariae exposure dose for both species. On average, metacarcariae size increased by 650% over 20 days. Our study highlights the importance of integrating field surveys, genetic tools, and experimental approaches to better understand the ecology of host-parasite interactions.

The human-snail transmission environment shapes long term schistosomiasis control outcomes: Implications for improving the accuracy of predictive modeling

INTRODUCTION

Schistosomiasis is a chronic parasitic trematode disease that affects over 240 million people worldwide. The Schistosoma lifecycle is complex, involving transmission via specific intermediate-host freshwater snails. Predictive mathematical models of Schistosoma transmission have often chosen to simplify or ignore the details of environmental human-snail interaction in their analyses. Schistosome transmission models now aim to provide better precision for policy planning of elimination of transmission. This heightens the importance of including the environmental complexity of vector-pathogen interaction in order to make more accurate projections.

METHODOLOGY AND PRINCIPAL FINDINGS

We propose a nonlinear snail force of infection (FOI) that takes into account an intermediate larval stage (miracidium) and snail biology. We focused, in particular, on the effects of snail force of infection (FOI) on the impact of mass drug administration (MDA) in human communities. The proposed (modified) model was compared to a conventional model in terms of their predictions. A longitudinal dataset generated in Kenya field studies was used for model calibration and validation. For each sample community, we calibrated modified and conventional model systems, then used them to model outcomes for a range of MDA regimens. In most cases, the modified model predicted more vigorous post-MDA rebound, with faster relapse to baseline levels of infection. The effect was pronounced in higher risk communities. When compared to observed data, only the modified system was able to successfully predict persistent rebound of Schistosoma infection.

CONCLUSION AND SIGNIFICANCE

The observed impact of varying location-specific snail inputs sheds light on the diverse MDA response patterns noted in operational research on schistosomiasis control, such as the recent SCORE project. Efficiency of human-to-snail transmission is likely to be much higher than predicted by standard models, which, in practice, will make local elimination by implementation of MDA alone highly unlikely, even over a multi-decade period.

Experimental investigation of alternative transmission functions: Quantitative evidence for the importance of nonlinear transmission dynamics in host-parasite systems

Understanding pathogen transmission is crucial for predicting and managing disease. Nonetheless, experimental comparisons of alternative functional forms of transmission remain rare, and those experiments that are conducted are often not designed to test the full range of possible forms.

To differentiate among 10 candidate transmission functions, we used a novel experimental design in which we independently varied four factors—duration of exposure, numbers of parasites, numbers of hosts and parasite density—in laboratory infection experiments.

We used interactions between amphibian hosts and trematode parasites as a model system and all candidate models incorporated parasite depletion. An additional manipulation involving anaesthesia addressed the effects of host behaviour on transmission form.

Across all experiments, nonlinear transmission forms involving either a power law or a negative binomial function were the best‐fitting models and consistently outperformed the linear density‐dependent and density‐independent functions. By testing previously published data for two other host–macroparasite systems, we also found support for the same nonlinear transmission forms.

Although manipulations of parasite density are common in transmission studies, the comprehensive set of variables tested in our experiments revealed that variation in density alone was least likely to differentiate among competing transmission functions. Across host–pathogen systems, nonlinear functions may often more accurately represent transmission dynamics and thus provide more realistic predictions for infection.