Strong and lasting impacts of past global warming on baleen whales and their prey.

Global warming is affecting the population dynamics and trophic interactions across a wide range of ecosystems and habitats. Translating these real-time effects into their long-term consequences remains a challenge. The rapid and extreme warming period that occurred after the Last Glacial Maximum (LGM) during the Pleistocene-Holocene transition (7-12 thousand years ago) provides an opportunity to gain insights into the long-term responses of natural populations to periods with global warming. The effects of this post-LGM warming period have been assessed in many terrestrial taxa, whereas insights into the impacts of rapid global warming on marine taxa remain limited, especially for megafauna. In order to understand how large-scale climate fluctuations during the post-LGM affected baleen whales and their prey, we conducted an extensive, large-scale analysis of the long-term effects of the post-LGM warming on abundance and inter-ocean connectivity in eight baleen whale and seven prey (fish and invertebrates) species across the Southern and the North Atlantic Ocean; two ocean basins that differ in key oceanographic features. The analysis was based upon 7032 mitochondrial DNA sequences as well as genome-wide DNA sequence variation in 100 individuals. The estimated temporal changes in genetic diversity during the last 30,000 years indicated that most baleen whale populations underwent post-LGM expansions in both ocean basins. The increase in baleen whale abundance during the Holocene was associated with simultaneous changes in their prey and climate. Highly correlated, synchronized and exponential increases in abundance in both baleen whales and their prey in the Southern Ocean were indicative of a dramatic increase in ocean productivity. In contrast, the demographic fluctuations observed in baleen whales and their prey in the North Atlantic Ocean were subtle, varying across taxa and time. Perhaps most important was the observation that the ocean-wide expansions and decreases in abundance that were initiated by the post-LGM global warming, continued for millennia after global temperatures stabilized, reflecting persistent, long-lasting impacts of global warming on marine fauna.

Whole genome resequencing reveals diagnostic markers for investigating global migration and hybridization between minke whale species.

BACKGROUND: In the marine environment, where there are few absolute physical barriers, contemporary contact between previously isolated species can occur across great distances, and in some cases, may be inter-oceanic. An example of this can be seen in the minke whale species complex. Antarctic minke whales are genetically and morphologically distinct from the common minke found in the north Atlantic and Pacific oceans, and the two species are estimated to have been isolated from each other for 5 million years or more. Recent atypical migrations from the southern to the northern hemisphere have been documented and fertile hybrids and back-crossed individuals between both species have also been identified. However, it is not known whether this represents a contemporary event, potentially driven by ecosystem changes in the Antarctic, or a sporadic occurrence happening over an evolutionary time-scale. We successfully used whole genome resequencing to identify a panel of diagnostic SNPs which now enable us address this evolutionary question. RESULTS: A large number of SNPs displaying fixed or nearly fixed allele frequency differences among the minke whale species were identified from the sequence data. Five panels of putatively diagnostic markers were established on a genotyping platform for validation of allele frequencies; two panels (26 and 24 SNPs) separating the two species of minke whale, and three panels (22, 23, and 24 SNPs) differentiating the three subspecies of common minke whale. The panels were validated against a set of reference samples, demonstrating the ability to accurately identify back-crossed whales up to three generations. CONCLUSIONS: This work has resulted in the development of a panel of novel diagnostic genetic markers to address inter-oceanic and global contact among the genetically isolated minke whale species and sub-species. These markers, including a globally relevant genetic reference data set for this species complex, are now openly available for researchers interested in identifying other potential whale hybrids in the world's oceans. The approach used here, combining whole genome resequencing and high-throughput genotyping, represents a universal approach to develop similar tools for other species and population complexes.

Geographical distribution of close kin in southern right whales on feeding grounds.

