Contemporary temperature-driven divergence in a Nordic freshwater fish under conditions commonly thought to hinder adaptation

Toward a universal measure of robustness across model organs and systems

The development of an individual must be capable of resisting the harmful effects of internal and external perturbations. This capacity, called robustness, can make the difference between normal variation and disease. Some systems and organs are more resilient in their capacity to correct the effects of internal disturbances such as mutations. Similarly, organs and organisms differ in their capacity to be resilient against external disturbances, such as changes in temperature. Furthermore, all developmental systems must be somewhat flexible to permit evolutionary change, and understanding robustness requires a comparative framework. Over the last decades, most research on developmental robustness has been focusing on specific model systems and organs. Hence, we lack tools that would allow cross-species and cross-organ comparisons. Here, we emphasize the need for a uniform framework to experimentally test and quantify robustness across study systems and suggest that the analysis of fluctuating asymmetry might be a powerful proxy to do so. Such a comparative framework will ultimately help to resolve why and how organs of the same and different species differ in their sensitivity to internal (e.g., mutations) and external (e.g., temperature) perturbations and at what level of biological organization buffering capacities exist and therefore create robustness of the developmental system.

Population structure and genome-wide evolutionary signatures reveal putative climate-driven habitat change and local adaptation in the large yellow croaker

The large yellow croaker (Larimichthys crocea) is one of the most economically valuable marine fish in China and is a notable species in ecological studies owing to a serious collapse of wild germplasm in the past few decades. The stock division and species distribution, which have important implications for ecological protection, germplasm recovery, and fishery resource management, have been debated since the 1960s. However, it is still uncertain even how many stocks exist in this species. To address this, we evaluated the fine-scale genetic structure of large yellow croaker populations distributed along the eastern and southern Chinese coastline based on 7.64 million SNP markers. Compared with the widely accepted stock boundaries proposed in the 1960s, our results revealed that a climate-driven habitat change probably occurred between the Naozhou (Nanhai) Stock and the Ming-Yuedong (Mindong) Stock. The boundary between these two stocks might have shifted northwards from the Pearl River Estuary to the northern area of the Taiwan Strait, accompanied by highly asymmetric introgression. In addition, we found divergent landscapes of natural selection between the stocks inhabiting northern and southern areas. The northern population exhibited highly agminated signatures of strong natural selection in genes related to developmental processes, whereas moderate and interspersed selective signatures were detected in many immune-related genes in the southern populations. These findings establish the stock status and genome-wide evolutionary landscapes of large yellow croaker, providing a basis for conservation, fisheries management and further evolutionary biology studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00165-2.

The terroir of the finch: How spatial and temporal variation shapes phenotypic traits in DARWIN'S finches

The term terroir is used in viticulture to emphasize how the biotic and abiotic characteristics of a local site influence grape physiology and thus the properties of wine. In ecology and evolution, such terroir (i.e., the effect of space or “site”) is expected to play an important role in shaping phenotypic traits. Just how important is the pure spatial effect of terroir (e.g., differences between sites that persist across years) in comparison to temporal variation (e.g., differences between years that persist across sites), and the interaction between space and time (e.g., differences between sites change across years)? We answer this question by analyzing beak and body traits of 4388 medium ground finches (Geospiza fortis) collected across 10 years at three locations in Galápagos. Analyses of variance indicated that phenotypic variation was mostly explained by site for beak size (η² = 0.42) and body size (η² = 0.43), with a smaller contribution for beak shape (η² = 0.05) and body shape (η² = 0.12), but still higher compared to year and site‐by‐year effects. As such, the effect of terroir seems to be very strong in Darwin's finches, notwithstanding the oft‐emphasized interannual variation. However, these results changed dramatically when we excluded data from Daphne Major, indicating that the strong effect of terroir was mostly driven by that particular population. These phenotypic results were largely paralleled in analyses of environmental variables (rainfall and vegetation indices) expected to shape terroir in this system. These findings affirm the evolutionary importance of terroir, while also revealing its dependence on other factors, such as geographical isolation.

Genotype-by-environment interactions drive the maintenance of genetic variation in a Salmo trutta L. hybrid zone

