Ecological adaptation in European eels is based on phenotypic plasticity

Putting Structural Variants Into Practice: The Role of Chromosomal Inversions in the Management of Marine Environments

Major threats to marine species and ecosystems include overfishing, invasive species, pollution and climate change. The changing climate not only imposes direct threats through the impacts of severe marine heatwaves, cyclones and ocean acidification but also complicates fisheries and invasive species management by driving species range shifts. The dynamic nature of these threats means that the future of our oceans will depend on the ability of species to adapt. This has led to calls for genetic interventions focussed on enhancing species' adaptive capacity, including translocations, restocking and selective breeding. Assessing the benefits and risks of such approaches requires an improved understanding of the genetic architecture of adaptive variation, not only in relation to climate-resilient phenotypes but also locally adapted populations and the fitness of hybrids. Large structural genetic variants such as chromosomal inversions play an important role in local adaptation by linking multiple adaptive loci. Consequently, inversions are likely to be particularly important when managing for adaptive capacity. However, under some circumstances, they also accumulate deleterious mutations, potentially increasing the risk of inbreeding depression. Genetic management that takes account of these dual roles on fitness is likely to be more effective at ensuring population persistence. We summarise evolutionary factors influencing adaptive and deleterious variation of inversions, review inversions found in marine taxa, and provide a framework to predict the consequences of ignoring inversions in key management scenarios. We conclude by describing practical methods to bridge the gap between evolutionary theory and practical application of inversions in conservation.

Isotope geolocation and population genomics in Vanessa cardui: Short- and long-distance migrants are genetically undifferentiated

The painted lady butterfly Vanessa cardui is renowned for its virtually cosmopolitan distribution and the remarkable long-distance migrations as part of its annual, multigenerational migratory cycle. In winter, V. cardui individuals inhabit breeding grounds north and south of the Sahara, suggesting distinct migratory behaviors within the species as individuals migrate southward from Europe in the autumn. However, the evolutionary and ecological factors shaping these differences in migratory behavior remain largely unexplored. Here, we performed whole-genome resequencing and analyzed the hydrogen and strontium isotopes of 40 V. cardui individuals simultaneously collected in the autumn from regions both north and south of the Sahara. Our investigation revealed two main migratory groups: (i) short-distance migrants, journeying from temperate Europe to the circum-Mediterranean region and (ii) long-distance migrants, originating from Europe, crossing the Mediterranean Sea and Sahara, and reaching West Africa, covering up to over 4,000 km. Despite these stark differences in migration distance, a genome-wide analysis revealed that short- and long-distance migrants belong to a single intercontinental panmictic population extending from northern Europe to sub-Saharan Africa. Contrary to common biogeographic patterns, the Sahara is not a catalyst for population structuring in this species. No significant genetic differentiation or signs of adaptation and selection were observed between the two migratory phenotypes. Nonetheless, two individuals, who were early arrivals to West Africa covering longer migration distances, exhibited some genetic differentiation. The lack of genetic structure between short- and long-distance migrants suggests that migration distance in V. cardui is a plastic response to environmental conditions.

Genomic Footprints of Hybridisation in North Atlantic Eels (Anguilla anguilla and A. rostrata)

Understanding interspecific introgressive hybridisation and the biological significance of introgressed variation remains an important goal in population genomics. European (Anguilla anguilla) and American eel (A. rostrata) represent a remarkable case of hybridisation. Both are panmictic and spawn in partial sympatry in the Sargasso Sea, occasionally producing viable, fertile hybrids, primarily found in Iceland. We studied introgressive hybridisation from American into European eel using whole-genome sequences of 78 individuals, including European, American and 21 putative hybrid eels. Previous studies using few genetic markers could not resolve whether hybridisation involved simple unidirectional backcrossing or a more complex hybrid swarm scenario. However, local ancestry inference along individual chromosomes revealed that Icelandic hybrids were primarily F1 or first-generation backcrosses towards European eel, with some showing more complex backcrossing. All European eels outside Iceland contained short chromosomal blocks from American eel, indicating a porous genome. We found no evidence for previously hypothesised geographical gradients of introgression in European eel outside Iceland. Several chromosomal regions showed high interspecific divergence, but haplotype blocks introgressed from American eel were identified both within and outside these regions. There was little correspondence between regions of high relative (FST) and absolute divergence (dXY), with the former reflecting selective sweeps within species or reduced recombination rather than barrier loci. A single genomic region showed evidence of repeated introgression from American into European eel under positive selection in both species. The study illustrates that species can maintain genetic integrity despite porous genomes and that standing variation in one species can potentially be available for future adaptive responses in the other species.

