Large-scale climate changes influence the geographic distribution of biodiversity. Many taxa have been reported to extend or reduce their geographic range, move poleward or displace other species. However, for closely related species that can hybridize in the natural environment, displacement is not the only effect of changes of environmental variables. Another option is subtler, hidden expansion, which can be found using genetic methods only. The marine blue mussels Mytilus are known to change their geographic distribution despite being sessile animals. In addition to natural dissemination at larval phase—enhanced by intentional or accidental introductions and rafting—they can spread through hybridization and introgression with local congeners, which can create mixed populations sustaining in environmental conditions that are marginal for pure taxa. The Mytilus species have a wide distribution in coastal regions of the Northern and Southern Hemisphere. In this study, we investigated the inter-regional genetic differentiation of the Mytilus species complex at 53 locations in the North Atlantic and adjacent Arctic waters and linked this genetic variability to key local environmental drivers. Of seventy-nine candidate single nucleotide polymorphisms (SNPs), all samples were successfully genotyped with a subset of 54 SNPs. There was a clear interregional separation of Mytilus species. However, all three Mytilus species hybridized in the contact area and created hybrid zones with mixed populations. Boosted regression trees (BRT) models showed that inter-regional variability was important in many allele models but did not prevail over variability in local environmental factors. Local environmental variables described over 40% of variability in about 30% of the allele frequencies of Mytilus spp. For the 30% of alleles, variability in their frequencies was only weakly coupled with local environmental conditions. For most studied alleles the linkages between environmental drivers and the genetic variability of Mytilus spp. were random in respect to “coding” and “non-coding” regions. An analysis of the subset of data involving functional genes only showed that two SNPs at Hsp70 and ATPase genes correlated with environmental variables. Total predictive ability of the highest performing models (r² between 0.550 and 0.801) were for alleles that discriminated most effectively M. trossulus from M. edulis and M. galloprovincialis, whereas the best performing allele model (BM101A) did the best at discriminating M. galloprovincialis from M. edulis and M. trossulus. Among the local environmental variables, salinity, water temperature, ice cover and chlorophyll a concentration were by far the greatest predictors, but their predictive performance varied among different allele models. In most cases changes in the allele frequencies along these environmental gradients were abrupt and occurred at a very narrow range of environmental variables. In general, regions of change in allele frequencies for M. trossulus occurred at 8–11 psu, 0–10 °C, 60%–70% of ice cover and 0–2 mg m⁻³ of chlorophyll a, M. edulis at 8–11 and 30–35 psu, 10–14 °C and 60%–70% of ice cover and for M. galloprovincialis at 30–35 psu, 14–20 °C.

Blue mussels of the genus Mytilus are important ecosystem engineers in intertidal and subtidal communities. The distribution of Mytilus mussels is influenced by a number of benthic and pelagic environmental variables (e.g., substratum type and availability, water movement, phytoplankton production, physical disturbance) as well as interactions between these variables. Because of its broad tolerance of environmental variation the Mediterranean species, Mytilus galloprovincialis, has the greatest ability of all blue mussels to colonise new geographic regions. Understanding how population genetic variation is related to, or caused by, environmental variation is important but has long been a challenge. The present study examined the genetic differentiation of native populations of M. galloprovincialis throughout its entire geographic range in the Mediterranean Sea, the Black Sea and the Sea of Azov using 53 single nucleotide polymorphisms (SNP loci). Mussels, in total 1004 individuals collected from 36 locations, were genotyped and combined with existing SNP data for mussels from 11 reference sites. Pairwise comparisons of F_ST values, correspondence analysis (CA) and STRUCTURE analysis all revealed four groups of populations: the Atlantic Ocean; the western Mediterranean; the Aegean Sea; and the Azov, Black and Marmara Seas. One population – from Algeria (Oran West) – was intermediate between the two main groups of the Mediterranean Sea and Atlantic Ocean. Seascape genetic analyses using GLM and DistLM analyses were employed to test site-specific genetic variation as a function of 13 environmental variables. The GLM identified five environmental variables that explained variation in site-specific F_ST values, whilst in the DistLM best-fit model only four were significant. These analyses suggest that a complex mix of environmental variables contribute to explaining genetic variation of M. galloprovincialis populations within the Mediterranean Sea, which most likely reflects the complex geological history of formation, isolation and reconnection among the regional sub-basins of the Sea.