Allopatric gene pools can evolve in different directions through adaptive and nonadaptive processes and are therefore a source of intraspecific diversity. The connection of these previously isolated gene pools through human intervention can lead to intraspecific diversity loss, through extirpation of native populations or hybridization. However, the mechanisms leading to these situations are not always explicitly documented and are thus rarely used to manage intraspecific diversity. In particular, genotype-by-environment (GxE) interactions can drive postzygotic reproductive isolation mechanisms that may result in a mosaic of diversity patterns, depending on the local environment. We test this hypothesis using a salmonid species (Salmo trutta) in the Mediterranean (MED) area, where intensive stocking from non-native Atlantic (ATL) origins has led to various outcomes of hybridization with the native MED lineage, going from MED resilience to total extirpation via full hybridization. We investigate patterns of offspring survival at egg stage in natural environments, based on parental genotypes in interaction with river temperature, to detect potential GxE interactions. Our results show a strong influence of maternal GxE interaction on embryonic survival, mediated by maternal effect through egg size, and a weak influence of paternal GxE interaction. In particular, when egg size is large and temperature is cold, the survival rate of offspring originating from MED females is three times higher than that of ATL females' offspring. Because river temperatures show contrast at small scale, this cold adaptation for MED females' offspring constitutes a potent postzygotic mechanism to explain small-scale spatial heterogeneity in diversity observed in MED areas where ATL fish have been stocked. It also indicates that management efforts could be specifically targeted at the environments that actively favor native intraspecific diversity through eco-evolutionary processes such as postzygotic selection.

Ambient temperature as a factor contributing to the developmental divergence in sympatric salmonids

Factors and mechanisms promoting resource-based radiation in animals still represent a main challenge to evolutionary biology. The modifications of phenotype tied with adaptive diversification may result from an environmentally related shift having occurred at the early stage of development. Here, we study the role of temperature dynamics on the reproductive sites in the early-life divergence and adaptive radiation of the salmonid fish Salvelinus malma dwelling in the Lake Kronotskoe basin (North-East Asia). Local sympatric charr ecomorphs demonstrate strict homing behaviour guiding the preordained distribution along tributaries and, hence, further development under different temperatures. We thoroughly assessed the annual temperature dynamics at the spawning grounds of each morph as compared to an ancestral anadromous morph. Then we carried out an experimental rearing of both under naturally diverging and uniformed temperatures. To compare the morphs' development under the dynamically changing temperatures, we have designed a method based on calculating the accumulated heat by the Arrhenius equation. The proposed equation shows a strong predictive power and, at the same time, is not bias-susceptible when the developmental temperature approximates 0°C. The temperature was found to significantly affect the charrs' early ontogeny, which underlies the divergence of developmental and growth rates between the morphs, as well as morph-specific ontogenetic adaptations to the spawning site's temperatures. As opposed to the endemic morphs from Lake Kronotskoe, the anadromous S. malma, being unexposed to selection оn highly specific reproduction conditions, showed a wide temperature tolerance, Our findings demonstrate that the hatch, onset timing of external feeding, and size dissimilarities between the sympatric morphs reveal themselves during the development under contrast temperatures. As a result of the observed developmental disparities, the morphs occupy specific definitive foraging niches in the lake.

Divergence in Thermal Physiology Could Contribute to Vertical Segregation in Intertidal Ecotypes of Littorina saxatilis

AbstractThermal stress is a potentially important selective agent in intertidal marine habitats, but the role that thermal tolerance might play in local adaptation across shore height has been underexplored. Northwest Spain is home to two morphologically distinct ecotypes of the periwinkle Littorina saxatilis, separated by shore height and subject to substantial differences in thermal stress exposure. However, despite other biotic and abiotic drivers of ecotype segregation being well studied, their thermal tolerance has not been previously characterized. We investigated thermal tolerance across multiple life history stages by employing the thermal death time (TDT) approach to determine (i) whether the two ecotypes differ in thermal tolerance and (ii) how any differences vary with life history stage. Adults of the two ecotypes differed in their thermal tolerance in line with their shore position: the upper-shore ecotype, which experiences more extreme temperatures, exhibited greater endurance of thermal stress compared with the lower-shore ecotype. This difference was most pronounced at the highest temperatures tested. The proximate physiological basis for these differences is unknown but likely due to a multifarious interaction of traits affecting different parts of the TDT curve. Differences in tolerance between ecotypes were less pronounced in early life history stages but increased with ontogeny, suggesting partial divergence of this trait during development. Thermal tolerance could potentially play an important role in maintaining population divergence and genetic segregation between the two ecotypes, since the increased thermal sensitivity of the lower-shore ecotype may limit its dispersal onto the upper shore and so restrict gene flow.

