Population genomics of the white-beaked dolphin (Lagenorhynchus albirostris): Implications for conservation amid climate-driven range shifts.

Climate change is rapidly affecting species distributions across the globe, particularly in the North Atlantic. For highly mobile and elusive cetaceans, the genetic data needed to understand population dynamics are often scarce. Cold-water obligate species such as the white-beaked dolphin (Lagenorhynchus albirostris) face pressures from habitat shifts due to rising sea surface temperatures in addition to other direct anthropogenic threats. Unravelling the genetic connectivity between white-beaked dolphins across their range is needed to understand the extent to which climate change and anthropogenic pressures may impact species-wide genetic diversity and identify ways to protect remaining habitat. We address this by performing a population genomic assessment of white-beaked dolphins using samples from much of their contemporary range. We show that the species displays significant population structure across the North Atlantic at multiple scales. Analysis of contemporary migration rates suggests a remarkably high connectivity between populations in the western North Atlantic, Iceland and the Barents Sea, while two regional populations in the North Sea and adjacent UK and Irish waters are highly differentiated from all other clades. Our results have important implications for the conservation of white-beaked dolphins by providing guidance for the delineation of more appropriate management units and highlighting the risk that local extirpation may have on species-wide genetic diversity. In a broader context, this study highlights the importance of understanding genetic structure of all species threatened with climate change-driven range shifts to assess the risk of loss of species-wide genetic diversity.

Dietary plasticity and broad North Atlantic origins inferred from bulk and amino acid-specific δ15N and δ13C favour killer whale range expansions into Arctic waters.

Killer whales (Orcinus orca) occur seasonally in the eastern Canadian Arctic (ECA), where their range expansion associated with declining sea ice have raised questions about the impacts of increasing killer whale predation pressure on Arctic-endemic prey. We assessed diet and distribution of ECA killer whales using bulk and compound-specific stable isotope analysis (CSIA) of amino acids (AA) of 54 skin biopsies collected from 2009 to 2020 around Baffin Island, Canada. Bulk ECA killer whale skin δ15N and δ13C values did not overlap with potential Arctic prey after adjustment for trophic discrimination, and instead reflected foraging history in the North Atlantic prior to their arrival in the ECA. Adjusted killer whale stable isotope (SI) values primarily overlapped with several species of North Atlantic baleen whales or tuna. Amino acid (AA)-specific δ15N values indicated the ECA killer whales fed primarily on marine mammals, having similar glutamic acid δ15N-phenylalanine δ15N (δ15NGlx-Phe) and threonine δ15N (δ15NThr) as mammal-eating killer whales from the eastern North Pacific (ENP) that served as a comparative framework. However, one ECA whale grouped with the fish-eating ENP ecotype based δ15NThr. Distinctive essential AA δ13C of ECA killer whale groups, along with bulk SI similarity to killer whales from different regions of the North Atlantic, indicates different populations converge in Arctic waters from a broad source area. Generalist diet and long-distance dispersal capacity favour range expansions, and integration of these insights will be critical for assessing ecological impacts of increasing killer whale predation pressure on Arctic-endemic species.

Quantitative fatty acid signature analysis reveals a high level of dietary specialization in killer whales across the North Atlantic.

Quantifying the diet composition of apex marine predators such as killer whales (Orcinus orca) is critical to assessing their food web impacts. Yet, with few exceptions, the feeding ecology of these apex predators remains poorly understood. Here, we use our newly validated quantitative fatty acid signature analysis (QFASA) approach on nearly 200 killer whales and over 900 potential prey to model their diets across the 5000 km span of the North Atlantic. Diet estimates show that killer whales mainly consume other whales in the western North Atlantic (Canadian Arctic, Eastern Canada), seals in the mid-North Atlantic (Greenland), and fish in the eastern North Atlantic (Iceland, Faroe Islands, Norway). Nonetheless, diet estimates also varied widely among individuals within most regions. This level of inter-individual feeding variation should be considered for future ecological studies focusing on killer whales in the North Atlantic and other oceans. These estimates reveal remarkable population- and individual-level variation in the trophic ecology of these killer whales, which can help to assess how their predation impacts community and ecosystem dynamics in changing North Atlantic marine ecosystems. This new approach provides researchers with an invaluable tool to study the feeding ecology of oceanic top predators.

