Short-Finned Pilot Whale Strandings Associated with Pilot Whale Morbillivirus, Brazil.

Cetacean morbillivirus (CeMV) causes illness and death in cetaceans worldwide; the CeMV strains circulating in the Southern Hemisphere are poorly known. We detected a pilot whale CeMV strain in 3 short-finned pilot whales (Globicephala macrorhynchus) stranded in Brazil during July-October 2020. Our results confirm this virus circulates in this species.

Numerical-modeling-based investigation of sound transmission and reception in the short-finned pilot whale (Globicephala macrorhynchus).

The sound-transmission, beam-formation, and sound-reception processes of a short-finned pilot whale (Globicephala macrorhynchus) were investigated using computed tomography (CT) scanning and numerical simulation. The results showed that sound propagations in the forehead were modulated by the upper jaw, air components, and soft tissues, which attributed to the beam formation in the external acoustic field. These structures owned different acoustic impedance and formed a multiphasic sound transmission system that can modulate sounds into a beam. The reception pathways composed of the solid mandible and acoustic fats in the lower head conducted sounds into the tympano-periotic complex. In the simulations, sounds were emitted in the forehead transmission system and propagated into water to interrogate a steel cylinder. The resulting echoes can be interpreted from multiple perspectives, including amplitude, waveform, and spectrum, to obtain the acoustic cues of the steel cylinder. By taking the short-finned pilot whale as an example, this study provides meaningful information to further deepen our understanding of biosonar system operations, and may expand sound-reception theory in odontocetes.

Echolocation click parameters of short-finned pilot whales (Globicephala macrorhynchus) in the wild.

Short-finned pilot whales (Globicephala macrorhynchus) are large, deep-diving predators with diverse foraging strategies, but little is known about their echolocation. To quantify the source properties of short-finned pilot whale clicks, we made 15 deployments off the coast of Tenerife of a deep-water hydrophone array consisting of seven autonomous time-synced hydrophone recorders (SoundTraps), enabling acoustic localization and quantification of click source parameters. Of 8185 recorded pilot whale clicks, 47 were classified as being recorded on-axis, with a mean peak-to-peak source level (SL) of 181 ± 7 dB re 1 µPa, a centroid frequency of 40 ± 4 kHz, and a duration of 57 ± 23 µs. A fit to a piston model yielded an estimated half-power (-3 dB) beam width of 13.7° [95% confidence interval (CI) 13.2°-14.5°] and a mean directivity index (DI) of 22.6 dB (95% CI 22.5-22.9 dB). These measured SLs and DIs are surprisingly low for a deep-diving toothed whale, suggesting we sampled the short-finned pilot whales in a context with little need for operating a long-range biosonar. The substantial spectral overlap with beaked whale clicks emitted in similar deep-water habitats implies that pilot whale clicks may constitute a common source of false detections in beaked whale passive acoustic monitoring efforts.

Propensity for Risk in Reproductive Strategy Affects Susceptibility to Anthropogenic Disturbance.

AbstractAnimals initiate, interrupt, or invest resources in reproduction in light of their physiology and the environment. The energetic risks entailed in an individual's reproductive strategy can influence the ability to cope with additional stressors, such as anthropogenic climate change and disturbance. To explore the trade-offs between internal state, external resource availability, and reproduction, we applied state-dependent life-history theory (SDLHT) to a dynamic energy budget (DEB) model for long-finned pilot whales (Globicephala melas). We investigated the reproductive strategies emerging from the interplay between fitness maximization and propensity to take energetic risks, as well as the resulting susceptibility of individual vital rates to disturbance. Without disturbance, facultative reproductive behavior from SDLHT and fixed rules in the DEB model led to comparable individual fitness. However, under disturbance, the reproductive strategies emerging from SDLHT increased vulnerability to energetic risks, resulting in lower fitness than fixed rules. These fragile strategies might therefore be unlikely to evolve in the first place. Heterogeneous resource availability favored more cautious (and thus more robust) strategies, particularly when knowledge of resource variation was accurate. Our results demonstrate that the assumptions regarding the dynamic trade-offs underlying an individual's decision-making can have important consequences for predicting the effects of anthropogenic stressors on wildlife populations.

Diving apart together: call propagation in diving long-finned pilot whales.

