Population dynamics of Baltic herring since the Viking Age revealed by ancient DNA and genomics

The Once and Future Fish: Assessing a Millennium of Atlantic Herring Exploitation Through Mixed-Stock Analysis and Ancient DNA

Small pelagic fish support profitable fisheries and are important for food security around the world. Yet, their sustainable management can be hindered by the indiscriminate impacts of simultaneous exploitation of fish from multiple distinct biological populations over extended periods of time. The quantification of such impacts is greatly facilitated by recently developed molecular tools-including diagnostic single nucleotide polymorphism (SNP) panels for mixed-stock analysis (MSA)-that can accurately detect the population identity of individual fish. However, the biological relevance of such tools over longer periods of time remains unknown. Here, we demonstrate that diagnostic SNP panels designed for contemporary MSA in Atlantic herring have a millennium-long biological relevance and applicability. We assign the population identity of ancient Atlantic herring specimens-obtained through famously profitable historic fisheries-up to 1300 years old from eight archaeological sites across Europe. Analyzing contemporary and ancient whole-genome data, we obtain evidence for the long-term mixed-stock exploitation of Atlantic herring. Despite such mixed-stock exploitation, we exclusively identify autumn-spawning herring amongst these archaeological remains, indicative of a specific biological availability or cultural preference for certain herring ecotypes in the past. Moreover, our results show that herring demographic patterns were relatively stable until the dramatic disruptions and stock collapses during the 20th century. We find small but significant reductions in genetic diversity over time, indicating long-term evolutionary consequences from 20th-century stock declines. The long-term applicability of diagnostic SNP panels underscores their biological relevance and cost-effective application for the genetic monitoring of herring stocks and highlights the utility of ancient DNA to obtain insights in herring ecology and population dynamics.

Genome-Wide Population Structure in a Marine Keystone Species, the European Flat Oyster (Ostrea edulis)

Ostrea edulis, the European flat oyster, was once a widespread economically and ecologically important marine species, but has suffered dramatic declines over the past two centuries. Consequently, there has been a surge in European restoration efforts, many of which focus on restocking as a conservation measure. In this study, we used whole-genome sequencing (WGS) data to investigate the population structure, demographic history, and patterns of local adaptation of O. edulis across its natural distribution with increased sampling densities at Scandinavian localities. Results revealed seven distinct genetic clusters, including previously undescribed complex population structure in Norway, and evidence for introgression between genetic clusters in Scandinavia. We detected large structural variants (SVs) on three pseudo-chromosomes. These megabase long regions were characterised by strong linkage disequilibrium and clear geographical differentiation, suggestive of chromosomal inversions potentially associated with local adaptation. The results indicated that genomic traces of past translocations of non-native O. edulis were still present in some individuals, but overall, we found limited evidence of major impacts of translocations on the scale of contemporary population structure. Our findings highlight the importance of considering population structure and signatures of selection in the design of effective conservation strategies to preserve and restore wild native European flat oyster populations, and we provide direct knowledge safeguarding sustainable mitigation actions in this important species.

Rapid speciation in the holopelagic ctenophore Mnemiopsis following glacial recession

Understanding how populations diverge is one of the oldest and most compelling questions in evolutionary biology. An in depth understanding of how this process operates in planktonic marine animals, where barriers for gene flow are seemingly absent, is critical to understanding the past, present, and future of ocean life. Mnemiopsis plays an important ecological role in its native habitat along the Atlantic coast of the Americas and is highly destructive in its non-native habitats in European waters. Although historical literature described three species of Mnemiopsis, the lack of stable morphological characters has led to the collapse of this group into a single species, Mnemiopsis leidyi. We generate high-quality reference genomes and use a whole-genome sequencing approach to reveal that there are two species of Mnemiopsis along its native range and show that historical divergence between the two species coincides with historical glacial melting. We define a hybridization zone between species and highlight that environmental sensing genes likely contribute to the invasive success of Mnemiopsis. Overall, this study provides insights into the fundamental question of how holopelagic species arise without clear barriers to gene flow and sheds light on the genomic mechanisms important for invasion success in a highly invasive species.