Genomic signatures of spatially divergent selection at clownfish range margins

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

Genomic differentiation and local adaptation on a microgeographic scale in a resident songbird

Elucidating forces capable of driving species diversification in the face of gene flow remains a key goal in evolutionary biology. Song sparrows, Melospiza melodia, occur as 25 subspecies in diverse habitats across North America, are among the continent's most widespread vertebrate species, and are exemplary of many highly variable species for which the conservation of locally adapted populations may be critical to their range-wide persistence. We focus here on six morphologically distinct subspecies resident in the San Francisco Bay region, including three salt-marsh endemics and three residents in upland and riparian habitats adjacent to the Bay. We used reduced-representation sequencing to generate 2,773 SNPs to explore genetic differentiation, spatial population structure, and demographic history. Clustering separated individuals from each of the six subspecies, indicating subtle differentiation at microgeographic scales. Evidence of limited gene flow and low nucleotide diversity across all six subspecies further supports a hypothesis of isolation among locally adapted populations. We suggest that natural selection for genotypes adapted to salt marsh environments and changes in demography over the past century have acted in concert to drive the patterns of diversification reported here. Our results offer evidence of microgeographic specialization in a highly polytypic bird species long discussed as a model of sympatric speciation and rapid adaptation, and they support the hypothesis that conserving locally adapted populations may be critical to the range-wide persistence of similarly highly variable species.

Temperature training improves transcriptional homeostasis after heat shock in juvenile Atlantic sturgeon (Acipenser oxyrinchus)

Exposure to high temperatures can lead to thermotolerance in fish, which is hypothesized to potentially improve post-release survival in species under restocking programs, like Atlantic sturgeon. The aim of this study was to determine whether Atlantic sturgeon juveniles exposed to a 4-week temperature treatment respond differently to a subsequent heat shock than juveniles exposed to heat shock for the first time (naive fish). Response to heat shock was assessed by mapping the liver transcriptome. In total, 838 unique contigs were differentially expressed between the trained and the control group (592 downregulated, 261 upregulated, and 15 down- or upregulated, depending on the condition), corresponding to genes involved in the response to heat, tissue damage, proteolysis, and metabolism. Temperature-trained fish showed 2-4-fold fewer dysregulated contigs than naive fish, indicating their ability to maintain and recover homeostasis faster. During heat shock, hspc1 was upregulated in both experimental groups, while hspa1 and dnaja4 were exclusively upregulated in the control. Overall, compensatory mechanisms were observed in addition to the heat shock response. Only two genes, fgg and apnl, were upregulated at nearly all timepoints in both groups. Peptidases were more strongly downregulated in control fish, which also showed a reduction in lipid metabolism during recovery. Keratins, pck1, gadd45ga, and gadd45gb were differentially expressed between trained and control fish, and due to their roles in tissue protection and ER stress reduction, they might be responsible for the maintenance of the transcriptional homeostasis observed in trained fish.

Cytosine methylation patterns suggest a role of methylation in plastic and adaptive responses to temperature in European grayling (Thymallus thymallus) populations

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

The role of common ancestry and gene flow in the evolution of human-directed play behaviour in dogs

Among-population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human-directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among-breed variance of human-directed play behaviour while also accounting for within-breed variance. We found that 75% of the among-breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human-directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within-species/among-population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among-breed variation of a behavioural phenotype.

Adaptations of early development to local spawning temperature in anadromous populations of pike (Esox lucius)

BACKGROUND