Group-living animals must communicate to stay in contact. In long-finned pilot whales, there is a trade-off between the benefits of foraging individually at depth and the formation of tight social groups at the surface. Using theoretical modelling and empirical data of tagged pairs within a group, we examined the potential of pilot whale social calls to reach dispersed group members during foraging periods. Both theoretical predictions and empirical data of tag pairs showed a potential for communication between diving and non-diving group members over separation distances up to 385 m (empirical) and 1800 m (theoretical). These distances match or exceed pilot whale dive depths recorded across populations. Call characteristics and environmental characteristics were analysed to investigate determinants of call detectability. Longer calls with a higher sound pressure level (SPL) that were received in a quieter environment were more often detected than their shorter, lower SPL counterparts within a noisier environment. In a noisier environment, calls were louder and had a lower peak frequency, indicating mechanisms for coping with varying conditions. However, the vulnerability of pilot whales to anthropogenic noise is still of concern as the ability to cope with increasing background noise may be limited. Our study shows that combining propagation modelling and actual tag recordings provides new insights into the communicative potential for social calls in orientation and reunion with group members for deep-diving pilot whales.

Host age predicts parasite occurrence, richness, and nested infracommunities in a pilot whale-helminth network.

Ecological data on marine mammal parasites represent an excellent opportunity to expand our understanding of host-parasite systems. In this study, we used a dataset of intestinal helminth parasites on 167 long-finned pilot whales Globicephala melas (Traill, 1809) from seven localities in the Faroe Islands to evaluate the extent to which the host's age and sex influence the occurrence, richness, and nested pattern of helminth parasites and the importance of individual hosts to the helminth community. We found positive effects of age on both the occurrence and richness of helminths. Older host individuals showed an ordered accumulation of parasites, as evidenced by the nested pattern in their composition. Males had a higher occurrence of parasites than females, but the richness of helminths did not differ between sexes. Our findings suggest that differences in host-parasite interactions in long-finned pilot whales result mainly from age-structured variations in biological and behavioral characteristics.

Oceanographic barriers, divergence, and admixture: Phylogeography and taxonomy of two putative subspecies of short-finned pilot whale.

Genomic phylogeography plays an important role in describing evolutionary processes and their geographic, ecological, or cultural drivers. These drivers are often poorly understood in marine environments, which have fewer obvious barriers to mixing than terrestrial environments. Taxonomic uncertainty of some taxa (e.g., cetaceans), due to the difficulty in obtaining morphological data, can hamper our understanding of these processes. One such taxon, the short-finned pilot whale, is recognized as a single global species but includes at least two distinct morphological forms described from stranding and drive hunting in Japan, the "Naisa" and "Shiho" forms. Using samples (n = 735) collected throughout their global range, we examine phylogeographic patterns of divergence by comparing mitogenomes and nuclear SNP loci. Our results suggest three types within the species: an Atlantic Ocean type, a western/central Pacific and Indian Ocean (Naisa) type, and an eastern Pacific Ocean and northern Japan (Shiho) type. mtDNA control region differentiation indicates these three types form two subspecies, separated by the East Pacific Barrier: Shiho short-finned pilot whale, in the eastern Pacific Ocean and northern Japan, and Naisa short-finned pilot whale, throughout the remainder of the species' distribution. Our data further indicate two diverging populations within the Naisa subspecies, in the Atlantic Ocean and western/central Pacific and Indian Oceans, separated by the Benguela Barrier off South Africa. This study reveals a process of divergence and speciation within a globally-distributed, mobile marine predator, and indicates the importance of the East Pacific Barrier to this evolutionary process.

A whale of a tale: A One Environmental Health approach to study metal pollution in the Sea of Cortez.

Marine metal pollution is an emerging concern for human, animal, and ecosystem health. We considered metal pollution in the Sea of Cortez, which is a relatively isolated sea rich in biodiversity. Here there are potentially significant anthropogenic inputs of pollution from agriculture and metal mining. We considered the levels of 23 heavy metals and selenium in seven distinct cetacean species found in the area. Our efforts considered two different periods of time: 1999 and 2016/17. We considered the metal levels in relation to (1) all species together across years, (2) differences between suborders Odontoceti and Mysticeti, (3) each species individually across years, and (4) gender differences for each of these comparisons. We further compared metal levels found in sperm whale skin samples collected during these voyages to a previous voyage in 1999, to assess changes in metal levels over a longer timescale. The metals Mg, Fe, Al, and Zn were found at the highest concentrations across all species and all years. For sperm whales, we observed decreased metal levels from 1999 to 2016/2017, except for iron (Fe), nickel (Ni), and chromium (Cr), which either increased or did not change during this time period. These results indicate a recent change in the metal input to the Sea of Cortez, which may indicate a decreased concern for human, animal, and ecosystem health for some metals, but raises concern for the genotoxic metals Cr and Ni. This work was supported by NIEHS grant ES016893 (J.P.W.) and numerous donors to the Wise Laboratory.

