An evolutionary epigenetic clock in plants.

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

Deep-sea predator niche segregation revealed by combined cetacean biologging and eDNA analysis of cephalopod prey.

Fundamental insight on predator-prey dynamics in the deep sea is hampered by a lack of combined data on hunting behavior and prey spectra. Deep-sea niche segregation may evolve when predators target specific prey communities, but this hypothesis remains untested. We combined environmental DNA (eDNA) metabarcoding with biologging to assess cephalopod community composition in the deep-sea foraging habitat of two top predator cetaceans. Risso's dolphin and Cuvier's beaked whale selectively targeted distinct epi/meso- and bathypelagic foraging zones, holding eDNA of 39 cephalopod taxa, including 22 known prey. Contrary to expectation, extensive taxonomic overlap in prey spectra between foraging zones indicated that predator niche segregation was not driven by prey community composition alone. Instead, intraspecific prey spectrum differences may drive differentiation for hunting fewer, more calorific, mature cephalopods in deeper waters. The novel combination of methods presented here holds great promise to disclose elusive deep-sea predator-prey systems, aiding in their protection.

Differential gene expression patterns related to lipid metabolism in response to ocean acidification in larvae and juveniles of Atlantic cod.

Elevated environmental carbon dioxide (pCO2) levels have been found to cause organ damage in the early life stages of different commercial fish species, including Atlantic cod (Gadus morhua). To illuminate the underlying mechanisms causing pathologies in the intestines, the kidney, the pancreas and the liver in response to elevated pCO2, we examined related gene expression patterns in Atlantic cod reared for two months under three different pCO2 regimes: 380 µatm (control), 1800 µatm (medium) and 4200 µatm (high). We extracted RNA from whole fish sampled during the larval (32 dph) and early juvenile stage (46 dph) for relative expression analysis of 18 different genes related to essential metabolic pathways. At 32 dph, larvae subjected to the medium treatment displayed an up-regulation of genes mainly associated with fatty acid and glycogen synthesis (GYS2, 6PGL, ACoA, CPTA1, FAS and PPAR1b). Larvae exposed to the high pCO2 treatment upregulated fewer but similar genes (6PGL, ACoA and PPAR1b,). These data suggest stress-induced alterations in the lipid and fatty acid metabolism and a disrupted lipid homeostasis in larvae, providing a mechanistic link to the findings of lipid droplet overload in the liver and organ pathologies. At 46 dph, no significant differences in gene expression were detected, confirming a higher resilience of juveniles in comparison to larvae when exposed to elevated pCO2 up to 4200 µatm.

Transcriptome profiling reveals exposure to predicted end-of-century ocean acidification as a stealth stressor for Atlantic cod larvae.

Ocean acidification (OA), a direct consequence of increasing atmospheric CO2 concentration dissolving in ocean waters, is impacting many fish species. Little is known about the molecular mechanisms underlying the observed physiological impacts in fish. We used RNAseq to characterize the transcriptome of 3 different larval stages of Atlantic cod (Gadus morhua) exposed to simulated OA at levels (1179 µatm CO2) representing end-of-century predictions compared to controls (503 µatm CO2), which were shown to induce tissue damage and elevated mortality in G. morhua. Only few genes were differentially expressed in 6 and 13 days-post-hatching (dph) (3 and 16 genes, respectively), during a period when maximal mortality as a response to elevated pCO2 occurred. At 36 dph, 1413 genes were differentially expressed, most likely caused by developmental asynchrony between the treatment groups, with individuals under OA growing faster. A target gene analysis revealed only few genes of the universal and well-defined cellular stress response to be differentially expressed. We thus suggest that predicted ocean acidification levels constitute a "stealth stress" for early Atlantic cod larvae, with a rapid breakdown of cellular homeostasis leading to organismal death that was missed even with an 8-fold replication implemented in this study.

Effects of parental acclimation and energy limitation in response to high CO2 exposure in Atlantic cod.

