- Research article
- Open Access
The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences
BMC Evolutionary Biology volume 18, Article number: 87 (2018)
Commercial whaling caused extensive demographic declines in many great whale species, including gray whales that were extirpated from the Atlantic Ocean and dramatically reduced in the Pacific Ocean. The Eastern Pacific gray whale has recovered since the 1982 ban on commercial whaling, but the Western Pacific gray whale—once considered possibly extinct—consists of only about 200 individuals and is considered critically endangered by some international authorities. Herein, we use whole-genome sequencing to investigate the demographic history of gray whales from the Pacific and use environmental niche modelling to make predictions about future gene flow.
Our sequencing efforts and habitat niche modelling indicate that: i) western gray whale effective population sizes have declined since the last glacial maximum; ii) contemporary gray whale genomes, both eastern and western, harbor less autosomal nucleotide diversity than most other marine mammals and megafauna; iii) the extent of inbreeding, as measured by autozygosity, is greater in the Western Pacific than in the Eastern Pacific populations; and iv) future climate change is expected to open new migratory routes for gray whales.
Our results indicate that gray whale genomes contain low nucleotide diversity and have been subject to both historical and recent inbreeding. Population sizes over the last million years likely peaked about 25,000 years before present and have declined since then. Our niche modelling suggests that novel migratory routes may develop within the next century and if so this could help retain overall genetic diversity, which is essential for adaption and successful recovery in light of global environmental change and past exploitation.
Widespread commercial whaling during the last two centuries unsustainably harvested many whale populations . Whale products such as oil, meat, blubber, and ambergris were commercially important and overharvesting greatly diminished many whale populations [2,3,4,5]. In 1982 the International Whaling Commission (IWC) instituted a moratorium on commercial whaling (https://iwc.int/commercial), and although some whale populations have since recovered to near their pre-whaling abundance, others remain compromised. Recent, anthropogenic bottlenecks due to commercial whaling can be contrasted with more ancient, natural bottlenecks often associated with climate and/or ecological change .
Great whales are important for marine ecosystems, as they facilitate nutrient transfer in the water column and stabilize ecosystems by increasing biodiversity . Whales are associated with areas of high primary productivity, and their sensitivity to environmental changes make them prime indicators of ecological perturbations . The marine ecosystem is rapidly changing due to anthropogenic impacts [8,9,10], most of which have unknown consequences for the future of marine environments and marine mammals [11, 12]. Scientists are just beginning to understand how large marine mammals have responded to past climatic cycles [13,14,15], and models predict that range and distribution patterns will shift towards the poles in the face of global warming .
One species severely affected by commercial whaling is the gray whale (Eschrichtius robustus). Gray whales were once common in the Northern hemisphere, but were extirpated from the Atlantic ocean by the early eighteenth century , potentially due to environmental change and/or by commercial whaling [2, 18]. Today, gray whales are found in the Eastern Pacific near the coast of North America and the Western Pacific near the coast of Asia (Fig. 1). There is evidence of gene flow between the two “stocks”, but there is also statistically detectable genetic differentiation between them [2, 19]. The eastern gray whale (EGW) population has been extensively studied, and post-whaling estimates based on genetic and ecological data indicate there are ~ 27,000 individuals [19,20,21]. In contrast, data on the western population is limited [22, 23]. Commercial whaling lasted considerably longer in the western Pacific , and today the western gray whale (WGW) is thought to be comprised of < 200 individuals and is listed as ‘critically endangered’ by the International Union for Conservation of Nature (IUCN) [25, 26].
During the late Pleistocene and Holocene (i.e., within the last ~ 100,000 years), the Northern hemisphere experienced massive changes to its marine ecosystems . Glacial periods led to ice cap oscillations that repeatedly opened and closed migration corridors [28,29,30], and fluctuations in water temperature and sea levels likely forced changes to habitats and feeding modes [2, 3]. Gray whale carrying capacities have been modelled based on shifts in feeding habitat during the last 120,000 years, and those data suggest that multiple demographic bottlenecks may have occurred . In addition to the trophic data, DNA sequences suggest that the EGW population has been subject to a genetic bottleneck during the last century . Although population fluctuations have not been investigated in the WGW, microsatellite and mitochondrial data suggest that the two populations have similar levels of neutral genetic diversity and thus may have similar long-term demographic histories [2, 19].
The ongoing reductions in the extent of sea ice provide gray whales with new potential migration routes and they may be shifting their range farther north in the Arctic . Gray whales have responded to climate changes by shifting the timing of their southbound migration . Because they annually migrate thousands of kilometres from their summer feeding grounds at high latitudes to their winter calving waters at lower latitudes, there may be opportunities for contemporary (i.e., within the last few dozen generations) gene flow between eastern and western Pacific populations (Fig. 1). There are indications of historical (referring to the Pleistocene and early Holocene) gene flow between Atlantic and Pacific gray whales , and—more recently—a satellite-tagged WGW has been tracked to an EGW wintering area near the Mexican coast . Furthermore, photographic identification has documented individual gray whales moving between the western Pacific (near Sakhalin Island, Russia) and the eastern Pacific . Collectively, these data suggest that the currently recognized WGW and EGW “populations” of this highly vagile species are not completely independent (i.e., gene flow is possible). Fortunately, population attributes such as historical demography, admixture (i.e., interbreeding between populations that have previously been isolated), and genetic diversity can now be addressed using whole genome sequences [33,34,35]. The ability to make population inferences from one or a few samples is especially important for rare species, where sampling efforts are often difficult, expensive, and should be minimized because of conservation concerns.
Herein, we employ genomic and computational techniques to infer population attributes of gray whales. The distribution of gray whales is largely disjunct today, but these geographic isolates were demographically and genetically connected in the past (as evidenced by the fact that they are recognized as a single species). Given the recent growth of the EGW population and ongoing climate change, there is reason to suspect that increased gene flow between EGW and WGW may occur in the future. We are interested in the long-term demographic trajectory of gray whales, both from a historical and a future perspective. Given the critically endangered status of the WGW, we were interested in comparing genomes of the WGW and the EGW to investigate levels of genetic diversity as a key component of adaptive potential. We used coalescent-based approaches to retrospectively gauge ancient admixture in gray whale genomes during the Pleistocene, and measures of autozygosity to directly assess inbreeding and search for signals of contemporary differentiation. Our habitat prediction models suggest that novel migratory routes may develop within the next century, which could influence the overall retention of genetic diversity in the species. This study presents the first genomic comparison of gray whales, and extends our insights into the molecular diversity and demographic history of this enigmatic species while contributing to our understanding of how our ocean’s great whales have responded to historical climate change.
Sampling, sequencing and SNP calling
We used previously published whole genome DNA sequences from DeWoody et al. . These sequences were derived from two gray whales sampled near Sakhalin Island, Russia, designated WGW1 (female) and WGW2 (male), and from one putative Eastern gray whale female (EGW) that was beached near Barrow, Alaska (Fig. 1). There is some uncertainty as to true population affinities of these individual gray whales. For example, WGW1 was biopsied near Sakhalin Island in the western Pacific but the same whale has been photographically identified in the eastern Pacific (Laguna San Ignacio; M. Scott, unpublished data). Nevertheless, we assigned geographical names to whales based on sampling locations in order to be comparable with previous genetic work on gray whales [2, 19, 37]. The data utilized herein consisted of 2x100bp paired-end (PE) libraries from each whale (~ 1 billion reads per individual; ~ 700 million high-quality reads per individual after quality-control; Additional file 1: Table S1). For detailed sampling and genome sequencing methodology see DeWoody et al. . For a summary of number of reads per individual and quality control effects, see Additional file 1: Table S1.
We used FASTQC v0.11.2 (www.bioinformatics.babraham.ac.uk/projects/fastqc) to generate summary statistics for the sequencing reads. TRIMMOMATIC v0.32  was used to remove adaptor sequences and trim low quality bases (< 20 Phred scores) from both the 5′ and 3′ end of each read. BWA v0.7.12  was used to map the PE reads to the published genome of the common minke whale (B. acutorostrata) (GenBank accession: SAMN02192642, ) using the ‘bwtsw’ function that indexes whole genomes, and the ‘mem’ function for mapping. PICARD-TOOLS v2.0.1 (http://broadinstitute.github.io/picard/.) was employed to mark duplicate reads. SAMTOOLS v1.3  was used for alignment manipulation. Local realignment, duplicate removal, and SNP variant calling were carried out with GATK v3.5  following ‘Best Practices protocol’ [43, 44]. Genotypes were called across all three samples together using the ‘gvcf’ option. We used a minimum base quality score of 20 (which corresponds to a base calling error rate of ~ 1% ) with a minimum mapping quality score of 20. In the downstream analyses, we only used SNPs with minimum 20× coverage, which should help minimize the number of heterozygotes falsely scored as homozygotes [46, 47]. Eight minke whale scaffolds are X-linked , and we removed gray whale reads that mapped to these scaffolds so they would not bias our downstream analyses. None of our gray whale reads mapped to Y-linked scaffolds because  reported none in their minke genome assembly.
Nucleotide diversity can be used to assess ancient admixture as well as contemporary differentiation [34, 49, 50]. We estimated observed heterozygosity for each individual, θgenome, based on the number of heterozygous sites / total number of sites where only sites with minimum 20× coverage were considered. We used θ values associated with each individual to independently estimate equilibrium effective population sizes (N e ) following θ = 4Neμ . To quantify differences in N e we compared θ among individuals, assuming that substitution rates do not vary appreciably across samples.
We directly quantified inbreeding levels by identifying the number and lengths of autosomal runs-of-homozygosity (ROHs) in each individual. A ROH is a genomic region that contains far less nucleotide variation than expected based on the genome–wide average for an individual . Under random mating, the length of ROH regions is expected to decrease with increasing number of generations to the ‘most recent common ancestor’ (MRCA) due to recombination and de novo mutations. In contrast, with inbreeding—as is often the case for critically endangered species—autozygosity is expected to increase over time, thus increasing the number and length of ROHs in the genome each generation. Analysis of ROH abundance and extent thus provides information on a population’s demographic history and on the genetic relationships among individuals .
We estimated four different ROH parameters: i) number of ROHs in each genome (NROH); ii) the mean length of ROHs (LROH); iii) the heterozygosity outside ROHs (θnoROH); and iv) the inbreeding coefficient FROH, the overall proportion of the genome contained in ROHs. We estimated θnoROH as the number of heterozygous SNPs / (total number of SNPs – SNPs in ROHs). When FROH is compared to θgenome, it quantifies the effect of inbreeding on overall levels of genomic variation. To compare our results directly to patterns of ROHs found in other species , we used PLINK v1.90b3.36  and defined ROHs as portions of the genome that spanned at least 20 homozygous sites allowing for a single heterozygous SNP (e.g., due to de novo mutation) and 1 missing SNP (e.g., a site with missing data) following Howrigan et al. . We searched the genomes for ROHs in consecutive 20 SNP sliding windows and, to facilitate detecting both short and long ROHs, we set the lower bound for ROHs to 1 kb. We used a Welch two-sample t-test to test for pairwise differences in LROH among individuals whereas pairwise ROH frequency distributions were compared among all three gray whales using a two-sample Kolmogorov–Smirnov test. All statistical tests were conducted in R .
Relatedness and population structure
We used PLINK to measure relatedness among individuals. Pairwise identical-by-state (IBS) comparisons were estimated based on the ratio of probabilities between a heterozygote–heterozygote site, p(HetHet) = 4p2q2, to the probability of a homozygote–homozygote site, p(HomHom) = p2q2. For each pair of individuals, the number of variable sites where they share no alleles (IBS = 0; e.g., discordant homozygotes AA/BB and BB/AA) are counted along with the number of sites where they share two alleles (IBS = 2) (e.g., heterozygotes AB/AB, BA/BA). On average, we expect this probability ratio to be 1:2 if the pair comes from a randomly mating population [55, 58]. A ‘HetHet’: HomHom’ ratio > 2 suggests that the individuals are more related than expected by chance, and a HetHet’: HomHom’ ratio < 2 suggests that the individuals have recent ancestry from different random mating populations. We used the ‘pairwise population concordance’ (PPC) test to evaluate if this probability ratio significantly deviated from the expected ratio under random mating, applying a significance level of 0.05, a minor allele frequency (MAF) of 0.01, and a minimum distance of 500 k base pairs between informative SNPs to limit the effects of linkage disequilibrium (LD).
To test for ancient admixture, we used the ABBA–BABA D-statistic test implemented in ANGSD v0.912 [34, 59]. The D-statistic tests for admixture between four individuals: two conspecific individuals (P1 and P2), a potential introgressor (P3), and an outgroup (O). At each polymorphic site in the genome the relationship among these four individuals and the topology of the species tree is compared. Sites that are inconsistent with the species tree are the sites where P2 shares a derived allele with P3 but not P1 (ABBA sites) or P1 shares derived sites with P3 but not P2 (BABA sites). An excess of either ABBA or BABA sites, compared to the sites supporting the species tree (i.e., AABB), is an indication of admixture between P2 and P3 or between P1 and P3, respectively. In the absence of ancient population structure, incomplete lineage sorting is the only process other than admixture that produces inconsistency with the species tree topology, but incomplete lineage sorting is expected to produce ABBA and BABA sites in an equal ratio [49, 50]. The D–test statistic evaluates the number (n) of ABBA and BABA sites (D = (nABBA - nBABA) / (nABBA + nBABA)) and D < 0 means that P1 is more closely related to P3 than to P2, whereas D > 0 indicates that P2 is more closely related to P3 than P1. The significance of the D test was evaluated with a Z-score, where |Z-scores| > 3 was used as the critical value for a significant test . As an outgroup, we used the common minke whale. The phylogenetic relationships among baleen whales are not completely resolved [60, 61], but the common minke whale is the closest relative with a published genome sequence . We used an LD block size of 10 Mb; increasing the block size (e.g., 20 Mb, 30 Mb) did not change the outcome of the ABBA-BABA test. We tested all scaffolds > 10 Mb in length in order to obtain a reliable Z-score. Admixture D-statistics were considered significant for |Z-scores| > 3.
Inference of demographic history
We used the PSMC’ mode implemented in MSMC [33, 63] to infer ancient demographic histories. Eleven scaffolds larger than 30 Mb in length, corresponding to a total of ~ 400 Mb, were used to improve the accuracy of inferring past recombination events [33, 64]. We ran the MSMC analysis for each individual separately using default settings; 20 iterations and averaging over 30 time segments. To quantify the variance in Ne we bootstrapped using the same MSMC settings. For each individual, 20 bootstrapped datasets were generated by randomly sampling 5 Mb sequences from each of the 11 scaffolds used to trace the mean Ne. Substitution rates—for both mitochondria and nuclear loci—are reportedly 8–10 fold slower in baleen whales than in other mammals [65, 66]. In order to convert θ to Ne over time, we applied an autosomal substitution rate of 4.8 × 10− 10 bp− 1 year− 1 (credibility interval (CI): 1.5 × 10− 10 – 10 × 10− 10) , and a generation time of 18.9 years which corresponds to the midpoint of estimated generation times which range between 15.5 and 22.3 years [68, 69]. MSMC runs were assessed for convergence using the R package CODA .
Prediction of suitable habitat
We used AQUAMAPS  to predict the relative probability of the future gray whale distribution across the Northern Hemisphere based on contemporary local conditions. Suitable habitat was based on occurrence records available via Ocean Biogeographic Information System (http://www.iobis.org) using the contemporary environmental envelope settings suggested by Alter et al.  (Additional file 1: Table S2), and future (year 2100) envelope settings from AQUAMAPS . We assumed that current environmental conditions are representative of the Holocene, as the Holocene climate has experienced relatively little variation compared to interglacial cycles .
We mapped, from each individual, high-quality PE reads from one eastern and two western gray whales to the minke whale genome. The mean depth of coverage per individual ranged from 27× to 30× (Additional file 1: Table S1), and this relatively deep coverage allowed us to assess nucleotide diversity with confidence. The level of genetic diversity represented by theta (θ) was lower in the individuals from the Western population (θ = 6.69 × 10− 4 and 6.64 × 10− 4) relative to the individual from the putative Eastern population θ = 8.00 × 10− 4 (Fig. 2). Thus, there is about a 1.2-fold difference in genetic diversity between East and West.
We found ROHs ranging from 1 to 559 Kb in length; few were longer than 300Kb (Fig. 3). Estimates of θgenome and θnoROH were lower in both WGWs than in the EGW (Table 1). The western individuals had fewer ROHs (WGW1: nROH = 188,012 and WGW2: nROH = 126,893) than the Eastern individual (nROH = 263,877), but their mean ROH length were significantly longer (WGW1: LROH = 11Kb and WGW2: LROH = 17Kb; both p = 2.2 × 10− 16) than in the eastern individual (EGW: LROH = 6Kb). ROHs covered a larger proportion of the western gray whale genomes (WGW1: TROH = 2.1 × 106 bp; WGW2: TROH = 2.2 × 106 bp) compared to the eastern gray whale (EGW: TROH = 1.6 × 106 bp) (Table 1). All individuals differed significantly from one another in LROH (p = 2.2 × 10− 16), and the Kolmogorov–Smirnov test showed that the ROH distributions were significantly different from one another (WGW1 & EGW D = 0.186, p = 2.2 × 10− 16; WGW2 & EGW D = 0.322, p = 2.2 × 10− 16, and WGW1 & WGW2 D = 0.142 p = 2.2 × 10− 16). This suggest that there are significant differences in genealogical histories between all individuals. Estimates of FROH were 0.088 (WGW1), 0.092 (WGW2), and 0.067 (EGW) and thus on average the WGWs were ~ 1.3 times as inbred as EGW.
Relatedness and population structure
To evaluate pairwise relatedness, we used the ‘HetHet’ to ‘HomHom’ ratios (where a ratio of 2.0 is expected for individuals from the same random mating population and a ratio > 2.0 suggests that the pair is more related to each other than expected based on chance alone). All pairwise comparisons yielded a ‘HetHet’ to ‘HomHom’ ratio ≥ 2 (Table 2), and the PPC test could not reject the null hypothesis: ‘HetHet’: HomHom’ ratio = 2 (Table 2). Thus this test is uninformative as the three individuals may or may not belong to the same gene pool.
The ABBA–BABA test revealed no significant support for ancient admixture (e.g., historical panmixia) between the Western and Eastern gray whales (Fig. 4; Additional file 1: Table S3). If ABBA and BABA patterns are equally common, then in theory D = 0 and the data are consistent with the tree. Deviations where D ≠ 0 can be due to: i) P3 exchanged genes with P1 or P2; ii) ancestral population (P1, P2 and P3’s founder) structure leading to discordant gene trees; or iii) P1 or P2 could have received genes from an unsampled ‘ghost’ population (Pg). The test is not influenced by demographic events assuming that P1, P2 and P3’s ancestral population was panmictic , which should be a reasonable assumption for gray whales .
Inference of long-term demographic history
We traced effective population size estimates over the last ~ 1,000,000 years using the PSMC’ method (Fig. 5; Additional file 1: Figure S1). The three individuals exhibit very consistent trajectories, indicating a step decline in Ne from Ne > 50,000 in the interval of ~ 1,000,000 years before present (YBP) until 100,000 YBP followed by a more stable period (~ 100,000–30,000 YBP) with Ne ~ 25,000 for both EGW and WGW populations. Prior to the LGM both populations increase in size to Ne ~ 45,000; hereafter a reduction in Ne to a population size of Ne~ 20,000 is observed in all three trajectories. These consistent results among individuals suggest there is relatively little noise in this PSMC’ analysis and that the trajectories themselves are likely a realistic representation of historical population dynamics. Furthermore, the most recent estimate of census population size (Nc) of the EGW is 27,000 . The concordance between Nc and recent Ne estimates (Fig. 5) suggests that the substitution rate we used (4.8 × 10− 10 bp− 1 year− 1) is a reasonable approximation of the true genome-wide substitution rate.
Predictions of suitable habitat
Our environmental niche modelling suggests that current habitat suitability is relatively high from Taiwan to Kamchatka through much of the Bering Sea and along the coast of North America to the Gulf of California (Fig. 1a). Currently marginal habitat, which is expected to improve in the future due to ongoing climate change, includes the Arctic and Chukchi Seas (Fig. 1b).
Anthropogenic factors are rapidly changing the global environment. We think that predictions regarding future biological impacts (e.g., species range shifts) are most informative when presented in a historical context. Genomic data have great potential in this regard as they can be used as a window to the past (e.g., the reconstruction of past demographic histories) and into the future (e.g., by identifying genes expected to face particular selection pressures, such as those related to thermoregulation). Using whole genome data from contemporary eastern and western gray whale populations, we quantified genetic diversity in gray whales and inferred key population attributes that bear on their evolution and conservation.
Genetic diversity and inbreeding
The genome-wide heterozygosity in gray whales is similar to the minke whale, but lower than other marine mammals—e.g., sperm whales (Physeter catodon), common bottlenose dolphin (Tursiops truncatus), killer whales (Orcinus orca), and manatees (Trichechus manatus latirostris)—and considerably lower than terrestrial megafauna (i.e., African elephant (Loxodonta africana), camels (Camelus bactrianus and C. ferus), white rhinoceros (Ceratotherium simum)) (Fig. 2). We expect that the variation in θ may be explained in part by differences in body size; larger animals have slower mutation rates, longer generation times, and produce fewer offspring—all factors that impact θ [74, 75]. Gray whales are the largest of the mammals surveyed here, which could partly explain their low genomic diversity, but population declines over the last ~ 20,000 years (Fig. 5) may also be a significant contributing factor.
Reduced genomic diversity is a concern as it constrains adaptive potential . We observed lower θgenome and θnoROH in western than eastern gray whales (Table 1), likely due to the smaller size of the western population compared to the eastern population. Small population sizes and reduced gene flow will lead to increased inbreeding that has the potential to reduce reproductive fitness due to homozygosity of deleterious recessive alleles and to reduced heterosis. The extent of ROHs in a genome is correlated with population size reductions and increased consanguinity [52, 53]. Our data indicate that, consistent with contemporary population sizes, ROHs significantly reduce overall nucleotide variation in the gray whale genome (Table 1). The timing and duration of bottlenecks are directly associated with the extent of ROHs; i.e., recent inbreeding leads to long ROHs whereas ancient inbreeding persists in the genome as shorter ROHs that have been disrupted by mutation and recombination [77, 78]. The eastern gray whale had more but shorter ROHs than the western gray whales (Table 1, Fig. 3). This is not surprising given that the eastern population is ~ 100× larger and has not experienced extensive recent inbreeding . In contrast, the western gray whale individuals had fewer but longer ROHs and a larger proportion of their genomes in ROHs (Table 1), a pattern that can be produced by a continuous small population size or a genetic bottleneck that persists for multiple generations . The small size of the western population (< 200 individuals) may not only have led to loss of genetic diversity, but also the loss of adaptive potential in the face of impending environmental change [8,9,10].
Relatedness, gene flow and geographical isolates
Gray whales are one of the most vagile species on earth; telemetry and photographic data indicate that some individuals annually move thousands of kilometres across the Pacific [22, 26]. This contemporary movement of individuals between eastern and western populations provides opportunities for gene flow. Furthermore, our niche modelling suggests that gray whales from the east and from the west could encounter the same suitable habitat (Fig. 1b). However, despite the potential overlap in suitable habitat and the known movement of individuals between the populations, their genomes significantly differ in terms of homozygosity (Fig. 3). Thus the ROH data are consistent with previous reports of population structure between eastern and western gray whales [2, 19]. However the PPC test could not reject the null hypothesis of random mating (Table 2; p = 1.00) and the relatedness analysis showed that the EGW was more closely related to both of the WGWs than expected by chance. These PPC and relatedness results are consistent with an earlier relatedness analysis based on 88 gene-associated SNPs, which found the EGW was no more or less related to the WGW population than expected on the basis of chance alone .
During the Pleistocene, climate-dependent dispersal occurred between the Pacific and Atlantic gray whale populations; prior to the last glacial period (110,000–11,700 YBP) and after the opening of the Bering Strait, gray whales migrated between the Pacific and Atlantic oceans . Analyses of mitochondrial sequences have documented haplotype sharing between the eastern and western populations, suggesting that recent maternal gene flow has occurred during the Holocene . In 2010, a Pacific gray whale was observed in the Mediterranean Sea, a sighting which produced speculation that climate-induced shrinking of the Arctic Sea ice may ultimately enable gray whales to recolonize the Atlantic . Thus, although gene flow between the western and eastern populations has no doubt occurred multiple times since the Pleistocene, the signal of contemporary genomic structure we detected between geographic populations is stronger than the signal of historical admixture (Fig. 4). For all D–test statistics we found the two WGWs to be more closely related to each other than either was to the EGW (Additional file 1: Table S3), although we could not reject the hypothesis that the individuals belong to the same randomly mating population (Table 2).
Dating and severity of population decline (s)
The dating of population size changes is inexact due to errors in the estimation of mutation rates and generation times, but genetic datasets are nevertheless often highly concordant with independent datasets (e.g., fossil evidence; [80, 81]). Demographic histories inferred from single whole genome sequences trace from the two haplotypes to their coalescence in the MRCA. This means that the most recent past is not well resolved, and if unphased haplotypes are used—as done in this study—this also affects deep (past) resolution . Thus, our PSMC’ analyses are unlikely to recover any Anthropocene population size changes associated with commercial whaling, as any genetic signal this may have left is much too recent for this method to detect. That said, trajectories of Ne over the last ~ 1,000,000 years are highly consistent with one another and suggest a similar demographic history in each lineage (Fig. 5). Pre-whaling eastern gray whale census population size (Nc) has been estimated at 96,000 (CI: 76,000–118,000) individuals based on nuclear microsatellites , whereas mitochondrial DNA sequences  yield Nc estimates of 100,670 (90% HPD: 59,940–111,550). These Nc estimates correspond to Ne of ~ 32,000 (CI: 25,000–39,000) for microsatellites, and Ne ~ 17,000 (90% HPD: 10,000–19,000) for mitochondrial DNA, which is similar to our post LGM Ne estimate ~ 20,000 (Fig. 5). These differences among studies may illustrate that using a subset of genomic markers does not accurately capture overall genomic diversity perhaps because of the ascertainment bias associated with the selection of highly polymorphic markers such as microsatellites .
Our data suggest an ancient population decline during previous ice ages and a more recent decline in the last ~ 25,000 years (Fig. 5). Glacial periods are often associated with population declines, and the large shifts in climate have impacted both terrestrial [83, 84] and marine mammals [85,86,87,88]. Taken together with the evidence for contemporary bottlenecks—occurring around the time of commercial whaling —these results support population models which indicate multiple bottlenecks have occurred in gray whales . Cumulatively, these bottlenecks may have contributed to the relative paucity of genetic diversity observed in gray whales (Fig. 2).
Western gray whales (WGWs)
We were particularly interested in tracing the demographic history and quantifying genetic diversity within the WGW because of its conservation status. We found that WGWs had increased autozygosity (higher FROH) and lower θgenome (Table 2) compared to the eastern gray whale, both of which would be expected in a small inbred population . However, despite having a more than 100–fold difference in census population size, the genomic differences were modest as Ne only differed 1.2 fold between the two geographic populations. The observed ROH patterns suggest that the western population has experienced population size reduction and an elevated level of inbreeding relative to the eastern individual (Fig. 3). These ROH patterns likely result from recent processes (e.g., inbreeding and drift) as opposed to a long-term small population size, which should be reflected in the Pleistocene Ne (Fig. 5). The small population size and low genetic diversity limit the potential evolutionary responses to future environmental change, and thus ongoing efforts to conserve the WGW are critical. Our samples sizes are large in terms of number of genetic loci, but small in terms of individual animals. Future studies will reveal whether the patterns we observe herein are indicative of the species as a whole.
Whole genome sequencing of cetaceans provides new insights into how these enigmatic animals have responded to past and ongoing changes in the marine environment. Herein, we present the first genome-scale study of gray whale demographic history. Our results show that gray whales from the eastern and western Pacific have low genetic diversity, that the past gray whale population (s) was much larger and experienced multiple declines since the Pleistocene, and that there is some evidence of geographic structuring between the populations. Ecological predictions for the year 2100 suggest the current habitat of gray whales in the Pacific Ocean is unlikely to decrease while their former habitat in the Atlantic Ocean could expand with global warming . Combined with decreasing sea ice cover in the Arctic, this expanding habitat could provide gray whales with opportunities to use alternative migration routes that could genetically bind east and west  but only time will tell how anthropogenic effects, genetic drift, inbreeding, and climate change will impact the population viability of gray whales over the long-term.
eastern gray whale
- F ROH :
International Union for Conservation of Nature
International Whaling Commission
Last Glacial Maximum
mean length of ROHs (LROH) minor allele frequency
most recent common ancestor
- N e :
effective population sizes
- N ROH :
number of ROHs
- P1; P2:
Pairwise population concordance test
western gray whale
years before present
- θ :
- θ genome :
- θ noROH :
heterozygosity outside ROHs
Reeves RR. Dolphins, whales and porpoises: 2002–2010 conservation action plan for the world's cetaceans, vol. 58. Gland: IUCN; 2003.
Alter SE, Meyer M, Post K, Czechowski P, Gravlund P, Gaines C, Rosenbaum HC, Kaschner K, Turvey ST, van der Plicht J, et al. Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100. Mol Ecol. 2015;24(7):1510–22.
Pyenson ND, Lindberg DR. What happened to gray whales during the Pleistocene? The ecological impact of sea-level change on benthic feeding areas in the North Pacific Ocean. PLoS One. 2011;6(7):e21295.
Fujiwara M, Caswell H. Demography of the endangered North Atlantic right whale. Nature. 2001;414:537.
Whitehead H, Christal J, Dufault S. Past and distant whaling and the rapid decline of sperm whales off the Galápagos Islands. Conserv Biol. 1997;11(6):1387–96.
Roman J, Estes JA, Morissette L, Smith C, Costa D, McCarthy J, Nation JB, Nicol S, Pershing A, Smetacek V. Whales as marine ecosystem engineers. Front Ecol Environ. 2014;12(7):377–85.
Moore SE. Marine mammals as ecosystem sentinels. J Mammal. 2008;89(3):534–40.
Lewison RL, Crowder LB, Wallace BP, Moore JE, Cox T, Zydelis R, McDonald S, DiMatteo A, Dunn DC, Kot CY. Global patterns of marine mammal, seabird, and sea turtle bycatch reveal taxa-specific and cumulative megafauna hotspots. PNAS. 2014;111(14):5271–6.
Schipper J, Chanson JS, Chiozza F, Cox NA, Hoffmann M, Katariya V, Lamoreux J, Rodrigues ASL, Stuart SN, Temple HJ, et al. The status of the World's land and marine mammals: diversity, threat, and knowledge. Science. 2008;322(5899):225–30.
Abram NJ, McGregor HV, Tierney JE, Evans MN, McKay NP, Kaufman DS, Thirumalai K. Consortium Pk: early onset of industrial-era warming across the oceans and continents. Nature. 2016;536(7617):411–8.
Jackson JBC. Ecological extinction and evolution in the brave new ocean. PNAS. 2008;105(Supplement 1):11458–65.
Doney SC, Ruckelshaus M, Duffy JE, Barry JP, Chan F, English CA, Galindo HM, Grebmeier JM, Hollowed AB, Knowlton N, et al. Climate change impacts on marine ecosystems. Annu Rev Mar Sci. 2012;4(1):11–37.
O’Corry-Crowe G. Climate change and the molecular ecology of Arctic marine mammals. Ecol Appl. 2008;18(sp2)
Phillips CD, Gelatt TS, Patton JC, Bickham JW. Phylogeography of Steller Sea lions: relationships among climate change, effective population size, and genetic diversity. J Mammal. 2011;92(5):1091–104.
Phillips CD, Hoffman JI, George JC, Suydam RS, Huebinger RM, Patton JC, Bickham JW. Molecular insights into the historic demography of bowhead whales: understanding the evolutionary basis of contemporary management practices. Ecol Evol. 2013;3(1):18–37.
Kaschner K, Tittensor DP, Ready J, Gerrodette T, Worm B. Current and future patterns of global marine mammal biodiversity. PLoS One. 2011;6(5):e19653.
Lindquist O. The North Atlantic gray whale (Escherichtius robustus): An historical outline based on Icelandic, Danish-Icelandic, English and Swedish sources dating from ca 1000 AD to 1792, vol. 1. Scotland: University of St. Andrews and Stirling; 2000.
Sumich J. E. Robustus the biology and human history of gray whales. Corvallis: whale cove marine education; 2014. p. 108–21.
Lang AR, Weller DW, Leduc RG, Burdin AM, Brownell RL Jr. Genetic differentiation between western and eastern (Eschrichtius robustus) gray whale populations using microsatellite markers: In: International Whaling Commision. US: University of California; 2010.
Alter SE, Newsome SD, Palumbi SR. Pre-whaling genetic diversity and population ecology in eastern Pacific gray whales: insights from ancient DNA and stable isotopes. PLoS One. 2012;7(5):e35039.
Durban JW, Weller DW, Perryman WL. Gray whale abundance estimates from shore-based counts off California in 2014/2015 and 2015/2016, vol. 4; 2017.
Mate BR, Ilyashenko VY, Bradford AL, Vertyankin VV, Tsidulko GA, Rozhnov VV, Irvine LM. Critically endangered western gray whales migrate to the eastern North Pacific. Biol Lett. 2015;11(4):20150071.
Bradford AL, Weller DW, Ivashchenko YV, Burdin AM, Brownell RL. Anthropogenic scarring of western gray whales (Eschrichtius robustus). Mar Mammal Sci. 2009;25(1):161–75.
Swartz SL, Taylor BL, Rugh DJ. Gray whale Eschrichtius robustus population and stock identity. Mammal Rev. 2006;36(1):66–84.
Eschrichtius robustus (western subpopulation). The IUCN Red List of Threatened Species 2008.
Weller DW, Klimek A, Bradford AL, Calambokidis J, Lang AR, Gisborne B, Burdin AM, Szaniszlo W, Urbán J, Unzueta AG-G. Movements of gray whales between the western and eastern North Pacific. Endanger Species Res. 2012;18(3):193–9.
Norris RD, Turner SK, Hull PM, Ridgwell A. Marine ecosystem responses to Cenozoic global change. Science. 2013;341(6145):492–8.
Darby DA, Polyak L, Bauch HA. Past glacial and interglacial conditions in the Arctic Ocean and marginal seas–a review. Prog Oceanogr. 2006;71(2):129–44.
McMahon BJ, Teeling EC, Höglund J. How and why should we implement genomics into conservation? Evol Appl. 2014;7(9):999–1007.
Hu A, Meehl GA, Otto-Bliesner BL, Waelbroeck C, Han W, Loutre M-F, Lambeck K, Mitrovica JX, Rosenbloom N. Influence of Bering Strait flow and North Atlantic circulation on glacial sea-level changes. Nat Geosci. 2010;3(2):118–21.
Moore SE, Huntington HP. Arctic marine mammals and climate change: impacts and resilience. Ecol Appl. 2008;18(sp2):S157–65.
Rugh DJ, Shelden KE, Schulman-Janiger A. Timing of the gray whale southbound migration. J Cetacean Res Manag. 2001;3:31–9.
Li H, Durbin R. Inference of human population history from individual whole-genome sequences. Nature. 2011;475(7357):493–U484.
Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, Genschoreck T, Webster T, Reich D. Ancient admixture in human history. Genetics. 2012;192(3):1065–93.
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo Z, Pool J, Xu X, Jiang H, Vinckenbosch N, Korneliussen T, et al. Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010;329:75–8.
DeWoody JA, Fernandez NB, Brüniche-Olsen A, Antonides JD, Doyle JM, Miguel PS, Westerman R, Vertyankin VV, Godard-Codding CAJ, Bickham JW. Characterization of the gray whale Eschrichtius robustus genome and a genotyping array based on single-nucleotidepolymorphisms in candidate genes. Biol Bull. 2017;232(3):186–97.
LeDuc RG, Weller DW, Hyde J, Burdin AM, Rosel PE, Brownell RL Jr, Wursig B, Dizon AE. Genetic differences between western and eastern gray whales (Eschrichtius robustus). J Cetac Res Manage. 2002;4(1):1–5.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26:589–95.
Yim H-S, Cho YS, Guang X, Kang SG, Jeong J-Y, Cha S-S, Oh H-M, Lee J-H, Yang EC, Kwon KK. Minke whale genome and aquatic adaptation in cetaceans. Nat Genet. 2014;46(1):88–92.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, Del Angel G, Rivas MA, Hanna M. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.
Auwera GA, Carneiro MO, Hartl C, Poplin R, del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J. From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43(11.10):1–33.
Nielsen R, Paul JS, Albrechtsen A, Song YS. Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011;12(6):443–51.
Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, Hall KP, Evers DJ, Barnes CL, Bignell HR, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456(7218):53–9.
Meynert AM, Ansari M, FitzPatrick DR, Taylor MS. Variant detection sensitivity and biases in whole genome and exome sequencing. BMC Bioinformatics. 2014;15(1):247.
Yim H-S, Cho YS, Guang X, Kang SG, Jeong J-Y, Cha S-S, Oh H-M, Lee J-H, Yang EC, Kwon KK, et al. Minke whale genome and aquatic adaptation in cetaceans. Nat Genet. 2013;46:88.
Durand EY, Patterson N, Reich D, Slatkin M. Testing for ancient admixture between closely related populations. Mol Biol Evol. 2011;28(8):2239–52.
Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH, et al. A draft sequence of the Neandertal genome. Science. 2010;328(5979):710–22.
Kimura M, Crow JF. The number of alleles that can be maintained in a finite population. Genetics. 1964;49(4):725–38.
McQuillan R, Leutenegger A-L, Abdel-Rahman R, Franklin CS, Pericic M, Barac-Lauc L, Smolej-Narancic N, Janicijevic B, Polasek O, Tenesa A. Runs of homozygosity in European populations. Am J Hum Genet. 2008;83(3):359–72.
Bosse M, Megens H-J, Madsen O, Paudel Y, Frantz LAF, Schook LB, Crooijmans RPMA, Groenen MAM. Regions of homozygosity in the porcine genome: consequence of demography and the recombination landscape. PLoS Genet. 2012;8(11):e1003100.
Brüniche-Olsen A, Kellner KK, Anderson CJ, DeWoody JA. Runs of homozygosity have utility in mammalian conservation and evolutionary studies. Evol Appl. in revision.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D, Maller J, Sklar P, de Bakker P, Daly M, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Howrigan DP, Simonson MA, Keller MC. Detecting autozygosity through runs of homozygosity: a comparison of three autozygosity detection algorithms. BMC Genomics. 2011;12:460.
Team RDC. R: A language and environment for statistical computing. Austria: Vienna; 2017.
Lee WC. Testing the genetic relation between two individuals using a panel of frequency-unknown single nucleotide polymorphisms. Ann Hum Genet. 2003;67(6):618–9.
Korneliussen TS, Albrechtsen A, Nielsen R. ANGSD: analysis of next generation sequencing data. BMC Bioinformatics. 2014;15(1):1–13.
Arnason U, Gullberg A, Janke A. Mitogenomic analyses provide new insights into cetacean origin and evolution. Gene. 2004;333:27–34.
Nikaido M, Hamilton H, Makino H, Sasaki T, Takahashi K, Goto M, Kanda N, Pastene LA, Okada N. Baleen whale phylogeny and a past extensive radiation event revealed by SINE insertion analysis. Mol Biol Evol. 2006;23(5):866–73.
Marx FG, Fordyce RE. Baleen boom and bust: a synthesis of mysticete phylogeny, diversity and disparity. R Soc Open Sci. 2015;2(4):140434.
Schiffels S, Durbin R. Inferring human population size and separation history from multiple genome sequences. Nat Genet. 2014;46(8):919–25.
Zhao S, Zheng P, Dong S, Zhan X, Wu Q, Guo X, Hu Y, He W, Zhang S, Fan W. Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation. Nat Genet. 2013;45(1):67–71.
Jackson JA, Baker CS, Vant M, Steel DJ, Medrano-González L, Palumbi SR. Big and slow: phylogenetic estimates of molecular evolution in baleen whales (suborder Mysticeti). Mol Biol Evol. 2009;26(11):2427–40.
Martin AP, Palumbi SR. Body size, metabolic rate, generation time, and the molecular clock. PNAS. 1993;90(9):4087–91.
Alter SE, Rynes E, Palumbi SR. DNA evidence for historic population size and past ecosystem impacts of gray whales. PNAS. 2007;104(38):15162–7.
DWW R, AADW R, Wolman AA. The life history and ecology of the gray whale (Eschrichtius robustus). Stillwater, Okla: American Society of Mammalogists; 1971.
Heppell SS, Caswell H, Crowder LB. Life histories and elasticity patterns: perturbation analysis for species with minimal demographic data. Ecology. 2000;81(3):654–65.
Plummer M, Best N, Cowles K, Vines K. CODA: convergence diagnosis and output analysis for MCMC. R news. 2006;6(1):7–11.
Ready J, Kaschner K, South AB, Eastwood PD, Rees T, Rius J, Agbayani E, Kullander S, Froese R. Predicting the distributions of marine organisms at the global scale. Ecol Model. 2010;221(3):467–78.
Kaschner K: Reviewed distribution maps for Eschrichtius robustus (gray whale), with modelled year 2100 native range map based on IPCC A2 emissions scenario. In., Aug. 2016 edn. www.aquamaps.org.
Haines A: Climate change 2001: the scientific basis. Contribution of working group 1 to the third assessment report of the intergovernmental panel on climate change. JT Houghton, Y Ding, DJ Griggs, M Noguer, PJ van der Winden, X Dai. Cambridge: Cambridge University Press, 2001, pp. 881,£ 34.95 (HB) ISBN: 0-21-01495-6;£ 90.00 (HB) ISBN: 0-521-80767-0. Int J Epidemiol 2003, 32(2):321–321.
Bromham L, Rambaut A, Harvey PH. Determinants of rate variation in mammalian DNA sequence evolution. J Mol Evol. 1996;43(6):610–21.
Chao L, Carr DE. The molecular clock and the relationship between population size and generation time. Evolution. 1993;47(2):688–90.
Frankham R. Genetics and extinction. Biol Conserv. 2005;126(2):131–40.
Prado-Martinez J, Sudmant PH, Kidd JM, Li H, Kelley JL, Lorente-Galdos B, Veeramah KR, Woerner AE, O’Connor TD, Santpere G, et al. Great ape genetic diversity and population history. Nature. 2013;499(7459):471–5.
Kirin M, McQuillan R, Franklin CS, Campbell H, McKeigue PM, Wilson JF. Genomic runs of homozygosity record population history and consanguinity. PLoS One. 2010;5(11):e13996.
Scheinin AP, Kerem D, MacLeod CD, Gazo M, Chicote CA, Castellote M. Gray whale (Eschrichtius robustus) in the Mediterranean Sea: anomalous event or early sign of climate-driven distribution change? Marine Biodiversity Records. 2011;4:e28.
Hu H, Petousi N, Glusman G, Yu Y, Bohlender R, Tashi T, Downie JM, Roach JC, Cole AM, Lorenzo FR, et al. Evolutionary history of Tibetans inferred from whole-genome sequencing. PLoS Genet. 2017;13(4):e1006675.
Fagundes NJR, Ray N, Beaumont M, Neuenschwander S, Salzano FM, Bonatto SL, Excoffier L. Statistical evaluation of alternative models of human evolution. Proc Natl Acad Sci. 2007;104(45):17614–9.
Väli Ü, Einarsson A, Waits L, Ellegren H. To what extent do microsatellite markers reflect genome-wide genetic diversity in natural populations? Mol Ecol. 2008;17(17):3808–17.
Shapiro B, Drummond AJ, Rambaut A, Wilson MC, Matheus PE, Sher AV, Pybus OG, Gilbert MT, Barnes I, Binladen J, et al. Rise and fall of the Beringian steppe bison. Science. 2004;306(5701):1561–5.
Blois JL, McGuire JL, Hadly EA. Small mammal diversity loss in response to late-Pleistocene climatic change. Nature. 2010;465(7299):771–4.
Kishida T. Population history of Antarctic and common minke whales inferred from individual whole-genome sequences. Marine Mammal Science. 2017;33(2):645–52.
Foote AD, Kaschner K, Schultze SE, Garilao C, Ho SYW, Post K, Higham TFG, Stokowska C, van der Es H, Embling CB, et al. Ancient DNA reveals that bowhead whale lineages survived late Pleistocene climate change and habitat shifts. Nat Commun. 2013;4:1677.
Moura AE, Janse van Rensburg C, Pilot M, Tehrani A, Best PB, Thornton M, Plön S, de Bruyn PJN, Worley KC, Gibbs RA, et al. Killer whale nuclear genome and mtDNA reveal widespread population bottleneck during the last glacial maximum. Mol Biol Evol. 2014;31(5):1121–31.
De Bruyn M, Hall BL, Chauke LF, Baroni C, Koch PL, Hoelzel AR. Rapid response of a marine mammal species to Holocene climate and habitat change. PLoS Genet. 2009;5(7):e1000554.
We thank P. A. Morin and two anonymous reviewers for providing constructive comments on an earlier version of the manuscript. A. Albrechtsen provided constructive comments on the analyses. K. Kellner assisted with developing R scripts. T. Rowles provided assistance with permits (National Marine Fisheries Service Office of Protected Resources’ Marine Mammal Health and Stranding Response Program permit 932-1905-MA-009526). A. Aziz, J. Dupont, M. Scott, and M. Swindoll provided support in obtaining the biopsies and associated metadata. The Institute of Ecology and Evolution of the Russian Academy of Sciences, the A.V. Zhirmunsky Institute of Marine Biology Far Eastern Branch, and Oregon State University provided invaluable support for the collection of the Western gray whale samples. C. George and R. Suydam (Department of Wildlife Management, North Slope Borough of Alaska) provided the Eastern gray whale sample.
This work was supported in part by Exxon Neftegas Limited and Sakhalin Energy Investment Company. The funding parties had no part in study design, data collection, analysis or interpreting the results. The content herein is solely the responsibility of the authors and does not necessarily represent the official views of the funding parties.
Availability of data and materials
Genome data were from Genbank acc. No. SRR5495100, SRR5495104, and SRR5495108.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Information on raw reads filtering statistics. Paired-end libraries were sequenced on an Illumina HiSeq 2500. Table S2. Environmental variables used in AQUAMAPS to generate maps of suitable habitat for gray whales during the Holocene. Table S3. The D–test statistic evaluates the number (n) of ABBA and BABA sites (D = (nABBA - nBABA) / (nABBA + nBABA)) and D < 0 means that P1 is more closely related to P3 than to P2, whereas D > 0 indicates that P2 is more closely related to P3 than P1. The significance of the D test was evaluated with a Z-score, where |Z-scores| > 3 was used as the critical value for a significant test. Figure S1. Inferred effective population sizes (Ne) over time. Estimates are averages based on 11 autosomal scaffolds larger than 30 Mb. A substitution rate of a) 10 × 10− 10 bp− 1 year− 1 and b) 1.5 × 10− 10 bp− 1 year− 1 were used. (DOCX 446 kb)
About this article
Cite this article
Brüniche-Olsen, A., Westerman, R., Kazmierczyk, Z. et al. The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences. BMC Evol Biol 18, 87 (2018) doi:10.1186/s12862-018-1204-3
- Demographic history