Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera).

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen.

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Natural variation of domestication-related genes contributed to latitudinal expansion and adaptation in soybean.

Soybean is a major source of protein and edible oil worldwide. Originating from the Huang-Huai-Hai region, which has a temperate climate, soybean has adapted to a wide latitudinal gradient across China. However, the genetic mechanisms responsible for the widespread latitudinal adaptation in soybean, as well as the genetic basis, adaptive differentiation, and evolutionary implications of theses natural alleles, are currently lacking in comprehensive understanding. In this study, we examined the genetic variations of fourteen major gene loci controlling flowering and maturity in 103 wild species, 1048 landraces, and 1747 cultivated species. We found that E1, E3, FT2a, J, Tof11, Tof16, and Tof18 were favoured during soybean improvement and selection, which explained 75.5% of the flowering time phenotypic variation. These genetic variation was significantly associated with differences in latitude via the LFMM algorithm. Haplotype network and geographic distribution analysis suggested that gene combinations were associated with flowering time diversity contributed to the expansion of soybean, with more HapA clustering together when soybean moved to latitudes beyond 35°N. The geographical evolution model was developed to accurately predict the suitable planting zone for soybean varieties. Collectively, by integrating knowledge from genomics and haplotype classification, it was revealed that distinct gene combinations improve the adaptation of cultivated soybeans to different latitudes. This study provides insight into the genetic basis underlying the environmental adaptation of soybean accessions, which could contribute to a better understanding of the domestication history of soybean and facilitate soybean climate-smart molecular breeding for various environments.

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.

Genome-wide identification and expression analysis revealed key transcription factors as potential regulators of high-temperature adaptation of Coriolopsis trogii.

Transcription factors (TFs) play a crucial role in gene expression, and studying them can lay the foundation for future research on the functional characterization of TFs involved in various biological processes. In this study, we conducted a genome-wide identification and analysis of TFs in the thermotolerant basidiomycete fungus, Coriolopsis trogii. The TF repertoire of C. trogii consisted of 350 TFs, with C2H2 and Zn2C6 being the largest TF families. When the mycelia of C. trogii were cultured on PDA and transferred from 25 to 35 °C, 14 TFs were up-regulated and 14 TFs were down-regulated. By analyzing RNA-seq data from mycelia cultured at different temperatures and under different carbon sources, we identified 22 TFs that were differentially expressed in more than three comparisons. Co-expression analysis revealed that seven differentially expressed TFs, including four Zn2C6s, one Hap4_Hap_bind, one HMG_box, and one Zinc_knuckle, showed significant correlation with 729 targeted genes. Overall, this study provides a comprehensive characterization of the TF family and systematically screens TFs involved in the high-temperature adaptation of C. trogii, laying the groundwork for further research into the specific roles of TFs in the heat tolerance mechanisms of filamentous fungi.

Discovering novel genomic regions explaining adaptation of bread wheat to conservation agriculture through GWAS.

To sustainably increase wheat yield to meet the growing world population's food demand in the face of climate change, Conservation Agriculture (CA) is a promising approach. Still, there is a lack of genomic studies investigating the genetic basis of crop adaptation to CA. To dissect the genetic architecture of 19 morpho-physiological traits that could be involved in the enhanced adaptation and performance of genotypes under CA, we performed GWAS to identify MTAs under four contrasting production regimes viz., conventional tillage timely sown (CTTS), conservation agriculture timely sown (CATS), conventional tillage late sown (CTLS) and conservation agriculture late sown (CALS) using an association panel of 183 advanced wheat breeding lines along with 5 checks. Traits like Phi2 (Quantum yield of photosystem II; CATS:0.37, CALS: 0.31), RC (Relative chlorophyll content; CATS:55.51, CALS: 54.47) and PS1 (Active photosystem I centers; CATS:2.45, CALS: 2.23) have higher mean values in CA compared to CT under both sowing times. GWAS identified 80 MTAs for the studied traits across four production environments. The phenotypic variation explained (PVE) by these QTNs ranged from 2.15 to 40.22%. Gene annotation provided highly informative SNPs associated with Phi2, NPQ (Quantum yield of non-photochemical quenching), PS1, and RC which were linked with genes that play crucial roles in the physiological adaptation under both CA and CT. A highly significant SNP AX94651261 (9.43% PVE) was identified to be associated with Phi2, while two SNP markers AX94730536 (30.90% PVE) and AX94683305 (16.99% PVE) were associated with NPQ. Identified QTNs upon validation can be used in marker-assisted breeding programs to develop CA adaptive genotypes.

Across two continents: The genomic basis of environmental adaptation in house mice (Mus musculus domesticus) from the Americas.

