Search

Many forest tree species have characteristics that make the study of their evolutionary ecology complex. For example, they are long-lived and thus have long generation times, and their often large, complex genomes have hampered establishing genomic resources. One way to tackle this challenge is to access multiple complementary data sources and analytical approaches when attempting to infer patterns of adaptive evolution. In the cover article of this issue of Molecular Ecology, Depardieu et al. (2021) combine large amounts of environmental, genomic, dendrochronological, and gene expression data in a common garden to explore the polygenic basis of drought resistance in white spruce (Picea glauca), a long-lived conifer. They identify candidate genes involved in growth and resistance to extreme drought events and show how multiple data sets may deliver complementary evidence to circumvent the manifold challenges of the research field.

Long generation times have been suggested to hamper rapid genetic adaptation of organisms to changing environmental conditions. We examined if environmental memory of the parental Scots pines (Pinus sylvestris L.) drive offspring survival and growth. We used seeds from trees growing under naturally dry conditions (control), irrigated trees (irrigated from 2003 to 2016), and formerly irrigated trees (“irrigation stop”; irrigated from 2003–2013; control condition since 2014). We performed two experiments, one under controlled greenhouse conditions and one at the experimental field site. In the greenhouse, the offspring from control trees exposed regularly to drought were more tolerant to hot–drought conditions than the offspring from irrigated trees and showed lower mortality even though there was no genetic difference. However, under optimal conditions (high water supply and full sunlight), these offspring showed lower growth and were outperformed by the offspring of the irrigated trees. This different offspring growth, with the offspring of the “irrigation‐stop” trees showing intermediate responses, points to the important role of transgenerational memory for the long‐term acclimation of trees. Such memory effects, however, may be overridden by climatic extremes during germination and early growth stages such as the European 2018 mega‐drought that impacted our field experiment., Environmental conditions (natural dry vs. irrigated) experienced by the mother trees of Pinus sylvestris significantly influenced the offspring's growth and survival, indicating that transgenerational environmental memory can shape trees' acclimation potential over the midterm to the long term.

International audience; There is growing concern that the implementation of political agreements on climate change and biodiversity will not be enough to protect forests in the short run and up to the end of the 21st century. As mitigation efforts are lagging behind self-imposed, reasonable targets, genetic diversity will have a large and significant part to play in the process of adapting forests to climate change. Genetic diversity, the raw material of evolution, can be used for adaptation by natural selection and artificial breeding, in naturally regenerated and plantation forests alike. The two-day scientific conference: “#rescueforests: Genetics to the rescue - Managing forests sustainably in a changing world”, addressed the genetic diversity of forests. More specifically, the conference was about natural as well as assisted adaptive processes, their spatial scale, from fine grain to landscape and ecoregions, and how much of the genome it involves. It also dealt with phenotypes and how much of their variation is determined by underlying genetic diversity. And finally, and perhaps most importantly, the conference emphasized the importance of conservation and sustainable use of this genetic diversity as a nature-based solution to adapt under the fast pace of climate change. The conference demonstrated how improved knowledge on genomic diversity and evolutionary mechanisms can help to rescue forests, either naturally or by means of management.

We present a conservation genetics tool kit, which offers two ready-to-use workflows for the routine application of genetic methods in conservation management. The workflows were optimized for work load and costs and are accompanied by an easy-to-read and richly illustrated manual with guidelines regarding sampling design, sampling of genetic material, necessary permits, laboratory methods, statistical analyses and documentation of results in a practice-oriented way. The manual also provides a detailed interpretation help for the implementation of the results in conservation management. One workflow deals with the identification of pond-breeding amphibians based on metabarcoding and environmental DNA (eDNA) from water samples. This workflow also discriminates the morphologically similar water frogs (Pelophylax sp.) and other closely related species (e.g. Triturus cristatus and T. carnifex). The second workflow studies connectivity among populations using microsatellite markers. Its statistical analyses encompass the detection of genetic groups and historical, recent and current dispersal and gene flow. Using the two workflows does not involve academic research institutes; they can be applied by environmental consultancies, laboratories from the private sector, governmental agencies or non-governmental organisations. These and additional conservation genetic workflows will hopefully foster the routine use of genetic methods in conservation management., Conservation Genetics Resources, 12 (4), ISSN:1877-7252, ISSN:1877-7260

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole‐genome pooled population sequencing to study genome‐wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over‐representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species‐specific.

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios. Measurement(s) growth ring • wood density Technology Type(s) measuring table • calculation Factor Type(s) tree species Sample Characteristic - Organism Betula pendula • Fagus sylvatica • Picea abies • Populus nigra • Pinus pinaster • Pinus sylvestris • Quercus petraea Sample Characteristic - Location Europe Published

