Host-parasite co-evolution and its genomic signature.

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

QTL study reveals candidate genes underlying host resistance in a Red Queen model system.

Specific interactions of host and parasite genotypes can lead to balancing selection, maintaining genetic diversity within populations. In order to understand the drivers of such specific coevolution, it is necessary to identify the molecular underpinnings of these genotypic interactions. Here, we investigate the genetic basis of resistance in the crustacean host, Daphnia magna, to attachment and subsequent infection by the bacterial parasite, Pasteuria ramosa. We discover a single locus with Mendelian segregation (3:1 ratio) with resistance being dominant, which we call the F locus. We use QTL analysis and fine mapping to localize the F locus to a 28.8-kb region in the host genome, adjacent to a known resistance supergene. We compare the 28.8-kb region in the two QTL parents to identify differences between host genotypes that are resistant versus susceptible to attachment and infection by the parasite. We identify 13 genes in the region, from which we highlight eight biological candidates for the F locus, based on presence/absence polymorphisms and differential gene expression. The top candidates include a fucosyltransferase gene that is only present in one of the two QTL parents, as well as several Cladoceran-specific genes belonging to a large family that is represented in multiple locations of the host genome. Fucosyltransferases have been linked to resistance in previous studies of Daphnia-Pasteuria and other host-parasite systems, suggesting that P. ramosa spore attachment could be mediated by changes in glycan structures on D. magna cuticle proteins. The Cladoceran-specific candidate genes suggest a resistance strategy that relies on gene duplication. Our results add a new locus to a growing genetic model of resistance in the D. magna-P. ramosa system. The identified candidate genes will be used in future functional genetic studies, with the ultimate aim to test for cycles of allele frequencies in natural populations.

Genome-Wide Allele Frequency Changes Reveal That Dynamic Metapopulations Evolve Differently.

