Disease hotspots or hot species? Infection dynamics in multi-host metacommunities controlled by species identity, not source location

Contribution of host species and pathogen clade to snake fungal disease hotspots in Europe

Infectious diseases are influenced by interactions between host and pathogen, and the number of infected hosts is rarely homogenous across the landscape. Areas with elevated pathogen prevalence can maintain a high force of infection and may indicate areas with disease impacts on host populations. However, isolating the ecological processes that result in increases in infection prevalence and intensity remains a challenge. Here we elucidate the contribution of pathogen clade and host species in disease hotspots caused by Ophidiomyces ophidiicola, the pathogen responsible for snake fungal disease, in 21 species of snakes infected with multiple pathogen strains across 10 countries in Europe. We found isolated areas of disease hotspots in a landscape where infections were otherwise low. O. ophidiicola clade had important effects on transmission, and areas with multiple pathogen clades had higher host infection prevalence. Snake species further influenced infection, with most positive detections coming from species within the Natrix genus. Our results suggest that both host and pathogen identity are essential components contributing to increased pathogen prevalence.

Energy input, habitat heterogeneity and host specificity drive avian haemosporidian diversity at continental scales

The correct identification of variables affecting parasite diversity and assemblage composition at different spatial scales is crucial for understanding how pathogen distribution responds to anthropogenic disturbance and climate change. Here, we used a database of avian haemosporidian parasites to test how the taxonomic and phylogenetic diversity and phylogenetic structure of the genera Plasmodium, Haemoproteus and Leucocytozoon from three zoogeographic regions are related to surrogate variables of Earth's energy input, habitat heterogeneity (climatic diversity, landscape heterogeneity, host richness and human disturbance) and ecological interactions (resource use), which was measured by a novel assemblage-level metric related to parasite niche overlap (degree of generalism). We found that different components of energy input explained variation in richness for each genus. We found that human disturbance influences the phylogenetic structure of Haemoproteus while the degree of generalism explained richness and phylogenetic structure of Plasmodium and Leucocytozoon genera. Furthermore, landscape attributes related to human disturbance (human footprint) can filter Haemoproteus assemblages by their phylogenetic relatedness. Finally, assembly processes related to resource use within parasite assemblages modify species richness and phylogenetic structure of Plasmodium and Leucocytozoon assemblages. Overall, our study highlighted the genus-specific patterns with the different components of Earth's energy budget, human disturbances and degree of generalism.

Pseudacris regilla metamorphs acquire resistance to Batrachochytrium dendrobatidis after exposure to the killed fungus

The pathogenic fungus Batrachochytrium dendrobatidis (Bd) is associated with drastic global amphibian declines. Prophylactic exposure to killed zoospores and the soluble chemicals they produce (Bd metabolites) can induce acquired resistance to Bd in adult Cuban treefrogs Osteopilus septentrionalis. Here, we exposed metamorphic frogs of a second species, the Pacific chorus frog Pseudacris regilla, to one of 2 prophylactic treatments prior to live Bd exposures: killed Bd zoospores with metabolites, killed zoospores alone, or a water control. Prior exposure to killed Bd zoospores with metabolites reduced Bd infection intensity in metamorphic Pacific chorus frogs by 60.4% compared to control frogs. Interestingly, Bd intensity in metamorphs previously exposed to killed zoospores alone did not differ in magnitude relative to the control metamorphs, nor to those treated with killed zoospores plus metabolites. Previous work indicated that Bd metabolites alone can induce acquired resistance in tadpoles, and so these findings together indicate that it is possible that the soluble Bd metabolites may contain immunomodulatory components that drive this resistance phenotype. Our results expand the generality of this prophylaxis work by identifying a second amphibian species (Pacific chorus frog) and an additional amphibian life stage (metamorphic frog) that can acquire resistance to Bd after metabolite exposure. This work increases hopes that a Bd-metabolite prophylaxis might be widely effective across amphibian species and life stages.

Wildlife susceptibility to infectious diseases at global scales

Disease transmission prediction across wildlife is crucial for risk assessment of emerging infectious diseases. Susceptibility of host species to pathogens is influenced by the geographic, environmental, and phylogenetic context of the specific system under study. We used machine learning to analyze how such variables influence pathogen incidence for multihost pathogen assemblages, including one of direct transmission (coronaviruses and bats) and two vector-borne systems (West Nile Virus [WNV] and birds, and malaria and birds). Here we show that this methodology is able to provide reliable global spatial susceptibility predictions for the studied host-pathogen systems, even when using a small amount of incidence information (i.e., [Formula: see text] of information in a database). We found that avian malaria was mostly affected by environmental factors and by an interaction between phylogeny and geography, and WNV susceptibility was mostly influenced by phylogeny and by the interaction between geographic and environmental distances, whereas coronavirus susceptibility was mostly affected by geography. This approach will help to direct surveillance and field efforts providing cost-effective decisions on where to invest limited resources.

