Negative frequency-dependent selection (NFDS) is an evolutionary mechanism suggested to govern host-parasite coevolution and the maintenance of genetic diversity at host resistance loci, such as the vertebrate MHC and R-genes in plants. Matching-allele interactions of hosts and parasites that prevent the emergence of host and parasite genotypes that are universally resistant and infective are a genetic mechanism predicted to underpin NFDS. The underlying genetics of matching-allele interactions are unknown even in host-parasite systems with empirical support for coevolution by NFDS, as is the case for the planktonic crustacean Daphnia magna and the bacterial pathogen Pasteuria ramosa. We fine-map one locus associated with D. magna resistance to P. ramosa and genetically characterize two haplotypes of the Pasteuria resistance (PR-) locus using de novo genome and transcriptome sequencing. Sequence comparison of PR-locus haplotypes finds dramatic structural polymorphisms between PR-locus haplotypes including a large portion of each haplotype being composed of non-homologous sequences resulting in haplotypes differing in size by 66 kb. The high divergence of PR-locus haplotypes suggest a history of multiple, diverse and repeated instances of structural mutation events and restricted recombination. Annotation of the haplotypes reveals striking differences in gene content. In particular, a group of glycosyltransferase genes that is present in the susceptible but absent in the resistant haplotype. Moreover, in natural populations, we find that the PR-locus polymorphism is associated with variation in resistance to different P. ramosa genotypes, pointing to the PR-locus polymorphism as being responsible for the matching-allele interactions that have been previously described for this system. Our results conclusively identify a genetic basis for the matching-allele interaction observed in a coevolving host-parasite system and provide a first insight into its molecular basis.

Dormancy is a common adaptation in invertebrates to survive harsh conditions. Triggered by environmental cues, populations produce resting eggs that allow them to survive temporally unsuitable conditions. Daphnia magna is a crustacean that reproduces by cyclical parthenogenesis, alternating between the production of asexual offspring and the sexual reproduction of diapausing eggs (ephippia). Prior to ephippia production, males (necessary to ensure ephippia fertilization) are produced parthenogenetically. Both the production of ephippia and the parthenogenetic production of males are induced by environmental factors. Here, we test the hypothesis that the induction of D. magna resting egg production shows a signature of local adaptation. We postulated that Daphnia from permanent ponds would produce fewer ephippia and males than Daphnia from intermittent ponds and that the frequency and season of habitat deterioration would correlate with the timing and amount of male and ephippia production. To test this, we quantified the production of males and ephippia in clonal D. magna populations in several different controlled environments. We found that the production of both ephippia and males varies strongly among populations in a way that suggests local adaptation. By performing quantitative trait locus mapping with parent clones from contrasting pond environments, we identified nonoverlapping genomic regions associated with male and ephippia production. As the traits are influenced by two different genomic regions, and both are necessary for successful resting egg production, we suggest that the genes for their induction co-evolve.

BackgroundAlthough Daphnia is increasingly recognized as a model for ecological genomics and biomedical research, there is, as of yet, no high-resolution genetic map for the genus. Such a map would provide an important tool for mapping phenotypes and assembling the genome. Here we estimate the genome size of Daphnia magna and describe the construction of an SNP array based linkage map. We then test the suitability of the map for life history and behavioural trait mapping. The two parent genotypes used to produce the map derived from D. magna populations with and without fish predation, respectively and are therefore expected to show divergent behaviour and life-histories.ResultsUsing flow cytometry we estimated the genome size of D. magna to be about 238 mb. We developed an SNP array tailored to type SNPs in a D. magna F2 panel and used it to construct a D. magna linkage map, which included 1,324 informative markers. The map produced ten linkage groups ranging from 108.9 to 203.6 cM, with an average distance between markers of 1.13 cM and a total map length of 1,483.6 cM (Kosambi corrected). The physical length per cM is estimated to be 160 kb. Mapping infertility genes, life history traits and behavioural traits on this map revealed several significant QTL peaks and showed a complex pattern of underlying genetics, with different traits showing strongly different genetic architectures.ConclusionsThe new linkage map of D. magna constructed here allowed us to characterize genetic differences among parent genotypes from populations with ecological differences. The QTL effect plots are partially consistent with our expectation of local adaptation under contrasting predation regimes. Furthermore, the new genetic map will be an important tool for the Daphnia research community and will contribute to the physical map of the D. magna genome project and the further mapping of phenotypic traits. The clones used to produce the linkage map are maintained in a stock collection and can be used for mapping QTLs of traits that show variance among the F2 clones.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-1033) contains supplementary material, which is available to authorized users.

Understanding the genetic architecture of host resistance is key for understanding the evolution of host-parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host-parasite interactions. In the QTL panel used here, Daphnia's resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host-parasite systems. Only the Pasteuria-Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium-Daphnia system remains unclear.