Replicated clines across environmental gradients can be strong evidence of adaptation. House mice (Mus musculus domesticus) were introduced to the Americas by European colonizers and are now widely distributed from Tierra del Fuego to Alaska. Multiple aspects of climate, such as temperature, vary predictably across latitude in the Americas. Past studies of North American populations across latitudinal gradients provided evidence of environmental adaptation in traits related to body size, metabolism, and behavior and identified candidate genes using selection scans. Here, we investigate genomic signals of environmental adaptation on a second continent, South America, and ask whether there is evidence of parallel adaptation across multiple latitudinal transects in the Americas. We first identified loci across the genome showing signatures of selection related to climatic variation in mice sampled across a latitudinal transect in South America, accounting for neutral population structure. Consistent with previous results, most candidate SNPs were in putatively regulatory regions. Genes that contained the most extreme outliers relate to traits such as body weight or size, metabolism, immunity, fat, eye function, and the cardiovascular system. We then compared these results with the results of analyses of published data from two transects in North America. While most candidate genes were unique to individual transects, we found significant overlap among candidate genes identified independently in the three transects. These genes are diverse, with functions relating to metabolism, immunity, cardiac function, and circadian rhythm, among others. We also found parallel shifts in allele frequency in candidate genes across latitudinal gradients. Finally, combining data from all three transects, we identified several genes associated with variation in body weight. Overall, our results provide strong evidence of shared responses to selection and identify genes that likely underlie recent environmental adaptation in house mice across North and South America.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States.

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

Deep metagenomic sequencing unveils novel SAR202 lineages and their vertical adaptation in the ocean.

SAR202 bacteria in the Chloroflexota phylum are abundant and widely distributed in the ocean. Their genome coding capacities indicate their potential roles in degrading complex and recalcitrant organic compounds in the ocean. However, our understanding of their genomic diversity, vertical distribution, and depth-related metabolisms is still limited by the number of assembled SAR202 genomes. In this study, we apply deep metagenomic sequencing (180 Gb per sample) to investigate microbial communities collected from six representative depths at the Bermuda Atlantic Time Series (BATS) station. We obtain 173 SAR202 metagenome-assembled genomes (MAGs). Intriguingly, 154 new species and 104 new genera are found based on these 173 SAR202 genomes. We add 12 new subgroups to the current SAR202 lineages. The vertical distribution of 20 SAR202 subgroups shows their niche partitioning in the euphotic, mesopelagic, and bathypelagic oceans, respectively. Deep-ocean SAR202 bacteria contain more genes and exhibit more metabolic potential for degrading complex organic substrates than those from the euphotic zone. With deep metagenomic sequencing, we uncover many new lineages of SAR202 bacteria and their potential functions which greatly deepen our understanding of their diversity, vertical profile, and contribution to the ocean's carbon cycling, especially in the deep ocean.

Evaluation of different bermudagrass germplasm at physiological and molecular level under shade along longitudinal and latitudinal gradients.

Responses of turfgrass to shade vary in individual species, and the degree and quality of low light; therefore, the selection of low light tolerant cultivars of turfgrass is important and beneficial for turf management rather than other practices. The stolons of thirteen bermudagrass genotypes were planted with two treatments and three replications of each treatment to establish for one month in the Yangzhou University Jiangsu China greenhouse. The established plants were transferred outside of the greenhouse, and 50% shading was applied to them with a black net. After 30 days of stress treatment, the morpho-physiological and biochemical analyses were performed. The expression of genes such as HEMA, HY5, PIF4, and Cu/ZnSOD was assessed. Cynodon dactylon is a C4, and perennial that grows as lawn grass and is used as forage. Based on different indicator measurements, the most shade-tolerant germplasm was L01 and L06 along the longitudes and L09 and L10 along the latitudes. At the same time, L02 and L08 were more susceptible, respectively. However, germplasm showed greater tolerance in higher latitudes while longitudinal plants showed less stress response. The current study aimed (1) to screen out the most shade-tolerant Cynodon dactylon genotype among 13 along longitudinal and latitudinal gradients in China. (2) to examine morpho-physiological indicators of different bermudagrassgenotypes; (3) to evaluate if and how differences in various indicators of bermudagrass correlated with geographic region. This study will significantly advance the use of Cynodon germplasm in breeding, genomics, management, nomenclature, and phylogeographical study. It will decisively define whether natural selection and migration can drive evolutionary responses for populations to adapt to their new environments effectively.

The loci of environmental adaptation in a model eukaryote.

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species.

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

Pan-Genome Provides Insights into Vibrio Evolution and Adaptation to Deep-Sea Hydrothermal Vents.

