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

Genetic Determination of a Cryptic Species in the Littoraria Genus With Whole-Genome Molecular Resolution

Recognizing cryptic species is crucial for understanding global biodiversity. The intertidal snail Littoraria flammea is potentially a cryptic species of L. melanostoma widely distributed in the Northwest Pacific. However, the evidence from traditional morphology and single genetic markers is inconsistent. Our study combined quantitative morphological and whole-genome molecular data to clarify the phylogenetic relationship of three species (L. flammea, L. aff. melanostoma, and L. melanostoma). Three-dimensional models of shells revealed significant differences in morphology between L. flammea and L. melanostoma. Neutral SNPs indicated that individuals of L. flammea and L. melanostoma were in different clusters. The ratio of interspecific F ST to intraspecific F ST between L. flammea and L. melanostoma (16) was much larger than the lowest ratio (2.31) in six published genera with cryptic species in gastropods. Non-neutral SNPs disclosed divergence in functional genes related to reproduction and protein binding. The morphological and phylogenetic analyses corroborated the transitional status of L. aff. melanostoma. These results confirmed that the L. flammea snails north of the Yangtze River Estuary is a cryptic species of L. melanostoma, and allopatric speciation occurs in the L. melanostoma complex.

Pushed waves, trailing edges, and extreme events: Eco-evolutionary dynamics of a geographic range shift in the owl limpet, Lottia gigantea

As climatic variation re-shapes global biodiversity, understanding eco-evolutionary feedbacks during species range shifts is of increasing importance. Theory on range expansions distinguishes between two different forms: "pulled" and "pushed" waves. Pulled waves occur when the source of the expansion comes from low-density peripheral populations, while pushed waves occur when recruitment to the expanding edge is supplied by high-density populations closer to the species' core. How extreme events shape pushed/pulled wave expansion events, as well as trailing-edge declines/contractions, remains largely unexplored. We examined eco-evolutionary responses of a marine invertebrate (the owl limpet, Lottia gigantea) that increased in abundance during the 2014-2016 marine heatwaves near the poleward edge of its geographic range in the northeastern Pacific. We used whole-genome sequencing from 19 populations across >11 degrees of latitude to characterize genomic variation, gene flow, and demographic histories across the species' range. We estimated present-day dispersal potential and past climatic stability to identify how contemporary and historical seascape features shape genomic characteristics. Consistent with expectations of a pushed wave, we found little genomic differentiation between core and leading-edge populations, and higher genomic diversity at range edges. A large and well-mixed population in the northern edge of the species' range is likely a result of ocean current anomalies increasing larval settlement and high-dispersal potential across biogeographic boundaries. Trailing-edge populations have higher differentiation from core populations, possibly driven by local selection and limited gene flow, as well as high genomic diversity likely as a result of climatic stability during the Last Glacial Maximum. Our findings suggest that extreme events can drive poleward range expansions that carry the adaptive potential of core populations, while also cautioning that trailing-edge extirpations may threaten unique evolutionary variation. This work highlights the importance of understanding how both trailing and leading edges respond to global change and extreme events.

Flight or fight: different strategies of intertidal periwinkle Littoraria sinensis coping with high temperature across populations

Intertidal organisms usually live near their upper thermal limits, and are vulnerable to future global warming. As a vital response to thermal stress, thermoregulatory strategy in physiological and behavioral performance is essential for organisms coping with thermal stress and surviving the changing world. To investigate the relationship between the thermoregulatory strategy and habitat temperature, in the present study, we comparatively investigated the thermal responsive strategy among different geographic populations of the supralittoral snail Littoraria sinensis by determining snails' cardiac function and behavioral performance. Our results indicated that populations inhabiting high ambient temperatures had higher sublethal temperatures (i.e. Arrhenius breakpoint temperatures, ABTs, the temperature at which the heart rate shapely decreases with further heating) and lethal temperatures (i.e. Flatline temperatures, FLTs, the temperature at which heart rate ceases), and behaved less actively (e.g. shorter moving distances and shorter moving time) in the face of high and rising temperatures-a physiological fight strategy. On the other hand, populations at relatively low ambient temperatures had relatively lower physiological upper thermal limits with lower ABTs and FLTs and moved more actively in the face of high and rising temperatures-a behavioral flight strategy. These results demonstrate that the thermoregulatory strategies of the snails are closely related to their habitat temperatures and are different among populations surviving divergent thermal environments.

