Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture.

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

Evolution at two time frames: ancient structural variants involved in post-glacial divergence of the European plaice (Pleuronectes platessa).

Changing environmental conditions can lead to population diversification through differential selection on standing genetic variation. Structural variant (SV) polymorphisms provide examples of ancient alleles that in time become associated with novel environmental gradients. The European plaice (Pleuronectes platessa) is a marine flatfish showing large allele-frequency differences at two putative SVs associated with environmental variation. In this study, we explored the contribution of these SVs to population structure across the North East Atlantic. We compared genome-wide population structure using sets of RAD-sequencing SNPs with the spatial structure of the SVs. We found that in contrast to the rest of the genome, the SVs were only weakly associated with an isolation-by-distance pattern. Indeed, both SVs showed important variation in haplogroup frequencies, with the same haplogroup increasing both along the salinity gradient of the Baltic Sea, and found in high frequency in the northern-range margin of the Atlantic. Phylogenetic analyses suggested that the SV alleles are much older than the age of the Baltic Sea itself. These results suggest that the SVs are older than the age of the environmental gradients with which they currently co-vary. Altogether, our results suggest that the plaice SVs were shaped by evolutionary processes occurring at two time frames, firstly following their origin, ancient spread and maintenance in the ancestral populations, and secondly related to their current association with more recently formed environmental gradients such as those found in the North Sea-Baltic Sea transition zone.

Genomic Signatures After Five Generations of Intensive Selective Breeding: Runs of Homozygosity and Genetic Diversity in Representative Domestic and Wild Populations of Turbot (Scophthalmus maximus).

Massive genotyping of single nucleotide polymorphisms (SNP) has opened opportunities for analyzing the way in which selection shapes genomes. Artificial or natural selection usually leaves genomic signatures associated with selective sweeps around the responsible locus. Strong selective sweeps are most often identified either by lower genetic diversity than the genomic average and/or islands of runs of homozygosity (ROHi). Here, we conducted an analysis of selective sweeps in turbot (Scophthalmus maximus) using two SNP datasets from a Northeastern Atlantic population (36 individuals) and a domestic broodstock (46 individuals). Twenty-six families (∼ 40 offspring per family) from this broodstock and three SNP datasets applying differing filtering criteria were used to adjust ROH calling parameters. The best-fitted genomic inbreeding estimate (FROH) was obtained by the sum of ROH longer than 1 Mb, called using a 21,615 SNP panel, a sliding window of 37 SNPs and one heterozygous SNP per window allowed. These parameters were used to obtain the ROHi distribution in the domestic and wild populations (49 and 0 ROHi, respectively). Regions with higher and lower genetic diversity within each population were obtained using sliding windows of 37 SNPs. Furthermore, those regions were mapped in the turbot genome against previously reported genetic markers associated with QTL (Quantitative Trait Loci) and outlier loci for domestic or natural selection to identify putative selective sweeps. Out of the 319 and 278 windows surpassing the suggestive pooled heterozygosity thresholds (ZHp) in the wild and domestic population, respectively, 78 and 54 were retained under more restrictive ZHp criteria. A total of 116 suggestive windows (representing 19 genomic regions) were linked to either QTL for production traits, or outliers for divergent or balancing selection. Twenty-four of them (representing 3 genomic regions) were retained under stricter ZHp thresholds. Eleven QTL/outlier markers were exclusively found in suggestive regions of the domestic broodstock, 7 in the wild population and one in both populations; one (broodstock) and two (wild) of those were found in significant regions retained under more restrictive ZHp criteria in the broodstock and the wild population, respectively. Genome mining and functional enrichment within regions associated with selective sweeps disclosed relevant genes and pathways related to aquaculture target traits, including growth and immune-related pathways, metabolism and response to hypoxia, which showcases how this genome atlas of genetic diversity can be a valuable resource to look for candidate genes related to natural or artificial selection in turbot populations.

Weak genetic structure despite strong genomic signal in lesser sandeel in the North Sea.

