Long-term isolation at a low effective population size greatly reduced genetic diversity in Gulf of California fin whales

Phylogenomics and Pervasive Genome-Wide Phylogenetic Discordance Among Fin Whales (Balaenoptera physalus)

Phylogenomics has the power to uncover complex phylogenetic scenarios across the genome. In most cases, no single topology is reflected across the entire genome as the phylogenetic signal differs among genomic regions due to processes, such as introgression and incomplete lineage sorting. Baleen whales are among the largest vertebrates on Earth with a high dispersal potential in a relatively unrestricted habitat, the oceans. The fin whale (Balaenoptera physalus) is one of the most enigmatic baleen whale species, currently divided into four subspecies. It has been a matter of debate whether phylogeographic patterns explain taxonomic variation in fin whales. Here we present a chromosome-level whole genome analysis of the phylogenetic relationships among fin whales from multiple ocean basins. First, we estimated concatenated and consensus phylogenies for both the mitochondrial and nuclear genomes. The consensus phylogenies based upon the autosomal genome uncovered monophyletic clades associated with each ocean basin, aligning with the current understanding of subspecies division. Nevertheless, discordances were detected in the phylogenies based on the Y chromosome, mitochondrial genome, autosomal genome and X chromosome. Furthermore, we detected signs of introgression and pervasive phylogenetic discordance across the autosomal genome. This complex phylogenetic scenario could be explained by a puzzle of introgressive events, not yet documented in fin whales. Similarly, incomplete lineage sorting and low phylogenetic signal could lead to such phylogenetic discordances. Our study reinforces the pitfalls of relying on concatenated or single locus phylogenies to determine taxonomic relationships below the species level by illustrating the underlying nuances that some phylogenetic approaches may fail to capture. We emphasize the significance of accurate taxonomic delineation in fin whales by exploring crucial information revealed through genome-wide assessments.

The genomic footprint of whaling and isolation in fin whale populations

Twentieth century industrial whaling pushed several species to the brink of extinction, with fin whales being the most impacted. However, a small, resident population in the Gulf of California was not targeted by whaling. Here, we analyzed 50 whole-genomes from the Eastern North Pacific (ENP) and Gulf of California (GOC) fin whale populations to investigate their demographic history and the genomic effects of natural and human-induced bottlenecks. We show that the two populations diverged ~16,000 years ago, after which the ENP population expanded and then suffered a 99% reduction in effective size during the whaling period. In contrast, the GOC population remained small and isolated, receiving less than one migrant per generation. However, this low level of migration has been crucial for maintaining its viability. Our study exposes the severity of whaling, emphasizes the importance of migration, and demonstrates the use of genome-based analyses and simulations to inform conservation strategies.

Fin whale acoustic populations present in New Zealand waters: Description of song types, occurrence and seasonality using passive acoustic monitoring

Southern fin whales (Balaenoptera physalus) are known to migrate from the Antarctic to mid-latitudes during winter for breeding, but the occurrence and distribution of this species is not well known in the waters around New Zealand. The 'doublet' calls are one of the main calls emitted specifically by fin whales and repeated in a regular pattern, which make the acoustic detection of these calls relevant to detect the presence of fin whales. Using a signal processing algorithm to detect 'doublet' calls emitted by fin whales, we studied the occurrence, characteristics and seasonality of these 'doublet' calls in two regions around New Zealand; Cook Strait in 2016/2017 and offshore Gisborne in 2014/2015. The call detection procedure consisted of binarization of the spectrogram and a cross-correlation between the binarized spectrogram and a template of binarized 'doublet' calls spectrogram. A binarization threshold for the data spectrograms and a cross correlation threshold were then determined through multiple trials on a training dataset and a Receiver Operating Characteristics (ROC) curve. Fin whale 'doublet' calls occurred on the east side of New Zealand's Cook Strait during austral winter, specifically in June 2017 and offshore Gisborne in June-August 2014. No 'doublet' calls were detected on the west side of Cook Strait. The 'doublet' calls' Inter-Note Interval (INI) was similar in both datasets. However, there was a difference in alternation of the mean frequency for both HF components of 'doublet' calls in Cook Strait and Gisborne. As the song types were compared with those previously described in the literature, our findings suggest that some fin whales wintering in New Zealand waters may be part of a broader 'acoustic population' whose range extends west to southern Australia and south to Antarctica.

Characterization of 25 new microsatellite markers for the fin whale (Balaenoptera physalus) and cross-species amplification in other cetaceans

Cetaceans are large mammals widely distributed on Earth. The fin whale, Balaenoptera physalus, is the second largest living animal. In the 20th century, commercial whaling reduced its global population by 70%, and in the Mediterranean Sea not only was their overall population depleted but the migration between the Mediterranean Sea and the Atlantic Ocean was reduced. Previous genetic studies identified isolation between these two regions, with a limited gene-flow between these adjacent populations based on nuclear and mitochondrial markers. However, only limited information exists for the Mediterranean population as genetic diversity and abundance trends are still unknown. In this study, 39 highly polymorphic microsatellite markers were tested, including 25 markers developed de novo together with 14 markers previously published. An average allelic diversity of 8.3 alleles per locus was reported, ranging from 3 to 15 alleles per locus, for B. physalus. Expected heterozygosity was variable among loci and ranged from 0.34 to 0.91. Only two markers in the new set were significantly deviant from the Hardy Weinberg equilibrium. Cross-species amplification was tested in four other cetacean species. A total of 27 markers were successfully amplified in the four species (Balaenoptera acutorostrata, Megaptera novaeangliae, Physeter macrocephalus and Globicephala melas). A multivariate analysis on the multilocus genotypes successfully discriminated the five species. This new set of microsatellite markers will not only provide a useful tool to identify and understand the genetic diversity and the evolution of the B. physalus population, but it will also be relevant for other cetacean species, and will allow further parentage analyses. Eventually, this new set of microsatellite markers will provide critical data that will shed light on important biological data within a conservation perspective.