Demographic Basis of Spatially Structured Fluctuations in a Threespine Stickleback Metapopulation.

AbstractUncovering the demographic basis of population fluctuations is a central goal of population biology. This is particularly challenging for spatially structured populations, which require disentangling synchrony in demographic rates from coupling via movement between locations. In this study, we fit a stage-structured metapopulation model to a 29-year time series of threespine stickleback abundance in the heterogeneous and productive Lake Mývatn, Iceland. The lake comprises two basins (North and South) connected by a channel through which the stickleback disperse. The model includes time-varying demographic rates, allowing us to assess the potential contributions of recruitment and survival, spatial coupling via movement, and demographic transience to the population's large fluctuations in abundance. Our analyses indicate that recruitment was only modestly synchronized between the two basins, whereas survival probabilities of adults were more strongly synchronized, contributing to cyclic fluctuations in the lake-wide population size with a period of approximately 6 years. The analyses further show that the two basins were coupled through movement, with the North Basin subsidizing the South Basin and playing a dominant role in driving the lake-wide dynamics. Our results show that cyclic fluctuations of a metapopulation can be explained in terms of the combined effects of synchronized demographic rates and spatial coupling.

Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback.

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole-genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h2  = 0.42-0.65), indicating adaptive potential for all traits. Environment-association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation - supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.

Spatial variation in the evolutionary potential and constraints of basal metabolic rate and body mass in a wild bird.

An organism's energy budget is strongly related to resource consumption, performance, and fitness. Hence, understanding the evolution of key energetic traits, such as basal metabolic rate (BMR), in natural populations is central for understanding life-history evolution and ecological processes. Here we used quantitative genetic analyses to study evolutionary potential of BMR in two insular populations of the house sparrow (Passer domesticus). We obtained measurements of BMR and body mass (Mb ) from 911 house sparrows on the islands of Leka and Vega along the coast of Norway. These two populations were the source populations for translocations to create an additional third, admixed 'common garden' population in 2012. With the use of a novel genetic group animal model concomitant with a genetically determined pedigree, we differentiate genetic and environmental sources of variation, thereby providing insight into the effects of spatial population structure on evolutionary potential. We found that the evolutionary potential of BMR was similar in the two source populations, whereas the Vega population had a somewhat higher evolutionary potential of Mb than the Leka population. BMR was genetically correlated with Mb in both populations, and the conditional evolutionary potential of BMR (independent of body mass) was 41% (Leka) and 53% (Vega) lower than unconditional estimates. Overall, our results show that there is potential for BMR to evolve independently of Mb , but that selection on BMR and/or Mb may have different evolutionary consequences in different populations of the same species.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus).

BACKGROUND: Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression. RESULTS: Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring. CONCLUSIONS: This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

Stream microbial communities and ecosystem functioning show complex responses to multiple stressors in wastewater.

Multiple anthropogenic drivers are changing ecosystems globally, with a disproportionate and intensifying impact on freshwater habitats. A major impact of urbanization are inputs from wastewater treatment plants (WWTPs). Initially designed to reduce eutrophication and improve water quality, WWTPs increasingly release a multitude of micropollutants (MPs; i.e., synthetic chemicals) and microbes (including antibiotic-resistant bacteria) to receiving environments. This pollution may have pervasive impacts on biodiversity and ecosystem services. Viewed through multiple lenses of macroecological and ecotoxicological theory, we combined field, flume, and laboratory experiments to determine the effects of wastewater (WW) on microbial communities and organic-matter processing using a standardized decomposition assay. First, we conducted a mensurative experiment sampling 60 locations above and below WWTP discharges in 20 Swiss streams. Microbial respiration and decomposition rates were positively influenced by WW inputs via warming and nutrient enrichment, but with a notable exception: WW decreased the activation energy of decomposition, indicating a "slowing" of this fundamental ecosystem process in response to temperature. Second, next-generation sequencing indicated that microbial community structure below WWTPs was altered, with significant compositional turnover, reduced richness, and evidence of negative MP influences. Third, a series of flume experiments confirmed that although diluted WW generally has positive influences on microbial-mediated processes, the negative effects of MPs are "masked" by nutrient enrichment. Finally, transplant experiments suggested that WW-borne microbes enhance decomposition rates. Taken together, our results affirm the multiple stressor paradigm by showing that different aspects of WW (warming, nutrients, microbes, and MPs) jointly influence ecosystem functioning in complex ways. Increased respiration rates below WWTPs potentially generate ecosystem "disservices" via greater carbon evasion from streams and rivers. However, toxic MP effects may fundamentally alter ecological scaling relationships, indicating the need for a rapprochement between ecotoxicological and macroecological perspectives.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development.

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life-stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non-linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life-stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