Eels' individual migratory behavior stems from a complex syndrome involving cognition, behavior, physiology, and life history

Variability within species is key for adaptability and biological evolution. To understand individualities in the context of animal movement, we focused on one of the most remarkable migrations-the journey of the endangered European eel from their birthplace in the Sargasso Sea to freshwater environments. Laboratory observations unveiled a continuum of diverse phenotypes of migrating eels: Some displayed a heightened tendency to swim against a constant water flow, while others a greater propensity to climb obstacles. Looking for the biological underpinnings of this migratory diversity, we characterized the eels' individual differences in traits of four key domains: life history, physiology, behavior, and cognition, among which we found significant variance and interconnectedness. Upon reducing this variance to its primary multivariate axes, we found that these predict the migratory types. Eels with 1) low exploration, high activity, low boldness, and high lateralization; 2) strong lateralization, enhanced quantitative abilities, short problem-solving time, high boldness, and low growth rates; or 3) enhanced problem-solving, reduced spatial learning, high cognitive flexibility, and shorter time to solve the cognitive tasks were more likely to display the climbing migratory type. Field sampling revealed how specific traits' combinations seemed to influence the distribution of eels once they begin to settle in the freshwater environment. Our study underscores the impressive diversity of individuals during this critical migration, emphasizing an intrinsic connection to multidomain trait variance. Preserving this diversity becomes paramount, as it likely contributes to the resilience and adaptability of endangered migratory species.

Genetic architecture of long-distance migration and population genomics of the endangered Japanese eel

The Japanese eel (Anguilla japonica), a flagship anguillid species for conservation, is known for its long-distance-oriented migration. However, our understanding of the genetic architecture underlying long-distance migration and population genomic characteristics of A. japonica is still limited. Here, we generated a high-quality chromosome-level genome assembly and conducted whole-genome resequencing of 218 individuals to explore these aspects. Strong signals of selection were found on genes involved in long-distance aerobic exercise and navigation, which might be associated with evolutionary adaptation to long-distance migrations. Low genetic diversity was detected, which might result from genetic drift associated with demographic declines. Both mitochondrial and nuclear genomic datasets supported the existence of a single panmictic population for Japanese eel, despite signals of single-generation selection. Candidate genes for local selection involved in functions like development and circadian rhythm. The findings can provide insights to adaptative evolution to long-distance migration and inform conservation efforts for A. japonica.

Limited Parallelism in Genetic Adaptation to Brackish Water Bodies in European Sprat and Atlantic Herring

The European sprat is a small plankton-feeding clupeid present in the northeastern Atlantic Ocean, in the Mediterranean Sea, and in the brackish Baltic Sea and Black Sea. This species is the target of a major fishery and, therefore, an accurate characterization of its genetic population structure is crucial to delineate proper stock assessments that aid ensuring the fishery's sustainability. Here, we present (i) a draft genome assembly, (ii) pooled whole genome sequencing of 19 population samples covering most of the species' distribution range, and (iii) the design and test of a single nucleotide polymorphism (SNP)-chip resource and use this to validate the population structure inferred from pooled sequencing. These approaches revealed, using the populations sampled here, three major groups of European sprat: Oceanic, Coastal, and Brackish with limited differentiation within groups even over wide geographical stretches. Genetic structure is largely driven by six large putative inversions that differentiate Oceanic and Brackish sprats, while Coastal populations display intermediate frequencies of haplotypes at each locus. Interestingly, populations from the Baltic and the Black Seas share similar frequencies of haplotypes at these putative inversions despite their distant geographic location. The closely related clupeids European sprat and Atlantic herring both show genetic adaptation to the brackish Baltic Sea, providing an opportunity to explore the extent of genetic parallelism. This analysis revealed limited parallelism because out of 125 independent loci detected in the Atlantic herring, three showed sharp signals of selection that overlapped between the two species and contained single genes such as PRLRA, which encodes the receptor for prolactin, a freshwater-adapting hormone in euryhaline species, and THRB, a receptor for thyroid hormones, important both for metabolic regulation and the development of red cone photoreceptors.

How Fish Population Genomics Can Promote Sustainable Fisheries: A Road Map