Evidence for discrimination between feeding sounds of familiar fish and unfamiliar mammal-eating killer whale ecotypes by long-finned pilot whales.

Killer whales (KW) may be predators or competitors of other cetaceans. Since their foraging behavior and acoustics differ among populations ('ecotypes'), we hypothesized that other cetaceans can eavesdrop on KW sounds and adjust their behavior according to the KW ecotype. We performed playback experiments on long-finned pilot whales (Globicephala melas) in Norway using familiar fish-eating KW sounds (fKW) simulating a sympatric population that might compete for foraging areas, unfamiliar mammal-eating KW sounds (mKW) simulating a potential predator threat, and two control sounds. We assessed behavioral responses using animal-borne multi-sensor tags and surface visual observations. Pilot whales barely changed behavior to a broadband noise (CTRL-), whereas they were attracted and exhibited spyhops to fKW, mKW, and to a repeated-tonal upsweep signal (CTRL+). Whales never stopped nor started feeding in response to fKW, whereas they reduced or stopped foraging to mKW and CTRL+. Moreover, pilot whales joined other subgroups in response to fKW and CTRL+, whereas they tightened individual spacing within group and reduced time at surface in response to mKW. Typical active intimidation behavior displayed to fKW might be an antipredator strategy to a known low-risk ecotype or alternatively a way of securing the habitat exploited by a heterospecific sympatric population. Cessation of feeding and more cohesive approach to mKW playbacks might reflect an antipredator behavior towards an unknown KW ecotype of potentially higher risk. We conclude that pilot whales are able to acoustically discriminate between familiar and unfamiliar KW ecotypes, enabling them to adjust their behavior according to the perceived disturbance type.

Bio-energetic modeling of medium-sized cetaceans shows high sensitivity to disturbance in seasons of low resource supply.

Understanding the full scope of human impact on wildlife populations requires a framework to assess the population-level repercussions of nonlethal disturbance. The Population Consequences of Disturbance (PCoD) framework provides such an approach, by linking the effects of disturbance on the behavior and physiology of individuals to their population-level consequences. Bio-energetic models have been used as implementations of PCoD, as these integrate the behavioral and physiological state of an individual with the state of the environment, to mediate between disturbance and biological significant changes in vital rates (survival, growth, and reproduction). To assess which levels of disturbance lead to adverse effects on population growth rate requires a bio-energetic model that covers the complete life cycle of the organism under study. In a density-independent setting, the expected lifetime reproductive output of a single female can then be used to predict the level of disturbance that leads to population decline. Here, we present such a model for a medium-sized cetacean, the long-finned pilot whale (Globicephala melas). Disturbance is modeled as a yearly recurrent period of no resource feeding for the pilot whale female and her calf. Short periods of disturbance lead to the pre-weaned death of the first one or more calves of the young female. Higher disturbance levels also affect survival of calves produced later in the life of the female, in addition to degrading female survival. The level of disturbance that leads to a negative population growth rate strongly depends on the available resources in the environment. This has important repercussion for the timing of disturbance if resource availability fluctuates seasonally. The model predicts that pilot whales can tolerate on average three times longer periods of disturbance in seasons of high resource availability, compared to disturbance happening when resources are low. Although our model is specifically parameterized for pilot whales, it provides useful insights into the general consequences of nonlethal disturbance. If appropriate data on life history and energetics are available, it can be used to provide management advice for specific species or populations.

Selective reactions to different killer whale call categories in two delphinid species.

The risk of predation is often invoked as an important factor influencing the evolution of social organization in cetaceans, but little direct information is available about how these aquatic mammals respond to predators or other perceived threats. We used controlled playback experiments to examine the behavioral responses of short-finned pilot whales (Globicephala macrorhynchus) off Cape Hatteras, NC, USA, and Risso's dolphins (Grampus griseus) off the coast of Southern California, USA, to the calls of a potential predator, mammal-eating killer whales. We transmitted calls of mammal-eating killer whales, conspecifics and baleen whales to 10 pilot whales and four Risso's dolphins equipped with multi-sensor archival acoustic recording tags (DTAGs). Only playbacks of killer whale calls resulted in significant changes in tagged animal heading. The strong responses observed in both species occurred only following exposure to a subset of killer whale calls, all of which contained multiple non-linear properties. This finding suggests that these structural features of killer whale calls convey information about predatory risk to pilot whales and Risso's dolphins. The observed responses differed between the two species; pilot whales approached the sound source while Risso's dolphins fled following playbacks. These divergent responses likely reflect differences in anti-predator response mediated by the social structure of the two species.