Ocean acidification (OA), the dissolution of excess anthropogenic carbon dioxide in ocean waters, is a potential stressor to many marine fish species. Whether species have the potential to acclimate and adapt to changes in the seawater carbonate chemistry is still largely unanswered. Simulation experiments across several generations are challenging for large commercially exploited species because of their long generation times. For Atlantic cod (Gadus morhua), we present first data on the effects of parental acclimation to elevated aquatic CO2 on larval survival, a fundamental parameter determining population recruitment. The parental generation in this study was exposed to either ambient or elevated aquatic CO2 levels simulating end-of-century OA levels (~1100 µatm CO2) for six weeks prior to spawning. Upon fully reciprocal exposure of the F1 generation, we quantified larval survival, combined with two larval feeding regimes in order to investigate the potential effect of energy limitation. We found a significant reduction in larval survival at elevated CO2 that was partly compensated by parental acclimation to the same CO2 exposure. Such compensation was only observed in the treatment with high food availability. This complex 3-way interaction indicates that surplus metabolic resources need to be available to allow a transgenerational alleviation response to ocean acidification.

Phylogeographic differentiation versus transcriptomic adaptation to warm temperatures in Zostera marina, a globally important seagrass.

Populations distributed across a broad thermal cline are instrumental in addressing adaptation to increasing temperatures under global warming. Using a space-for-time substitution design, we tested for parallel adaptation to warm temperatures along two independent thermal clines in Zostera marina, the most widely distributed seagrass in the temperate Northern Hemisphere. A North-South pair of populations was sampled along the European and North American coasts and exposed to a simulated heatwave in a common-garden mesocosm. Transcriptomic responses under control, heat stress and recovery were recorded in 99 RNAseq libraries with ~13 000 uniquely annotated, expressed genes. We corrected for phylogenetic differentiation among populations to discriminate neutral from adaptive differentiation. The two southern populations recovered faster from heat stress and showed parallel transcriptomic differentiation, as compared with northern populations. Among 2389 differentially expressed genes, 21 exceeded neutral expectations and were likely involved in parallel adaptation to warm temperatures. However, the strongest differentiation following phylogenetic correction was between the three Atlantic populations and the Mediterranean population with 128 of 4711 differentially expressed genes exceeding neutral expectations. Although adaptation to warm temperatures is expected to reduce sensitivity to heatwaves, the continued resistance of seagrass to further anthropogenic stresses may be impaired by heat-induced downregulation of genes related to photosynthesis, pathogen defence and stress tolerance.

Bateman's principle and immunity in a sex-role reversed pipefish.

In diverse animal species, from insects to mammals, females display a more efficient immune defence than males. Bateman's principle posits that males maximize their fitness by increasing mating frequency whereas females gain fitness benefits by maximizing their lifespan. As a longer lifespan requires a more efficient immune system, these implications of Bateman's principle may explain widespread immune dimorphism among animals. Because in most extant animals, the provisioning of eggs and a higher parental investment are attributes of the female sex, sex-role reversed species provide a unique opportunity to assess whether or not immune dimorphism depends on life history and not on sex per se. In the broad-nosed pipefish Syngnathus typhle, males brood and nourish the eggs in a ventral pouch and thus invest more into reproduction than females. We found males to have a more active immune response both in field data from four populations and also in an experiment under controlled laboratory conditions. This applied to different measures of immunocompetence using innate as well as adaptive immune system traits. We further determined the specificity of immune response initiation after a fully factorial primary and secondary exposure to a common marine pathogen Vibrio spp. Males not only had a more active but also a more specific immune defence than females. Our results thus indeed suggest that the sex that invests more into the offspring has the stronger immune defence.

Dr. Zompo: an online data repository for Zostera marina and Posidonia oceanica ESTs.