This study investigated the close kinship structure of southern right whales on feeding grounds during austral summer seasons. The study was based on biopsy samples of 171 individual whales, which were genotyped with 14 microsatellite DNA loci. Kinship was investigated by using the LOD (Log Odds) score, a relatedness index for a pair of genotypes. Based on a cut-off point of LODPO > 6, which was chosen to balance false positives and negatives, a total of 28 dyads were inferred. Among these, 25 were classified as parent-offspring pairs. Additional genetic (mitochondrial DNA haplotypes) and biological (estimated body length, sex) data were used to provide additional information on the inferred close kin pairs. The elapsed time between sampling varied from 0 (close kin detected in the same austral summer season) to 17 years. All the kin pairs occurred within the Antarctic Indo sector (85°-135°E) and no pair occurred between whales within and outside of this sector. Six pairs were between individuals in high (Antarctic) and lower latitudes. Results of the present analysis on kinship are consistent with the views that whales in the Indo sector of the Antarctic are related with the breeding ground in Southwest Australia, and that whales from this population can occupy different feeding grounds. The present study has the potential to contribute to the conservation of the southern right whales through the monitoring of important population parameters such as population sizes and growth rate, in addition to assist the interpretation of stock structure derived from standard population genetic analyses.

Hybrids between common and Antarctic minke whales are fertile and can back-cross.

BACKGROUND: Minke whales are separated into two genetically distinct species: the Antarctic minke whale found in the southern hemisphere, and the common minke whale which is cosmopolitan. The common minke whale is further divided into three allopatric sub-species found in the North Pacific, southern hemisphere, and the North Atlantic. Here, we aimed to identify the genetic ancestry of a pregnant female minke whale captured in the North Atlantic in 2010, and her fetus, using data from the mtDNA control region, 11 microsatellite loci and a sex determining marker. RESULTS: All statistical parameters demonstrated that the mother was a hybrid displaying maternal and paternal contribution from North Atlantic common and Antarctic minke whales respectively. Her female fetus displayed greater genetic similarity to North Atlantic common minke whales than herself, strongly suggesting that the hybrid mother had paired with a North Atlantic common minke whale. CONCLUSION: This study clearly demonstrates, for the first time, that hybrids between minke whale species may be fertile, and that they can back-cross. Whether contact between these species represents a contemporary event linked with documented recent changes in the Antarctic ecosystem, or has occurred at a low frequency over many years, remains open.

Molecular cloning of urea transporters from the kidneys of baleen and toothed whales.

Urea transport in the kidney is important for the production of concentrated urine. This process is mediated by urea transporters (UTs) encoded by two genes, UT-A (Slc14a2) and UT-B (Slc14a1). Our previous study demonstrated that cetaceans produce highly concentrated urine than terrestrial mammals, and that baleen whales showed higher concentrations of urinary urea than sperm whales. Therefore, we hypothesized that cetaceans have unique actions of UTs to maintain fluid homeostasis in marine habitat. Kidney samples of common minke (Balaenoptera acutorostrata), sei (B. borealis), Bryde's (B. brydei) and sperm whales (Physeter macrocephalus) were obtained to determine the nucleotide sequences of mRNAs encoding UT. The sequences of 2.5-kb cDNAs encode 397-amino acid proteins, which are 90-94% identical to the mammalian UT-A2s. Two putative glycosylation sites are conserved between the whales and the terrestrial mammals, whereas consensus sites for protein kinases are not completely conserved; only a single protein kinase A consensus site was identified in the whale UT-A2s. Two protein kinase C consensus sites are present in the baleen whale UT-A2s, however, a single protein kinase C consensus site was identified in the sperm whale UT-A2. These different phosphorylation sites of whale UT-A2s may result in the high concentrations of urinary urea in whales, by reflecting their urea permeability.

Phylogenetic relationships among cetaceans revealed by Y-chromosome sequences.