The combination of highly resolved climatic and genomic data allows assessment of putative maladaptation of populations to climate change and can identify high-risk populations. Now, a study that accounts for migration and dispersal shows high maladaptation of a North American tree species in the northern and eastern distribution range.

International audience; During the last years, simple sequence repeats (SSRs, also known as microsatellites) and single-nucleotide polymorphisms (SNPs) have become the most popular molecular markers for describing neutral genetic variation in populations of a wide range of organisms. However, only a limited number of studies has focused on comparing the performance of these two types of markers for describing the underlying genetic structure of wild populations. Moreover, none of these studies targeted fungi, the group of organisms with one of the most complex reproductive strategies. We evaluated the utility of SSRs and SNPs for inferring the neutral genetic structure of Armillaria cepistipes (basidiomycetes) at different spatial scales. For that, 407 samples were collected across a small (150 km(2)) area in the Ukrainian Carpathians and a large (41 000 km(2)) area in the Swiss Alps. All isolates were analyzed at 17 SSR loci distributed throughout the whole genome and at 24 SNP loci located in different single-copy conserved genes. The two markers showed different patterns of structure within the two spatial scales studied. The multi-allelic SSR markers seemed to be best suited for detecting genetic structure in indigenous fungal populations at a rather small spatial scale (radius of ~50-100 km). The pattern observed at SNP markers rather reflected ancient divergence of distant (~1000 km) populations that in addition are separated by mountain ranges. Despite these differences, both marker types were suitable for detecting the weak genetic structure of the two A. cepistipes populations investigated.

It is generally accepted that the spatial distribution of neutral genetic diversity within a species' native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome-wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, FST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (He ) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.

Silver fir (Abies alba Mill.) is a keystone conifer of European montane forest ecosystems that has experienced large fluctuations in population size during during the Quaternary and, more recently, due to land-use change. To forecast the species’ future distribution and survival, it is important to investigate the genetic basis of adaptation to environmental change, notably to extreme events. For this purpose, we here provide a first draft genome assembly and annotation of the silver fir genome, established through a community-based initiative. DNA obtained from haploid megagametophyte and diploid needle tissue was used to construct and sequence Illumina paired-end and mate-pair libraries, respectively, to high depth. The assembled A. alba genome sequence accounted for over 37 million scaffolds corresponding to 18.16 Gb, with a scaffold N50 of 14,051 bp. Despite the fragmented nature of the assembly, a total of 50,757 full-length genes were functionally annotated in the nuclear genome. The chloroplast genome was also assembled into a single scaffold (120,908 bp) that shows a high collinearity with both the A. koreana and A. sibirica complete chloroplast genomes. This first genome assembly of silver fir is an important genomic resource that is now publicly available in support of a new generation of research. By genome-enabling this important conifer, this resource will open the gate for new research and more precise genetic monitoring of European silver fir forests.

Numerous landscape genomic studies have identified single-nucleotide polymorphisms (SNPs) and genes potentially involved in local adaptation. Rarely, it has been explicitly evaluated whether these environmental associations also hold true beyond the populations studied. We tested whether putatively adaptive SNPs in Arabidopsis halleri (Brassicaceae), characterized in a previous study investigating local adaptation to a highly heterogeneous environment, show the same environmental associations in an independent, geographically enlarged set of 18 populations. We analysed new SNP data of 444 plants with the same methodology (partial Mantel tests, PMTs) as in the original study and additionally with a latent factor mixed model (LFMM) approach. Of the 74 candidate SNPs, 41% (PMTs) and 51% (LFMM) were associated with environmental factors in the independent data set. However, only 5% (PMTs) and 15% (LFMM) of the associations showed the same environment–allele relationships as in the original study. In total, we found 11 genes (31%) containing the same association in the original and independent data set. These can be considered prime candidate genes for environmental adaptation at a broader geographical scale. Our results suggest that selection pressures in highly heterogeneous alpine environments vary locally and signatures of selection are likely to be population-specific. Thus, genotype-by-environment interactions underlying adaptation are more heterogeneous and complex than is often assumed, which might represent a problem when testing for adaptation at specific loci.

