Mariann Árnyasi

The Duroc-derived reference genome Sscrofa11.1 has provided a critical foundation for pig genomics, providing a high-quality reference genome for accurate variant detection and comparative genomics but does not capture breed-specific variation. Here, we present a near-complete, gap-free genome assembly for the Landrace pig (Landrace_v1, GCA_963921485.1), spanning all 20 chromosomes and totaling 2.6 Gb, including 176 Mb of sequence absent from Sscrofa11.1. Comparative analyses with recently published high-quality pig genomes reveal a conserved centromere organization across breeds, accompanied by substantial variation in repeat composition and length, and identify a pig specific pattern of telomere variant repeats across eight pig breeds. The improved resolution of repetitive regions in Landrace_v1 enables more complete reconstruction of complex gene families, including olfactory receptors, and uncovers structural variation at the KIT proto-oncogene receptor tyrosine kinase locus not represented in the Duroc reference. Together, these findings highlight the limitations of single-reference genomes and demonstrate the value of breed-specific assemblies for capturing genomic diversity and improving downstream analyses.

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.

The phylum Bryozoa is an understudied, yet commonly occurring, globally distributed bilaterian metazoan organismal group. They have a colonial lifestyle and an evolutionary history that spans at least 480 million years but likely longer. Despite their contentious phylogenetic affinities among metazoans, disproportionately few genomic investigations have been performed thus far. Here, we describe the first chromosome-level genome assembly of an individual Flustra foliacea colony belonging to the order Cheilostomatida, collected in southern Norway. The haplotype-resolved assembly of F. foliacea contains two pseudo-haplotypes spanning 956 megabases and 880 megabases, respectively. Both assemblies are highly complete both in terms of scaffolding (>90% of sequences placed in eight autosomal chromosomal pseudomolecules), and gene content (BUSCO completeness scores >90%). We also present gene and repeat annotations of the two assemblies. A comparison of our newly sequenced F. foliacea with five previously published bryozoan genomes supports the hypothesis that the group has undergone extensive genome rearrangements. This includes multiple chromosomal fusions in F. foliacea since their split with other cheilostome bryozoans. These fusions were enriched with long terminal repeat retrotransposons, highlighting the complex interplay between genome organization and genomic repeats. Our study contributes to a deeper understanding of bryozoan genome evolution and the role of repeats in metazoan genome organization.

Norwegian Red (NR) cattle are the main dairy breed in Norway, bred according to a broad breeding goal including health and fertility since the 1970s. Genomic studies on NR cattle have relied on the public Hereford reference, thus increasing the risk of missing or misrepresenting NR breed-specific variation. Moreover, the Hereford reference is a pseudohaploid assembly, representing homologous chromosomes in a collapsed manner, which results in loss of haplotype-specific alleles and misrepresentation of complex variants. To develop more refined NR specific resources, we utilised long-read sequencing (PacBio HiFi + ONT) and trio-binning to construct six new haplotype-resolved assemblies representing NR genomes. These six NR2025 assemblies show high completeness (BUSCO: 95.82-98.11%) and contiguity (N50: 73.8-88.5 Mb) and are accurately phased (hamming error rate: 0.46-2.52%). Most autosomes have been assembled into the acrocentric centromere, and mapping of bovine satellite sequences reveal distinct organisational patterns of different satellite units across this highly repetitive region. The NR2025 assemblies provide a valuable resource for identification of novel variants and haplotypes in the NR population, which will enable more accurate association studies of genotype-to-phenotype relations and genomic predictions in NR cattle, ultimately enhancing the efficacy of selection for desirable traits.

The phylum Bryozoa is an understudied, yet commonly-occurring, globally distributed bilaterian metazoan organismal group. They have a colonial lifestyle and an evolutionary history that spans at least 480 million years but likely longer. Despite their contentious phylogenetic affinities among metazoans, disproportionately little genomic investigations have been performed thus far. Here, we describe the first chromosome-level genome assembly of an individual Flustra foliacea colony belonging to the order Cheilostomatida, collected in southern Norway. The haplotype-resolved assembly of F. foliacea contains two pseudo-haplotypes spanning 956 megabases and 880 megabases, respectively. Both assemblies are highly complete both in terms of scaffolding (>90% of sequences placed in 8 autosomal chromosomal pseudomolecules), and gene content (BUSCO completeness scores >90%). We also present gene and repeat annotations of the two assemblies. A comparison of our newly sequenced F. foliacea with five previously published bryozoan genomes supports the hypothesis that the group has undergone extensive genome rearrangements. This includes multiple chromosomal fusions in F. foliacea s ince their split with other cheilostome bryozoans. These fusions were enriched with long terminal repeat (LTR) retrotransposons, highlighting the complex interplay between genome organization and genomic repeats. Our study contributes to a deeper understanding of bryozoan genome evolution and the role of repeats in metazoan genome organization.

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.

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon ( Salmo salar ), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (< 50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and at smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F _ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

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.

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.

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.