Cryptic lineages of a common alpine mayfly show strong life-history divergence.

Understanding ecological divergence of morphologically similar but genetically distinct species - previously considered as a single morphospecies - is of key importance in evolutionary ecology and conservation biology. Despite their morphological similarity, cryptic species may have evolved distinct adaptations. If such ecological divergence is unaccounted for, any predictions about their responses to environmental change and biodiversity loss may be biased. We used spatio-temporally replicated field surveys of larval cohort structure and population genetic analyses (using nuclear microsatellite markers) to test for life-history divergence between two cryptic lineages of the alpine mayfly Baetis alpinus in the Swiss Alps. We found that the more widespread and abundant cryptic lineage represents a 'generalist' with at least two cohorts per year, whereas the less abundant lineage is restricted to higher elevations and represents a 'specialist' with a single cohort per year. Importantly, our results indicate partial temporal segregation in reproductive periods between these lineages, potentially facilitating local coexistence and reproductive isolation. Taken together, our findings emphasize the need for a taxonomic revision: widespread and apparently generalist morphospecies can hide cryptic lineages with much narrower ecological niches and distribution ranges.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps.

BACKGROUND: Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation. RESULTS: We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages. CONCLUSIONS: Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

Molecular phenotyping of maternally mediated parallel adaptive divergence within Rana arvalis and Rana temporaria.

When similar selection acts on the same traits in multiple species or populations, parallel evolution can result in similar phenotypic changes, yet the underlying molecular architecture of parallel phenotypic divergence can be variable. Maternal effects can influence evolution at ecological timescales and facilitate local adaptation, but their contribution to parallel adaptive divergence is unclear. In this study, we (i) tested for variation in embryonic acid tolerance in a common garden experiment and (ii) used molecular phenotyping of egg coats to investigate the molecular basis of maternally mediated parallel adaptive divergence in two amphibian species (Rana arvalis and Rana temporaria). Our results on three R. arvalis and two R. temporaria populations show that adaptive divergence in embryonic acid tolerance is mediated via maternally derived egg coats in both species. We find extensive polymorphism in egg jelly coat glycoproteins within both species and that acid-tolerant clutches have more negatively charged egg jelly - indicating that the glycosylation status of the jelly coat proteins is under divergent selection in acidified environments, likely due to its impact on jelly water balance. Overall, these data provide evidence for parallel mechanisms of adaptive divergence in two species. Our study highlights the importance of studying intraspecific molecular variation in egg coats and, specifically, their glycoproteins, to increase understanding of underlying forces maintaining variation in jelly coats.

Evolution of egg coats: linking molecular biology and ecology.

One central goal of evolutionary biology is to explain how biological diversity emerges and is maintained in nature. Given the complexity of the phenotype and the multifaceted nature of inheritance, modern evolutionary ecological studies rely heavily on the use of molecular tools. Here, we show how molecular tools help to gain insight into the role of egg coats (i.e. the extracellular structures surrounding eggs and embryos) in evolutionary diversification. Egg coats are maternally derived structures that have many biological functions from mediating fertilization to protecting the embryo from environmental hazards. They show great molecular, structural and functional diversity across species, but intraspecific variability and the role of ecology in egg coat evolution have largely been overlooked. Given that much of the variation that influences egg coat function is ultimately determined by their molecular phenotype, cutting-edge molecular tools (e.g. proteomics, glycomics and transcriptomics), combined with functional assays, are needed for rigorous inferences on their evolutionary ecology. Here, we identify key research areas and highlight emerging molecular techniques that can increase our understanding of the role of egg coats in the evolution of biological diversity, from adaptation to speciation.

Mechanistic basis of adaptive maternal effects: egg jelly water balance mediates embryonic adaptation to acidity in Rana arvalis.

