Genetic structure and host–parasite co‐divergence: evidence for trait‐specific local adaptation

Host–parasite co‐evolution can lead to genetic differentiation among isolated host–parasite populations and local adaptation between parasites and their hosts. However, tests of local adaptation rarely consider multiple fitness‐related traits although focus on a single component of fitness can be misleading. Here, we concomitantly examined genetic structure and co‐divergence patterns of the trematode Coitocaecum parvum and its crustacean host Paracalliope fluviatilis among isolated populations using the mitochondrial cytochrome oxidase I gene (COI). We then performed experimental cross‐infections between two genetically divergent host–parasite populations. Both hosts and parasites displayed genetic differentiation among populations, although genetic structure was less pronounced in the parasite. Data also supported a co‐divergence scenario between C. parvum and P. fluviatilis potentially related to local co‐adaptation. Results from cross‐infections indicated that some parasite lineages seemed to be locally adapted to their sympatric (home) hosts in which they achieved higher infection and survival rates than in allopatric (away) amphipods. However, local, intrinsic host and parasite characteristics (host behavioural or immunological resistance to infections, parasite infectivity or growth rate) also influenced patterns of host–parasite interactions. For example, overall host vulnerability to C. parvum varied between populations, regardless of parasite origin (local vs. foreign), potentially swamping apparent local co‐adaptation effects. Furthermore, local adaptation effects seemed trait specific; different components of parasite fitness (infection and survival rates, growth) responded differently to cross‐infections. Overall, data show that genetic differentiation is not inevitably coupled with local adaptation, and that the latter must be interpreted with caution in a multi‐trait context.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Linkage disequilibrium clustering‐based approach for association mapping with tightly linked genomewide data

Genomewide association studies (GWAS) aim to identify genetic markers strongly associated with quantitative traits by utilizing linkage disequilibrium (LD) between candidate genes and markers. However, because of LD between nearby genetic markers, the standard GWAS approaches typically detect a number of correlated SNPs covering long genomic regions, making corrections for multiple testing overly conservative. Additionally, the high dimensionality of modern GWAS data poses considerable challenges for GWAS procedures such as permutation tests, which are computationally intensive. We propose a cluster‐based GWAS approach that first divides the genome into many large nonoverlapping windows and uses linkage disequilibrium network analysis in combination with principal component (PC) analysis as dimensional reduction tools to summarize the SNP data to independent PCs within clusters of loci connected by high LD. We then introduce single‐ and multilocus models that can efficiently conduct the association tests on such high‐dimensional data. The methods can be adapted to different model structures and used to analyse samples collected from the wild or from biparental F₂ populations, which are commonly used in ecological genetics mapping studies. We demonstrate the performance of our approaches with two publicly available data sets from a plant (Arabidopsis thaliana) and a fish (Pungitius pungitius), as well as with simulated data.     

Genetic variation within and between populations of hermaphroditic <Emphasis Type="Italic">Bulinus truncatus</Emphasis> tetraploid freshwater snails of the Albertine Rift,East Africa

Genetic variability within and among Bulinus truncatus of the Albertine Rift freshwater bodies were assessed to investigate the degree of inbreeding and gene flow in the snail populations. The effect of ploidy on the genetic structuring of B. truncatus is also described. We characterized the genetic structure of seven B. truncatus populations from Lake Albert, Lake Kivu, and Katosho swamp in Tanzania using five polymorphic microsatellite loci. Genetic differentiation was quantified using pairwise FST values and Nei’s standard genetic distances. Different alleles were observed across all loci and genetic diversity was low although it varied greatly across populations; observed heterozygosity was, however, higher than the expected heterozygosity in three of the populations studied. Significant heterozygote deficiencies were observed coupled with significant linkage disequilibria in five populations for all the five loci examined in this study. We found significant genetic differentiation among the seven freshwater bodies; private alleles were observed across all loci indicating restricted or absence of gene flow between populations. Limited snail dispersal and the reproductive biology of B. truncatus are the major forces shaping the genetic variation observed. Low genetic variation within B. truncatus populations exposes them to a high parasite infection risk as predicted in the Red Queen hypothesis.     