In the wake of climate change many environments will be exposed to increased and more variable temperatures. Knowledge about how species and populations respond to altered temperature regimes is therefore important to improve projections of how ecosystems will be affected by global warming, and to aid management. We conducted a common garden, split-brood temperature gradient (4.5 °C, 9.7 °C and 12.3 °C) experiment to study the effects of temperature in two populations (10 families from each population) of anadromous pike (Esox lucius) that normally experience different temperatures during spawning. Four offspring performance measures (hatching success, day degrees until hatching, fry survival, and fry body length) were compared between populations and among families.

RESULTS

Temperature affected all performance measures in a population-specific manner. Low temperature had a positive effect on the Harfjärden population and a negative effect on the Lervik population. Further, the effects of temperature differed among families within populations.

CONCLUSIONS

The population-specific responses to temperature indicate genetic differentiation in developmental plasticity between populations, and may reflect an adaptation to low temperature during early fry development in Harfjärden, where the stream leading up to the wetland dries out relatively early in the spring, forcing individuals to spawn early. The family-specific responses to temperature treatment indicate presence of genetic variation for developmental plasticity (G x E) within both populations. Protecting between- and within-population genetic variation for developmental plasticity and high temperature-related adaptive potential of early life history traits will be key to long-term viability and persistence in the face of continued climate change.

Comparative transcriptome analysis reveals potential evolutionary differences in adaptation of temperature and body shape among four Percidae species

Considering the divergent temperature habitats and morphological traits of four Percidae species: yellow perch (Perca flavescens), Eurasian perch (Perca fluviatilis), pike perch (Sander lucioperca), and ruffe (Gymnocephalus cernua), we stepped into the transcriptome level to discover genes and mechanisms that drive adaptation to different temperature environments and evolution in body shape. Based on 93,566 to 181,246 annotated unigenes of the four species, we identified 1,117 one-to-one orthologous genes and subsequently constructed the phylogenetic trees that are consistent with previous studies. Together with the tree, the ratios of nonsynonymous to synonymous substitutions presented decreased evolutionary rates from the D. rerio branch to the sub-branch clustered by P. flavescens and P. fluviatilis. The specific 93 fast-evolving genes and 57 positively selected genes in P. flavescens, compared with 22 shared fast-evolving genes among P. fluviatilis, G. cernua, and S. lucioperca, showed an intrinsic foundation that ensure its adaptation to the warmer Great Lakes and farther south, especially in functional terms like “Cul4-RING E3 ubiquitin ligase complex.” Meanwhile, the specific 78 fast-evolving genes and 41 positively selected genes in S. lucioperca drew a clear picture of how it evolved to a large and elongated body with camera-type eyes and muscle strength so that it could occupy the highest position in the food web. Overall, our results uncover genetic basis that support evolutionary adaptation of temperature and body shape in four Percid species, and could furthermore assist studies on environmental adaptation in fishes.

Carryover effects from natal habitat type upon competitive ability lead to trait divergence or source-sink dynamics

Local adaptation to rare habitats is difficult due to gene flow, but can occur if the habitat has higher productivity. Differences in offspring phenotypes have attracted little attention in this context. We model a scenario where the rarer habitat improves offspring's later competitive ability - a carryover effect that operates on top of local adaptation to one or the other habitat type. Assuming localised dispersal, so the offspring tend to settle in similar habitat to the natal type, the superior competitive ability of offspring remaining in the rarer habitat hampers immigration from the majority habitat. This initiates a positive feedback between local adaptation and trait divergence, which can thereafter be reinforced by coevolution with dispersal traits that match ecotype to habitat type. Rarity strengthens selection on dispersal traits and promotes linkage disequilibrium between locally adapted traits and ecotype-habitat matching dispersal. We propose that carryover effects may initiate isolation by ecology.

Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a QST - FST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater-dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site-specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half-sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family-specific heritable plasticity that could facilitate adaptive change.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

Mismatch between fishway operation and timing of fish movements: a risk for cascading effects in partial migration systems