Connaître en détails la composition du régime alimentaire des grands prédateurs marins tels que les orques (Orcinus orca) est primordial afin d'évaluer leurs impacts sur les écosystèmes. Pourtant, à quelques exceptions près, l'écologie alimentaire de ces super-prédateurs reste mal comprise. Ici, nous utilisons notre nouvelle approche d'analyse quantitative des signatures d'acides gras (QFASA) sur près de 200 orques et plus de 900 proies potentielles pour modéliser leur régime alimentaire à travers l'Atlantique Nord. Les estimations de leurs régimes alimentaires montrent que les orques consomment principalement d'autres baleines dans l'ouest de l'Atlantique Nord (Arctique canadien, Est du Canada), des phoques dans le milieu de l'Atlantique Nord (Groenland) et des poissons dans l'est de l'Atlantique Nord (Islande, îles Féroé, Norvège). Néanmoins, ces estimations varient considérablement d'un individu à l'autre dans la plupart des régions. Cette variation alimentaire importante entre les individus doit être prise en compte dans les futures études écologiques qui s'intéressent aux orques de l'Atlantique Nord et d'ailleurs. Ces estimations révèlent des variations remarquables dans l'écologie trophique des orques tant au niveau des population que de l'individu, ce qui peut aider à évaluer l'impact de leur prédation sur la dynamique des communautés et des écosystèmes dans un contexte de changements climatiques en l'Atlantique Nord. Cette nouvelle approche fournit aux chercheurs un outil inestimable pour étudier l'écologie alimentaire des super-prédateurs océaniques.

Varying Diet Composition Causes Striking Differences in Legacy and Emerging Contaminant Concentrations in Killer Whales across the North Atlantic.

Lipophilic persistent organic pollutants (POPs) tend to biomagnify in food chains, resulting in higher concentrations in species such as killer whales (Orcinus orca) feeding on marine mammals compared to those consuming fish. Advancements in dietary studies include the use of quantitative fatty acid signature analysis (QFASA) and differentiation of feeding habits within and between populations of North Atlantic (NA) killer whales. This comprehensive study assessed the concentrations of legacy and emerging POPs in 162 killer whales from across the NA. We report significantly higher mean levels of polychlorinated biphenyls (PCBs), organochlorine pesticides, and flame retardants in Western NA killer whales compared to those of Eastern NA conspecifics. Mean ∑PCBs ranged from ∼100 mg/kg lipid weight (lw) in the Western NA (Canadian Arctic, Eastern Canada) to ∼50 mg/kg lw in the mid-NA (Greenland, Iceland) to ∼10 mg/kg lw in the Eastern NA (Norway, Faroe Islands). The observed variations in contaminant levels were strongly correlated with diet composition across locations (inferred from QFASA), emphasizing that diet and not environmental variation in contaminant concentrations among locations is crucial in assessing contaminant-associated health risks in killer whales. These findings highlight the urgency for implementing enhanced measures to safely dispose of POP-contaminated waste, prevent further environmental contamination, and mitigate the release of newer and potentially harmful contaminants.

Genomic insights into the historical and contemporary demographics of the grey reef shark.

