Epigenetic age estimation of wild mice using faecal samples.

Age is a key parameter in population ecology, with a myriad of biological processes changing with age as organisms develop in early life then later senesce. As age is often hard to accurately measure with non-lethal methods, epigenetic methods of age estimation (epigenetic clocks) have become a popular tool in animal ecology and are often developed or calibrated using captive animals of known age. However, studies typically rely on invasive blood or tissue samples, which limit their application in more sensitive or elusive species. Moreover, few studies have directly assessed how methylation patterns and epigenetic age estimates compare across environmental contexts (e.g. captive or laboratory-based vs. wild animals). Here, we built a targeted epigenetic clock from laboratory house mice (strain C57BL/6, Mus musculus) using DNA from non-invasive faecal samples, and then used it to estimate age in a population of wild mice (Mus musculus domesticus) of unknown age. This laboratory mouse-derived epigenetic clock accurately predicted adult wild mice to be older than juveniles and showed that wild mice typically increased in epigenetic age over time, but with wide variation in epigenetic ageing rate among individuals. Our results also suggested that, for a given body mass, wild mice had higher methylation across targeted CpG sites than laboratory mice (and consistently higher epigenetic age estimates as a result), even among the smallest, juvenile mice. This suggests wild and laboratory mice may display different CpG methylation levels from very early in life and indicates caution is needed when developing epigenetic clocks on laboratory animals and applying them in the wild.

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.

Daphnia magna egg piRNA cluster expression profiles change as mothers age.

BACKGROUND: PiRNAs prevent transposable elements wreaking havoc on the germline genome. Changes in piRNA expression over the lifetime of an individual may impact on ageing through continued suppression, or release, of transposable element expression. We identified piRNA producing clusters in the genome of Daphnia magna by a combination of bioinformatic methods, and then contrasted their expression between parthenogenetically produced eggs representing maternally-deposited germline piRNAs of young (having their 1st clutch) and old (having their 5th clutch) mothers. Results from eggs were compared to cluster expression in three generations of adults. RESULTS: As for other arthropods, D. magna encodes long uni-directionally transcribed non-coding RNAs consisting of fragmented transposable elements which account for most piRNAs expressed. Egg tissues showed extensive differences between clutches from young mothers and those from old mothers, with 578 and 686 piRNA clusters upregulated, respectively. Most log fold-change differences for significant clusters were modest, however. When considering only highly expressed clusters, there was a bias towards 1st clutch eggs at 41 upregulated versus eight clusters in the eggs from older mothers. F0 generation differences between young and old mothers were fewer than eggs, as 179 clusters were up-regulated in young versus 170 old mothers. This dropped to 31 versus 22 piRNA clusters when comparing adults in the F1 generation, and no differences were detected in the F3 generation. Inter-generational losses of differential piRNA cluster were similar to that observed for D. magna micro-RNA expression. CONCLUSIONS: Little overlap in differentially expressed clusters was found between adults containing mixed somatic and germline (ovary) tissues and germ-line representing eggs. A cluster encompassing a Tudor domain containing gene important in the piRNA pathway was upregulated in the eggs from old mothers. We hypothesise that regulation of this gene could form part of a feedback loop that reduces piRNA pathway activity explaining the reduced number of highly-expressed clusters in eggs from old mothers.

Pathogen Dynamics across the Diversity of Aging.

AbstractReproduction, mortality, and immune function often change with age but do not invariably deteriorate. Across the tree of life, there is extensive variation in age-specific performance and changes to key life-history traits. These changes occur on a spectrum from classic senescence, where performance declines with age, to juvenescence, where performance improves with age. Reproduction, mortality, and immune function are also important factors influencing the spread of infectious disease, yet there exists no comprehensive investigation into how the aging spectrum of these traits impacts epidemics. We used a model laboratory infection system to compile an aging profile of a single organism, including traits directly linked to pathogen susceptibility and those that should indirectly alter pathogen transmission by influencing demography. We then developed generalizable epidemiological models demonstrating that different patterns of aging produce dramatically different transmission landscapes: in many cases, aging can reduce the probability of epidemics, but it can also promote severity. This work provides context and tools for use across taxa by empiricists, demographers, and epidemiologists, advancing our ability to accurately predict factors contributing to epidemics or the potential repercussions of senescence manipulation.