Habitat fragmentation is a growing problem worldwide. Particularly in river systems, numerous dams and weirs hamper the movement of a wide variety of species. With the aim to preserve connectivity for fish, many barriers in river systems are equipped with fishways (also called fish passages or fish ladders). However, few fishways provide full connectivity. Here we hypothesized that restricted seasonal opening times of fishways can importantly reduce their effectiveness by interfering with the timing of fish migration, for both spring- and autumn-spawning species. We empirically tested our hypothesis, and discuss the possible eco-evolutionary consequences of affected migration timing. We analyzed movements of two salmonid fishes, spring-spawning European grayling (Thymallus thymallus) and autumn-spawning brown trout (Salmo trutta), in Norway's two largest river systems. We compared their timing of upstream passage through four fishways collected over 28 years with the timing of fish movements in unfragmented river sections as monitored by radiotelemetry. Confirming our hypothesis, late opening of fishways delayed the migration of European grayling in spring, and early closure of fishways blocked migration for brown trout on their way to spawning locations during late autumn. We show in a theoretical framework how restricted opening times of fishways can induce shifts from migratory to resident behavior in potamodromous partial migration systems, and propose that this can induce density-dependent effects among fish accumulating in lower regions of rivers. Hence, fragmentation may not only directly affect the migratory individuals in the population, but may also have effects that cascade downstream and alter circumstances for resident fish. Fishway functionality is inadequate if there is a mismatch between natural fish movements and fishway opening times in the same river system, with ecological and possibly evolutionary consequences for fish populations.

Thermal plasticity in farmed, wild and hybrid Atlantic salmon during early development: has domestication caused divergence in low temperature tolerance?

BACKGROUND

In the past three decades, millions of domesticated Atlantic salmon Salmo salar L. have escaped from farms into the wild. Their offspring display reduced survival in the natural environment, which demonstrates that gene-flow is likely to have a negative effect on wild populations. However, inter-population differences in introgression of farmed salmon have been observed, and the underlying ecological mechanisms remain enigmatic. We hypothesised that domestication-driven divergence in tolerance to low temperatures during early development may contribute to lower survival of farmed salmon offspring in the wild, which in turn, may influence patterns of introgression among populations exposed to different temperature regimes. We reared the offspring of 35 families of wild, farmed and hybrid origin at three temperatures (3.9, 5.6 and 12°C) from the onset of exogenous feeding and throughout their first summer. Thermal reaction norms for growth and survival were investigated along the gradient.

RESULTS

The main results of this study, which is based upon the analysis of juvenile salmon from five wild strains, two farmed strains and two hybrid strains, can be summarised as; (i) salmon of all origins were able to successfully initiate feeding at all temperatures and similar survival reaction norms were detected in all strains across the temperature gradient; (ii) deviating growth reaction norms were detected between strains, although this result was most likely due to an overall lack of growth in the lower temperature treatments.

CONCLUSIONS

This study revealed no evidence of domesticated-driven divergence in low temperature tolerance in Atlantic salmon during early development. Although the potential interaction between low temperature and other river-specific factors cannot be excluded, our results indicate that the reduced survival of farmed offspring in the wild is not explained by farmed salmon displaying impaired abilities to initiate feeding at low temperatures. We therefore suggest that the observed inter-population patterns of introgression are not low-temperature driven and that other ecological or biological factors may explain why detection of farmed salmon in wild rivers is not synonymous with introgression. In general, our results support the literature indicating that phenotypic plasticity instead of thermal adaption has been selected for in Atlantic salmon.

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

Genetic parentage analysis confirms a polygynandrous breeding system in the European grayling (Thymallus thymallus)

Knowing the breeding system of a species is important in order to understand individual variation in reproductive success. Large variation in reproductive success and thus reproductive skew strongly impacts on the effective number of breeders and thus the long-term effective population size (Ne). Fishes, in particular species belonging to the salmonid family, exhibit a wide diversity of breeding systems. In general, however, breeding systems are rarely studied in detail in the wild. Here we examine the breeding system of the spring-spawning European grayling Thymallus thymallus from a small Norwegian stream using parentage assignment based on the genotyping of 19 polymorphic microsatellite loci. In total 895 individual grayling fry and 154 mature grayling (57 females and 97 males) were genotyped. A total of 466 offspring were assigned a father, a mother, or a parent pair with a confidence of 90% or higher. Successfully reproducing males had on average 11.9 ± 13.3 (SD) offspring with on average 2.1 ± 1.2 partners, whereas successful females had on average 9.5 ± 12.8 offspring and 2.3 ± 1.5 partners. Parents with more partners also produced more offspring. Thus the grayling breeding system within this small stream revealed a polygynandrous breeding system, similar to what has been observed for many other salmonid fish species. The present study thus unambiguously corroborates a polygynadrous breeding system in the European grayling. This knowledge is critical for managing populations of this species, which has suffered significant local population declines throughout its range over the last several decades.