As ecosystem engineers, seagrasses are angiosperms of paramount ecological importance in shallow shoreline habitats around the globe. Furthermore, the ancestors of independent seagrass lineages have secondarily returned into the sea in separate, independent evolutionary events. Thus, understanding the molecular adaptation of this clade not only makes significant contributions to the field of ecology, but also to principles of parallel evolution as well. With the use of Dr. Zompo, the first interactive seagrass sequence database presented here, new insights into the molecular adaptation of marine environments can be inferred. The database is based on a total of 14 597 ESTs obtained from two seagrass species, Zostera marina and Posidonia oceanica, which have been processed, assembled and comprehensively annotated. Dr. Zompo provides experimentalists with a broad foundation to build experiments and consider challenges associated with the investigation of this class of non-domesticated monocotyledon systems. Our database, based on the Ruby on Rails framework, is rich in features including the retrieval of experimentally determined heat-responsive transcripts, mining for molecular markers (SSRs and SNPs), and weighted key word searches that allow access to annotation gathered on several levels including Pfam domains, GeneOntology and KEGG pathways. Well established plant genome sites such as The Arabidopsis Information Resource (TAIR) and the Rice Genome Annotation Project are interfaced by Dr. Zompo. With this project, we have initialized a valuable resource for plant biologists in general and the seagrass community in particular. The database is expected to grow together with more data to come in the near future, particularly with the recent initiation of the Zostera genome sequencing project.The Dr. Zompo database is available at http://drzompo.uni-muenster.de/

Small effective population sizes in two planktonic freshwater copepod species (Eudiaptomus) with apparently large census sizes.

In small planktonic organisms, large census sizes (N(c)) suggest large effective population sizes (N(e)), but reliable estimates are rare. Here, we present N(e)/N(c) ratios for two freshwater copepod species (Eudiaptomus sp.) using temporal samples of multilocus microsatellite genotypes and a pseudo-likelihood approach. N(e)/N(c) ratios were very small in both Eudiaptomus species (10(-7)-10(-8)). Although we hypothesized that the species producing resting eggs (E. graciloides) had a larger N(e) than the other (E. gracilis), estimates were not statistically different (E. graciloides: N(e) = 672.7, CI: 276-1949; E. gracilis: N(e) = 1027.4, CI: 449-2495), suggesting that the propagule bank of E. graciloides had no detectable influence on N(e).

Individual MHC class I and MHC class IIB diversities are associated with male and female reproductive traits in the three-spined stickleback.

Genes of the major histocompatibility complex (MHC) are indispensable for pathogen defence in vertebrates. With wild-caught three-spined sticklebacks (Gasterosteus aculeatus) we conducted the first study to relate individual reproductive parameters to both MHC class I and II diversities. An optimal MHC class IIB diversity was found for male nest quality. However, male breeding colouration was most intense at a maximal MHC class I diversity. One MHC class I allele was associated with male redness. Similarly, one MHC class IIB allele was associated with continuous rather than early female reproduction, possibly extending the reproductive period. Both alleles occurred more frequently with increasing individual allele diversity. We suggest that if an allele is currently not part of the optimum, it had not been propagated by choosy females. The parasite against which this allele provides resistance is therefore unlikely to have been predominant the previous year - a step to negative frequency-dependent selection.

Genetic variation in MHC class II expression and interactions with MHC sequence polymorphism in three-spined sticklebacks.

Genes of the major histocompatibility complex (MHC) have been studied for several decades because of their pronounced allelic polymorphism. Structural allelic polymorphism is, however, not the only source of variability subjected to natural selection. Genetic variation may also exist in gene expression patterns. Here, we show that in a natural population of three-spined sticklebacks (Gasterosteus aculeatus) the expression of MHC class IIB genes was positively correlated with parasite load, which indicates increased immune activation of the MHC when infections are frequent. To experimentally study MHC expression, we used laboratory-bred sticklebacks that were exposed to three naturally occurring species of parasite. We found strong differences in MHC class IIB expression patterns among fish families, which were consistent over two generations, thus demonstrating a genetic component. The average number of MHC class IIB sequence variants within families was negatively correlated to the MHC expression level suggesting compensatory up-regulation in fish with a low (i.e. suboptimal) MHC sequence variability. The observed differences among families and the negative correlation with individual sequence diversity imply that MHC expression is evolutionary relevant for the onset and control of the immune response in natural populations.