Host-parasite Interactions between a Copepod (Pharodes tortugensis) and Small Reef-associated Gobies (Coryphopterus) in the British Virgin Islands

The effects of parasitic copepods on free-living hosts are infrequently documented, and the copepod Pharodes tortugensis has remained virtually unstudied since described. For the first time, we document its host range in the British Virgin Islands (BVI), the prevalence and intensity of infections on wild hosts, and its impacts on host morphology and performance. Infections were observed on four benthic gobies in the BVI (Coryphopterus glaucofraenum, C. venezuelae, C. dicrus and C. eidolon) but not on other host species previously reported from other parts of the western Atlantic. Infected gobies were widespread in the BVI (detected at 33 of 52 sites, prevalence from 1-25%) but extremely rare elsewhere in the Caribbean (detected at 2 of 16 sites, prevalence < 0.006%). As is typical of macroparasite infections, P. tortugensis was over-dispersed in BVI host populations (mean intensity = 4.7, range = 1-17). Infections were most common in juvenile and female hosts, and rarely found in larger male hosts. The copepods attach in the branchial chamber of the goby; female copepods show high attachment fidelity to the ventral surface of the chamber, while male copepods attached most often to the first two gill arches and in the branchial chamber adjacent to the female. Infections caused substantial damage to the host's branchial chamber and gill filaments. Parasitized gobies also had larger livers and smaller gonads than unparasitized individuals of similar length. The changes in organ mass of infected gobies were not sizeable enough to affect total body mass, and host condition (the body-length vs. body-mass relationship) was similar for gobies with and without infections. Parasitized gobies were, however, significantly smaller in body mass at a given age, reflecting slower overall growth. Effects of P. tortugensis on individual hosts were broadly similar to those of other parasitic copepods that infect fish gills and, for unknown reasons, the BVI appears to be a persistent hotspot of infections on these goby hosts.

High prevalence does not necessarily equal maintenance species: Avoiding biased claims of disease reservoirs when using surveillance data

1. Many pathogens of public health and conservation concern persist in host communities. Identifying candidate maintenance and reservoir species is therefore a central component of disease management. The term maintenance species implies that if all species but the putative maintenance species were removed, then the pathogen would still persist. In the absence of field manipulations, this statement inherently requires a causal or mechanistic model to assess. 2. However, we lack a systematic understanding of i) how often conclusions are made about maintenance and reservoir species without reference to mechanistic models ii) what types of biases may be associated with these conclusions and iii) how explicitly invoking causal or mechanistic modeling can help ameliorate these biases. Filling these knowledge gaps is critical for robust inference on pathogen persistence and spillover in multihost parasite systems, with clear implications for human and wildlife health. 3. To address these gaps, we performed a literature review on the evidence previous studies have used to make claims regarding maintenance or reservoir species. We then developed a multihost-parasite model to explore and demonstrate common biases that could arise when inferring maintenance potential from observational prevalence data. Finally, we developed new theory to show how model-driven inference of maintenance species can minimize and eliminate emergent biases. 4. In our review, we found that 83% of studies used some form of observational prevalence data to draw conclusions on maintenance potential and only 6% of these studies combined observational data with mechanistic modeling. Using our model, we demonstrate how the community, spatial, and temporal context of observational data can lead to substantial biases in inferences of maintenance potential. Importantly, our theory identifies that model-driven inference of maintenance species elucidates other streams of observational data that can be leveraged to correct these biases. 5. Model-driven inference is an essential, yet underused, component of multidisciplinary studies that make inference on host reservoir and maintenance species. Better integration of wildlife disease surveillance and mechanistic models is necessary to improve the robustness and reproducibility of our conclusions regarding maintenance and reservoir species.

Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis)

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd), termed Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the site of origin for the common ancestor of modern Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swabs and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs (Rana catesbeiana) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential site of origin for Bd-GPL.

Ranavirus Amplification in Low-Diversity Amphibian Communities

In an era where emerging infectious diseases are a serious threat to biodiversity, epidemiological patterns need to be identified, particularly the complex mechanisms driving the dynamics of multi-host pathogens in natural communities. Many amphibian species have faced unprecedented population declines associated with diseases. Yet, specific processes shaping host-pathogen relationships within and among communities for amphibian pathogens such as ranaviruses (RV) remain poorly understood. To address this gap, we conducted a comprehensive study of RV in low-diversity amphibian communities in north-western Canada to assess the effects of biotic factors (species identity, species richness, abundance) and abiotic factors (conductivity, pH) on the pathogen prevalence and viral loads. Across 2 years and 18 sites, with communities of up to three hosts (wood frog, Rana sylvatica; boreal chorus frog, Pseudacris maculata; Canadian toad, Anaxyrus hemiophrys), we observed that RV prevalence nearly doubled with each additional species in a community, suggesting an amplification effect in aquatic, as well as terrestrial life-history stages. Infection intensity among infected wood frogs and boreal chorus frogs also significantly increased with an increase in species richness. Interestingly, we did not observe any effects of host abundance or abiotic factors, highlighting the importance of including host identity and species richness when investigating multi-host pathogens. Ultimately, only such a comprehensive approach can improve our understanding of complex and often highly context-dependent host-pathogen interactions.