This study delves into the genomic features of 10 Vibrio strains collected from deep-sea hydrothermal vents in the Pacific Ocean, providing insights into their evolutionary history and ecological adaptations. Through sequencing and pan-genome analysis involving 141 Vibrio species, we found that deep-sea strains exhibit larger genomes with unique gene distributions, suggesting adaptation to the vent environment. The phylogenomic reconstruction of the investigated isolates revealed the presence of 2 main clades: The first is monophyletic, consisting exclusively of Vibrio alginolyticus, while the second forms a monophyletic clade comprising both Vibrio antiquarius and Vibrio diabolicus species, which were previously isolated from deep-sea vents. All strains carry virulence and antibiotic resistance genes related to those found in human pathogenic Vibrio species which may play a wider ecological role other than host infection in these environments. In addition, functional genomic analysis identified genes potentially related to deep-sea survival and stress response, alongside candidate genes encoding for novel antimicrobial agents. Ultimately, the pan-genome we generated represents a valuable resource for future studies investigating the taxonomy, evolution, and ecology of Vibrio species.

Genomics-based identification of a cold adapted clade in Deinococcus.

BACKGROUND: Microbes in the cold polar and alpine environments play a critical role in feedbacks that amplify the effects of climate change. Defining the cold adapted ecotype is one of the prerequisites for understanding the response of polar and alpine microbes to climate change. RESULTS: Here, we analysed 85 high-quality, de-duplicated genomes of Deinococcus, which can survive in a variety of harsh environments. By leveraging genomic and phenotypic traits with reverse ecology, we defined a cold adapted clade from eight Deinococcus strains isolated from Arctic, Antarctic and high alpine environments. Genome-wide optimization in amino acid composition and regulation and signalling enable the cold adapted clade to produce CO2 from organic matter and boost the bioavailability of mineral nitrogen. CONCLUSIONS: Based primarily on in silico genomic analysis, we defined a potential cold adapted clade in Deinococcus and provided an updated view of the genomic traits and metabolic potential of Deinococcus. Our study would facilitate the understanding of microbial processes in the cold polar and alpine environments.

Innate and germline immune memory: specificity and heritability of the ancient immune mechanisms for adaptation and survival.

The immune memory is one of the defensive strategies developed by both unicellular and multicellular organisms for ensuring their integrity and functionality. While the immune memory of the vertebrate adaptive immune system (based on somatic recombination) is antigen-specific, encompassing the generation of memory T and B cells that only recognize/react to a specific antigen epitope, the capacity of vertebrate innate cells to remember past events is a mostly non-specific mechanism of adaptation. This "innate memory" can be considered as germline-encoded because its effector tools (such as innate receptors) do not need somatic recombination for being active. Also, in several organisms the memory-related information is integrated in the genome of germline cells and can be transmitted to the progeny for several generations, but it can also be erased depending on the environmental conditions. Overall, depending on the organism, its environment and its living habits, innate immune memory appears to be a mechanism for achieving better protection and survival against repeated exposure to microbes/stressful agents present in the same environment or occurring in the same anatomical district, able to adapt to changes in the environmental cues. The anatomical and functional complexity of the organism and its lifespan drive the generation of different immune memory mechanisms, for optimal adaptation to changes in the living/environmental conditions. The concept of innate immunity being non-specific needs to be revisited, as a wealth of evidence suggests a significant degree of specificity both in the primary immune reaction and in the ensuing memory-like responses. This is clearly evident in invertebrate metazoans, in which distinct scenarios can be observed, with both non-specific (immune enhancement) or specific (immune priming) memory-like responses. In the case of mammals, there is evidence that some degree of specificity can be attained in different situations, for instance as organ-specific protection rather than microorganism-specific reaction. Thus, depending on the challenges and conditions, innate memory can be non-specific or specific, can be integrated in the germline and transmitted to the progeny or be short-lived, thereby representing an exceptionally plastic mechanism of defensive adaptation for ensuring individual and species survival.

Chromosome-level genome assembly of a cliff plant Taihangia rupestris var. ciliata provides insights into its adaptation and demographic history.

BACKGROUND: Cliffs are recognized as one of the most challenging environments for plants, characterized by harsh conditions such as drought, infertile soil, and steep terrain. However, they surprisingly host ancient and diverse plant communities and play a crucial role in protecting biodiversity. The Taihang Mountains, which act as a natural boundary in eastern China, support a rich variety of plant species, including many unique to cliff habitats. However, it is little known how cliff plants adapt to harsh habitats and the demographic history in this region. RESULTS: To better understand the demographic history and adaptation of cliff plants in this area, we analyzed the chromosome-level genome of a representative cliff plant, T. rupestris var. ciliata, which has a genome size of 769.5 Mb, with a scaffold N50 of 104.92 Mb. The rapid expansion of transposable elements may have contributed to the increasing genome and its ability to adapt to unique and challenging cliff habitats. Comparative analysis of the genome evolution between Taihangia and non-cliff plants in Rosaceae revealed a significant expansion of gene families associated with oxidative phosphorylation, which is likely a response to the abiotic stresses faced by cliff plants. This expansion may explain the long-term adaptation of Taihangia to harsh cliff environments. The effective population size of the two varieties has continuously decreased due to climatic fluctuations during the Quaternary period. Furthermore, significant differences in gene expression between the two varieties may explain the varied leaf phenotypes and adaptations to harsh conditions in different natural distributions. CONCLUSION: Our study highlights the extraordinary adaptation of T. rupestris var. ciliata, shedding light on the evolution of cliff plants worldwide.

