Morphological and molecular characterization of Stomachicola muraenesocis Yamaguti, 1934 (Digenea: Hemiuridae) from the daggertooth pike conger Muraenesox cinereus (Forsskål).

Hemiurid digeneans conspecific with Stomachicola muraenesocis Yamaguti, 1934 (the type species of the genus Stomachicola Yamaguti, 1934) were collected from the stomach of the daggertooth pike conger Muraenesox cinereus (Forsskål) off the Persian Gulf of Iran. This study aimed to provide a detailed characterization of Stom. muraenesocis, including measurements, illustrations and scanning electron microscopy (s.e.m.) representations. Comparisons with the original and previous descriptions revealed morphological and metrical variations in several features (i.e. body size and shape, arrangement of reproductive organs, soma to ecsoma length ratio, position of genital opening, number of vitelline tubules and extension of uterine coils) between Stom. muraenesocis from different hosts and localities. This study presents the first molecular sequence data associated with the small (18S) and large (28S) subunit nuclear ribosomal RNA genes (rDNA) for Stom. muraenesocis. Phylogenetic analyses of the 18S dataset placed Stom. muraenesocis as sister lineage to a clade formed of a group of species of Lecithaster Lühe, 1901 (Lecithasteridae Odhner, 1905). In contrast, phylogenetic analyses based on the 28S consistently recovered a sister relationship between Stom. muraenesocis and representatives of the Hemiuridae Looss, 1899. Further comprehensive phylogenetically based classification in light of morphology and taxonomic history of the Hemiuridae and Lecithasteridae is required to infer phylogenetic affinities and historical biogeography of Stomachicola. A comprehensive list of previously reported species of Stomachicola together with their associated hosts, localities and morphometric data is provided.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments.

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

Fate of MHCII in salmonids following 4WGD.

Major histocompatibility complex (MHC) genes are key players in the adaptive immunity providing a defense against invading pathogens. Although the basic structures are similar when comparing mammalian and teleost MHC class II (MHCII) molecules, there are also clear-cut differences. Based on structural requirements, the teleosts non-classical MHCII molecules do not comply with a function similar to the human HLA-DM and HLA-DO, i.e., assisting in peptide loading and editing of classical MHCII molecules. We have previously studied the evolution of teleost class II genes identifying various lineages and tracing their phylogenetic occurrence back to ancient ray-finned fishes. We found no syntenic MHCII regions shared between cyprinids, salmonids, and neoteleosts, suggesting regional instabilities. Salmonids have experienced a unique whole genome duplication 94 million years ago, providing them with the opportunity to experiment with gene duplicates. Many salmonid genomes have recently become available, and here we set out to investigate how MHCII has evolved in salmonids using Northern pike as a diploid sister phyla, that split from the salmonid lineage prior to the fourth whole genome duplication (4WGD) event. We identified 120 MHCII genes in pike and salmonids, ranging from 11 to 20 genes per species analyzed where DB-group genes had the most expansions. Comparing the MHC of Northern pike with that of Atlantic salmon and other salmonids species provides a tale of gene loss, translocations, and genome rearrangements.

Environmental DNA detection of endangered and invasive species in Kejimkujik National Park and Historic Site.

The use of environmental DNA (eDNA) allows the early detection of aquatic species at low densities (e.g., elusive and invasive species), which otherwise could be challenging to monitor using conventional techniques. Here, we assess the ability of eDNA sampling to detect the presence or absence of one species at risk (Blanding's turtle) and two invasive species (chain pickerel and smallmouth bass) in Kejimkujik National Park and National Historic Site, Nova Scotia, where the aquatic system is highly acidic and rich in organic compounds. Five replicates of 1 L water samples were taken per sampling site. Water filtration and eDNA extractions were performed on-site, while qPCR reactions were performed in the laboratory using species-specific assays. Samples were treated with an inhibition removal kit and analyzed pre- and post-inhibition removal. Despite the low pH and PCR inhibitors in water samples, our results showed positive eDNA detections in almost all expected positive sites (except in one site for Blanding's turtle). Detections of the target species were also observed at sites where their presence was previously unknown. Our study supports the advantage of eDNA to monitor species at low densities, revealing new distributions or recently invaded areas. We also demonstrate how eDNA can directly instruct management strategies in Kejimkujik.