Environmental stress, such as acidification, can challenge persistence of natural populations and act as a powerful evolutionary force at ecological time scales. The ecological and evolutionary responses of natural populations to environmental stress at early life-stages are often mediated via maternal effects. During early life-stages, maternal effects commonly arise from egg coats (the extracellular structures surrounding the embryo), but the role of egg coats has rarely been studied in the context of adaptation to environmental stress. Previous studies on the moor frog Rana arvalis found that the egg coat mediated adaptive divergence along an acidification gradient in embryonic acid stress tolerance. However, the exact mechanisms underlying these adaptive maternal effects remain unknown. Here, we investigated the role of water balance and charge state (zeta potential) of egg jelly coats in embryonic adaptation to acid stress in three populations of R. arvalis. We found that acidic pH causes severe water loss in the egg jelly coat, but that jelly coats from an acid-adapted population retained more water than jelly coats from populations not adapted to acidity. Moreover, embryonic acid tolerance (survival at pH 4.0) correlated with both water loss and charge state of the jelly, indicating that negatively charged glycans influence jelly water balance and contribute to embryonic adaptation to acidity. These results indicate that egg coats can harbor extensive intra-specific variation, probably facilitated in part via strong selection on water balance and glycosylation status of egg jelly coats. These findings shed light on the molecular mechanisms of environmental stress tolerance and adaptive maternal effects.

Multifarious selection through environmental change: acidity and predator-mediated adaptive divergence in the moor frog (Rana arvalis).

Environmental change can simultaneously cause abiotic stress and alter biological communities, yet adaptation of natural populations to co-changing environmental factors is poorly understood. We studied adaptation to acid and predator stress in six moor frog (Rana arvalis) populations along an acidification gradient, where abundance of invertebrate predators increases with increasing acidity of R. arvalis breeding ponds. First, we quantified divergence among the populations in anti-predator traits (behaviour and morphology) at different rearing conditions in the laboratory (factorial combinations of acid or neutral pH and the presence or the absence of a caged predator). Second, we evaluated relative fitness (survival) of the populations by exposing tadpoles from the different rearing conditions to predation by free-ranging dragonfly larvae. We found that morphological defences (relative tail depth) as well as survival of tadpoles under predation increased with increasing pond acidity (under most experimental conditions). Tail depth and larval size mediated survival differences among populations, but the contribution of trait divergence to survival was strongly dependent on prior rearing conditions. Our results indicate that R. arvalis populations are adapted to the elevated predator pressure in acidified ponds and emphasize the importance of multifarious selection via both direct (here: pH) and indirect (here: predators) environmental changes.

Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus.

BACKGROUND: A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure. RESULTS: These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect. CONCLUSIONS: Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Effects of anthropogenic stress on hosts and their microbiomes: Treated wastewater alters performance and gut microbiome of a key detritivore (Asellus aquaticus).

Human activity is a major driver of ecological and evolutionary change in wild populations and can have diverse effects on eukaryotic organisms as well as on environmental and host-associated microbial communities. Although host-microbiome interactions can be a major determinant of host fitness, few studies consider the joint responses of hosts and their microbiomes to anthropogenic changes. In freshwater ecosystems, wastewater is a widespread anthropogenic stressor that represents a multifarious environmental perturbation. Here, we experimentally tested the impact of treated wastewater on a keystone host (the freshwater isopod Asellus aquaticus) and its gut microbiome. We used a semi-natural flume experiment, in combination with 16S rRNA amplicon sequencing, to assess how different concentrations (0%, 30%, and 80%) of nonfiltered wastewater (i.e. with chemical toxicants, nutrients, organic particles, and microbes) versus ultrafiltered wastewater (i.e. only dissolved pollutants and nutrients) affected host survival, growth, and food consumption as well as mid- and hindgut bacterial community composition and diversity. Our results show that while host survival was not affected by the treatments, host growth increased and host feeding rate decreased with nonfiltered wastewater - potentially indicating that A. aquaticus fed on organic matter and microbes available in nonfiltered wastewater. Furthermore, even though the midgut microbiome (diversity and composition) was not affected by any of our treatments, nonfiltered wastewater influenced bacterial composition (but not diversity) in the hindgut. Ultrafiltered wastewater, on the other hand, affected both community composition and bacterial diversity in the hindgut, an effect that in our system differed between sexes. While the functional consequences of microbiome changes and their sex specificity are yet to be tested, our results indicate that different components of multifactorial stressors (i.e. different constituents of wastewater) can affect hosts and their microbiome in distinct (even opposing) manners and have a substantial impact on eco-evolutionary responses to anthropogenic stressors.

