Dissecting the genetic architecture of a stepwise infection process

Acclimation to warmer temperatures can protect host populations from both further heat stress and the potential invasion of pathogens

Thermal acclimation can provide an essential buffer against heat stress for host populations, while acting simultaneously on various life-history traits that determine population growth. In turn, the ability of a pathogen to invade a host population is intimately linked to these changes via the supply of new susceptible hosts, as well as the impact of warming on its immediate infection dynamics. Acclimation therefore has consequences for hosts and pathogens that extend beyond simply coping with heat stress-governing both population growth trajectories and, as a result, an inherent propensity for a disease outbreak to occur. The impact of thermal acclimation on heat tolerances, however, is rarely considered simultaneously with metrics of both host and pathogen population growth, and ultimately fitness. Using the host Daphnia magna and its bacterial pathogen, we investigated how thermal acclimation impacts host and pathogen performance at both the individual and population scales. We first tested the effect of maternal and direct thermal acclimation on the life-history traits of infected and uninfected individuals, such as heat tolerance, fecundity, and lifespan, as well as pathogen infection success and spore production. We then predicted the effects of each acclimation treatment on rates of host and pathogen population increase by deriving a host's intrinsic growth rate (rm) and a pathogen's basic reproductive number (R0). We found that direct acclimation to warming enhanced a host's heat tolerance and rate of population growth, despite a decline in life-history traits such as lifetime fecundity and lifespan. In contrast, pathogen performance was consistently worse under warming, with within-host pathogen success, and ultimately the potential for disease spread, severely hampered at higher temperatures. Our results suggest that hosts could benefit more from warming than their pathogens, but only by linking multiple individual traits to population processes can the full impact of higher temperatures on host and pathogen population dynamics be realised.

Is there a sicker sex? Dose relationships modify male-female differences in infection prevalence

Throughout the animal kingdom, there are striking differences in the propensity of one sex or the other to become infected. However, precisely when we should expect males or females to be the sicker sex remains unclear. A major barrier to answering this question is that very few studies have considered how the susceptibility of males and females changes across the full range of pathogen doses encountered in nature. Without quantifying this 'dose-susceptibility' relationship, we have likely underestimated the scope for sex differences to arise. Here, we use the Daphnia magnia-Pasteuria ramosa system to reveal that sex differences in susceptibility are entirely dose-dependent, with pathogens having a higher probability of successfully establishing an infection in mature males at low doses, but mature females at high doses. The scope for male-female differences to emerge is therefore much greater than previously appreciated-extending to sex differences in the upper limits to infection success, per-propagule infectivity risks and density-dependent pathogen behaviour. Applying this expanded scope across the animal kingdom will help us understand when and why a sicker sex emerges, and the implications for diseases in nature-where sex ratios, age structure and pathogen densities vary drastically.

Pharmaceutical pollution alters the cost of bacterial infection and its relationship to pathogen load

The relationship between pathogen proliferation and the cost of infection experienced by a host drives the ecology and evolution of host-pathogen dynamics. While environmental factors can shape this relationship, there is currently limited knowledge on the consequences of emerging contaminants, such as pharmaceutical pollutants, on the relationship between a pathogen's growth within the host and the damage it causes, termed its virulence. Here, we investigated how exposure to fluoxetine (Prozac), a commonly detected psychoactive pollutant, could alter this key relationship using the water flea Daphnia magna and its bacterial pathogen Pasteuria ramosa as a model system. Across a variety of fluoxetine concentrations, we found that fluoxetine shaped the damage a pathogen caused, such as the reduction in fecundity or intrinsic growth experienced by infected individuals, but with minimal change in average pathogen spore loads. Instead, fluoxetine modified the relationship between the degree of pathogen proliferation and its virulence, with both the strength of this trade-off and the component of host fitness most affected varying by fluoxetine concentration and host genotype. Our study underscores the potential for pharmaceutical pollution to modify the virulence of an invading pathogen, as well as the fundamental trade-off between host and pathogen fitness, even at the trace amounts increasingly found in natural waterways.

Can pharmaceutical pollution alter the spread of infectious disease? A case study using fluoxetine