Two important characteristics of metapopulations are extinction-(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation-selection-drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations.

Unraveling coevolutionary dynamics using ecological genomics.

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Multispecies coexistence in fragmented landscapes.

Spatial dynamics have long been recognized as an important driver of biodiversity. However, our understanding of species' coexistence under realistic landscape configurations has been limited by lack of adequate analytical tools. To fill this gap, we develop a spatially explicit metacommunity model of multiple competing species and derive analytical criteria for their coexistence in fragmented heterogeneous landscapes. Specifically, we propose measures of niche and fitness differences for metacommunities, which clarify how spatial dynamics and habitat configuration interact with local competition to determine coexistence of species. We parameterize our model with a Bayesian approach using a 36-y time-series dataset of three Daphnia species in a rockpool metacommunity covering >500 patches. Our results illustrate the emergence of interspecific variation in extinction and recolonization processes, including their dependencies on habitat size and environmental temperature. We find that such interspecific variation contributes to the coexistence of Daphnia species by reducing fitness differences and increasing niche differences. Additionally, our parameterized model allows separating the effects of habitat destruction and temperature change on species extinction. By integrating coexistence theory and metacommunity theory, our study provides platforms to increase our understanding of species' coexistence in fragmented heterogeneous landscapes and the response of biodiversity to environmental changes.

The cuticle proteome of a planktonic crustacean.

The cuticles of arthropods provide an interface between the organism and its environment. Thus, the cuticle's structure influences how the organism responds to and interacts with its surroundings. Here, we used label-free quantification proteomics to provide a proteome of the moulted cuticle of the aquatic crustacean Daphnia magna, which has long been a prominent subject of studies on ecology, evolution, and developmental biology. We detected a total of 278 high-confidence proteins. Using protein sequence domain and functional enrichment analyses, we identified chitin-binding structural proteins and chitin-modifying enzymes as the most abundant protein groups in the cuticle proteome. Structural cuticular protein families showed a similar distribution to those found in other arthropods and indicated proteins responsible for the soft and flexible structure of the Daphnia cuticle. Finally, cuticle protein genes were also clustered as tandem gene arrays in the D. magna genome. The cuticle proteome presented here will be a valuable resource to the Daphnia research community, informing genome annotations and investigations on diverse topics such as the genetic basis of interactions with predators and parasites.

Uncovering the Genomic Basis of Infection Through Co-genomic Sequencing of Hosts and Parasites.

Understanding the genomic basis of infectious disease is a fundamental objective in co-evolutionary theory with relevance to healthcare, agriculture, and epidemiology. Models of host-parasite co-evolution often assume that infection requires specific combinations of host and parasite genotypes. Co-evolving host and parasite loci are, therefore, expected to show associations that reflect an underlying infection/resistance allele matrix, yet little evidence for such genome-to-genome interactions has been observed among natural populations. We conducted a study to search for this genomic signature across 258 linked host (Daphnia magna) and parasite (Pasteuria ramosa) genomes. Our results show a clear signal of genomic association between multiple epistatically interacting loci in the host genome, and a family of genes encoding for collagen-like protein in the parasite genome. These findings are supported by laboratory-based infection trials, which show strong correspondence between phenotype and genotype at the identified loci. Our study provides clear genomic evidence of antagonistic co-evolution among wild populations.

Sperm competition favours intermediate sperm size in a hermaphrodite1.

Sperm competition is a potent mechanism of postcopulatory sexual selection that has been found to shape reproductive morphologies and behaviours in promiscuous animals. Especially sperm size has been argued to evolve in response to sperm competition through its effect on sperm longevity, sperm motility, the ability to displace competing sperm, and ultimately fertilization success. Additionally, sperm size has been observed to co-evolve with female reproductive morphology. Theoretical work predicts that sperm competition may select for longer sperm but may also favour shorter sperm if sperm size trades-off with number. In this study, we studied the relationship between sperm size and postmating success in the free-living flatworm, Macrostomum lignano. Specifically, we used inbred isolines of M. lignano that varied in sperm size to investigate how sperm size translated into the ability of worms to transfer and deposit sperm in a mating partner. Our results revealed a hump-shaped relationship with individuals producing sperm of intermediate size having the highest sperm competitiveness. This finding broadens our understanding of the evolution of sperm morphology by providing empirical support for stabilizing selection on sperm size under sperm competition.

Evolutionary change in flight-to-light response in urban moths comes with changes in wing morphology.

Moths and other insects are attracted by artificial light sources. This flight-to-light behaviour disrupts their general activity focused on finding resources, such as mating partners, and increases predation risk. It thus has substantial fitness costs. In illuminated urban areas, spindle ermine moths Yponomeuta cagnagella were reported to have evolved a reduced flight-to-light response. Yet, the specific mechanism remained unknown, and was hypothesized to involve either changes in visual perception or general flight ability or overall mobility traits. Here, we test whether spindle ermine moths from urban and rural populations-with known differences in flight-to-light responses-differ in flight-related morphological traits. Urban individuals were found to have on average smaller wings than rural moths, which in turn correlated with a lower probability of being attracted to an artificial light source. Our finding supports the reduced mobility hypothesis, which states that reduced mobility in urban areas is associated with specific morphological changes in the flight apparatus.

Genetic Drift Shapes the Evolution of a Highly Dynamic Metapopulation.

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

Population-Genomic Analysis Identifies a Low Rate of Global Adaptive Fixation in the Proteins of the Cyclical Parthenogen Daphnia magna.

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative 'arms race' genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution.

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the "Gymnomera"). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.].

The limits of stress-tolerance for zooplankton resting stages in freshwater ponds.

In seasonal environments, many organisms evolve strategies such as diapause to survive stressful periods. Understanding the link between habitat stability and diapause strategy can help predict a population's survival in a changing world. Indeed, resting stages may be an important way freshwater organisms can survive periods of drought or freezing, and as the frequency and extent of drought or freezing vary strongly among habitats and are predicted to change with climate change, it raises questions about how organisms cope with, and survive, environmental stress. Using Daphnia magna as a model system, we tested the ability of resting stages from different populations to cope with stress during diapause. The combination of elevated temperatures and wet conditions during diapause shows to prevent hatching altogether. In contrast, hatching is relatively higher after a dry and warm diapause, but declines with rising temperatures, while time to hatch increases. Resting stages produced by populations from summer-dry habitats perform slightly, but consistently, better at higher temperatures and dryness, supporting the local adaptation hypothesis. A higher trehalose content in resting eggs from summer-dry habitat might explain such pattern. Considering that temperatures and summer droughts are projected to increase in upcoming years, it is fundamental to know how resting stages resist stressful conditions so as to predict and protect the ecological functioning of freshwater ecosystems.

Balancing Selection for Pathogen Resistance Reveals an Intercontinental Signature of Red Queen Coevolution.