Context dependent variation in corticosterone and phenotypic divergence of Rana arvalis populations along an acidification gradient.

BACKGROUND: Physiological processes, as immediate responses to the environment, are important mechanisms of phenotypic plasticity and can influence evolution at ecological time scales. In stressful environments, physiological stress responses of individuals are initiated and integrated via the release of hormones, such as corticosterone (CORT). In vertebrates, CORT influences energy metabolism and resource allocation to multiple fitness traits (e.g. growth and morphology) and can be an important mediator of rapid adaptation to environmental stress, such as acidification. The moor frog, Rana arvalis, shows adaptive divergence in larval life-histories and predator defense traits along an acidification gradient in Sweden. Here we take a first step to understanding the role of CORT in this adaptive divergence. We conducted a fully factorial laboratory experiment and reared tadpoles from three populations (one acidic, one neutral and one intermediate pH origin) in two pH treatments (Acid versus Neutral pH) from hatching to metamorphosis. We tested how the populations differ in tadpole CORT profiles and how CORT is associated with tadpole life-history and morphological traits. RESULTS: We found clear differences among the populations in CORT profiles across different developmental stages, but only weak effects of pH treatment on CORT. Tadpoles from the acid origin population had, on average, lower CORT levels than tadpoles from the neutral origin population, and the intermediate pH origin population had intermediate CORT levels. Overall, tadpoles with higher CORT levels developed faster and had shorter and shallower tails, as well as shallower tail muscles. CONCLUSIONS: Our common garden results indicate among population divergence in CORT levels, likely reflecting acidification mediated divergent selection on tadpole physiology, concomitant to selection on larval life-histories and morphology. However, CORT levels were highly environmental context dependent. Jointly these results indicate a potential role for CORT as a mediator of multi-trait divergence along environmental stress gradients in natural populations. At the same time, the population level differences and high context dependency in CORT levels suggest that snapshot assessment of CORT in nature may not be reliable bioindicators of stress.

Short-term responses of Rana arvalis tadpoles to pH and predator stress: adaptive divergence in behavioural and physiological plasticity?

Environmental stress is a major driver of ecological and evolutionary processes in nature. To cope with stress, organisms can adjust through phenotypic plasticity and/or adapt through genetic change. Here, we compared short-term behavioural (activity) and physiological (corticosterone levels, CORT) responses of Rana arvalis tadpoles from two divergent populations (acid origin, AOP, versus neutral origin, NOP) to acid and predator stress. Tadpoles were initially reared in benign conditions at pH 7 and then exposed to a combination of two pH (acid versus neutral) and two predator cue (predator cue versus no predator cue) treatments. We assessed behavioural activity within the first 15 min, and tissue CORT within 8 and 24 h of stress exposure. Both AOP and NOP tadpoles reduced their activity in acidic pH, but the response to the predator cue differed between the populations: AOP tadpoles increased whereas NOP tadpoles decreased their activity. The AOP and NOP tadpoles differed also in their CORT responses, with AOP being more responsive (CORT levels of NOP tadpoles did not differ statistically across treatments). After 8 h exposure, AOP tadpoles had elevated CORT levels in the acid-predator cue treatment and after 24 h exposure they had elevated CORT levels in all three stress treatments (relative to the benign neutral-no-cue treatment). These results suggest that adaptation to environmental acidification in R. arvalis is mediated, in part, via behavioural and hormonal plasticity.

Wastewater constituents impact biofilm microbial community in receiving streams.