Maintenance of genetic diversity in marine fishes targeted by commercial fishing is a grand challenge for the future. Most of these species are abundant and therefore important for marine ecosystems and food security. Here, we present a road map of how population genomics can promote sustainable fisheries. In these species, the development of reference genomes and whole genome sequencing is key, because genetic differentiation at neutral loci is usually low due to large population sizes and gene flow. First, baseline allele frequencies representing genetically differentiated populations within species must be established. These can then be used to accurately determine the composition of mixed samples, forming the basis for population demographic analysis to inform sustainably set fish quotas. SNP-chip analysis is a cost-effective method for determining baseline allele frequencies and for population identification in mixed samples. Finally, we describe how genetic marker analysis can transform stock identification and management.

Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

Effects of the invasive swim bladder parasite Anguillicola crassus on health and condition indicators in the European eel

Parasites negatively affect biological processes within their hosts, which may alter for example health, growth, and reproductive ability. Non-native invasive parasites, in particular, may have large effects on the endemic hosts, given that the hosts lack evolved specific defences against such parasites. The swim bladder nematode Anguillicola crassus, an invasive parasite originating from Asia, is found in the European eel (Anguilla anguilla, L. 1758), since the 1980s. We investigated whether A. crassus affected several indicators related to health of the European eel (spleen- and liver size, body fat content and relative condition). Our results indicate that during the continental residency of the eels, infection by A. crassus had no major negative impacts on the investigated health indicators at the generally low infection intensities present in this study (median 2-3 visible parasites). Given that many of the adult eels were found to have swim bladder damage, concerns about their spawning migration through deeper oceanic environments can still be raised. To allow further investigations, we suggest that quantification of swim bladder damage should be implemented in eel-monitoring programs. Compared to other parasite pressure parameters, swim bladder damage provides additional information about past infections and future problems.

The genomic basis and environmental correlates of local adaptation in the Atlantic horse mackerel (Trachurus trachurus)

Understanding how populations adapt to their environment is increasingly important to prevent biodiversity loss due to overexploitation and climate change. Here we studied the population structure and genetic basis of local adaptation of Atlantic horse mackerel, a commercially and ecologically important marine fish that has one of the widest distributions in the eastern Atlantic. We analyzed whole-genome sequencing and environmental data of samples collected from the North Sea to North Africa and the western Mediterranean Sea. Our genomic approach indicated low population structure with a major split between the Mediterranean Sea and the Atlantic Ocean and between locations north and south of mid-Portugal. Populations from the North Sea are the most genetically distinct in the Atlantic. We discovered that most population structure patterns are driven by a few highly differentiated putatively adaptive loci. Seven loci discriminate the North Sea, two the Mediterranean Sea, and a large putative inversion (9.9 Mb) on chromosome 21 underlines the north-south divide and distinguishes North Africa. A genome-environment association analysis indicates that mean seawater temperature and temperature range, or factors correlated to them, are likely the main environmental drivers of local adaptation. Our genomic data broadly support the current stock divisions, but highlight areas of potential mixing, which require further investigation. Moreover, we demonstrate that as few as 17 highly informative SNPs can genetically discriminate the North Sea and North African samples from neighboring populations. Our study highlights the importance of both, life history and climate-related selective pressures in shaping population structure patterns in marine fish. It also supports that chromosomal rearrangements play a key role in local adaptation with gene flow. This study provides the basis for more accurate delineation of the horse mackerel stocks and paves the way for improving stock assessments.

Adaptive, maladaptive, neutral, or absent plasticity: Hidden caveats of reaction norms

Adaptive phenotypic plasticity may improve the response of individuals when faced with new environmental conditions. Typically, empirical evidence for plasticity is based on phenotypic reaction norms obtained in reciprocal transplant experiments. In such experiments, individuals from their native environment are transplanted into a different environment, and a number of trait values, potentially implicated in individuals' response to the new environment, are measured. However, the interpretations of reaction norms may differ depending on the nature of the assessed traits, which may not be known beforehand. For example, for traits that contribute to local adaptation, adaptive plasticity implies nonzero slopes of reaction norms. By contrast, for traits that are correlated to fitness, high tolerance to different environments (possibly due to adaptive plasticity in traits that contribute to adaptation) may, instead, result in flat reaction norms. Here we investigate reaction norms for adaptive versus fitness-correlated traits and how they may affect the conclusions regarding the contribution of plasticity. To this end, we first simulate range expansion along an environmental gradient where plasticity evolves to different values locally and then perform reciprocal transplant experiments in silico. We show that reaction norms alone cannot inform us whether the assessed trait exhibits locally adaptive, maladaptive, neutral, or no plasticity, without any additional knowledge of the traits assessed and species' biology. We use the insights from the model to analyse and interpret empirical data from reciprocal transplant experiments involving the marine isopod Idotea balthica sampled from two geographical locations with different salinities, concluding that the low-salinity population likely has reduced adaptive plasticity relative to the high-salinity population. Overall, we conclude that, when interpreting results from reciprocal transplant experiments, it is necessary to consider whether traits assessed are locally adaptive with respect to the environmental variable accounted for in the experiments or correlated to fitness.