Despite decreasing sequencing costs, whole-genome sequencing for population-based genome scans for selection is still prohibitively expensive for organisms with large genomes. Moreover, the repetitive nature of large genomes often represents a challenge in bioinformatic and downstream analyses. Here we use in-depth transcriptome sequencing to design probes for exome capture in Swiss stone pine (Pinus cembra), a conifer with an estimated genome size of 29.3 Gbp and no reference genome available. We successfully applied around 55,000 self-designed probes, targeting 25,000 contigs, to DNA pools of seven populations from the Swiss Alps and identified > 140,000 SNPs in around 13,000 contigs. The probes performed equally well in pools of the closely related species Pinus sibirica; in both species, more than 70% of the targeted contigs were sequenced at a depth ≥ 40x, i.e. the number of haplotypes in the pool. However, a thorough analysis of individually sequenced P. cembra samples indicated that a majority of the contigs (63%) represented multi-copy genes. We therefore removed paralogous contigs based on heterozygote excess and deviation from allele balance. Without putatively paralogous contigs, allele frequencies of population pools represented accurate estimates of individually determined allele frequencies. Using population genetic and landscape genomic methods, we show that inferences of neutral and adaptive genetic variation may be biased when not accounting for such multi-copy genes. Future studies should therefore put more emphasis on identifying paralogous loci, which will be facilitated by the establishment of additional high-quality reference genomes.

When plants adapt to local environments, strong signatures of selection are expected in the genome, particularly in high-stress environments such as trace metal element enriched (metalliferous) soils. Using Arabidopsis halleri, a model species for metal homeostasis and adaptation to extreme environments, we identifid genes, gene variants, and pathways that are associated with soil properties and may thus contribute to adaptation to high concentrations of trace metal elements. We analysed whole-genome Pool-seq data from two metallicolous (from metalliferous soils) and two non-metallicolous populations (in total 119 individuals) and associated allele frequencies of the identified single-nucleotide polymorphisms (SNPs) with soil variables measured on site. Additionally, we accounted for polygenic adaptation by searching for gene pathways showing enrichment of signatures of selection. Out of >2.5 million SNPs, we identified 57 SNPs in 19 genes that were significantly associated with soil variables and are members of three enriched pathways. At least three of these candidate genes and pathways are involved in transmembrane transport and/or associated with responses to various stresses such as oxidative stress. We conclude that both allocation and detoxification processes play a crucial role in A. halleri for coping with these unfavourable conditions., Scientific Reports, 8, ISSN:2045-2322

Parasite communities of fishes are known to respond directly to the abiotic environment of the host, for example, to water quality and water temperature. Biotic factors are also important as they affect the exposure profile through heterogeneities in parasite distribution in the environment. Parasites in a particular environment may pose a strong selection on fish. For example, ecological differences in selection by parasites have been hypothesized to facilitate evolutionary differentiation of freshwater fish morphs specializing on different food types. However, as parasites may also respond directly to abiotic environment the parasite risk does not depend only on biotic features of the host environment. It is possible that different morphs experience specific selection gradients by parasites but it is not clear how consistent the selection is when abiotic factors change. We examined parasite pressure in sympatric morphs of threespine stickleback (Gasterosteus aculeatus) across a temperature gradient in two large Icelandic lakes, Myvatn and Thingvallavatn. Habitat‐specific temperature gradients in these lakes are opposite. Myvatn lava rock morph lives in a warm environment, while the mud morph lives in the cold. In Thingvallavatn, the lava rock morph lives in a cold environment and the mud morph in a warm habitat. We found more parasites in fish living in higher temperature in both lakes, independent of the fish morph, and this pattern was similar for the two dominating parasite taxa, trematodes and cestodes. However, at the same time, we also found higher parasite abundance in a third morph living in deep cold–water habitat in Thingvallavatn compared to the cold‐water lava morph, indicating strong effect of habitat‐specific biotic factors. Our results suggest complex interactions between water temperature and biotic factors in determining the parasite community structure, a pattern that may have implications for differentiation of stickleback morphs. ISSN:2045-7758

Testing how populations are locally adapted and predicting their response to their future environment is of key importance in view of climate change. Landscape genomics is a powerful approach to investigate genes and environmental factors involved in local adaptation. In a pooled amplicon sequencing approach of 94 genes in 71 populations, we tested if >3'500 single nucleotide polymorphisms (SNPs) in the three most common oak species in Switzerland (Quercus petraea, Q. pubescens, Q. robur) show an association with abiotic factors related to local topography, historical climate, and soil characteristics. In the analysis including all species, the most frequently associated environmental factors were those best describing the habitats of the species. In the species-specific analyses, the most important environmental factors and associated SNPs greatly differed among species. However, we identified one SNP and seven genes that were associated to the same environmental factor across all species. We finally used regressions of allele frequencies of the most strongly associated SNPs along environmental gradients to predict the risk of non-adaptedness (RONA), which represents the average change in allele frequency at climate-associated loci theoretically required to match future climatic conditions. RONA is considerable for some populations and species (up to 48% in single populations) and strongly differs among species. Given the long generation time of oaks, some of the required allele frequency changes might not be realistic to achieve based on standing genetic variation. Hence, future adaptedness requires gene flow or planting of individuals carrying beneficial alleles from habitats currently matching future climatic conditions. This article is protected by copyright. All rights reserved.

