Long-term stability in the circumpolar foraging range of a Southern Ocean predator between the eras of whaling and rapid climate change.

Assessing environmental changes in Southern Ocean ecosystems is difficult due to its remoteness and data sparsity. Monitoring marine predators that respond rapidly to environmental variation may enable us to track anthropogenic effects on ecosystems. Yet, many long-term datasets of marine predators are incomplete because they are spatially constrained and/or track ecosystems already modified by industrial fishing and whaling in the latter half of the 20th century. Here, we assess the contemporary offshore distribution of a wide-ranging marine predator, the southern right whale (SRW, Eubalaena australis), that forages on copepods and krill from ~30°S to the Antarctic ice edge (>60°S). We analyzed carbon and nitrogen isotope values of 1,002 skin samples from six genetically distinct SRW populations using a customized assignment approach that accounts for temporal and spatial variation in the Southern Ocean phytoplankton isoscape. Over the past three decades, SRWs increased their use of mid-latitude foraging grounds in the south Atlantic and southwest (SW) Indian oceans in the late austral summer and autumn and slightly increased their use of high-latitude (>60°S) foraging grounds in the SW Pacific, coincident with observed changes in prey distribution and abundance on a circumpolar scale. Comparing foraging assignments with whaling records since the 18th century showed remarkable stability in use of mid-latitude foraging areas. We attribute this consistency across four centuries to the physical stability of ocean fronts and resulting productivity in mid-latitude ecosystems of the Southern Ocean compared with polar regions that may be more influenced by recent climate change.

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability.

Pleistocene environmental changes are generally assumed to have dramatically affected species' demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120-14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome-wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present-day levels during the Last Glacial Maximum (LGM; 26-19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species' demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species' demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

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.

The population genomic structure of green turtles (Chelonia mydas) suggests a warm-water corridor for tropical marine fauna between the Atlantic and Indian oceans during the last interglacial.

The occasional westward transport of warm water of the Agulhas Current, "Agulhas leakage", around southern Africa has been suggested to facilitate tropical marine connectivity between the Atlantic and Indian oceans, but the "Agulhas leakage" hypothesis does not explain the signatures of eastward gene flow observed in many tropical marine fauna. We investigated an alternative hypothesis: the establishment of a warm-water corridor during comparatively warm interglacial periods. The "warm-water corridor" hypothesis was investigated by studying the population genomic structure of Atlantic and Southwest Indian Ocean green turtles (N = 27) using 12,035 genome-wide single nucleotide polymorphisms (SNPs) obtained via ddRAD sequencing. Model-based and multivariate clustering suggested a hierarchical population structure with two main Atlantic and Southwest Indian Ocean clusters, and a Caribbean and East Atlantic sub-cluster nested within the Atlantic cluster. Coalescent-based model selection supported a model where Southwest Indian Ocean and Caribbean populations diverged from the East Atlantic population during the transition from the last interglacial period (130-115 thousand years ago; kya) to the last glacial period (115-90 kya). The onset of the last glaciation appeared to isolate Atlantic and Southwest Indian Ocean green turtles into three refugia, which subsequently came into secondary contact in the Caribbean and Southwest Indian Ocean when global temperatures increased after the Last Glacial Maximum. Our findings support the establishment of a warm-water corridor facilitating tropical marine connectivity between the Atlantic and Southwest Indian Ocean during warm interglacials.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae).

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

Decadal shift in foraging strategy of a migratory southern ocean predator.

