Mitogenome phylogeographic analysis of a planktonic crustacean.

Phylogeography places population genetics in an explicitly spatial context, and in doing so attempts to reconstruct the historical and contemporary evolutionary processes acting across a species range through space and time. Here we present the phylogeographical structure of Daphnia magna as determined for full mitochondrial genomes from samples of 60 populations throughout much of the species known range, including Europe, the Middle East, and Asia. Contrary to previous analyses, the present analysis of the mitochondrial genome reveals coarse-grained (continental scale) evidence for spatial structure, and in particular a deep split between Western Eurasia and East Asian D. magna lineages. In contrast to previous analyses with nuclear genetic markers, our mitogenomic analysis reveals much less structure within lineages. We quantify divergence between species using the full mitochondrial genome sequence of a closely related species, D. similis. The distribution of European and Middle Eastern genetic diversity is consistent with a rapid demographic expansion following the end of the most recent ice age about 10,000 years before present. By estimating species wide distributions of dN/dS in mtDNA, we provide evidence that the effectiveness of purifying selection on protein coding genes in the mitochondrial genome of coastal rock pool populations, which have pronounced extinction-colonization dynamics, is reduced compared to larger and more stable non-rock pool populations. The present study adds important insights into the evolutionary history of a widely used model organism in ecology, evolution and ecotoxicology, and highlights the utility of phylogeographic analysis of organellar genomes to understand evolutionary processes.

Transcriptome profiling during a natural host-parasite interaction.

BACKGROUND: Infection outcome in some coevolving host-pathogens is characterised by host-pathogen genetic interactions, where particular host genotypes are susceptible only to a subset of pathogen genotypes. To identify candidate genes responsible for the infection status of the host, we exposed a Daphnia magna host genotype to two bacterial strains of Pasteuria ramosa, one of which results in infection, while the other does not. At three time points (four, eight and 12 h) post pathogen exposure, we sequenced the complete transcriptome of the hosts using RNA-Seq (Illumina). RESULTS: We observed a rapid and transient response to pathogen treatment. Specifically, at the four-hour time point, eight genes were differentially expressed. At the eight-hour time point, a single gene was differentially expressed in the resistant combination only, and no genes were differentially expressed at the 12-h time point. CONCLUSIONS: We found that pathogen-associated transcriptional activity is greatest soon after exposure. Genome-wide resistant combinations were more likely to show upregulation of genes, while susceptible combinations were more likely to be downregulated, relative to controls. Our results also provide several novel candidate genes that may play a pivotal role in determining infection outcomes.

Novel insights into the insect trancriptome response to a natural DNA virus.

BACKGROUND: Little is known about invertebrate responses to DNA viruses. Here, we infect a commercially important pest moth species Plodia interpunctella with its naturally infecting DNA virus. We sequenced, assembled and annotated the complete transcriptome of the moth, and a partial transcriptome of the virus. We then tested for differential gene expression between moths that were exposed to the virus and controls. RESULTS: We found 51 genes that were differentially expressed in moths exposed to a DNA baculovirus compared to controls. Gene set enrichment analysis revealed that cuticle proteins were significantly overrepresented in this group of genes. Interestingly, 6 of the 7 differentially expressed cuticle proteins were downregulated, suggesting that baculoviruses are able to manipulate its host's response. In fact, an additional 29 of the 51 genes were also downregulated in exposed compared with control animals, including a gram-negative binding protein. In contrast, genes involved in transposable element movement were upregulated after infection. CONCLUSIONS: We present the first experiment to measure genome-wide gene expression in an insect after infection with a natural DNA virus. Our results indicate that cuticle proteins might be key genes underpinning the response to DNA viruses. Furthermore, the large proportion of genes that were downregulated after viral exposure suggests that this virus is actively manipulating the insect immune response. Finally, it appears that transposable element activity might increase during viral invasion. Combined, these results provide much needed host candidate genes that respond to DNA viral invaders.

The development of pathogen resistance in Daphnia magna: implications for disease spread in age-structured populations.