Islands within an island: repeated adaptive divergence in a single population

Physical barriers to gene flow were once viewed as prerequisites for adaptive evolutionary divergence. However, a growing body of theoretical and empirical work suggests that divergence can proceed within a single population. Here we document genetic structure and spatially replicated patterns of phenotypic divergence within a bird species endemic to 250 km(2) Santa Cruz Island, California, USA. Island scrub-jays (Aphelocoma insularis) in three separate stands of pine habitat had longer, shallower bills than jays in oak habitat, a pattern that mirrors adaptive differences between allopatric populations of the species' mainland congener. Variation in both bill measurements was heritable, and island scrub-jays mated nonrandomly with respect to bill morphology. The population was not panmictic; instead, we found a continuous pattern of isolation by distance across the east-west axis of the island, as well as a subtle genetic discontinuity across the boundary between the largest pine stand and adjacent oak habitat. The ecological factors that appear to have facilitated adaptive differentiation at such a fine scale--environmental heterogeneity and localized dispersal--are ubiquitous in nature. These findings support recent arguments that microgeographic patterns of adaptive divergence may be more common than currently appreciated, even in mobile taxonomic groups like birds.

Transcriptomic response to heat stress among ecologically divergent populations of redband trout

Background

As ectothermic organisms have evolved to differing aquatic climates, the molecular basis of thermal adaptation is a key area of research. In this study, we tested for differential transcriptional response of ecologically divergent populations of redband trout (Oncorhynchus mykiss gairdneri) that have evolved in desert and montane climates. Each pure strain and their F1 cross were reared in a common garden environment and exposed over four weeks to diel water temperatures that were similar to those experienced in desert climates within the species’ range. Gill tissues were collected from the three strains of fish (desert, montane, F1 crosses) at the peak of heat stress and tested for mRNA expression differences across the transcriptome with RNA-seq.

Results

Strong differences in transcriptomic response to heat stress were observed across strains confirming that fish from desert environments have evolved diverse mechanisms to cope with stressful environments. As expected, a large number of total transcripts (12,814) were differentially expressed in the study (FDR ≤ 0.05) with 2310 transcripts in common for all three strains, but the desert strain had a larger number of unique differentially expressed transcripts (2875) than the montane (1982) or the F1 (2355) strain. Strongly differentiated genes (>4 fold change and FDR ≤ 0.05) were particularly abundant in the desert strain (824 unique contigs) relative to the other two strains (montane = 58; F1 = 192).

Conclusions

This study demonstrated patterns of acclimation (i.e., phenotypic plasticity) within strains and evolutionary adaptation among strains in numerous genes throughout the transcriptome. Key stress response genes such as molecular chaperones (i.e., heat shock proteins) had adaptive patterns of gene expression among strains, but also a much higher number of metabolic and cellular process genes were differentially expressed in the desert strain demonstrating these biological pathways are critical for thermal adaptation to warm aquatic climates. The results of this study further elucidate the molecular basis for thermal adaptation in aquatic ecosystems and extend the potential for identifying genes that may be critical for adaptation to changing climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1246-5) contains supplementary material, which is available to authorized users.

AFLPs and mitochondrial haplotypes reveal local adaptation to extreme thermal environments in a freshwater gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.

Gene pleiotropy constrains gene expression changes in fish adapted to different thermal conditions

Understanding the factors that shape the evolution of gene expression is a central goal in biology, but the molecular mechanisms behind this remain controversial. A related major goal is ascertaining how such factors may affect the adaptive potential of a species or population. Here we demonstrate that temperature-driven gene expression changes in fish adapted to differing thermal environments are constrained by the level of gene pleiotropy estimated by either the number of protein interactions or gene biological processes. Genes with low pleiotropy levels were the main drivers of both plastic and evolutionary global expression profile changes, while highly pleiotropic genes had limited expression response to temperature treatment. Our study provides critical insights into the molecular mechanisms by which natural populations can adapt to changing environments. In addition to having important implications for climate change adaptation, these results suggest that gene pleiotropy should be considered more carefully when interpreting expression profiling data.