Simulated climate change,epidemic size,and host evolution across host–parasite populations

Climate change is causing warmer and more variable temperatures as well as physical flux in natural populations, which will affect the ecology and evolution of infectious disease epidemics. Using replicate seminatural populations of a coevolving freshwater invertebrate‐parasite system (host: Daphnia magna, parasite: Pasteuria ramosa), we quantified the effects of ambient temperature and population mixing (physical flux within populations) on epidemic size and population health. Each population was seeded with an identical suite of host genotypes and dose of parasite transmission spores. Biologically reasonable increases in environmental temperature caused larger epidemics, and population mixing reduced overall epidemic size. Mixing also had a detrimental effect on host populations independent of disease. Epidemics drove parasite‐mediated selection, leading to a loss of host genetic diversity, and mixed populations experienced greater evolution due to genetic drift over the season. These findings further our understanding of how diversity loss will reduce the host populations’ capacity to respond to changes in selection, therefore stymying adaptation to further environmental change.     

Experimental evolution of local parasite maladaptation

Sign and magnitude of local adaptation in host–parasite systems may vary with ecological, epidemiological or genetic parameters. To investigate the role of host genetic background, we established long‐term experimental populations of different genotypes of the protozoan Paramecium caudatum, infected with the bacterial parasite Holospora undulata. We observed the evolution of an overall pattern of parasite local maladaptation for infectivity, indicating a general coevolutionary disadvantage of this parasite. Maladaptation extended to host populations with the same genetic background, similar to extending from the local to a higher regional level in natural populations. Patterns for virulence were qualitatively similar, but with less statistical support. A nonsignificant correlation with levels of (mal)adaptation for infectivity suggests independent evolution of these traits. Our results indicate similar (co)evolutionary trajectories in populations with different genetic backgrounds. Nonetheless, the correlated clines of genetic distance and parasite performance illustrate how genetic background can shape spatial gradients of local adaptation.     

THE GEOGRAPHY OF COEVOLUTION: COMPARATIVE POPULATION STRUCTURES FOR A SNAIL AND ITS TREMATODE PARASITE

Gene flow and the genetic structure of host and parasite populations are critical to the coevolutionary process, including the conditions under which antagonistic coevolution favors sexual reproduction. Here we compare the genetic structures of different populations of a freshwater New Zealand snail (Potamopyrgus antipodarum) with its trematode parasite (Microphallus sp.) using allozyme frequency data. Allozyme variation among snail populations was found to be highly structured among lakes; but for the parasite there was little allozyme structure among lake populations, suggesting much higher levels of parasite gene flow. The overall pattern of variation was confirmed with principal component analysis, which also showed that the organization of genetic differentiation for the snail (but not the parasite) was strongly related to the geographic arrangement of lakes. Some snail populations from different sides of the Alps near mountain passes were more similar to each other than to other snail populations on the same side of the Alps. Furthermore, genetic distances among parasite populations were correlated with the genetic distances among host populations, and genetic distances among both host and parasite populations were correlated with “stepping-stone” distances among lakes. Hence, the host snail and its trematode parasite seem to be dispersing to adjacent lakes in a stepping-stone fashion, although parasite dispersal among lakes is clearly greater. High parasite gene flow should help to continuously reintroduce genetic diversity within local populations where strong selection might otherwise isolate “host races.” Parasite gene flow can thereby facilitate the coevolutionary (Red Queen) dynamics that confer an advantage to sexual reproduction by restoring lost genetic variation.     

Genetic diversity,population structure and sex‐biased dispersal in three co‐evolving species

Genetic diversity and spatial structure of populations are important for antagonistic coevolution. We investigated genetic variation and population structure of three closely related European ant species: the social parasite Harpagoxenus sublaevis and its two host species Leptothorax acervorum and Leptothorax muscorum. We sampled populations in 12 countries and analysed eight microsatellite loci and an mtDNA sequence. We found high levels of genetic variation in all three species, only slightly less variation in the host L. muscorum. Using a newly introduced measure of differentiation (Jost’s D_est ), we detected strong population structuring in all species and less male‐biased dispersal than previously thought. We found no phylogeographic patterns that could give information on post‐glacial colonization routes – northern populations are as variable as more southern populations. We conclude that conditions for Thompson’s geographic mosaic of coevolution are ideal in this system: all three species show ample genetic variation and strong population structure.     