Human activity is changing global environments at an unprecedented rate, imposing new ecological and evolutionary ramifications on wildlife dynamics, including host-parasite interactions. Here we investigate how an emerging concern of modern human activity, pharmaceutical pollution, influences the spread of disease in a population, using the water flea Daphnia magna and the bacterial pathogen Pasteuria ramosa as a model system. We found that exposure to different concentrations of fluoxetine-a widely prescribed psychoactive drug and widespread contaminant of aquatic ecosystems-affected the severity of disease experienced by an individual in a non-monotonic manner. The direction and magnitude of any effect, however, varied with both the infection outcome measured and the genotype of the pathogen. By contrast, the characteristics of unexposed animals, and thus the growth and density of susceptible hosts, were robust to fluoxetine. Using our data to parameterize an epidemiological model, we show that fluoxetine is unlikely to lead to a net increase or decrease in the likelihood of an infectious disease outbreak, as measured by a pathogen's transmission rate or basic reproductive number. Instead, any given pathogen genotype may experience a twofold change in likely fitness, but often in opposing directions. Our study demonstrates that changes in pharmaceutical pollution give rise to complex genotype-by-environment interactions in its influence of disease dynamics, with repercussions on pathogen genetic diversity and evolution. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

QTL study reveals candidate genes underlying host resistance in a Red Queen model system

Specific interactions of host and parasite genotypes can lead to balancing selection, maintaining genetic diversity within populations. In order to understand the drivers of such specific coevolution, it is necessary to identify the molecular underpinnings of these genotypic interactions. Here, we investigate the genetic basis of resistance in the crustacean host, Daphnia magna, to attachment and subsequent infection by the bacterial parasite, Pasteuria ramosa. We discover a single locus with Mendelian segregation (3:1 ratio) with resistance being dominant, which we call the F locus. We use QTL analysis and fine mapping to localize the F locus to a 28.8-kb region in the host genome, adjacent to a known resistance supergene. We compare the 28.8-kb region in the two QTL parents to identify differences between host genotypes that are resistant versus susceptible to attachment and infection by the parasite. We identify 13 genes in the region, from which we highlight eight biological candidates for the F locus, based on presence/absence polymorphisms and differential gene expression. The top candidates include a fucosyltransferase gene that is only present in one of the two QTL parents, as well as several Cladoceran-specific genes belonging to a large family that is represented in multiple locations of the host genome. Fucosyltransferases have been linked to resistance in previous studies of Daphnia-Pasteuria and other host-parasite systems, suggesting that P. ramosa spore attachment could be mediated by changes in glycan structures on D. magna cuticle proteins. The Cladoceran-specific candidate genes suggest a resistance strategy that relies on gene duplication. Our results add a new locus to a growing genetic model of resistance in the D. magna-P. ramosa system. The identified candidate genes will be used in future functional genetic studies, with the ultimate aim to test for cycles of allele frequencies in natural populations.

Genetic variation and microbiota in bumble bees cross-infected by different strains of C. bombi

The bumblebee Bombus terrestris is commonly infected by a trypanosomatid gut parasite Crithidia bombi. This system shows a striking degree of genetic specificity where host genotypes are susceptible to different genotypes of parasite. To a degree, variation in host gene expression underlies these differences, however, the effects of standing genetic variation has not yet been explored. Here we report on an extensive experiment where workers of twenty colonies of B. terrestris were each infected by one of twenty strains of C. bombi. To elucidate the host's genetic bases of susceptibility to infection (measured as infection intensity), we used a low-coverage (~2 x) genome-wide association study (GWAS), based on angsd, and a standard high-coverage (~15x) GWAS (with a reduced set from a 8 x 8 interaction matrix, selected from the full set of twenty). The results from the low-coverage approach remained ambiguous. The high-coverage approach suggested potentially relevant genetic variation in cell surface and adhesion processes. In particular, mucin, a surface mucoglycoprotein, potentially affecting parasite binding to the host gut epithelia, emerged as a candidate. Sequencing the gut microbial community of the same bees showed that the abundance of bacterial taxa, such as Gilliamella, Snodgrassella, or Lactobacillus, differed between 'susceptible' and 'resistant' microbiota, in line with earlier studies. Our study suggests that the constitutive microbiota and binding processes at the cell surface are candidates to affect infection intensity after the first response (captured by gene expression) has run its course. We also note that a low-coverage approach may not be powerful enough to analyse such complex traits. Furthermore, testing large interactions matrices (as with the full 20 x 20 combinations) for the effect of interaction terms on infection intensity seems to blur the specific host x parasite interaction effects, likely because the outcome of an infection is a highly non-linear process dominated by variation in individually different pathways of host defence (immune) responses.

Population-level variation in parasite resistance due to differences in immune initiation and rate of response