Analyses of genetic diversity can shed light on both the origins of biodiversity hotspots, as well as the conservation status of species that are impacted by human activities. With these objectives, we assembled a genomic dataset of 14,935 single nucleotide polymorphisms from 513 grey reef sharks (Carcharhinus amblyrhynchos) sampled across 17 locations in the tropical Indo-Pacific. We analysed geographic variation in genetic diversity, estimated ancient and contemporary effective population size (Ne) across sampling locations (using coalescent and linkage disequilibrium methods) and modelled the history of gene flow between the Coral Triangle and the Coral Sea. Genetic diversity decreased with distance away from the Coral Triangle and north-western Australia, implying that C. amblyrhynchos may have originated in this region. Increases in Ne were detected across almost all sampling locations 40,000-90,000 generations ago (approximately 0.6-1.5 mya, given an estimated generation time of 16.4 years), suggesting a range expansion around this time. More recent, secondary increases in Ne were inferred for the Misool and North Great Barrier Reef sampling locations, but joint modelling did not clarify whether these were due to population growth, migration, or both. Despite the greater genetic diversity and ancient Ne observed at sites around Australia and the Coral Triangle, remote reefs around north-western New Caledonia had the highest contemporary Ne, demonstrating the importance of using multiple population size assessment methods. This study provides insight into both the past and present demographics of C. amblyrhynchos and contributes to our understanding of evolution in marine biodiversity hotspots.

Mercury and selenium concentrations, and selenium:mercury molar ratios in small cetaceans taken off St. Vincent, West Indies.