Immunity in vertebrates is well established to develop with time, but the ontogeny of defence in invertebrates is markedly less studied. Yet, age-specific capacity for defence against pathogens, coupled with age structure in populations, has widespread implications for disease spread. Thus, we sought to determine the susceptibility of hosts of different ages in an experimental invertebrate host-pathogen system. In a series of experiments, we show that the ability of Daphnia magna to resist its natural bacterial pathogen Pasteuria ramosa changes with host age. Clonal differences make it difficult to draw general conclusions, but the majority of observations indicate that resistance increases early in the life of D. magna, consistent with the idea that the defence system develops with time. Immediately following this, at about the time when a daphnid would be most heavily investing in reproduction, resistance tends to decline. Because many ecological factors influence the age structure of Daphnia populations, our results highlight a broad mechanism by which ecological context can affect disease epidemiology. We also show that a previously observed protective effect of restricted maternal food persists throughout the entire juvenile period, and that the protective effect of prior treatment with a small dose of the pathogen ('priming') persists for 7 days, observations that reinforce the idea that immunity in D. magna can change over time. Together, our experiments lead us to conclude that invertebrate defence capabilities have an ontogeny that merits consideration with respect to both their immune systems and the epidemic spread of infection.

Daphnia magna shows reduced infection upon secondary exposure to a pathogen.

Previous pathogen exposure is an important predictor of the probability of becoming infected. This is deeply understood for vertebrate hosts, and increasingly so for invertebrate hosts. Here, we test if an initial pathogen exposure changes the infection outcome to a secondary pathogen exposure in the natural host-pathogen system Daphnia magna and Pasteuria ramosa. Hosts were initially exposed to an infective pathogen strain, a non-infective pathogen strain or a control. The same hosts underwent a second exposure, this time to an infective pathogen strain, either immediately after the initial encounter or 48 h later. We observed that an initial encounter with a pathogen always conferred protection against infection compared with controls.

Immune genes undergo more adaptive evolution than non-immune system genes in Daphnia pulex.

BACKGROUND: Understanding which parts of the genome have been most influenced by adaptive evolution remains an unsolved puzzle. Some evidence suggests that selection has the greatest impact on regions of the genome that interact with other evolving genomes, including loci that are involved in host-parasite co-evolutionary processes. In this study, we used a population genetic approach to test this hypothesis by comparing DNA sequences of 30 putative immune system genes in the crustacean Daphnia pulex with 24 non-immune system genes. RESULTS: In support of the hypothesis, results from a multilocus extension of the McDonald-Kreitman (MK) test indicate that immune system genes as a class have experienced more adaptive evolution than non-immune system genes. However, not all immune system genes show evidence of adaptive evolution. Additionally, we apply single locus MK tests and calculate population genetic parameters at all loci in order to characterize the mode of selection (directional versus balancing) in the genes that show the greatest deviation from neutral evolution. CONCLUSIONS: Our data are consistent with the hypothesis that immune system genes undergo more adaptive evolution than non-immune system genes, possibly as a result of host-parasite arms races. The results of these analyses highlight several candidate loci undergoing adaptive evolution that could be targeted in future studies.

Length variation in 18S rRNA expansion segment 43/e4 of Daphnia obtusa: ancient or recurring polymorphism?

Expansion segments in ribosomal DNA (rDNA) can show length variation at the level of the individual, yet our understanding of the evolutionary forces shaping this variation is incomplete. Previous studies of expansion segment 43/e4 of the 18S rRNA gene in Daphnia obtusa have examined this variation in six individuals; however, it is not known if the variation documented at this locus is representative of variation across the species' geographic range. Furthermore, it is unclear whether length variants found in multiple individuals share common ancestry, or were generated de novo through recombination. We quantified expansion segment length variant frequencies in 134 individual D. obtusa from 33 populations at 15 sites across the species range in the US, and used a phylogeographic approach to determine whether recombination continues to add to the standing crop of variation at this locus. We identified seven length variants across the sampling range, which spans almost 3000 km. Based on the phylogeographic distribution of length variants in the expansion segment, we conclude that they are shared ancient polymorphisms that have persisted despite the operation of molecular mechanisms that cause the concerted evolution of multigene families such as rDNA.

An ancient immunity gene duplication in Daphnia magna: RNA expression and sequence analysis of two nitric oxide synthase genes.

