Molecular insights into the historic demography of bowhead whales: understanding the evolutionary basis of contemporary management practices

Does Effective Population Size Govern Evolutionary Differences in Telomere Length?

Lengths of telomeres vary by an order of magnitude across mammalian species. Similarly, age- and sex-standardized telomere lengths differ by up to 1 kb (14%) across human populations. How to explain these differences? Telomeres play a central role in senescence and aging, and genes that affect telomere length are likely under weak selection (i.e. telomere length is a trait that is subject to nearly neutral evolution). Importantly, natural selection is more effective in large populations than in small populations. Here, we propose that observed differences in telomere length across species and populations are largely due to differences in effective population sizes. In this perspective, we present preliminary evolutionary genetic evidence supporting this hypothesis and highlight the need for more data.

Eco-evolutionary factors that influence its demographic oscillations in Prochilodus costatus (Actinopterygii: Characiformes) populations evidenced through a genetic spatial-temporal evaluation

The human activity impact on wild animal populations is indicated by eco-evolutionary and demographic processes, along with their survival and capacity to evolve; consequently, such data can contribute toward enhancing genetic-based conservation programs. In this context, knowledge on the life-history and the eco-evolutionary processes is required to understand extant patterns of population structure in Prochilodus costatus a Neotropical migratory fish that has been threatened due to loss and fragmentation of its natural habitat since 1960s promoted by the expansion of hydroelectric power plant construction programs. This study evaluated the eco-evolutionary parameters that cause oscillations in the demography and structure of P. costatus populations. An integrated approach was used, including temporal and spatial sampling, next-generation sequencing of eight microsatellite loci, multivariate genetic analysis, and demographic life-history reconstruction. The results provided evidence of the complex interplay of ecological-evolutionary and human-interference events on the life history of this species in the upper basin. In particular, spawning wave behavior might have ecological triggers resulting in an overlapping of distinct genetic generations, and arising distinct migratory and nonmigratory genetic patterns living in the same area. An abrupt decrease in the effective population size of the P. costatus populations in the recent past (1960-80) was likely driven by environment fragmentation promoted by the construction of the Três Marias hydropower dam. The low allelic diversity that resulted from this event is still detected today; thus, active stocking programs are not effective at expanding the genetic diversity of this species in the river basin. Finally, this study highlights the importance of using mixed methods to understand spatial and temporal variation in genetic structure for effective mitigation and conservation programs for threatened species that are directly affected by human actions.

High genomic diversity in the endangered East Greenland Svalbard Barents Sea stock of bowhead whales (Balaena mysticetus)

The East Greenland-Svalbard-Barents Sea (EGSB) bowhead whale stock (Balaena mysticetus) was hunted to near extinction and remains Endangered on the International Union of Conservation of Nature Red List. The intense, temporally extensive hunting pressure may have left the population vulnerable to other perturbations, such as environmental change. However, the lack of genomic baseline data renders it difficult to evaluate the impacts of various potential stressors on this stock. Twelve EGSB bowhead whales sampled in 2017/2018 were re-sequenced and mapped to a previously published draft genome. All individuals were unrelated and void of significant signs of inbreeding, with similar observed and expected homo- and heterozygosity levels. Despite the small population size, mean autosome-wide heterozygosity was 0.00102, which is higher than that of most mammals for which comparable estimates are calculated using the same parameters, and three times higher than a conspecific individual from the Eastern-Canada-West-Greenland bowhead whale stock. Demographic history analyses indicated a continual decrease of N_e from ca. 1.5 million to ca. 250,000 years ago, followed by a slight increase until ca. 100,000 years ago, followed by a rapid decrease in N_e between 50,000 and 10,000 years ago. These estimates are lower than previously suggested based on mitochondrial DNA, but suggested demographic patterns over time are similar.

Mitogenomics and the genetic differentiation of contemporary Balaena mysticetus (Cetacea) from Svalbard