The Y chromosome has recently come into the spotlight as a new and efficient genetic marker for tracing paternal lineages. We reconstructed cetacean phylogeny using a 1.7-kbp fragment of the non-recombining Y chromosome (NRY), including the SRY gene and a flanking non-coding region. The topology of the Y-chromosome tree is robust to various methods of analysis and exhibits high branch-support values, possibly due to the absence of recombination, small effective population size, and low homoplasy. The Y-chromosome tree indicates monophyly of each suborder, Mysticeti and Odontoceti, with high branch support values (BS> or =86%; PP> or =98%). In the Odontoceti clade, three superfamilies, Physeteroidea, Ziphioidea, and Delphinoidea, diverged soon after the split between Mysticeti and Odontoceti. Our analysis allows resolution of this rapid radiation and indicates that Physeteroidea is basal in the Odontoceti clade (BS, 99%; PP, 100%; MBS, 61%). The major split within the superfamily Delphinoidea is between the Delphinidae clade and the Monodontidae+ Phocoenidae clade. The phylogenetic relationships among delphinid species are ambiguous, probably because of the rapid radiation of this family. In the Mysticeti clade, the first major split is between Balaenidae and Balaenopteridae; within Balaenopteridae, a Balaenoptera acutorostrata+B. bonaerensis (minke whales) clade forms a sister clade with the other balaenopterid species. Megaptera novaeangliae is nested within Balaenoptera, making the latter paraphyletic. The low homoplasy exhibited by the Y-chromosome data presented here suggests that an extended data set incorporating longer sequences would provide better resolution of cetacean lower-level pylogeny.

Sequence variation in the Tbx4 gene in marine mammals.

The amino-acid sequences of the T-domain region of the Tbx4 gene, which is required for hindlimb development, are 100% identical in humans and mice. Cetaceans have lost most of their hindlimb structure, although hindlimb buds are present in very early cetacean embryos. To examine whether the Tbx4 gene has the same function in cetaceans as in other mammals, we analyzed Tbx4 sequences from cetaceans, dugong, artiodactyls and marine carnivores. A total of 39 primers were designed using human and dog Tbx4 nucleotide sequences. Exons 3, 4, 5, 6, 7, and 8 of the Tbx4 genes from cetaceans, artiodactyls, and marine carnivores were sequenced. Non-synonymous substitution sites were detected in the T-domain regions from some cetacean species, but were not detected in those from artiodactyls, the dugong, or the carnivores. The C-terminal regions contained a number of non-synonymous substitutions. Although some indels were present, they were in groups of three nucleotides and therefore did not cause frame shifts. The dN/dS values for the T-domain and C-terminal regions of the cetacean and artiodactylous Tbx4 genes were much lower than 1, indicating that the Tbx4 gene maintains it function in cetaceans, although full expression leading to hindlimb development is suppressed.

Genetic variation of the MHC DQB locus in the finless porpoise (Neophocaena phocaenoides).

The Major Histocompatibility Complex (MHC) is a large multigene coding for glycoproteins that play a key role in the initiation of immune responses in vertebrates. The exon 2 region of the MHC DQB locus was analyzed using 160 finless porpoises from 5 populations in Japanese waters. The 5 populations were based on a previous mitochondrial DNA control region analysis, which showed distinct geographical separation. Eight DQB alleles were detected, and the geographical distribution of the alleles indicated that most of them are shared among the populations. Heterozygosity of the DQB alleles in each population ranged from 0.55 to 0.78, and for all 5 populations was 0.78. Low MHC variability is not a common feature in marine mammals, but the finless porpoise populations inhabiting coastal waters had a relatively high MHC heterozygosity. Balancing selection in the MHC DQB alleles of the finless porpoise was indicated by the higher rate of nonsynonymous than synonymous substitutions for PBR; however, an excess of hetrozygotes compared to expectation was not observed. This suggests that the MHC DQB locus in the finless porpoise may have been under balancing selection for a long evolutionary time period, and is influenced by genetic drift beyond the effect of balancing selection for short time periods in small local populations.

Plasma and urine levels of electrolytes, urea and steroid hormones involved in osmoregulation of cetaceans.