Microbial life in natural biofilms is dominated by prokaryotes and microscopic eukaryotes living in dense association. In stream ecosystems, microbial biofilms influence primary production, elemental cycles, food web interactions as well as water quality. Understanding how biofilm communities respond to anthropogenic impacts, such as wastewater treatment plant (WWTP) effluent, is important given the key role of biofilms in stream ecosystem function. Here, we implemented 16S and 18S rRNA gene sequencing of stream biofilms upstream (US) and downstream (DS) of WWTP effluents in four Swiss streams to test how bacterial and eukaryotic communities respond to wastewater constituents. Stream biofilm composition was strongly affected by geographic location - particularly for bacteria. However, the abundance of certain microbial community members was related to micropollutants in the wastewater - among bacteria, micropollutant-associated members were found e.g. in Alphaproteobacteria, and among eukaryotes e.g. in Bacillariophyta (algal diatoms). This study corroborates several previously characterized responses (e.g. as seen in diatoms), but also reveals previously unknown community responses - such as seen in Alphaproteobacteria. This study advances our understanding of the ecological impact of the current wastewater treatment practices and provides information about potential new marker organisms to assess ecological change in stream biofilms.

Variation in egg size and offspring phenotype among and within seven Arctic charr morphs.

Maternal effects have the potential to alter early developmental processes of offspring and contribute to adaptive diversification. Egg size is a major contributor to offspring phenotype, which can influence developmental trajectories and potential resource use. However, to what extent intraspecific variation in egg size facilitates evolution of resource polymorphism is poorly understood. We studied multiple resource morphs of Icelandic Arctic charr, ranging from an anadromous morph-with a phenotype similar to the proposed ancestral phenotype-to sympatric morphs that vary in their degree of phenotypic divergence from the ancestral anadromous morph. We characterized variation in egg size and tested whether egg size influenced offspring phenotype at early life stages (i.e., timing of- and size at- hatching and first feeding [FF]). We predicted that egg size would differ among morphs and be less variable as morphs diverge away from the ancestral anadromous phenotype. We also predicted that egg size would correlate with offspring size and developmental timing. We found morphs had different egg size, developmental timing, and size at hatching and FF. Egg size increased as phenotypic proximity to the ancestral anadromous phenotype decreased, with larger eggs generally giving rise to larger offspring, especially at FF, but egg size had no effect on developmental rate. The interaction between egg size and the environment may have a profound impact on offspring fitness, where the resulting differences in early life-history traits may act to initiate and/or maintain resource morphs diversification.

Adaptive divergence of the moor frog (Rana arvalis) along an acidification gradient.

BACKGROUND: Environmental stress can result in strong ecological and evolutionary effects on natural populations, but to what extent it drives adaptive divergence of natural populations is little explored. We used common garden experiments to study adaptive divergence in embryonic and larval fitness traits (embryonic survival, larval growth, and age and size at metamorphosis) in eight moor frog, Rana arvalis, populations inhabiting an acidification gradient (breeding pond pH 4.0 to 7.5) in southwestern Sweden. Embryos were raised until hatching at three (pH 4.0, 4.3 and 7.5) and larvae until metamorphosis at two (pH 4.3 and 7.5) pH treatments. To get insight into the putative selective agents along this environmental gradient, we measured relevant abiotic and biotic environmental variables from each breeding pond, and used linear models to test for phenotype-environment correlations. RESULTS: We found that acid origin populations had higher embryonic and larval acid tolerance (survival and larval period were less negatively affected by low pH), higher larval growth but slower larval development rates, and metamorphosed at a larger size. The phenotype-environment correlations revealed that divergence in embryonic acid tolerance and metamorphic size correlated most strongly with breeding pond pH, whereas divergence in larval period and larval growth correlated most strongly with latitude and predator density, respectively. CONCLUSION: Our results suggest that R. arvalis has diverged in response to pH mediated selection along this acidification gradient. However, as latitude and pH were closely spatially correlated in this study, further studies are needed to disentangle the specific agents of natural selection along acidification gradients. Our study highlights the need to consider the multiple interacting selective forces that drive adaptive divergence of natural populations along environmental stress gradients.

On biological evolution and environmental solutions.

Drawing insights from multiple disciplines is essential for finding integrative solutions that are required to tackle complex environmental problems. Human activities are causing unprecedented influence on global ecosystems, culminating in the loss of species and fundamental changes in the selective environments of organisms across the tree of life. Our collective understanding about biological evolution can help identify and mitigate many of the environmental problems in the Anthropocene. To this end, we propose a stronger integration of environmental sciences with evolutionary biology.