Abstract
Ancient DNA (aDNA) approaches have been successfully used to infer the long-term impacts of climate change, domestication, and human exploitation in a range of terrestrial species. Nonetheless, studies investigating such impacts using aDNA in marine species are rare. Atlantic cod (Gadus morhua), is an economically important species that has experienced dramatic census population declines during the last century. Here, we investigated 48 ancient mitogenomes from historical specimens obtained from a range of archeological excavations in northern Europe dated up to 6,500 BCE. We compare these mitogenomes to those of 496 modern conspecifics sampled across the North Atlantic Ocean and adjacent seas. Our results confirm earlier observations of high levels of mitogenomic variation and a lack of mutation-drift equilibrium—suggestive of population expansion. Furthermore, our temporal comparison yields no evidence of measurable mitogenomic changes through time. Instead, our results indicate that mitogenomic variation in Atlantic cod reflects past demographic processes driven by major historical events (such as oscillations in sea level) and subsequent gene flow rather than contemporary fluctuations in stock abundance. Our results indicate that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on mitogenomic diversity in a wide-spread marine species with high gene flow such as Atlantic cod. These observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and the importance of genetic variation in implementing conservation strategies. Instead, these observations imply that any measures toward the demographic recovery of Atlantic cod in the eastern Atlantic, will not be constrained by recent loss of historical mitogenomic variation.
Introduction
Continuous human activities and a changing climate have influenced terrestrial and marine ecosystems for millennia (Venter et al., 2016; Rodrigues et al., 2019; Mitchell and Rawlence, 2021), impacting the evolutionary potential and population demography of a range of species (Seersholm et al., 2018). Ancient mitochondrial DNA (mtDNA) has been widely used to understand long-term genomic consequences of such impacts (Shapiro et al., 2004; Nyström et al., 2006; Stiller et al., 2010; Paijmans et al., 2013; Fortes and Paijmans, 2015; Casas-Marce et al., 2017). Nonetheless, most ancient mtDNA studies have focused on terrestrial species, and studies that investigate the impacts of long-term human activities and/or climatic variation on fish, using whole genome sequencing approaches, are relatively rare. Long-term commercial fisheries—covering many centuries—have contributed to the decline of economically and ecologically important marine species (Exadactylos et al., 2007; Pinnegar and Engelhard, 2008; Barrett, 2019). The consequences of intensive fishing in recent times may be difficult to assess as this requires an understanding of historical population dynamics (Selim et al., 2016). The analysis of long-term biological and demographic fluctuations can therefore help to improve guidelines for sustainable fisheries management and optimal conservation measures (Barrett, 2019). In order to provide a long-term perspective on fishing exploitation impacts, the use of archeological evidence, such as fish bone remains, is essential for those periods for which little or no historical data are available. Recent developments in whole genome aDNA methods now allow the inference of demographic histories and the estimation of genetic fluctuations over time from fishbone samples (Oosting et al., 2019; Ferrari et al., 2021). Such combined molecular analyses of historical and modern samples can potentially provide an understanding of the association between human-environmental impact and population declines (Hofman et al., 2015).
Several studies have shown the utility of temporal mtDNA analyses in the marine environment. For instance, ancient mitogenomes have investigated impacts of climate and hunting on the Atlantic walrus (Star et al., 2018; Keighley et al., 2019; Barrett et al., 2020), narwhals (Louis et al., 2020), and the extinct great auk (Thomas et al., 2019). In fish, such studies remain limited to partial mitogenome data. For example, a shift in sturgeon species distributions was detected during the Holocene in the North East Atlantic based on CytB amplicon data (Nikulina and Schmölcke, 2016). Moreover, impacts of habitat destruction and human activities during the 1800s were associated with a reduction of the mtDNA diversity of Chinook salmon from the Columbia River in the 12S and control region by comparing ancient and modern samples (Johnson et al., 2018). Similarly, impacts of human exploitation and climate oscillations were associated with losses of haplotypic CytB variation in Atlantic cod during the 15th to 16th centuries in Iceland (Olafsdottir et al., 2014). In contrast, comparable levels of ancient mtDNA genetic diversity were found between ancient and modern samples of herring specimens, despite continuous human exploitation (Speller et al., 2012). Notwithstanding these examples, human-environmental impacts and population declines remain unclear for a wide range of marine species and populations.
Atlantic cod (Gadus morhua L. 1758) is a benthopelagic predatory fish with high reproductive rates and with a fundamental ecological role in marine ecosystems (Barth et al., 2017; Edvardsson et al., 2019). It has been one of the most exploited fish species in the North Atlantic Ocean (Carr et al., 1995; Árnason et al., 2000; Nicholls et al., 2021). The distribution of this species extends through the cold waters of North America, across the continental shelves of Greenland and Iceland, and northern Europe (Lait et al., 2018). Relatively large population sizes have been characteristic throughout its entire distribution even during the expansion of long-distance fish trading during the 12th to 13th centuries in the eastern Atlantic and at the beginning of the 16th century in the western Atlantic (Barrett et al., 2004, 2011; Orton et al., 2014; Castañeda et al., 2020). However, intensive fishing activities during the 20th century (Mieszkowska et al., 2009; Jonsson et al., 2016; Brattey et al., 2018) resulted in the severe depletion of several stocks, for instance the North Sea stock, which was decimated from annual landings of 354,000 to 50,000 tons during this period (Bannister, 2004). In addition to past human exploitation, climatic events like the Little Ice Age—a cooling period that varied regionally in timing and duration but occurred between ca. 1300–1850 CE—may have caused large declines between the sixteenth and 17th centuries (Edvardsson et al., 2019).
The genomic consequences of such population dynamics and declines in Atlantic cod remain unclear. Based on partial and whole mtDNA data, Atlantic cod populations between the western and eastern Atlantic Ocean show significant structure (Árnason, 2004; Jørgensen et al., 2018; Lait et al., 2018), whereas low to no mtDNA differentiation has been found across a wide range of eastern Atlantic locations (Carr et al., 1995; Árnason and Palsson, 1996; Árnason et al., 1998, 2000; Sigurgíslason and Árnason, 2003). Here, we compared modern and ancient Atlantic cod mitogenomes—dated up to 6500 BCE—from different fishing locations in northern Europe. We evaluated whether Atlantic cod in the eastern Atlantic have experienced any loss of genetic variation, analyzed long term patterns of effective population size, and related any observed decline to the impact of commercial fisheries or climate change.
Materials and Methods
Sample Collection
Ancient samples of Atlantic cod (n = 48) were obtained from 11 excavation sites (Figure 1 and Supplementary Table 1) and were stored dry and unfrozen. The specimens were all supplied by the relevant archeological organizations, or sampled with permission on their premises. The shipment of Atlantic cod bones does not require CITES or other wildlife regulation permits for transport or analysis. Where practicable, only a subsample of bone was employed for the aDNA research, leaving material for other studies. Dating of the samples (Supplementary Table 1) was based on archeological context. Ancient samples were morphologically and genetically identified as Atlantic cod. A total of 472 available modern mitogenomes were obtained from Jørgensen et al. (2018), Lait et al. (2018), and Barth et al. (2019). Novel mtDNA sequence data from modern specimens sampled in 2016 in Orkney, United Kingdom (n = 24) were also included (Figure 1 and Supplementary Table 2). The collection of the Orkney specimens complied with the Nagoya Protocol and Convention on Biological Diversity, which the United Kingdom signed up to in 2016. All specimens were deceased when the fin clip was collected.
FIGURE 1
DNA Extraction, Amplification and Sequencing
DNA extraction and library preparation from ancient samples were performed in the aDNA laboratory at the University of Oslo under rigorous measures (Cooper and Poinar, 2000; Gilbert et al., 2005). All ancient samples were processed with the same DNA extraction and library protocols according to Ferrari et al. (2021). In short, bones were UV-treated for 10 min per side and pulverized using a stainless-steel mortar (Gondek et al., 2018). Per specimen, two aliquots containing between 150 and 200 mg of bone powder were used as starting material for DNA extraction. Double-indexed blunt-end sequencing libraries were built from 15 to 16 μl of DNA extract using the Meyer-Kircher protocol (Meyer and Kircher, 2010; Kircher et al., 2012) with the modifications listed in Schroeder et al. (2015) and the single-tube (BEST) protocol (Carøe et al., 2018) with the modifications described in Mak et al. (2017). Sequencing reads were processed using PALEOMIX v1.2.13 (Schubert et al., 2014). Trimming of residual adapter contamination, filtering and collapse of reads was done using AdapterRemoval v.2.1.7 (Lindgreen, 2012). Sequencing reads shorter than 25 bp were discarded. Mapping of remaining reads was performed against the Atlantic cod GadMor3.0 nuclear genome (RefSeq assembly accession GCF_902167405.1; Star et al., 2011; Tørresen et al., 2017) and mitochondrial genome (Johansen and Bakke, 1996) using BWA v.0.7.12 (Li and Durbin, 2009) with the aln algorithm, disabled seeding and minimum quality score of 25. The resulting BAM files were indexed with samtools v.1.9 (Li et al., 2009) and DNA postmortem damage assessed using MapDamage v.2.0.9 (Jónsson et al., 2013). DNA from modern Orkney samples were extracted using a DNeasy Blood & Tissue kit (Qiagen). Libraries were assembled with a TrueSeq DNA PCR-Free Preparation Kit and sequenced on an Illumina HiSeq 2,500. Modern alignment—including Orkney and Barth et al. (2019) samples, and the outgroup Alaska pollock (Gadus chalcogrammus; Malmstrøm et al., 2016) – was carried out using BWA v.0.7.12 with the mem algorithm, and a minimum quality score of 25.
Mitogenomic Analysis
Variant calling was performed using GATK v.4.1.4. (McKenna et al., 2010) simultaneously in all ancient, modern Orkney and Barth et al. (2019) samples, including the outgroup. gVCF files were created for each sample using GATK HaplotypeCaller (with ploidy set to 1). Individual genotypes were combined in one file using GATK CombineGVCFs and GenotypeGVCFs. Filtering was performed using bcftools v.1.9. (Li et al., 2009) and vcftools v.0.1.16. (Danecek et al., 2011) with the following thresholds: FS < 60.0, SOR < 4, MQ > 30.0, QD > 2.0, SnpGap = 10, minGQ = 15, minDP = 3, remove indels = yes, biallelic loci = yes, meanDP < 30 and read depth > 3. Consensus sequences were built using bcftools consensus and aligned using MAFFT v.7.429 (Katoh and Standley, 2013). Available modern mitogenomes obtained from Jørgensen et al. (2018) and Lait et al. (2018) were manually inspected using MEGA v.7 (Kumar et al., 2016) to set as missing the consistent nucleotide differences (between 50 and 100%; Supplementary Table 3) between their Illumina, Sanger and/or Roche 454 technologies with the Illumina sequenced mitogenomes in this study. Control region and half of the tRNA-Pro sequences from all the mitogenomes were excluded from further analyses as these two regions were not fully complete (i.e., 15,696–15,815 positions) after aligning sequences obtained from Jørgensen et al. (2018) and Lait et al. (2018) with the sequences presented in this study and Barth et al. (2019) samples. Thus, all sequences analyzed had 15,695 bp in length. Validated SNPs were annotated as transversion and/or transition using SNP-sites (Page et al., 2016). Checked and modified modern sequences (Jørgensen et al., 2018; Lait et al., 2018) were added and aligned to our multi-fasta alignment using MAFFT v.7.429. Unique sequences were identified with IQTREE v.1.6.12 (Nguyen et al., 2015). File formats required for different software and/or packages were obtained with seqinr and ape (i.e., nexus format; Paradis and Schliep, 2019; Charif et al., 2020), and phylotools (i.e., phylip format; Zhang et al., 2017) packages implemented in R.
Different sample combinations were used to compare the genetic diversity of the ancient samples to those of the modern conspecifics. Given the low spatial structure in the eastern Atlantic region (Árnason and Palsson, 1996; Árnason et al., 1998; Sigurgíslason and Árnason, 2003) and lack of consistent spatial structure amongst specimens (Supplementary Figures 3, 5, 6), all 48 ancient samples were compared as a single group to modern samples grouped into larger marine locations (according to their geographical proximity or ecotype; Figure 1). In addition, a comparison of subsets of multiple specimens from two archeological locations (Quoygrew and Haithabu) for which a more specific temporal pair from the same geographical region could be identified, was performed (Supplementary Table 1). Quoygrew specimens were locally sourced (Harland and Barrett, 2012; Star et al., 2017). Therefore, modern specimens sampled in the same area (i.e., modern Orkney) provide a logical, spatially consistent temporal comparison. However, specimens from Haithabu, were sourced from northern Norway (Star et al., 2017), and belonged to the North East Arctic ecotype. For these traded specimens, the North East Arctic ecotypes provide a spatially relevant temporal comparison, rather than North Sea or western Baltic specimens.
Haplotype (h) and nucleotide diversities (π), number of haplotypes (Nh) and number of polymorphic sites (S) were calculated using DnaSP v.6 (Rozas et al., 2017). To allow direct comparison with earlier CytB results (Árnason, 2004; Olafsdottir et al., 2014; Jørgensen et al., 2018), specific CytB haplotypes based on 250 bp gene fragment as previously reported by Árnason (2004) were identified using MEGA v.7. Demographic histories were determined by Tajima’s D (TD) and Fu’s F (F) neutrality in DnaSP v.6. A different number of specimens were obtained for ancient and modern locations. We corrected for such differences in sample size by randomly downsampling the modern specimens for each of the temporal comparisons (North East Arctic and Orkney) using 1,000 bootstrap replicates. A 95% confidence interval of the genetic parameters; genetic variation (π) and patterns of population demography (TD and F) was calculated from these 1,000 bootstrap replicates that were sampled using a without replacement approach with the sample function implemented in R (R Core Team, 2020) and the fasta.sample function in the FastaUtils package also in R (Salazar, 2020). For the bootstrapping test, π, TD and F from temporally spaced modern locations were re-calculated with the pegas (Paradis, 2010) and PopGenome (Pfeifer et al., 2020) packages implemented in R. Relationships among ancient and modern samples were visualized for whole mitogenome and CytB sequence data, by constructing a mitochondrial haplotype-genealogy graph using Fitchi (Matschiner, 2016) with the ML-based phylogenetic tree obtained with IQTREE v.1.6.12 as input.
Population Dynamics and Demographic Reconstruction
An analysis of molecular variance (AMOVA, 1,000,000 permutations) and population pairwise genetic distances (ΦST) were obtained in Arlequin v.3.5 (Excoffier and Lischer, 2010), to determine the distribution of variation between marine locations and temporally spaced locations. Divergence and coalescent analyses were based on unique sequences only (n = 525 sequences including the outgroup). Substitution model selection for unique sequences was performed using PHYML v.3.1 (Guindon et al., 2010) as implemented in JMODELTEST v.2.1.10 (Guindon and Gascuel, 2003; Darriba et al., 2012). Model selection was determined on the following partitions: 1st, 2nd, and 3rd codons from protein coding regions, rRNAs and tRNAs. Best-fitting models were selected according to the Aikaike Information Criterion (AIC; Supplementary Table 4). Based on these results, phylogenetic estimates were obtained using BEAST v.2.6.3 (Bouckaert et al., 2019).
Bayesian settings for all phylogenetic analyses included two sets of partitions: coding region and non-coding region. Three independent runs to test for chain convergence were run under the Coalescent Constant Population Tree Prior. Tip ages (ancient and modern dates) were included for each set of runs (Supplementary Tables 1, 2). Sample dates for ancient specimens were rounded to a midpoint date—from a given range—where necessary. To achieve high effective sample sizes (ESS = ≥ 200), chain lengths were run 800,000,000 under a substitution rate of 1.14 × 10–8 substitution/site/year as per Lait (2016) assuming a GTR + I (for coding regions) and TIM1 + I (for non-coding regions) models of evolution and a strict clock. Tracer v.1.71 (Rambaut et al., 2018) was used to check for convergence of MCMC and to ensure sufficient sampling. Consensus trees were obtained using TreeAnnotator v.2.6.2—implemented in BEAST v.2.6.3—after 10% burn-in. Final phylogenetic trees were viewed and edited in FigTree v.1.4.4.
Finally, a Coalescent Bayesian Skyline (CBS) analysis was completed to reconstruct the demographic history—including female effective population size (Ne) —of Atlantic cod through time. To assess any confounding effect of past or contemporary population structure (Heller et al., 2013), we analyzed demographic history using 6 different data sets (excluding the outgroup): (I) all 524 sequences, (II) 476 modern sequences (excluding 48 ancient samples), (III) 273 sequences (excluding clades associated with most western Atlantic and Baltic Sea samples), (IV) 368 sequences (excluding the clade associated with most Baltic Sea samples), (V) 429 sequences (excluding clades associated with most western Atlantic samples) and (VI) 48 ancient sequences (excluding all modern samples). The specific clades that were excluded in III, IV and V can be found in Supplementary Figure 4. We used the same MCMC sampling procedure described before with 3 independent runs reaching convergence at high effective sample sizes (ESS = ≥ 200). Chain lengths were run 800,000,000 for data sets I, II and V with a number of bPopSize and bGroupSize of 10; while chain length for data sets III and IV were run 500,000,000 and 50,000 for data set VI with a number of bPopSize and bGroupSize of 5.
Results
Mitogenomic Variation
Sequencing reads from all ancient specimens showed the expected patterns of DNA fragmentation and deamination rates that were consistent with those of authentic aDNA (Supplementary Figure 1). We obtained 48 mitogenomes with at least 3-fold average coverage. We also obtained mitogenomes for 24 modern Orkney specimens (Supplementary Table 2). A total of 2135 SNPs (∼13% of mitogenome positions) were identified among all 545 samples – including the outgroup species Alaska pollock –: 1219 SNPs corresponded to informative sites and 916 SNPs were singletons (Supplementary Table 5).
Nucleotide diversity (π) between modern locations ranged between 0.002 and 0.003 (Table 1) and π of ancient samples did not vary from the values obtained in modern locations. The temporal comparison of specific sites (Quoygrew-Orkney and Haithabu-North East Arctic), showed limited significant differences between genetic statistics of temporally spaced ancient and modern locations (Supplementary Table 6 and Supplementary Figure 2), where Haithabu has significantly lower π and higher F compared to the North East Arctic (Supplementary Figure 2).
TABLE 1
| Location | Code | N | h | Nh | S | π | TD | F | |
| Modern | western Atlantic | WA | 124 | 1.000 | 124 | 759 | 0.002 | –2.652* | –5.668* |
| Baffin Island | BI | 18 | 0.791 | 7 | 68 | 0.002 | 1.570 | 1.688 | |
| western Baltic | WB | 43 | 0.996 | 40 | 221 | 0.002 | –1.719 | –2.992* | |
| eastern Baltic | EB | 36 | 1.000 | 36 | 249 | 0.002 | –1.888* | –3.305* | |
| Tvedestrand (fjord) | TV | 37 | 0.982 | 31 | 256 | 0.002 | –1.570 | –2.757* | |
| North Sea | NS | 99 | 0.999 | 96 | 678 | 0.003 | –2.350* | –4.421* | |
| North East Arctic | NE | 97 | 0.999 | 92 | 716 | 0.003 | –2.408* | –4.695* | |
| Norwegian coast | NC | 41 | 1.000 | 41 | 377 | 0.002 | –2.187* | –3.694* | |
| Ancient | Ancient | AN | 48 | 0.998 | 46 | 364 | 0.002 | –2.283* | –3.906* |
Estimates of genetic diversity statistics for Atlantic cod at whole mitogenomes from different marine locations or ecotypes in the North Atlantic (see Supplementary Tables 1, 2, 5).
Significant values are indicated with * (p ≤ 0.01). N, sample size; h, haplotype diversity; Nh, number of haplotypes; S, number of polymorphic sites; π, nucleotide diversity; TD, Tajima’s D; F, Fu’s F. Locations with 1 sample are excluded for genetic analysis (i.e., Irish Sea = IS).
Neutrality tests showed significant negative values for all Tajima’s D (TD) and F statistics in most locations, except for the western location Baffin Island, and the eastern locations Tvedestrand fjord and western Baltic (Table 1). Overall, there were 486 haplotypes —- including the outgroup—across all 545 samples, of which only 26 were shared between individuals (Figure 2 and Supplementary Table 7). Ancient CytB variation consisted of 7 different haplotypes, including four main haplotypes (A, C, D, and E) previously identified in modern mtDNA studies (Árnason, 2004; Jørgensen et al., 2018). Two novel variations of existing CytB haplotypes were found in western Baltic (haplotype ED) and North Sea (haplotype LI), while another 2 novel variations of existing CytB haplotypes were found among ancient samples (haplotypes LJ and TI). The most prevalent ancient haplotypes were A and E (∼ 40 and 38%, respectively, Supplementary Tables 1, 7), which were also commonly found in modern samples (Supplementary Table 8). The haplotype genealogy for whole mitogenome and CytB sequence data showed an extensive distribution of ancient samples across marine locations (Figure 2 and Supplementary Figure 4). Limited geographic mitogenome structure was observed, except for elevated divergence between western Atlantic and eastern Atlantic locations, and between locations in the western and eastern Baltic Sea and other eastern Atlantic locations (Figures 2B,C and Supplementary Figures 4B,C). A star-like topology is observed for the whole mitogenome and CytB genealogies (Figure 2 and Supplementary Figure 4).
FIGURE 2
Demographic Patterns and Population Structure
The AMOVA assigned 7.58% of the variation between marine locations (including ancient samples as a single group) while 91.47% of the variation was represented between individuals (ΦCT = 0.076, p ≤ 0.001; ΦST = 0.085, p ≤ 0.000). Pairwise ΦST values (Figure 3 and Supplementary Tables 9, 10) showed significant differentiation levels between all ancient samples and western Atlantic, Baffin Island, western Baltic, eastern Baltic and Tvedestrand fjord. Ancient samples showed higher differentiation when compared to western Atlantic (ΦST = 0.117), and Baffin Island (ΦST = 0.192) in comparison to other eastern Atlantic locations. Among modern samples, western Atlantic, Baffin Island, western Baltic and eastern Baltic showed significant ΦST values when compared to all other locations (Supplementary Tables 9, 10). ΦST values were not significant between North Sea, North East Arctic, Norwegian coast and Ancient samples. Pairwise ΦST values between temporal spaced locations also showed no significant differentiation (Quoygrew and modern Orkney: ΦST = 0.000; p = 0.807; and Haithabu and North East Arctic: ΦST = 0.000; p = 0.456).
FIGURE 3
The time-calibrated Bayesian phylogeny for ancient and modern Atlantic cod samples resulted in 2 main clades with an estimated divergence from the most recent common ancestor at 220 kya (95% highest posterior density (HPD) = 194,780—249,980 kya; Figure 4). The first clade, which is not further divided, includes mitogenomes from 6 different widely scattered localities. The second clade was composed by 16 subclades with posterior probability > 0.8, with divergence times of ca. 100 kya. Clades and subclades in the phylogeny were not geographically structured, with the exception of most samples from western Atlantic, and most samples from western and eastern Baltic, which clustered together (Figures 2, 4).
FIGURE 4
The Bayesian skyline analysis using different subsets of the data revealed a consistent pattern of step-wise population expansions followed by periods of constant population size. Expansions around 150, 50, and 10 kya are present in most subsets (Figure 5). A population expansion of Atlantic cod was identified ca. 50 kya in all subsets. The most recent expansion identified (around 10 kya), is only present in data sets that include clades with most Baltic Sea specimens (Figures 5A,B,E). Despite such differences, all analyses agree with a high and increasing female effective population size (Ne) of Atlantic cod (Ne = ca. 1,000,000—10,000,000) during the last ca. 100 kya, with highest estimates of Ne during the last few millennia (Figure 5).
FIGURE 5
Discussion
Here, we compared modern and ancient mtDNA diversity in Atlantic cod to investigate whether observed historical and contemporaneous census population declines (Hutchinson et al., 2003; Hylen et al., 2008; Limburg et al., 2008; Bartolino et al., 2012; Jonsson et al., 2016; Brattey et al., 2018) have had mitogenomic consequences. The temporal comparison of 48 ancient specimens to 496 modern conspecifics did not reveal consistent significant mitogenomic changes or measurable effective genetic population declines through time. Below, we discuss reasons why such genomic impacts may not be observed.
First, mitogenomic variation is high in modern Atlantic cod and is characterized by limited genetic differentiation between populations and incomplete lineage sorting over large spatial scales across its range in the North Atlantic (Jørgensen et al., 2018; Lait et al., 2018). Low observed genetic differentiation (ΦST) between Tvedestrand fjord and other Norwegian coastal locations, as well as between the North Sea, the North East Arctic and the Norwegian coast confirm this lack of geographic structuring over large parts of the eastern Atlantic (Figure 3). Indeed, the non-significant differentiation of all ancient samples with modern North Sea, North East Arctic and Norwegian coast is fully consistent with their presumed geographical origin and highlights the long-term lack of mtDNA structure in this region. Non-significant ΦST values between the Norwegian coastal locations and Tvedestrand fjord indicate possible recent migration of fish between such coastal communities and more restricted fjord populations (Knutsen et al., 2011). Compared to many terrestrial ecosystems, where populations can often be isolated by physical barriers—which restrain interbreeding and dispersal—(Hauser and Carvalho, 2008; Exadactylos et al., 2019), in marine ecosystems the absence of physical barriers promotes larger panmictic populations and Atlantic cod is no exception (Berg et al., 2016, 2017; Sodeland et al., 2016; Barth et al., 2017). Thus, a combination of low spatial resolution of mtDNA data as a result of continuous gene flow and connectivity may mask any local temporal erosion of mitogenomic diversity (Welch et al., 2012) in Atlantic cod.
Second, we determined high long-term estimates of effective population size (Ne = ca. 1,000,000–10,000,000;Figure 5), which is in agreement with earlier observations in Atlantic cod (Hardie et al., 2006; Therkildsen et al., 2010; Pinsky et al., 2021). Estimates of Ne can remain high in economically important fish species, even if their populations have experienced a large biomass decline (Hauser and Carvalho, 2008) since it takes hundreds of generations (i.e., depending on the generation time of the species; Amos and Balmford, 2001; Frankham et al., 2002) for the actual population numbers and breeding populations to be reflected in Ne (Hauser and Carvalho, 2008). In fact, simulations have shown that a population with theoretical Ne of 100 (which is several orders of magnitude lower than observed in Atlantic cod) would retain 75% of heterozygosity after 57 generations (Frankham et al., 2002; Welch et al., 2012). Given that such population declines take a very long time to lead to measurable genomic consequences, mtDNA—as a single locus—will have limited power to record such changes in populations of high Ne (Allentoft et al., 2014; Johnson et al., 2018; Thomas et al., 2019; Spencer, 2020). The absence of significant genetic changes in this study is consistent with the absence of such changes in genome-wide data using historical samples of Atlantic cod from the western and eastern Atlantic (Pinsky et al., 2021) and with the absence of such changes in mitogenomic data from other taxa that have similarly high estimates of Ne as Atlantic cod, such as the Pacific herring (Speller et al., 2012; Moss et al., 2016), the Hawaiian petrel (Welch et al., 2012) and even extinct species such as the New Zealand moa (Allentoft et al., 2014), the passenger pigeon (Murray et al., 2017) and the great auk (Thomas et al., 2019).
In contrast, temporal losses of mitogenomic diversity and/or declines in Ne have been reported in species that have suffered population fragmentation (e.g., resulting in small effective population sizes) or that have experienced limited connectivity, such as the steppe bison (Shapiro et al., 2004), the Scandinavian arctic fox (Nyström et al., 2006), cave bears (Stiller et al., 2010), the Iberian lynx (Casas-Marce et al., 2017), the Iberian salmon (Consuegra et al., 2002) and the common bream (Ciesielski and Makowiecki, 2005). Interestingly, a loss of haplotypic variation has been identified—using CytB sequence data—for a single period (i.e., 15th to 16th centuries, out of 6 temporal periods investigated) in an Icelandic population of Atlantic cod (Olafsdottir et al., 2014). There are two potential explanations for this discrepancy. First, nearly all substitutions that comprise the CytB haplotypes can be affected by post-mortem deamination (i.e., they consist of C > T and G > A substitutions). Most of the ancient sequences (90%) investigated in Olafsdottir et al. (2014) were obtained in a single round of PCR without evaluation of such post-mortem deamination. Therefore, such bias due to post-mortem damage cannot be excluded. Second, our sampling does not include many specimens from Iceland (Figure 1 and Supplementary Table 1), and it remains possible that—with 156 samples—a local effect has been observed in Olafsdottir et al. (2014), which we do not detect in our data.
Third, we do not observe major novel mtDNA lineages in the ancient data, nor observe a significant loss of such lineages over time. Instead, the majority of Atlantic cod mtDNA lineages observed in ancient and modern samples today have originated ca. 100–150 kya (Figure 4), during a period of population expansion (Figure 5). Therefore, the gain of such lineages—and associated population expansions—in Atlantic cod is more likely caused by changes in abundance driven by major historical climatic events such as eustatic oscillations in sea level, and the interglacial and warming periods experienced during the last glacial maximum ca. 23,000 kya (Bigg et al., 2008) and the Wisconsinan (ca. 110–120 kya) and Illinoian (ca. 200–130 kya) glaciations (Gibbard and Van Kolfschoten, 2005) as described by Lait et al. (2018). For instance, we only observe the most recent population expansion ca. 10 kya (Figures 5A,B,E) when including those mtDNA clades which are strongly associated with the Baltic Sea. The timing of this expansion is in agreement with the development of the Baltic Sea (ca. 7,000–8,000 years; Ojaveer et al., 2010; Wenne et al., 2020) which has led to genetically distinct Atlantic cod populations that have adapted to local environmental conditions (i.e., salinity and temperature; Johannesson and Andre, 2006; Berg et al., 2015; Wenne et al., 2020). Therefore, the observed changes in Ne reflect past population demography rather than recent and contemporary demographic changes (Lombal et al., 2020).
It is clear from zooarcheological evidence that Atlantic cod has periodically experienced intense exploitation in the distant past, particularly around the North Sea and the Baltic Sea (Barrett et al., 1999; Enghoff, 1999; Olson and Walther, 2007; Orton et al., 2011). This fishing pressure became even greater in the 19th and 20th centuries (e.g., Thurstan et al., 2010). Landings of Atlantic cod exceeded 4,000,000 tons during 1960–1990s in the North Atlantic Ocean (Shelton and Morgan, 2014). In particular, landings surpassed 600,000 tons in Iceland by ca. 1930s (Drinkwater, 2006), 354,000 tons in the North Sea during ca. 1970s (Bannister, 2004), 200–400,000 tons in the eastern Baltic during 1960–1990s (MacKenzie et al., 2002), 650,000 tons in North East Arctic between 1937 and 1938 up to 800–1,200,000 tons in ca. 1950s (Sætersdal and Hylen, 1964; Hylen, 2002). Such high levels of exploitation led to major reductions in present abundances of most Atlantic cod populations (i.e., Food and Agriculture Organization [Fao], 2020-2021a,b). Nonetheless, for the reasons discussed above, our results indicate that such population declines of Atlantic cod did not lead to a detectable impact on the mtDNA genome on the time scale we investigated here.
Taken together, our results highlight that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on the ancient mitogenomic diversity of a wide-spread marine species with high gene flow such as Atlantic cod. Future ancient DNA studies should consider the inclusion of nuclear genomic data and extensive sampling on a local scale—considering a temporal comparison of specimens from the same geographical region—to assess the effects of climate and human exploitation with greater statistical power. Finally, our observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and they do not oppose information about the important implications of genetic variation in evolutionary biology, ecology and conservation biology. Instead, our observations suggest that conservation management measures aimed toward the demographic recovery of Atlantic cod in the eastern Atlantic, if achievable by conservation management measures, will not be constrained by recent loss of historical mitogenomic variation.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/ena, PRJEB42959.
Author contributions
BS and JB: conceptualization, project design, and supervision. GF and LM-G: laboratory work and data curation. LM-G with input from GF, TO, LA, BS, and JB: formal analysis. HB and MS: analytical advice. JB, IY, IJ, JH, RN, DO, SH-D, US, AH, and RB: ancient Atlantic cod specimens. FB: modern Orkney specimens. JB, IY, IJ, JH, RN, DO, SH-D, US, AH, and RB: archeological context information. LM-G and BS: data visualization. JB, BS, KJ, SJ, and DO: funding acquisition. LM-G and BS with input from JB: writing—original draft. All authors writing—review and editing.
Funding
This work was supported by Research Council of Norway projects “Catching the Past” (262777) and “The Aqua Genome Project” (221734), Leverhulme Trust Project MRF-2013-065 and the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie grant agreement No. 813383.
Acknowledgments
We thank Agata T. Gondek-Wyrozemska for help with processing the ancient specimens and Francis Neat for providing modern Orkney specimens. We also thank Oliver Kersten and Michael Matschiner for advice during analyses. Finally, we thank M. Skage, S. Kollias and A. Tooming-Klunderud at the Norwegian Sequencing Centre for sequencing and processing of samples. Analyses were performed on the SAGA Cluster using the resources and assistance from the SIGMA2 Metacenter, the Norwegian National Infrastructure for High Performance Computing and Data Storage.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2021.671281/full#supplementary-material
References
1
AllentoftM. E.HellerR.OskamC. L.LorenzenE. D.HaleM. L.GilbertM. T. P.et al (2014). Extinct New Zealand megafauna were not in decline before human colonization.Proc. Natl. Acad. Sci.U. S. A.1114922–4927. 10.1073/pnas.1314972111
2
AmosW.BalmfordA. (2001). When does conservation genetics matter?Heredity87257–265. 10.1046/j.1365-2540.2001.00940.x
3
ÁrnasonE. (2004). Mitochondrial cytochrome b DNA variation in the high-fecundity Atlantic cod: trans-Atlantic clines and shallow gene genealogy.Genetics1661871–1885. 10.1093/genetics/166.4.1871
4
ÁrnasonE.PalssonS. (1996). Mitochondrial cytochrome b DNA sequence variation of Atlantic cod Gadus morhua, from Norway.Mol. Ecol.5715–724. 10.1111/j.1365-294x.1996.tb00368.x
5
ÁrnasonE.PetersenP. H.KristinssonK.SigurgíslasonH.PálssonS. (2000). Mitochondrial cytochrome b DNA sequence variation of Atlantic cod from Iceland and Greenland.J. Fish Biol.56409–430. 10.1006/jfbi.1999.1167
6
ÁrnasonE.PetersenP. H.PálssonS. (1998). Mitochondrial cytochrome b DNA sequence variation of Atlantic cod, Gadus morhua, from the Baltic and the White Seas.Hereditas12937–43. 10.1111/j.1601-5223.1998.00037.x
7
BannisterR. C. A. (2004). “The rise and fall of cod (Gadus morhua, L.) in the North Sea,” inManagement of Shared Fish Stocks, edsPayneA. I. L.O’BrienC. M.RogersS. I. (Oxford: Blackwell Publishing), 316–338. 10.1002/9780470999936.ch19
8
BarrettJ. H. (2019). An environmental (pre) history of European fishing: past and future archaeological contributions to sustainable fisheries.J. Fish Biol.941033–1044. 10.1111/jfb.13929
9
BarrettJ. H.BoessenkoolS.KnealeC. J.O’ConnellT. C.StarB. (2020). Ecological globalisation, serial depletion and the medieval trade of walrus rostra.Quat. Sci. Rev.229:106122. 10.1016/j.quascirev.2019.106122
10
BarrettJ. H.LockerA. M.RobertsC. M. (2004). ‘Dark Age Economics’ revisited: the English fish bone evidence AD 600-1600.Antiquity78618–636. 10.1017/s0003598x00113262
11
BarrettJ. H.NicholsonR. A.Cerón-CarrascoR. (1999). Archaeo-ichthyological Evidence for Long-term Socioeconomic Trends in Northern Scotland: 3500 BC to AD 1500.J. Archaeol. Sci.26353–388. 10.1006/jasc.1998.0336
12
BarrettJ. H.OrtonD.JohnstoneC.HarlandJ.Van NeerW.ErvynckA.et al (2011). Interpreting the expansion of sea fishing in medieval Europe using stable isotope analysis of archaeological cod bones.J. Archaeol. Sci.381516–1524. 10.1016/j.jas.2011.02.017
13
BarthJ. M.BergP. R.JonssonP. R.BonanomiS.CorellH.Hemmer−HansenJ.et al (2017). Genome architecture enables local adaptation of Atlantic cod despite high connectivity.Mol. Ecol.264452–4466. 10.1111/mec.14207
14
BarthJ. M. I.Villegas-RíosD.FreitasC.MolandE.StarB.AndréC.et al (2019). Disentangling structural genomic and behavioural barriers in a sea of connectivity.Mol. Ecol.281394–1411. 10.1111/mec.15010
15
BartolinoV.CardinaleM.SvedängH.Linderholm, HansW.CasiniM.et al (2012). Historical spatiotemporal dynamics of eastern North Sea cod.Can. J. Fish. Aquat. Sci.69833–841. 10.1139/f2012-028
16
BergP. R.JentoftS.StarB.RingK. H.KnutsenH.LienS.et al (2015). Adaptation to Low Salinity Promotes Genomic Divergence in Atlantic Cod (Gadus morhua L.).Genome Biol. Evol.71644–1663. 10.1093/gbe/evv093
17
BergP. R.StarB.PampoulieC.BradburyI. R.BentzenP.HutchingsJ. A.et al (2017). Trans-oceanic genomic divergence of Atlantic cod ecotypes is associated with large inversions.Heredity119418–428. 10.1038/hdy.2017.54
18
BergP. R.StarB.PampoulieC.SodelandM.BarthJ. M. I.KnutsenH.et al (2016). Three chromosomal rearrangements promote genomic divergence between migratory and stationary ecotypes of Atlantic cod.Sci. Rep.6:23246. 10.1038/srep23246
19
BiggG. R.CunninghamC. W.OttersenG.PogsonG. H.WadleyM. R.WilliamsonP. (2008). Ice-age survival of Atlantic cod: agreement between palaeoecology models and genetics.Proc. R. Soc. B Biol. Sci.275163–173. 10.1098/rspb.2007.1153
20
BouckaertR.VaughanT. G.Barido-SottaniJ.DuchêneS.FourmentM.GavryushkinaA.et al (2019). BEAST 2.5: an advanced software platform for Bayesian evolutionary analysis.PLoS Comput. Biol.15:e1006650. 10.1371/journal.pcbi.1006650
21
BratteyJ.CadiganN.DwyerK. S.HealeyB. P.IngsD. W.LeeE. M.et al (2018). Assessment of the Northern Cod (Gadus morhua) stock in NAFO Divisions 2J3KL in 2016.DFO Can. Sci. Advis. Sec. (CSAS) Res. Doc. 2018/01810720–24.
22
CarøeC.GopalakrishnanS.VinnerL.MakS. S.SindingM. H. S.SamaniegoJ. A.et al (2018). Single−tube library preparation for degraded DNA.Methods Ecol. Evol.9410–419. 10.1111/2041-210x.12871
23
CarrS.SnellenA.HowseK.WroblewskiJ. (1995). Mitochondrial DNA sequence variation and genetic stock structure of Atlantic cod (Gadus morhua) from bay and offshore locations on the Newfoundland continental shelf.Mol. Ecol.479–88. 10.1111/j.1365-294x.1995.tb00194.x
24
Casas-MarceM.MarmesatE.SorianoL.Martínez-CruzB.Lucena-PerezM.NoceteF.et al (2017). Spatiotemporal dynamics of genetic variation in the Iberian lynx along its path to extinction reconstructed with ancient DNA.Mol. Biol. Evol.342893–2907. 10.1093/molbev/msx222
25
CastañedaR. A.BurliukC. M. M.CasselmanJ. M.CookeS. J.DunmallK. M.ForbesL. S.et al (2020). A Brief History of Fisheries in Canada.Fisheries45303–318. 10.1002/fsh.10449
26
CharifD.LobryJ. R.NecsuleaA.PalmeiraL.PenelS.PerriereG.et al (2020). Package ‘seqinr. URL: http://seqinr.r-forge.r-project.org/
27
CiesielskiS.MakowieckiD. (2005). Ancient and modern mitochondrial haplotypes of common bream (Abramis brama L.) in Poland.Ecol. Freshw. Fish14278–282. 10.1111/j.1600-0633.2005.00097.x
28
ConsuegraS.García, de LeánizC.SerdioA.González MoralesM.StrausL.et al (2002). Mitochondrial DNA variation in Pleistocene and modern Atlantic salmon from the Iberian glacial refugium.Mol. Ecol.112037–2048. 10.1046/j.1365-294x.2002.01592.x
29
CooperA.PoinarH. N. (2000). Ancient DNA: do it right or not at all.Science2891139–1139.
30
R Core Team. (2020). R: A language and environment for statistical computing.Vienna, Austria: R Foundation for Statistical Computing.
31
DanecekP.AutonA.AbecasisG.AlbersC. A.BanksE.DePristoM. A.et al (2011). The variant call format and VCFtools.Bioinformatics272156–2158. 10.1093/bioinformatics/btr330
32
DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and parallel computing.Nat. Methods9772–772. 10.1038/nmeth.2109
33
DrinkwaterK. F. (2006). The regime shift of the 1920s and 1930s in the North Atlantic.Prog. Oceanogr.68134–151. 10.1016/j.pocean.2006.02.011
34
EdvardssonR.PattersonW. P.BárðarsonH.TimsicS.ÓlafsdóttirG. Á (2019). Change in Atlantic cod migrations and adaptability of early land-based fishers to severe climate variation in the North Atlantic.Quat. Res.20191–11. 10.1017/qua.2018.147
35
EnghoffI. B. (1999). Fishing in the Baltic Region from the 5th century BC to the 16th century AD: evidence from Fish Bones.Archaeofauna841–85.
36
ExadactylosA.RigbyM. J.GeffenA. J.ThorpeJ. P. (2007). Conservation aspects of natural populations and captive-bred stocks of turbot (Scophthalmus maximus) and Dover sole (Solea solea) using estimates of genetic diversity.ICES J. Mar. Sci.641173–1181. 10.1093/icesjms/fsm086
37
ExadactylosA.VafidisD.TsigenopoulosC. S.GkafasG. A. (2019). High Connectivity of the White Seabream (Diplodus sargus, L. 1758) in the Aegean Sea, Eastern Mediterranean Basin.Animals9:979. 10.3390/ani9110979
38
ExcoffierL.LischerH. E. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows.Mol. Ecol. Resour.10564–567. 10.1111/j.1755-0998.2010.02847.x
39
FerrariG.CuevasA.Gondek-WyrozemskaA. T.BallantyneR.KerstenO.PálsdóttirA. H.et al (2021). The preservation of ancient DNA in archaeological fish bone.J. Archaeol. Sci.126:105317. 10.1016/j.jas.2020.105317
40
Food and Agriculture Organization [Fao]. (2020-2021a). International Council for the Exploration of the Sea (ICES). ICES Advice 2019. Cod - Western English Channel, Bristol Channel, Celtic Sea and Southwest of Ireland. FIRMS Reports. In: Fisheries and Resources Monitoring System (FIRMS) [online].Rome: Food and Agriculture Organization10.1016/j.jas.2020.105317
41
Food and Agriculture Organization [Fao]. (2020-2021b). International Council for the Exploration of the Sea (ICES). ICES Advice 2020. Cod - Baltic Sea (Eastern part). FIRMS Reports. In: Fisheries and Resources Monitoring System (FIRMS) [online].Rome: Food and Agriculture Organization.
42
FortesG. G.PaijmansJ. L. (2015). “Analysis of whole mitogenomes from ancient samples,” inWhole genome amplification, (Heidelberg: Springer), 179–195. 10.1007/978-1-4939-2990-0_13
43
FrankhamR.BallouJ. D.BriscoeD. A. (2002). “Chapter 10: Loss of genetic diversity in small populations,” inIntroduction to conservation genetics, ed.McInnessK. H. (Cambridge, UK: Cambridge University Press), 227–253. 10.1017/cbo9780511808999.013
44
GibbardP.Van KolfschotenT. (2005). “The Pleistocene and Holocene Epochs,” inA Geologic Time Scale 2004, edsSmithA. G.GradsteinF. M.OggJ. G. (Cambridge: Cambridge University Press), 441–452. 10.1017/cbo9780511536045.023
45
GilbertM. T. P.BandeltH.-J.HofreiterM.BarnesI. (2005). Assessing ancient DNA studies.Trends Ecol. Evol.20541–544. 10.1016/j.tree.2005.07.005
46
GondekA. T.BoessenkoolS.StarB. (2018). A stainless-steel mortar, pestle and sleeve design for the efficient fragmentation of ancient bone.BioTechniques64266–269. 10.2144/btn-2018-0008
47
GuindonS.DufayardJ.-F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: assessing the Performance of PhyML 3.0.Syst. Biol.59307–321. 10.1093/sysbio/syq010
48
GuindonS.GascuelO. (2003). A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.Syst. Biol.52696–704. 10.1080/10635150390235520
49
HardieD. C.GillettR. M.HutchingsJ. A. (2006). The effects of isolation and colonization history on the genetic structure of marine-relict populations of Atlantic cod (Gadus morhua) in the Canadian Arctic.Can. J. Fish. Aquat. Sci.631830–1839. 10.1139/f06-085
50
HarlandJ.BarrettJ. H. (2012). “Chapter 7: The Maritime Economy, Fish Bone,” inBeing an Islander: Production and identity at Quoygrew, Orkney, AD 900-1600, ed.BarrettJ. H. (Cambridge: McDonald Institute for Archaeological Research), 115–138.
51
HauserL.CarvalhoG. R. (2008). Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts.Fish Fish.9333–362. 10.1111/j.1467-2979.2008.00299.x
52
HellerR.ChikhiL.SiegismundH. R. (2013). The Confounding Effect of Population Structure on Bayesian Skyline Plot Inferences of Demographic History.PLoS One8:e62992. 10.1371/journal.pone.0062992
53
HofmanC. A.RickT. C.FleischerR. C.MaldonadoJ. E. (2015). Conservation archaeogenomics: ancient DNA and biodiversity in the Anthropocene.Trends Ecol. Evol.30540–549. 10.1016/j.tree.2015.06.008
54
HutchinsonW. F.OosterhoutC. V.RogersS. I.CarvalhoG. R. (2003). Temporal analysis of archived samples indicates marked genetic changes in declining North Sea cod (Gadus morhua).Proc. R. Soc. Lon. B Biol. Sci.2702125–2132. 10.1098/rspb.2003.2493
55
HylenA. (2002). Fluctuations in abundance of Northeast Arctic cod during the 20th century.ICES Mar. Sci. Symp.215543–550.
56
HylenA.NakkenO.NedreaasK. (2008). Northeast Arctic cod: fisheries, life history, stock fluctuations and management. Norwegian spring-spawning herring and Northeast Arctic cod.Trondheim: Tapir Academic Press, 83–118.
57
JohannessonK.AndreC. (2006). INVITED REVIEW: life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea.Mol. Ecol.152013–2029. 10.1111/j.1365-294x.2006.02919.x
58
JohansenS.BakkeI. (1996). The complete mitochondrial DNA sequence of Atlantic cod (Gadus morhua): relevance to taxonomic studies among codfishes.Mol. Mar. Biol. Biotechnol.5203–214.
59
JohnsonB. M.KempB. M.ThorgaardG. H. (2018). Increased mitochondrial DNA diversity in ancient Columbia River basin Chinook salmon Oncorhynchus tshawytscha.PLoS One13:e0190059. 10.1371/journal.pone.0190059
60
JónssonH.GinolhacA.SchubertM.JohnsonP. L.OrlandoL. (2013). mapDamage2. 0: fast approximate Bayesian estimates of ancient DNA damage parameters.Bioinformatics291682–1684. 10.1093/bioinformatics/btt193
61
JonssonP. R.CorellH.AndréC.SvedängH.MoksnesP. O. (2016). Recent decline in cod stocks in the North Sea–Skagerrak–Kattegat shifts the sources of larval supply.Fish. Oceanogr.25210–228. 10.1111/fog.12146
62
JørgensenT. E.KarlsenB. O.EmblemÅBreinesR.AndreassenM.RoungeT. B.et al (2018). Mitochondrial genome variation of Atlantic cod.BMC Res. Notes11:397. 10.1186/s13104-018-3506-3
63
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Mol. Biol. Evol.30772–780. 10.1093/molbev/mst010
64
KeighleyX.PálssonS.EinarssonB. F.PetersenA.Fernández-CollM.JordanP.et al (2019). Disappearance of Icelandic walruses coincided with Norse settlement.Mol. Biol. Evol.362656–2667. 10.1093/molbev/msz196
65
KircherM.SawyerS.MeyerM. (2012). Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform.Nucleic Acids Res.40:e3. 10.1093/nar/gkr771
66
KnutsenH.OlsenE. M.JordeP. E.EspelandS. H.AndréC.StensethN. C. (2011). Are low but statistically significant levels of genetic differentiation in marine fishes ’biologically meaningful’? A case study of coastal Atlantic cod.Mol Ecol.20768–783. 10.1111/j.1365-294X.2010.04979.x
67
KumarS.StecherG.TamuraK. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets.Mol. Biol. Evol.331870–1874. 10.1093/molbev/msw054
68
LaitL. A. (2016). A mitogenomic study of four at-risk maring fish species: Atlantic wolffish, spotted wolffish, northern wolffish, and Atlantic cod, with special emphasis on the waters off Newfoundland and Labrador.Ph.D. thesis, Canada: Memorial University of Newfoundland.
69
LaitL. A.MarshallH. D.CarrS. M. (2018). Phylogeographic mitogenomics of Atlantic cod Gadus morhua: variation in and among trans−Atlantic, trans−Laurentian, Northern cod, and landlocked fjord populations.Ecol. Evol.86420–6437. 10.1002/ece3.3873
70
LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform.Bioinformatics251754–1760. 10.1093/bioinformatics/btp324
71
LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al (2009). The sequence alignment/map format and SAMtools.Bioinformatics252078–2079. 10.1093/bioinformatics/btp352
72
LimburgK. E.WaltherY.HongB.OlsonC.StoråJ. (2008). Prehistoric versus modern Baltic Sea cod fisheries: selectivity across the millennia.Proc. R. Soc. B Biol. Sci.2752659–2665. 10.1098/rspb.2008.0711
73
LindgreenS. (2012). AdapterRemoval: easy cleaning of next-generation sequencing reads.BMC Res. Notes5:337. 10.1186/1756-0500-5-337
74
LombalA. J.O’dwyerJ. E.FriesenV.WoehlerE. J.BurridgeC. P. (2020). Identifying mechanisms of genetic differentiation among populations in vagile species: historical factors dominate genetic differentiation in seabirds.Biol. Rev.95625–651. 10.1111/brv.12580
75
LouisM.SkovrindM.Samaniego CastruitaJ. A.GarilaoC.KaschnerK.GopalakrishnanS.et al (2020). Influence of past climate change on phylogeography and demographic history of narwhals.Monodon monoceros.Proc. R. Soc. B287:20192964. 10.1098/rspb.2019.2964
76
MacKenzieB. R.AlheitJ.ConleyD. J.HolmP.KinzeC. C. (2002). Ecological hypotheses for a historical reconstruction of upper trophic level biomass in the Baltic Sea and Skagerrak.Can. J. Fish. Aqua. Sci.59173–190. 10.1139/f01-201
77
MakS. S. T.GopalakrishnanS.CarøeC.GengC.LiuS.SindingM. S.et al (2017). Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing.Gigascience61–13. 10.1093/gigascience/gix049
78
MalmstrømM.MatschinerM.TørresenO. K.StarB.SnipenL. G.HansenT. F.et al (2016). Evolution of the immune system influences speciation rates in teleost fishes.Nat. Genet.481204–1210. 10.1038/ng.3645
79
MatschinerM. (2016). Fitchi: haplotype genealogy graphs based on the Fitch algorithm.Bioinformatics321250–1252. 10.1093/bioinformatics/btv717
80
McKennaA.HannaM.BanksE.SivachenkoA.CibulskisK.KernytskyA.et al (2010). The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.Genome Res.201297–1303. 10.1101/gr.107524.110
81
MeyerM.KircherM. (2010). Illumina sequencing library preparation for highly multiplexed target capture and sequencing.Cold Spring Harb. Protoc.2010:db.rot5448.
82
MieszkowskaN.GennerM. J.HawkinsS. J.SimsD. W. (2009). “Chapter 3 Effects of Climate Change and Commercial Fishing on Atlantic Cod Gadus morhua,” inAdvances in Marine Biology, ed.SimsD. W. (Cambridge: Academic Press), 213–273. 10.1016/s0065-2881(09)56003-8
83
MitchellK. J.RawlenceN. J. (2021). Examining Natural History through the Lens of Palaeogenomics.Trends Ecol. Evol.36258–267. 10.1016/j.tree.2020.10.005
84
MossM. L.RodriguesA. T.SpellerC. F.YangD. Y. (2016). The historical ecology of Pacific herring: tracing Alaska Native use of a forage fish.J. Archaeol. Sci. Rep.8504–512. 10.1016/j.jasrep.2015.10.005
85
MurrayG. G.SoaresA. E.NovakB. J.SchaeferN. K.CahillJ. A.BakerA. J.et al (2017). Natural selection shaped the rise and fall of passenger pigeon genomic diversity.Science358951–954. 10.1126/science.aao0960
86
NguyenL.-T.SchmidtH. A.Von HaeselerA.MinhB. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies.Mol. Biol. Evol.32268–274. 10.1093/molbev/msu300
87
NichollsJ.AllaireB.HolmP. (2021). The Capacity Trend Method: a new approach for enumerating the Newfoundland cod fisheries (1675–1790).Histori. Methods J. Quan. Interdiscip. History20211–14. 10.1080/01615440.2020.1853643
88
NikulinaE. A.SchmölckeU. (2016). Reconstruction of the historical distribution of sturgeons (Acipenseridae) in the eastern North Atlantic based on ancient DNA and bone morphology of archaeological remains: implications for conservation and restoration programmes.Divers. Distrib.221036–1044. 10.1111/ddi.12461
89
NyströmV.AngerbjörnA.DalénL. (2006). Genetic consequences of a demographic bottleneck in the Scandinavian arctic fox.Oikos11484–94. 10.1111/j.2006.0030-1299.14701.x
90
OjaveerH.JaanusA.MacKenzieB. R.MartinG.OleninS.RadziejewskaT.et al (2010). Status of Biodiversity in the Baltic Sea.PLoS One5:e12467. 10.1371/journal.pone.0012467
91
OlafsdottirG. A.WestfallK. M.EdvardssonR.PalssonS. (2014). Historical DNA reveals the demographic history of Atlantic cod (Gadus morhua) in medieval and early modern Iceland.Proc. Biol. Sci.281:20132976. 10.1098/rspb.2013.2976
92
OlsonC.WaltherY. (2007). Neolithic cod (Gadus morhua) and herring (Clupea harengus) fisheries in the Baltic Sea, in the light of fine-mesh sieving: a comparative study of subfossil fishbone from the late Stone Age sites at Ajvide, Gotland, Sweden and Jettböle, Åland, Finland.Environ. Archaeol.12175–185. 10.1179/174963107x226435
93
OostingT.StarB.BarrettJ. H.WellenreutherM.RitchieP. A.RawlenceN. J. (2019). Unlocking the potential of ancient fish DNA in the genomic era.Evol. Appl.121513–1522. 10.1111/eva.12811
94
OrtonD. C.MakowieckiD.de RooT.JohnstoneC.HarlandJ.JonssonL.et al (2011). Stable Isotope Evidence for Late Medieval (14th–15th C) Origins of the Eastern Baltic Cod (Gadus morhua) Fishery.PLoS One6:e27568. 10.1371/journal.pone.0027568
95
OrtonD. C.MorrisJ.LockerA.BarrettJ. H. (2014). Fish for the city: meta-analysis of archaeological cod remains and the growth of London’s northern trade.Antiquity88516–530. 10.1017/S0003598X00101152
96
PageA. J.TaylorB.DelaneyA. J.SoaresJ.SeemannT.KeaneJ. A.et al (2016). SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments.Microb. Genom.2:e000056. 10.1099/mgen.0.000056
97
PaijmansJ. L.GilbertM. T. P.HofreiterM. (2013). Mitogenomic analyses from ancient DNA.Mol. Phylogenet. Evol.69404–416. 10.1016/j.ympev.2012.06.002
98
ParadisE. (2010). pegas: an R package for population genetics with an integrated–modular approach.Bioinformatics26419–420. 10.1093/bioinformatics/btp696
99
ParadisE.SchliepK. (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R.Bioinformatics35526–528. 10.1093/bioinformatics/bty633
100
PfeiferB.WittelsbuergerU.LiH.HandsakerB.PfeiferM. B. (2020). Package ‘PopGenome’.URL: https://CRAN.R-project.org/package=PopGenome
101
PinnegarJ. K.EngelhardG. H. (2008). The ‘shifting baseline’phenomenon: a global perspective.Rev. Fish Biol. Fish.181–16. 10.1007/s11160-007-9058-6
102
PinskyM. L.EikesetA. M.HelmersonC.BradburyI. R.BentzenP.MorrisC.et al (2021). Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic.Proc. Natl. Acad. Sci. U. S. A.118:e2025453118. 10.1073/pnas.2025453118
103
RambautA.DrummondA. J.XieD.BaeleG.SuchardM. A. (2018). Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7.Syst. Biol.67901–904. 10.1093/sysbio/syy032
104
RodriguesA. S.MonsarratS.CharpentierA.BrooksT. M.HoffmannM.ReevesR.et al (2019). Unshifting the baseline: a framework for documenting historical population changes and assessing long-term anthropogenic impacts.Philos. Trans. R. Soc. B374:20190220. 10.1098/rstb.2019.0220
105
RozasJ.Ferrer-MataA.Sánchez-DelBarrioJ. C.Guirao-RicoS.LibradoP.Ramos-OnsinsS. E.et al (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets.Mol. Biol. Evol.343299–3302. 10.1093/molbev/msx248
106
SalazarG. (2020). FastaUtils: Utilities for DNA/RNA sequence processing. URL: https://github.com/GuillemSalazar/FastaUtils.git
107
SchroederH.Ávila-ArcosM. C.MalaspinasA.-S.PoznikG. D.Sandoval-VelascoM.CarpenterM. L.et al (2015). Genome-wide ancestry of 17th-century enslaved Africans from the Caribbean.Proc. Natl. Acad. Sci. U. S. A.1123669–3673. 10.1073/pnas.1421784112
108
SchubertM.ErminiL.Der SarkissianC.JónssonH.GinolhacA.SchaeferR.et al (2014). Characterization of ancient and modern genomes by SNP detection and phylogenomic and metagenomic analysis using PALEOMIX.Nat. Protoc.9:1056. 10.1038/nprot.2014.063
109
SeersholmF. V.ColeT. L.GrealyA.RawlenceN. J.GreigK.KnappM.et al (2018). Subsistence practices, past biodiversity, and anthropogenic impacts revealed by New Zealand-wide ancient DNA survey.Proc. Natl. Acad. Sci. U. S. A.115:7771. 10.1073/pnas.1803573115
110
SelimS. A.BlanchardJ. L.BedfordJ.WebbT. J. (2016). Direct and indirect effects of climate and fishing on changes in coastal ecosystem services: a historical perspective from the North Sea.Reg. Environ. Change16341–351. 10.1007/s10113-014-0635-7
111
ShapiroB.DrummondA. J.RambautA.WilsonM. C.MatheusP. E.SherA. V.et al (2004). Rise and fall of the Beringian steppe bison.Science3061561–1565. 10.1126/science.1101074
112
SheltonP. A.MorganM. J. (2014). Impact of maximum sustainable yield-based fisheries management frameworks on rebuilding North Atlantic cod stocks.J. Northwest Atlant. Fish. Sci.4615–25. 10.2960/j.v46.m697
113
SigurgíslasonH.ÁrnasonE. (2003). Extent of mitochondrial DNA sequence variation in Atlantic cod from the Faroe Islands: a resolution of gene genealogy.Heredity91557–564. 10.1038/sj.hdy.6800361
114
SodelandM.JordeP. E.LienS.JentoftS.BergP. R.GroveH.et al (2016). “Islands of Divergence” in the Atlantic Cod Genome Represent Polymorphic Chromosomal Rearrangements.Genome Biol. Evol.81012–1022. 10.1093/gbe/evw057
115
SpellerC. F.HauserL.LepofskyD.MooreJ.RodriguesA. T.MossM. L.et al (2012). High Potential for Using DNA from Ancient Herring Bones to Inform Modern Fisheries Management and Conservation.PLoS One7:e51122. 10.1371/journal.pone.0051122
116
SpencerH. G. (2020). Beyond Equilibria: the Neglected Role of History in Ecology and Evolution.Quart. Rev. Biol.95311–321. 10.1086/711803
117
StarB.BarrettJ. H.GondekA. T.BoessenkoolS. (2018). Ancient DNA reveals the chronology of walrus ivory trade from Norse Greenland.Proc. R. Soc. B Biol. Sci.285:20180978. 10.1098/rspb.2018.0978
118
StarB.BoessenkoolS.GondekA. T.NikulinaE. A.HufthammerA. K.PampoulieC.et al (2017). Ancient DNA reveals the Arctic origin of Viking Age cod from Haithabu. Germany.Proc. Natl. Acad. Sci. U. S. A.1149152–9157. 10.1073/pnas.1710186114
119
StarB.NederbragtA. J.JentoftS.GrimholtU.MalmstrømM.GregersT. F.et al (2011). The genome sequence of Atlantic cod reveals a unique immune system.Nature477207–210.
120
StillerM.BaryshnikovG.BocherensH.Grandal, d’AngladeA.HilpertB.et al (2010). Withering Away—25,000 Years of Genetic Decline Preceded Cave Bear Extinction.Mol. Biol. Evol.27975–978. 10.1093/molbev/msq083
121
SætersdalG.HylenA. (1964). The decline of the skrei fisheries: a review of the landing statistics 1866-1957 and an evaluation of the effects of the postwarincrease in the total exploitation of the arctic cod.Fiskeridirektoratets havforskningsinstitutt1356–69.
122
TherkildsenN. O.NielsenE. E.SwainD. P.PedersenJ. S. (2010). Large effective population size and temporal genetic stability in Atlantic cod (Gadus morhua) in the southern Gulf of St. Lawrence.Can. J. Fish. Aquat. Sci.671585–1595. 10.1139/f10-084
123
ThomasJ. E.CarvalhoG. R.HaileJ.RawlenceN. J.MartinM. D.HoS. Y. W.et al (2019). Demographic reconstruction from ancient DNA supports rapid extinction of the great auk.Elife8:e47509. 10.7554/eLife.47509
124
ThurstanR. H.BrockingtonS.RobertsC. M. (2010). The effects of 118 years of industrial fishing on UK bottom trawl fisheries.Nat. Commun.1:15. 10.1038/ncomms1013
125
TørresenO. K.StarB.JentoftS.ReinarW. B.GroveH.MillerJ. R.et al (2017). An improved genome assembly uncovers prolific tandem repeats in Atlantic cod.BMC Genomics18:95. 10.1186/s12864-016-3448-x
126
VenterO.SandersonE. W.MagrachA.AllanJ. R.BeherJ.JonesK. R.et al (2016). Sixteen years of change in the global terrestrial human footprint and implications for biodiversity conservation.Nat. Commun.7:12558. 10.1038/ncomms12558
127
WelchA. J.WileyA. E.JamesH. F.OstromP. H.StaffordJr.FleischerR. C. (2012). Ancient DNA reveals genetic stability despite demographic decline: 3,000 years of population history in the endemic Hawaiian petrel.Mol. Biol. Evol.293729–3740. 10.1093/molbev/mss185
128
WenneR.BernaśR.KijewskaA.Poćwierz-KotusA.StrandJ.PetereitC.et al (2020). SNP genotyping reveals substructuring in weakly differentiated populations of Atlantic cod (Gadus morhua) from diverse environments in the Baltic Sea.Sci. Rep.101–15.
129
ZhangJ.PeiN.MiX.ZhangM. J. (2017). Package ‘phylotools’. dimension12. URL: https://github.com/helixcn/phylotools
Summary
Keywords
population structure, fisheries, human exploitation, phylogenomics, population expansion, demographic history
Citation
Martínez-García L, Ferrari G, Oosting T, Ballantyne R, van der Jagt I, Ystgaard I, Harland J, Nicholson R, Hamilton-Dyer S, Baalsrud HT, Brieuc MSO, Atmore LM, Burns F, Schmölcke U, Jakobsen KS, Jentoft S, Orton D, Hufthammer AK, Barrett JH and Star B (2021) Historical Demographic Processes Dominate Genetic Variation in Ancient Atlantic Cod Mitogenomes. Front. Ecol. Evol. 9:671281. doi: 10.3389/fevo.2021.671281
Received
26 February 2021
Accepted
06 May 2021
Published
04 June 2021
Volume
9 - 2021
Edited by
Melissa T. R. Hawkins, Smithsonian Institution, United States
Reviewed by
Martin Taylor, University of East Anglia, United Kingdom; Athanasios Exadactylos, University of Thessaly, Greece
Updates
Copyright
© 2021 Martínez-García, Ferrari, Oosting, Ballantyne, van der Jagt, Ystgaard, Harland, Nicholson, Hamilton-Dyer, Baalsrud, Brieuc, Atmore, Burns, Schmölcke, Jakobsen, Jentoft, Orton, Hufthammer, Barrett and Star.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lourdes Martínez-García, l.m.garcia@ibv.uio.noBastiaan Star, bastiaan.star@ibv.uio.no
†ORCID: Lourdes Martínez-García, orcid.org/0000-0002-1582-3611; Giada Ferrari, orcid.org/0000-0002-0850-1518; Tom Oosting, orcid.org/0000-0002-7031-0747; Rachel Ballantyne, orcid.org/0000-0002-6506-3163; Helle Tessand Baalsrud, orcid.org/0000-0002-4161-3247; Marine Servane Ono Brieuc, orcid.org/0000-0001-8601-2122; Lane M. Atmore, orcid.org/0000-0002-8903-8149; Kjetill S. Jakobsen, orcid.org/0000-0002-8861-5397; Sissel Jentoft, orcid.org/0000-0001-8707-531X; David Orton, orcid.org/0000-0003-4069-8004; James H. Barrett, orcid.org/0000-0002-6683-9891; Bastiaan Star, orcid.org/0000-0003-0235-9810
This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Ecology and Evolution
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.