GENETIC VARIATION IN ASTERIONELLA FORMOSA (BACILLARIOPHYCEAE): IS IT LINKED TO FREQUENT EPIDEMICS OF HOST‐SPECIFIC PARASITIC FUNGI?1

Understanding of the genetic basis for susceptibility and resistance is still lacking for most aquatic host–parasite systems, for instance, for phytoplankton and their fungal parasites. Fungal parasites can have significant effects on phytoplankton populations, mainly through their ability to decimate algal host populations during epidemics. We used random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) analysis to study levels of genetic variation within a population of the freshwater diatom Asterionella formosa Hassall in relation to parasitism by the obligate, host‐specific, fungal parasite Zygorhizidium planktonicum Canter. The level of genetic variation within the A. formosa population in Lake Maarsseveen, The Netherlands was found to be high despite the presumed absence or very low frequency of sexual reproduction in this species, the limited gene flow, and the severity of parasite attack that would purge the population from susceptible genotypes. RAPD analysis revealed four distinct banding patterns, with 3 of 21 markers (14%) being polymorphic. In AFLP analysis, every single isolate of A. formosa showed a unique banding pattern, and 120 of the 210 AFLP markers (57%) were found to be polymorphic. Furthermore, character compatibility analysis revealed that sexual reproduction may be one of the mechanisms that generates and maintains genetic variation in the A. formosa population in Lake Maarsseveen. The presence of genetic variation in A. formosa was reflected in infection experiments, which showed that genetically different A. formosa strains differed in their susceptibility to various Z. planktonicum strains and that parasite strains differed in their ability to infect particular host strains.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

THE CELLULAR IMMUNE RESPONSE OF DAPHNIA MAGNA UNDER HOST–PARASITE GENETIC VARIATION AND VARIATION IN INITIAL DOSE

In invertebrate–parasite systems, the likelihood of infection following parasite exposure is often dependent on the specific combination of host and parasite genotypes (termed genetic specificity). Genetic specificity can maintain diversity in host and parasite populations and is a major component of the Red Queen hypothesis. However, invertebrate immune systems are thought to only distinguish between broad classes of parasite. Using a natural host–parasite system with a well‐established pattern of genetic specificity, the crustacean Daphnia magna and its bacterial parasite Pasteuria ramosa, we found that only hosts from susceptible host–parasite genetic combinations mounted a cellular response following exposure to the parasite. These data are compatible with the hypothesis that genetic specificity is attributable to barrier defenses at the site of infection (the gut), and that the systemic immune response is general, reporting the number of parasite spores entering the hemocoel. Further supporting this, we found that larger cellular responses occurred at higher initial parasite doses. By studying the natural infection route, where parasites must pass barrier defenses before interacting with systemic immune responses, these data shed light on which components of invertebrate defense underlie genetic specificity.     

Within‐ and among‐population variation in infectivity,latency and spore production in a host–pathogen system

In spatially structured populations, host–parasite coevolutionary potential depends on the distribution of genetic variation within and among populations. Inoculation experiments using the plant, Silene latifolia, and its fungal pathogen, Microbotryum violaceum, revealed little overall differentiation in infectivity/resistance, latency or spore production among host or pathogen populations. Within populations, fungal strains had similar means, but varied in performance across plant populations. Variation in resistance among seed families indicates the potential for parasite‐mediated selection, whereas there was little evidence for local pathogen genotype × plant genotype interactions assumed by most theoretical coevolution models. Lower spore production on sympatric than allopatric hosts confirmed local fungal maladaptation already observed for infectivity. Correlations between infectivity and latency or spore production suggest a common mechanism for variation in these traits. Our results suggest low variation available to this pathogen for tracking its coevolving host. This may be caused by random drift, breeding system or migration characteristic of metapopulation dynamics.     