Rapid anthropogenic environmental change is expected to impact a host of ecological parameters in Southern Ocean ecosystems. Of critical concern are the consequences of these changes on the range of species that show fidelity to migratory destinations, as philopatry is hypothesized to help or hinder adaptation to climate change depending on the circumstances. Many baleen whales show philopatry to feeding grounds and are also capital breeders that meet migratory and reproductive costs through seasonal energy intake. Southern right whales (Eubalaena australis, SRWs) are capital breeders that have a strong relationship between reproductive output and foraging success. The population dynamics of South Africa's population of SRWs are characterized by two distinct periods: the 1990s, a period of high calving rates; and the late 2010s, a period associated with lowered calving rates. Here we use analyses of stable carbon (δ13 C) and nitrogen (δ15 N) isotope values from SRW biopsy samples (n = 122) collected during these two distinct periods to investigate foraging ecology of the South African population of SRWs over a time period coincident with the demographic shift. We show that South African SRWs underwent a dramatic northward shift, and diversification, in foraging strategy from 1990s to 2010s. Bayesian mixing model results suggest that during the 1990s, South African SRWs foraged on prey isotopically similar to South Georgia/Islas Georgias del Sur krill. In contrast, in the 2010s, South African SRWs foraged on prey isotopically consistent with the waters of the Subtropical Convergence, Polar Front and Marion Island. We hypothesize that this shift represents a response to changes in preferred habitat or prey, for example, the decrease in abundance and southward range contraction of Antarctic krill. By linking reproductive decline to changing foraging strategies for the first time in SRWs, we show that altering foraging strategies may not be sufficient to adapt to a changing ocean.

Fin whale (Balaenoptera physalus) mitogenomics: A cautionary tale of defining sub-species from mitochondrial sequence monophyly.

The advent of massive parallel sequencing technologies has resulted in an increase of studies based upon complete mitochondrial genome DNA sequences that revisit the taxonomic status within and among species. Spatially distinct monophyly in such mitogenomic genealogies, i.e., the sharing of a recent common ancestor among con-specific samples collected in the same region has been viewed as evidence for subspecies. Several recent studies in cetaceans have employed this criterion to suggest subsequent intraspecific taxonomic revisions. We reason that employing intra-specific, spatially distinct monophyly at non-recombining, clonally inherited genomes is an unsatisfactory criterion for defining subspecies based upon theoretical (genetic drift) and practical (sampling effort) arguments. This point was illustrated by a re-analysis of a global mitogenomic assessment of fin whales, Balaenoptera physalus spp., published by Archer et al. (2013), which proposed to further subdivide the Northern Hemisphere fin whale subspecies, B. p. physalus. The proposed revision was based upon the detection of spatially distinct monophyly among North Atlantic and North Pacific fin whales in a genealogy based upon complete mitochondrial genome DNA sequences. The extended analysis conducted in this study (1676 mitochondrial control region, 162 complete mitochondrial genome DNA sequences and 20 microsatellite loci genotyped in 380 samples) revealed that the apparent monophyly among North Atlantic fin whales reported by Archer et al. (2013) to be due to low sample sizes. In conclusion, defining sub-species from monophyly (i.e., the absence of para- or polyphyly) can lead to erroneous conclusions due to relatively "trivial" aspects, such as sampling. Basic population genetic processes (i.e., genetic drift and migration) also affect the time to the most recent common ancestor and hence the probability that individuals in a sample are monophyletic.

Evaluation of fecal hormones for noninvasive research on reproduction and stress in humpback whales (Megaptera novaeangliae).

Fecal hormone analysis shows high potential for noninvasive assessment of population-level patterns in stress and reproduction of marine mammals. However, the marine environment presents unique challenges for fecal sample collection. Data are still lacking on collection methodology and assay validations for most species, particularly for those mysticete whales that have variable diets. In this study we tested collection techniques for fecal samples of free-swimming humpback whales (Megaptera novaeangliae), and validated immunoassays for five steroid and thyroid hormones. Resulting data were used for preliminary physiological validations, i.e., comparisons to independently confirmed sex and reproductive state. Pregnant females had significantly higher fecal progestins and glucocorticoids than did other demographic categories of whales. Two possible cases of previously undetected pregnancies were noted. Males had significantly higher fecal testosterone metabolites than nonpregnant females. Fecal glucocorticoids were significantly elevated in pregnant females and mature males compared to nonpregnant females. Calf fecal samples had elevated concentrations of all fecal hormones. Fecal thyroid hormones showed a significant seasonal decline from spring to summer. Though sample sizes were small, and sampling was necessarily opportunistic, these patterns indicate that noninvasive fecal hormone analysis may facilitate studies of reproduction, stress and potentially energetics in humpback whales.

Order within chaos: potential migratory strategies and individual associations in fin whales feeding off Iceland.

