Standing genetic variation fuels rapid adaptation to ocean acidification

Evaluation of potential bioaccumulation of Bisphenol A in the mussel Mytilus galloprovincialis

Bisphenol A (BPA) is an endocrine disruptor, widespread into the marine environment, which can be accumulated in marine organisms, representing a serious threat to human health, even if it is present at low concentrations. This study presents a new methodology for detecting BPA in the edible mussel Mytilus galloprovincialis which, compared to the existing ones, has several advantages. It is cheaper, faster, independent of the amount of organism material considered, and can be used even for the analysis of very small samples: LOD of 0.13 μg g-1 was obtained extracting 50 mg of tissue. The detection of BPA was obtained by means of a tailored method able to spot the contaminant spiked or bioaccumulated by living mussels at different concentrations and time points. Under environmental conditions inedible mussels (not depurated) bioaccumulated up to 0.6 μg g-1 (d.w.). Our methodology was applied to in vivo experiments whereas the concentration of BPA detected in mussels' tissues was found to depend on the quantity administered. This methodology may provide the baseline for future research and improve environmental regulation and risk assessment protocols.

Seasonal genetic variation and genetic structure of Spodoptera exigua in Liaoning Province, Northeast China: insights from 11 years of microsatellite data

Background

The beet armyworm (BAW), Spodoptera exigua, is a destructive migratory pest worldwide that has caused severe economic losses in China's major crop-producing regions. To control this pest effectively, it is crucial to investigate its seasonal genetic variation and population genetic structure in northern China.

Methods

In this study, we used eight nuclear microsatellite loci to investigate the seasonal genetic variation and genetic structure of BAW in Shenyang, Liaoning Province, Northeast China, from 2012-2022, collected from a single location on Welsh onion.

Results

Microsatellite data revealed moderate levels of genetic variation among 50 seasonal populations of BAW sampled from 2012-2022, along with significant genetic differentiation among these populations. Neighbor-joining dendrograms, STRUCTURE analysis, and principal coordinate analysis (PCoA) revealed two genetically distinct groups: the SY2012-2018 group and the SY2019-2022 group. Our results revealed seasonal variation in the genetic subconstruction at this location, which may be related to the presence of different migratory individuals throughout the year. Accordingly, our unique insights into the population genetics of BAW will contribute to the development of effective management strategies for this migratory pest.

Directed Evolution of Silicatein Reveals Biomineralization Synergism between Protein Sequences

Biomineralization is a green synthesis route for a variety of metal nanoparticles. Silicatein is a biomineralization protein originally found in marine sponge Tethya aurantia that converts inorganic precursors to metal oxide nanoparticles. In this work, we investigate the popular catalytic triad hypothesis and implement directed evolution with the aim to improve the solubility and kinetics of silicatein to enable increased nanoparticle synthesis. Site-directed mutagenesis with catalytic triad residues did not abolish biomineralization activity, aligning with the results seen in one previous study. Recombinant production of silicatein and mutants in Escherichia coli following library generation and a survival screen yielded several mutant proteins with augmented biomineralization activity. Sequence analysis of these mutant proteins reveals multiple sequences within a single cell that contribute to enhanced biomineralization. Combined with the sequence analysis of silicateins from different marine sponges, these results suggest the protein is permissive to wide sequence variations and that multiple protein sequences act synergistically for enhanced biomineralization.

Shallow-water mussels (Mytilus galloprovincialis) adapt to deep-sea environment through transcriptomic and metagenomic insights

Recent studies have unveiled the deep sea as a rich biosphere, populated by species descended from shallow-water ancestors post-mass extinctions. Research on genomic evolution and microbial symbiosis has shed light on how these species thrive in extreme deep-sea conditions. However, early adaptation stages, particularly the roles of conserved genes and symbiotic microbes, remain inadequately understood. This study examined transcriptomic and microbiome changes in shallow-water mussels Mytilus galloprovincialis exposed to deep-sea conditions at the Site-F cold seep in the South China Sea. Results reveal complex gene expression adjustments in stress response, immune defense, homeostasis, and energy metabolism pathways during adaptation. After 10 days of deep-sea exposure, shallow-water mussels and their microbial communities closely resembled those of native deep-sea mussels, demonstrating host and microbiome convergence in response to adaptive shifts. Notably, methanotrophic bacteria, key symbionts in native deep-sea mussels, emerged as a dominant group in the exposed mussels. Host genes involved in immune recognition and endocytosis correlated significantly with the abundance of these bacteria. Overall, our analyses provide insights into adaptive transcriptional regulation and microbiome dynamics of mussels in deep-sea environments, highlighting the roles of conserved genes and microbial community shifts in adapting to extreme environments.

Genomic signals of adaptation to a natural CO2 gradient over a striking microgeographic scale

Our study explores genomic signs of adaptation in A. lixula to different water pH conditions. To achieve this, we analysed the genomics variation of A. lixula individuals living across a natural pH gradient in Canary Islands, Spain. We use a 2b-RADseq protocol with 74 samples from sites with varying pH levels (from 7.3 to 7.9 during low tide) and included a control site. We identified 14,883 SNPs, with 432 identified as candidate SNPs under selection to pH variations through redundancy analysis. While all SNPs indicated genomic homogeneity, the 432 candidate SNPs under selection displayed genomic differences among sites and along the pH gradient. Out of these 432 loci, 17 were annotated using published A. lixula transcriptomes, involved in biological functions such as growth. Therefore, our findings suggest local adaptation in A. lixula populations to acidification in CO2 vents, even over short distances of 75 m, underscoring their potential resistance to future Ocean Acidification.

Experimental evolution suggests rapid assembly of the 'selfing syndrome' from standing variation in Mimulus guttatus