Sandeels are an ecologically important group of fishes; they are a key part of the food chain serving as food for marine mammals, seabirds and fish. Sandeels are further targeted by a large industrial fishery, which has led to concern about ecosystem effects. In the North Sea, the lesser sandeel Ammodytes marinus is by far the most prevalent species of sandeel in the fishery. Management of sandeel in the North Sea plus the Kattegat is currently divided into seven geographical areas, based on subtle differences in demography, population dynamics and results from simulations of larval dispersal. However, little is known about the underlying genetic population structure. In this study, we used 2,522 SNPs derived from restriction site-associated DNA sequencing (RADseq) typed in 429 fish representing four main sandeel management areas. Our main results showed (a) a lack of a clear spatially defined genetic structure across the majority of genetic markers and (b) the existence of a group of at least 13 SNPs under strong linkage disequilibrium which together separate North Sea sandeel into three haplotype clusters, suggestive of one or more structural variants in the genome. Analyses of the spatial distribution of these putative structural variants suggest at least partial reproductive isolation of sandeel in the western management area along the Scottish coast, supporting a separate management. Our results highlight the importance of the application of a large number of markers to be able to detect weak patterns of differentiation. This study contributes to increasing the genetic knowledge of this important exploited species, and results can be used to improve our understanding of population dynamics and stock structure.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genetic analyses reveal complex dynamics within a marine fish management area.

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Stable coexistence of genetically divergent Atlantic cod ecotypes at multiple spatial scales.

Coexistence in the same habitat of closely related yet genetically different populations is a phenomenon that challenges our understanding of local population structure and adaptation. Identifying the underlying mechanisms for such coexistence can yield new insight into adaptive evolution, diversification and the potential for organisms to adapt and persist in response to a changing environment. Recent studies have documented cryptic, sympatric populations of Atlantic cod (Gadus morhua) in coastal areas. We analysed genetic origin of 6,483 individual cod sampled annually over 14 years from 125 locations along the Norwegian Skagerrak coast and document stable coexistence of two genetically divergent Atlantic cod ecotypes throughout the study area and study period. A "fjord" ecotype dominated in numbers deep inside fjords while a "North Sea" ecotype was the only type found in offshore North Sea. Both ecotypes coexisted in similar proportions throughout coastal habitats at all spatial scales. The size-at-age of the North Sea ecotype on average exceeded that of the fjord ecotype by 20% in length and 80% in weight across all habitats. Different growth and size among individuals of the two types might be one of several ecologically significant variables that allow for stable coexistence of closely related populations within the same habitat. Management plans, biodiversity initiatives and other mitigation strategies that do not account for the mixture of species ecotypes are unlikely to meet objectives related to the sustainability of fish and fisheries.

Range-wide genomic data synthesis reveals transatlantic vicariance and secondary contact in Atlantic cod.

Recent advances in genetic and genomic analysis have greatly improved our understanding of spatial population structure in marine species. However, studies addressing phylogeographic patterns at oceanic spatial scales remain rare. In Atlantic cod (Gadus morhua), existing range-wide examinations suggest significant transatlantic divergence, although the fine-scale contemporary distribution of populations and potential for secondary contact are largely unresolved. Here, we explore transatlantic phylogeography in Atlantic cod using a data-synthesis approach, integrating multiple genome-wide single-nucleotide polymorphism (SNP) datasets representative of different regions to create a single range-wide dataset containing 1,494 individuals from 54 locations and genotyped at 796 common loci. Our analysis highlights significant transatlantic divergence and supports the hypothesis of westward post-glacial colonization of Greenland from the East Atlantic. Accordingly, our analysis suggests the presence of transatlantic secondary contact off eastern North America and supports existing perspectives on the phylogeographic history of Atlantic cod with an unprecedented combination of genetic and geographic resolution. Moreover, we demonstrate the utility of integrating distinct SNP databases of high comparability.

Genome architecture enables local adaptation of Atlantic cod despite high connectivity.