Full mitochondrial genomes were assembled for 12 recently sampled animals from the Svalbard bowhead whale (Balaena mysticetus) stock via high-throughput sequencing data, facilitating analysis of the demographic history of the population for the first time. The Svalbard population has retained noticeable amounts of mitochondrial genome diversity despite extreme historical harvest levels. Haplotype and nucleotide diversities were similar to those estimated earlier for other bowhead whale populations. The reconstructed demographic history was in accordance with a boom–bust scenario, combining a slight Pleistocene population growth 25 000–35 000 years ago and a Holocene decline. Employing a mutation rate of 3.418 × 10–8 substitutions per site per year, the time to the most recent common ancestor for the mitochondrial genomes of the contemporary Svalbard bowhead whales was estimated to be 68 782 (54 353–83 216) years before the present. Based on 370 bp fragments of the D-loop region, significant genetic differentiation was detected between all extant bowhead whale populations across the circumpolar Arctic. Thus, the Svalbard bowhead whales can be regarded as a population with its own genetic legacy.

Evolutionary impacts differ between two exploited populations of northern bottlenose whale (Hyperoodon ampullatus)

Interpretation of conservation status should be informed by an appreciation of genetic diversity, past demography, and overall trends in population size, which contribute to a species' evolutionary potential and resilience to genetic risks. Low genetic diversity can be symptomatic of rapid demographic declines and impose genetic risks to populations, but can also be maintained by natural processes. The northern bottlenose whale Hyperoodon ampullatus has the lowest known mitochondrial diversity of any cetacean and was intensely whaled in the Northwest Atlantic over the last century, but whether exploitation imposed genetic risks that could limit recovery is unknown. We sequenced full mitogenomes and genotyped 37 novel microsatellites for 128 individuals from known areas of abundance in the Scotian Shelf, Northern and Southern Labrador, Davis Strait, and Iceland, and a newly discovered group off Newfoundland. Despite low diversity and shared haplotypes across all regions, both markers supported the Endangered Scotian Shelf population as distinct from the combined northern regions. The genetic affinity of Newfoundland was uncertain, suggesting an area of mixing with no clear population distinction for the region. Demographic reconstruction using mitogenomes suggests that the northern region underwent population expansion following the last glacial maximum, but for the peripheral Scotian Shelf population, a stable demographic trend was followed by a drastic decline over a temporal scale consistent with increasing human activity in the Northwest Atlantic. Low connectivity between the Scotian Shelf and the rest of the Atlantic likely compounded the impact of intensive whaling for this species, potentially imposing genetic risks affecting recovery of this population. We highlight how the combination of historical environmental conditions and modern exploitation of this species has had very different evolutionary impacts on structured populations of northern bottlenose whales across the western North Atlantic.

The inference of gray whale (Eschrichtius robustus) historical population attributes from whole-genome sequences

BACKGROUND

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.

RESULTS

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.

CONCLUSION

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.

Incorporating non-equilibrium dynamics into demographic history inferences of a migratory marine species

Understanding how dispersal and gene flow link geographically separated the populations over evolutionary history is challenging, particularly in migratory marine species. In southern right whales (SRWs, Eubalaena australis), patterns of genetic diversity are likely influenced by the glacial climate cycle and recent history of whaling. Here we use a dataset of mitochondrial DNA (mtDNA) sequences (n = 1327) and nuclear markers (17 microsatellite loci, n = 222) from major wintering grounds to investigate circumpolar population structure, historical demography and effective population size. Analyses of nuclear genetic variation identify two population clusters that correspond to the South Atlantic and Indo-Pacific ocean basins that have similar effective breeder estimates. In contrast, all wintering grounds show significant differentiation for mtDNA, but no sex-biased dispersal was detected using the microsatellite genotypes. An approximate Bayesian computation (ABC) approach with microsatellite markers compared the scenarios with gene flow through time, or isolation and secondary contact between ocean basins, while modelling declines in abundance linked to whaling. Secondary-contact scenarios yield the highest posterior probabilities, implying that populations in different ocean basins were largely isolated and came into secondary contact within the last 25,000 years, but the role of whaling in changes in genetic diversity and gene flow over recent generations could not be resolved. We hypothesise that these findings are driven by factors that promote isolation, such as female philopatry, and factors that could promote dispersal, such as oceanographic changes. These findings highlight the application of ABC approaches to infer the connectivity in mobile species with complex population histories and, currently, low levels of differentiation.