BACKGROUND: The life cycle of most baleen whales involves annual migrations from low-latitude breeding grounds to high latitude feeding grounds. In most species, these migrations are traditionally considered to be carried out according to information acquired through vertical social learning during the first months of life and made individually. However, some recent studies have suggested a more complex scenario, particularly for the species of the Balaenoptera genus. METHODS: Here, we studied the variation of δ15N and δ13C values along the growth axis of the baleen plate from 24 fin whales feeding off western Iceland to delve into their pattern of movements and to identify potential associations between individuals. The segment of baleen plate analyzed informed about at least two complete migratory cycles. We performed cluster analyses through two different methodologies and, whenever possible, we genotyped 20 microsatellite loci to determine potential existence of kinship. RESULTS: Results of the of δ15N and δ13C values agree with a dispersion strategy in the winter breeding grounds. However, and despite the overall large variability, several pairs or groups of individuals with no kinship showed highly similar isotopic patterns for two consecutive years for both δ15N and δ13C values. CONCLUSIONS: Our results suggest that, notably, some whales without kinship share the same migratory regime and destinations. We hypothesize that this could reflect either: (i) the sharing of particularly beneficial migratory regimes, and/or (ii) long-term association between individuals.

Inferring recent historic abundance from current genetic diversity.

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

Wild pedigrees inform mutation rates and historic abundance in baleen whales.

Phylogeny-based estimates suggesting a low germline mutation rate (µ) in baleen whales have influenced research ranging from assessments of whaling impacts to evolutionary cancer biology. We estimated µ directly from pedigrees in four baleen whale species for both the mitochondrial control region and nuclear genome. The results suggest values higher than those obtained through phylogeny-based estimates and similar to pedigree-based values for primates and toothed whales. Applying our estimate of µ reduces previous genetic-based estimates of preexploitation whale abundance by 86% and suggests that µ cannot explain low cancer rates in gigantic mammals. Our study shows that it is feasible to estimate µ directly from pedigrees in natural populations, with wide-ranging implications for ecological and evolutionary research.

Empirical evaluation of humpback whale telomere length estimates; quality control and factors causing variability in the singleplex and multiplex qPCR methods.

BACKGROUND: Telomeres, the protective cap of chromosomes, have emerged as powerful markers of biological age and life history in model and non-model species. The qPCR method for telomere length estimation is one of the most common methods for telomere length estimation, but has received recent critique for being too error-prone and yielding unreliable results. This critique coincides with an increasing awareness of the potentials and limitations of the qPCR technique in general and the proposal of a general set of guidelines (MIQE) for standardization of experimental, analytical, and reporting steps of qPCR. In order to evaluate the utility of the qPCR method for telomere length estimation in non-model species, we carried out four different qPCR assays directed at humpback whale telomeres, and subsequently performed a rigorous quality control to evaluate the performance of each assay. RESULTS: Performance differed substantially among assays and only one assay was found useful for telomere length estimation in humpback whales. The most notable factors causing these inter-assay differences were primer design and choice of using singleplex or multiplex assays. Inferred amplification efficiencies differed by up to 40% depending on assay and quantification method, however this variation only affected telomere length estimates in the worst performing assays. CONCLUSION: Our results suggest that seemingly well performing qPCR assays may contain biases that will only be detected by extensive quality control. Moreover, we show that the qPCR method for telomere length estimation can be highly precise and accurate, and thus suitable for telomere measurement in non-model species, if effort is devoted to optimization at all experimental and analytical steps. We conclude by highlighting a set of quality controls which may serve for further standardization of the qPCR method for telomere length estimation, and discuss some of the factors that may cause variation in qPCR experiments.

How many genetic markers to tag an individual? An empirical assessment of false matching rates among close relatives.