The factors that shape the evolution of gene expression and their role in adaptation are poorly understood. Here, Papakostas et al. show that gene pleiotropy constrains protein expression evolution in freshwater salmonids adapted to different temperatures.

Evidence for adaptive divergence of thermal responses among Bemisia tabaci populations from tropical Colombia following a recent invasion

There is an increasing evidence that populations of ectotherms can diverge genetically in response to different climatic conditions, both within their native range and (in the case of invasive species) in their new range. Here, we test for such divergence in invasive whitefly Bemisia tabaci populations in tropical Colombia, by considering heritable variation within and between populations in survival and fecundity under temperature stress, and by comparing population differences with patterns established from putatively neutral microsatellite markers. We detected significant differences among populations linked to mean temperature (for survival) and temperature variation (for fecundity) in local environments. A QST  - FST analysis indicated that phenotypic divergence was often larger than neutral expectations (QST  > FST ). Particularly, for survival after a sublethal heat shock, this divergence remained linked to the local mean temperature after controlling for neutral divergence. These findings point to rapid adaptation in invasive whitefly likely to contribute to its success as a pest species. Ongoing evolutionary divergence also provides challenges in predicting the likely impact of Bemisia in invaded regions.

Plastic and evolutionary responses to climate change in fish

The physical and ecological 'fingerprints' of anthropogenic climate change over the past century are now well documented in many environments and taxa. We reviewed the evidence for phenotypic responses to recent climate change in fish. Changes in the timing of migration and reproduction, age at maturity, age at juvenile migration, growth, survival and fecundity were associated primarily with changes in temperature. Although these traits can evolve rapidly, only two studies attributed phenotypic changes formally to evolutionary mechanisms. The correlation-based methods most frequently employed point largely to 'fine-grained' population responses to environmental variability (i.e. rapid phenotypic changes relative to generation time), consistent with plastic mechanisms. Ultimately, many species will likely adapt to long-term warming trends overlaid on natural climate oscillations. Considering the strong plasticity in all traits studied, we recommend development and expanded use of methods capable of detecting evolutionary change, such as the long term study of selection coefficients and temporal shifts in reaction norms, and increased attention to forecasting adaptive change in response to the synergistic interactions of the multiple selection pressures likely to be associated with climate change.

Testing for local adaptation in brown trout using reciprocal transplants

BACKGROUND

Local adaptation can drive the divergence of populations but identification of the traits under selection remains a major challenge in evolutionary biology. Reciprocal transplant experiments are ideal tests of local adaptation, yet rarely used for higher vertebrates because of the mobility and potential invasiveness of non-native organisms. Here, we reciprocally transplanted 2500 brown trout (Salmo trutta) embryos from five populations to investigate local adaptation in early life history traits. Embryos were bred in a full-factorial design and raised in natural riverbeds until emergence. Customized egg capsules were used to simulate the natural redd environment and allowed tracking the fate of every individual until retrieval. We predicted that 1) within sites, native populations would outperform non-natives, and 2) across sites, populations would show higher performance at 'home' compared to 'away' sites.

RESULTS

There was no evidence for local adaptation but we found large differences in survival and hatching rates between sites, indicative of considerable variation in habitat quality. Survival was generally high across all populations (55% ± 3%), but ranged from 4% to 89% between sites. Average hatching rate was 25% ± 3% across populations ranging from 0% to 62% between sites.

CONCLUSION

This study provides rare empirical data on variation in early life history traits in a population network of a salmonid, and large-scale breeding and transplantation experiments like ours provide powerful tests for local adaptation. Despite the recently reported genetic and morphological differences between the populations in our study area, local adaptation at the embryo level is small, non-existent, or confined to ecological conditions that our experiment could not capture.

Sex-biased genetic component distribution among populations: additive genetic and maternal contributions to phenotypic differences among populations of Chinook salmon