Ecological and evolutionary changes are likely to occur rapidly when outcrossing populations experience pollinator loss. However, the number and identify of plant traits that will respond to this form of selection, as well as the overall predictability of evolutionary responses, remain unclear. We experimentally evolved 20 large replicate populations of Mimulus guttatus for 10 generations under three treatments: pure outcrossing, mixed mating (10% outcrossing) and pure selfing. These populations were founded from the same genetically diverse and outcrossing natural population. After 10 generations, all measured traits evolved with flower size, phenology, and reproductive traits diverging consistently among mating system treatments. Autogamy increased dramatically in the selfing treatment, but the magnitude of adaptation only becomes clear once inbreeding depression is factored out. Selfing treatment plants evolved reduced stigma-anther separation, and also exhibited declines in flower size and per-flower reproductive capacity. Flower size also declined in selfing populations but this was driven mainly by inbreeding depression and cannot be attributed to adaptation towards the selfing syndrome. Generally, the mixed mating populations evolved trait values intermediate to the fully selfing and outcrossing populations. Overall, our experimental treatments reiterated differences that have been documented in interspecific comparisons between selfing and outcrossing species pairs. Given that such contrasts involve species separated by thousands or even millions of generations, it is noteworthy that large evolutionary responses were obtained from genetic variation segregating within a single natural population.

Natural acidified marine systems: Lessons and predictions

Natural acidified marine systems (ASs) are environments with relatively low pH levels due to natural causes such as volcanic activity, geochemical reactions, and biological processes. These systems act as natural laboratories for the study of the effects of ocean acidification, allowing for the observation of long-term ecological and evolutionary responses. Understanding these systems is crucial for predicting the effects of anthropogenic ocean acidification (OA) on marine ecosystems. There are 23 ASs in which scientific research has shown significant parallelisms in their results worldwide, such as the disappearance of calcareous organisms and the loss of species with key ecological functions under OA conditions. Future research should emphasize continuous collaboration among teams, as well as public access to oceanographic and biological data along with the monitoring of environmental variables at each AS. To preserve these areas, it is imperative to employ non-destructive methods and protect them as human heritage sites.

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Environmental epigenetics: Exploring phenotypic plasticity and transgenerational adaptation in fish

Epigenetics plays a vital role in the interaction between living organisms and their environment by regulating biological functions and phenotypic plasticity. Considering that most aquaculture activities take place in open or natural habitats that are vulnerable to environmental changes. Promising findings from recent research conducted on various aquaculture species have provided preliminary evidence suggesting a link between epigenetic mechanisms and economically valuable characteristics. Environmental stressors, including climate changes (thermal stress, hypoxia, and water salinity), anthropogenic impacts such as (pesticides, crude oil pollution, nutritional impacts, and heavy metal) and abiotic factors (infectious diseases), can directly trigger epigenetic modifications in fish. While experiments have confirmed that many epigenetic alterations caused by environmental factors have plastic responses, some can be permanently integrated into the genome through genetic integration and promoting rapid transgenerational adaptation in fish. These environmental factors might cause irregular DNA methylation patterns in genes related to many biological events leading to organs dysfunction by inducing alterations in genes related to oxidative stress or apoptosis. Moreover, these environmental issues alter DNA/histone methylation leading to decreased reproductive competence. This review emphasizes the importance of understanding the effects of environmentally relevant issues on the epigenetic regulation of phenotypic variations in fish. The goal is to expand our knowledge of how epigenetics can either facilitate or hinder species' adaptation to these adverse conditions. Furthermore, this review outlines the areas that warrant further investigation in understanding epigenetic reactions to various environmental issues.

Chromosomal-scale genome assembly of the Mediterranean mussel Mytilus galloprovincialis

The Mediterranean mussel, Mytilus galloprovincialis, is a significant marine bivalve species that has ecological and economic importance. This species is robustly resilient and highly invasive. Despite the scientific and commercial interest in studying its biology and aquaculture, there remains a need for a high-quality, chromosome-scale reference genome. In this study, we have assembled a high-quality chromosome-scale reference genome for M. galloprovincialis. The total length of our reference genome is 1.41 Gb, with a scaffold N50 sequence length of 96.9 Mb. BUSCO analysis revealed a 97.5% completeness based on complete BUSCOs. Compared to the four other available M. galloprovincialis assemblies, the assembly described here is dramatically improved in both contiguity and completeness. This new reference genome will greatly contribute to a deeper understanding of the resilience and invasiveness of M. galloprovincialis.

Predicting species invasiveness with genomic data: Is genomic offset related to establishment probability?

Predicting the risk of establishment and spread of populations outside their native range represents a major challenge in evolutionary biology. Various methods have recently been developed to estimate population (mal)adaptation to a new environment with genomic data via so-called Genomic Offset (GO) statistics. These approaches are particularly promising for studying invasive species but have still rarely been used in this context. Here, we evaluated the relationship between GO and the establishment probability of a population in a new environment using both in silico and empirical data. First, we designed invasion simulations to evaluate the ability to predict establishment probability of two GO computation methods (Geometric GO and Gradient Forest) under several conditions. Additionally, we aimed to evaluate the interpretability of absolute Geometric GO values, which theoretically represent the adaptive genetic distance between populations from distinct environments. Second, utilizing public empirical data from the crop pest species Bactrocera tryoni, a fruit fly native from Northern Australia, we computed GO between "source" populations and a diverse range of locations within invaded areas. This practical application of GO within the context of a biological invasion underscores its potential in providing insights and guiding recommendations for future invasion risk assessment. Overall, our results suggest that GO statistics represent good predictors of the establishment probability and may thus inform invasion risk, although the influence of several factors on prediction performance (e.g., propagule pressure or admixture) will need further investigation.

Adaptive colonization across a parasitism-mutualism gradient