Combining genome‐wide association study and FST‐based approaches to identify targets of Borrelia‐mediated selection in natural rodent hosts

Recent advances in high‐throughput sequencing technologies provide opportunities to gain novel insights into the genetic basis of phenotypic trait variation. Yet to date, progress in our understanding of genotype–phenotype associations in nonmodel organisms in general and natural vertebrate populations in particular has been hampered by small sample sizes typically available for wildlife populations and a resulting lack of statistical power, as well as a limited ability to control for false‐positive signals. Here we propose to combine a genome‐wide association study (GWAS) and F_ST‐based approach with population‐level replication to partly overcome these limitations. We present a case study in which we used this approach in combination with genotyping‐by‐sequencing (GBS) single nucleotide polymorphism (SNP) data to identify genomic regions associated with Borrelia afzelii resistance or susceptibility in the natural rodent host of this Lyme disease‐causing spirochete, the bank vole (Myodes glareolus). Using this combined approach we identified four consensus SNPs located in exonic regions of the genes Slc26a4, Tns3, Wscd1 and Espnl, which were significantly associated with the voles’ Borrelia infectious status within and across populations. Functional links between host responses to bacterial infections and most of these genes have previously been demonstrated in other rodent systems, making them promising new candidates for the study of evolutionary host responses to Borrelia emergence. Our approach is applicable to other systems and may facilitate the identification of genetic variants underlying disease resistance or susceptibility, as well as other ecologically relevant traits, in wildlife populations.     

Prioritization of SNPs for genome-wide association studies using an interaction model of genetic variation,gene expression,and trait variation

The identification of true causal loci to unravel the statistical evidence of genotype-phenotype correlations and the biological relevance of selected single-nucleotide polymorphisms (SNPs) is a challenging issue in genome-wide association studies (GWAS). Here, we introduced a novel method for the prioritization of SNPs based on p-values from GWAS. The method uses functional evidence from populations, including phenotype-associated gene expressions. Based on the concept of genetic interactions, such as perturbation of gene expression by genetic variation, phenotype and gene expression related SNPs were prioritized by adjusting the p-values of SNPs. We applied our method to GWAS data related to drug-induced cytotoxicity. Then, we prioritized loci that potentially play a role in druginduced cytotoxicity. By generating an interaction model, our approach allowed us not only to identify causal loci, but also to find intermediate nodes that regulate the flow of information among causal loci, perturbed gene expression, and resulting phenotypic variation.     

Higher parasite resistance in Daphnia populations with recent epidemics

Natural populations often show genetic variation in parasite resistance, forming the basis for evolutionary response to selection imposed by parasitism. We investigated whether previous epidemics selected for higher resistance to novel parasite isolates in a Daphnia galeata–microparasite system by comparing susceptibility of host clones from populations with varying epidemic history. We manipulated resource availability to evaluate whether diet influences Daphnia susceptibility as epidemics are common in nutrient‐rich lakes. Exposing clones from 10 lakes under two food treatments to an allopatric protozoan parasite, we found that Daphnia originating from lakes (mainly nutrient rich) with previous epidemics better resist infection. Despite this result, there was a tendency of higher susceptibility in the low food treatment, suggesting that higher resistance of clones from populations with epidemic background is not directly caused by lake nutrient level. Rather, our results imply that host populations respond to parasite‐mediated selection by evolving higher parasite resistance.     

Local parasite pressures and host genotype modulate epigenetic diversity in a mixed‐mating fish

Parasite‐mediated selection is one of the main drivers of genetic variation in natural populations. The persistence of long‐term self‐fertilization, however, challenges the notion that low genetic variation and inbreeding compromise the host's ability to respond to pathogens. DNA methylation represents a potential mechanism for generating additional adaptive variation under low genetic diversity. We compared genetic diversity (microsatellites and AFLPs), variation in DNA methylation (MS‐AFLPs), and parasite loads in three populations of Kryptolebias hermaphroditus, a predomintanly self‐fertilizing fish, to analyze the potential adaptive value of DNA methylation in relation to genetic diversity and parasite loads. We found strong genetic population structuring, as well as differences in parasite loads and methylation levels among sampling sites and selfing lineages. Globally, the interaction between parasites and inbreeding with selfing lineages influenced DNA methylation, but parasites seemed more important in determining methylation levels at the local scale.     