Genetic identification of individuals is now commonplace, enabling the application of tagging methods to elusive species or species that cannot be tagged by traditional methods. A key aspect is determining the number of loci required to ensure that different individuals have non-matching multi-locus genotypes. Closely related individuals are of particular concern because of elevated matching probabilities caused by their recent co-ancestry. This issue may be addressed by increasing the number of loci to a level where full siblings (the relatedness category with the highest matching probability) are expected to have non-matching multi-locus genotypes. However, increasing the number of loci to meet this "full-sib criterion" greatly increases the laboratory effort, which in turn may increase the genotyping error rate resulting in an upward-biased mark-recapture estimate of abundance as recaptures are missed due to genotyping errors. We assessed the contribution of false matches from close relatives among 425 maternally related humpback whales, each genotyped at 20 microsatellite loci. We observed a very low (0.5-4%) contribution to falsely matching samples from pairs of first-order relatives (i.e., parent and offspring or full siblings). The main contribution to falsely matching individuals from close relatives originated from second-order relatives (e.g., half siblings), which was estimated at 9%. In our study, the total number of observed matches agreed well with expectations based upon the matching probability estimated for unrelated individuals, suggesting that the full-sib criterion is overly conservative, and would have required a 280% relative increase in effort. We suggest that, under most circumstances, the overall contribution to falsely matching samples from close relatives is likely to be low, and hence applying the full-sib criterion is unnecessary. In those cases where close relatives may present a significant issue, such as unrepresentative sampling, we propose three different genotyping strategies requiring only a modest increase in effort, which will greatly reduce the number of false matches due to the presence of related individuals.

Individual and population dietary specialization decline in fin whales during a period of ecosystem shift.

This study sought to estimate the effect of an anthropogenic and climate-driven change in prey availability on the degree of individual and population specialization of a large marine predator, the fin whale (Balaenoptera physalus). We examined skin biopsies from 99 fin whales sampled in the St. Lawrence Estuary (Canada) over a nine year period (1998-2006) during which environmental change was documented. We analyzed stable isotope ratios in skin and fatty acid signatures in blubber samples of whales, as well as in seven potential prey species, and diet was quantitatively assessed using Bayesian isotopic models. An abrupt change in fin whale dietary niche coincided with a decrease in biomass of their predominant prey, Arctic krill (Thysanoessa spp.). This dietary niche widening toward generalist diets occurred in nearly 60% of sampled individuals. The fin whale population, typically composed of specialists of either krill or lipid-rich pelagic fishes, shifted toward one composed either of krill specialists or true generalists feeding on various zooplankton and fish prey. This change likely reduced intraspecific competition. In the context of the current "Atlantification" of northern water masses, our findings emphasize the importance of considering individual-specific foraging tactics and not only population or group average responses when assessing population resilience or when implementing conservation measures.

Genetic analyses of historic and modern marbled murrelets suggest decoupling of migration and gene flow after habitat fragmentation.

The dispersal of individuals among fragmented populations is generally thought to prevent genetic and demographic isolation, and ultimately reduce extinction risk. In this study, we show that a century of reduction in coastal old-growth forests, as well as a number of other environmental factors, has probably resulted in the genetic divergence of marbled murrelets (Brachyramphus marmoratus) in central California, despite the fact that 7 per cent of modern-sampled murrelets in this population were classified as migrants using genetic assignment tests. Genetic differentiation appears to persist because individuals dispersing from northern populations contributed relatively few young to the central California population, as indicated by the fact that migrants were much less likely to be members of parent-offspring pairs than residents (10.5% versus 45.4%). Moreover, a recent 1.4 per cent annual increase in the proportion of migrants in central California, without appreciable reproduction, may have masked an underlying decline in the resident population without resulting in demographic rescue. Our results emphasize the need to understand the behaviour of migrants and the extent to which they contribute offspring in order to determine whether dispersal results in gene flow and prevents declines in resident populations.

Characterizing dispersal patterns in a threatened seabird with limited genetic structure.

Genetic assignment methods provide an appealing approach for characterizing dispersal patterns on ecological time scales, but require sufficient genetic differentiation to accurately identify migrants and a large enough sample size of migrants to, for example, compare dispersal between sexes or age classes. We demonstrate that assignment methods can be rigorously used to characterize dispersal patterns in a marbled murrelet (Brachyramphus marmoratus) population from central California that numbers approximately 600 individuals and is only moderately differentiated (F(ST) approximately 0.03) from larger populations to the north. We used coalescent simulations to select a significance level that resulted in a low and approximately equal expected number of type I and II errors and then used this significance level to identify a population of origin for 589 individuals genotyped at 13 microsatellite loci. The proportion of migrants in central California was greatest during winter when 83% of individuals were classified as migrants compared to lower proportions during the breeding (6%) and post-breeding (8%) seasons. Dispersal was also biased toward young and female individuals, as is typical in birds. Migrants were rarely members of parent-offspring pairs, suggesting that they contributed few young to the central California population. A greater number of migrants than expected under equilibrium conditions, a lack of individuals with mixed ancestry, and a small number of potential source populations (two), likely allowed us to use assignment methods to rigorously characterize dispersal patterns for a population that was larger and less differentiated than typically thought required for the identification of migrants.

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