Effective population size and the genetic consequences of commercial whaling on the humpback whales (Megaptera novaeangliae) from Southwestern Atlantic Ocean

Genotypes of 10 microsatellite loci of 420 humpback whales from the Southwestern Atlantic Ocean population were used to estimate for the first time its contemporary effective (N_e) and census (N_c) population sizes and to test the genetic effect of commercial whaling. The results are in agreement with our previous studies that found high genetic diversity for this breeding population. Using an approximate Bayesian computation approach, the scenario of constant N_e was significantly supported over scenarios with moderate to strong size changes during the commercial whaling period. The previous generation N_c (N_e multiplied by 3.6), which should corresponds to the years between around 1980 and 1990, was estimated between ~2,600 and 6,800 whales (point estimate ~4,000), and is broadly compatible with the recent abundance surveys extrapolated to the past using a growth rate of 7.4% per annum. The long-term N_c in the constant scenario (point estimate ~15,000) was broadly compatible (considering the confidence interval) with pre-whaling catch records estimates (point estimate ~25,000). Overall, our results shown that the Southwestern Atlantic Ocean humpback whale population is genetically very diverse and resisted well to the strong population reduction during commercial whaling.

The eastern migratory caribou: the role of genetic introgression in ecotype evolution

Understanding the evolutionary history of contemporary animal groups is essential for conservation and management of endangered species like caribou (Rangifer tarandus). In central Canada, the ranges of two caribou subspecies (barren-ground/woodland caribou) and two woodland caribou ecotypes (boreal/eastern migratory) overlap. Our objectives were to reconstruct the evolutionary history of the eastern migratory ecotype and to assess the potential role of introgression in ecotype evolution. STRUCTURE analyses identified five higher order groups (i.e. three boreal caribou populations, eastern migratory ecotype and barren-ground). The evolutionary history of the eastern migratory ecotype was best explained by an early genetic introgression from barren-ground into a woodland caribou lineage during the Late Pleistocene and subsequent divergence of the eastern migratory ecotype during the Holocene. These results are consistent with the retreat of the Laurentide ice sheet and the colonization of the Hudson Bay coastal areas subsequent to the establishment of forest tundra vegetation approximately 7000 years ago. This historical reconstruction of the eastern migratory ecotype further supports its current classification as a conservation unit, specifically a Designatable Unit, under Canada's Species at Risk Act. These findings have implications for other sub-specific contact zones for caribou and other North American species in conservation unit delineation.

The transcriptome of the bowhead whale Balaena mysticetus reveals adaptations of the longest-lived mammal

Mammals vary dramatically in lifespan, by at least two-orders of magnitude, but the molecular basis for this difference remains largely unknown. The bowhead whale Balaena mysticetus is the longest-lived mammal known, with an estimated maximal lifespan in excess of two hundred years. It is also one of the two largest animals and the most cold-adapted baleen whale species. Here, we report the first genome-wide gene expression analyses of the bowhead whale, based on the de novo assembly of its transcriptome. Bowhead whale or cetacean-specific changes in gene expression were identified in the liver, kidney and heart, and complemented with analyses of positively selected genes. Changes associated with altered insulin signaling and other gene expression patterns could help explain the remarkable longevity of bowhead whales as well as their adaptation to a lipid-rich diet. The data also reveal parallels in candidate longevity adaptations of the bowhead whale, naked mole rat and Brandt's bat. The bowhead whale transcriptome is a valuable resource for the study of this remarkable animal, including the evolution of longevity and its important correlates such as resistance to cancer and other diseases.