Survey of Helminth Parasites in Populations of <Emphasis Type="Italic">Alouatta palliata mexicana</Emphasis> and <Emphasis Type="Italic">A. pigra</Emphasis> in Continuous and in Fragmented Habitat in Southern Mexico

The 2 howler species that occur in southern Mexico, Alouatta palliata mexicana and Alouatta pigra are endangered, mainly as a result of habitat loss and fragmentation from human activity. Little is known about the gastrointestinal parasite communities affecting their populations, and lack of baseline information for populations of howler species in continuous forest habitats, makes evaluations of gastrointestinal parasite prevalence in populations in fragmented landscapes difficult. We report the results of a one-time broad survey of gastrointestinal parasites in fecal samples of individuals from several demographically stable populations of Alouatta palliata mexicana and A. pigra existing in continuous and/or protected forests. We further report similar data for populations of both species in human-fragmented landscapes. We detected 6 parasites for each howler monkey species, but only 3 of them (Trematode I, Controrchis biliophilus, Trypanoxyuris sp.) were common to both species. While parasitic prevalence in populations of both howler species was, in general, higher in the fragmented habitat than in continuous and/or protected forests. The difference is only marginally significant in Alouatta pigra. Some parasites (Coccidia and Strongylid) only appeared in populations in fragmented landscapes. Preliminary data suggest that adult males tended to have higher parasite prevalence values than those of adult females in both howler species. Parasite prevalence is associated to average group size, but not to population density in Alouatta pigra.     

Parasite infection reflects host genetic diversity among non-native populations of pumpkinseed sunfish in Europe

Species introductions often coincide with loss of genetic diversity and natural enemies. Anthropogenic translocation of the North-American pumpkinseed Lepomis gibbosus (L., 1758) (Centrarchidae) and its further spread have resulted in recent species establishment in most European countries. This study determines genetic differentiation of non-native European pumpkinseed populations and identifies how their genetic structure relates to the distribution and abundance of parasite species. Microsatellite analysis indicated presence of three genetic lineages, which were well supported by discriminant analysis based on parasite abundance data. The first lineage clustered pumpkinseed populations from northern and southern France and showed high allelic richness, heterozygosity and parasite richness. The second included populations along the “Southern invasion corridor” connecting the rivers Rhine, Main and Danube. The fish exhibited low to high genetic and parasite diversity and generally high parasite abundance. The third lineage clustered populations with low genetic and parasite diversity, located in Portuguese reservoirs and water bodies along the upper Elbe. Parasite species richness was significantly associated with host microsatellite heterozygosity and allelic richness, a trend partially affected by richness of North-American parasites. Furthermore, our results indicate that parasite community composition may serve as a useful biological tool to discriminate non-native fish populations and their inter-relationships.

     

Large‐scale GWAS in sorghum reveals common genetic control of grain size among cereals

Grain size is a key yield component of cereal crops and a major quality attribute. It is determined by a genotype’s genetic potential and its capacity to fill the grains. This study aims to dissect the genetic architecture of grain size in sorghum. An integrated genome‐wide association study (GWAS) was conducted using a diversity panel (n = 837) and a BC‐NAM population (n = 1421). To isolate genetic effects associated with genetic potential of grain size, rather than the genotype’s capacity to fill the grains, a treatment of removing half of the panicle was imposed during flowering. Extensive and highly heritable variation in grain size was observed in both populations in 5 field trials, and 81 grain size QTL were identified in subsequent GWAS. These QTL were enriched for orthologues of known grain size genes in rice and maize, and had significant overlap with SNPs associated with grain size in rice and maize, supporting common genetic control of this trait among cereals. Grain size genes with opposite effect on grain number were less likely to overlap with the grain size QTL from this study, indicating the treatment facilitated identification of genetic regions related to the genetic potential of grain size. These results enhance understanding of the genetic architecture of grain size in cereal, and pave the way for exploration of underlying molecular mechanisms and manipulation of this trait in breeding practices.