Daphnia magna modifies its gene expression extensively in response to caloric restriction revealing a novel effect on haemoglobin isoform preference.

Caloric restriction (CR) produces clear phenotypic effects within and between generations of the model crustacean Daphnia magna. We have previously established that micro-RNAs and cytosine methylation change in response to CR in this organism, and we demonstrate here that CR has a dramatic effect on gene expression. Over 6,000 genes were differentially expressed between CR and well-fed D. magna, with a bias towards up-regulation of genes under caloric restriction. We identified a highly expressed haemoglobin gene that responds to CR by changing isoform proportions. Specifically, a transcript containing three haem-binding erythrocruorin domains was strongly down-regulated under CR in favour of transcripts containing fewer or no such domains. This change in the haemoglobin mix is similar to the response to hypoxia in Daphnia, which is mediated through the transcription factor hypoxia-inducible factor 1, and ultimately the mTOR signalling pathway. This is the first report of a role for haemoglobin in the response to CR. We also observed high absolute expression of superoxide dismutase (SOD) in normally fed individuals, which contrasts with observations of high SOD levels under CR in other taxa. However, key differentially expressed genes, like SOD, were not targeted by differentially expressed micro-RNAs. Whether the link between haemoglobin and CR occurs in other organisms, or is related to the aquatic lifestyle, remains to be tested. It suggests that one response to CR may be to simply transport less oxygen and lower respiration.

Genome-wide methylation is modified by caloric restriction in Daphnia magna.

BACKGROUND: The degradation of epigenetic control with age is associated with progressive diseases of ageing, including cancers, immunodeficiency and diabetes. Reduced caloric intake slows the effects of ageing and age-related disease in vertebrates and invertebrates, a process potentially mediated by the impact of caloric restriction on epigenetic factors such as DNA methylation. We used whole genome bisulphite sequencing to study how DNA methylation patterns change with diet in a small invertebrate, the crustacean Daphnia magna. Daphnia show the classic response of longer life under caloric restriction (CR), and they reproduce clonally, which permits the study of epigenetic changes in the absence of genetic variation. RESULTS: Global cytosine followed by guanine (CpG) methylation was 0.7-0.9%, and there was no difference in overall methylation levels between normal and calorie restricted replicates. However, 333 differentially methylated regions (DMRs) were evident between the normally fed and CR replicates post-filtering. Of these 65% were hypomethylated in the CR group, and 35% were hypermethylated in the CR group. CONCLUSIONS: Our results demonstrate an effect of CR on the genome-wide methylation profile. This adds to a growing body of research in Daphnia magna that demonstrate an epigenomic response to environmental stimuli. Specifically, gene Ontology (GO) term enrichment of genes associated with hyper and hypo-methylated DMRs showed significant enrichment for methylation and acyl-CoA dehydrogenase activity, which are linked to current understanding of their roles in CR in invertebrate model organisms.

Caging and Uncaging Genetics.

It is important for biology to understand if observations made in highly reductionist laboratory settings generalise to harsh and noisy natural environments in which genetic variation is sorted to produce adaptation. But what do we learn by studying, in the laboratory, a genetically diverse population that mirrors the wild? What is the best design for studying genetic variation? When should we consider it at all? The right experimental approach depends on what you want to know.

Daphnia magna microRNAs respond to nutritional stress and ageing but are not transgenerational.