Catch-related and genetic outcome of adult northern pike Esox lucius stocking in a large river system.

Genetic introgression from stocked adult northern pike Esox lucius to a wild self-recruiting population was detected in a large river system and some stocked E. lucius survived up to two spawning seasons and dispersed over several kilometres in the river. Moreover, the catch rate of stocked E. lucius by anglers was low (9.6%), hence suggesting that the efficiency of stocking activity is questionable.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.).

BACKGROUND: Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes. RESULTS: The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30-38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10-30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides. CONCLUSIONS: Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

The axon degeneration gene SARM1 is evolutionarily distinct from other TIR domain-containing proteins.

Many forms of neurodegenerative disease are characterized by Wallerian degeneration, an active program of axonal destruction. Recently, the important player which enacts Wallerian degeneration was discovered, the multidomain protein SARM1. Since the SARM1 protein has classically been thought of as an innate immune molecule, its role in Wallerian degeneration has raised questions on the evolutionary forces acting on it. Here, we synthesize a picture of SARM1's evolution through various organisms by examining the molecular and genetic changes of SARM1 and the genes around it. Using proteins that possess domains homologous to SARM1, we established distances and Ka/Ks values through 5671 pairwise species-species comparisons. We demonstrate that SARM1 diverged across species in a pattern similar to other SAM domain-containing proteins. This is surprising, because it was expected that SARM1 would behave more like its TIR domain relatives. Going along with this divorce from TIR, we also noted that SARM1's TIR is under stronger purifying selection than the rest of the TIR domain-containing proteins (remaining highly conserved). In addition, SARM1's synteny analysis reveals that the surrounding gene cluster is highly conserved, functioning as a potential nexus of gene functionality across species. Taken together, SARM1 demonstrates a unique evolutionary pattern, separate from the TIR domain protein family.

Species richness in the African pike genus Hepsetus: a perfect match between genetics and morphology.

In this study, morphological differences were found that corroborate earlier results that showed the existence of six species within the previously monospecific African pike genus Hepsetus. Additional genetic data (coI, mtDNA and rag1, nDNA) confirm the morphology-based species delineations. Deep genetic divergences imply a relatively old age for diversification within the genus. An identification key for the six species is provided in the present study.

Individual heterogeneity in life histories and eco-evolutionary dynamics.

Individual heterogeneity in life history shapes eco-evolutionary processes, and unobserved heterogeneity can affect demographic outputs characterising life history and population dynamical properties. Demographic frameworks like matrix models or integral projection models represent powerful approaches to disentangle mechanisms linking individual life histories and population-level processes. Recent developments have provided important steps towards their application to study eco-evolutionary dynamics, but so far individual heterogeneity has largely been ignored. Here, we present a general demographic framework that incorporates individual heterogeneity in a flexible way, by separating static and dynamic traits (discrete or continuous). First, we apply the framework to derive the consequences of ignoring heterogeneity for a range of widely used demographic outputs. A general conclusion is that besides the long-term growth rate lambda, all parameters can be affected. Second, we discuss how the framework can help advance current demographic models of eco-evolutionary dynamics, by incorporating individual heterogeneity. For both applications numerical examples are provided, including an empirical example for pike. For instance, we demonstrate that predicted demographic responses to climate warming can be reversed by increased heritability. We discuss how applications of this demographic framework incorporating individual heterogeneity can help answer key biological questions that require a detailed understanding of eco-evolutionary dynamics.

Evolutionary divergence of adult body size and juvenile growth in sympatric subpopulations of a top predator in aquatic ecosystems.