Species-level classification of beaked whale echolocation signals detected in the northern Gulf of Mexico.

This study presents and evaluates several methods for automated species-level classification of echolocation clicks from three beaked whale species recorded in the northern Gulf of Mexico. The species included are Cuvier's and Gervais' beaked whales, as well as an unknown species denoted Beaked Whale Gulf. An optimal feature set for discriminating the three click types while also separating detected clicks from unidentified delphinids was determined using supervised step-wise discriminant analysis. Linear and quadratic discriminant analyses both achieved error rates below 1% with three features, determined by tenfold cross validation. The waveform fractal dimension was found to be a highly ranked feature among standard spectral and temporal parameters. The top-ranking features were Higuchi's fractal dimension, spectral centroid, Katz's fractal dimension, and -10 dB duration. Six clustering routines, including four popular network-based algorithms, were also evaluated as unsupervised classification methods using the selected feature set. False positive rates of 0.001 and 0.024 were achieved by Chinese Whispers and spectral clustering, respectively, across 200 randomized trials. However, Chinese Whispers clustering yielded larger false negative rates. Spectral clustering was further tested on clicks from encounters of beaked, sperm, and pilot whales in the Tongue of the Ocean, Bahamas.

Familial social structure and socially driven genetic differentiation in Hawaiian short-finned pilot whales.

Social structure can have a significant impact on divergence and evolution within species, especially in the marine environment, which has few environmental boundaries to dispersal. On the other hand, genetic structure can affect social structure in many species, through an individual preference towards associating with relatives. One social species, the short-finned pilot whale (Globicephala macrorhynchus), has been shown to live in stable social groups for periods of at least a decade. Using mitochondrial control sequences from 242 individuals and single nucleotide polymorphisms from 106 individuals, we examine population structure among geographic and social groups of short-finned pilot whales in the Hawaiian Islands, and test for links between social and genetic structure. Our results show that there are at least two geographic populations in the Hawaiian Islands: a Main Hawaiian Islands (MHI) population and a Northwestern Hawaiian Islands/Pelagic population (FST and ΦST p < .001), as well as an eastern MHI community and a western MHI community (FST p = .009). We find genetically driven social structure, or high relatedness among social units and clusters (p < .001), and a positive relationship between relatedness and association between individuals (p < .0001). Further, socially organized clusters are genetically distinct, indicating that social structure drives genetic divergence within the population, likely through restricted mate selection (FST p = .05). This genetic divergence among social groups can make the species less resilient to anthropogenic or ecological disturbance. Conservation of this species therefore depends on understanding links among social structure, genetic structure and ecological variability within the species.

High diving metabolic rate indicated by high-speed transit to depth in negatively buoyant long-finned pilot whales.

To maximize foraging duration at depth, diving mammals are expected to use the lowest cost optimal speed during descent and ascent transit and to minimize the cost of transport by achieving neutral buoyancy. Here, we outfitted 18 deep-diving long-finned pilot whales with multi-sensor data loggers and found indications that their diving strategy is associated with higher costs than those of other deep-diving toothed whales. Theoretical models predict that optimal speed is proportional to (basal metabolic rate/drag)1/3 and therefore to body mass0.05 The transit speed of tagged animals (2.7±0.3 m s-1) was substantially higher than the optimal speed predicted from body mass (1.4-1.7 m s-1). According to the theoretical models, this choice of high transit speed, given a similar drag coefficient (median, 0.0035) to that in other cetaceans, indicated greater basal metabolic costs during diving than for other cetaceans. This could explain the comparatively short duration (8.9±1.5 min) of their deep dives (maximum depth, 444±85 m). Hydrodynamic gliding models indicated negative buoyancy of tissue body density (1038.8±1.6 kg m-3, ±95% credible interval, CI) and similar diving gas volume (34.6±0.6 ml kg-1, ±95% CI) to those in other deep-diving toothed whales. High diving metabolic rate and costly negative buoyancy imply a 'spend more, gain more' strategy of long-finned pilot whales, differing from that in other deep-diving toothed whales, which limits the costs of locomotion during foraging. We also found that net buoyancy affected the optimal speed: high transit speeds gradually decreased during ascent as the whales approached neutral buoyancy owing to gas expansion.