Cetaceans are well adapted to their hyperosmotic environment by properly developed osmoregulatory ability. A question here is how they regulate water and mineral balances in marine habitats. In the present study, we determined blood and urine levels of various chemicals involved in osmoregulation, compared them with those in artiodactyls, and characterized the values in the whales. Blood and urine samples obtained from baleen whales of common minke (Balaenoptera acutorostrata), sei (B. borealis), and Bryde's whales (B. brydei), and toothed whales of sperm whales (Physeter macrocephalus) were analyzed for osmolality, major electrolytes, urea, steroid hormones and glucose. The urine osmolality and Na(+) concentrations in the cetaceans were much higher than those in the cattle. Furthermore, the cetaceans had 5 to 11-fold urea in plasma than the cattle, and 2 to 4-fold urea in urine. There were no significant difference in the plasma concentrations of corticosteroids between the cetaceans and the cattle. The present results indicate that the osmoregulatory parameters seem to be not affected by the reproductive stage and sex steroid hormones. The concentrations of urea in plasma and urine of the baleen whales were higher than those of the sperm whales, indicating a possibility that their osmoregulatory mechanisms may be correlated to their feeding habits. The present results suggest that cetaceans have unique osmoregulatory mechanisms by which they excrete strongly hypertonic urine to maintain fluid homeostasis in marine habitats.

Investigating population genetic structure in a highly mobile marine organism: the minke whale Balaenoptera acutorostrata acutorostrata in the North East Atlantic.

Inferring the number of genetically distinct populations and their levels of connectivity is of key importance for the sustainable management and conservation of wildlife. This represents an extra challenge in the marine environment where there are few physical barriers to gene-flow, and populations may overlap in time and space. Several studies have investigated the population genetic structure within the North Atlantic minke whale with contrasting results. In order to address this issue, we analyzed ten microsatellite loci and 331 bp of the mitochondrial D-loop on 2990 whales sampled in the North East Atlantic in the period 2004 and 2007-2011. The primary findings were: (1) No spatial or temporal genetic differentiations were observed for either class of genetic marker. (2) mtDNA identified three distinct mitochondrial lineages without any underlying geographical pattern. (3) Nuclear markers showed evidence of a single panmictic population in the NE Atlantic according STRUCTURE's highest average likelihood found at K = 1. (4) When K = 2 was accepted, based on the Evanno's test, whales were divided into two more or less equally sized groups that showed significant genetic differentiation between them but without any sign of underlying geographic pattern. However, mtDNA for these individuals did not corroborate the differentiation. (5) In order to further evaluate the potential for cryptic structuring, a set of 100 in silico generated panmictic populations was examined using the same procedures as above showing genetic differentiation between two artificially divided groups, similar to the aforementioned observations. This demonstrates that clustering methods may spuriously reveal cryptic genetic structure. Based upon these data, we find no evidence to support the existence of spatial or cryptic population genetic structure of minke whales within the NE Atlantic. However, in order to conclusively evaluate population structure within this highly mobile species, more markers will be required.

Migration of Antarctic minke whales to the Arctic.

The Antarctic minke whale (Balaenoptera bonaerensis), and the common minke whale found in the North Atlantic (Balaenoptera acutorostrata acutorostrata), undertake synchronized seasonal migrations to feeding areas at their respective poles during spring, and to the tropics in the autumn where they overwinter. Differences in the timing of seasons between hemispheres prevent these species from mixing. Here, based upon analysis of mitochondrial and microsatellite DNA profiles, we report the observation of a single B. bonaerensis in 1996, and a hybrid with maternal contribution from B. bonaerensis in 2007, in the Arctic Northeast Atlantic. Paternal contribution was not conclusively resolved. This is the first documentation of B. bonaerensis north of the tropics, and, the first documentation of hybridization between minke whale species.

Radiation and speciation of pelagic organisms during periods of global warming: the case of the common minke whale, Balaenoptera acutorostrata.