Rapid and transient evolution of local adaptation to seasonal host fruits in an invasive pest fly

Both local adaptation and adaptive phenotypic plasticity can influence the match between phenotypic traits and local environmental conditions. Theory predicts that environments stable for multiple generations promote local adaptation, whereas highly heterogeneous environments favor adaptive phenotypic plasticity. However, when environments have periods of stability mixed with heterogeneity, the relative importance of local adaptation and adaptive phenotypic plasticity is unclear. Here, we used Drosophila suzukii as a model system to evaluate the relative influence of genetic and plastic effects on the match of populations to environments with periods of stability from three to four generations. This invasive pest insect can develop within different fruits, and persists throughout the year in a given location on a succession of distinct host fruits, each one being available for only a few generations. Using reciprocal common environment experiments of natural D. suzukii populations collected from cherry, strawberry, and blackberry, we found that both oviposition preference and offspring performance were higher on medium made with the fruit from which the population originated than on media made with alternative fruits. This pattern, which remained after two generations in the laboratory, was analyzed using a statistical method we developed to quantify the contributions of local adaptation and adaptive plasticity in determining fitness. Altogether, we found that genetic effects (local adaptation) dominate over plastic effects (adaptive phenotypic plasticity). Our study demonstrates that spatially and temporally variable selection does not prevent the rapid evolution of local adaptation in natural populations. The speed and strength of adaptation may be facilitated by several mechanisms including a large effective population size and strong selective pressures imposed by host plants.

Morphological and physiological variation of soybean seedlings in response to shade

Soybean (Glycine max) is a legume species that is widely used in intercropping. Quantitative analyses of plasticity and genetic differences in soybean would improve the selection and breeding of soybean in intercropping. Here, we used data of 20 varieties from one year artificial shading experiment and one year intercropping experiment to characterize the morphological and physiological traits of soybean seedlings grown under shade and full sun light conditions. Our results showed that shade significantly decreased biomass, leaf area, stem diameter, fraction of dry mass in petiole, leaf mass per unit area, chlorophyll a/b ratio, net photosynthetic rate per unit area at PAR of 500 μmol m-2 s-1 and 1,200 μmol m-2 s-1 of soybean seedling, but significantly increased plant height, fraction of dry mass in stem and chlorophyll content. Light × variety interaction was significant for all measured traits, light effect contributed more than variety effect. The biomass of soybean seedlings was positively correlated with leaf area and stem diameter under both shade and full sunlight conditions, but not correlated with plant height and net photosynthetic rate. The top five (62.75% variation explained) most important explanatory variables of plasticity of biomass were that the plasticity of leaf area, leaf area ratio, leaflet area, plant height and chlorophyll content, whose total weight were 1, 0.9, 0.3, 0.2, 0.19, respectively. The plasticity of biomass was positively correlated with plasticity of leaf area and leaflet area but significant negative correlated with plasticity of plant height. The principal component one account for 42.45% variation explain. A cluster analysis further indicated that soybean cultivars were classified into three groups and cultivars; Jiandebaimaodou, Gongdou 2, and Guixia 3 with the maximum plasticity of biomass. These results suggest that for soybean seedlings grown under shade increasing the capacity for light interception by larger leaf area is more vital than light searching (plant height) and light conversion (photosynthetic rate).

Genome-wide methylation in the panmictic European eel (Anguilla anguilla)

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, likely reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggests interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight hox clusters, and the results suggest important roles of methylation at hox genes for adaptive processes.

Adaptation insights from comparative transcriptome analysis of two Opisthopappus species in the Taihang mountains