Temporal Shifts in Poly- and Perfluoroalkyl Substances (PFASs) in North Atlantic Pilot Whales Indicate Large Contribution of Atmospheric Precursors.

Poly- and perfluoroalkyl substances (PFASs) are persistent, bioaccumulative anthropogenic compounds associated with adverse health impacts on humans and wildlife. PFAS production changed in North America and Europe around the year 2000, but impacts on wildlife appear to vary across species and location. Unlike other mammal species, cetaceans lack the enzyme for transforming an important intermediate precursor (perfluorooctane sulfonamide: FOSA), into a prevalent compound in most wildlife (perfluorooctanesulfonate: PFOS). Thus, their tissue burden differentiates these two compounds while other mammals contain PFOS from both direct exposure and precursor degradation. Here we report temporal trends in 15 PFASs measured in muscle from juvenile male North Atlantic pilot whales (Globicephala melas) harvested between 1986 and 2013. FOSA accounted for a peak of 84% of the 15 PFASs around 2000 but declined to 34% in recent years. PFOS and long-chained PFCAs (C9-C13) increased significantly over the whole period (2.8% yr-1 to 8.3% yr-1), but FOSA declined by 13% yr-1 after 2006. Results from FOSA partitioning and bioaccumulation modeling forced by changes in atmospheric inputs reasonably capture magnitudes and temporal patterns in FOSA concentrations measured in pilot whales. Rapid changes in atmospheric FOSA in polar and subpolar regions around 2000 helps to explain large declines in PFOS exposure for species that metabolize FOSA, including seafood consuming human populations. This work reinforces the importance of accounting for biological exposures to PFAS precursors.

Steroid hormones and persistent organic pollutants in plasma from North-eastern Atlantic pilot whales.

Persistent organic pollutants (POPs) are known to have endocrine disruptive effects, interfering with endogenous steroid hormones. The present study examined nine steroid hormones and their relationships with the concentrations of selected POPs in pilot whales (Globicephala melas) from the Faroe Islands, NE Atlantic. The different steroids were detected in 15 to all of the 26 individuals. High concentrations of progesterone (83.3-211.7pmol/g) and pregnenolone (PRE; 4.68-5.69pmol/g) were found in three adult females indicating that they were pregnant or ovulating. High androgen concentrations in two of the males reflected that one was adult and that one (possibly) had reached puberty. In males a significant positive and strong correlation between body length and testosterone (TS) levels was identified. Furthermore, positive and significant correlations were found between 4-OH-CB107/4'-OH-CB108 and 17ß-estradiol in males. In adult females significant positive correlations were identified between PRE and CB149 and t-nonachlor, between estrone and CB138, -149, -187 and p,p'-DDE, between androstenedione and CB187, and between TS and CB-99 and -153. Although relationships between the POPs and the steroid hormones reported herein are not evidence of cause-effect relationships, the positive correlations between steroids and POPs, particularly in females, suggest that POPs may have some endocrine disrupting effects on the steroid homeostasis in this species.

Vocal repertoire of long-finned pilot whales (Globicephala melas) in northern Norway.

The knowledge of the vocal repertoire of pilot whales is very limited. In this paper, the vocal repertoire of long-finned pilot whales recorded during different encounters in the Vestfjord in northern Norway between November 2006 and August 2010 are described. Sounds were analysed using two different methods: (1) an observer-based audio-visual inspection of FFT-derived spectrograms, with which, besides a general variety of clicks, buzzes, nonharmonic sounds, and whistles, 129 different distinct call types and 25 subtypes were distinguished. These call types included pulsed calls and discrete structured whistles varying from simple to highly complex structures composed of several segments and elements. In addition, ultrasonic whistles previously not described for pilot whales were found. In addition to the diversity of single calls, call sequences consisting of repetitions and combinations of specific call types were recorded and (2) a parametric approach that permitted the confirmation of the high variability in pilot whale call structures was developed. It is concluded that the pilot whale vocal repertoire is among the most complex for the mammalian species and the high structural variability, along with call repetitions and combinations, require a closer investigation to judge their importance for vocal communication.

Acoustic properties of a short-finned pilot whale head with insight into temperature influence on tissues' sound velocity.