Effects of an alternative host on the prevalence and infection intensity of a bumble bee parasite

Several bee parasites are transmitted through flowers, and some of them can infect multiple host species. Given the shared use of flowers by bee species, parasites can potentially encounter multiple host species, which could affect the evolution of parasite virulence. We used the trypanosomatid parasite Crithidia bombi and its host, the common eastern bumble bee (Bombus impatiens), to explore the effect of infecting an alternative host, the alfalfa leaf-cutter bee (Megachile rotundata), on parasite infectivity and ability to replicate. We conducted a serial passage experiment on primary and alternative hosts, assessing infectivity and intensity of infection during five passes. Parasite cells from each pass through the alternative host were also used to infect a group of primary hosts. We found that serial passes through the alternative host increased infectivity, but there was no effect on intensity of infection. Interestingly, both the probability and intensity of infection on the primary host increased after serial passage through the alternative host. This increase in intensity of infection could be due to maladaptation after selection of new C. bombi strains has occurred in the alternative host. This study suggests that host switching has the potential to affect the adaptation of bee parasites to their hosts.

Genetic variation of Batrachochytrium dendrobatidis is linked to skin bacterial diversity in the Pacific treefrog Hyliola regilla (hypochondriaca)

Symbiotic bacterial communities are crucial to combating infections and contribute to host health. The amphibian skin microbiome plays an important role in protecting their hosts against pathogens such as Batrachochytrium dendrobatidis (Bd), one of the causative agents of chytridiomycosis, which is responsible for dramatic amphibian population declines worldwide. Although symbiotic skin bacteria are known to inhibit Bd growth, an understanding of the relationship between Bd genetic variability, environmental conditions, and skin bacterial communities is limited. Therefore, we examined the associations between Bd infection load, Bd genetic diversity and skin bacterial communities in five populations of Hyliola regilla (hypochondriaca) from environmentally contrasting sites in Baja California, Mexico. We observed differences in Bd genetics and infection load among sites and environments. Genetic analysis of Bd isolates revealed patterns of spatial structure corresponding to the five sites sampled. Amphibian skin bacterial diversity and community structure differed among environments and sites. Bacterial community composition was correlated with Bd genetic differences and infection load, with specific bacterial taxa enriched on infected and un-infected frogs. Our results indicate that skin-associated bacteria and Bd strains likely interact on the host skin, with consequences for microbial community structure and Bd infection intensity.

Community-level prevalence of a forest pathogen, not individual-level disease risk, declines with tree diversity

Understanding why diversity sometimes limits disease is essential for managing outbreaks; however, mechanisms underlying this 'dilution effect' remain poorly understood. Negative diversity-disease relationships have previously been detected in plant communities impacted by an emerging forest disease, sudden oak death. We used this focal system to empirically evaluate whether these relationships were driven by dilution mechanisms that reduce transmission risk for individuals or from the fact that disease was averaged across the host community. We integrated laboratory competence measurements with plant community and symptom data from a large forest monitoring network. Richness increased disease risk for bay laurel trees, dismissing possible dilution mechanisms. Nonetheless, richness was negatively associated with community-level disease prevalence because the disease was aggregated among hosts that vary in disease susceptibility. Aggregating observations (which is surprisingly common in other dilution effect studies) can lead to misinterpretations of dilution mechanisms and bias towards a negative diversity-disease relationship.

Dynamic of Mycobacterium avium subspecies paratuberculosis infection in a domestic-wildlife interface: Domestic sheep and guanaco as reservoir community

Natural herbivore populations have experienced uninterrupted pressures from direct and evident domestic-wildlife interactions and competition, to indirect or less obvious ones such as pathogen transmission. Thus, pathogen spillover between wild and domestic animals is a constant concern because the domestic-wildlife interface represents the ecological frontier in which pathogen transmission takes place in both directions. In Patagonian steppe communities, extensive sheep ranching and guanaco (Lama guanicoe) populations coexist, and guanaco have shown to be infected by pathogens such as Mycobacterium avium subspecies paratuberculosis (MAP) likely transmitted from livestock. MAP causes chronic enteritis and affects mostly domestic ruminants. We evaluated MAP prevalence and pathogen shedding in both species' faeces collected in non-shared and shared sites according to presence/absence of sheep and guanaco along a year, in four different seasons (autumn, winter, and spring 2018, and summer 2019). Our results indicate that MAP circulates in both sheep and guanaco populations with self-sustained transmission; however, both species differ in their levels of competence. We detected higher pathogen shedding in sites occupied by sheep, suggesting that sheep populations may be the main source of infection for susceptible animals due to their large numbers which drive MAP dynamics.

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.