Evolutionary theory predicts that different selective regimes may contribute to divergent evolution of body size and growth rate among populations, but most studies have focused on allopatric populations. Here, we studied five sympatric subpopulations of anadromous northern pike (Esox lucius) in the Baltic Sea subjected to allopatric habitats for a short period of their life cycle due to homing behavior. We report differences in adult body size among subpopulations that were in part due to variation in growth rate. Body size of emigrating juveniles also differed among subpopulations, and differences remained when individuals were reared in a common environment, thus indicating evolutionary divergence among subpopulations. Furthermore, a QST-FST comparison indicated that differences had evolved due to divergent selection rather than genetic drift, possibly in response to differences in selective mortality among spawning habitats during the allopatric life stage. Adult and juvenile size were negatively correlated across subpopulations, and reconstruction of growth trajectories of adult fishes suggested that body size differences developed gradually and became accentuated throughout the first years of life. These results represent rare evidence that sympatric subpopulations can evolve differences in key life-history traits despite being subjected to allopatric habitats during only a very short fraction of their life.

Loss of genetic variation and ancestral sex determination system in North American northern pike characterized by whole-genome resequencing.

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = -1.07 and -1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby, the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby(-) males in Alaska and amhby(+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby).

Parvalbumin Gene: A Valuable Marker for Pike Authentication and Allergen Risk Assessment.

Fish from the pike (Esox) genus are valued in gastronomy for their superior meat quality. However, they can cause allergic reactions in sensitive consumers. This work aimed to fill the gap in the detection of pike allergens using molecular-biological techniques. New, fast, and accurate loop-mediated isothermal amplification (LAMP) and real-time PCR (qPCR) assays were designed to detect pike DNA using the parvalbumin gene as a marker. LAMP was assessed by electrophoresis, SYBR green optical detection, and real-time fluorescence detection. The latter was the most sensitive, detecting as little as 0.78 ng of pike DNA; the qPCR detection limit was 0.1 ng. The LAMP analysis took 20-70 min, which is significantly faster than qPCR. The study provides reliable detection and quantification of the parvalbumin gene in both fresh and processed samples and further highlights the versatility of the use of the parvalbumin gene for the authentication of food products and consumer protection via refined allergen risk assessment that is independent of the type of tissue or food processing method used.

An evolutionary perspective on Elovl5 fatty acid elongase: comparison of Northern pike and duplicated paralogs from Atlantic salmon.

BACKGROUND: The ability to produce physiologically critical LC-PUFA from dietary fatty acids differs greatly among teleost species, and is dependent on the possession and expression of fatty acyl desaturase and elongase genes. Atlantic salmon, as a result of a recently duplicated genome, have more of these enzymes than other fish. Recent phylogenetic studies show that Northern pike represents the closest extant relative of the preduplicated ancestral salmonid. Here we characterise a pike fatty acyl elongase, elovl5, and compare it to Atlantic salmon elovl5a and elovl5b duplicates. RESULTS: Phylogenetic analyses show that Atlantic salmon paralogs are evolving symmetrically, and they have been retained in the genome by purifying selection. Heterologous expression in yeast showed that Northern pike Elovl5 activity is indistinguishable from that of the salmon paralogs, efficiently elongating C18 and C20 substrates. However, in contrast to salmon, pike elovl5 was predominantly expressed in brain with negligible expression in liver and intestine. CONCLUSIONS: We suggest that the predominant expression of Elovl5b in salmon liver and Elovl5a in salmon intestine is an adaptation, enabled by genome duplication, to a diet rich in terrestrial invertebrates which are relatively poor in LC-PUFA. Pike have retained an ancestral expression profile which supports the maintenance of PUFA in the brain but, due to a highly piscivorous LC-PUFA-rich diet, is not required in liver and intestine. Thus, the characterisation of elovl5 in Northern pike provides insights into the evolutionary divergence of duplicated genes, and the ecological adaptations of salmonids which have enabled colonisation of nutrient poor freshwaters.

Loss of gut microbial diversity in the cultured, agastric fish, Mexican pike silverside (Chirostoma estor: Atherinopsidae).