Tracing 600 years of long-distance Atlantic cod trade in medieval and post-medieval Oslo using stable isotopes and ancient DNA

Marine resources have been important for the survival and economic development of coastal human communities across northern Europe for millennia. Knowledge of the origin of such historic resources can provide key insights into fishing practices and the spatial extent of trade networks. Here, we combine ancient DNA and stable isotopes (δ13C, δ15N, non-exchangeable δ2H and δ34S) to investigate the geographical origin of archaeological cod remains in Oslo from the eleventh to seventeenth centuries CE. Our findings provide genetic evidence that Atlantic cod was obtained from different geographical populations, including a variety of distant-water populations like northern Norway and possibly Iceland. Evidence for such long-distance cod trade is already observed from the eleventh century, contrasting with archaeological and historical evidence from Britain and other areas of Continental Europe around the North and Baltic Seas, where such trade increased during the thirteenth to fourteenth centuries. The genomic assignments of specimens to different populations coincide with significantly different δ13C values between those same specimens, indicating that multiple Atlantic cod populations living in different environments were exploited. This research provides novel information about the exploitation timeline of specific Atlantic cod stocks and highlights the utility of combining ancient DNA (aDNA) methods and stable isotope analysis to describe the development of medieval and post-medieval marine fisheries.

Genetic Stock Identification Reveals Mismatches Between Management Areas and Population Genetic Structure in a Migratory Pelagic Fish

Sustainable fisheries management is important for the continued harvest of the world's marine resources, especially as they are increasingly challenged by a range of climatic and anthropogenic factors. One of the pillars of sustainable fisheries management is the accurate identification of the biological units, i.e., populations. Here, we developed and implemented a genetic baseline for Atlantic herring harvested in the Norwegian offshore fisheries to investigate the validity of the current management boundaries. This was achieved by genotyping > 15,000 herring from the northern European seas, including samples of all the known populations in the region, with a panel of population-informative SNPs mined from existing genomic resources. The final genetic baseline consisted of ~1000 herring from 12 genetically distinct populations. We thereafter used the baseline to investigate mixed catches from the North and Norwegian Seas, revealing that each management area consisted of multiple populations, as previously suspected. However, substantial numbers (up to 50% or more within a sample) of herring were found outside of their expected management areas, e.g., North Sea autumn-spawning herring north of 62° N (average = 19.2%), Norwegian spring-spawning herring south of 62° N (average = 13.5%), and western Baltic spring-spawning herring outside their assumed distribution area in the North Sea (average = 20.0%). Based upon these extensive observations, we conclude that the assessment and management areas currently in place for herring in this region need adjustments to reflect the populations present. Furthermore, we suggest that for migratory species, such as herring, a paradigm shift from using static geographic stock boundaries towards spatial dynamic boundaries is needed to meet the requirements of future sustainable management regimes.

Archaeogenomic analysis of Chesapeake Atlantic sturgeon illustrates shaping of its populations in recovery from severe overexploitation

Atlantic sturgeon (Acipenser oxyrinchus ssp. oxyrinchus) has been a food resource in North America for millennia. However, industrial-scale fishing activities following the establishment of European colonies led to multiple collapses of sturgeon stocks, driving populations such as those in the Chesapeake area close to extinction. While recent conservation efforts have been successful in restoring census numbers, little is known regarding genomic consequences of the population bottleneck. Here, we characterize its effect on present-day population structuring and genomic diversity in James River populations. To establish a pre-collapse baseline, we collected genomic data from archaeological remains from Middle Woodland Maycock's Point (c. 200-900 CE), as well as Jamestown and Williamsburg colonial sites. Demographic analysis of recovered mitogenomes reveals a historical collapse in effective population size, also reflected in diminished present-day mitogenomic diversity and structure. We infer that James River fall- and spring-spawning populations likely took shape in recent years of population recovery, where genetic drift enhanced the degree of population structure. The mismatch of mitogenomic lineages to geographical-seasonal groupings implies that despite their homing instinct and differential adaptation manifested as season-specific behaviour, colonization of new rivers has been a key ecological strategy for Atlantic sturgeon over evolutionary timescales.

Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species

Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.