How a complex life cycle can improve a parasite's sex life.

How complex life cycles of parasites are maintained is still a fascinating and unresolved topic. Complex life cycles using three host species, free-living stages, asexual and sexual reproduction are widespread in parasitic helminths. For such life cycles, we propose here that maintaining a second intermediate host in the life cycle can be advantageous for the individual parasite to increase the intermixture of different clones and therefore decrease the risk of matings between genetically identical individuals in the definitive host. Using microsatellite markers, we show that clone mixing occurs from the first to the second intermediate host in natural populations of the eye-fluke Diplostomum pseudospathaceum. Most individuals released by the first intermediate host belonged to one clone. In contrast, the second intermediate host was infected with a diverse array of mostly unique parasite genotypes. The proposed advantage of increased parasite clone intermixture may be a novel selection pressure favouring the maintenance of complex life cycles.

Local genetic structure in a clonal dioecious angiosperm.

We used seven microsatellite loci to characterize genetic structure and clonal architecture at three different spatial scales (from meters to centimetres) of a Cymodocea nodosa population. C. nodosa exhibits both sexual reproduction and vegetative propagation by rhizome elongation. Seeds remain buried in the sediment nearby the mother plant in a dormant stage until germination. Seed dispersal potential is therefore expected to be extremely restricted. High clonal diversity (up to 67% of distinct genotypes) and a highly intermingled configuration of genets at different spatial scales were found. No significant differences in genetic structure were found among the three spatial scales, indicating that genetic diversity is evenly distributed along the meadow. Autocorrelation analyses of kinship estimates confirmed the absence of spatial clumping of genets at small spatial scale and the expectations of a very restricted seed dispersal (observed dispersal range 1-21 m) in this species.

Evolutionary divergence and possible incipient speciation in post-glacial populations of a cosmopolitan aquatic plant.

Habitat configuration is expected to have a major influence on genetic exchange and evolutionary divergence among populations. Aquatic organisms occur in two fundamentally different habitat types, the sea and freshwater lakes, making them excellent models to study the contrasting effects of continuity vs. isolation on genetic divergence. We compared the divergence in post-glacial populations of a cosmopolitan aquatic plant, the pondweed Potamogeton pectinatus that simultaneously occurs in freshwater lakes and coastal marine sites. Relative levels of gene flow were inferred in 12 lake and 14 Baltic Sea populations in northern Germany using nine highly polymorphic microsatellite markers developed for P. pectinatus. We found highly significant isolation-by-distance in both habitat types (P < 0.001). Genetic differentiation increased approximately 2.5-times faster among freshwater populations compared with those from the Baltic Sea. As different levels of genetic drift or population history cannot explain these differences, higher population connectivity in the sea relative to freshwater populations is the most likely source of contrasting evolutionary divergence. These findings are consistent with the notion that freshwater angiosperms are more conducive to allopatric speciation than their life-history counterparts in the sea, the relative species poor seagrasses. Surprisingly, population pairs from different habitat types revealed almost maximal genetic divergence expected for complete reproductive isolation, regardless of their respective geographical distance. Hence, the barrier to gene flow between lake and sea habitat types cannot be due to dispersal limitation. We may thus have identified a case of rapid incipient speciation in post-glacial populations of a widespread aquatic plant.

Flexible mating: cross-pollination affects sex-expression in a marine clonal plant.