Adaptive colonization is a process wherein a colonizing population exhibits an adaptive change in response to a novel environment, which may be critical to its establishment. To date, theoretical models of adaptive colonization have been based on single-species introductions. However, given their pervasiveness, symbionts will frequently be co-introduced with their hosts to novel areas. We present an individual-based model to investigate adaptive colonization by hosts and their symbionts across a parasite-mutualist continuum. The host must adapt in order to establish itself in the novel habitat, and the symbiont must adapt to track evolutionary change in the host. First, we classify the qualitative shifts in the outcome that can potentially be driven by non-neutral effects of the symbiont-host interaction into three main types: parasite-driven co-extinction, parasite release, and mutualistic facilitation. Second, we provide a detailed description of a specific example for each type of shift. Third, we disentangle how the interplay between symbiont transmissibility, host migration, and selection strength determines: (a) which type of shift is more likely to occur and (b) the size of the interaction effects necessary to produce it. Overall, we demonstrate the crucial role of host and symbiont dispersal scales in shaping the impacts of parasitism and mutualism on adaptive colonization.

Environment-induced heritable variations are common in Arabidopsis thaliana

Parental or ancestral environments can induce heritable phenotypic changes, but whether such environment-induced heritable changes are a common phenomenon remains unexplored. Here, we subject 14 genotypes of Arabidopsis thaliana to 10 different environmental treatments and observe phenotypic and genome-wide gene expression changes over four successive generations. We find that all treatments caused heritable phenotypic and gene expression changes, with a substantial proportion stably transmitted over all observed generations. Intriguingly, the susceptibility of a genotype to environmental inductions could be predicted based on the transposon abundance in the genome. Our study thus challenges the classic view that the environment only participates in the selection of heritable variation and suggests that the environment can play a significant role in generating of heritable variations.

Transcriptomic Response of the Atlantic Surfclam (Spisula solidissima) to Acute Heat Stress

There is clear evidence that the oceans are warming due to anthropogenic climate change, and the northeastern coast of USA contains some of the fastest warming areas. This warming is projected to continue with serious biological and social ramifications for fisheries and aquaculture. One species particularly vulnerable to warming is the Atlantic surfclam (Spisula solidissima). The surfclam is a critically important species, linking marine food webs and supporting a productive, lucrative, and sustainable fishery. The surfclam is also emerging as an attractive candidate for aquaculture diversification, but the warming of shallow coastal farms threatens the expansion of surfclam aquaculture. Little is known about the adaptive potential of surfclams to cope with ocean warming. In this study, the surfclam transcriptome under heat stress was examined. Two groups of surfclams were subjected to heat stress to assess how artificial selection may alter gene expression. One group of clams had been selected for greater heat tolerance (HS) and the other was composed of random control clams (RC). After a 6-h exposure to 16 or 29 °C, gill transcriptome expression profiles of the four temperature/group combinations were determined by RNA sequencing and compared. When surfclams experienced heat stress, they exhibited upregulation of heat shock proteins (HSPs), inhibitors of apoptosis (IAPs), and other stress-response related genes. RC clams differentially expressed 1.7 times more genes than HS clams, yet HS clams had a stronger response of key stress response genes, including HSPs, IAPs, and genes involved with mitigating oxidative stress. The findings imply that the HS clams have a more effective response to heat stress after undergoing the initial selection event due to genetic differences created by the selection, epigenetic memory of the first heat shock, or both. This work provides insights into how surfclams adapt to heat stress and should inform future breeding programs that attempt to breed surfclam for greater heat tolerance, and ultimately bring greater resiliency to shellfish farms.

Ancestors' Gift: Parental Early Exposure to the Environmentally Realistic Pesticide Mixture Drives Offspring Phenotype in a Larger Extent Than Direct Exposure in the Pacific Oyster, Crassostrea gigas

Marine organisms are threatened by the presence of pesticides in coastal waters. Among them, the Pacific oyster is one of the most studied invertebrates in marine ecotoxicology where numerous studies highlighted the multiscale impacts of pesticides. In the past few years, a growing body of literature has reported the epigenetic outcomes of xenobiotics. Because DNA methylation is an epigenetic mark implicated in organism development and is meiotically heritable, it raises the question of the multigenerational implications of xenobiotic-induced epigenetic alterations. Therefore, we performed a multigenerational exposure to an environmentally relevant mixture of 18 pesticides (nominal sum concentration: 2.85 μg·L-1) during embryo-larval stages (0-48 hpf) of a second generation (F1) for which parents where already exposed or not in F0. Gene expression, DNA methylation, and physiological end points were assessed throughout the life cycle of individuals. Overall, the multigenerational effect has a greater influence on the phenotype than the exposure itself. Thus, multigenerational phenotypic effects were observed: individuals descending from exposed parents exhibited lower epinephrine-induced metamorphosis and field survival rates. At the molecular level, RNA-seq and Methyl-seq data analyses performed in gastrula embryos and metamorphosis-competent pediveliger (MCP) larvae revealed a clear F0 treatment-dependent discrimination. Some genes implicated into shell secretion and immunity exhibited F1:F0 treatment interaction patterns (e.g., Calm and Myd88). Those results suggest that low chronic environmental pesticide contamination can alter organisms beyond the individual scale level and have long-term adaptive implications.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

Projection of current and future distribution of adaptive genetic units in an alpine ungulate

Climate projections predict major changes in alpine environments by the end of the 21st century. To avoid climate-induced maladaptation and extinction, many animal populations will either need to move to more suitable habitats or adapt in situ to novel conditions. Since populations of a species exhibit genetic variation related to local adaptation, it is important to incorporate this variation into predictive models to help assess the ability of the species to survive climate change. Here, we evaluate how the adaptive genetic variation of a mountain ungulate-the Northern chamois (Rupicapra rupicapra)-could be impacted by future global warming. Based on genotype-environment association analyses of 429 chamois using a ddRAD sequencing approach, we identified genetic variation associated with climatic gradients across the European Alps. We then delineated adaptive genetic units and projected the optimal distribution of these adaptive groups in the future. Our results suggest the presence of local adaptation to climate in Northern chamois with similar genetic adaptive responses in geographically distant but climatically similar populations. Furthermore, our results predict that future climatic changes will modify the Northern chamois adaptive landscape considerably, with various degrees of maladaptation risk.