The link between long-term host-parasite coevolution and genetic diversity is key to understanding genetic epidemiology and the evolution of resistance. The model of Red Queen host-parasite coevolution posits that high genetic diversity is maintained when rare host resistance variants have a selective advantage, which is believed to be the mechanistic basis for the extraordinarily high levels of diversity at disease-related genes such as the major histocompatibility complex in jawed vertebrates and R-genes in plants. The parasites that drive long-term coevolution are, however, often elusive. Here we present evidence for long-term balancing selection at the phenotypic (variation in resistance) and genomic (resistance locus) level in a particular host-parasite system: the planktonic crustacean Daphnia magna and the bacterium Pasteuria ramosa. The host shows widespread polymorphisms for pathogen resistance regardless of geographic distance, even though there is a clear genome-wide pattern of isolation by distance at other sites. In the genomic region of a previously identified resistance supergene, we observed consistent molecular signals of balancing selection, including higher genetic diversity, older coalescence times, and lower differentiation between populations, which set this region apart from the rest of the genome. We propose that specific long-term coevolution by negative-frequency-dependent selection drives this elevated diversity at the host's resistance loci on an intercontinental scale and provide an example of a direct link between the host's resistance to a virulent pathogen and the large-scale diversity of its underlying genes.

A Two-Locus System with Strong Epistasis Underlies Rapid Parasite-Mediated Evolution of Host Resistance.

Parasites are a major evolutionary force, driving adaptive responses in host populations. Although the link between phenotypic response to parasite-mediated natural selection and the underlying genetic architecture often remains obscure, this link is crucial for understanding the evolution of resistance and predicting associated allele frequency changes in the population. To close this gap, we monitored the response to selection during epidemics of a virulent bacterial pathogen, Pasteuria ramosa, in a natural host population of Daphnia magna. Across two epidemics, we observed a strong increase in the proportion of resistant phenotypes as the epidemics progressed. Field and laboratory experiments confirmed that this increase in resistance was caused by selection from the local parasite. Using a genome-wide association study, we built a genetic model in which two genomic regions with dominance and epistasis control resistance polymorphism in the host. We verified this model by selfing host genotypes with different resistance phenotypes and scoring their F1 for segregation of resistance and associated genetic markers. Such epistatic effects with strong fitness consequences in host-parasite coevolution are believed to be crucial in the Red Queen model for the evolution of genetic recombination.

Demographic history shapes genomic variation in an intracellular parasite with a wide geographical distribution.

Analysing variation in a species' genomic diversity can provide insights into its historical demography, biogeography and population structure, and thus its ecology and evolution. Although such studies are rarely undertaken for parasites, they can be highly revealing because of the parasite's co-evolutionary relationships with hosts. Modes of reproduction and transmission are thought to be strong determinants of genomic diversity for parasites and vary widely among microsporidia (fungal-related intracellular parasites), which are known to have high intraspecific genetic diversity and interspecific variation in genome architecture. Here we explore genomic variation in the microsporidium Hamiltosporidium, a parasite of the freshwater crustacean Daphnia magna, looking especially at which factors contribute to nucleotide variation. Genomic samples from 18 Eurasian populations and a new, long-read-based reference genome were used to determine the roles that reproduction mode, transmission mode and geography play in determining population structure and demographic history. We demonstrate two main Hamiltosporidium tvaerminnensis lineages and a pattern of isolation-by-distance, but note an absence of congruence between these two parasite lineages and the two Eurasian host lineages. We suggest a comparatively recent parasite spread through Northern Eurasian host populations after a change from vertical to mixed-mode transmission and the loss of sexual reproduction. While gaining knowledge about the ecology and evolution of this focal parasite, we also identify common features that shape variation in genomic diversity for many parasites, such as distinct modes of reproduction and the intertwining of host-parasite demographies.

A new microsporidian parasite, Ordospora pajunii sp. nov (Ordosporidae), of Daphnia longispina highlights the value of genomic data for delineating species boundaries.

