Host-pathogen coevolution increases genetic variation in susceptibility to infection
Abstract
It is common to find considerable genetic variation in susceptibility to infection in natural populations. We have investigated whether natural selection increases this variation by testing whether host populations show more genetic variation in susceptibility to pathogens that they naturally encounter than novel pathogens. In a large cross-infection experiment involving four species of Drosophila and four host-specific viruses, we always found greater genetic variation in susceptibility to viruses that had coevolved with their host. We went on to examine the genetic architecture of resistance in one host species, finding that there are more major-effect genetic variants in coevolved host-pathogen interactions. We conclude that selection by pathogens has increased genetic variation in host susceptibility, and much of this effect is caused by the occurrence of major-effect resistance polymorphisms within populations.
https://doi.org/10.7554/eLife.46440.001Introduction
From bacteria to plants and insects to humans, it is common to find considerable genetic variation in susceptibility to infection in natural populations (Chapman and Hill, 2012; Bangham et al., 2008a; Hammond-Kosack and Jones, 1997; Lazzaro et al., 2004). This variation in susceptibility can determine the impact of disease on health and economic output (Cooke and Hill, 2001; King and Lively, 2012; Alonso-Blanco and Méndez-Vigo, 2014; Burgner et al., 2006). In nature and breeding programs, it determines the ability of populations to evolve resistance to infection. Insect populations, like those of other organisms, typically contain considerable genetic variation in susceptibility to infection (Bangham et al., 2008a; Lazzaro et al., 2004; Tinsley et al., 2006; Obbard and Dudas, 2014), and provide a convenient laboratory model in which to investigate basic questions about how this variation is maintained (Magwire et al., 2012). Within vector species like mosquitoes, resistant genotypes are less likely to transmit pathogens, and this has the potential to reduce disease in vertebrate populations (Beerntsen et al., 2000). Where pathogens are contributing the decline of beneficial species like pollinators, high levels of genetic variation may allow populations to recover (Maori et al., 2007). Understanding the origins of genetic variation in susceptibility is therefore a fundamental question in infectious disease biology.
As pathogens are harmful, natural selection is expected to favour resistant host genotypes. Directional selection on standing genetic variation will drive alleles to fixation, removing variants from the population (Falconer, 1960; Falconer and Mackay, 1996; Blows and Hoffmann, 2005). However, as directional selection also increases the frequency of new mutations that change the trait in the direction of selection, at equilibrium it is expected to have no effect on levels of standing genetic variation (relative to mutation-drift balance; Hill, 1982). However, selection mediated by pathogens may be different. Coevolution with pathogens can result in the maintenance of both resistant and susceptible alleles by negative frequency dependent selection (Woolhouse et al., 2002; Haldane, 1949). Similarly, when infection prevalence exhibits geographical or temporal variation, selection can maintain genetic variation, especially if pleiotropic costs to resistance provide an advantage to susceptible individuals when infection is rare (Nuismer et al., 2003; Thompson, 1999; Koskella, 2018). Even when there is simple directional selection on alleles that increase resistance, the direction of selection by pathogens may frequently change so populations may not be at equilibrium. If selection favours rare alleles – such as new mutations – directional selection can transiently increase genetic variation during their spread through the population (Barton and Turelli, 1987; Bangham et al., 2007; Magwire et al., 2011).
As part of a whole genome association study, we have previously estimated levels of genetic variation in the susceptibility of D. melanogaster to four different viruses (Magwire et al., 2012). We found that there was more genetic variation in susceptibility to the two viruses that were isolated from D. melanogaster than the two viruses from other insect species. Furthermore, in each of these naturally coevolved host-pathogen associations we detected a single major-effect polymorphism affecting resistance. This led us to propose that coevolution had increased genetic variation in susceptibility due to the presence of major-effect resistance polymorphisms. However, this conclusion remains anecdotal. First, aside from two sigma viruses, the viruses we used were mostly very distantly related, so their biology may differ for many reasons. Second, our association study had low statistical power, so conclusions about the genetics were based on just two genes. Finally, and most importantly, the link between coevolutionary history and genetic variation is based on a single host and four viruses, and so could arise by chance. For this reason, here we return to this question and formally test whether a history of coevolution alters the amount and nature of genetic variation.
To examine how selection by a pathogen affects levels of genetic variation we used a natural host-virus system; Drosophila and sigma viruses (Longdon and Jiggins, 2012; Longdon et al., 2012). Sigma viruses are a clade of insect RNA viruses with negative-sense genomes in the family Rhabdoviridae (Longdon et al., 2017; Longdon et al., 2015a; Longdon et al., 2010; Longdon et al., 2011a). They are vertically transmitted through eggs and sperm, and each sigma virus infects a single host species, simplifying studies of coevolution (Longdon et al., 2017; Longdon et al., 2011a). In Drosophila melanogaster, despite the virus causing little adult mortality, infection reduces host fitness by approximately 25% (Wilfert and Jiggins, 2013; Yampolsky et al., 1999). As prevalence in wild populations is typically around 10% (Wilfert and Jiggins, 2013), there is the necessary selective pressure for resistance to evolve (Magwire et al., 2011; Cao et al., 2016). In D. melanogaster, three major-effect resistance alleles have been identified (Bangham et al., 2008a; Magwire et al., 2012; Bangham et al., 2007; Magwire et al., 2011; Cao et al., 2016; Bangham et al., 2008b; Wayne et al., 1996). There has been a recent sweep of D. melanogaster sigma virus (DMelSV) genotypes that are able to overcome one of these host resistance genes (Fleuriet, 1988; Wilfert and Jiggins, 2010a; Wilfert and Jiggins, 2014). Given the power of Drosophila genetics, this system is an excellent model of a coevolutionary arms race between hosts and pathogens.
Sigma viruses offer a novel way to test how coevolution with a pathogen alters the amount of genetic variation in host susceptibility. As sigma viruses are vertically transmitted, we can be certain of which hosts and viruses are naturally coevolving and which are not. In this study we have used four species of Drosophila (D. affinis, D. immigrans, D. melanogaster and D. obscura) that shared a common ancestor approximately 40 million years ago (Obbard et al., 2012; Tamura et al., 2004), and their natural sigma viruses (DAffSV, DImmSV, DMelSV and DObsSV, which have amino acid identities of <55% in the most conserved gene) (Longdon et al., 2015a; Longdon et al., 2010; Longdon et al., 2011b). Despite their vertical mode of transmission, the phylogenies of the viruses and their hosts are incongruent, suggesting they have jumped between host species during their evolution (Longdon et al., 2011b). To test whether selection by viruses increases the amount of genetic variation in host susceptibility we have compared the viral load of endemic viruses that naturally infect each of the four host species to non-endemic viruses. We then examined how selection by these viruses has altered the genetic architecture of resistance by mapping loci that confer resistance to endemic and non-endemic viruses in D. melanogaster.
Results
Genetic variation in susceptibility to infection is greatest in coevolved host-virus associations
To test whether selection by viruses has increased genetic variation in susceptibility to infection, we compared endemic host-virus associations with novel associations that have no history of coevolution. We used four different species of Drosophila, each of which is naturally host to a different sigma virus (Longdon and Jiggins, 2012; Longdon et al., 2012; Longdon et al., 2017; Longdon et al., 2010; Longdon et al., 2011a). We collected four species from the wild and created genetically diverse populations in the laboratory. Using flies from these populations we crossed single males to single females to create full-sib families. The progeny of these crosses were then injected with either the virus isolated from that species, or a virus isolated from one of the other species (Figure 1). Studies on DMelSV have shown that loci that reduce loads when the virus is injected also reduce infection rates in both the lab and field (Bangham et al., 2008a; Longdon et al., 2012; Bregliano, 1970; Brun and Plus, 1980; Ohanessian-Guillemain, 1963; Wilfert and Jiggins, 2010b). As infection is costly, this is expected to increase host fitness. Fifteen days post infection we extracted RNA from the flies and measured viral load by quantitative RT-PCR. The differences between the viral load of different families allowed us to estimate the genetic variance (VG) in viral load – a measure of how much viral resistance varies in the population due to genetic as opposed to environmental causes. In total we infected 52,592 flies and measured the viral load in 4295 biological replicates (a vial containing a mean of 12 flies) across 1436 full-sib families (details of sample sizes in Additional Methods).
Within populations of all four species, we found significantly greater genetic variance in susceptibility to the sigma virus that naturally infects that species compared to viruses from other species (Figure 2; Figure 2—figure supplement 1; Supplementary file 1 and 2). These different variances reflect considerable differences in the mean viral load between families (Figure 2A). For example, when families of D. obscura in the 2nd and 98th percentile were compared, there was a 1294 fold difference between the viral loads of the coevolved virus in the families (Figure 2A). In contrast, for the non-coevolved viruses there was a 27 fold difference in DMelSV loads and a 19 fold difference in DAffSV loads (see Figure 2A for statistics). The data in Figure 2 is zero-centred to allow comparison of the variances, but across the four species, there was no consistent difference between the mean viral load in coevolved versus non-coevolved associations (i.e. the coevolved virus does not always replicate to higher levels suggesting this is not an artefact of simply replicating poorly in a host – see Figure 2—figure supplement 1). Additionally, there was no correlation between the genetic variance in viral load and the mean viral load (Figure 2—figure supplement 1, Spearman’s correlation: = -0.38, S=296, P=0.22).
Major-effect genetic variants that are known to provide resistance to DMelSV do not protect against other viruses
To examine whether the genetic basis of resistance to coevolved and non-coevolved viruses was different we estimated the genetic correlations (rg) in their viral loads. These represent the proportion of genetic variance in viral load between pairs of viruses that shares the same genetic causes. In D. melanogaster these were 0.40 for DMelSV-DAffSV, (95% CIs: 0.20, 0.61) and 0.25 for DMelSV-DObsSV (95% CIs: −0.01,0.47). Similar results were obtained using the D. melanogaster mapping population described below (Supplementary file 3). In the other species our estimates sometimes had wide credible intervals, but the genetic correlations between coevolved and non-coevolved viruses were mostly below 0.5 (Supplementary file 3; note correlations are frequently low between pairs of non-endemic viruses too). Therefore if natural selection increases genetic variation in susceptibility to a natural pathogen, there is expected to be a smaller effect on non-coevolved viruses.
In D. melanogaster a substantial proportion of the genetic variance in susceptibility to DMelSV is explained by major-effect variants in the genes CHKov1 and p62 (also known as Ref(2)P) (Magwire et al., 2012; Bangham et al., 2007; Magwire et al., 2011; Wayne et al., 1996; Contamine et al., 1989). The resistant allele in each of these genes has arisen recently by mutation and been driven up in frequency by natural selection, presumably due to the presence of DMelSV in natural populations (Bangham et al., 2007; Magwire et al., 2011). Therefore, if these genetic variants confer resistance to DMelSV but not the other sigma viruses, then this may explain the differences in genetic variance that we observed.
To examine whether CHKov1 or p62 contributed to the differences in genetic variance we observed in D. melanogaster, we genotyped the parents of the full sib families for variants that confer resistance (Supplementary file 4). Assuming the effects of the resistant alleles are additive, we estimated that the load of the coevolved virus DMelSV was more than halved in homozygous CHKov1 resistant flies compared to susceptible flies (reduction in log2 viral load = 1.2, 95% CI = 0.6, 1.8). In contrast we found no significant effect of this gene on loads of the non-coevolved viruses (DAffSV = −0.2, 95% CI = 0.2,–0.8; DObsSV = −0.4, 95% CI = −0.04, 1.0). The resistant allele of p62 was present at such a low frequency (1.5%) in the population that we lacked statistical power to investigate its effects.
To confirm these results we infected 1869 flies from 32 inbred D. melanogaster lines (Mackay et al., 2012) that had known CHKov1 or p62 genotypes. The effect of these genes was greater on the naturally occurring virus than the viruses from other species (Figure 3; effect of genotype on DMelSV load: F2,28 = 13.2, p=0.00001; DAffSV: F2,29 = 4.9, p=0.01; DObsSV: F2,29 = 5.7, p=0.01).
There are a greater number of major-effect variants in coevolved host-virus associations
To investigate how coevolution shapes the genetics of resistance, we mapped loci controlling resistance using a D. melanogaster advanced intercross population (the DSPR panel [King et al., 2012]). This population samples genetic variation in a small number of genotypes from around the world (the experiments above sampled many genotypes from a single location). It was founded by allowing two sets of 8 inbred founder lines to interbreed for 50 generations, then creating recombinant inbred lines (RILs) whose genomes are a fine-scale mosaic of the original founder genomes. We used 377 RILs from these populations, which have up to 15 alleles of each gene (one founder line is shared between the two populations). We infected 15,916 flies across 1362 biological replicates (a vial containing a mean of 12 flies) with DMelSV, DAffSV or DObsSV and measured viral load as above (see Materials and methods).
We first estimated the genetic variance in viral load within our mapping population. The results recapitulated what we had found above in a natural population of flies — there was considerably more genetic variation in susceptibility to the coevolved virus than the non-coevolved viruses (Figure 4A, filled circles). Therefore, our earlier result from a single population holds when sampling flies from across six continents, although the magnitude of the effect is considerably greater in this mapping population.
To examine the genetic basis of virus resistance, we looked for associations between genotype and viral load across the genome (Figure 4B). In the coevolved association (DMelSV) we identified seven QTL associated with resistance, compared to one that affects DObsSV and none affecting DAffSV (Supplementary file 5; this excludes one DMelSV QTL on the X chromosome that did not remain significant after accounting for the other QTL). The QTL affecting DObsSV also has a significant effect on DMelSV. One of the QTL corresponded to p62 (2L 53 cM). The susceptible allele of CHKov1 was not present in the fly lines assayed.
To examine the effect that the QTL have on viral load, we first split the founder alleles into a resistant class and a susceptible class (see Materials and methods) and then estimated the difference in viral load between the functionally distinct alleles. Six of the seven QTL resulted in greater reductions in the load of the coevolved virus (DMelSV) than the viruses isolated from other species (Figure 4C). There were only two cases where there was substantial cross-resistance to multiple viruses—3R 49 cM confers strong resistance to DMelSV and weak resistance to DObsSV, while 3L 46 cM confers weak resistance to DMelSV and strong resistance to DObsSV.
Together, this modest number of loci with substantial effects on resistance explains most of the high genetic variance in resistance to the coevolved virus (Figure 4A, filled versus open circles). Individually, resistant alleles cause an approximate 3–7 fold reduction in viral load (Figure 4C), and together they explain 59% of the genetic variance in susceptibility to DMelSV, 77% for DObsSV and 3% for DAffSV (Figure 4A, filled versus open circles). However, even after accounting for these genes there remains a significantly higher genetic variance in the viral load of the coevolved virus (Figure 4A, open circles, non-overlapping 95% CI).
Discussion
We have found greater genetic variation in susceptibility to viruses that naturally infect Drosophila compared to viruses that do not, suggesting that selection by these pathogens has acted to increase the amount of genetic variation in susceptibility. This effect was largely caused by a modest number of major-effect genes that explain over half of the genetic variance in resistance.
As the genetic variants in the genes p62 (ref(2)P) and CHKov1 that confer resistance to DMelSV have been identified, this has previously allowed us to use patterns of DNA sequence variation to infer how selection has acted on resistance in D. melanogaster. In both these genes the resistant alleles have arisen relatively recently by mutation and natural selection has pushed them rapidly up in frequency, leaving a characteristic signature of elevated linkage disequilibrium and low genetic diversity around the variant causing resistance (Bangham et al., 2007; Magwire et al., 2011; Wayne et al., 1996). There is no indication of negative frequency dependent selection, and these polymorphisms appear to have arisen from partial selective sweeps (Bangham et al., 2007; Magwire et al., 2011).
The most parsimonious explanation of these observations is that there has been directional selection favouring resistance alleles (although this type of data cannot rule out negative frequency dependent selection as is predicted by models of coevolution). At equilibrium, directional selection on a trait is not expected to affect its genetic variance (relative to a population under mutation-drift balance; Hill, 1982). However, the genetic variance will transiently increase if the variants under selection are initially at low frequency (Barton and Turelli, 1987), as was the case for both p62 and CHKov1 (Bangham et al., 2007; Magwire et al., 2011). A particular feature of pathogens is that the direction of selection is likely to continually change as new pathogens appear in populations or existing pathogens evolve to overcome host defences. For example, in France and Germany in the 1980s, DMelSV evolved to largely overcome the effects of the resistant allele of p62 (Fleuriet and Periquet, 1993; Fleuriet and Sperlich, 1992). Similarly, DImmSV has swept through European populations of D. immigrans in the last ~16 years and DObsSV through UK populations of D. obscura in the last ~11 years (Longdon et al., 2017; Longdon et al., 2011a). If selection by pathogens continually changes and resistance evolves from new mutations, then this may cause a sustained increase in genetic variance in susceptibility to infection.
A key question is whether the increased genetic variation that we see in coevolved Drosophila–sigma virus interactions will hold for coevolved pathogens more generally. Theory suggests that a critical factor determining levels of genetic variation is whether resistance is costly to evolve, as this can result in the maintenance of variation by negative frequency dependent selection (Antonovics and Thrall, 1994; Boots et al., 2014). In humans this has been proposed as an explanation of why there is less genetic variation in susceptibility to pathogens that are effectively controlled by the adaptive immune response, as these resistance mechanisms may be less costly (Baker and Antonovics, 2012). However, it seems unlikely that virus resistance in Drosophila is costly, as experiments have failed to detect costs of DCV resistance (Faria et al., 2015) despite costs of parasitoid and bacterial resistance being repeatedly detected (McGonigle et al., 2017; McKean et al., 2008; Ye et al., 2009). Sigma viruses are also extreme host specialists, so evolutionary changes in resistance will tend to alter pathogen prevalence and so the strength of selection. These epidemiological feedbacks are predicted to frequently increase genetic diversity (Boots et al., 2014; Best et al., 2009; Boots and Haraguchi, 1999). However, in D. melanogaster we see large amounts of genetic variation in susceptibility to viruses that have a broader host range than sigma viruses (DCV and Kallithea virus that infect D. melanogaster and D. simulans in the wild) (Magwire et al., 2012; Palmer et al., 2018; Christian, 1987; Webster et al., 2015). Therefore, it seems unlikely that our conclusions will be a quirk of the sigma virus system. In an analogous study of rust resistance in wild flax plants, sympatric (putatively coevolved) populations had fewer partial resistances than allopatric populations, suggesting more major gene effects, even though overall there was somewhat less genetic variation in susceptibility to sympatric fungal pathogens (Antonovics et al., 2011).
Quantitative traits are typically controlled by a very large number of genetic variants, each of which tends to have a very small effect (Shi et al., 2016; Yang et al., 2010). However, susceptibility to sigma viruses has a simpler genetic basis, with seven polymorphisms explaining over half the genetic variance. This confirms our previous work in D. melanogaster showing a simple genetic basis of virus resistance (Magwire et al., 2012; Cogni et al., 2016). As these genetic variants mostly only affect the naturally occurring pathogen of D. melanogaster, our results suggest that not only is selection by pathogens increasing the genetic variance but it is also altering the genetic architecture of resistance by introducing major-effect variants into the population. One explanation for this observation is that most quantitative traits are under stabilising selection, so major effect variants will tend to be deleterious and removed by selection (Gibson, 2018). In contrast, selection by pathogens likely changes through time and populations may be far from their optimal level of resistance. If this is the case, Fisher’s geometric model predicts that major effect variants will be favoured by directional selection (Fisher, 1930). Alternatively, the coevolution of hosts and parasites can favour discrete susceptible and resistant hosts (Boots et al., 2014), and at the genetic level this may result in major-effect variants (although this theory does not explicitly address the underlying genetics). The simple genetics may also be driven by mutation—for many traits major-effect mutations that increase fitness may be extremely rare, but this may not be the case for virus resistance. For example a single (loss of function) mutation may prevent a virus binding to a host receptor or utilising other parts of the host cellular machinery, and so confer strong resistance.
Regardless of its causes, it may be common that susceptibility to infectious disease has a simple genetic basis. In humans, Hill (2012) advocated the view that susceptibility to infectious disease is qualitatively different from other traits and has a much simpler genetic basis (Hill, 2012). In Drosophila, resistance to DCV and parasitoid wasps both have a simple genetic basis (Magwire et al., 2012; Cao et al., 2017; Orr and Irving, 1997). In plants, major-effect polymorphisms in R genes are commonplace (Hammond-Kosack and Jones, 1997) and in a plant-fungi system genotype-by-genotype interactions explain a larger proportion of the total variance in sympatric (more coevolved) associations (Antonovics et al., 2011). In contrast, studies of bacterial resistance in Drosophila have typically used pathogens that are unlikely to have any history of coevolution, and have found a polygenic basis to resistance (Bou Sleiman et al., 2015; Hotson and Schneider, 2015; Howick and Lazzaro, 2017; Wang et al., 2017). In these studies the polymorphism with the largest effect was found against the only natural D. melanogaster pathogen tested (a polymorphism in Diptericin detected using Providencia rettgeri infection), anecdotally supporting the patterns seen here (Howick and Lazzaro, 2017; Unckless et al., 2015).
A major source of emerging infectious disease is pathogens jumping into novel hosts where they have no co-evolutionary history (Longdon et al., 2014; Parrish et al., 2008). Our results suggest that when a pathogen infects a novel host species, there may be far less genetic variation in susceptibility among individuals than is normally the case. This may create a ‘monoculture effect’ (King and Lively, 2012; Lively, 2010; Ostfeld and Keesing, 2012), which could leave populations vulnerable to epidemics of pathogens that have previously circulated in other host species. Longer term, low levels of pre-standing genetic variation may slow down the rate at which the new host can evolve resistance to a new pathogen.
In conclusion, we have demonstrated that selection by pathogens has increased the amount of genetic variation in host susceptibility. We find resistance has a simple underlying genetic architecture and is largely controlled by major effect resistance loci.
Materials and methods
Virus extraction and infection
Request a detailed protocolWe extracted the sigma viruses DAffSV, DImmSV, DMelSV and DObsSV from infected stocks of D. affinis (line: NC10), D. immigrans (line: DA2), D. melanogaster (line: E320 Ex) and D. obscura (line: 10A) respectively (Longdon et al., 2017; Longdon et al., 2015a; Longdon et al., 2010; Longdon et al., 2011a; Brun and Plus, 1980). These infected lines were collected from the wild between 2007–2012 (all from the UK, bar D.affinis which was collected in the USA). Infected fly stocks were checked for infection by exposing the flies to 100% CO2 at 12°C for 15mins then paralysed flies were collected 30mins later. The DImmSV infected line does not show CO2 sensitivity and so was confirmed to have a high level of infection using RT-PCR. Infected flies were frozen at −80°C and later homogenised in ringers solution (2.5 μl per fly) and centrifuged at 13,000 g for 10 min at 4°C. The supernatant was collected, 2% v/v FBS was added then virus solutions were aliquoted and stored at −80°C.
Experimental design
Request a detailed protocolWe set up a common garden experiment to measure genetic variation in susceptibility to natural and non-natural viruses across four host species. In a fractional factorial experiment each species was infected it with its own virus, as well as two viruses that do not infect that host species (see Figure 1). All of the viruses replicate in all hosts (with the exception of DImmSV in D. melanogaster that showed limited evidence of replication – this combination was not used in the experiment). All fly stocks used (see additional methods for stock details) were tested for existing sigma virus infection using RT-PCR over two generations. For all species we collected flies from the wild and we used a full-sib mating design. The progeny of these crosses were infected by injecting them with 69 nl of the viruses intrathoracically and measuring viral RNA loads 15 days post infection, as in Longdon et al. (2011c). This time point was selected as RNA viral load tends to plateau from around day 15 post infection and there is no mortality from infection in this period. The specifics for each host species, including sample sizes, are detailed in the additional methods.
Known resistance genes in D. melanogaster
Request a detailed protocolWe genotyped parents of each D. melanogaster full sib family from the experiment above for two resistance alleles that are known to confer protection against DMelSV; p62 (Ref(2)P) and CHKov1. We genotyped parental flies using PCR assays that produce different sized products depending on whether flies carry resistant or susceptible alleles. Information on these PCRs and primer sequences can be found in Supplementary file 4. We then calculated the number of resistance alleles in each family by summing the number of alleles from both mothers and fathers. We produced genotype information for 230 of the 255 families.
Another resistance allele has been identified in the gene Ge-1 (Cao et al., 2016). However, this allele has been found to occur at a low frequency in wild populations. We genotyped 184 parental flies from our experiment (parental flies for some families could not be collected) and found the resistant allele was not present, suggesting it is rare or absent.
We further examined the effect of alleles known to affect susceptibility to DMelSV on all three viruses. Firstly, we infected 32 lines from the Drosophila Genetic Reference Panel (DGRP) (Mackay et al., 2012) that were susceptible for both p62 and CHKov1 (n = 11 lines), were resistant for CHKov1 only (n = 13 lines), or were resistant for both genes (n = 8 lines) with DMelSV, DAffSV and DObsSV. No lines in the panel were resistant for p62 and susceptible for CHKov1. We infected a mean of 18 flies per line (range = 3–22).
Mapping resistance genes in D. melanogaster
Request a detailed protocolWe used 377 DSPR lines (154 from panel A and 223 from panel B, http://FlyRILs.org [King et al., 2012; Long et al., 2014]), kindly provided by S.J. Macdonald, University of Kansas) to carry out a Quantitative Trait Locus (QTL) study to examine the genetic basis of resistance to DAffSV, DMelSV and DObsSV in D. melanogaster.
Three females and three males from each DSPR line were placed into yeasted cornmeal vials and allowed to lay for 3–4 days, at 25°C. Male offspring were collected at 0–4 days post-eclosion, and placed at 18°C for 4–6 days. Flies were then injected with DMelSV, DAffSV or DObsSV as described above. Injected males were maintained on unyeasted cornmeal at 18°C, and frozen on day 15 post-infection as above.
Injections were carried out over 13 weeks. Each day of injection a mean of 47 unique lines (range 20–60) and 51 replicate vials were injected with 1–3 different viruses. In total we assayed 377 DSPR lines (108 lines had two biological replicates). Each replicate vial contained a mean of 12 flies (range 1–22). In total, 15,916 flies were injected across both panels of DSPR lines. We injected 319 lines with all three viruses, 38 with 2 viruses and 20 with one virus. The order of injection of lines and of viruses, was randomised across injection days. Independent biological replicates were injected on different days. Panel A and Panel B lines were assayed in two overlapping blocks.
Measuring viral load
Request a detailed protocolWe 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). RNA was extracted from flies homogenised in Trizol and reverse transcribed with 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 2 min followed by 40 cycles of: 95°C for 5 s followed by 60°C for 30 s. Two qRT-PCR reactions (technical replicates) were carried out per sample with both the viral and endogenous control primers. Each qRT-PCR plate contained three standard samples, and all experimental samples were split across plates in a blocked design. A linear model was used to correct the cycle threshold (Ct) values for differences between qRT-PCR plates. Primer sequences are in Supplementary file 6.
To estimate viral load, we calculated ΔCt as the difference between the qRT-PCR cycle thresholds of the virus and the endogenous control. Viral load calculated without using the endogenous control is strongly correlated to ΔCt for all species.
Statistical analysis full-sib experiments
Request a detailed protocolWe used a linear mixed model to examine the amount of genetic variation in susceptibility to the different viruses. We used a trivariate model with the load of the three viruses as the response variable. For each species the model was structured as:
Where is the log2 viral load of the ith biological replicate of full-sib family f infected with virus v. are the fixed effects, with being the mean viral load of each virus. u are the random effects for full-sib families (f) and for the day of injection (d), e are the residuals. By assuming that all the genetic variation in the population is additive (Hill et al., 2008), we estimated the genetic variance (VG) of the viral load as twice the between-family variance (Falconer, 1960; Falconer and Mackay, 1996). Both empirical data and theory suggest additive genetic variation makes up large proportion of the total genetic variance (Hill et al., 2008).
In addition, for D. melanogaster we ran a further model that included the additional fixed effects and that are the linear effects of the CHKov1 and p62 (Ref(2)P) resistance alleles. We assumed these genes had additive effects, and modelled their effects simply as the proportion of resistant alleles in a family (if one parent was heterozygous and the other homozygous susceptible, the value is 0.25).
The model was fitted using the MCMCglmm package in R (Hadfield, 2010). The random effects (and residuals) are assumed to be multivariate normal with zero mean and covariance structure V ⊗ I. I is an identity matrix, and V a matrix of estimated variances and covariances. For the random effects V is a 3 × 3 covariance matrix describing the variances for each virus and the covariances between them. The off-diagonal elements of V for the residual were set to zero because the covariances between traits at these levels are not estimable by design.
Diffuse independent normal priors were placed on the fixed effects (means of zero and variances of 108). Parameter expanded priors were placed on the covariance 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 130 million iterations with a burn-in of 30 million iterations and a thinning interval of 100,000.
We confirmed the results were not sensitive to the choice of prior by also fitting models with inverse-Wishart and flat priors for the variance covariance matrices (described in Longdon et al., 2011c), as well as fitting the models by REML in ASReml R (Gilmour et al., 2002). These analyses all gave qualitatively similar results (data not shown).
Statistical analysis of QTL experiment
Request a detailed protocolBased on genotyping data, the probability that each Recombinant Inbred Line (RIL) in the DSPR panel was derived from each of the eight founder lines has been estimated at 10 kB intervals across the genome (King et al., 2012). To identify QTL affecting viral load we first calculated the mean viral load (ΔCt) across the biological replicates of each RIL. We then regressed the mean viral load against the eight genotype probabilities and calculated logarithm of odds (LOD) scores using the DSPRqtl package in R (King et al., 2012). These LOD scores were calculated separately for DSPR Panel A and Panel B, and then summed at each genomic location. To obtain a significance threshold, we permuted our mean viral load estimates across the RILs within each panel, repeated the analysis above and recorded the highest LOD score across the entire genome. This process was repeated 1000 times to obtain a null distribution of the maximum LOD score.
To estimate the effect of each QTL we assumed that there was a single genetic variant affecting viral load, so the founder alleles could be assigned to two functionally distinct allelic classes. First, we regressed the mean viral load against the genotype probabilities (as described above), resulting in estimates of the mean viral load of each founder allele in the dataset. The two DSPR panels had one founder line (line 8) in common. For QTL where the line eight allele was present in both panels, this analysis included data from both panels, and ‘panel’ was included as a fixed effect in the analysis. When this was not the case, we analysed only data from the panel where the QTL was most significant. We then ranked the founder alleles by viral load estimate, and split this ranked list into all possible groups of two alleles. For each split, the genotype probabilities in the first group of founder alleles were summed. We then regressed mean viral load against each of these combined genotype probabilities. The regression model with the highest likelihood was taken as the most likely classification into allelic classes. The effect size of the QTL was then estimated from this model.
To estimate the genetic variance in viral load within the DSPR panels we modified the model described in Equation 1 as follows. is the log2 viral load of the ith biological replicate of each RIL f infected with virus v. There was a single fixed effect, , of the panel the line is from. u is the random effect for each RIL (f). The day of injection (d) was omitted. As all the RILs are homozygous, we estimated the genetic variance in viral load (VG) as half the between-RIL variance. This assumes all the genetic variation is additive.
To estimate the proportion of the genetic variance that is explained by the QTL we identified, we repeated this analysis but included the 7 QTL we identified as fixed effects in the model. Each QTL was included by estimating the probability that each line carried the resistant allele of the QTL and adding this as a fixed effect to the model. The between-RIL variance then allowed us to estimate the genetic variance in viral load after removing the effects of the QTL.
Additional methods
Request a detailed protocolAll lines were screened for their retrospective sigma virus over two generations by RT-PCR, and infected isofemale lines discarded prior to the experiment.
Drosophila melanogaster
Request a detailed protocolWe created an outcrossed population by combining 150 isofemale lines of D. melanogaster (collected in Accra, Ghana (5.593,–0.188) in 2014) in a population cage. The population was maintained throughout the experiment with a large population size (~1500–2000 flies), with eggs collected from the population cage used to set up each subsequent generation. All rearing was carried out on cornmeal medium (recipe below) sprinkled with live yeast (‘yeasted’) at 25°C.
Virgin flies were collected daily from bottles set up at a controlled egg density. Full-sib families were set up using crosses of single male and female virgins placed in the same vial and aged for 3 days. Each of these families was tipped onto fresh food daily for 5 days to create replicate vials. After 5 days the adult flies were frozen for later genotyping. 12 days after laying, male offspring were collected from each replicate vial and split into two vials of cornmeal medium without any yeast on the surface (‘unyeasted’) and placed at 18°C.
After 5 days these flies were injected with 69 nl of virus extract intra-abdominally using a Nanoject II micro-injector (Drummond scientific). Injected flies were kept at 18°C and tipped onto fresh unyeasted cornmeal every 5 days, before being homogenised in Trizol (Invitrogen) and frozen at −80°C on day 15 post injection for later RNA extraction and qRT-PCR.
Injections were carried out over 25 overlapping blocks. Each block consisted of 10 families, and each day, two vials per family were injected with two different viruses. Each replicate vial contained a mean of 14 flies (range 3–28 flies). In total we measured 255 families over 1567 biological replicates. We aimed to carry out a minimum of 2 replicates of each virus per family, but where possible we carried out 3 or four replicates (92 virus-family combinations had one replicate, 555 had two replicates, 82 had 3 replicates and 25 had four replicates.). 248 families had replicates for all three viruses. Blocks were staggered to overlap with at least 20 families being infected on any one day. The order families were injected in was randomised, and the order the different viruses were injected was blocked across days.
D. immigrans
Request a detailed protocol92 D. immigrans lines were collected from Madingley, Cambridge, UK (52.225, 0.043) in 2012 and 2015. Full-sib families were set up using crosses between the 92 isofemale lines of Drosophila immigrans. Flies were reared on malt food (recipe below) at 18°C. Crosses were between different isofemale lines (i.e. excluding reciprocal crosses) and maximising the number of lines used. Families were established from 2 to 4 day old single female and male virgin flies placed in the same vial for 7 days. These crosses were tipped onto fresh food every 7 days to generate replicate vials of each family. Eclosed males were collected 27–34 days after initial egg laying and injected with DImmSV, DMelSV or DObsSV 1–3 days post-collection, then maintained and frozen on day 15 post-infection as above.
Injections were carried out over 18 overlapping blocks. Each block consisted of an average of 19 families and 46 replicate vials. Each replicate vial contained a mean of 14 flies (range: 4–26). In total we assayed 341 families over 812 biological replicates. We aimed to have a minimum of 2 replicates per virus per family (235 virus-family combinations had one replicate, 270 had two replicates, 11 had 3 replicates and 1 had four replicates). 140 families had replicates across two different viruses and 18 families had replicates for all three viruses. Blocks were staggered to overlap with a mean of 39 families being infected on any one day. The order families were injected in was randomised, and the order the different viruses were injected was blocked across days.
D. affinis
Request a detailed protocolSite | N |
---|---|
Athens, Georgia, USA, (33.946,–83.384) in 2012 | 13 |
Great Smokey Mountain National Park, Gatlinburg, USA (35.698,–83.613) in 2015 | 23 |
Rochester, New York, USA (43.135,–77.599) in 2012 | 4 |
Full-sib families were set up using crosses between 40 isofemale lines of Drosophila affinis (see above for collection details) collected in the U.S. Flies were reared on malt food (recipe below) at 18°C. Crosses were between different isofemale lines (i.e. excluding reciprocal crosses) and maximising the number of lines used. Families were established from 6 day old single female and male virgin flies placed in the same vial for 7 days. These crosses were tipped onto fresh food every 7 days to generate replicate vials of each family. Eclosed males were collected 35–42 days after initial egg laying and then injected with DAffSV, DImmSV or DMelSV 1–3 days post-collection, then maintained and frozen on day 15 post-infection as above.
Injections were carried out over 27 overlapping blocks. Each block consisted of an average of 19 families and 28 replicate vials. Each replicate vial contained a mean of 11 flies (range: 3–23). In total we assayed 520 families over 1003 biological replicates. We aimed to have a minimum of 2 replicates per virus per family (336 virus-family combinations had one replicate, 286 had 2 replicates and 30 had three replicates). 109 families had replicates across two different viruses and 12 families had replicates for all three viruses. Blocks were staggered to overlap with a mean of 23 families being infected on any one day. The order families were injected in was randomised, and the order the different viruses were injected was blocked across days.
D. obscura
Request a detailed protocolSite | N |
---|---|
Derbyshire Site A, UK (52.978,–1.440) in 2012 | 4 |
Derbyshire Site C, UK (52.903,–1.374), in 2012 | 1 |
Les Gorges du Chambon, France (45.622, 0.555) in 2012 | 1 |
Madingley, Cambridge, UK (52.226, 0.046) in 2014 | 15 |
D. obscura were collected in the United Kingdom and France (see above). Males and females were separated, and females were placed in vials to establish isofemale lines. Full-sib families were set up using crosses between 21 isofemale lines of Drosophila obscura collected in the UK. Flies were reared on banana food (recipe below) at 18°C. Crosses were between different isofemale lines (i.e. excluding reciprocal crosses) and maximising the number of lines used. Families were established from 6 day old single female and male virgin flies placed in the same vial for 7 days. These crosses were tipped onto fresh food every 7 days to generate replicate vials of each family. Eclosed males were collected 35–42 days after initial egg laying and then injected with DAffSV, DMelSV, or DObsSV 1–3 days post-collection, then maintained and frozen on day 15 post-infection as above.
Injections were carried out over 25 overlapping blocks. Each block consisted of a mean of 16 families with a mean of 76 vials being injected each day for 12 days. Each replicate vial contained a mean of 8 flies (range: 1–15). In total we assayed 320 families over 913 biological replicates. We aimed to have a minimum of 2 replicates per virus per family (126 virus-family combinations had one replicate, 314 had 2 replicates and 49 had 3 replicates and 3 had four replicates). 94 families had replicates across two different viruses and 39 families had replicates for all three viruses. Blocks were staggered to overlap with at least 10 families being infected on any one day. The order families were injected in was randomised, and the order the different viruses were injected was blocked across days.
Sample size estimation
Request a detailed protocolThe number of full-sib families required for estimating genetic variance in susceptibility was determined by simulation using previous estimates of genetic variation to DMelSV in D. melanogaster (Magwire et al., 2012). After carrying out the full-sib experiment in D. melanogaster, we then down-sampled this data to calculate the minimum number of families required to provide accurate estimates for the other species. Sample sizes for the DSPR experiment were based on previous data (Faria et al., 2015).
Food recipes
Request a detailed protocolBanana:
Mixture 1
1000 ml water
30 g yeast
10 g agar
Mixture 2
20 ml Nipagin
150 g pureed banana
50 g corn syrup
30 g malt powder
Bring mixture one to the boil for 3–4 min, whisk constantly. Add to Mixture two to Mixture 1. Whisk constantly and simmer for 5 min.
Cornmeal:
1200 ml water
13 g agar
105 g dextrose
105 g maize
23 g yeast
Combine and bring to a boil for 5mins, cool to 70°C before adding 35 ml Nipagin (10%)
Malt:
1000 ml water
10 g agar
60 g semolina
20 g yeast
80 g malt extract
Combine and bring to a boil for 5mins, cool to 70°C and then add 14 ml Nipagin (10%) and 5 ml propionic acid.
Data availability
Request a detailed protocolDatasets and R code for estimating the amount of genetic variation in susceptibility https://doi.org/10.6084/m9.figshare.6743339
DGRP dataset https://doi.org/10.6084/m9.figshare.6743354
DSPR dataset and R code https://doi.org/10.6084/m9.figshare.7195751
Data availability
All data generated or analysed during this study are available at Figshare: Datasets and R code for estimating the amount of genetic variation in susceptibility https://doi.org/10.6084/m9.figshare.6743339; DGRP dataset https://doi.org/10.6084/m9.figshare.6743354; DSPR dataset and R code https://doi.org/10.6084/m9.figshare.7195751.
-
FigshareSusceptibility of different Drosophila melanogaster DGRP lines to three sigma viruses.https://doi.org/10.6084/m9.figshare.6743354
-
FigshareNatural selection by pathogens increases genetic variation in susceptibility to infection.https://doi.org/10.6084/m9.figshare.6743339
-
FigshareSusceptibility of Drosophila melanogaster DSPR lines to three sigma viruses.https://doi.org/10.6084/m9.figshare.7195751
References
-
Genetic architecture of naturally occurring quantitative traits in plants: an updated synthesisCurrent Opinion in Plant Biology 18:37–43.https://doi.org/10.1016/j.pbi.2014.01.002
-
Cost of resistance and the maintenance of Genetic-Polymorphism in Host-Pathogen systemsProceedings of the Royal Society B-Biological Sciences 257:105–110.https://doi.org/10.1098/rspb.1994.0101
-
The age and evolution of an antiviral resistance mutation in Drosophila melanogasterProceedings of the Royal Society B: Biological Sciences 274:2027–2034.https://doi.org/10.1098/rspb.2007.0611
-
Genetics of mosquito vector competenceMicrobiology and Molecular Biology Reviews 64:115–137.https://doi.org/10.1128/MMBR.64.1.115-137.2000
-
The implications of coevolutionary dynamics to host-parasite interactionsThe American Naturalist 173:779–791.https://doi.org/10.1086/598494
-
The evolution of costly resistance in Host-Parasite systemsThe American Naturalist 153:359–370.https://doi.org/10.1086/303181
-
Study of infection of germ line in female Drosophila infected with sigma virus. 2. evidence of a correspondance between ovarian cysts with increased virus yield and stabilized progenyAnnales De L Institut Pasteur 119:685–704.
-
BookThe viruses of DrosophilaIn: Ashburner M, Wright T. R. F, editors. The Genetics and Biology of Drosophila. New York: Academic Press. pp. 625–702.
-
Genetic susceptibility to infectious diseases: big is beautiful, but will bigger be even better?The Lancet Infectious Diseases 6:653–663.https://doi.org/10.1016/S1473-3099(06)70601-6
-
Human genetic susceptibility to infectious diseaseNature Reviews Genetics 13:175–188.https://doi.org/10.1038/nrg3114
-
ThesisStudies of Drosophila C and a Viruses in Australian Populations of Drosophila MelanogasterAustralian National University.
-
The genetic architecture of resistance to virus infection in DrosophilaMolecular Ecology 25:5228–5241.https://doi.org/10.1111/mec.13769
-
Genetic resistance to viral infection: the molecular cloning of a Drosophila gene that restricts infection by the rhabdovirus sigmaGenetics 123:525–533.
-
Genetics of susceptibility to human infectious diseaseNature Reviews Genetics 2:967–977.https://doi.org/10.1038/35103577
-
Evolution of the Drosophila melanogaster-sigma virus system in a natural population from tübingenTAG. Theoretical and Applied Genetics. Theoretische Und Angewandte Genetik 85:186–189.https://doi.org/10.1007/BF00222858
-
Population genetics and GWAS: a primerPLOS Biology 16:e2005485.https://doi.org/10.1371/journal.pbio.2005485
-
SoftwareASReml User Guide Release 1.0Vsni.
-
MCMC methods for Multi-Response generalized Linear mixed models: the MCMCglmm R packageJournal of Statistical Software 33:1–22.https://doi.org/10.18637/jss.v033.i02
-
Plant disease resistance genesAnnual Review of Plant Physiology and Plant Molecular Biology 48:575–607.https://doi.org/10.1146/annurev.arplant.48.1.575
-
Evolution, revolution and heresy in the genetics of infectious disease susceptibilityPhilosophical Transactions of the Royal Society B: Biological Sciences 367:840–849.https://doi.org/10.1098/rstb.2011.0275
-
Drosophila melanogaster natural variation affects growth dynamics of infecting listeria monocytogenesG3: Genes|Genomes|Genetics 5:2593–2600.https://doi.org/10.1534/g3.115.022558
-
The effect of host genetic diversity on disease spreadThe American Naturalist 175:E149–E152.https://doi.org/10.1086/652430
-
Sigma viruses from three species of Drosophila form a major new clade in the rhabdovirus phylogenyProceedings of the Royal Society B: Biological Sciences 277:35–44.https://doi.org/10.1098/rspb.2009.1472
-
The evolution and genetics of virus host shiftsPLOS Pathogens 10:e1004395.https://doi.org/10.1371/journal.ppat.1004395
-
Vertically transmitted rhabdoviruses are found across three insect families and have dynamic interactions with their hostsProceedings of the Royal Society B: Biological Sciences 284:20162381.https://doi.org/10.1098/rspb.2016.2381
-
Vertically transmitted viral endosymbionts of insects: do sigma viruses walk alone?Proceedings of the Royal Society B: Biological Sciences 279:3889–3898.https://doi.org/10.1098/rspb.2012.1208
-
The evolutionary costs of immunological maintenance and deploymentBMC Evolutionary Biology 8:76.https://doi.org/10.1186/1471-2148-8-76
-
Coevolution between hosts and parasites with partially overlapping geographic rangesJournal of Evolutionary Biology 16:1337–1345.https://doi.org/10.1046/j.1420-9101.2003.00609.x
-
Estimating divergence dates and substitution rates in the Drosophila phylogenyMolecular Biology and Evolution 29:3459–3473.https://doi.org/10.1093/molbev/mss150
-
The genetics of host-virus coevolution in invertebratesCurrent Opinion in Virology 8:73–78.https://doi.org/10.1016/j.coviro.2014.07.002
-
Etude de facteurs genetiques controlant les relations Du virus sigma et de la drosophile son hoteAnnales De Genetique 5:G1–G64.
-
Effects of host diversity on infectious diseaseAnnual Review of Ecology, Evolution, and Systematics 43:157–182.https://doi.org/10.1146/annurev-ecolsys-102710-145022
-
Cross-species virus transmission and the emergence of new epidemic diseasesMicrobiology and Molecular Biology Reviews 72:457–470.https://doi.org/10.1128/MMBR.00004-08
-
Contrasting the genetic architecture of 30 complex traits from summary association dataThe American Journal of Human Genetics 99:139–153.https://doi.org/10.1016/j.ajhg.2016.05.013
-
Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocksMolecular Biology and Evolution 21:36–44.https://doi.org/10.1093/molbev/msg236
-
Specific hypotheses on the geographic mosaic of coevolutionThe American Naturalist 153:S1–S14.https://doi.org/10.1086/303208
-
Molecular population genetics of ref(2)P, a locus which confers viral resistance in DrosophilaMolecular Biology and Evolution 13:191–199.https://doi.org/10.1093/oxfordjournals.molbev.a025555
-
Host-parasite coevolution: genetic variation in a virus population and the interaction with a host geneJournal of Evolutionary Biology 23:1447–1455.https://doi.org/10.1111/j.1420-9101.2010.02002.x
-
Rapid accumulation of a vertically transmitted parasite triggered by relaxation of natural selection among hostsEvolutionary Ecology Research 1:581–589.
Article and author information
Author details
Funding
NERC Environmental Bioinformatics Centre (NE/L004232/1)
- Elizabeth ML Duxbury
- Jonathan P Day
- Francis M Jiggins
- Ben Longdon
Wellcome (Sir Henry Dale Fellowship (Grant Number 109356/Z/15/Z))
- Ben Longdon
H2020 European Research Council (281668)
- Elizabeth ML Duxbury
Royal Society (Sir Henry Dale Fellowship (Grant Number 109356/Z/15/Z))
- Ben Longdon
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Many thanks to: Alastair Wilson and Jarrod Hadfield for useful advice and discussion; Stuart Macdonald for providing DSPR lines, Trudy Mackay for providing the DGRP fly lines and Kelly Dyer with help collecting D. affinis; Darren Obbard for providing photographs of Drosophila species; Camille Bonneaud, Katherine Roberts, Ryan Imrie and the Unckless lab group for constructive comments on this work. Many thanks to Brian Lazzaro, Janis Antonovics, Bruno Lemaître and one anonymous reviewer for constructive comments that greatly improved the manuscript.
Copyright
© 2019, Duxbury et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 8,462
- views
-
- 757
- downloads
-
- 50
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation by distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60km2 per generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.
-
- Evolutionary Biology
The majority of highly polymorphic genes are related to immune functions and with over 100 alleles within a population, genes of the major histocompatibility complex (MHC) are the most polymorphic loci in vertebrates. How such extraordinary polymorphism arose and is maintained is controversial. One possibility is heterozygote advantage (HA), which can in principle maintain any number of alleles, but biologically explicit models based on this mechanism have so far failed to reliably predict the coexistence of significantly more than ten alleles. We here present an eco-evolutionary model showing that evolution can result in the emergence and maintenance of more than 100 alleles under HA if the following two assumptions are fulfilled: first, pathogens are lethal in the absence of an appropriate immune defence; second, the effect of pathogens depends on host condition, with hosts in poorer condition being affected more strongly. Thus, our results show that HA can be a more potent force in explaining the extraordinary polymorphism found at MHC loci than currently recognized.