Acoustic properties of odontocete head tissues, including sound velocity, density, and acoustic impedance, are important parameters to understand dynamics of its echolocation. In this paper, acoustic properties of head tissues from a freshly dead short-finned pilot whale (Globicephala macrorhynchus) were reconstructed using computed tomography (CT) and ultrasound. The animal's forehead soft tissues were cut into 188 ordered samples. Sound velocity, density, and acoustic impedance of each sample were either directly measured or calculated by formula, and Hounsfield Unit values (HUs) were obtained from CT scanning. According to relationships between HUs and sound velocity, HUs and density, as well as HUs and acoustic impedance, distributions of acoustic properties in the head were reconstructed. The inner core in the melon with low-sound velocity and low-density is an evidence for its potential function of sound focusing. The increase in acoustic impedance of forehead tissues from inner core to outer layer may be important for the acoustic impedance matching between the outer layer tissue and seawater. In addition, temperature dependence of sound velocity in soft tissues was also examined. The results provide a guide to the simulation of the sound emission of the short-finned pilot whale.

Comparative Chromosome Map and Heterochromatin Features of the Gray Whale Karyotype (Cetacea).

Cetacean karyotypes possess exceptionally stable diploid numbers and highly conserved chromosomes. To date, only toothed whales (Odontoceti) have been analyzed by comparative chromosome painting. Here, we studied the karyotype of a representative of baleen whales, the gray whale (Eschrichtius robustus, Mysticeti), by Zoo-FISH with dromedary camel and human chromosome-specific probes. We confirmed a high degree of karyotype conservation and found an identical order of syntenic segments in both branches of cetaceans. Yet, whale chromosomes harbor variable heterochromatic regions constituting up to a third of the genome due to the presence of several types of repeats. To investigate the cause of this variability, several classes of repeated DNA sequences were mapped onto chromosomes of whale species from both Mysticeti and Odontoceti. We uncovered extensive intrapopulation variability in the size of heterochromatic blocks present in homologous chromosomes among 3 individuals of the gray whale by 2-step differential chromosome staining. We show that some of the heteromorphisms observed in the gray whale karyotype are due to distinct amplification of a complex of common cetacean repeat and heavy satellite repeat on homologous autosomes. Furthermore, we demonstrate localization of the telomeric repeat in the heterochromatin of both gray and pilot whale (Globicephala melas, Odontoceti). Heterochromatic blocks in the pilot whale represent a composite of telomeric and common repeats, while heavy satellite repeat is lacking in the toothed whale consistent with previous studies.

Relationships between concentrations of selected organohalogen contaminants and thyroid hormones and vitamins A, E and D in Faroese pilot whales.

Pilot whales (Globicephala melas) from the Faroe Islands, North-East Atlantic, have high body concentrations of organohalogenated compounds (OHCs), such as polychlorinated biphenyls (PCBs), organochlorinated pesticides (OCPs) and brominated flame retardants (BFRs). The aim of the present study was to examine if and to what extent blood plasma and liver concentrations of several groups of these OHCs are related to concentrations of relevant nutritional and hormonal biomarkers in pilot whales. Thyroid hormones (THs: total and free thyroxine and total and free triiodothyronine) and vitamin A (retinol), D (25-hydroxyvitamin D3) and E (α-tocopherol) were analysed in plasma (n=27) and vitamin A (total vitamin A, retinol and retinyl palmitate) and E (α- and γ-tocopherol) were analysed in liver (n=37) of Faroe Island pilot whales. Correlative relationships between the biomarkers and OHC concentrations previously analysed in the same tissues in these individuals were studied. The TH concentrations in plasma were significantly higher in juveniles than in adults. Vitamin D concentrations in plasma and α- and γ-tocopherol in liver were higher in adults than in juveniles. Multivariate statistical modelling showed that the age and sex influenced the relationship between biomarkers and OHCs. Some significant positive relationships were found between OHCs and thyroid hormone concentrations in the youngest juveniles (p<0.05). In plasma of juvenile whales α-tocopherol was also positively correlated with all the OHCs (p<0.05). Only few significant correlations were found between single OHCs and retinol and vitamin D in plasma within the age groups. There were significant negative relationships between hepatic PBDE concentrations and retinol (BDE-47) and γ-tocopherol (BDE-49, -47, -100, -99, -153) in liver. The relationships between OHCs and THs or vitamins suggest that in pilot whales OHCs seem to have minor effects on TH and vitamin concentrations.