Closely related populations often differ in resistance to a given parasite, as measured by infection success or failure. Yet, the immunological mechanisms of these evolved differences are rarely specified. Does resistance evolve via changes to the host's ability to recognize that an infection exists, actuate an effective immune response, or attenuate that response? We tested whether each of these phases of the host response contributed to threespine sticklebacks' recently evolved resistance to their tapeworm Schistocephalus solidus. Although marine stickleback and some susceptible lake fish permit fast-growing tapeworms, other lake populations are resistant and suppress tapeworm growth via a fibrosis response. We subjected lab-raised fish from three populations (susceptible marine "ancestors," a susceptible lake population, and a resistant lake population) to a novel immune challenge using an injection of (1) a saline control, (2) alum, a generalized pro-inflammatory adjuvant that causes fibrosis, (3) a tapeworm protein extract, or (4) a combination of alum and tapeworm protein. With enough time, all three populations generated a robust fibrosis response to the alum treatments. Yet, only the resistant population exhibited a fibrosis response to the tapeworm protein alone. Thus, these populations differed in their ability to respond to the tapeworm protein but shared an intact fibrosis pathway. The resistant population also initiated fibrosis faster in response to alum, and was able to attenuate fibrosis, unlike the susceptible populations' slow but longer lasting response to alum. As fibrosis has pathological side effects that reduce fecundity, the faster recovery by the resistant population may reflect an adaptation to mitigate the costs of immunity. Broadly, our results confirm that parasite detection and immune initiation, activation speed, and immune attenuation simultaneously contribute to the evolution of parasite resistance and adaptations to infection in natural populations.

Infection phenotypes of a coevolving parasite are highly diverse, structured, and specific

Understanding how diversity is maintained in natural populations is a major goal of evolutionary biology. In coevolving hosts and parasites, negative frequency-dependent selection is one mechanism predicted to maintain genetic variation. While much is known about host diversity, parasite diversity remains understudied in coevolutionary research. Here, we survey natural diversity in a bacterial parasite by characterizing infection phenotypes for over 50 isolates in relation to 12 genotypes of their host, Daphnia magna. We find striking phenotypic variation among parasite isolates, and we discover the parasite can infect its host through at least five different attachment sites. Variation in attachment success at each site is explained to varying degrees by host and parasite genotypes. A spatial correlation analysis showed that infectivity of different isolates does not correlate with geographic distance, meaning isolates from widespread populations are equally able to infect the host. Overall, our results reveal that infection phenotypes of this parasite are highly diverse. Our results are consistent with the prediction that under Red Queen coevolutionary dynamics both the host and the parasite should show high genetic diversity for traits of functional importance in their interactions.

Host Controls of Within-Host Disease Dynamics: Insight from an Invertebrate System

AbstractWithin-host processes (representing the entry, establishment, growth, and development of a parasite inside its host) may play a key role in parasite transmission but remain challenging to observe and quantify. We develop a general model for measuring host defenses and within-host disease dynamics. Our stochastic model breaks the infection process down into the stages of parasite exposure, entry, and establishment and provides associated probabilities for a host's ability to resist infections with barriers and clear internal infections. We tested our model on Daphnia dentifera and the parasitic fungus Metschnikowia bicuspidata and found that when faced with identical levels of parasite exposure, Daphnia patent (transmitting) infections depended on the strength of internal clearance. Applying a Gillespie algorithm to the model-estimated probabilities allowed us to visualize within-host dynamics, within which signatures of host defense could be clearly observed. We also found that early within-host stages were the most vulnerable to internal clearance, suggesting that hosts have a limited window during which recovery can occur. Our study demonstrates how pairing longitudinal infection data with a simple model can reveal new insight into within-host dynamics and mechanisms of host defense. Our model and methodological approach may be a powerful tool for exploring these properties in understudied host-parasite interactions.

A Two-Locus System with Strong Epistasis Underlies Rapid Parasite-Mediated Evolution of Host Resistance

Parasites are a major evolutionary force, driving adaptive responses in host populations. Although the link between phenotypic response to parasite-mediated natural selection and the underlying genetic architecture often remains obscure, this link is crucial for understanding the evolution of resistance and predicting associated allele frequency changes in the population. To close this gap, we monitored the response to selection during epidemics of a virulent bacterial pathogen, Pasteuria ramosa, in a natural host population of Daphnia magna. Across two epidemics, we observed a strong increase in the proportion of resistant phenotypes as the epidemics progressed. Field and laboratory experiments confirmed that this increase in resistance was caused by selection from the local parasite. Using a genome-wide association study, we built a genetic model in which two genomic regions with dominance and epistasis control resistance polymorphism in the host. We verified this model by selfing host genotypes with different resistance phenotypes and scoring their F1 for segregation of resistance and associated genetic markers. Such epistatic effects with strong fitness consequences in host–parasite coevolution are believed to be crucial in the Red Queen model for the evolution of genetic recombination.

Among-Strain Variation in Resistance of Paramecium caudatum to the Endonuclear Parasite Holospora undulata: Geographic and Lineage-Specific Patterns