International audience; Abstract• Key messageWe show that joint multivariate analyses of leaf morphological characters and molecular-genetic markers improve the taxonomic assignment in hybridizing European white oaks. However, model-based approaches using genetic data alone represent straightforward alternatives to laborious, detailed morphological assessments.• ContextIn European white oaks, species delimitation is debated because of large overlap of morphological characteristics likely due to hybridization.• AimsWe tested whether joint multivariate analyses of leaf morphology and molecular markers improve the identification of three oak species (Quercus petraea, Quercus pubescens, Quercus robur) compared to approaches using morphological or genetic variables only.• MethodsWe assessed 13 leaf morphological characters and applied eight nuclear microsatellite markers in almost 1400 trees of 71 oak populations across Switzerland. We performed two multivariate approaches with three variable sets (morphology, genetics, combined) and assessed their performance in separating the taxa. We also compared the taxon assignment to a model-based clustering approach (Structure) based on genetic data alone.• ResultsA joint use of morphological and genetic variables led to an improved taxon assignment. Whereas Q. robur could clearly be separated from the two other taxa, there was a certain overlap between Q. petraea and Q. pubescens. The Structure clustering led to the same taxon assignment in 85 % of the individuals.• ConclusionIt is important to consider both morphological and genetic properties in morphologically similar and hybridizing species. However, it might be more efficient to concentrate only on genetic markers than on time-consuming morphological assessments.

Understanding genetic specificity in factors determining the outcome of host-parasite interactions is especially important as it contributes to parasite epidemiology, virulence, and maintenance of genetic variation. Such specificity, however, is still generally poorly understood. We examined genetic specificity in interactions among coinfecting parasites. In natural populations, individual hosts are often simultaneously infected by multiple parasite species and genotypes that interact. Such interactions could maintain genetic variation in parasite populations if they are genetically specific so that the relative fitness of parasite genotypes varies across host individuals depending on (1) the presence/absence of coinfections and/or (2) the genetic composition of the coinfecting parasite community. We tested these predictions using clones of fish eye flukes Diplostomum pseudospathaceum and Diplostomum gasterostei. We found that interactions among parasites had a strong genetic basis and that this modified genetic variation in infection success of D. pseudospathaceum between single and multiple infections as well as across multiply infected host individuals depending on the genetic identity of the coinfecting D. gasterostei. The relative magnitude of these effects, however, depended on the exposure dose, suggesting that ecological factors can modify genetic interactions between parasites.

Co-infecting parasite genotypes typically compete for host resources limiting their fitness. The intensity of such competition depends on whether parasites are reproducing in a host, or using it primarily as a transmission vehicle while not multiplying in host tissues (referred to as ‘competition hypothesis’). Alternatively, simultaneous attack and co-infection by several parasite genotypes might facilitate parasite infection because such a diverse attack could present an additional challenge to host immune defence (referred to as ‘facilitation hypothesis’). We tested the competition hypothesis by comparing the production of transmission stages (cercariae) from snails infected with one or two genotypes of the trematode Diplostomum pseudospathaceum . We found that cercarial production did not differ between the two groups of snails, suggesting lower per genotype production in double infections, and competition for host resources. Second, we tested the facilitation hypothesis by comparing parasite infection success on fishes (proportion of parasites establishing in the host) using cercariae originating from single-infected snails, double-infected snails and artificial mixtures of the single genotypes. In both cases, we found higher infection success when fishes were challenged with two parasite genotypes instead of one, supporting the facilitation hypothesis. Our results suggest that constraints defining the success of multiple genotype infections in parasites with multiple host life cycles include both between-genotype resource competition in the host and performance of host immune defences against a diverse parasite challenge.