The impacts of thermal heterogeneity across microhabitats on post-settlement selection of intertidal mussels

Rapid genetic selection is critical for allowing natural populations to adapt to different thermal environments such as those that occur across intertidal microhabitats with high degrees of thermal heterogeneity. To address the question of how thermal regimes influence selection and adaptation in the intertidal black mussel Mytilisepta virgata, we continuously recorded environmental temperatures in both tidal pools and emergent rock microhabitats and then assessed genetic differentiation, gene expression patterns, RNA editing level, and cardiac performance. Our results showed that the subpopulations in the tidal pool and on emergent rocks had different genetic structures and exhibited different physiological and molecular responses to high-temperature stress. These results indicate that environmental heterogeneity across microhabitats is important for driving genetic differentiation and shed light on the importance of post-settlement selection for adaptively modifying the genetic composition and thermal responses of these intertidal mussels.

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

CO2 enrichment and excess nitrogen supply synergistically increase toxicity of marine dinoflagellate Alexandrium minutum

Discharges of CO2 and nutrients by anthropogenic activities have notable contributions to CO2 enrichment and eutrophication in coastal systems. Following our previous study that toxic dinoflagellate Alexandrium minutum will increase their growth rate and cellular toxicity under elevated levels of CO2, we further examined the joint effects of CO2 enrichment and excess nitrogen supply through a 29-day experiment under three CO2 levels (400, 800 and 1200 ppm) with a high N/P ratio of 80. It was found that the two factors have synergistical effects in promoting the increase of cellular toxicity. There were remarkable increases in toxin biosynthesis of paralytic shellfish toxins when both intracellular and extracellular toxins were considered. Under the joint impacts of CO2 enrichment and excess nitrogen supply, the apparent transformation from gonyautoxins2/3 to gonyautoxins1/4, with much higher toxicity and lower rate of release, is likely to be another major factor accounting for the increasing toxicity. The increasing growth rate and cellular toxicity of A. minutum under the scenarios with elevated concentrations of both CO2 and nitrogen in coastal systems in the future will increase the risks associated with such toxic algal blooms.

Earth beyond six of nine planetary boundaries

This planetary boundaries framework update finds that six of the nine boundaries are transgressed, suggesting that Earth is now well outside of the safe operating space for humanity. Ocean acidification is close to being breached, while aerosol loading regionally exceeds the boundary. Stratospheric ozone levels have slightly recovered. The transgression level has increased for all boundaries earlier identified as overstepped. As primary production drives Earth system biosphere functions, human appropriation of net primary production is proposed as a control variable for functional biosphere integrity. This boundary is also transgressed. Earth system modeling of different levels of the transgression of the climate and land system change boundaries illustrates that these anthropogenic impacts on Earth system must be considered in a systemic context.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Genomic signatures suggesting adaptation to ocean acidification in a coral holobiont from volcanic CO2 seeps

Ocean acidification, caused by anthropogenic CO₂ emissions, is predicted to have major consequences for reef-building corals, jeopardizing the scaffolding of the most biodiverse marine habitats. However, whether corals can adapt to ocean acidification and how remains unclear. We addressed these questions by re-examining transcriptome and genome data of Acropora millepora coral holobionts from volcanic CO₂ seeps with end-of-century pH levels. We show that adaptation to ocean acidification is a wholistic process involving the three main compartments of the coral holobiont. We identified 441 coral host candidate adaptive genes involved in calcification, response to acidification, and symbiosis; population genetic differentiation in dinoflagellate photosymbionts; and consistent transcriptional microbiome activity despite microbial community shifts. Coral holobionts from natural analogues to future ocean conditions harbor beneficial genetic variants with far-reaching rapid adaptation potential. In the face of climate change, these populations require immediate conservation strategies as they could become key to coral reef survival.

Mixed strain pathogen populations accelerate the evolution of antibiotic resistance in patients

Antibiotic resistance poses a global health threat, but the within-host drivers of resistance remain poorly understood. Pathogen populations are often assumed to be clonal within hosts, and resistance is thought to emerge due to selection for de novo variants. Here we show that mixed strain populations are common in the opportunistic pathogen P. aeruginosa. Crucially, resistance evolves rapidly in patients colonized by multiple strains through selection for pre-existing resistant strains. In contrast, resistance evolves sporadically in patients colonized by single strains due to selection for novel resistance mutations. However, strong trade-offs between resistance and growth rate occur in mixed strain populations, suggesting that within-host diversity can also drive the loss of resistance in the absence of antibiotic treatment. In summary, we show that the within-host diversity of pathogen populations plays a key role in shaping the emergence of resistance in response to treatment.

Shedding light on variation in reproductive success through studies of population genetic structure in a Southeast Pacific Coast mussel

Phylogeography often focuses on the spatial dimension of genetic diversity, rarely including the temporal dynamics occurring interannually among local populations, which can provide insight into past variations in reproductive success. Currently, there is an intense aquaculture industry of Mytilus spp. on the Southeast Pacific Coast which depends entirely on the spat released by natural populations forming a relevant and sensitive social-ecological system. Temporal and spatial spat variability from natural mussel beds could be related to interannual reproductive dynamics with variable reproductive success and recruitment, which leave genetic signatures. Temporal and spatial genetic structure was evaluated in six natural beds in the Southeast Pacific (from 39°25'S to 43°07'S) on the most abundant and widespread Mytilus lineage detected, Mytilus cf. chilensis, in 4 consecutive years. Analyses included data from >180 individuals per year, with a total of 751 (mitochondrial COI) and 747 (nuclear H1) individuals, respectively. Overall, both markers showed high haplotype diversity and low spatial and temporal genetic differentiation. Likely, the high dispersal capacity of Mytilus cf. chilensis maintains population homogeneity and prevents diversity erosion. The slight differences in genetic variance of COI were better explained by differences among sites (space), and conversely, the H1 genetic variance was better explained by interannual (temporal) comparisons, which could explain temporal variability in spat availability. This study highlights the important insights achieved with the evaluation of both temporal and spatial population genetic structures in marine species with high reproductive output, which can condition the success and sustainability of the relevant social-ecological system.