Resistance is a key determinant in interactions between hosts and their parasites. Understanding the amount and distribution of variation in this trait between strains can provide insights into (co)evolutionary processes and their potential to shape patterns of diversity in natural populations. Using controlled inoculation in experimental mass cultures, we investigated the quantitative variation in resistance to the bacterial parasite Holospora undulata across a worldwide collection of strains of its ciliate host Paramecium caudatum. We combined the observed variation with available information on the phylogeny and biogeography of the strains. We found substantial variation in resistance among strains, with upper-bound values of broad-sense heritability >0.5 (intraclass correlation coefficients). Strain estimates of resistance were repeatable between laboratories and ranged from total resistance to near-complete susceptibility. Early (1 week post inoculation) measurements provided higher estimates of resistance heritability than did later measurements (2-3 weeks), possibly due to diverging epidemiological dynamics in replicate cultures of the same strains. Genetic distance (based on a neutral marker) was positively correlated with the difference in resistance phenotype between strains (r = 0.45), essentially reflecting differences between highly divergent clades (haplogroups) within the host species. Haplogroup A strains, mostly European, were less resistant to the parasite (49% infection prevalence) than non-European haplogroup B strains (28%). At a smaller geographical scale (within Europe), strains that are geographically closer to the parasite origin (Southern Germany) were more susceptible to infection than those from further away. These patterns are consistent with a picture of local parasite adaptation. Our study demonstrates ample natural variation in resistance on which selection can act and hints at symbiont adaptation producing signatures in geographic and lineage-specific patterns of resistance in this model system.

Pathogen exposure reduces sexual dimorphism in a host's upper thermal limits

The climate is warming at an unprecedented rate, pushing many species toward and beyond the upper temperatures at which they can survive. Global change is also leading to dramatic shifts in the distribution of pathogens. As a result, upper thermal limits and susceptibility to infection should be key determinants of whether populations continue to persist, or instead go extinct. Within a population, however, individuals vary in both their resistance to both heat stress and infection, and their contributions to vital growth rates. No more so is this true than for males and females. Each sex often varies in their response to pathogen exposure, thermal tolerances, and particularly their influence on population growth, owing to the higher parental investment that females typically make in their offspring. To date, the interplay between host sex, infection, and upper thermal limits has been neglected. Here, we explore the response of male and female Daphnia to bacterial infection and static heat stress. We find that female Daphnia, when uninfected, are much more resistant to static heat stress than males, but that infection negates any advantage that females are afforded. We discuss how the capacity of a population to cope with multiple stressors may be underestimated unless both sexes are considered simultaneously.

Neutral genomic signatures of host-parasite coevolution

Background

Coevolution is a selective process of reciprocal adaptation in hosts and parasites or in mutualistic symbionts. Classic population genetics theory predicts the signatures of selection at the interacting loci of both species, but not the neutral genome-wide polymorphism patterns. To bridge this gap, we build an eco-evolutionary model, where neutral genomic changes over time are driven by a single selected locus in hosts and parasites via a simple biallelic gene-for-gene or matching-allele interaction. This coevolutionary process may lead to cyclic changes in the sizes of the interacting populations.

Results

We investigate if and when these changes can be observed in the site frequency spectrum of neutral polymorphisms from host and parasite full genome data. We show that changes of the host population size are too smooth to be observable in its polymorphism pattern over the course of time. Conversely, the parasite population may undergo a series of strong bottlenecks occurring on a slower relative time scale, which may lead to observable changes in a time series sample. We also extend our results to cases with 1) several parasites per host accelerating relative time, and 2) multiple parasite generations per host generation slowing down rescaled time.

Conclusions

Our results show that time series sampling of host and parasite populations with full genome data are crucial to understand if and how coevolution occurs. This model provides therefore a framework to interpret and draw inference from genome-wide polymorphism data of interacting species.

Independent effects on cellular and humoral immune responses underlie genotype-by-genotype interactions between Drosophila and parasitoids

It is common to find abundant genetic variation in host resistance and parasite infectivity within populations, with the outcome of infection frequently depending on genotype-specific interactions. Underlying these effects are complex immune defenses that are under the control of both host and parasite genes. We have found extensive variation in Drosophila melanogaster's immune response against the parasitoid wasp Leptopilina boulardi. Some aspects of the immune response, such as phenoloxidase activity, are predominantly affected by the host genotype. Some, such as upregulation of the complement-like protein Tep1, are controlled by the parasite genotype. Others, like the differentiation of immune cells called lamellocytes, depend on the specific combination of host and parasite genotypes. These observations illustrate how the outcome of infection depends on independent genetic effects on different aspects of host immunity. As parasite-killing results from the concerted action of different components of the immune response, these observations provide a physiological mechanism to generate phenomena like epistasis and genotype-interactions that underlie models of coevolution.