Speciation is a complex and continuous process that makes the delineation of species boundaries a challenging task in particular in species with little morphological differentiation, such as parasites. In this case, the use of genomic data is often necessary, such as for the intracellular Microsporidian parasites. Here, we characterize the genome of a gut parasite of the cladoceran Daphnia longispina (isolate FI-F-10), which we propose as a new species within the genus Ordospora: Ordospora pajunii sp. nov (Ordosporidae). FI-F-10 closest relative, Ordospora colligata is only found in D. magna. Both microsporidian species share several morphological features. Although it is not possible to estimate divergence times for Microsporidia due to the lack of fossil records and accelerated evolutionary rates, we base our proposal on the phylogenomic and genomic distances between both microsporidian lineages. The phylogenomic reconstruction shows that FI-F-10 forms an early diverging branch basal to the cluster that contains all known O. colligata strains. Whole-genome comparisons show that FI-F-10 presents a greater divergence at the sequence level than observed among O. colligata strains, and its genomic average nucleotide identity (ANI) values against O. colligata are beyond the intraspecific range previously established for yeast and prokaryotes. Our data confirm that the ANI metrics are useful for fine genetic divergence calibration across Microsporidia taxa. In combination with phylogenetic and ecological data, genome-based metrics provide a powerful approach to delimitate species boundaries.

Trehalose provisioning in Daphnia resting stages reflects local adaptation to the harshness of diapause conditions.

Environmental fluctuations often select for adaptations such as diapause states, allowing species to outlive harsh conditions. The natural sugar trehalose which provides both cryo- and desiccation-protection, has been found in diapause stages of diverse taxa. Here, we hypothesize that trehalose deposition in resting stages is a locally adapted trait, with higher concentrations produced in harsher habitats. We used resting stages, produced under standardized conditions, by 37 genotypes of Daphnia magna collected from Western Palaearctic habitats varying in their propensity to dry in summer and freeze in winter. Resting eggs produced by D. magna from populations from summer-dry habitats showed significantly higher trehalose than those from summer-wet habitats, suggesting that trehalose has a protective function during desiccation. By contrast, winter-freezing did not explain variation in trehalose content. Adaptations to droughts are important, as summer dryness of water bodies is foreseen to increase with ongoing climate change.

Reduced structural connectivity in the corpus callosum in patients with anorexia nervosa.

OBJECTIVE: Previous diffusion tensor imaging studies reported a reduced fractional anisotropy in the body of the corpus callosum in patients with anorexia nervosa, which may indicate impaired white matter integrity in interhemispheric connections. The aim of the current study was to investigate whether structural connectivity is affected in patients with anorexia nervosa. METHOD: To this end, we compared the number of streamlines (a model of the white matter fibre tracts) and the total volume filled by these streamlines in different subsections of the corpus callosum in 33 women with and 33 without anorexia nervosa as well as in 20 recovered individuals. RESULTS: The volume of streamlines in the anterior and mid-anterior subsection of the corpus callosum was reduced in women with, but not in those who had recovered from anorexia nervosa. No differences in number of streamlines was detected in the corpus callosum between patients with anorexia nervosa, healthy controls and recovered patients. CONCLUSIONS: Alterations of the corpus callosum have been repeatedly reported in anorexia nervosa. Since the recovered group did not differ from the healthy control group, the reported alterations in acute patients appear to represent a state and not a trait marker.

Genome-Wide Association Analysis Identifies a Genetic Basis of Infectivity in a Model Bacterial Pathogen.

Knowledge of the genetic architecture of pathogen infectivity and host resistance is essential for a mechanistic understanding of coevolutionary processes, yet the genetic basis of these interacting traits remains unknown for most host-pathogen systems. We used a comparative genomic approach to explore the genetic basis of infectivity in Pasteuria ramosa, a Gram-positive bacterial pathogen of planktonic crustaceans that has been established as a model for studies of Red Queen host-pathogen coevolution. We sequenced the genomes of a geographically, phenotypically, and genetically diverse collection of P. ramosa strains and performed a genome-wide association study to identify genetic correlates of infection phenotype. We found multiple polymorphisms within a single gene, Pcl7, that correlate perfectly with one common and widespread infection phenotype. We then confirmed this perfect association via Sanger sequencing in a large and diverse sample set of P. ramosa clones. Pcl7 codes for a collagen-like protein, a class of adhesion proteins known or suspected to be involved in the infection mechanisms of a number of important bacterial pathogens. Consistent with expectations under Red Queen coevolution, sequence variation of Pcl7 shows evidence of balancing selection, including extraordinarily high diversity and absence of geographic structure. Based on structural homology with a collagen-like protein of Bacillus anthracis, we propose a hypothesis for the structure of Pcl7 and the physical location of the phenotype-associated polymorphisms. Our results offer strong evidence for a gene governing infectivity and provide a molecular basis for further study of Red Queen dynamics in this model host-pathogen system.