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.

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.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s)

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Mitochondrial genomes reveal mid-Pleistocene population divergence, and post-glacial expansion, in Australasian snapper (Chrysophrys auratus)

Glacial cycles play important roles in determining the phylogeographic structure of terrestrial species, however, relatively little is known about their impacts on the distribution of marine biota. This study utilised modern (n = 350) and ancient (n = 26) mitochondrial genomes from Australasian snapper (Chrysophrys auratus) sampled in New Zealand to assess their demographic and phylogeographic history. We also tested for changes in genetic diversity using the up to 750-year-old mitochondrial genomes from pre-European archaeological sites to assess the potential impacts of human exploitation. Nucleotide diversity and haplotype diversity was high (π = 0.005, h = 0.972). There was no significant change in nucleotide diversity over the last 750 years (p = 0.343), with no detectable loss of diversity as a result of indigenous and industrial-scale fishing activity. While there was no evidence for contemporary population structure (AMOVA, p = 0.764), phylogeographic analyses identified two distinct mitochondrial clades that diverged approximately 650,000 years ago during the mid-Pleistocene, suggesting the species experienced barriers to gene flow when sea levels dropped over 120 m during previous glacial maxima. An exponential population increase was also observed around 8000 years ago consistent with a post-glacial expansion, which was likely facilitated by increased ocean temperatures and rising sea levels. This study demonstrates that glacial cycles likely played an important role in the demographic history of C. auratus and adds to our growing understanding of how dynamic climatic changes have influenced the evolution of coastal marine species.

Ancient DNA evidence for the ecological globalization of cod fishing in medieval and post-medieval Europe

Understanding the historical emergence and growth of long-range fisheries can provide fundamental insights into the timing of ecological impacts and the development of coastal communities during the last millennium. Whole-genome sequencing approaches can improve such understanding by determining the origin of archaeological fish specimens that may have been obtained from historic trade or distant water. Here, we used genome-wide data to individually infer the biological source of 37 ancient Atlantic cod specimens (ca 1050-1950 CE) from England and Spain. Our findings provide novel genetic evidence that eleventh- to twelfth-century specimens from London were predominantly obtained from nearby populations, while thirteenth- to fourteenth-century specimens were derived from distant sources. Our results further suggest that Icelandic cod was indeed exported to London earlier than previously reported. Our observations confirm the chronology and geography of the trans-Atlantic cod trade from Newfoundland to Spain starting by the early sixteenth century. Our findings demonstrate the utility of whole-genome sequencing and ancient DNA approaches to describe the globalization of marine fisheries and increase our understanding regarding the extent of the North Atlantic fish trade and long-range fisheries in medieval and early modern times.

Ancient DNA from the extinct New Zealand grayling (Prototroctes oxyrhynchus) reveals evidence for Miocene marine dispersal

The evolutionary history of Southern Hemisphere graylings (Retropinnidae) in New Zealand (NZ), including their relationship to the Australian grayling, is poorly understood. The NZ grayling (Prototroctes oxyrhynchus) is the only known fish in NZ to have gone extinct since human arrival there. Despite its historical abundance, only 23 wet and dried, formalin-fixed specimens exist in museums. We used high-throughput DNA sequencing to generate mitogenomes from formalin-fixed P. oxyrhynchus specimens, and analysed these in a temporal phylogenetic framework of retropinnids and osmerids. We recovered a strong sister-relationship between NZ and Australian grayling (P. mareana), with a common ancestor ~13.8 Mya [95% highest posterior density (HPD): 6.1–23.2 Mya], after the height of Oligocene marine inundation in NZ. Our temporal phylogenetic analysis suggests a single marine dispersal between NZ and Australia, although the direction of dispersal is equivocal, followed by divergence into genetically and morphologically distinguishable species through isolation by distance. This study provides further insights into the possible extinction drivers of the NZ grayling, informs discussion regarding reintroduction of Prototroctes to NZ and highlights how advances in palaeogenetics can be used to test evolutionary hypotheses in fish, which, until relatively recently, have been comparatively neglected in ancient-DNA research.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.