Predicting evolution in experimental range expansions of an aquatic model system

Predicting range expansion dynamics is an important goal of both fundamental and applied research in conservation and global change biology. However, this is challenging if ecological and evolutionary processes occur on the same time scale. Using the freshwater ciliate Paramecium caudatum, we combined experimental evolution and mathematical modeling to assess the predictability of evolutionary change during range expansions. In the experiment, we followed ecological dynamics and trait evolution in independently replicated microcosm populations in range core and front treatments, where episodes of natural dispersal alternated with periods of population growth. These eco-evolutionary conditions were recreated in a predictive mathematical model, parametrized with dispersal and growth data of the 20 founder strains in the experiment. We found that short-term evolution was driven by selection for increased dispersal in the front treatment and general selection for higher growth rates in all treatments. There was a good quantitative match between predicted and observed trait changes. Phenotypic divergence was further mirrored by genetic divergence between range core and front treatments. In each treatment, we found the repeated fixation of the same cytochrome c oxidase I (COI) marker genotype, carried by strains that also were the most likely winners in our model. Long-term evolution in the experimental range front lines resulted in the emergence of a dispersal syndrome, namely a competition-colonization trade-off. Altogether, both model and experiment highlight the potential importance of dispersal evolution as a driver of range expansions. Thus, evolution at range fronts may follow predictable trajectories, at least for simple scenarios, and predicting these dynamics may be possible from knowledge of few key parameters.

Sea of opportunities: marine genomics in an era of global environmental change

Overexploitation of natural resources and pollution of seas, acidification of the ocean, and rising temperatures all contribute to the destruction of marine habitats and, in 2015, the protection of the ocean became one of the UN Sustainable Development Goals (SDG 14: Life Below Water). This collection aims to highlight the molecular genetic changes currently happening in marine organisms.

Transgenerational endocrine disruptor effects of cadmium in zebrafish and contribution of standing epigenetic variation to adaptation

Evidence has emerged that environmentally-induced epigenetic changes can have long-lasting effects on gene transcription across generations. These recent findings highlight the need to investigate the transgenerational impacts of pollutants to assess their long term effects on populations. In this study, we investigated the transgenerational effect of cadmium on zebrafish across 4 generations. A first whole methylome approach carried out on fish of the first two generations led us to focus our investigations on the estradiol receptor alpha gene (esr1). We observed a sex-dependent transgenerational inheritance of Cd-induced DNA methylation changes up to the last generation. These changes were associated with single nucleotide polymorphisms (SNPs) that were themselves at the origin of the creation or deletion of methylation sites. Thus, Cd-induced genetic selection gave rise to DNA methylation changes. We also analyzed the transcription level of various sections of esr1 as well as estrogen responsive genes. While Cd triggered transgenerational disorders, Cd-induced epigenetic changes in esr1 contributed to the rapid transgenerational adaptation of fish to Cd. Our results provide insight into the processes underpinning rapid adaptation and highlight the need to maintain genetic diversity within natural populations to bolster the resilience of species faced with the global environmental changes.

Using conservation genetics to prioritise management options for an endangered songbird

Genetic data can be highly informative for answering questions relevant to practical conservation efforts, but remain one of the most neglected aspects of species recovery plans. Framing genetic questions with reference to practical and tractable conservation objectives can help bypass this limitation of the application of genetics in conservation. Using a single-nucleotide polymorphism dataset from reduced-representation sequencing (DArTSeq), we conducted a genetic assessment of remnant populations of the endangered forty-spotted pardalote (Pardalotus quadragintus), a songbird endemic to Tasmania, Australia. Our objectives were to inform strategies for the conservation of genetic diversity in the species and estimate effective population sizes and patterns of inter-population movement to identify management units relevant to population conservation and habitat restoration. We show population genetic structure and identify two small populations on mainland Tasmania as 'satellites' of larger Bruny Island populations connected by migration. Our data identify management units for conservation objectives relating to genetic diversity and habitat restoration. Although our results do not indicate the immediate need to genetically manage populations, the small effective population sizes we estimated for some populations indicate that they are vulnerable to genetic drift, highlighting the urgent need to implement habitat restoration to increase population size and to conduct genetic monitoring. We discuss how our genetic assessment can be used to inform management interventions for the forty-spotted pardalote and show that by assessing contemporary genetic aspects, valuable information for conservation planning and decision-making can be produced to guide actions that account for genetic diversity and increase chances of recovery in species of conservation concern.

RNAi Silencing of the Biomineralization Gene Perlucin Impairs Oyster Ability to Cope with Ocean Acidification

Calcifying marine organisms, including the eastern oyster (Crassostrea virginica), are vulnerable to ocean acidification (OA) because it is more difficult to precipitate calcium carbonate (CaCO₃). Previous investigations of the molecular mechanisms associated with resilience to OA in C. virginica demonstrated significant differences in single nucleotide polymorphism and gene expression profiles among oysters reared under ambient and OA conditions. Converged evidence generated by both of these approaches highlighted the role of genes related to biomineralization, including perlucins. Here, gene silencing via RNA interference (RNAi) was used to evaluate the protective role of a perlucin gene under OA stress. Larvae were exposed to short dicer-substrate small interfering RNA (DsiRNA-perlucin) to silence the target gene or to one of two control treatments (control DsiRNA or seawater) before cultivation under OA (pH ~7.3) or ambient (pH ~8.2) conditions. Two transfection experiments were performed in parallel, one during fertilization and one during early larval development (6 h post-fertilization), before larval viability, size, development, and shell mineralization were monitored. Silenced oysters under acidification stress were the smallest, had shell abnormalities, and had significantly reduced shell mineralization, thereby suggesting that perlucin significantly helps larvae mitigate the effects of OA.