Impact of population structure in the estimation of recent historical effective population size by the software GONE

BACKGROUND

Effective population size (Ne) is a crucial parameter in conservation genetics and animal breeding. A recent method, implemented by the software GONE, has been shown to be rather accurate in estimating recent historical changes in Ne from a single sample of individuals. However, GONE estimations assume that the population being studied has remained isolated for a period of time, that is, without migration or confluence of other populations. If this occurs, the estimates of Ne can be heavily biased. In this paper, we evaluate the impact of migration and admixture on the estimates of historical Ne provided by GONE through a series of computer simulations considering several scenarios: (a) the mixture of two or more ancestral populations; (b) subpopulations that continuously exchange individuals through migration; (c) populations receiving migrants from a large source; and (d) populations with balanced systems of chromosomal inversions, which also generate genetic structure.

RESULTS

Our results indicate that the estimates of historical Ne provided by GONE may be substantially biased when there has been a recent mixture of populations that were previously separated for a long period of time. Similarly, biases may occur when the rate of continued migration between populations is low, or when chromosomal inversions are present at high frequencies. However, some biases due to population structuring can be eliminated by conducting population structure analyses and restricting the estimation to the differentiated groups. In addition, disregarding the genomic regions that are involved in inversions can also remove biases in the estimates of Ne.

CONCLUSIONS

Different kinds of deviations from isolation and panmixia of the populations can generate biases in the recent historical estimates of Ne. Therefore, estimation of past demography could benefit from performing population structure analyses beforehand, by mitigating the impact of these biases on historical Ne estimates.

Ancient DNA is preserved in fish fossils from tropical lake sediments

Tropical freshwater lakes are well known for their high biodiversity, and particularly the East African Great Lakes are renowned for their adaptive radiation of cichlid fishes. While comparative phylogenetic analyses of extant species flocks have revealed patterns and processes of their diversification, little is known about evolutionary trajectories within lineages, the impacts of environmental drivers, or the scope and nature of now-extinct diversity. Time-structured palaeodata from geologically young fossil records, such as fossil counts and particularly ancient DNA (aDNA) data, would help fill this large knowledge gap. High ambient temperatures can be detrimental to the preservation of DNA, but refined methodology now allows data generation even from very poorly preserved samples. Here, we show for the first time that fish fossils from tropical lake sediments yield endogenous aDNA. Despite generally low endogenous content and high sample dropout, the application of high-throughput sequencing and, in some cases, sequence capture allowed taxonomic assignment and phylogenetic placement of 17% of analysed fish fossils to family or tribe level, including remains which are up to 2700 years old or weigh less than 1 mg. The relationship between aDNA degradation and the thermal age of samples is similar to that described for terrestrial samples from cold environments when adjusted for elevated temperature. Success rates and aDNA preservation differed between the investigated lakes Chala, Kivu and Victoria, possibly caused by differences in bottom water oxygenation. Our study demonstrates that the sediment records of tropical lakes can preserve genetic information on rapidly diversifying fish taxa over time scales of millennia.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.

Understanding and applying biological resilience, from genes to ecosystems

The natural world is under unprecedented and accelerating pressure. Much work on understanding resilience to local and global environmental change has, so far, focussed on ecosystems. However, understanding a system's behaviour requires knowledge of its component parts and their interactions. Here we call for increased efforts to understand 'biological resilience', or the processes that enable components across biological levels, from genes to communities, to resist or recover from perturbations. Although ecologists and evolutionary biologists have the tool-boxes to examine form and function, efforts to integrate this knowledge across biological levels and take advantage of big data (e.g. ecological and genomic) are only just beginning. We argue that combining eco-evolutionary knowledge with ecosystem-level concepts of resilience will provide the mechanistic basis necessary to improve management of human, natural and agricultural ecosystems, and outline some of the challenges in achieving an understanding of biological resilience.