Opisthopappus is a major wild source of Asteraceae with resistance to cold and drought. Two species of this genus (Opisthopappus taihangensis and O. longilobus) have been employed as model systems to address the evolutionary history of perennial herb biomes in the Taihang Mountains of China. However, further studies on the adaptive divergence processes of these two species are currently impeded by the lack of genomic resources. To elucidate the molecular mechanisms involved, a comparative analysis of these two species was conducted. Among the identified transcription factors, the bHLH members were most prevalent, which exhibited significantly different expression levels in the terpenoid metabolic pathway. O. longilobus showed higher level of expression than did O. taihangensis in terms of terpenes biosynthesis and metabolism, particularly monoterpenoids and diterpenoids. Analyses of the positive selection genes (PSGs) identified from O. taihangensis and O. longilobus revealed that 1203 genes were related to adaptative divergence, which were under rapid evolution and/or have signs of positive selection. Differential expressions of PSG occurred primarily in the mitochondrial electron transport, starch degradation, secondary metabolism, as well as nucleotide synthesis and S-metabolism pathway processes. Several PSGs were obviously differentially expressed in terpenes biosynthesis that might result in the fragrances divergence between O. longilobus and O. taihangensis, which would provide insights into adaptation of the two species to different environments that characterized by sub-humid warm temperate and temperate continental monsoon climates. The comparative analysis for these two species in Opisthopappus not only revealed how the divergence occurred from molecular perspective, but also provided novel insights into how differential adaptations occurred in Taihang Mountains.

The role of phenotypic plasticity in the establishment of range margins

It has been argued that adaptive phenotypic plasticity may facilitate range expansions over spatially and temporally variable environments. However, plasticity may induce fitness costs. This may hinder the evolution of plasticity. Earlier modelling studies examined the role of plasticity during range expansions of populations with fixed genetic variance. However, genetic variance evolves in natural populations. This may critically alter model outcomes. We ask: how does the capacity for plasticity in populations with evolving genetic variance alter range margins that populations without the capacity for plasticity are expected to attain? We answered this question using computer simulations and analytical approximations. We found a critical plasticity cost above which the capacity for plasticity has no impact on the expected range of the population. Below the critical cost, by contrast, plasticity facilitates range expansion, extending the range in comparison to that expected for populations without plasticity. We further found that populations may evolve plasticity to buffer temporal environmental fluctuations, but only when the plasticity cost is below the critical cost. Thus, the cost of plasticity is a key factor involved in range expansions of populations with the potential to express plastic response in the adaptive trait. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Comparative genomics and signatures of selection in North Atlantic eels

Comparative genomic approaches can identify putative private and shared signatures of selection. We performed a comparative genomic study of North Atlantic eels, European eel (Anguilla Anguilla) and American eel (A. rostrata). The two sister species are nearly undistinguishable at the phenotypic level and despite a wide non-overlapping continental distribution, they spawn in partial sympatry in the Sargasso Sea. Taking advantage of the newly assembled and annotated genome, we used genome wide RAD sequencing data of 359 individuals retrieved from Sequence Nucleotide Archive and state-of-the-art statistic tests to identify putative genomic signatures of selection in North Atlantic eels. First, using the F_ST and XP-EHH methods, we detected apparent islands of divergence on a total of 7 chromosomes, particularly on chromosomes 6 and 10. Gene ontology analyses suggested candidate genes mainly related to energy production, development and regulation, which could reflect strong selection on traits related to eel migration and larval duration time. Gene effect prediction using SNPeff showed a high number of SNPs in noncoding regions, pointing to a possible regulatory role. Second, using the iHS method we detected shared regions under selection on a total of 11 chromosomes. Several hypotheses might account for the detection of shared islands of selection in North Atlantic eels, including parallel evolution due to adaptation to similar environments and introgression. Future comparative genomic studies will be needed to further clarify the causes and consequences of introgression, including the directionality of these introgression events.

The changing times of Europe's largest remaining commercially harvested population of eel Anguilla anguilla L

This study quantifies the processes involved in regulating the European eel population of Lough Neagh, a lake in Northern Ireland. The relationship between glass eel input and silver eel output for the 1923-1997 cohorts was best described by a Beverton-Holt stock recruitment model. Glass eel input time series was not complete and was thus derived from the relationship between catches elsewhere in Europe and Lough Neagh, together with the addition of stocked glass eel. Silver eel output was the sum of silver eel escapement, catch and yellow eel catch converted to silver eel equivalents. Natural mortality increased with glass eel density, ranging from 0.017 to 0.142 year^-1 . The mean carrying capacity increased from ≈3.25 M silver eels (≈26 kg ha^-1 ) for the 1923-1943 cohorts to ≈5.0 M (≈40 kg ha^-1 ) for the 1948-1971 cohorts before regressing back to ≈3.25 M. The total silver eel output was highest during the late 1970s/early 1980s at 35-45 kg ha^-1  year^-1 and lowest during the early years of the 20th century and is currently at 10-15 kg ha^-1  year^-1 . The findings are discussed in relation to (a) the ecological changes that have occurred within the lough, associated with eutrophication and the introduction of roach (Rutilus rutilus L.), and (b) the decline of the wider European eel stock across its distribution range. The findings from this study have relevance for the wider management of the European eel stock.