Many functionally hermaphroditic plants have evolved mechanisms to reduce interference between the sex functions and to optimize reproductive output. In addition to physical mechanisms such as the spatial (herkogamy) and temporal (dichogamy) separation of male and female functions, plasticity in sex expression by means of mate-recognition (flexible mating) could be important in plants with variable access to cross-pollen. This applies particularly to clonal plants because of their modular growth form. We experimentally tested for the effects of pollen source and vegetative neighbourhood on instantaneous sex ratio and seed production in the self-compatible clonal marine angiosperm Zostera marina L. To this end, we exposed the (monoecious) flowering shoots to self and cross-pollen and to neighbourhoods of their own and a mix of foreign vegetative shoots. Flowering shoots that had been exposed to cross-pollen showed (1) a significantly lower female/male ratio at peak flowering, evidence for mate-recognition, and (2) a significantly higher seed set by the end of the season. Both effects were independent of the genetic composition of their vegetative neighbourhood. The results suggest that Z. marina maintains a cryptic self-incompatibility system not previously described for angiosperms with sub-aqueous pollination. In Z. marina, and possibly other self-compatible clonal plant species, mate-recognition could be a means of increasing the out-crossing probability for flowering shoots with central positions within their clone.

Multiple parasites are driving major histocompatibility complex polymorphism in the wild.

Parasite mediated selection may result in arms races between host defence and parasite virulence. In particular, simultaneous infections from multiple parasite species should cause diversification (i.e. balancing selection) in resistance genes both at the population and the individual level. Here, we tested these ideas in highly polymorphic major histocompatibility complex (MHC) genes from three-spined sticklebacks (Gasterosteus aculeatus L.). In eight natural populations, parasite diversity (15 different species), and MHC class IIB diversity varied strongly between habitat types (lakes vs. rivers vs. estuaries) with lowest values in rivers. Partial correlation analysis revealed an influence of parasite diversity on MHC class IIB variation whereas general genetic diversity assessed at seven microsatellite loci was not significantly correlated with parasite diversity. Within individual fish, intermediate, rather than maximal allele numbers were associated with minimal parasite load, supporting theoretical models of self-reactive T-cell elimination. The optimal individual diversity matched those values female fish try to achieve in their offspring by mate choice. We thus present correlative evidence supporting the 'allele counting' strategy for optimizing the immunocompetence in stickleback offspring.

Genetic neighbourhood of clone structures in eelgrass meadows quantified by spatial autocorrelation of microsatellite markers.

Limited dispersal distances in plant populations frequently cause local genetic structure, which can be quantified by spatial autocorrelation. In clonal plants, three levels of spatial organization can contribute to positive autocorrelation; namely, the neighbourhood of (a) ramets, (b) clone fragments and (c) entire clones. Here we use data from an exhaustive sampling scheme on a clonal plant to measure the contribution of the neighbourhoods of each distinct clonal structure to total spatial autocorrelation. Four plots (256 grid points each) within dense meadows of the marine clonal plant Zostera marina (eelgrass) were sampled for clone structure with nine microsatellite markers ( approximately 80 alleles). We found significant coancestry (f(ij)), at all three levels of spatial organization, even when not allowing for joins between samples of identical genets. In addition, absolute values of f(ij) and the maximum distance with significant positive f(ij) decreased with the progressive exclusion of joins between alike genotypes. The neighbourhood of this clonal plant thus consists of three levels of organization, which are reflected in different kinship structures. Each of these kinship structures may affect the level of biparental inbreeding and the physical distance between flowering shoots and their outcrossing neighbourhood. These results also emphasize the notion that spatial autocorrelation crucially depends on the scale and intensity of sampling.

Inbreeding depression influences genet size distribution in a marine angiosperm.