An approach frequently used to demonstrate a genetic basis for population-level phenotypic differences is to employ common garden rearing designs, where observed differences are assumed to be attributable to primarily additive genetic effects. Here, in two common garden experiments, we employed factorial breeding designs between wild and domestic, and among wild populations of Chinook salmon (Oncorhynchus tshawytscha). We measured the contribution of additive (V(A)) and maternal (V(M)) effects to the observed population differences for 17 life history and fitness-related traits. Our results show that, in general, maternal effects contribute more to phenotypic differences among populations than additive genetic effects. These results suggest that maternal effects are important in population phenotypic differentiation and also signify that the inclusion of the maternal source of variation is critical when employing models to test population differences in salmon, such as in local adaptation studies.

Contemporary divergence in early life history in grayling (Thymallus thymallus)

Background

Following colonization of new habitats and subsequent selection, adaptation to environmental conditions might be expected to be rapid. In a mountain lake in Norway, Lesjaskogsvatnet, more than 20 distinct spawning demes of grayling have been established since the lake was colonized, some 20-25 generations ago. The demes spawn in tributaries consistently exhibiting either colder or warmer temperature conditions during spawning in spring and subsequent early development during early summer. In order to explore the degree of temperature-related divergence in early development, a multi-temperature common-garden experiment was performed on embryos from four different demes experiencing different spring temperatures.

Results

Early developmental characters were measured to test if individuals from the four demes respond differently to the treatment temperatures. There was clear evidence of among-deme differences (genotype - environment interactions) in larval growth and yolk-to-body-size conversion efficiency. Under the cold treatment regime, larval growth rates were highest for individuals belonging to cold streams. Individuals from warm streams had the highest yolk-consumption rate under cold conditions. As a consequence, yolk-to-body-mass conversion efficiency was highest for cold-deme individuals under cold conditions. As we observed response parallelism between individuals from demes belonging to similar thermal groups for these traits, some of the differentiation seems likely to result from local adaptation

Conclusion

The observed differences in length at age during early larval development most likely have a genetic component, even though both directional and random processes are likely to have influenced evolutionary change in the demes under study.

Strong gene flow and lack of stable population structure in the face of rapid adaptation to local temperature in a spring-spawning salmonid, the European grayling (Thymallus thymallus)

Gene flow has the potential to both constrain and facilitate adaptation to local environmental conditions. The early stages of population divergence can be unstable because of fluctuating levels of gene flow. Investigating temporal variation in gene flow during the initial stages of population divergence can therefore provide insights to the role of gene flow in adaptive evolution. Since the recent colonization of Lake Lesjaskogsvatnet in Norway by European grayling (Thymallus thymallus), local populations have been established in over 20 tributaries. Multiple founder events appear to have resulted in reduced neutral variation. Nevertheless, there is evidence for local adaptation in early life-history traits to different temperature regimes. In this study, microsatellite data from almost a decade of sampling were assessed to infer population structuring and its temporal stability. Several alternative analyses indicated that spatial variation explained 2-3 times more of the divergence in the system than temporal variation. Over all samples and years, there was a significant correlation between genetic and geographic distance. However, decomposed pairwise regression analysis revealed differing patterns of genetic structure among local populations and indicated that migration outweighs genetic drift in the majority of populations. In addition, isolation by distance was observable in only three of the six years, and signals of population bottlenecks were observed in the majority of samples. Combined, the results suggest that habitat-specific adaptation in this system has preceded the development of consistent population substructuring in the face of high levels of gene flow from divergent environments.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

A climate for contemporary evolution

A new study of divergence in freshwater fish provides strong evidence of rapid, temperature-mediated adaptation. This study is particularly important in the ongoing debate over the extent and significance of evolutionary response to climate change because divergence has occurred in relatively few generations in spite of ongoing gene flow and in the aftermath of a significant genetic bottleneck, factors that have previously been considered obstacles to evolution. Climate change may thus be more likely to foster contemporary evolutionary responses than has been anticipated, and I argue here for the importance of investigating their possible occurrence.

See Research article: http://www.biomedcentral.com/1471-2148/10/350/abstract