Teleost fish are the most diverse group of extant vertebrates and have varied digestive anatomical structures and strategies, suggesting they also possess an array of different host-microbiota interactions. Differences in fish gut microbiota have been shown to affect host development, the process of gut colonization, and the outcomes of gene-environment or immune system-microbiota interactions. There is generally a lack of studies on the digestive mechanisms and microbiota of agastric short-intestine fish however, meaning that we do not understand how changes in gut microbial diversity might influence the health of these types of fish. To help fill these gaps in knowledge, we decided to study the Mexican pike silverside (Chirostoma estor) which has a simplified alimentary canal (agastric, short-intestine, 0.7 gut relative length) to observe the diversity and metabolic potential of its intestinal microbiota. We characterized gut microbial populations using high-throughput sequencing of the V3 region in bacterial 16S rRNA genes while searching for population shifts resulting associated with fish development in different environments and cultivation methods. Microbiota samples were taken from the digesta, anterior and posterior intestine (the three different intestinal components) of fish that grew wild in a lake, that were cultivated in indoor tanks, or that were raised in outdoor ponds. Gut microbial diversity was significantly higher in wild fish than in cultivated fish, suggesting a loss of diversity when fish are raised in controlled environments. The most abundant phyla observed in these experiments were Firmicutes and Proteobacteria, particularly of the genera Mycoplasma, Staphylococcus, Spiroplasma, and Aeromonas. Of the 14,161 OTUs observed in this experiment, 133 were found in all groups, and 17 of these, belonging to Acinetobacter, Aeromonas, Pseudomonas, and Spiroplasma genera, were found in all samples suggesting the existence of a core C. estor microbiome. Functional metagenomic prediction of bacterial ecological functions using PICRUSt2 suggested that different intestinal components select for functionally distinct microbial populations with variation in pathways related to the metabolism of amino acids, vitamins, cofactors, and energy. Our results provide, for the first time, information on the bacterial populations present in an agastric, short-gut teleost with commercial potential and show that controlled cultivation of this fish reduces the diversity of its intestinal microbiota.

The rise and fall of the ancient northern pike master sex-determining gene.

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome.

BACKGROUND: Salmonids are one of the most intensely studied fish, in part due to their economic and environmental importance, and in part due to a recent whole genome duplication in the common ancestor of salmonids. This duplication greatly impacts species diversification, functional specialization, and adaptation. Extensive new genomic resources have recently become available for Atlantic salmon (Salmo salar), but documentation of allelic versus duplicate reference genes remains a major uncertainty in the complete characterization of its genome and its evolution. RESULTS: From existing expressed sequence tag (EST) resources and three new full-length cDNA libraries, 9,057 reference quality full-length gene insert clones were identified for Atlantic salmon. A further 1,365 reference full-length clones were annotated from 29,221 northern pike (Esox lucius) ESTs. Pairwise dN/dS comparisons within each of 408 sets of duplicated salmon genes using northern pike as a diploid out-group show asymmetric relaxation of selection on salmon duplicates. CONCLUSIONS: 9,057 full-length reference genes were characterized in S. salar and can be used to identify alleles and gene family members. Comparisons of duplicated genes show that while purifying selection is the predominant force acting on both duplicates, consistent with retention of functionality in both copies, some relaxation of pressure on gene duplicates can be identified. In addition, there is evidence that evolution has acted asymmetrically on paralogs, allowing one of the pair to diverge at a faster rate.

Environmental exposure of northern pike to a primary wastewater effluent: Impact on the lipidomic profile and lipid metabolism.