How do populations of highly mobile species inhabiting open environments become reproductively isolated and evolve into new species? We test the hypothesis that elevated ocean-surface temperatures can facilitate allopatry among pelagic populations and thus promote speciation. Oceanographic modelling has shown that increasing surface temperatures cause localization and reduction of upwelling, leading to fragmentation of feeding areas critical to pelagic species. We test our hypothesis by genetic analyses of populations of two closely related baleen whales, the Antarctic minke whale (Balaenoptera bonaerensis) and common minke whale (Balaenoptera acutorostrata) whose current distributions and migration patterns extent are largely determined by areas of consistent upwelling with high primary production. Phylogeographic and population genetic analyses of mitochondrial DNA control-region nucleotide sequences collected from 467 whales sampled in four different ocean basins were employed to infer the evolutionary relationship among populations of B. acutorostrata by rooting an intraspecific phylogeny with a population of B. bonaerensis. Our findings suggest that the two species diverged in the Southern Hemisphere less than 5 million years ago (Ma). This estimate places the speciation event during a period of extended global warming in the Pliocene. We propose that elevated ocean temperatures in the period facilitated allopatric speciation by disrupting the continuous belt of upwelling maintained by the Antarctic Circumpolar Current. Our analyses revealed that the current populations of B. acutorostrata likely diverged after the Pliocene some 1.5 Ma when global temperatures had decreased and presumably coinciding with the re-establishment of the polar-equatorial temperature gradient that ultimately drives upwelling. In most population samples, we detected genetic signatures of exponential population expansions, consistent with the notion of increasing carrying capacity after the Pliocene. Our hypothesis that prolonged periods of global warming facilitate speciation in pelagic marine species that depend on upwelling should be tested by comparative analyses in other pelagic species.

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Baleen whale phylogeny and a past extensive radiation event revealed by SINE insertion analysis.

Baleen whales (suborder Mysticeti) comprise 11 extant species that are classified into four families. Although several phylogenetic hypotheses about these taxa have been proposed, their phylogenetic relationships remain confused. We addressed this problem using short interspersed repetitive element (SINE) insertion data, which now are regarded as almost ideal shared, derived characters at the molecular level. We reconstructed the phylogenetic relationships of baleen whales by characterizing 36 informative SINE loci. One of the intriguing conclusions is that balaenopterids and eschrichtiids radiated very rapidly during a very short evolutionary period. During this period, speciation occurred in balaenopterids and eschrichtiids while newly inserted SINE loci remains polymorphic. Later on, these SINEs were sorted incompletely into each lineage. Thus, there are now inconsistencies among species regarding the presence or absence of a given SINE. This is in sharp contrast to the phylogeny of toothed whales, for which no SINE inconsistencies have been found. Furthermore, we found monophyletic groupings between humpback and fin whales as well as between (sei+Bryde's) whales and blue whales, both of which have not previously been recognized. The comprehensive SINE insertion data, together with the mitochondrial DNA phylogeny that was recently completed (Sasaki, T., M. Nikaido, H. Healy et al. 2005. Mitochondrial phylogenetics and evolution of mysticete whales. Syst. Biol. 56:77-90; Rychel, A. L., T. W. Reeder, and A. Berta. 2004. Phylogeny of mysticete whales based on mitochondrial and nuclear data. Mol. Phylogenet. Evol. 32:892-901), provide a nearly complete picture of the evolutionary history of baleen whales.

Mitochondrial phylogenetics and evolution of mysticete whales.

The phylogenetic relationships among baleen whales (Order: Cetacea) remain uncertain despite extensive research in cetacean molecular phylogenetics and a potential morphological sample size of over 2 million animals harvested. Questions remain regarding the number of species and the monophyly of genera, as well as higher order relationships. Here, we approach mysticete phylogeny with complete mitochondrial genome sequence analysis. We determined complete mtDNA sequences of 10 extant Mysticeti species, inferred their phylogenetic relationships, and estimated node divergence times. The mtDNA sequence analysis concurs with previous molecular studies in the ordering of the principal branches, with Balaenidae (right whales) as sister to all other mysticetes base, followed by Neobalaenidae (pygmy right whale), Eschrichtiidae (gray whale), and finally Balaenopteridae (rorquals + humpback whale). The mtDNA analysis further suggests that four lineages exist within the clade of Eschrichtiidae + Balaenopteridae, including a sister relationship between the humpback and fin whales, and a monophyletic group formed by the blue, sei, and Bryde's whales, each of which represents a newly recognized phylogenetic relationship in Mysticeti. We also estimated the divergence times of all extant mysticete species, accounting for evolutionary rate heterogeneity among lineages. When the mtDNA divergence estimates are compared with the mysticete fossil record, several lineages have molecular divergence estimates strikingly older than indicated by paleontological data. We suggest this discrepancy reflects both a large amount of ancestral polymorphism and long generation times of ancestral baleen whale populations.