Genome-wide sequencing reveals geographical variations in the thermal adaptation of an aquaculture species with frequent seedling introductions

Climate change and human activity are essential factors affecting marine biodiversity and aquaculture, and understanding the impacts of human activities on the genetic structure to increasing high temperatures is crucial for sustainable aquaculture and marine biodiversity conservation. As a commercially important bivalve, the Manila clam Ruditapes philippinarum is widely distributed along the coast of China, and it has been frequently introduced from Fujian Province, China, to other regions for aquaculture. In this study, we collected four populations of Manila clams from different areas to evaluate their thermal tolerance by measuring cardiac performance and genetic variations using whole-genome resequencing. The upper thermal limits of the clams showed high variations within and among populations. Different populations displayed divergent genetic compositions, and the admixed population was partly derived from the Zhangzhou population in Fujian Province, implying a complex genomic landscape under the influence of local genetic sources and human introductions. Multiple single nucleotide polymorphisms (SNPs) were associated with the cardiac functional traits, and some of these SNPs can affect the codon usage and the structural stability of the resulting protein. This study shed light on the importance of establishing long-term ecological and genetic monitoring programs at the local level to enhance resilience to future climate change.

Analysis of Genetic Diversity and Structure of Eight Populations of Nerita yoldii along the Coast of China Based on Mitochondrial COI Gene

Nerita yoldii is a euryhaline species commonly found in the intertidal zone. To investigate the genetic diversity of 233 N. yoldii individuals from eight locations along the coast of China, we utilized the mitochondrial COI gene as a molecular marker. A total of 34 haplotypes were detected, exhibiting a mean haplotype diversity (Hd) of 0.5915 and a mean nucleotide diversity (Pi) of 0.0025, indicating high levels of genetic diversity among all populations. An analysis of molecular variance (AMOVA) indicated that the primary source of genetic variation occurs within populations. In addition, neutral tests and mismatch analyses suggested that N. yoldii populations may have experienced bottleneck events. Moderate genetic differentiation was observed between Xiapu and other populations, excluding the Taizhou population, and may be attributed to the ocean currents. Intensively studying the genetic variation and population structure of N. yoldii populations contributes to understanding the current population genetics of N. yoldii in the coastal regions of China. This not only provides a reference for the study of other organisms in the same region but also lays the foundation for the systematic evolution of the Neritidae family.

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.

Roles of multi-level temperature-adaptive responses and microhabitat variation in establishing distributions of intertidal species

How intertidal species survive their harsh environment and how best to evaluate and forecast range shifts in species distribution are two important and closely related questions for intertidal ecologists and global change biologists. Adaptive variation in responses of organisms to environmental change across all levels of biological organization - from behavior to molecular systems - is of key importance in setting distribution patterns, yet studies often neglect the interactions of diverse types of biological variation (e.g. differences in thermal optima owing to genetic and acclimation-induced effects) with environmental variation, notably at the scale of microhabitats. Intertidal species have to cope with extreme and frequently changing thermal stress, and have shown high variation in thermal sensitivities and adaptive responses at different levels of biological organization. Here, I review the physiological and biochemical adaptations of intertidal species to environmental temperature on multiple spatial and temporal scales. With fine-scale datasets for the thermal limits of individuals and for environmental temperature variation at the microhabitat scale, we can map the thermal sensitivity for each individual in different microhabitats, and then scale up the thermal sensitivity analysis to the population level and, finally, to the species level by incorporating physiological traits into species distribution models. These more refined mechanistic models that include consideration of physiological variations have higher predictive power than models that neglect these variations, and they will be crucial to answering the questions posed above concerning adaptive mechanisms and the roles they play in governing distribution patterns in a rapidly changing world.

Evaluation of species thermal sensitivity with individual-based physiological performance

Ignoring intraspecific variations can prevent us from accurately assessing species' thermal sensitivity to global warming. Individual-based physiological performance provides a feasible solution to depict species' thermal sensitivity using a bottom-up approach. We measured the cardiac performance of intertidal bivalves (1159 individuals from multiple populations of six bivalves), determined the upper thermal limit of each individual, calculated the proportions of individuals suffering sublethal/lethal heat stress, and mapped sensitive regions to high temperatures. Results showed that high inter-individual variations of physiological performance existed in levels of populations and species, and species' thermal sensitivity was positively related to the intraspecific variations of heat tolerance. This bottom-up approach scaled up from individual, population to species emphasizes the importance of individual-based physiology performance in assessing thermal sensitivity across different hierarchical levels and enables better evaluating and forecasting of species responses to global warming.