Lipids play important roles in growth, reproduction, locomotion, and migration of fish. Municipal effluents, which are complex mixtures of biological and chemical compounds including flame retardants, have been shown to alter lipid metabletabolism in environmentally and experimentally exposed fish. Down-regulation of several genes coding for fatty acid metabolism enzymes has previously been reported in male northern pike (Esox lucius) collected in the St. Lawrence River (QC, Canada) downstream of a major primary wastewater treatment plant (WWTP) point of discharge. The main objective of this study was to investigate the effects of exposure to the Montreal's WWTP effluent on the lipidomic profile (i.e., fatty acids, acylcarnitines, and phospholipids) as well as the transcription of genes related to lipid metabolism in the liver of northern pike collected upstream and downstream of this WWTP effluent. Halogenated flame retardant concentrations were also determined in pike liver and used as markers of exposure to this effluent. Greater concentrations of saturated and monounsaturated lysophosphatidylcholines (LPCs) and lower concentrations of polyunsaturated LPCs were determined in the liver of pike collected downstream of the WWTP compared to those collected upstream. Lower mRNA levels of peroxisome proliferator-activated receptor alpha (pparα), a major regulator of lipid metabolism, were also measured in pike exposed to Montreal's WWTP effluent. In addition, the relative contributions (%) of LPC 18:2 and LPC14:0, compounds used as markers of obesity and inflammation, were significantly correlated with halogenated flame retardant concentrations and fish girth. Results of the present study suggest that chronic environmental exposure to a primary WWTP effluent can modulate the transcription of genes related to lipid metabolism, and hence affect the hepatic phospholipid composition of pike from the St. Lawrence River.

Harvest-induced disruptive selection increases variance in fitness-related traits.

The form of Darwinian selection has important ecological and management implications. Negative effects of harvesting are often ascribed to size truncation (i.e. strictly directional selection against large individuals) and resultant decrease in trait variability, which depresses capacity to buffer environmental change, hinders evolutionary rebound and ultimately impairs population recovery. However, the exact form of harvest-induced selection is generally unknown and the effects of harvest on trait variability remain unexplored. Here we use unique data from the Windermere (UK) long-term ecological experiment to show in a top predator (pike, Esox lucius) that the fishery does not induce size truncation but disruptive (diversifying) selection, and does not decrease but rather increases variability in pike somatic growth rate and size at age. This result is supported by complementary modelling approaches removing the effects of catch selectivity, selection prior to the catch and environmental variation. Therefore, fishing most likely increased genetic variability for somatic growth in pike and presumably favoured an observed rapid evolutionary rebound after fishery relaxation. Inference about the mechanisms through which harvesting negatively affects population numbers and recovery should systematically be based on a measure of the exact form of selection. From a management perspective, disruptive harvesting necessitates combining a preservation of large individuals with moderate exploitation rates, and thus provides a comprehensive tool for sustainable exploitation of natural resources.

Four decades of opposing natural and human-induced artificial selection acting on Windermere pike (Esox lucius).

The ability of natural selection to drive local adaptation has been appreciated ever since Darwin. Whether human impacts can impede the adaptive process has received less attention. We tested this hypothesis by quantifying natural selection and harvest selection acting on a freshwater fish (pike) over four decades. Across the time series, directional natural selection tended to favour large individuals whereas the fishery targeted large individuals. Moreover, non-linear natural selection tended to favour intermediate sized fish whereas the fishery targeted intermediate sized fish because the smallest and largest individuals were often not captured. Thus, our results unequivocally demonstrate that natural selection and fishery selection often acted in opposite directions within this natural system. Moreover, the two selective factors combined to produce reduced fitness overall and stronger stabilizing selection relative to natural selection acting alone. The long-term ramifications of such human-induced modifications to adaptive landscapes are currently unknown and certainly warrant further investigation.

Causes and consequences of intra-specific variation in vertebral number.

Intraspecific variation in vertebral number is taxonomically widespread. Much scientific attention has been directed towards understanding patterns of variation in vertebral number among individuals and between populations, particularly across large spatial scales and in structured environments. However, the relative role of genes, plasticity, selection, and drift as drivers of individual variation and population differentiation remains unknown for most systems. Here, we report on patterns, causes and consequences of variation in vertebral number among and within sympatric subpopulations of pike (Esox lucius). Vertebral number differed among subpopulations, and common garden experiments indicated that this reflected genetic differences. A QST-FST comparison suggested that population differences represented local adaptations driven by divergent selection. Associations with fitness traits further indicated that vertebral counts were influenced both by stabilizing and directional selection within populations. Overall, our study enhances the understanding of adaptive variation, which is critical for the maintenance of intraspecific diversity and species conservation.