Deep resilience: An evolutionary perspective on calcification in an age of ocean acidification

The success of today's calcifying organisms in tomorrow's oceans depends, in part, on the resilience of their skeletons to ocean acidification. To the extent this statement is true there is reason to have hope. Many marine calcifiers demonstrate resilience when exposed to environments that mimic near-term ocean acidification. The fossil record similarly suggests that resilience in skeletons has increased dramatically over geologic time. This "deep resilience" is seen in the long-term stability of skeletal chemistry, as well as a decreasing correlation between skeletal mineralogy and extinction risk over time. Such resilience over geologic timescales is often attributed to genetic canalization-the hardening of genetic pathways due to the evolution of increasingly complex regulatory systems. But paradoxically, our current knowledge on biomineralization genetics suggests an opposing trend, where genes are co-opted and shuffled at an evolutionarily rapid pace. In this paper we consider two possible mechanisms driving deep resilience in skeletons that fall outside of genetic canalization: microbial co-regulation and macroevolutionary trends in skeleton structure. The mechanisms driving deep resilience should be considered when creating risk assessments for marine organisms facing ocean acidification and provide a wealth of research avenues to explore.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria

The increasing concentration of CO₂ in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO₂ and reared under acidified (pH ~ 7.5, pCO₂ ~ 1200 ppm) or control (pH ~ 7.9, pCO₂ ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

Selective effects of a short transient environmental fluctuation on a natural population

Natural populations experience continuous and often transient changes of environmental conditions. These in turn may result in fluctuating selection pressures leading to variable demographic and evolutionary population responses. Rapid adaptation as short-term response to a sudden environmental change has in several cases been attributed to polygenic traits, but the underlying genomic dynamics and architecture are poorly understood. In this study, we took advantage of a natural experiment in an insect population of the non-biting midge Chironomus riparius by monitoring genome-wide allele frequencies before and after a cold snap event. Whole genome pooled sequencing of time series samples revealed 10 selected haplotypes carrying ancient polymorphisms, partially with signatures of balancing selection. By constantly cold exposing genetically variable individuals in the laboratory, we could demonstrate with whole genome resequencing (i) that among the survivors, the same alleles rose in frequency as in the wild, and (ii) that the identified variants additively predicted fitness (survival time) of its bearers. Finally, by simultaneously sequencing the genome and the transcriptome of cold exposed individuals we could tentatively link some of the selected SNPs to the cis- and trans-regulation of genes and pathways known to be involved in cold response of insects, such as cytochrome P450 and fatty acid metabolism. Altogether, our results shed light on the strength and speed of selection in natural populations and the genomic architecture of its underlying polygenic trait. Population genomic time series data thus appear as promising tool for measuring the selective tracking of fluctuating selection in natural populations.

Within-population variation in body size plasticity in response to combined nutritional and thermal stress is partially independent from variation in development time

Ongoing climate change has forced animals to face changing thermal and nutritional environments. Animals can adjust to such combinations of stressors via plasticity. Body size is a key trait influencing organismal fitness, and plasticity in this trait in response to nutritional and thermal conditions varies among genetically diverse, locally adapted populations. The standing genetic variation within a population can also influence the extent of body size plasticity. We generated near-isogenic lines from a newly collected population of Drosophila melanogaster at the mid-point of east coast Australia and assayed body size for all lines in combinations of thermal and nutritional stress. We found that isogenic lines showed distinct underlying patterns of body size plasticity in response to temperature and nutrition that were often different from the overall population response. We then tested whether plasticity in development time could explain, and therefore regulate, variation in body size to these combinations of environmental conditions. We selected five genotypes that showed the greatest variation in response to combined thermal and nutritional stress and assessed the correlation between response of developmental time and body size. While we found significant genetic variation in development time plasticity, it was a poor predictor of body size among genotypes. Our results therefore suggest that multiple developmental pathways could generate genetic variation in body size plasticity. Our study emphasizes the need to better understand genetic variation in plasticity within a population, which will help determine the potential for populations to adapt to ongoing environmental change.

Genetic diversity goals and targets have improved, but remain insufficient for clear implementation of the post-2020 global biodiversity framework

Genetic diversity among and within populations of all species is necessary for people and nature to survive and thrive in a changing world. Over the past three years, commitments for conserving genetic diversity have become more ambitious and specific under the Convention on Biological Diversity's (CBD) draft post-2020 global biodiversity framework (GBF). This Perspective article comments on how goals and targets of the GBF have evolved, the improvements that are still needed, lessons learned from this process, and connections between goals and targets and the actions and reporting that will be needed to maintain, protect, manage and monitor genetic diversity. It is possible and necessary that the GBF strives to maintain genetic diversity within and among populations of all species, to restore genetic connectivity, and to develop national genetic conservation strategies, and to report on these using proposed, feasible indicators.

Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes

Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.

Monitoring genetic diversity with new indicators applied to an alpine freshwater top predator

Genetic diversity is the basis for population adaptation and long-term survival, yet rarely considered in biodiversity monitoring. One key issue is the need for useful and straightforward indicators of genetic diversity. We monitored genetic diversity over 40 years (1970-2010) in metapopulations of brown trout (Salmo trutta) inhabiting 27 small mountain lakes representing 10 lake systems in central Sweden using >1200 fish per time point. We tested six newly proposed indicators; three were designed for broad, international use in the UN Convention on Biological Diversity (CBD) and are currently applied in several countries. The other three were recently elaborated for national use by a Swedish science-management effort and applied for the first time here. The Swedish indicators use molecular genetic data to monitor genetic diversity within and between populations (indicators ΔH and ΔF_ST , respectively) and assess the effective population size (N_e -indicator). We identified 29 genetically distinct populations, all retained over time. Twelve of the 27 lakes harboured more than one population indicating that brown trout biodiversity hidden as cryptic, sympatric populations are more common than recognized. The N_e indicator showed values below the threshold (N_e  ≤ 500) in 20 populations with five showing N_e  < 100. Statistically significant genetic diversity reductions occurred in several populations. Metapopulation structure appears to buffer against diversity loss; applying the indicators to metapopulations suggest mostly acceptable rates of change in all but one system. The CBD indicators agreed with the Swedish ones but provided less detail. All these indicators are appropriate for managers to initiate monitoring of genetic biodiversity.

Increased genetic diversity loss and genetic differentiation in a model marine diatom adapted to ocean warming compared to high CO2

Although high CO2 and warming could act interactively on marine phytoplankton, little is known about the molecular basis for this interaction on an evolutionary scale. Here we explored the adaptation to high CO2 in combination with warming in a model marine diatom Phaeodactylum tricornutum. Whole-genome re-sequencing identifies, in comparison to populations grown under control conditions, a larger genetic diversity loss and a higher genetic differentiation in the populations adapted for 2 years to warming than in those adapted to high CO2. However, this diversity loss was less under high CO2 combined with warming, suggesting that the evolution driven by warming was constrained by high CO2. By integrating genomics, transcriptomics, and physiological data, we found that the underlying molecular basis for this constraint is associated with the expression of genes involved in some key metabolic pathways or biological processes, such as the glyoxylate pathway, amino acid and fatty acid metabolism, and diel variability. Our results shed new light on the evolutionary responses of marine phytoplankton to multiple environmental changes in the context of global change and provide new insights into the molecular basis underpinning interactions among those multiple drivers.

Genetic and epigenetic interplay allows rapid transgenerational adaptation to metal pollution in zebrafish

Despite still being a matter of debate, there is growing evidence that pollutant-induced epigenetic changes can be propagated across generations. Whereas such modifications could have long-lasting effects on organisms and even on population, environmentally relevant data from long-term exposure combined with follow-up through multiple generations remain scarce for non-mammalian species. We performed a transgenerational experiment comprising four successive generations of zebrafish. Only fish from the first generation were exposed to an environmentally realistic concentration of cadmium (Cd). Using a whole methylome analysis, we first identified the DNA regions that were differentially methylated in response to Cd exposure and common to fish of the first two generations. Among them, we then focused our investigations on the exon 3 (ex3) of the cep19 gene. We indeed recorded transgenerational growth disorders in Cd-exposed fish, and a mutation in this exon is known to cause morbid obesity in mammals. Its methylation level was thus determined in zebrafish from all the four generations by means of a targeted and base resolution method. We observed a transgenerational inheritance of Cd-induced DNA methylation changes up to the fourth generation. However, these changes were closely associated with genetic variations, mainly a single nucleotide polymorphism. This single nucleotide polymorphism was itself at the origin of the creation or deletion of a methylation site and deeply impacted the methylation level of neighboring methylation sites. Cd-induced epigenetic changes were associated with different mRNA transcripts and an improved condition of Cd fish. Our results emphasize a tight relationship between genetic and epigenetic mechanisms and suggest that their interplay and pre-existing diversity can allow rapid adaptation to anthropogenic environmental changes.

Variations in phenotypic plasticity in a cosmopolitan copepod species across latitudinal hydrographic gradients

Studies assessing latitudinal variations in habitat conditions and phenotypic plasticity among populations yield evidence of the mechanisms governing differentiation in the potential to adapt to current/future habitat changes. The cosmopolitan copepod species Acartia tonsa thrives across ocean clines delimiting Seasonal (30–40° S) and Permanent (10–30° S) Upwelling coastal provinces established during the middle–late Pliocene (3.6–1.8 Ma) alongshore the South East Pacific (SEP), nowadays exhibiting contrasting variability features related to several ocean drivers (temperature, salinity, pH, and food availability). Latitudinal variation across the range of environmental conditions of the coastal provinces can contribute toward shaping divergent A. tonsa’s phenotypes, for example, through specific patterns of phenotypic plasticity in morphological and physiological traits and tolerance to environmental drivers. With the aim of contributing to the understanding of these adaptive processes in a relatively little studied oceanic region, here we compared the expression of parental (i.e., adult size, egg production, and ingestion rate) and offspring (i.e., egg size) traits in relation to variation in environmental habitat conditions across different cohorts of two distant (&gt; 15° latitude) A. tonsa populations inhabiting estuarine and upwelling habitats located in the Seasonal and Permanent Upwelling province, respectively. Mean conditions and ranges of variability in the habitat conditions and phenotypic plasticity of parental and offspring traits within and among cohorts of A. tonsa populations varied significantly across the different examined regions (i.e., Seasonal vs. Permanent). We also found significant differences in the coupling of habitat variability and trait expression, suggesting that the differences in trait expressions might be related to habitat variability. The phenotypic divergence was translated to cohort-related patterns of trait trade-offs regulating reproduction and tolerance of egg production efficiency that can jointly determine the level of plasticity, genetic structure, or local adaptation. The current findings provide novel evidence of how divergent phenotypes might sustain A. tonsa populations across variable coastal provinces of the SEP.

A Transcriptomic Analysis of Phenotypic Plasticity in Crassostrea virginica Larvae under Experimental Acidification

Ocean acidification (OA) is a major threat to marine calcifiers, and little is known regarding acclimation to OA in bivalves. This study combined physiological assays with next-generation sequencing to assess the potential for recovery from and acclimation to OA in the eastern oyster (Crassostrea virginica) and identify molecular mechanisms associated with resilience. In a reciprocal transplant experiment, larvae transplanted from elevated pCO2 (~1400 ppm) to ambient pCO2 (~350 ppm) demonstrated significantly lower mortality and larger size post-transplant than oysters remaining under elevated pCO2 and had similar mortality compared to those remaining in ambient conditions. The recovery after transplantation to ambient conditions demonstrates the ability for larvae to rebound and suggests phenotypic plasticity and acclimation. Transcriptomic analysis supported this hypothesis as genes were differentially regulated under OA stress. Transcriptomic profiles of transplanted and non-transplanted larvae terminating in the same final pCO2 converged, further supporting the idea that acclimation underlies resilience. The functions of differentially expressed genes included cell differentiation, development, biomineralization, ion exchange, and immunity. Results suggest acclimation as a mode of resilience to OA. In addition, the identification of genes associated with resilience can serve as a valuable resource for the aquaculture industry, as these could enable marker-assisted selection of OA-resilient stocks.

Molecular basis of ocean acidification sensitivity and adaptation in Mytilus galloprovincialis

Predicting the potential for species adaption to climate change is challenged by the need to identify the physiological mechanisms that underpin species vulnerability. Here, we investigated the sensitivity to ocean acidification in marine mussels during early development, and specifically the trochophore stage. Using RNA and DNA sequencing and in situ RNA hybridization, we identified developmental processes associated with abnormal development and rapid adaptation to low pH. Trochophores exposed to low pH seawater exhibited 43 differentially expressed genes. Gene annotation and in situ hybridization of differentially expressed genes point to pH sensitivity of (1) shell field development and (2) cellular stress response. Five genes within these two processes exhibited shifts in allele frequencies indicative of a potential for rapid adaptation. This case study contributes direct evidence that protecting species’ existing genetic diversity is a critical management action to facilitate species resilience to climate change.

•

Marine mussel larval development and genetic adaptation in low pH seawater

•

RNA and DNA responses reveal impacts on shell field development and cell stress

•

Five genes exhibited both physiological sensitivity and long-term adaptive potential

•

Conserving standing genetic variation could bolster resilience to global change

Acclimatory gene expression of primed clams enhances robustness to elevated pCO2

Sub-lethal exposure to environmental challenges may enhance ability to cope with chronic or repeated change, a process known as priming. In a previous study, pre-exposure to seawater enriched with pCO₂ improved growth and reduced antioxidant capacity of juvenile Pacific geoduck Panopea generosa, suggesting that transcriptional shifts may drive phenotypic modifications post-priming. To this end, juvenile clams were sampled and TagSeq gene expression data analyzed after 1) a 110-day acclimation under ambient (921 μatm, naïve) and moderately-elevated pCO₂ (2870 μatm, pre-exposed); then following 2) a second 7-day exposure to three pCO₂ treatments (ambient: 754 μatm; moderately-elevated: 2750 μatm; severely-elevated: 4940 μatm), a 7-day return to ambient pCO₂ , and a third 7-day exposure to two pCO₂ treatments (ambient: 967 μatm; moderately-elevated: 3030 μatm). Pre-exposed geoducks frontloaded genes for stress and apoptosis/innate immune response, homeostatic processes, protein degradation, and transcriptional modifiers. Pre-exposed geoducks were also responsive to subsequent encounters, with gene sets enriched for mitochondrial recycling and immune defense under elevated pCO₂ and energy metabolism and biosynthesis under ambient recovery. In contrast, gene sets with higher expression in naïve clams were enriched for fatty-acid degradation and glutathione components, suggesting naïve clams could be depleting endogenous fuels, with unsustainable energetic requirements if changes in carbonate chemistry persist. Collectively, our transcriptomic data indicates pCO₂ priming during post-larval periods could, via gene expression regulation, enhance robustness in bivalves to environmental change. Such priming approaches may be beneficial for aquaculture, as seafood demand intensifies concurrent with increasing climate change in marine systems.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Genome-wide signatures of synergistic epistasis during parallel adaptation in a Baltic Sea copepod