Coping with harsh heat environments: molecular adaptation of metabolic depression in the intertidal snail Echinolittorina radiata

Harsh thermal environments in the rocky intertidal zone pose serious physiological and molecular challenges to the inhabitants. Metabolic depression is regarded as an energy-conserving feature of intertidal species. To understand the molecular mechanism of metabolic depression, we investigated physiological and transcriptomic responses in the intertidal snail Echinolittorina radiata. The metabolic rate and expression of most genes were insensitive to temperatures ranging from 33 to 45 °C and then increased with further heating to 52 °C. Different from other genes, the genes involved in heat shock response (HSR) and oxidative stress response (OSR) (e.g., genes encoding heat shock protein 70 (HSP70) and cytochrome P450 protein (CYP450)) kept upregulating during metabolic depression. These high levels of HSR and OSR genes should be important for surviving the harsh thermal environments on the rocky shore. In the population experiencing more frequent moderate heat events, the depression breadth was larger, and the change in magnitude of upregulation was insensitive for HSR genes (e.g., HSP70s) but heat-sensitive for OSR genes (e.g., CYP450s) at the temperature of 37 to 45 °C. These findings indicate that both the thermal sensitivity of HSR and OSR genes and the insensitivity of metabolic genes are crucial for surviving extreme intertidal environments, and different populations of the same species rely on various physiological mechanisms to differing extents to deal with heat stress. The cellular stress response is not a "one size fits all" response across populations largely depending on local thermal regimes.

Single Nucleotide Polymorphism Markers with Applications in Conservation and Exploitation of Aquatic Natural Populations

An increasing number of aquatic species have been studied for genetic polymorphism, which extends the knowledge on their natural populations. One type of high-resolution molecular marker suitable for studying the genetic diversity of large numbers of individuals is single nucleotide polymorphism (SNP). This review is an attempt to show the range of applications of SNPs in studies of natural populations of aquatic animals. In recent years, SNPs have been used in the genetic analysis of wild and enhanced fish and invertebrate populations in natural habitats, exploited migratory species in the oceans, migratory anadromous and freshwater fish and demersal species. SNPs have been used for the identification of species and their hybrids in natural environments, to study the genetic consequences of restocking for conservation purposes and the negative effects on natural populations of fish accidentally escaping from culture. SNPs are very useful for identifying genomic regions correlated with phenotypic variants relevant for wildlife protection, management and aquaculture. Experimental size-selective catches of populations created in tanks have caused evolutionary changes in life cycles of fishes. The research results have been discussed to clarify whether the fish populations in natural conditions can undergo changes due to selective harvesting targeting the fastest-growing fishes.

Thermal heterogeneity is an important factor for maintaining the genetic differentiation pattern of the pelagic barnacle Lepas anatifera in the northwest Pacific

Macrobenthic invertebrates are ubiquitously distributed in the epipelagic zone of the open ocean. Yet, our understanding of their genetic structure patterns remains poorly understood. Investigating the genetic differentiation patterns of pelagic Lepas anatifera and clarifying the potential roles of temperature maintaining this pattern are crucial for our understanding of the distribution and biodiversity of pelagic macrobenthos. In the present study, mitochondrial cytochrome oxidase subunit I (mtDNA COI) from three South China Sea (SCS) populations and six Kuroshio Extension (KE) region populations of L. anatifera sampled from fixed buoys and genome-wide SNPs from a subset of populations (two SCS populations and four KE region populations) were sequenced and analyzed for investigating the genetic pattern of the pelagic barnacle. Water temperature was different among sampling sites; in other words, the water temperature decreased with latitude increases, and the water temperature on the surface was higher than in the subsurface. Our result showed that three lineages with clear genetic differentiation were found in different geographical locations and depths based on mtDNA COI, all SNPs, neutral SNPs, and outlier SNPs. Lineage 1 and lineage 2 were dominant in the subsurface populations and surface populations from the KE region, respectively. Lineage 3 was dominant in the SCS populations. Historical events during the Pliocene epoch shaped the differentiation of the three lineages, while, nowadays, temperature heterogeneity maintains the current genetic pattern of L. anatifera in the northwest Pacific. The subsurface populations were genetically isolated from the surface populations in the Kuroshio Extension (KE) region, implying small-scale vertical thermal heterogeneity was also an important factor maintaining the genetic differentiation pattern of the pelagic species.