Maternal effects, where the performance of offspring is determined by the condition of their mother, are widespread and may in some cases be adaptive. The crustacean Daphnia magna shows strong maternal effects: offspring size at birth and other proxies for fitness are altered when their mothers are older or when mothers have experienced dietary restriction. The mechanisms for this transgenerational transmission of maternal experience are unknown, but could include changes in epigenetic patterning. MicroRNAs (miRNAs) are regulators of gene expression that have been shown to play roles in intergenerational information transfer, and here, we test whether miRNAs are involved in D. magna maternal effects. We found that miRNAs were differentially expressed in mothers of different ages or nutritional state. We then examined miRNA expression in their eggs, their adult daughters and great granddaughters, which did not experience any treatments. The maternal (treatment) generation exhibited differential expression of miRNAs, as did their eggs, but this was reduced in adult daughters and lost by great granddaughters. Thus, miRNAs are a component of maternal provisioning, but do not appear to be the cause of transgenerational responses under these experimental conditions. MicroRNAs may act in tandem with egg provisioning (e.g., with carbohydrates or fats), and possibly other small RNAs or epigenetic modifications.

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.

Disease spread in age structured populations with maternal age effects.

Fundamental ecological processes, such as extrinsic mortality, determine population age structure. This influences disease spread when individuals of different ages differ in susceptibility or when maternal age determines offspring susceptibility. We show that Daphnia magna offspring born to young mothers are more susceptible than those born to older mothers, and consider this alongside previous observations that susceptibility declines with age in this system. We used a susceptible-infected compartmental model to investigate how age-specific susceptibility and maternal age effects on offspring susceptibility interact with demographic factors affecting disease spread. Our results show a scenario where an increase in extrinsic mortality drives an increase in transmission potential. Thus, we identify a realistic context in which age effects and maternal effects produce conditions favouring disease transmission.

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.

Maternal effects on offspring consumption can stabilize fluctuating predator-prey systems.

Maternal effects, where the conditions experienced by mothers affect the phenotype of their offspring, are widespread in nature and have the potential to influence population dynamics. However, they are very rarely included in models of population dynamics. Here, we investigate a recently discovered maternal effect, where maternal food availability affects the feeding rate of offspring so that well-fed mothers produce fast-feeding offspring. To understand how this maternal effect influences population dynamics, we explore novel predator-prey models where the consumption rate of predators is modified by changes in maternal prey availability. We address the 'paradox of enrichment', a theoretical prediction that nutrient enrichment destabilizes populations, leading to cycling behaviour and an increased risk of extinction, which has proved difficult to confirm in the wild. Our models show that enriched populations can be stabilized by maternal effects on feeding rate, thus presenting an intriguing potential explanation for the general absence of 'paradox of enrichment' behaviour in natural populations. This stabilizing influence should also reduce a population's risk of extinction and vulnerability to harvesting.

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.

Maternal food quantity affects offspring feeding rate in Daphnia magna.

Maternal effects have wide-ranging effects on life-history traits. Here, using the crustacean Daphnia magna, we document a new effect: maternal food quantity affects offspring feeding rate, with low quantities of food triggering mothers to produce slow-feeding offspring. Such a change in the rate of resource acquisition has broad implications for population growth or dynamics and for interactions with, for instance, predators and parasites. This maternal effect can also explain the previously puzzling situation that the offspring of well-fed mothers, despite being smaller, grow and reproduce better than the offspring of food-starved mothers. As an additional source of variation in resource acquisition, this maternal effect may also influence relationships between life-history traits, i.e. trade-offs, and thus constraints on adaptation. Maternal nutrition has long-lasting effects on health and particularly diet-related traits in humans; finding an effect of maternal nutrition on offspring feeding rate in Daphnia highlights the utility of this organism as a powerful experimental model for exploring the relationship between maternal diet and offspring fitness.

Host nutrition alters the variance in parasite transmission potential.