The impacts and interactions between hatchery-bred fish and wild fish populations has been a topic of active investigation in recent decades. In some instances, the benefits of stocking can be overshadowed by negative effects such as genetic introgression with natural populations, loss of genetic diversity, and dilution of local adaptations. Methods that facilitate the identification of stocked fish enable us to estimate not only the effectiveness of stocking but also the level of natural reproduction and the degree of hybridization. The longest Baltic river, the Vistula, also has the second highest discharge. Historically, it hosted numerous populations of the anadromous form of brown trout (sea trout); however, dam construction has since interfered with and reduced spawning migration to a rate that is much lower than before. Reduced spawning has resulted in a population collapse and a negative flow-on effect on commercial catches. In response, Poland (along with many other Baltic countries) initiated an intensive stocking program which continues today and which sees the average annual release of 700,000 smolts. As a consequence, today’s main-river and inshore catches come from stock-enhanced populations. High-throughput single-nucleotide polymorphism (SNP) genotyping was performed on samples of sea trout from southern Baltic populations; results suggest that a significant portion of the sea trout catches in the Vistula mouth region have direct hatchery origin and indicate the presence of Pomeranian specimens. SNP loci identified as outliers indicate a potential selection pressure that may be related with effects of hatchery breeding and mixing with natural populations. The brown trout SNP array applied in this study showed high effectiveness not only for population differentiation, but more importantly, it emerged as a sensitive tool to provide evidence of detection selection.

Comparative genomic studies in birds have revealed that bird genomes are relatively repeat-poor and stable in terms of karyotype, size, and gene synteny/collinearity compared to other vertebrates. One notable exception is the owls, with cytogenetic studies demonstrating large variations in karyotypes and the evolution of unusual centromeric satellite repeats in some species. However, there has so far not been an investigation into genome architecture and repeat landscape of owls. Here, we present a chromosome-level genome assembly for the snowy owl ( Bubo scandiacus ). We find that the repeat DNA content in the relatively large snowy owl genome (1.6 Gb) is among the highest reported for any bird genome to date (28.34% compared to an average of ∼10% in other birds). The bulk of the snowy owl genomic repeat landscape consists of centromeric satellite DNA, which appears to have originated from an endogenous retrovirus (ERV1). Using gene collinearity analyses we show that the position of these evolutionary new centromeres (ECNs) are not homologous with chicken centromeres, and are located in regions with collinearity breaks to other bird genomes due to chromosomal rearrangements. Our results support rapid transposable element-driven evolution of lineage-specific centromeres, which could have played a role in reproductive isolation and speciation of the snowy owl.

Comparative genomic studies in birds have revealed that bird genomes are relatively repeat-poor and stable in terms of karyotype, size, and gene synteny/collinearity compared to other vertebrates. One notable exception is the owls, with cytogenetic studies demonstrating large variations in karyotypes and the evolution of unusual centromeric satellite repeats in some species. However, there has so far not been an investigation into genome architecture and repeat landscape of owls. Here, we present a chromosome-level genome assembly for the snowy owl (Bubo scandiacus). We find that the repeat DNA content in the relatively large snowy owl genome (1.6 Gb) is among the highest reported for any bird genome to date (28.34% compared to an average of ~ 10% in other birds). The bulk of the snowy owl genomic repeat landscape consists of extensive centromeric satellite arrays that form large contiguous blocks across multiple chromosomes. These arrays represent one of the largest centromeric satellite expansions reported in birds and are consistent with rapid turnover and amplification of centromeric repeats. Using gene collinearity analyses we show that the position of the centromeres are not homologous with chicken centromeres, and are located in regions with collinearity breaks to other bird genomes due to chromosomal rearrangements. Our results support rapid satellite turnover at lineage-specific centromeres, which could have played a role in centromeric drive, reproductive isolation and speciation of the snowy owl. The extensive centromeric satellite arrays uncovered in this genome highlight centromere evolution as a potentially important driver of genome evolution in owls.

Currently available genome assemblies for Atlantic cod (Gadus morhua) have been constructed using DNA from fish belonging to the Northeast Arctic Cod (NEAC) population; a migratory population feeding in the cold Barents Sea. These assemblies have been crucial for the development of genetic markers which have been used to study population differentiation and adaptive evolution in Atlantic cod, pinpointing four discrete islands of genomic divergence located on linkage groups 1, 2, 7 and 12. In this paper, we present a high-quality reference genome from a male Atlantic cod representing a southern population inhabiting the Celtic sea. Structurally, the genome assembly (gadMor_Celtic) was produced from long-read nanopore data and has a combined contig size of 686 Mb with a N50 of 10 Mb. Integrating contigs with genetic linkage mapping information enabled us to construct 23 chromosome sequences which mapped with high confidence to the latest NEAC population assembly (gadMor3) and allowed us to characterize in detail large chromosomal inversions on linkage groups 1, 2, 7 and 12. In most cases, inversion breakpoints could be located within single nanopore contigs. Our results suggest the presence of inversions in Celtic cod on linkage groups 6, 11 and 21, although these remain to be confirmed. Further, we identified a specific repetitive element that is relatively enriched at predicted centromeric regions. Our gadMor_Celtic assembly provides a resource representing a ‘southern’ cod population which is complementary to the existing ‘northern’ population based genome assemblies and represents the first step towards developing pan-genomic resources for Atlantic cod.