Although inbreeding depression is a major genetic phenomena influencing individual fitness, it is difficult to measure in wild populations. An alternative approach is to correlate heterozygosity, measured using highly polymorphic markers, with a fitness-correlated trait. In clonal plants, genet size is predicted to be fitness correlated. Here we test the prediction that the genet size distribution of the marine clonal plant Zostera marina (eelgrass) is influenced by inbreeding depression. We used nine polymorphic microsatellite markers to access the fine scale clonal structure and to measure individual heterozygosity within 4 plots (each corresponds to 256 m2, sampled at 1-m intervals) in two populations along the German Baltic Coast. The same plots were also sampled for flowering and vegetative shoots to obtain estimates for sexual reproductive output at the level of the genetic individual. We found substantial differences in the genet size distribution between the two populations that may be explained by different disturbance frequency. In both populations, clone size was significantly positively correlated with the total number of flowering shoots, indicating that larger clones have a higher reproductive output. Individual heterozygosity was significantly positively associated with clone size. The effect was much stronger in Falkenstein (low disturbance) than in Maasholm (high disturbance). The results indicate that in a low disturbance population the relatively outbred clones occupy a higher proportion of the available space, possibly because they outcompete relatively inbred neighbours.

Habitat differentiation vs. isolation-by-distance: the genetic population structure of Elymus athericus in European salt marshes.

We investigated genetic differentiation among populations of the clonal grass Elymus athericus, a common salt-marsh species occurring along the Wadden Sea coast of Europe. While E. athericus traditionally occurs in the high salt marsh, it recently also invaded lower parts of the marsh. In one of the first analyses of the genetic population structure in salt-marsh species, we were interested in population differentiation through isolation-by-distance, and among strongly divergent habitats (low and high marsh) in this wind- and water-dispersed species. High and low marsh habitats were sampled at six sites throughout the Wadden Sea. Based on reciprocal transplantation experiments conducted earlier revealing lower survival of foreign genotypes we predicted reduced gene flow among habitats. Accordingly, an analysis with polymorphic cross-species microsatellite primers revealed significant genetic differentiation between high and low marsh habitats already on a very small scale (< 100 m), while isolation-by-distance was present only on larger scales (60-443 km). In an analysis of molecular variance we found that 14% of the genetic variance could be explained by the differentiation between habitats, as compared to only 8.9% to geographical (isolation-by-distance) effects among six sites 2.5-443 km distant from each other. This suggests that markedly different selection regimes between these habitats, in particular intraspecific competition and herbivory, result in habitat adaptation and restricted gene flow over distances as small as 80 m. Hence, the genetic population structure of plant species can only be understood when considering geographical and selection-mediated restrictions to gene flow simultaneously.

Rapid genetic divergence in postglacial populations of threespine stickleback (Gasterosteus aculeatus): the role of habitat type, drainage and geographical proximity.

The goal of this study was to evaluate the role of common ancestry, and of geographical or reproductive isolation, in genetic divergence in populations of threespine sticklebacks (Gasterosteus aculeatus). Using seven DNA microsatellite loci we compared the effects of habitat type, drainage system and geographical proximity on genetic distance among 16 populations situated in an area in Schleswig-Holstein (Germany) that became deglaciated approximately 12 000 years ago. Stickleback population structure correlated only weakly with drainage system, whereas the primary divergence was among habitat types. Phylogenetic analysis revealed that lake (n = 7) and river (n = 5) populations formed two distinct clades (Cavalli-Sforza's and Edwards' chord distance, 82-100% bootstrap support) at approximately equal genetic distances to a third clade, comprising putative estuarine (n = 4) ancestors. Allele frequencies in lake and river populations represented different subsets of the genetically more diverse estuarine populations. In nested amovas approximately twice the genetic variance was distributed among lake vs. river vs. estuarine populations as compared with the combined effects of drainage system and geographical distance. Limited gene flow between habitat types must have been established after postglacial colonization, suggesting ecological hybrid inferiority or behavioural mating barriers between ecotypes. Within estuarine and lake populations, population differentiation followed an isolation-by-distance model. Given the high observed heterozygosities within the 16 study populations (HO = 0.65-0.87), the mean divergence between lake and river population pairs (FST = 0.18 +/- 0.007) would be reached after 300-6000 generations in a stepwise mutation model, depending on the size of N(e). This demonstrates both the utility of hypervariable microsatellites for detecting recent population divergences and the danger of operating at temporal or spatial scales which are beyond their resolution.