The Causes and Consequences of Changes in Virulence following Pathogen Host Shifts
Many emerging infectious diseases are the result of a host shift, with the pathogen jumping into the new host from another species. Virulence—the harm a pathogen does to its host—can be extremely high following a host shift (for example HIV, SARs and Ebola), while other host shifts may go undetected as they cause few symptoms in the new host. We have found that variation in virulence following host shifts can be extremely large and were highly predictable from the host phylogeny, with hosts clustering together in distinct clades displaying high or low virulence. These changes in virulence result from changes in viral load, and therefore the transmission potential of the virus. This suggests there is no clear rule to predict whether a pathogen will be virulent in a novel host. However, it does suggest a simple rule of thumb may be that if a pathogen causes high levels of virulence in any given host species, it will typically cause similar levels of virulence in closely related hosts.
Published in the journal:
The Causes and Consequences of Changes in Virulence following Pathogen Host Shifts. PLoS Pathog 11(3): e32767. doi:10.1371/journal.ppat.1004728
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1004728
Summary
Many emerging infectious diseases are the result of a host shift, with the pathogen jumping into the new host from another species. Virulence—the harm a pathogen does to its host—can be extremely high following a host shift (for example HIV, SARs and Ebola), while other host shifts may go undetected as they cause few symptoms in the new host. We have found that variation in virulence following host shifts can be extremely large and were highly predictable from the host phylogeny, with hosts clustering together in distinct clades displaying high or low virulence. These changes in virulence result from changes in viral load, and therefore the transmission potential of the virus. This suggests there is no clear rule to predict whether a pathogen will be virulent in a novel host. However, it does suggest a simple rule of thumb may be that if a pathogen causes high levels of virulence in any given host species, it will typically cause similar levels of virulence in closely related hosts.
Introduction
Virulence—which we define as the harm a pathogen does to its host—can sometimes dramatically increase when pathogens shift to infect new host species, resulting in new and devastating outbreaks and epidemics [1,2]. For example, in bats, Ebola virus appears to be largely asymptomatic, but it is frequently fatal when it crosses the species barrier into humans and other primates [3,4]. Similarly Henipaviruses appear to be non-pathogenic in pteropoid bats, but can cause high levels of mortality in livestock and humans [5,6].
Despite the importance of these virulence changes, we have little understanding of what determines changes in virulence following a host shift. Historically, it was thought that novel host-parasite associations result in high levels of virulence and that long-term host-parasite interactions lead to the pathogen evolving towards avirulence [7]. However, despite several high profile examples of host shifts resulting in extremely high virulence (see above), there may be a large ascertainment bias in detecting novel infections following hosts shifts, with only the most virulent likely to be detected and many benign infections going unnoticed [8]. Likewise, it has been shown that long term coevolution does not necessarily lead to avirulence [9–12].
To predict how virulence will change following a host shift, we need to understand what causes virulence. Virulence is sometimes thought to be a direct consequence of pathogen replication, with greater levels of replication causing greater amounts of damage to the host [13]. However, virulence and parasite replication can be decoupled [14], especially if hosts mount an inappropriate immune response against novel pathogens [15]. Similarly, virulence can result from infection of a tissue type that has no adaptive value to the pathogen [14,16]. For example, virulence in bacterial meningitis is a result of infection of the central nervous system, and has no relation to transmission potential, or infection of other tissues, where it causes no symptoms [16].
Novel host-parasite interactions have not been under direct selection [17], so the level of pathogen virulence is likely to be maladaptive not only for the host but also the pathogen. This contrasts with natural associations where pathogens have been selected to optimise their virulence in order to maximise their between-host transmission potential [9], and hosts have been selected to reduce the harm caused by their pathogens (to ‘tolerate’ infection) [13,18–22], although one-sided host resistance may evolve in cases where the pathogen-causes a dead end infection [17]. The lack of adaptation in novel host-pathogen associations may explain dead-end chains of transmission seen in some emerging diseases [23]. Here maladaptive levels of virulence may mean the host has less chance of transmitting the pathogen [24], potentially preventing it from becoming established in the population. The classic example of the myxoma virus highlights how virulence may be maladaptive following a host shift. For example, the myxoma virus is relatively benign in South American Sylvilagus rabbits, but it was transferred to European rabbits (Oryctolagus cuniculus) as a biological control agent where it initially had case-mortality rates as high as 99.8% in Australia and drove many populations to near extinction [20,25]. However, following its introduction to European rabbit populations in Australia, case-mortality rates dropped rapidly. This was due to the spread of attenuated virus strains that had higher rates of transmission due to infected hosts surviving, and therefore transmitting the virus, for longer [20,25].
Whilst we know little about how virulence changes following a host shift, the ability of pathogens to infect novel hosts has been more widely studied [26–28]. A number of studies have revealed that the host phylogeny is an important factor in determining the susceptibility of novel hosts [29–35]. The first effect of the host phylogeny is that species closely related to the natural host of a pathogen tend to be more susceptible. This is thought to be because parasites evolve specialised adaptations to their natural host, such as binding to host receptors, avoiding immune responses or utilising host resources, and these break down if the environment provided by the novel host is too different [28]. Secondly, independent from the genetic distance from the pathogen’s natural host, closely-related groups of hosts may have similar levels of susceptibility. Such patterns could arise due to the loss or gain of immune or cellular components in different lineages (e.g. [36,37]), resulting in the host phylogeny being a patchwork of host clades varying in their susceptibility.
We have used a large cross-infection experiment to understand how virulence changes following a host shift. We infected 48 species of Drosophilidae with Drosophila C virus (DCV) and measured virulence, the change in viral load, and the transmission potential of the virus. DCV is a positive sense RNA virus in the family Discistroviridae that was isolated from D. melanogaster and naturally infects D. melanogaster and D. simulans in the wild [38]. The full host range of DCV in the wild is unknown, but it can infect other species of Drosophilidae in the laboratory [39] and can replicate when injected into other dipterans and a moth [40]. DCV causes a reduction in metabolic rate, lowers activity levels, causes intestinal obstruction, lowers the pH of the haemolymph and can ultimately result in death [41–44]. DCV is targeted by antiviral RNAi and a JAK-STAT pathway dependent immune response [45,46]. Here, we have investigated how the host phylogeny determines variation in virulence between species, and the relationships between virulence, viral load and the transmission potential of the virus.
Methods
Virus production
DCV was produced in Schneider’s Drosophila line 2 (DL2) cells [47] as described in [48]. The DCV strain used was isolated from D. melanogaster collected in Charolles, France [49]. Cells were cultured at 26.5°C in Schneider’s Drosophila Medium with 10% Fetal Bovine Serum, 100 U/ml penicillin and 100 μg/ml streptomycin (all Invitrogen, UK). Cells were then freeze-thawed twice to lyse cells and centrifuged at 4000g for 10 minutes at 4°C to remove any cellular components or bacteria. The resulting virus suspension was aliquoted and frozen at -80°C. Uninfected cell culture for control sham inoculations was produced by growing DL2 cells as for virus production but DCV was replaced with Drosophila Ringer’s solution [50]. To calculate infectivity of the virus, serial dilutions of virus from 10-1 to 10-12 were carried out in Schneider’s medium, and each dilution was added to 8 wells of a plate of DL2 cells. After 7 days the wells were examined and classed as “infected” when cell death and cytopathic effects were clearly visible. The Tissue Culture Infective Dose 50 (TCID50) was calculated by the Reed-Muench end-point method [51].
Inoculating flies
48 species of flies which share a common ancestor ~20–50mya were used in this study [52]. All fly stocks were reared at 22°C and 70% relative humidity. Stocks of each fly species were kept in 250ml bottles at staggered ages, and each day freshly eclosed flies were sexed, females were removed, and males were placed on cornmeal medium for 2 days before inoculation. The food medium used for rearing and details of the fly stocks and food recipes used can be found in Table A in S1 Text.
To inoculate flies, a 0.0125 mm diameter stainless steel needle (26002–10, Fine Science Tools, CA, USA) was bent ~0.25 mm from the end, dipped in DCV (TCID50 = 4.64×109) or control solution (uninfected cell culture medium), and the bent part of the needle pricked into the pleural suture on the thorax of flies. Flies were then placed into vials of cornmeal medium and kept at 22°C and 70% relative humidity. These conditions were chosen as they were suitable to maintain all 48 species and carry out infection assays.
As a measure of virulence we recorded mortality after infection. The number of dead flies was counted each day for 20 days and flies were transferred onto fresh medium every 3 days to minimise mortality unrelated to infection. To measure the change in viral load, half of the flies were snap frozen immediately after inoculation in liquid nitrogen as a reference sample to control for relative dose, and the rest were kept for 2 days before being snap frozen in liquid nitrogen. The day 2 time-point was chosen based on time-course data for 10 host species (Figure A in S1 Text), with viral load beginning to plateau after this time. Frozen flies were then homogenised in Trizol reagent (Invitrogen) and stored at -80°C for later RNA extractions.
The mortality inoculations were carried out over a period of 3 days, with the aim of completing a control and virus treatment biological replicate for each fly species each day; i.e. each fly species was included in 2 blocks each day of the experiment. Treatment (virus or control) and the order in which fly species were inoculated was randomized between blocks. The inoculations for measuring the change in viral load were carried out over 6 days, with each species being inoculated each day. Treatment (frozen immediately or on day 2 post infection) and the order the fly species were infected was randomized each day.
In total we measured survival or the change in viral load in 12, 276 flies. In the mortality treatment we measured survival in a mean of 22 flies per replicate (range across species means = 7–25 flies). Out of the 48 species, 41 had 6 biological replicates, 2 had 5 biological replicates, 2 had 4 biological replicates and 3 had 3 biological replicates for the mortality treatments. On average we quantified viral load in a pool of 23 flies per replicate (range across species means = 10–25 flies) with the aim of producing 3 pairs of day 0 and day 2 biological replicates. Out of the 48 species, 45 had 6 biological replicates and 3 had 4 biological replicates for the viral load measurements.
Other factors
Fly stocks were tested for Wolbachia bacterial endosymbionts using PCR primers that amplify the wsp gene [53] prior to the experiment. Wolbachia have been shown to provide resistance to DCV [47,54]) and so only Wolbachia-free species were used; 43 species were naturally Wolbachia free, 5 species were derived from antibiotic treated lines. Species that had a pre-existing DCV infection were also excluded from the experiment.
We also checked that the results were not affected by differences in body size between the species. Wing length is commonly used as a body size measure in Drosophila and strongly correlates with thorax length [55,56]. Wings were removed from ethanol-stored flies collected at the start of the experiment, and photographed under a dissecting microscope. On average we measured wings from 34 individuals per species (range = 13–47). The length of the IV longitudinal vein from the tip of the proximal segment to where the distal segment joins vein V [57] was measured (relative to a standard measurement) using ImageJ software (v1.48) [58].
Measuring the change in viral load
We measured the change in RNA viral load using qRT-PCR. The viral RNA load was expressed relative to the endogenous control housekeeping gene RpL32 (Rp49). We sequenced RpL32 for each species and designed specific RpL32 primers for each species in two conserved regions as described in [30]. We found the DCV primers we used amplified multiple products for D. tropicalis, so a different primer pair was used for this species (as we are measuring the change in relative viral load this should not affect the results). Primer efficiencies were calculated using a dilution series and were close to 100% (DCV: 104%; alternative DCV primers for D. tropicalis: 97%; RpL32: mean = 106%; range of 98–112%). DCV primers were DCV qPCR 599F (5’-GACACTGCCTTTGATTAG-3’) and DCV qPCR 733R (5’-CCCTCTGGGAACTAAATG-3’)[48], the alternative D. tropicalis DCV primers were DCV qPCR 3477F (5’-TTCTTGGTTAGGTCGATTCTTTT-3’) and DCV qPCR 3611R (5’-AATTCTTCGGCTCCAGCTTC-3’).
Total RNA was extracted from Trizol homogenised flies, reverse-transcribed with Promega GoScript reverse transcriptase (Promega) and random hexamer primers, and then diluted 1:10 with nuclease free water. The qRT-PCR was performed on an Applied Biosystems StepOnePlus system using Sensifast Hi-Rox Sybr kit (Bioline) with the following PCR cycle: 95°C for 2min followed by 40 cycles of: 95°C for 5 sec followed by 60°C for 30 sec. Two qRT-PCR reactions (technical replicates) were carried out per sample with both the viral and endogenous control primers. Each qRT-PCR plate contained four standard samples, and all experimental samples were split across plates in a randomised block design. A linear model was used to correct the cycle threshold (Ct) values for differences between qRT-PCR plates. To estimate the change in viral load, we first calculated ΔCt as the difference between the cycle thresholds of the DCV qRT-PCR and the endogenous control. The viral load of day 2 flies relative to day 0 flies was then calculated as 2-ΔΔCt, where ΔΔCt = ΔCtday0 –ΔCtday2, where ΔCtday0 and ΔCtday2 are a pair of ΔCt values from a day 0 biological replicate and a day 2 biological replicate for a particular species. We note that calculating the change in viral load without the use of the endogenous control gene gives equivalent results (Spearman’s correlation between viral load calculated with and without endogenous control = 0.97, P<0.001).
Host phylogeny
The host phylogeny was inferred using the COI, COII, 28S rDNA, Adh, SOD, Amyrel and RpL32 genes. We downloaded all the available sequences from Genbank, and attempted to sequence COI, COII, 28S rDNA, Adh and Amyrel in those species from which they were missing, as described in [30]. This resulted in sequence for all species for COI, COII and RpL32 and partial coverage for the other genes (34 out of 336 species-locus combinations were missing, S1 Table). The sequences of each gene were aligned using ClustalW (alignments and Genbank accession numbers for sequences used are available in an online repository http://dx.doi.org/10.6084/m9.figshare.1112749 and as S1 Dataset). To reconstruct the phylogeny we used BEAST (v1.8.0) [59] as this allows construction of an ultrametric (time-based) tree using a relaxed molecular clock model. Using ultrametric trees assumes that rates of evolutionary change (virulence and viral load) are proportional to time (rather than proportional to the rate of sequence evolution in the genes used to construct the phylogeny). The genes were partitioned into 3 groups each with their own substitution and molecular clock models. The three partitions were: mitochondrial (COI, COII); ribosomal (28S); and nuclear (Adh, SOD, Amyrel, RpL32). Each of the partitions used an HKY substitution model [60] (which allows transitions and transversions to occur at different rates) with a gamma distribution of rate variation with 4 categories and estimated base frequencies. Additionally the mitochondrial and nuclear data sets were partitioned into codon positions 1+2 and 3, with unlinked substitution rates and base frequencies across codon positions. Empirical studies suggest that HKY models with codon partitions are a good fit for most protein-coding data sets [61]. A random starting tree was used, with a relaxed uncorrelated lognormal molecular clock and we used no external temporal information, so all dates are relative to the root age. The tree-shape prior was set to a speciation-extinction (birth-death) process. The BEAST analysis was run twice for 100 million MCMC generations sampled every 1000 steps. The two MCMC processes were examined using the program Tracer (v1.6) [62] to ensure convergence and adequate sampling. The two runs were combined using the programme LogCombiner (v1.8). A maximum clade credibility tree was created following a 10% burnin using TreeAnnotator (v1.8) and visualised using FigTree (v. 1.3.1) [63] (Tree nexus file http://dx.doi.org/10.6084/m9.figshare.1112749).
Statistical analysis
We used a phylogenetic mixed model to examine the effects of host species relatedness [64–66] using the MCMCglmm R package ([67], R Foundation for Statistical Computing, Vienna, Austria). We used a trivariate formulation with mortality of the controls, mortality of the virus infected flies and viral load as response variables. Viral load was treated as Gaussian, whereas the mortality data were treated as longitudinal binomial data as in event history analysis [68]. The daily binomial counts are the number of flies that died and survived between successive days, with days where both counts are zero (i.e. all flies in the vial died the previous day) omitted as they are uninformative. The mortality model is equivalent to a proportional odds survival analysis (with censoring) [69]. Mortality and viral load data are provided in S2 Dataset.
The model structure is similar to the trivariate model outlined in [30]:
where ηthij is the jth linear predictor for the ith biological replicate of host species h for trait t. βv For the viral load trait, j always equals 1, but for the mortality traits j takes on values 1 to n, where n is either 20 days (the length of time over which mortality was recorded) or the number of days it took for all flies in a vial to die if this occurred in less than 20 days.β are fixed effects with β1 being the intercepts for each trait and β2 and β3 the linear effects of wing size and genetic distance from D. melanogaster respectively. For the two mortality traits, separate linear and quadratic age effects were also included for each treatment, as was a linear effect of days since the flies were transferred onto fresh food medium: δ is one if the trait is a mortality trait and zero otherwise
u are random effects for vial (v—mortality traits only), phylogenetic species effects (p) and non-phylogenetic species effects (s) and e are residuals. The random effects (and residuals) are assumed to be multivariate normal with zero mean and covariance structure V ⊗ A for phylogenetic effects or V ⊗ I otherwise. I is an identity matrix, A the phylogenetic relatedness matrix and V a matrix of estimated variances and covariances. For the phylogenetic and non-phylogenetic species effects V is a 3x3 covariance matrix describing the (phylogenetic and non-phylogenetic) inter-specific variances in each trait and the inter-specific covariances between them. For the residuals V is also a 3x3 matrix, but for the vial effects that are only defined for the two mortality traits, V is a 2x2 covariance matrix. The off-diagonal elements of V for the residual and vial effects were set to zero because the covariances between traits at these levels are not estimable by design.
The linear predictor defines the expected value of the response E[ythij]/wthij = g-1 (ηthij) where g-1 is the inverse-link function, and wthij is a weight. For the viral load trait g-1 is the identity function and the weights are one. For the mortality traits g-1 is the inverse logit function (exp(ηthij))/(1+ exp(ηthij)) and the weights are the number of flies alive at the end of the previous observation (or for j = 1 the number of flies initially present).
Diffuse independent normal priors were placed on the fixed effects (means of zero and variances of 108). Parameter expanded priors were placed on the (co)variance matrices resulting in scaled multivariate F distributions which have the property that the marginal distributions for the variances are scaled (by 1000) F1,1. The exceptions were the residual variances for which an inverse-gamma prior was used with shape and scale equal to 0.001. The MCMC chain was ran for 1300000 iterations with a burn-in of 300000 and a thinning interval of 500.
We confirmed the results were not sensitive to the choice of prior by also fitting models with inverse-Wishart and flat priors for the matrices describing the phylogenetic and non-phylogenetic species effects (described in [30]), which gave qualitatively similar results (data not shown). We also ran models with only phylogenetic or non-phylogenetic species effects, which gave equivalent results in all instances (data not shown). To confirm our results were robust to different analysis methods we also analysed the mortality data by calculating the proportion of flies dead on each day of the experiment, taking the mean across all the days (which will be equivalent to the area under the survival curve), and arcsine square-root transforming this number. We then used this together with the change in viral load as responses in a trivariate model as described above. This gave very similar results to the proportional odds model (data not shown). Finally, we also ran models where the logit-linear dependency varies across species (variable age-dependent mortality), and these models also gave very similar results (data not shown).
Transmission of DCV
In order to understand how DCV is transmitted, and how viral load might relate to transmission, we performed two additional transmission experiments. In the first experiment a single infected adult fly of the focal species was housed with adult sentinel flies (from a D. melanogaster isogenetic line) and the infection status of the sentinels (adult-adult transmission) and their offspring (adult+parent-offspring transmission) assayed. In the second experiment corpses of single adult flies (of the focal species) that died through DCV infection were housed with sentinel eggs, and the infection status of the emergent adult sentinels assayed (dead adult-embryo/larvae transmission). In both cases, we examined the transmission rates and tested for a correlation between transmission success and viral load from the previous experiment, which was measured 2 days after inoculation. To ensure that this viral load was correlated with the viral load in corpses of flies that died from a DCV infection, we also measured viral load in corpses for a subset of species. Full methods (S1 Text) and data (S1 and S3 Datasets) can be found in the Supporting Information.
Results
Host shifts result in large changes in virulence
To investigate the effect of switching hosts on virulence, we inoculated 2967 flies from 48 species of Drosophilidae with DCV and recorded mortality (virulence) for 20 days. We found dramatic differences in the virulence of DCV when it infects different host species (Fig. 1). In some species it appears to be a largely benign infection, with the virus causing no increase in mortality over 20 days. At the other extreme, ~90% of D. persimilis and D. pseudoobscura flies are killed by day 4 post-infection.
To confirm that these differences reflect the virulence of the virus rather than intrinsic differences in the survivorship of the different species, alongside the virus infections we also inoculated 2993 flies from the 48 species with a control solution. There was far less mortality in the controls than the virus infected flies (range of mortality for controls: 0–18% and infected: 0–89%, see Figs. 1 and 2). Nonetheless, there was significant inter-specific variation (the sum of both phylogenetic and non-phylogenetic variation) in control mortality, which accounted for 43% of the between-vial variation in mortality (95% CI = 15%-74%), but this was not correlated with mortality in infected flies (inter-specific correlation = -0.08, CI = -0.53, 0.34). We also found no significant relationship between wing size and the mortality of infected flies (-0.44, 95% CI = -2.56, 1.78).
Virulence following host shifts is determined by the host phylogeny
To examine how the host phylogeny influences changes in virulence, we reconstructed the phylogeny of our 48 Drosophilidae species using seven genes. The resulting tree is broadly consistent with previous studies of these taxa [70], with the close phylogenetic relationships being generally well supported (Fig. 2).
Because DCV was isolated from and naturally infects D. melanogaster [38], we examined whether virulence was affected by genetic distance from this host species. We found there was no significant change in virulence with distance from D. melanogaster (Fig. 2; slope = -0.76, 95% CI = -3.37, 1.62).
Despite there being no effect of genetic distance from D. melanogaster on virulence, there is a striking pattern of high- and low-virulence hosts clustering together on the phylogeny (Fig. 2). We used a phylogenetic mixed model to partition the inter-specific variance in mortality into that which is explained by a Brownian motion model of evolution on the host phylogeny (vp), and a species-specific component independent of the phylogeny (vs). We then calculated the proportion of between-species variance in the treated group that is explained by the host phylogeny (vp/(vp + vs), which is related to Pagel’s lambda or phylogenetic heritability [65,66,71]. 75% of the between-species variation in virulence is explained by the host phylogeny (95% CI = 48%-97%).
The effect of the host phylogeny on virulence can be clearly seen by using this model to reconstruct mortality rates across the phylogeny on different days post infection (Fig. 3 and video http://dx.doi.org/10.6084/m9.figshare.1192890). Early in the course of infection significant mortality only occurs in a single small clade of species (obscura group), but as time progresses other clades start to experience high mortality, and by day 16 most of the survivors belong to a single clade of species (virilis group) where virulence is relatively low (Fig. 3).
Virulence following host shifts is strongly correlated to viral load
To measure viral load we inoculated 6316 flies from all 48 species and measured the change in viral RNA load by qRT-PCR (Figure B in S1 Text). Similar to virulence, there are extremely large differences in viral load between species, with the total inter-specific variance explaining 91% of the variance in viral load. Indeed, we found roughly a billion times the amount of viral RNA in the most susceptible host species compared to the most resistant species (Fig. 2). This was unrelated to the size of the flies, with no significant effect of wing size on viral load (slope = 0.15, 95% CI = -6.04, 6.27),
Species with high viral loads tend to cluster together on the phylogeny (Fig. 4). This was confirmed by fitting a phylogenetic mixed model, where 67% (95% CI = 33%-93%) of the between-species variation in viral load is explained by the host phylogeny. As was the case for virulence, we found there was no significant change in viral load with distance from D. melanogaster, the species from which the virus was isolated (slope = -2.01, 95% CI = -7.90, 7.08).
Variation in virulence could either be a consequence of different pathogen loads, or due to similar pathogen loads causing different levels of harm to their hosts. To separate these hypotheses, we examined the relationship between virulence and viral load. We found strong positive phylogenetic and non-phylogenetic correlations between virulence and viral load (Fig. 4, phylogenetic correlation = 0.76, 95% CI = 0.44, 0.99; non-phylogenetic correlation = 0.47, 95% CI = 0.05, 0.83). There is no correlation between viral load and the mortality that occurred in the mock-infected control flies (inter-specific correlation = 0.06, 95% CI = -0.38, 0.52).
Virulence and transmission rate in novel hosts
We investigated the efficiency of different routes of transmission by allowing uninfected sentinel D. melanogaster to come into contact with infected flies (see supplementary information). We then tested pools of sentinel flies for the presence of DCV. In our assay we could not detect any adult to larval transmission in the 8 species we tested, despite these species all having high viral loads (n = 21 biological replicates). This suggests shedding of the virus is not a major source of infection for larvae. Next, we measured transmission between adult flies, and detected transmission in 11% of vials (n = 107 biological replicates across 38 species). Finally, we tested whether transmission occurred from the corpses of dead flies to embryos or larvae, and found that this was a highly efficient transmission route, with 80% of vials of sentinel flies testing positive for DCV infection (n = 44 biological replicates across 15 species).
Because transmission from dead flies occurs at a high rate, the viral load in dead flies may determine the transmission rate of the virus. Therefore, we checked that our measurements of viral load on day 2 post-infection reflected viral loads on the day of death. We infected 15 species of fly, collected their corpses on the day that they died, and measured their viral load. We found that the viral load on the day of death is tightly correlated to the day 2 viral load (Figure C in S1 Text, inter-specific correlation: 0.84, 95% CI = 0.62, 0.99).
Discussion
We have observed dramatic variation in virulence when we infected different host species with a viral pathogen. We found both increases and decreases in virulence compared to D. melanogaster, the host from which the virus was isolated. These changes in virulence were explained by the host phylogeny, with very similar levels of virulence shown by closely related hosts (Fig. 2). Virulence was tightly coupled to viral load (Fig. 4), suggesting the amount of harm caused to a host is a result of virus accumulation in a host, and there is little interspecific variation in tolerance. This in turn means that changes in virulence following host shifts may be accompanied by changes in the transmission rate of the virus.
Large shifts in virulence following pathogen host shifts can result in devastating epidemics, such as myxomatosis in rabbits [19] or Ebola shifting from bats to humans [3]. These changes in virulence often appear random and unpredictable consequences of the pathogen being poorly adapted to its new host [8] or the host being poorly adapted to the pathogen [17]. Our results show that the changes in virulence can be explained by the host phylogeny, with the virus causing similar levels of mortality in closely related species. However, those clades of hosts where virulence is high are scattered across the phylogeny (Fig. 3), so host shifts across large genetic distances may sometimes result in virulent infections and other times benign associations.
DCV was isolated from D. melanogaster where its virulence is at an intermediate level, and transfer to different species can both increase and decrease virulence. Despite seeing dramatic changes in virulence, there is no tendency for it to decline with genetic distance from the natural host. This contrasts with studies that have investigated infection success in novel hosts and found that it tends to decline in species less closely related to the natural host [29–35]. However, it is difficult to compare our results to these studies as there has been no systematic survey of the host range of DCV, and it may infect species outside of the melanogaster subgroup in nature.
The patterns that we see will depend on the taxonomic scale we are looking at—outside of insects it is unlikely that DCV can infect any hosts [40]. The species we looked at shared a common ancestor 40 million years ago. Combined with the rapid generation time of insects, this makes our results relevant to the taxonomic distances over which hosts shifts normally occur. For example, mammals shared a common ancestor 95 million years ago [72], but have far slower generation times and lower rates of molecular evolution [73,74]
The physiological and molecular reasons why DCV virulence is so variable across Drosophila species are a matter of speculation. A virus’s ability to successfully infect a novel host will depend on its ability to bind to host cells, replicate in suitable tissues and avoid or suppress the host immune response. The ability of viruses to successfully bind to host receptors is important for a successful host shift. For example, changes to the viral capsid of parvoviruses that allow the binding of host cells are essential for them to successfully host shift from cats to dogs [75]. In D. melanogaster RNAi is known to be an important antiviral immune response [46]; to counter this host immune response, DCV has evolved a suppressor of RNAi that binds dsRNA and inhibits the production of siRNAs [76]. Therefore, the patterns observed may be due to species varying in their ability to produce an antiviral RNAi response, or due to the viral suppressor of RNAi being more or less efficient in different host species, as has been observed for another Drosophila virus [77]. In D. melanogaster, certain alleles of a restriction factor called pastrel confer DCV resistance [78], and similar factors may be responsible for the differences in virulence observed between species. The patterns observed could also be due to certain host clades being assayed at sub-optimal environmental conditions (e.g. temperature) that affect immunity or other physiological processes. Whilst we used hosts without Wolbachia endosymbionts [47,54], other components of the host microbiota or unknown viral infections may also affect resistance. Regardless of the cause, these changes have happened in the common ancestor of certain clades, resulting in all the descendant species in different clades having similar viral loads and levels of virulence, resulting in the strong phylogenetic patterns we observe (Fig. 3).
Evolutionary theories of virulence commonly assume a direct link between virulence and pathogen loads, but this assumption has been criticised because virulence is often an indirect effect of the pathogen, such as an inappropriate immune response ([8,79], reviewed in [7]). However, in this system we found virulence to be very tightly correlated with viral load. This suggests virulence is a direct consequence of the damage caused by high parasite burdens, although it could also result from viral virulence factors resulting in the high viral loads.
Because virulence is associated with high viral loads, changes in virulence following host shifts are likely to be linked to changes in the transmission rate of the virus. Indeed, it is difficult to conceive that the million-fold range in viral loads that we observed would not affect transmission, and our limited data on transmission tentatively supports this link (see supplementary information). Therefore, host shifts may result in maladaptive viral loads and low rates of transmission in the new host.
When virulence is coupled to pathogen load it can sometimes create a trade-off—higher viral loads result in a higher rate of transmission per unit time but the host dies faster and transmits for less time—meaning that intermediate levels of virulence may be optimal for the pathogen to maximise its fitness [9]. In our case we found high rates of transmission of DCV from the corpses of dead flies to larvae, suggesting killing the host could actually aid transmission. If the majority of transmission occurs from the corpses of dead flies in nature there will be a strong mechanistic coupling between transmission and virulence [80], suggesting the most successful strategy for DCV to evolve will be to maximise its replication rate at any cost to the host. However, the true optimal virulence for DCV remains a matter of speculation, as in nature it will be affected by host population structure [81], environmental conditions [82] and ecological interactions [83]. Regardless, it is clear that host shifts may result in DCV having maladaptive levels of virulence that result in low rates of onward transmission, and this may prevent the pathogen from becoming established in that species.
In conclusion, we have found that host shifts can result in both large increases and decreases in virulence. The fact that virulence is associated with high pathogen burdens, which in turn may lead to higher levels of transmission, suggests that avirulent host shifts may be unsuccessful. We found virulence was largely explained by the host phylogeny, with clades of closely related hosts displaying similar levels of virulence. While this study suggests there is no clear rule to predict whether a pathogen will be virulent in a novel host, it does suggest that if a pathogen causes high levels of virulence in any given host species, it will typically cause similar levels of virulence in closely related hosts.
Supporting Information
Zdroje
1. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, et al. (2008) Global trends in emerging infectious diseases. Nature 451: 990–993. doi: 10.1038/nature06536 18288193
2. Woolhouse ME, Haydon DT, Antia R (2005) Emerging pathogens: the epidemiology and evolution of species jumps. Trends Ecol Evol 20: 238–244. 16701375
3. Leroy EM, Kumulungui B, Pourrut X, Rouquet P, Hassanin A, et al. (2005) Fruit bats as reservoirs of Ebola virus. Nature 438: 575–576. 16319873
4. Swanepoel R, Leman PA, Burt FJ, Zachariades NA, Braack LEO, et al. (1996) Experimental inoculation of plants and animals with Ebola virus. Emerging Infectious Diseases 2: 321–325. 8969248
5. Chua KB, Bellini WJ, Rota PA, Harcourt BH, Tamin A, et al. (2000) Nipah virus: a recently emergent deadly paramyxovirus. Science 288: 1432–1435. 10827955
6. Middleton DJ, Morrissy CJ, van der Heide BM, Russell GM, Braun MA, et al. (2007) Experimental Nipah virus infection in pteropid bats (Pteropus poliocephalus). J Comp Pathol 136: 266–272. 17498518
7. Alizon S, Hurford A, Mideo N, Van Baalen M (2009) Virulence evolution and the trade-off hypothesis: history, current state of affairs and the future. J Evol Biol 22: 245–259. doi: 10.1111/j.1420-9101.2008.01658.x 19196383
8. Weiss RA (2002) Virulence and pathogenesis. Trends Microbiol 10: 314–317. 12110209
9. Anderson RM, May RM (1982) Coevolution of hosts and parasites. Parasitology 85 (Pt 2): 411–426. 6755367
10. Anderson RM, May RM (1991) Infectious Diseases of humans; dynamics and control. Oxford: Oxford University Press.
11. Read AF (1994) The evolution of virulence. Trends Microbiol 2: 73–76. 8156274
12. Toft CA, Karter AJ (1990) Parasite-host coevolution. Trends Ecol Evol 5: 326–329. doi: 10.1016/0169-5347(90)90179-H 21232384
13. Ebert D, Bull JJ (2008) The evolution and expression of virulence. In: Stearns SC, Koella JC, editors. Evolution in Health and Disease. 2nd ed. Oxford: Oxford University Press. pp. 153–167.
14. Margolis E, Levin BR (2008) The evolution of bacteria-host interactions: virulence and the immune over-response. In: Gutirrez JQ, Baquero F, editors. Introduction to the evolutionary biology of bacterial and fungal pathogens. Washington D.C.: ASM Press.
15. Graham AL, Allen JE, Read AF (2005) Evolutionary causes and consequences of immunopathology. Annual Review of Ecology Evolution and Systematics 36: 373–397.
16. Levin BR, Bull JJ (1994) Short-sighted evolution and the virulence of pathogenic microorganisms. Trends Microbiol 2: 76–81. 8156275
17. Antonovics J, Boots M, Ebert D, Koskella B, Poss M, et al. (2013) The origin of specificity by means of natural selection: evolved and nonhost resistance in host-pathogen interactions. Evolution 67: 1–9. doi: 10.1111/j.1558-5646.2012.01793.x 23289557
18. Jensen KH, Little TJ, Skorping A, Ebert D (2006) Empirical support for optimal virulence in a castrating parasite. PLoS Biol 4: e197. doi: 10.1371/annotation/6f580f9f-d724-433c-9e12-f402fac28829 16719563
19. Marshall ID, Fenner F (1958) Studies in the epidemiology of infectious myxomatosis of rabbits. V. Changes in the innate resistance of Australian wild rabbits exposed to myxomatosis. J Hyg (Lond) 56: 288–302. 13563871
20. Kerr PJ (2012) Myxomatosis in Australia and Europe: a model for emerging infectious diseases. Antiviral Res 93: 387–415. doi: 10.1016/j.antiviral.2012.01.009 22333483
21. de Roode JC, Yates AJ, Altizer S (2008) Virulence-transmission trade-offs and population divergence in virulence in a naturally occurring butterfly parasite. Proc Natl Acad Sci U S A 105: 7489–7494. doi: 10.1073/pnas.0710909105 18492806
22. Fraser C, Hollingsworth TD, Chapman R, de Wolf F, Hanage WP (2007) Variation in HIV-1 set-point viral load: Epidemiological analysis and an evolutionary hypothesis. Proceedings of the National Academy of Sciences of the United States of America 104: 17441–17446. 17954909
23. Blumberg S, Lloyd-Smith JO (2013) Inference of R(0) and transmission heterogeneity from the size distribution of stuttering chains. PLoS Comput Biol 9: e1002993. doi: 10.1371/journal.pcbi.1002993 23658504
24. Andre JB, Hochberg ME (2005) Virulence evolution in emerging infectious diseases. Evolution 59: 1406–1412. 16153027
25. Fenner F, Ratcliffe FN (1965) Myxomatosis. Cambridge: Cambridge University press.
26. Engelstadter J, Hurst GD (2006) The dynamics of parasite incidence across host species. Evolutionary Ecology 20: 603–616.
27. Waxman D, Weinert LA, Welch JJ (2014) Inferring host range dynamics from comparative data: the protozoan parasites of new world monkeys. Am Nat 184: 65–74. doi: 10.1086/676589 24921601
28. Longdon B, Brockhurst MA, Russell CA, Welch JJ, Jiggins FM (2014) The Evolution and Genetics of Virus Host Shifts. PLoS Pathog 10: e1004395. doi: 10.1371/journal.ppat.1004395 25375777
29. Faria NR, Suchard MA, Rambaut A, Streicker DG, Lemey P (2013) Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraints. Philos Trans R Soc Lond B Biol Sci 368: 20120196. doi: 10.1098/rstb.2012.0196 23382420
30. Longdon B, Hadfield JD, Webster CL, Obbard DJ, Jiggins FM (2011) Host phylogeny determines viral persistence and replication in novel hosts. PLoS Pathogens 7: e1002260. doi: 10.1371/journal.ppat.1002260 21966271
31. Streicker DG, Turmelle AS, Vonhof MJ, Kuzmin IV, McCracken GF, et al. (2010) Host Phylogeny Constrains Cross-Species Emergence and Establishment of Rabies Virus in Bats. Science 329: 676–679. doi: 10.1126/science.1188836 20689015
32. de Vienne DM, Hood ME, Giraud T (2009) Phylogenetic determinants of potential host shifts in fungal pathogens. Journal of Evolutionary Biology 22: 2532–2541. doi: 10.1111/j.1420-9101.2009.01878.x 19878406
33. Gilbert GS, Webb CO (2007) Phylogenetic signal in plant pathogen-host range. Proceedings of the National Academy of Sciences of the United States of America 104: 4979–4983. 17360396
34. Perlman SJ, Jaenike J (2003) Infection success in novel hosts: An experimental and phylogenetic study of Drosophila-parasitic nematodes. Evolution 57: 544–557. 12703944
35. Tinsley MC, Majerus MEN (2007) Small steps or giant leaps for male-killers? Phylogenetic constraints to male-killer host shifts. Bmc Evolutionary Biology 7.
36. Jiggins FM, Kim KW (2005) The evolution of antifungal peptides in Drosophila. Genetics 171: 1847–1859. 16157672
37. Salazar-Jaramillo L, Paspati A, van de Zande L, Vermeulen CJ, Schwander T, et al. (2014) Evolution of a cellular immune response in Drosophila: a phenotypic and genomic comparative analysis. Genome Biology and Evolution.
38. Christian PD (1987) Studies of Drosophila C and A viruses in Australian populations of Drosophila melanogaster: Australian National University. 305 p.
39. Kapun M, Nolte V, Flatt T, Schlotterer C (2010) Host Range and Specificity of the Drosophila C Virus. Plos One 5: e12421. doi: 10.1371/journal.pone.0012421 20865043
40. Jousset FX (1976) Host Range of Drosophila-Melanogaster C Virus among Diptera and Lepidoptera. Annales De Microbiologie A127: 529–&. 823856
41. Chtarbanova S (2011) Tissue-specific pathologies induced by two RNA viruses in Drosophila melanogaster: University of Strasbourg.
42. Arnold PA, Johnson KN, White CR (2013) Physiological and metabolic consequences of viral infection in Drosophila melanogaster. J Exp Biol 216: 3350–3357. doi: 10.1242/jeb.088138 23685974
43. Chtarbanova S, Lamiable O, Lee KZ, Galiana D, Troxler L, et al. (2014) Drosophila C virus systemic infection leads to intestinal obstruction. J Virol 88: 14057–14069. doi: 10.1128/JVI.02320-14 25253354
44. Ferreira AG, Naylor H, Esteves SS, Pais IS, Martins NE, et al. (2014) The toll-dorsal pathway is required for resistance to viral oral infection in Drosophila. PLoS Pathog 10: e1004507. doi: 10.1371/journal.ppat.1004507 25473839
45. Dostert C, Jouanguy E, Irving P, Troxler L, Galiana-Arnoux D, et al. (2005) The Jak-STAT signaling pathway is required but not sufficient for the antiviral response of drosophila. Nat Immunol 6: 946–953. 16086017
46. Kemp C, Mueller S, Goto A, Barbier V, Paro S, et al. (2013) Broad RNA interference-mediated antiviral immunity and virus-specific inducible responses in Drosophila. J Immunol 190: 650–658. doi: 10.4049/jimmunol.1102486 23255357
47. Teixeira L, Ferreira A, Ashburner M (2008) The Bacterial Symbiont Wolbachia Induces Resistance to RNA Viral Infections in Drosophila melanogaster. Plos Biology 6: 2753–2763.
48. Longdon B, Cao C, Martinez J, Jiggins FM (2013) Previous Exposure to an RNA Virus Does Not Protect against Subsequent Infection in Drosophila melanogaster. Plos One 8: e73833. doi: 10.1371/journal.pone.0073833 24040086
49. Jousset FX, Plus N, Croizier G, Thomas M (1972) [Existence in Drosophila of 2 groups of picornavirus with different biological and serological properties]. C R Acad Sci Hebd Seances Acad Sci D 275: 3043–3046. 4631976
50. Sullivan W, Ashburner M, Hawley S (2000) Drosophila Protocols. New York: Cold Spring Harbor Laboratory Press.
51. Reed LJ, Muench H (1938) A simple method of estimating fifty per cent endpoints. The American Journal of Hygiene 27: 493–497.
52. Obbard DJ, Maclennan J, Kim K-W, Rambaut A, O’Grady PM, et al. (2012) Estimating divergence dates and substitution rates in the Drosophila phylogeny. Molecular Biology and Evolution 29: 3459–3473. doi: 10.1093/molbev/mss150 22683811
53. Zhou W, Rousset F, O'Neil S (1998) Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences. Proc Biol Sci 265: 509–515. 9569669
54. Hedges LM, Brownlie JC, O'Neill SL, Johnson KN (2008) Wolbachia and Virus Protection in Insects. Science 322: 702–702. doi: 10.1126/science.1162418 18974344
55. Huey RB, Moreteau B, Moreteau JC, Gibert P, Gilchrist GW, et al. (2006) Sexual size dimorphism in a Drosophila clade, the D-obscura group. Zoology 109: 318–330. 16978850
56. Sokoloff A (1966) Morphological Variation in Natural and Experimental Populations of Drosophila Pseudoobscura and Drosophila Persimilis. Evolution 20: 49–71.
57. Gilchrist GW, Huey RB, Serra L (2001) Rapid evolution of wing size clines in Drosophila subobscura. Genetica 112–113: 273–286.
58. Rasband WS (1997–2011) ImageJ, U. S. National Institutes of Health, Bethesda,Maryland, USA. Available: http://imagej.nih.gov/ij/. v1.43u ed.
59. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. Bmc Evolutionary Biology 7: 214. 17996036
60. Hasegawa M, Kishino H, Yano TA (1985) DATING OF THE HUMAN APE SPLITTING BY A MOLECULAR CLOCK OF MITOCHONDRIAL-DNA. Journal of Molecular Evolution 22: 160–174. 3934395
61. Shapiro B, Rambaut A, Drummond AJ (2006) Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol Biol Evol 23: 7–9. 16177232
62. Rambaut A, Drummond AJ (2007) Tracer v14, Available from http://beast.bio.ed.ac.uk/Tracer
63. Rambaut A (2011) FigTree. v1.3 ed.
64. Hadfield JD, Nakagawa S (2010) General quantitative genetic methods for comparative biology: phylogenies, taxonomies and multi-trait models for continuous and categorical characters. Journal of Evolutionary Biology 23: 494–508. doi: 10.1111/j.1420-9101.2009.01915.x 20070460
65. Housworth EA, Martins EP, Lynch M (2004) The phylogenetic mixed model. American Naturalist 163: 84–96. 14767838
66. Lynch M (1991) METHODS FOR THE ANALYSIS OF COMPARATIVE DATA IN EVOLUTIONARY BIOLOGY. Evolution 45: 1065–1080.
67. Hadfield JD (2010) MCMC Methods for Multi-Response Generalized Linear Mixed Models: The MCMCglmm R Package. Journal of Statistical Software 33: 1–22. 20808728
68. Aalen OO, Borgan Ø, Gjessing HK (2008) Survival and Event History Analysis: A Process Point of View. New York: Springer.
69. Bennett S (1983) Log-Logistic Regression-Models for Survival-Data. Applied Statistics-Journal of the Royal Statistical Society Series C 32: 165–171.
70. van der Linde K, Houle D, Spicer GS, Steppan SJ (2010) A supermatrix-based molecular phylogeny of the family Drosophilidae. Genetics Research 92: 25–38. doi: 10.1017/S001667231000008X 20433773
71. Pagel M (1999) Inferring the historical patterns of biological evolution. Nature 401: 877–884. 10553904
72. Murphy WJ, Pringle TH, Crider TA, Springer MS, Miller W (2007) Using genomic data to unravel the root of the placental mammal phylogeny. Genome Research 17: 413–421. 17322288
73. Welch J, Bininda-Emonds O, Bromham L (2008) Correlates of substitution rate variation in mammalian protein-coding sequences. Bmc Evolutionary Biology 8: 53. doi: 10.1186/1471-2148-8-53 18284663
74. Thomas JA, Welch JJ, Lanfear R, Bromham L (2010) A Generation Time Effect on the Rate of Molecular Evolution in Invertebrates. Molecular Biology and Evolution 27: 1173–1180. doi: 10.1093/molbev/msq009 20083649
75. Truyen U, Evermann JF, Vieler E, Parrish CR (1996) Evolution of canine parvovirus involved loss and gain of feline host range. Virology 215: 186–189. 8560765
76. van Rij RP, Saleh MC, Berry B, Foo C, Houk A, et al. (2006) The RNA silencing endonuclease Argonaute 2 mediates specific antiviral immunity in Drosophila melanogaster. Genes & Development 20: 2985–2995.
77. van Mierlo JT, Overheul GJ, Obadia B, van Cleef KW, Webster CL, et al. (2014) Novel Drosophila viruses encode host-specific suppressors of RNAi. PLoS Pathog 10: e1004256. doi: 10.1371/journal.ppat.1004256 25032815
78. Magwire MM, Fabian DK, Schweyen H, Cao C, Longdon B, et al. (2012) Genome-wide association studies reveal a simple genetic basis of resistance to naturally coevolving viruses in Drosophila melanogaster. Plos Genetics 8: e1003057. doi: 10.1371/journal.pgen.1003057 23166512
79. Ebert D, Bull JJ (2003) Challenging the trade-off model for the evolution of virulence: is virulence management feasible? Trends in Microbiology 11: 15–20. 12526850
80. Frank SA, Schmid-Hempel P (2008) Mechanisms of pathogenesis and the evolution of parasite virulence. J Evol Biol 21: 396–404. doi: 10.1111/j.1420-9101.2007.01480.x 18179516
81. Boots M, Sasaki A (1999) ‘Small worlds’ and the evolution of virulence: infection occurs locally and at a distance. 1933–1938 p.
82. Vale PF, Choisy M, Little TJ (2013) Host nutrition alters the variance in parasite transmission potential. Biol Lett 9: 20121145. doi: 10.1098/rsbl.2012.1145 23407498
83. Martinez J, Longdon B, Bauer S, Chan YS, Miller WJ, et al. (2014) Symbionts commonly provide broad spectrum resistance to viruses in insects: a comparative analysis of Wolbachia strains. PLoS Pathog 10: e1004369. doi: 10.1371/journal.ppat.1004369 25233341
84. Russo CA, Takezaki N, Nei M (1995) Molecular phylogeny and divergence times of drosophilid species. Mol Biol Evol 12: 391–404. 7739381
Štítky
Hygiena a epidemiológia Infekčné lekárstvo LaboratóriumČlánok vyšiel v časopise
PLOS Pathogens
2015 Číslo 3
- Parazitičtí červi v terapii Crohnovy choroby a dalších zánětlivých autoimunitních onemocnění
- Očkování proti virové hemoragické horečce Ebola experimentální vakcínou rVSVDG-ZEBOV-GP
- Koronavirus hýbe světem: Víte jak se chránit a jak postupovat v případě podezření?
Najčítanejšie v tomto čísle
- Bacterial Immune Evasion through Manipulation of Host Inhibitory Immune Signaling
- Antimicrobial-Induced DNA Damage and Genomic Instability in Microbial Pathogens
- Is Antigenic Sin Always “Original?” Re-examining the Evidence Regarding Circulation of a Human H1 Influenza Virus Immediately Prior to the 1918 Spanish Flu
- An 18 kDa Scaffold Protein Is Critical for Biofilm Formation