Adaptation to local conditions is a fundamental process in evolution; however, mechanisms maintaining local adaptation despite high gene flow are still poorly understood. Marine ecosystems provide a wide array of diverse habitats that frequently promote ecological adaptation even in species characterized by strong levels of gene flow. As one example, populations of the marine fish Atlantic cod (Gadus morhua) are highly connected due to immense dispersal capabilities but nevertheless show local adaptation in several key traits. By combining population genomic analyses based on 12K single nucleotide polymorphisms with larval dispersal patterns inferred using a biophysical ocean model, we show that Atlantic cod individuals residing in sheltered estuarine habitats of Scandinavian fjords mainly belong to offshore oceanic populations with considerable connectivity between these diverse ecosystems. Nevertheless, we also find evidence for discrete fjord populations that are genetically differentiated from offshore populations, indicative of local adaptation, the degree of which appears to be influenced by connectivity. Analyses of the genomic architecture reveal a significant overrepresentation of a large ~5 Mb chromosomal rearrangement in fjord cod, previously proposed to comprise genes critical for the survival at low salinities. This suggests that despite considerable connectivity with offshore populations, local adaptation to fjord environments may be enabled by suppression of recombination in the rearranged region. Our study provides new insights into the potential of local adaptation in high gene flow species within fine geographical scales and highlights the importance of genome architecture in analyses of ecological adaptation.

Genetic diversity and connectivity within Mytilus spp. in the subarctic and Arctic.

Climate changes in the Arctic are predicted to alter distributions of marine species. However, such changes are difficult to quantify because information on present species distribution and the genetic variation within species is lacking or poorly examined. Blue mussels, Mytilus spp., are ecosystem engineers in the coastal zone globally. To improve knowledge of distribution and genetic structure of the Mytilus edulis complex in the Arctic, we analyzed 81 SNPs in 534 Mytilus spp. individuals sampled at 13 sites to provide baseline data for distribution and genetic variation of Mytilus mussels in the European Arctic. Mytilus edulis was the most abundant species found with a clear genetic split between populations in Greenland and the Eastern Atlantic. Surprisingly, analyses revealed the presence of Mytilus trossulus in high Arctic NW Greenland (77°N) and Mytilus galloprovincialis or their hybrids in SW Greenland, Svalbard, and the Pechora Sea. Furthermore, a high degree of hybridization and introgression between species was observed. Our study highlights the importance of distinguishing between congener species, which can display local adaptation and suggests that information on dispersal routes and barriers is essential for accurate predictions of regional susceptibility to range expansions or invasions of boreal species in the Arctic.

Historical DNA documents long-distance natal homing in marine fish.

The occurrence of natal homing in marine fish remains a fundamental question in fish ecology as its unequivocal demonstration requires tracking of individuals from fertilization to reproduction. Here, we provide evidence of long-distance natal homing (>1000 km) over more than 60 years in Atlantic cod (Gadus morhua), through genetic analysis of archived samples from marked and recaptured individuals. Using a high differentiation single-nucleotide polymorphism assay, we demonstrate that the vast majority of cod tagged in West Greenland and recaptured on Icelandic spawning grounds belonged to the Iceland offshore population, strongly supporting a hypothesis of homing. The high degree of natal fidelity observed provides the evolutionary settings for development of locally adapted populations in marine fish and emphasize the need to consider portfolio effects in marine fisheries management strategies.

Conserved structure and expression of hsp70 paralogs in teleost fishes.

The cytosolic 70KDa heat shock proteins (Hsp70s) are widely used as biomarkers of environmental stress in ecological and toxicological studies in fish. Here we analyze teleost genome sequences to show that two genes encoding inducible hsp70s (hsp70-1 and hsp70-2) are likely present in all teleost fish. Phylogenetic and synteny analyses indicate that hsp70-1 and hsp70-2 are distinct paralogs that originated prior to the diversification of the teleosts. The promoters of both genes contain a TATA box and conserved heat shock elements (HSEs), but unlike mammalian HSP70s, both genes contain an intron in the 5' UTR. The hsp70-2 gene has undergone tandem duplication in several species. In addition, many other teleost genome assemblies have multiple copies of hsp70-2 present on separate, small, genomic scaffolds. To verify that these represent poorly assembled tandem duplicates, we cloned the genomic region surrounding hsp70-2 in Fundulus heteroclitus and showed that the hsp70-2 gene copies that are on separate scaffolds in the genome assembly are arranged as tandem duplicates. Real-time quantitative PCR of F. heteroclitus genomic DNA indicates that four copies of the hsp70-2 gene are likely present in the F. heteroclitus genome. Comparison of expression patterns in F. heteroclitus and Gasterosteus aculeatus demonstrates that hsp70-2 has a higher fold increase than hsp70-1 following heat shock in gill but not in muscle tissue, revealing a conserved difference in expression patterns between isoforms and tissues. These data indicate that ecological and toxicological studies using hsp70 as a biomarker in teleosts should take this complexity into account.