Phenotyping revealed tolerance traits and genotypes for acidity and aluminum toxicity in European Vicia faba L.

Soil acidity is a global issue; soils with pH <4.5 are widespread in Europe. This acidity adversely affects nutrient availability to plants; pH levels <5.0 lead to aluminum (Al3+) toxicity, a significant problem that hinders root growth and nutrient uptake in faba bean (Vicia faba L.) and its symbiotic relationship with Rhizobium. However, little is known about the specific traits and tolerant genotypes among the European faba beans. This study aimed to identify response traits associated with tolerance to root zone acidity and Al3+ toxicity and potentially tolerant genotypes for future breeding efforts. Germplasm survey was conducted using 165 genotypes in a greenhouse aquaponics system. Data on the root and shoot systems were collected. Subsequently, 12 genotypes were selected for further phenotyping in peat medium, where data on physiological and morphological parameters were recorded along with biochemical responses in four selected genotypes. In the germplasm survey, about 30% of genotypes showed tolerance to acidity and approximately 10% exhibited tolerance to Al3+, while 7% showed tolerance to both. The phenotyping experiment indicated diverse morphological and physiological responses among treatments and genotypes. Acid and Al3+ increased proline concentration. Interaction between genotype and environment was observed for ascorbate peroxidase activity, malondialdehyde, and proline concentrations. Genomic markers associated with acidity and acid+Al3+-toxicity tolerances were identified using GWAS analysis. Four faba bean genotypes with varying levels of tolerance to acidity and Al3+ toxicity were identified.

Biotic interactions promote local adaptation to soil in plants.

Although different ecological factors shape adaptative evolution in natural habitats, we know little about how their interactions impact local adaptation. Here we used eight generations of experimental evolution with outcrossing Brassica rapa plants as a model system, in eight treatment groups that varied in soil type, herbivory (with/without aphids), and pollination mode (hand- or bumblebee-pollination), to study how biotic interactions affect local adaptation to soil. First, we show that several plant traits evolved in response to biotic interactions in a soil-specific way. Second, using a reciprocal transplant experiment, we demonstrate that significant local adaptation to soil-type evolved in the "number of open flowers", a trait used as a fitness proxy, but only in plants that evolved with herbivory and bee pollination. Whole genome re-sequencing of experimental lines revealed that biotic interactions caused a 10-fold increase in the number of SNPs across the genome with significant allele frequency change, and that alleles with opposite allele frequency change in different soil types (antagonistic pleiotropy) were most common in plants with an evolutionary history of herbivory and bee pollination. Our results demonstrate that the interaction with mutualists and antagonists can facilitate local adaptation to soil type through antagonistic pleiotropy.

Uncovering the genetic diversity and adaptability of Butuo Black Sheep through whole-genome re-sequencing.

The Butuo Black Sheep (BBS) is well-known for its ability to thrive at high altitudes, resist diseases, and produce premium-quality meat. Nonetheless, there is insufficient data regarding its genetic diversity and population-specific Single nucleotide polymorphisms (SNPs). This paper centers on the genetic diversity of (BBS). The investigation conducted a whole-genome resequencing of 33 BBS individuals to recognize distinct SNPs exclusive to BBS. The inquiry utilized bioinformatic analysis to identify and explain SNPs and pinpoint crucial mutation sites. The findings reveal that reproductive-related genes (GHR, FSHR, PGR, BMPR1B, FST, ESR1), lipid-related genes (PPARGC1A, STAT6, DGAT1, ACACA, LPL), and protein-related genes (CSN2, LALBA, CSN1S1, CSN1S2) were identified as hub genes. Functional enrichment analysis showed that genes associated with reproduction, immunity, inflammation, hypoxia, PI3K-Akt, and AMPK signaling pathways were present. This research suggests that the unique ability of BBS to adapt to low oxygen levels in the plateau environment may be owing to mutations in a variety of genes. This study provides valuable insights into the genetic makeup of BBS and its potential implications for breeding and conservation efforts. The genes and SPNs identified in this study could serve as molecular markers for BBS.

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

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