This study measured the concentration of total mercury (THg) and selenium (Se), and calculated the Se:Hg molar ratios in the muscle, blubber, liver, and kidney of small cetaceans (false killer whale, Pseudorca crassidens; killer whale, Orcinus orca; Risso's dolphin, Grampus griseus; short-finned pilot whale, Globicephala macrorhynchus; and dolphins of the genus Stenella) taken for human consumption off St. Vincent, West Indies. Overall, 122 samples were analyzed; mean THg concentrations (µg/g dry weight) were highest in the liver (730), followed by the kidney (274), muscle (76.4), and blubber (4.57). To explain variability in muscle THg concentrations, carbon (δ13C) and nitrogen (δ15N) stable isotope ratios were analyzed to explore differences in dietary carbon source and relative trophic position, respectively, among species. There was no relationship between δ15N and THg concentration, but there was a positive relationship between δ13C and THg concentration. On average for each species, the Se:Hg molar ratios were >1 in blubber and <1 in muscle. All liver samples and the majority of kidney, muscle, and blubber samples exceeded the FAO/WHO human consumption advisory level of 1 µg/g wet weight. Based on our estimations, consuming only 6.6 g of muscle a week would exceed the MeHg provisional tolerable weekly intake of 1.6 µg MeHg/kg body weight/week for a 60 kg person. Given the high THg concentration in these cetaceans and the frequency at which these tissues are consumed, this is a potential human health issue that warrants further investigation.

Inter-individual differences in ontogenetic trophic shifts among three marine predators.

Ontogenetic niche shifts are widespread. However, individual differences in size at birth, morphology, sex, and personalities can cause variability in behavior. As such, inherent inter-individual differences within populations may lead to context-dependent changes in behavior with animal body size, which is of concern for understanding population dynamics and optimizing ecological monitoring. Using stable carbon and nitrogen isotope values from concurrently sampled tissues, we quantified the direction and magnitude of intraspecific variation in trophic shifts among three shark species, and how these changed with body size: spurdogs (Squalus spp.) in deep-sea habitats off La Réunion, bull sharks (Carcharhinus leucas) in estuarine habitats of the Florida Everglades, and blacktip reef sharks (Carcharhinus melanopterus) in coral reef ecosystems of Moorea, French Polynesia. Intraspecific variation in trophic shifts was limited among spurdogs, and decreased with body size, while bull sharks exhibited greater individual differences in trophic shifts, but also decreased in variability through ontogeny. In contrast, blacktip reef sharks exhibited increased intraspecific variation in trophic interactions with body size. Variability in trophic interactions and ontogenetic shifts are known to be associated with changes in energetic requirements, but can vary with ecological context. Our results suggest that environmental stability may affect variability within populations, and ecosystems with greater spatial and/or temporal variability in environmental conditions, and those with more diverse food webs may facilitate greater individual differences in trophic interactions, and thus ontogenetic trophic shifts. In light of concerns over environmental disturbance, elucidating the contexts that promote or dampen phenotypic variability is invaluable for predicting population- and community-level responses to environmental changes.

Effect of body length, trophic position and habitat use on mercury concentrations of sharks from contrasted ecosystems in the southwestern Indian Ocean.

The non-essential metal mercury (Hg) can have deleterious effects on health of organisms, and tends to bioaccumulate with age in long-lived organisms and to biomagnify along food webs. Because elasmobranchs are fished for human consumption and their Hg levels are frequently above the maximum Hg concentration recommended for fish consumption, understanding the drivers of Hg concentration is of considerable interest. Total Hg concentrations were analysed in muscle tissues of 14 shark and 2 batoid species (n = 339 individuals) sampled across multiple habitats (coastal, open ocean and bathyal) in the southwestern Indian Ocean. Stable isotope ratios of carbon (δ13C) and nitrogen (δ15N) were analysed to assess whether relative trophic position and foraging habitats affected Hg concentrations. Hg concentrations increased with δ15N and body length, highlighting the mechanisms of bioaccumulation and biomagnification in relation with the trophic position and size of the individuals. Habitats where elasmobranchs were collected also affected their Hg concentrations. Bathyal sharks had high Hg concentrations that were almost similar to those of oceanic species, despite their lower relative trophic position. Higher bioavailability of Hg due to its enhanced methylation in deeper waters was considered as the most likely explanation for this result. These results highlight that multiple factors contribute to mercury accumulation in elasmobranchs.

Effects of anticoagulants on stable-isotope values (δ13 C and δ15 N) of shark blood components.

The effects of anticoagulant EDTA and sodium heparin (SH) on stable carbon δ13 C and nitrogen δ15 N isotopic values of red blood cells (RBC) and blood plasma in juvenile blacktip reef sharks Carcharhinus melanopterus were analysed. Plasma preserved with anticoagulants was not isotopically distinct from plasma stored in no-additive control tubes but RBC δ15 N values exhibited small enrichments when preserved with EDTA and SH. Results suggest EDTA and SH are viable anticoagulants for stable isotopic analyses of blood fractions but further studies are advised to validate results.

Using unmanned aerial vehicle (UAV) surveys and image analysis in the study of large surface-associated marine species: a case study on reef sharks Carcharhinus melanopterus shoaling behaviour.

A novel image analysis-based technique applied to unmanned aerial vehicle (UAV) survey data is described to detect and locate individual free-ranging sharks within aggregations. The method allows rapid collection of data and quantification of fine-scale swimming and collective patterns of sharks. We demonstrate the usefulness of this technique in a small-scale case study exploring the shoaling tendencies of blacktip reef sharks Carcharhinus melanopterus in a large lagoon within Moorea, French Polynesia. Using our approach, we found that C. melanopterus displayed increased alignment with shoal companions when distributed over a sandflat where they are regularly fed for ecotourism purposes as compared with when they shoaled in a deeper adjacent channel. Our case study highlights the potential of a relatively low-cost method that combines UAV survey data and image analysis to detect differences in shoaling patterns of free-ranging sharks in shallow habitats. This approach offers an alternative to current techniques commonly used in controlled settings that require time-consuming post-processing effort.

Contrasted accumulation patterns of persistent organic pollutants and mercury in sympatric tropical dolphins from the south-western Indian Ocean.

Due to their high trophic position and long life span, small cetaceans are considered as suitable bioindicators to monitor the presence of contaminants in marine ecosystems. Here, we document the contamination with persistent organic pollutants (POPs) and total mercury (T-Hg) of spinner (Stenella longirostris, n =21) and Indo-Pacific bottlenose dolphins (Tursiops aduncus, n=32) sampled from the coastal waters of La Réunion (south-western Indian Ocean). In addition, seven co-occurring teleost fish species were sampled and analyzed as well. Blubber samples from living dolphins and muscle from teleosts were analyzed for polychlorinated biphenyls (PCBs), DDT and metabolites (DDTs), chlordanes (CHLs), hexachlorocyclohexanes (HCHs), hexachlorobenzene (HCB), and polybrominated diphenyl ethers (PBDEs). Methoxylated PBDEs (MeO-PBDEs), reported as having a natural origin, were also analyzed. T-Hg levels were measured in blubber and skin biopsies of the two dolphin species. Stable isotopes δ(13)C and δ(15)N values were determined in skin of the dolphins and in the muscle of teleosts. For PCBs, HCHs and T-Hg, concentrations were significantly higher in T. aduncus than in S. longirostris. For other POP levels, intra-species variability was high. MeO-PBDEs were the dominant compounds (55% of the total POPs) in S. longirostris, while PCBs dominated (50% contribution) in T. aduncus. Other contaminants showed similar profiles between the two species. Given the different patterns of POPs and T-Hg contamination and the δ(15)N values observed among analyzed teleosts, dietary and foraging habitat preferences most likely explain the contrasted contaminant profiles observed in the two dolphin species. Levels of each class of contaminants were significantly higher in males than females. Despite their spatial and temporal overlap in the waters of La Réunion, S. longirostris and T. aduncus are differently exposed to contaminant accumulation.

Epidemiology of lobomycosis-like disease in bottlenose dolphins Tursiops spp. from South America and southern Africa.

We report on the epidemiology of lobomycosis-like disease (LLD), a cutaneous disorder evoking lobomycosis, in 658 common bottlenose dolphins Tursiops truncatus from South America and 94 Indo-Pacific bottlenose dolphins T. aduncus from southern Africa. Photographs and stranding records of 387 inshore residents, 60 inshore non-residents and 305 specimens of undetermined origin (inshore and offshore) were examined for the presence of LLD lesions from 2004 to 2015. Seventeen residents, 3 non-residents and 1 inshore dolphin of unknown residence status were positive. LLD lesions appeared as single or multiple, light grey to whitish nodules and plaques that may ulcerate and increase in size over time. Among resident dolphins, prevalence varied significantly among 4 communities, being low in Posorja (2.35%, n = 85), Ecuador, and high in Salinas, Ecuador (16.7%, n = 18), and Laguna, Brazil (14.3%, n = 42). LLD prevalence increased in 36 T. truncatus from Laguna from 5.6% in 2007-2009 to 13.9% in 2013-2014, albeit not significantly. The disease has persisted for years in dolphins from Mayotte, Laguna, Salinas, the Sanquianga National Park and Bahía Málaga (Colombia) but vanished from the Tramandaí Estuary and the Mampituba River (Brazil). The geographical range of LLD has expanded in Brazil, South Africa and Ecuador, in areas that have been regularly surveyed for 10 to 35 yr. Two of the 21 LLD-affected dolphins were found dead with extensive lesions in southern Brazil, and 2 others disappeared, and presumably died, in Ecuador. These observations stress the need for targeted epidemiological, histological and molecular studies of LLD in dolphins, especially in the Southern Hemisphere.

Legacy persistent organic pollutants among multiple cetacean species in the Northwest Atlantic.

The historical contamination of eastern Canadian shelf waters remains an ongoing concern, predominantly stemming from anthropogenic discharges in the Great Lakes region. Although legacy persistent organic pollutants (POPs) were banned decades ago, it remains unclear whether their concentrations have sufficiently decreased to safer levels in cetaceans that feed in the continental shelf waters of the northwestern Atlantic. This study compares polychlorinated biphenyl (PCB) and organochlorine pesticide (OC) accumulation in six cetacean species sampled in the Northwest Atlantic from 2015 to 2022. We assessed the influence of relative trophic level and foraging habitat preferences on POP accumulations among species using stable isotopes and fatty acids as dietary tracers. We further identified the species most susceptible to the effects of these contaminants. Killer whales (Orcinus orca) exhibited the highest PCB (∼100 mg/kg lw) and OC concentrations, followed by other odontocetes, with lowest concentrations in mysticetes. Stable isotope analysis revealed an unexpected lack of correlation between δ15N values and contaminant levels. However, there was a positive correlation between δ13C values and POP concentrations. Cetaceans foraging on pelagic prey species, as indicated by elevated proportions of the FA markers 22:1n11 and 20:1n9, had lower contaminant loads compared to cetaceans with benthic/coastal FA signatures. PCB and DDT (dichlorodiphenyltrichloroethane) concentrations are lower now in most cetacean species than in the 1980s and 1990s, likely due to regulatory measures restricting their production and use. Although current PCB concentrations for most species are under the thresholds for high risks of immune and reproductive failure, concentrations in killer whales exceed all established toxicity thresholds, underscoring the need for further action to reduce sources of these contaminants to the continental shelf waters of the northwestern Atlantic.

Agonistic interactions initiated by adult bottlenose dolphins on Antillean manatee calves in the Caribbean Sea.

The dynamics and drivers of inter-species interactions in the wild are poorly understood, particularly those involving social animal species. Inter-species interactions between cetaceans and sirenians have rarely been documented and investigated. Here, we report 10 cases of interaction initiated by adult bottlenose dolphins (Tursiops truncatus) towards Antillean manatee (Trichechus manatus manatus). Interactions were documented through behavioral observations in the wild (n = 7) and from the examination of orphaned calves (i.e., tooth rake marks on their body; n = 4) that entered a rehabilitation facility, one individual both observed interacting with dolphins and found stranded with bite marks. Bottlenose dolphins were observed interacting with orphan manatee calves and with mother-calf pairs, exhibiting agonistic behavior (n = 2), affiliative or neutral behaviors (n = 1), but the behavioral contexts of these interactions remain unclear in most cases (n = 7). Information on stranded individuals was collected from four calves (of 13 examined calves) recovered in poor condition with bottlenose dolphin tooth rakes and bite wounds on their bodies, one of which died. Injury from bite wounds varied in extent and severity, ranging from superficial scratches leaving rake marks to deep lacerations. Our findings suggest the regular occurrence of agonistic behaviors initiated by adult bottlenose dolphins and directed toward manatee calves. However, the drivers of these interactions remain unknown and need to be further investigated.

Spatial Variation in Mercury Accumulation in Bottlenose Dolphins (Tursiops spp.) in Southeastern U.S.A.

Bottlenose dolphins (Tursiops spp.) inhabit bays, sounds, and estuaries (BSEs) throughout the southeast region of the U.S.A. and are sentinel species for human and ecosystem-level health. Dolphins are vulnerable to the bioaccumulation of contaminants through the coastal food chain because they are high-level predators. Currently, there is limited information on the spatial dynamics of mercury accumulation in these dolphins. Total mercury (THg) was measured in dolphin skin from multiple populations across the U.S. Southeast Atlantic and Gulf of Mexico coasts, and the influence of geographic origin, sex, and age class was investigated. Mercury varied significantly among sampling sites and was greatest in dolphins in St. Joseph Bay, Florida Everglades, and Choctawhatchee Bay (14,193 ng/g ± 2196 ng/g, 10,916 ng/g ± 1532 ng/g, and 7333 ng/g ± 1405 ng/g wet mass (wm), respectively) and lowest in dolphins in Charleston and Skidaway River Estuary (509 ng/g ± 32.1 ng/g and 530 ng/g ± 58.4 ng/g wm, respectively). Spatial mercury patterns were consistent regardless of sex or age class. Bottlenose dolphin mercury exposure can effectively represent regional trends and reflect large-scale atmospheric mercury input and local biogeochemical processes. As a sentinel species, the bottlenose dolphin data presented here can direct future studies to evaluate mercury exposure to human residents in St. Joseph Bay, Choctawhatchee Bay, and Florida Coastal Everglades, as well as additional sites with similar geographical, oceanographic, or anthropogenic parameters. These data may also inform state and federal authorities that establish fish consumption advisories to determine if residents in these locales are at heightened risk for mercury toxicity.

Directed conservation of the world's reef sharks and rays.

Many shark populations are in decline around the world, with severe ecological and economic consequences. Fisheries management and marine protected areas (MPAs) have both been heralded as solutions. However, the effectiveness of MPAs alone is questionable, particularly for globally threatened sharks and rays ('elasmobranchs'), with little known about how fisheries management and MPAs interact to conserve these species. Here we use a dedicated global survey of coral reef elasmobranchs to assess 66 fully protected areas embedded within a range of fisheries management regimes across 36 countries. We show that conservation benefits were primarily for reef-associated sharks, which were twice as abundant in fully protected areas compared with areas open to fishing. Conservation benefits were greatest in large protected areas that incorporate distinct reefs. However, the same benefits were not evident for rays or wide-ranging sharks that are both economically and ecologically important while also threatened with extinction. We show that conservation benefits from fully protected areas are close to doubled when embedded within areas of effective fisheries management, highlighting the importance of a mixed management approach of both effective fisheries management and well-designed fully protected areas to conserve tropical elasmobranch assemblages globally.

Habitat partitioning, co-occurrence patterns, and mixed-species group formation in sympatric delphinids.

Numerous species have been reported to form mixed-species groups, however, little is known about the interplay between niche partitioning and mixed-species group formation. Furthermore, it is often unclear whether species come together by chance due to overlapping habitat preferences, by shared attraction to resources, or by attraction between them. We assessed habitat partitioning, co-occurrence patterns, and mixed-species group formation of sympatric Australian humpback (Sousa sahulensis) and Indo-Pacific bottlenose dolphins (Tursiops aduncus) around the North West Cape, Western Australia, with a joint species distribution model and temporal analyses of sighting data. Australian humpback dolphins preferred shallower and more nearshore waters than Indo-Pacific bottlenose dolphins, yet these species co-occurred more often than expected by chance given shared responses to environmental variables. Indo-Pacific bottlenose dolphins were sighted more often than Australian humpback dolphins during the afternoon, however, we did not find any temporal patterns in the occurrence of mixed-species groups. We propose that the positive association in the species' occurrence indicates the active formation of mixed-species groups. By evaluating habitat partitioning and co-occurrence patterns, this study provides direction for future work which should proceed to investigate the benefits that these species may gain from grouping with each other.

Critically endangered Rice's whales (Balaenoptera ricei) selectively feed on high-quality prey in the Gulf of Mexico.

Determining the drivers of prey selection in marine predators is critical when investigating ecosystem structure and function. The newly recognized Rice's whale (Balaenoptera ricei) is one of the most critically endangered large whales in the world and endemic to the industrialized Gulf of Mexico. Here, we investigated the drivers of resource selection by Rice's whales in relation to prey availability and energy density. Bayesian stable isotope (δ13C, δ15N) mixing models suggest that Rice's whales feed primarily on a schooling fish, Ariomma bondi (66.8% relative contribution). Prey selection using the Chesson's index revealed that active prey selection was found to be positive for three out of the four potential prey identified in the mixing model. A low degree of overlap between prey availability and diet inferred from the mixing model (Pianka Index: 0.333) suggests that prey abundance is not the primary driver of prey selection. Energy density data suggest that prey selection may be primarily driven by the energy content. Results from this study indicate that Rice's whales are selective predators consuming schooling prey with the highest energy content. Environmental changes in the region have the potential to influence prey species that would make them less available to Rice's whales.