The Gulf of California, Mexico is home to many cetacean species, including a presumed resident population of fin whales, Balaenoptera physalus. Past studies reported very low levels of genetic diversity among Gulf of California fin whales and a significant level of genetic differentiation from con-specifics in the eastern North Pacific. The aim of the present study was to assess the degree and timing of the isolation of Gulf of California fin whales in a population genetic analysis of 18 nuclear microsatellite genotypes from 402 samples and 565 mitochondrial control region DNA sequences (including mitochondrial sequences retrieved from NCBI). The analyses revealed that the Gulf of California fin whale population was founded ~2.3 thousand years ago and has since remained at a low effective population size (~360) and isolated from the eastern North Pacific (Nem between 0.89-1.4). The low effective population size and high degree of isolation implied that Gulf of California fin whales are vulnerable to the negative effects of genetic drift, human-caused mortality and habitat change.

Decline in abundance and apparent survival rates of fin whales (Balaenoptera physalus) in the northern Gulf of St. Lawrence.

Estimates of abundance and survivorship provide quantifiable measures to monitor populations and to define and understand their conservation status. This study investigated changes in abundance and survival rates of fin whales (Balaenoptera physalus) in the northern Gulf of St. Lawrence in the context of anthropogenic pressures and changing environmental conditions. A long-term data set, consisting of 35 years of photo-identification surveys and comprising more than 5,000 identifications of 507 individuals, formed the basis of this mark-recapture study. Based on model selection using corrected Akaike Information Criterion, the most parsimonious Cormack-Jolly-Seber model included a linear temporal trend in noncalf apparent survival rates with a sharp decline in the last 5 years of the study and a median survival rate of 0.946 (95% confidence interval (CI) 0.910-0.967). To account for capture heterogeneity due to divergent patterns of site fidelity, agglomerative hierarchical cluster analysis was employed to categorize individuals based on their annual and survey site fidelity indices. However, the negative trend in survivorship remained and was corroborated by a significant decline in the estimated super-population size from 335 (95% CI 321-348) individuals in 2004-2010 to 291 (95% CI 270-312) individuals in 2010-2016. Concurrently, a negative trend was estimated in recruitment to the population, supported by a sharp decrease in the number of observed calves. Ship strikes and changes in prey availability are potential drivers of the observed decline in fin whale abundance. The combination of clustering methods with mark-recapture represents a flexible way to investigate the effects of site fidelity on demographic variables and is broadly applicable to other individual-based studies.

Stable isotopes reveal winter feeding in different habitats in blue, fin and sei whales migrating through the Azores.

Knowing the migratory movements and behaviour of baleen whales is fundamental to understanding their ecology. We compared δ15N and δ13C values in the skin of blue (Balaenoptera musculus), fin (Balaenoptera physalus) and sei (Balaenoptera borealis) whales sighted in the Azores in spring with the values of potential prey from different regions within the North Atlantic using Bayesian mixing models to investigate their trophic ecology and migration patterns. Fin whale δ15N values were higher than those recorded in blue and sei whales, reflecting feeding at higher trophic levels. Whales' skin δ15N and δ13C values did not reflect prey from high-latitude summer foraging grounds; instead mixing models identified tropical or subtropical regions as the most likely feeding areas for all species during winter and spring. Yet, differences in δ13C values among whale species suggest use of different regions within this range. Blue and sei whales primarily used resources from the Northwest African upwelling and pelagic tropical/subtropical regions, while fin whales fed off Iberia. However, determining feeding habitats from stable isotope values remains difficult. In conclusion, winter feeding appears common among North Atlantic blue, fin and sei whales, and may play a crucial role in determining their winter distribution. A better understanding of winter feeding behaviour is therefore fundamental for the effective conservation of these species.