The role of epistasis in driving adaptation has remained an unresolved problem dating back to the Evolutionary Synthesis. In particular, whether epistatic interactions among genes could promote parallel evolution remains unexplored. To address this problem, we employ an Evolve and Resequence (E&R) experiment, using the copepod Eurytemora affinis, to elucidate the evolutionary genomic response to rapid salinity decline. Rapid declines in coastal salinity at high latitudes are a predicted consequence of global climate change. Based on time-resolved pooled whole-genome sequencing, we uncover a remarkably parallel, polygenic response across ten replicate selection lines, with 79.4% of selected alleles shared between lines by the tenth generation of natural selection. Using extensive computer simulations of our experiment conditions, we find that this polygenic parallelism is consistent with positive synergistic epistasis among alleles, far more so than other mechanisms tested. Our study provides experimental and theoretical support for a novel mechanism promoting repeatable polygenic adaptation, a phenomenon that may be common for selection on complex physiological traits.

WOS, Steinberg DK, Swadling KM, Tarling GA, Thorpe SE, Veytia D, Ward P, Weldrick CK, Yang G. Status, Change, and Futures of Zooplankton in the Southern Ocean

In the Southern Ocean, several zooplankton taxonomic groups, euphausiids, copepods, salps and pteropods, are notable because of their biomass and abundance and their roles in maintaining food webs and ecosystem structure and function, including the provision of globally important ecosystem services. These groups are consumers of microbes, primary and secondary producers, and are prey for fishes, cephalopods, seabirds, and marine mammals. In providing the link between microbes, primary production, and higher trophic levels these taxa influence energy flows, biological production and biomass, biogeochemical cycles, carbon flux and food web interactions thereby modulating the structure and functioning of ecosystems. Additionally, Antarctic krill (Euphausia superba) and various fish species are harvested by international fisheries. Global and local drivers of change are expected to affect the dynamics of key zooplankton species, which may have potentially profound and wide-ranging implications for Southern Ocean ecosystems and the services they provide. Here we assess the current understanding of the dominant metazoan zooplankton within the Southern Ocean, including Antarctic krill and other key euphausiid, copepod, salp and pteropod species. We provide a systematic overview of observed and potential future responses of these taxa to a changing Southern Ocean and the functional relationships by which drivers may impact them. To support future ecosystem assessments and conservation and management strategies, we also identify priorities for Southern Ocean zooplankton research.