NO (nitric oxide) is a highly reactive free radical gas thought to play a major role in the invertebrate immune response by harming pathogens and limiting their growth. Here we report on studies of nitric oxide synthase (NOS) genes in the crustacean Daphnia, one of the few non-insect arthropod models used to study host-pathogen interactions. While the NOS gene is found as a single copy in other invertebrates, we found two copies (NOS1 and NOS2), which a phylogenetic reconstruction showed to be the result of an ancient duplication event. Both genes bear features commonly found in invertebrate NOS, however, the two genes differ in their rate of evolution, intraspecific polymorphism and expression level. We tested whether the more rapid evolution of NOS2 could be due to positive selection, but found the rate of amino-acid substitutions between Daphnia species to be compatible with a neutral model. To associate NOS or NO activity with infection, we performed infection experiments with Daphnia magna and one of its natural pathogens (the bacterium Pasteuria ramosa). In one set of experimental infections, we supplemented D. magna with L-arginine, the NOS substrate, or with L-NAME, a NOS antagonist, and found this to result in lower and higher infection levels, respectively, which is at least compatible with the notion that NO may aid defence against Pasteuria. A second experiment indicated that NOS transcription does not increase following exposure to Pasteuria. Thus, the function of NOS in Daphnia immunity remains uncertain, but the pattern of gene duplication and subsequent divergence suggests evolution via neo- or subfunctionalization.

The components of the Daphnia pulex immune system as revealed by complete genome sequencing.

BACKGROUND: Branchiopod crustaceans in the genus Daphnia are key model organisms for investigating interactions between genes and the environment. One major theme of research on Daphnia species has been the evolution of resistance to pathogens and parasites, but lack of knowledge of the Daphnia immune system has limited the study of immune responses. Here we provide a survey of the immune-related genome of D. pulex, derived from the newly completed genome sequence. Genes likely to be involved in innate immune responses were identified by comparison to homologues from other arthropods. For each candidate, the gene model was refined, and we conducted an analysis of sequence divergence from homologues from other taxa. RESULTS AND CONCLUSION: We found that some immune pathways, in particular the TOLL pathway, are fairly well conserved between insects and Daphnia, while other elements, in particular antimicrobial peptides, could not be recovered from the genome sequence. We also found considerable variation in gene family copy number when comparing Daphnia to insects and present phylogenetic analyses to shed light on the evolution of a range of conserved immune gene families.

Nucleotide polymorphism and within-gene recombination in Daphnia magna and D. pulex, two cyclical parthenogens.

Theory predicts that partially asexual organisms may make the "best of both worlds": for the most part, they avoid the costs of sexual reproduction, while still benefiting from an enhanced efficiency of selection compared to obligately asexual organisms. There is, however, little empirical data on partially asexual organisms to test this prediction. Here we examine patterns of nucleotide diversity at eight nuclear loci in continentwide samples of two species of cyclically parthenogenetic Daphnia to assess the effect of partial asexual reproduction on effective population size and amount of recombination. Both species have high nucleotide diversities and show abundant evidence for recombination, yielding large estimates of effective population sizes (300,000-600,000). This suggests that selection will act efficiently even on mutations with small selection coefficients. Divergence between the two species is less than one-tenth of previous estimates, which were derived using a mitochondrial molecular clock. As the two species investigated are among the most distantly related species of the genus, this suggests that the genus Daphnia may be considerably younger than previously thought. Daphnia has recently received increased attention because it is being developed as a model organism for ecological and evolutionary genomics. Our results confirm the attractiveness of Daphnia as a model organism, because the high nucleotide diversity and low linkage disequilibrium suggest that fine-scale mapping of genes affecting phenotypes through association studies should be feasible.

Rates of recombination in the ribosomal DNA of apomictically propagated Daphnia obtusa lines.

Ribosomal (r)DNA undergoes concerted evolution, the mechanisms of which are unequal crossing over and gene conversion. Despite the fundamental importance of these mechanisms to the evolution of rDNA, their rates have been estimated only in a few model species. We estimated recombination rate in rDNA by quantifying the relative frequency of intraindividual length variants in an expansion segment of the 18S rRNA gene of the cladoceran crustacean, Daphnia obtusa, in four apomictically propagated lines. We also used quantitative PCR to estimate rDNA copy number. The apomictic lines were sampled every 5 generations for 90 generations, and we considered each significant change in the frequency distribution of length variants between time intervals to be the result of a recombination event. Using this method, we calculated the recombination rate for this region to be 0.02-0.06 events/generation on the basis of three different estimates of rDNA copy number. In addition, we observed substantial changes in rDNA copy number within and between lines. Estimates of haploid copy number varied from 53 to 233, with a mean of 150. We also measured the relative frequency of length variants in 30 lines at generations 5, 50, and 90. Although length variant frequencies changed significantly within and between lines, the overall average frequency of each length variant did not change significantly between the three generations sampled, suggesting that there is little or no bias in the direction of change due to recombination.

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa.

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.