Archived DNA reveals fisheries and climate induced collapse of a major fishery.

Fishing and climate change impact the demography of marine fishes, but it is generally ignored that many species are made up of genetically distinct locally adapted populations that may show idiosyncratic responses to environmental and anthropogenic pressures. Here, we track 80 years of Atlantic cod (Gadus morhua) population dynamics in West Greenland using DNA from archived otoliths in combination with fish population and niche based modeling. We document how the interacting effects of climate change and high fishing pressure lead to dramatic spatiotemporal changes in the proportions and abundance of different genetic populations, and eventually drove the cod fishery to a collapse in the early 1970s. Our results highlight the relevance of fisheries management at the level of genetic populations under future scenarios of climate change.

Population genomics of marine fishes: next-generation prospects and challenges.

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

The use of archived tags in retrospective genetic analysis of fish.

Collections of historical tissue samples from fish (e.g. scales and otoliths) stored in museums and fisheries institutions are precious sources of DNA for conducting retrospective genetic analysis. However, in some cases, only external tags used for documentation of spatial dynamics of fish populations have been preserved. Here, we test the usefulness of fish tags as a source of DNA for genetic analysis. We extract DNA from historical tags from cod collected in Greenlandic waters between 1950 and 1968. We show that the quantity and quality of DNA recovered from tags is comparable to DNA from archived otoliths from the same individuals. Surprisingly, levels of cross-contamination do not seem to be significantly higher in DNA from external (tag) than internal (otolith) sources. Our study therefore demonstrates that historical tags can be a highly valuable source of DNA for retrospective genetic analysis of fish.

Spatiotemporal SNP analysis reveals pronounced biocomplexity at the northern range margin of Atlantic cod Gadus morhua.

Accurate prediction of species distribution shifts in the face of climate change requires a sound understanding of population diversity and local adaptations. Previous modeling has suggested that global warming will lead to increased abundance of Atlantic cod (Gadus morhua) in the ocean around Greenland, but the dynamics of earlier abundance fluctuations are not well understood. We applied a retrospective spatiotemporal population genomics approach to examine the temporal stability of cod population structure in this region and to search for signatures of divergent selection over a 78-year period spanning major demographic changes. Analyzing >900 gene-associated single nucleotide polymorphisms in 847 individuals, we identified four genetically distinct groups that exhibited varying spatial distributions with considerable overlap and mixture. The genetic composition had remained stable over decades at some spawning grounds, whereas complete population replacement was evident at others. Observations of elevated differentiation in certain genomic regions are consistent with adaptive divergence between the groups, indicating that they may respond differently to environmental variation. Significantly increased temporal changes at a subset of loci also suggest that adaptation may be ongoing. These findings illustrate the power of spatiotemporal population genomics for revealing biocomplexity in both space and time and for informing future fisheries management and conservation efforts.

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod.

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80-year period. Screening of >1000 gene-associated single-nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.

A genomic island linked to ecotype divergence in Atlantic cod.

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

Illegal, Unreported and Unregulated fishing has had a major role in the overexploitation of global fish populations. In response, international regulations have been imposed and many fisheries have been 'eco-certified' by consumer organizations, but methods for independent control of catch certificates and eco-labels are urgently needed. Here we show that, by using gene-associated single nucleotide polymorphisms, individual marine fish can be assigned back to population of origin with unprecedented high levels of precision. By applying high differentiation single nucleotide polymorphism assays, in four commercial marine fish, on a pan-European scale, we find 93-100% of individuals could be correctly assigned to origin in policy-driven case studies. We show how case-targeted single nucleotide polymorphism assays can be created and forensically validated, using a centrally maintained and publicly available database. Our results demonstrate how application of gene-associated markers will likely revolutionize origin assignment and become highly valuable tools for fighting illegal fishing and mislabelling worldwide.