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.

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.

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.

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.

DNA methylation differs extensively between strains of the same geographical origin and changes with age in Daphnia magna.

BACKGROUND: Patterns of methylation influence lifespan, but methylation and lifespan may also depend on diet, or differ between genotypes. Prior to this study, interactions between diet and genotype have not been explored together to determine their influence on methylation. The invertebrate Daphnia magna is an excellent choice for testing the epigenetic response to the environment: parthenogenetic offspring are identical to their siblings (making for powerful genetic comparisons), they are relatively short lived and have well-characterised inter-strain life-history trait differences. We performed a survival analysis in response to caloric restriction and then undertook a 47-replicate experiment testing the DNA methylation response to ageing and caloric restriction of two strains of D. magna. RESULTS: Methylated cytosines (CpGs) were most prevalent in exons two to five of gene bodies. One strain exhibited a significantly increased lifespan in response to caloric restriction, but there was no effect of food-level CpG methylation status. Inter-strain differences dominated the methylation experiment with over 15,000 differently methylated CpGs. One gene, Me31b, was hypermethylated extensively in one strain and is a key regulator of embryonic expression. Sixty-one CpGs were differentially methylated between young and old individuals, including multiple CpGs within the histone H3 gene, which were hypermethylated in old individuals. Across all age-related CpGs, we identified a set that are highly correlated with chronological age. CONCLUSIONS: Methylated cytosines are concentrated in early exons of gene sequences indicative of a directed, non-random, process despite the low overall DNA methylation percentage in this species. We identify no effect of caloric restriction on DNA methylation, contrary to our previous results, and established impacts of caloric restriction on phenotype and gene expression. We propose our approach here is more robust in invertebrates given genome-wide CpG distributions. For both strain and ageing, a single gene emerges as differentially methylated that for each factor could have widespread phenotypic effects. Our data showed the potential for an epigenetic clock at a subset of age positions, which is exciting but requires confirmation.

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.

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.

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.

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.

Bigger is better: changes in body size explain a maternal effect of food on offspring disease resistance.

Maternal effects triggered by changes in the environment (e.g., nutrition or crowding) can influence the outcome of offspring-parasite interactions, with fitness consequences for the host and parasite. Outside of the classic example of antibody transfer in vertebrates, proximate mechanisms have been little studied, and thus, the adaptive significance of maternal effects on infection is not well resolved. We sought to determine why food-stressed mothers give birth to offspring that show a low rate of infection when the crustacean Daphnia magna is exposed to an orally infective bacterial pathogen. These more-resistant offspring are also larger at birth and feed at a lower rate. Thus, reduced disease resistance could result from slow-feeding offspring ingesting fewer bacterial spores or because their larger size allows for greater immune investment. To distinguish between these theories, we performed an experiment in which we measured body size, feeding rate, and susceptibility, and were able to show that body size is the primary mechanism causing altered susceptibility: Larger Daphnia were less likely to become infected. Contrary to our predictions, there was also a trend that fast-feeding Daphnia were less likely to become infected. Thus, our results explain how a maternal environmental effect can alter offspring disease resistance (though body size), and highlight the potential complexity of relationship between feeding rate and susceptibility in a host that encounters a parasite whilst feeding.

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 transcriptome by RNA-Seq across 12 environmental stressors.

The full exploration of gene-environment interactions requires model organisms with well-characterized ecological interactions in their natural environment, manipulability in the laboratory and genomic tools. The waterflea Daphnia magna is an established ecological and toxicological model species, central to the food webs of freshwater lentic habitats and sentinel for water quality. Its tractability and cyclic parthenogenetic life-cycle are ideal to investigate links between genes and the environment. Capitalizing on this unique model system, the STRESSFLEA consortium generated a comprehensive RNA-Seq data set by exposing two inbred genotypes of D. magna and a recombinant cross of these genotypes to a range of environmental perturbations. Gene models were constructed from the transcriptome data and mapped onto the draft genome of D. magna using EvidentialGene. The transcriptome data generated here, together with the available draft genome sequence of D. magna and a high-density genetic map will be a key asset for future investigations in environmental genomics.

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.

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.

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.

Effects of juvenile host density and food availability on adult immune response, parasite resistance and virulence in a Daphnia-parasite system.

Host density can increase infection rates and reduce host fitness as increasing population density enhances the risk of becoming infected either through increased encounter rate or because host condition may decline. Conceivably, potential hosts could take high host density as a cue to up-regulate their defence systems. However, as host density usually covaries with food availability, it is difficult to examine the importance of host density in isolation. Thus, we performed two full-factorial experiments that varied juvenile densities of Daphnia magna (a freshwater crustacean) and food availability independently. We also included a simulated high-density treatment, where juvenile experimental animals were kept in filtered media that previously maintained Daphnia at high-density. Upon reaching adulthood, we exposed the Daphnia to their sterilizing bacterial parasite, Pasteuria ramosa, and examined how the juvenile treatments influenced the likelihood and severity of infection (Experiment I) and host immune investment (Experiment II). Neither juvenile density nor food treatments affected the likelihood of infection; however, well-fed hosts that were well-fed as juveniles produced more offspring prior to sterilization than their less well-fed counterparts. By contrast, parasite growth was independent of host juvenile resources or host density. Parasite-exposed hosts had a greater number of circulating haemocytes than controls (i.e., there was a cellular immune response), but the magnitude of immune response was not mediated by food availability or host density. These results suggest that density dependent effects on disease arise primarily through correlated changes in food availability: low food could limit parasitism and potentially curtail epidemics by reducing both the host's and parasite's reproduction as both depend on the same food.

The bacterial parasite Pasteuria ramosa is not killed if it fails to infect: implications for coevolution.

Strong selection on parasites, as well as on hosts, is crucial for fueling coevolutionary dynamics. Selection will be especially strong if parasites that encounter resistant hosts are destroyed and diluted from the local environment. We tested whether spores of the bacterial parasite Pasteuria ramosa were passed through the gut (the route of infection) of their host, Daphnia magna, and whether passaged spores remained viable for a "second chance" at infecting a new host. In particular, we tested if this viability (estimated via infectivity) depended on host genotype, whether or not the genotype was susceptible, and on initial parasite dose. Our results show that Pasteuria spores generally remain viable after passage through both susceptible and resistant Daphnia. Furthermore, these spores remained infectious even after being frozen for several weeks. If parasites can get a second chance at infecting hosts in the wild, selection for infection success in the first instance will be reduced. This could also weaken reciprocal selection on hosts and slow the coevolutionary process.