The environmental conditions experienced by hosts are known to affect their mean parasite transmission potential. How different conditions may affect the variance of transmission potential has received less attention, but is an important question for disease management, especially if specific ecological contexts are more likely to foster a few extremely infectious hosts. Using the obligate-killing bacterium Pasteuria ramosa and its crustacean host Daphnia magna, we analysed how host nutrition affected the variance of individual parasite loads, and, therefore, transmission potential. Under low food, individual parasite loads showed similar mean and variance, following a Poisson distribution. By contrast, among well-nourished hosts, parasite loads were right-skewed and overdispersed, following a negative binomial distribution. Abundant food may, therefore, yield individuals causing potentially more transmission than the population average. Measuring both the mean and variance of individual parasite loads in controlled experimental infections may offer a useful way of revealing risk factors for potential highly infectious hosts.

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.

Recent and recurrent selective sweeps of the antiviral RNAi gene Argonaute-2 in three species of Drosophila.

Antagonistic host-parasite interactions can drive rapid adaptive evolution in genes of the immune system, and such arms races may be an important force shaping polymorphism in the genome. The RNA interference pathway gene Argonaute-2 (AGO2) is a key component of antiviral defense in Drosophila, and we have previously shown that genes in this pathway experience unusually high rates of adaptive substitution. Here we study patterns of genetic variation in a 100-kbp region around AGO2 in three different species of Drosophila. Our data suggest that recent independent selective sweeps in AGO2 have reduced genetic variation across a region of more than 50 kbp in Drosophila melanogaster, D. simulans, and D. yakuba, and we estimate that selection has fixed adaptive substitutions in this gene every 30-100 thousand years. The strongest signal of recent selection is evident in D. simulans, where we estimate that the most recent selective sweep involved an allele with a selective advantage of the order of 0.5-1% and occurred roughly 13-60 Kya. To evaluate the potential consequences of the recent substitutions on the structure and function of AGO2, we used fold-recognition and homology-based modeling to derive a structural model for the Drosophila protein, and this suggests that recent substitutions in D. simulans are overrepresented at the protein surface. In summary, our results show that selection by parasites can consistently target the same genes in multiple species, resulting in areas of the genome that have markedly reduced genetic diversity.

The coevolution of virulence: tolerance in perspective.

Coevolutionary interactions, such as those between host and parasite, predator and prey, or plant and pollinator, evolve subject to the genes of both interactors. It is clear, for example, that the evolution of pollination strategies can only be understood with knowledge of both the pollinator and the pollinated. Studies of the evolution of virulence, the reduction in host fitness due to infection, have nonetheless tended to focus on parasite evolution. Host-centric approaches have also been proposed--for example, under the rubric of "tolerance", the ability of hosts to minimize virulence without necessarily minimizing parasite density. Within the tolerance framework, however, there is room for more comprehensive measures of host fitness traits, and for fuller consideration of the consequences of coevolution. For example, the evolution of tolerance can result in changed selection on parasite populations, which should provoke parasite evolution despite the fact that tolerance is not directly antagonistic to parasite fitness. As a result, consideration of the potential for parasite counter-adaptation to host tolerance--whether evolved or medially manipulated--is essential to the emergence of a cohesive theory of biotic partnerships and robust disease control strategies.

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.

Epidemiological, evolutionary, and coevolutionary implications of context-dependent parasitism.

Abstract Victims of infection are expected to suffer increasingly as parasite population growth increases. Yet, under some conditions, faster-growing parasites do not appear to cause more damage, and infections can be quite tolerable. We studied these conditions by assessing how the relationship between parasite population growth and host health is sensitive to environmental variation. In experimental infections of the crustacean Daphnia magna and its bacterial parasite Pasteuria ramosa, we show how easily an interaction can shift from a severe interaction, that is, when host fitness declines substantially with each unit of parasite growth, to a tolerable relationship by changing only simple environmental variables: temperature and food availability. We explored the evolutionary and epidemiological implications of such a shift by modeling pathogen evolution and disease spread under different levels of infection severity and found that environmental shifts that promote tolerance ultimately result in populations harboring more parasitized individuals. We also find that the opportunity for selection, as indicated by the variance around traits, varied considerably with the environmental treatment. Thus, our results suggest two mechanisms that could underlie coevolutionary hotspots and coldspots: spatial variation in tolerance and spatial variation in the opportunity for selection.