Due to their particular way of life, dispersal of parasites is often mediated by their host's biology. Dispersal distance is relevant for parasites because high degree of dispersal leads to high gene flow, which counters the rate of parasite local adaptation in the host populations. Parasites with complex life cycles need to exploit sequentially more than one host species to complete their life cycle. Most trematode parasites have such complex life cycles involving invertebrate and vertebrate hosts. The spatial scales of invertebrate and vertebrate host populations are often different, which may decrease the probability that the parasite cycles locally in the intermediate host population. We used neutral microsatellite markers to determine genetic structure in Diplostomum pseudospathaceum parasites collected from local populations of freshwater snails (Lymnaea stagnalis). D. pseudospathaceum is a trematode that has two intermediate hosts (snail and fish) and a highly motile definitive host (bird). We found that the parasite population infecting the local snail populations showed no genetic structure over a large geographic range (over 300 km). We also did not detect evidence for isolation by distance in the parasite. We conclude that dispersal in the motile definitive host is likely to prevent emergence of local population genetic structure in the parasite. Our results suggest that parasite dispersal in the definitive host may limit local cycling of the parasites in the intermediate host populations.

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool-Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using FST-based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo-climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo-climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment., Molecular Ecology, 22 (22), ISSN:0962-1083, ISSN:1365-294X

Characterizing genetic variation in parasite transmission traits and its contribution to parasite vigor is essential for understanding the evolution of parasite life‐history traits. We measured genetic variation in output, activity, survival, and infection success of clonal transmission stages (cercaria larvae) of a complex life cycle parasite (Diplostomum pseudospathaceum). We further tested if variation in host nutritional stage had an effect on these traits by keeping hosts on limited or ad libitum diet. The traits we measured were highly variable among parasite genotypes indicating significant genetic variation in these life‐history traits. Traits were also phenotypically variable, for example, there was significant variation in the measured traits over time within each genotype. However, host nutritional stage had no effect on the parasite traits suggesting that a short‐term reduction in host resources was not limiting the cercarial output or performance. Overall, these results suggest significant interclonal and phenotypic variation in parasite transmission traits that are not affected by host nutritional status., Ecology and Evolution, 3 (7), ISSN:2045-7758

In nature, parasites commonly share hosts with other conspecific parasite genotypes. While adult parasites typically show aggregated distribution in their final hosts, aggregation of clonal parasite genotypes in intermediate hosts, such as those of trematodes in molluscs, is not generally known. However, infection of a host by multiple parasite genotypes has significant implications for evolution of virulence and host-parasite coevolution. Aggregated distribution of the clonal stages can increase host mortality and reduce larval output of each infecting genotype through interclonal competition, and therefore have significant implications for parasite epidemiology. The aim of this study was (i) to find out how common multiple genotype infections (MGIs) are in aquatic snails serving as intermediate hosts for different trematode species; (ii) to find out whether the prevalence of infection could be used to predict MGI frequencies and (iii) to use the relationship to infer whether MGIs aggregate in molluscan hosts. We determined the prevalence of trematode (Diplostomum pseudospathaceum) infections and the frequency of MGIs in snail (Lymnaea stagnalis) host populations as well as compiled corresponding literature data from a range of snail-trematode systems. We used Bayesian simulations to explore the relationship between prevalence of infection and MGI frequency in these data, and tested whether genotypes aggregate in snails by comparing the simulated relations with null model (Poisson and demographic Poisson) expectations. Our results show that MGIs are common in aquatic snails with up to 90% of the infected snails carrying MGIs. Parasite prevalence is a good predictor of MGI frequencies at a range of prevailing natural prevalences of infection (0-50%). The frequency of MGIs was higher than expected by both null models, indicating parasite aggregation at genotype level. These findings are in sharp contrast with the absence of multiple infections in snails at level of trematode species, suggesting that co-infections by multiple species and multiple genotypes of one species are controlled by different biological processes. Aggregation of MGIs in snail hosts appears to be widespread across different snail-trematode systems.

Infections by multiple parasite species are common in nature and have important consequences for between species interactions and coevolutionary dynamics with the host populations. For example ecological and evolutionary factors underlying the structure of parasite communities determine the range of hosts a parasite can infect and set the basis for both evolution of host defences and parasite virulence as well as management of diseases. Studies investigating these factors have been facilitated in the recent past by genetic methods which surmount difficulties of traditional morphological taxonomy in identifying individual parasite species. Here we take a step further and present a novel methodological approach to analyze parasite communities as a whole. We determined the relative frequencies of interspecific SNP alleles by pyrosequencing naturally pooled samples of closely related Diplostomum spp. trematodes infecting eye lenses of freshwater fish. Pyrosequencing allowed us to use naturally pooled community samples (lenses) to increase the sample size and statistical power without sequencing single parasite specimens. In a case study we applied this method to analyze the community structure of the eye flukes across fish host species of a freshwater system in Finland. We found that the fish species harboured significantly different parasite communities and that multiple species infections were common. Our study provides an example of how quantitative pyrosequencing can be used to answer evolutionary and ecological questions in natural communities of parasites. (C) 2011 Elsevier B.V. All rights reserved.