A review of the biology, fisheries and conservation of the whale shark Rhincodon typus

Persistent Organic Pollutants in Whale Shark (Rhincodon typus) Skin Biopsies from Bahía de Los Ángeles, Mexico

The whale shark (Rhincodon typus) is a filter-feeding organism that can be considered a sentinel species, and Bahía de los Ángeles (BLA) in the Gulf of California is an important sighting site for these elasmobranchs. This filter-feeding organism can be considered a pollutant sampler from the marine environment. Persistent organic pollutants are toxic compounds with high mobility and environmental persistence, bioaccumulation and trophic transfer. Among these are polycyclic aromatic hydrocarbons (PAHs) and organochlorine pesticides (OCPs). The present work aimed to determine concentrations of PAHs and OCPs in whale shark skin biopsies, collected in 2021 at BLA. Mean detected levels of PAHs and OCPs were 279.4 ng/g dw (dry weight) and 1478.1 ng/g dw, respectively. Analysis of similarities between the ordered sizes (4.2-7.6 m) and the concentrations of PAHs and OCPs indicated no significant differences. Individual PAHs detected indicate pyrogenic and petrogenic sources; the presence of pesticides at levels higher than those of hydrocarbons may be related to agricultural activity in the areas surrounding the Baja California peninsula. This study is the first report of PAH levels in R. typus for the Gulf of California and Mexico.

The biology and ecology of the basking shark: A review

Here we review the literature on the basking shark (Cetorhinus maximus, Gunnerus, 1765), well known as the second largest extant shark (and fish) species globally. Previous reviews were published by Kunzlik in 1988 and Sims in 2008, but in the last 15 years modern electronic and DNA sequencing technologies have resulted in considerable advances in our knowledge of the species' behaviour and ecology. Basking sharks are planktivores and under appropriate conditions spend prolonged periods at the ocean surface feeding on copepod prey that primarily make up their diet, the behaviour that gave rise to their common name. In general, they are migratory and move into higher latitude waters during the summer months, when loose surface-feeding aggregations may form at favoured sites, the best known of which at present occur at hotspots on the west coasts of Britain and Ireland. The species is found circumglobally in temperate waters, but they are also now known on occasion to migrate at depth between northern and southern hemispheres, as well as across oceans within the northern hemisphere. In the past basking shark were more abundant across much of their range, but, consequent on targeted fisheries and in some places intentional eradication, became everywhere scarce, with recent population recovery in the north-east Atlantic being the result of protective measures initiated in the 1990s. Despite their charismatic nature, some of their most fundamental biological processes including copulation, gestation and birth remain largely unknown, due to their migratory and often deep-water lifestyle. In contrast, the deployment of small-scale archival and satellite tags has revealed the details of both broadscale migratory movements and horizontal and vertical foraging behaviours. Recent genetic studies support evidence suggesting a degree of site fidelity in relation to seasonal feeding grounds, which likely explains why in the past local populations have collapsed following periods of intensive fishing. Other recent research using aerial drones and towed cameras has revealed within loose feeding aggregations elements of social behaviour that may have a courtship function as well as enhance feeding efficiency.

Minke whale feeding rate limitations suggest constraints on the minimum body size for engulfment filtration feeding

Bulk filter feeding has enabled gigantism throughout evolutionary history. The largest animals, extant rorqual whales, utilize intermittent engulfment filtration feeding (lunge feeding), which increases in efficiency with body size, enabling their gigantism. The smallest extant rorquals (7-10 m minke whales), however, still exhibit short-term foraging efficiencies several times greater than smaller non-filter-feeding cetaceans, raising the question of why smaller animals do not utilize this foraging modality. We collected 437 h of bio-logging data from 23 Antarctic minke whales (Balaenoptera bonaerensis) to test the relationship of feeding rates (λf) to body size. Here, we show that while ultra-high nighttime λf (mean ± s.d.: 165 ± 40 lunges h-1; max: 236 lunges h-1; mean depth: 28 ± 46 m) were indistinguishable from predictions from observations of larger species, daytime λf (mean depth: 72 ± 72 m) were only 25-40% of predicted rates. Both λf were near the maxima allowed by calculated biomechanical, physiological and environmental constraints, but these temporal constraints meant that maximum λf was below the expected λf for animals smaller than ~5 m-the length of weaned minke whales. Our findings suggest that minimum size for specific filter-feeding body plans may relate broadly to temporal restrictions on filtration rate and have implications for the evolution of filter feeding.

Species associated with whale sharks Rhincodontypus (Orectolobiformes, Rhincodontidae) in the Galapagos Archipelago

Whale sharks Rhincodontypus frequently appear to interact or associate with other species, which vary depending on the community structure and the demographic of the whale sharks at each location globally. Here, we present the species sighted frequently around whale sharks in the Galapagos Archipelago and reported by dive guides and scientists and also in earlier publications. These associated species include cetacean species: bottlenose dolphins Tursiopstruncatus, other shark species: silky sharks Carcharhinusfalciformis, Galapagos sharks Carcharhinusgalapagensis, scalloped hammerhead sharks Sphyrnalewini, tiger sharks Galeocerdocuvier and teleost fish species: remoras Remora remora, yellowfin tuna Thunnusalbacares, almaco jacks Seriolarivoliana and black jacks Caranxlugubris. The recording of interspecies associations and interactions may lead to better understanding of the natural history of whale sharks and can show important symbiotic relationships or interdependence between different species.

To what extent are filter feeder elasmobranchs exposed to marine pollution? A systematic review

Filter feeding elasmobranchs may be considered as biological indicators of marine pollution, despite most of these species are under some degree of extinction risk. Among threats to this taxonomic group, marine pollution might represent an additional concern for their survival. In this review, a comprehensive systematic search of scientific literature on pollutants in filter feeding elasmobranchs was conducted to evaluate the bioaccumulation patterns, and risk for human consumers. We found that, despite an increasing trend in the number of published studies, the geographical coverage is still very limited and most of the studies focused solely on trace elements (70.8%). Among sharks, Rhincodon typus was the most represented species (66.7%), while Mobula mobular the most studied ray species (41.7%). Comparing the levels of pollutants in filter feeders between ocean basins, this review highlighted that Hg, As and Cd levels are mostly higher in those areas affected by both strong natural and anthropogenic source of emissions, such as the Indian Ocean. With regards to OCs, ΣPCB levels in muscle of C. maximus were between 4.3 and 50.5 μg kg^-1 ww, highlighting a persistent contamination of PCB in the Mediterranean Sea. Some species exceeded the maximum allowable limits for foodstuff consumption for As, Cd and Pb. A total of 77.8% of the analyzed species exceeded the Environmental Quality Standards for Hg, while they were always below the EQS_biota for HCB, PBDEs, PFOS and DDT. Given their feeding mechanism that continuously samples the marine environment, further investigations are urgently needed to determine not only the extent of contaminant exposure in different hotspot locations but also the risks posed to the elasmobranch health.

Global collision-risk hotspots of marine traffic and the world's largest fish, the whale shark

Marine traffic is increasing globally yet collisions with endangered megafauna such as whales, sea turtles, and planktivorous sharks go largely undetected or unreported. Collisions leading to mortality can have population-level consequences for endangered species. Hence, identifying simultaneous space use of megafauna and shipping throughout ranges may reveal as-yet-unknown spatial targets requiring conservation. However, global studies tracking megafauna and shipping occurrences are lacking. Here we combine satellite-tracked movements of the whale shark, Rhincodon typus, and vessel activity to show that 92% of sharks’ horizontal space use and nearly 50% of vertical space use overlap with persistent large vessel (>300 gross tons) traffic. Collision-risk estimates correlated with reported whale shark mortality from ship strikes, indicating higher mortality in areas with greatest overlap. Hotspots of potential collision risk were evident in all major oceans, predominantly from overlap with cargo and tanker vessels, and were concentrated in gulf regions, where dense traffic co-occurred with seasonal shark movements. Nearly a third of whale shark hotspots overlapped with the highest collision-risk areas, with the last known locations of tracked sharks coinciding with busier shipping routes more often than expected. Depth-recording tags provided evidence for sinking, likely dead, whale sharks, suggesting substantial “cryptic” lethal ship strikes are possible, which could explain why whale shark population declines continue despite international protection and low fishing-induced mortality. Mitigation measures to reduce ship-strike risk should be considered to conserve this species and other ocean giants that are likely experiencing similar impacts from growing global vessel traffic.

Microbiome structure in large pelagic sharks with distinct feeding ecologies

Background

Sharks play essential roles in ocean food webs and human culture, but also face population declines worldwide due to human activity. The relationship between sharks and the microbes on and in the shark body is unclear, despite research on other animals showing the microbiome as intertwined with host physiology, immunity, and ecology. Research on shark-microbe interactions faces the significant challenge of sampling the largest and most elusive shark species. We leveraged a unique sampling infrastructure to compare the microbiomes of two apex predators, the white (Carcharodon carcharias) and tiger shark (Galeocerdo cuvier), to those of the filter-feeding whale shark (Rhincodon typus), allowing us to explore the effects of feeding mode on intestinal microbiome diversity and metabolic function, and environmental exposure on the diversity of microbes external to the body (on the skin, gill).

Results

The fecal microbiomes of white and whale sharks were highly similar in taxonomic and gene category composition despite differences in host feeding mode and diet. Fecal microbiomes from these species were also taxon-poor compared to those of many other vertebrates and were more similar to those of predatory teleost fishes and toothed whales than to those of filter-feeding baleen whales. In contrast, microbiomes of external body niches were taxon-rich and significantly influenced by diversity in the water column microbiome.

Conclusions

These results suggest complex roles for host identity, diet, and environmental exposure in structuring the shark microbiome and identify a small, but conserved, number of intestinal microbial taxa as potential contributors to shark physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00168-x.

Assessing the exposure risk of large pelagic fish to oil spills scenarios in the deep waters of the Gulf of Mexico

Exposure risk is assessed based on modeling suitable habitat of large pelagic fish and oil spill scenarios originating at three wells located in the western GM's deep waters. Since the fate of the oil depends on the oceanographic conditions present during the accident, as well as the magnitude and duration of the spill, which are not known a priori, the scenarios used are a statistical representation of the area in which oil spilled from the well could be found, given all possible outcomes. The ecological vulnerability assessment identified a subset of bony fish with low-medium vulnerability and elasmobranchs with medium-high vulnerability. The oiling probability and exposure risk of both bony fish and elasmobranchs hotspots vary by well analyzed. Thus, these results provide essential information for a risk management plan for the assessed species and others with economic or conservation importance distributed in the GM and worldwide.

Pieces in a global puzzle: Population genetics at two whale shark aggregations in the western Indian Ocean

The whale shark Rhincodon typus is found throughout the world's tropical and warm-temperate ocean basins. Despite their broad physical distribution, research on the species has been concentrated at a few aggregation sites. Comparing DNA sequences from sharks at different sites can provide a demographically neutral understanding of the whale shark's global ecology. Here, we created genetic profiles for 84 whale sharks from the Saudi Arabian Red Sea and 72 individuals from the coast of Tanzania using a combination of microsatellite and mitochondrial sequences. These two sites, separated by approximately 4500 km (shortest over-water distance), exhibit markedly different population demographics and behavioral ecologies. Eleven microsatellite DNA markers revealed that the two aggregation sites have similar levels of allelic richness and appear to be derived from the same source population. We sequenced the mitochondrial control region to produce multiple global haplotype networks (based on different alignment methodologies) that were broadly similar to each other in terms of population structure but suggested different demographic histories. Data from both microsatellite and mitochondrial markers demonstrated the stability of genetic diversity within the Saudi Arabian aggregation site throughout the sampling period. These results contrast previously measured declines in diversity at Ningaloo Reef, Western Australia. Mapping the geographic distribution of whale shark lineages provides insight into the species' connectivity and can be used to direct management efforts at both local and global scales. Similarly, understanding historical fluctuations in whale shark abundance provides a baseline by which to assess current trends. Continued development of new sequencing methods and the incorporation of genomic data could lead to considerable advances in the scientific understanding of whale shark population ecology and corresponding improvements to conservation policy.

Global shifts in species richness have shaped carpet shark evolution

Background

The evolutionary processes that shape patterns of species richness in marine ecosystems are complex and may differ between organismal groups. There has been considerable interest in understanding the evolutionary processes that led to marine species richness being concentrated in specific geographical locations. In this study we focus on the evolutionary history of a group of small-to-medium sized sharks known as carpet sharks. While a few carpet shark species are widespread, the majority of carpet shark species richness is contained within a biodiversity hotspot at the boundary of the Indian and Pacific oceans. We address the significance of this biodiversity hotspot in carpet shark evolution and speciation by leveraging a rich fossil record and molecular phylogenetics to examine the prehistoric distribution of carpet sharks.

Results

We find that carpet sharks species richness was greatest in shallow seas connected to the Atlantic Ocean during the Late Cretaceous, but that there was a subsequent loss of biodiversity in Atlantic waters. Fossil evidence from sites in close geographic proximity to the current center of carpet shark diversity are generally restricted to younger geologic strata.

Conclusions

From this data we conclude that (1) center of carpet shark biodiversity has shifted during the last 100 million years, (2) carpet sharks have repeatedly dispersed to nascent habitat (including to their current center of diversity), and (3) the current center of carpet shark biodiversity conserves lineages that have been extirpated from this prehistoric range and is a source of new carpet shark species. Our findings provide insights into the roles of marine biodiversity hotspots for higher-tropic level predators and the methods applied here can be used for additional studies of shark evolution.

Checklist of marine and estuarine fishes from the AlaskaYukon Border, Beaufort Sea, to Cabo San Lucas, Mexico

This paper is a checklist of the fishes that have been documented, through both published and unpublished sources, in marine and estuarine waters, and out 200 miles, from the United States-Canadian border on the Beaufort Sea to Cabo San Lucas, Mexico. A minimum of 241 families and 1,644 species are known within this range, including both native and nonnative species. For each of these species, we include maximum size, geographic and depth ranges, whether it is native or nonnative, as well as a brief mention of any taxonomic issues.

General descriptions of the dermis structure of a juvenile whale shark Rhincodon typus from the Gulf of California

The aim of this study is to provide preliminary observations on the microanatomy of Rhincodon typus skin using histology and electron microscopy analyses. Skin biopsies were obtained from a deceased juvenile male shark (548 cm total length) stranded in La Paz, Mexico, during February 2018. The results of this study evidenced the basic structure of the dermal denticles in the epidermis of the trunk of the shark, as well as the composition of the connective tissue in the hypodermis. Histological images of the hypodermis showed a high concentration of collagen fibres, formed by a large number of fine and wavy fibres of compact shape and little intercellular substance.

Individual haplotyping of whale sharks from seawater environmental DNA

Population genetic data can provide valuable information on the demography of a species. For rare and elusive marine megafauna, samples for generating the data are traditionally obtained from tissue biopsies, which can be logistically difficult and expensive to collect and require invasive sampling techniques. Analysis of environmental DNA (eDNA) offers an alternative, minimally invasive approach to provide important genetic information. Although eDNA approaches have been studied extensively for species detection and biodiversity monitoring in metabarcoding studies, the potential for the technique to address population-level questions remains largely unexplored. Here, we applied "eDNA haplotyping" to obtain estimates of the intraspecific genetic diversity of a whale shark (Rhincodon typus) aggregation at Ningaloo reef, Australia. Over 2 weeks, we collected seawater samples directly behind individual sharks prior to taking a tissue biopsy sample from the same animal. Our data showed a 100% match between mtDNA sequences recovered in the eDNA and tissue sample for all 28 individuals sampled. In the seawater samples, >97% of all reads were assigned to six dominant haplotypes, and a clear dominant signal (~99% of sample reads) was recovered in each sample. Our study demonstrates accurate individual-level haplotyping from seawater eDNA. When DNA from one individual clearly dominates each eDNA sample, it provides many of the same opportunities for population genetic analyses as a tissue sample, potentially removing the need for tissue sampling. Our results show that eDNA approaches for population-level analyses have the potential to supply critical demographic data for the conservation and management of marine megafauna.

Regional movements of satellite-tagged whale sharks Rhincodon typus in the Gulf of Aden

To gain insight into whale shark (Rhincodon typus) movement patterns in the Western Indian Ocean, we deployed eight pop-up satellite tags at an aggregation site in the Arta Bay region of the Gulf of Tadjoura, Djibouti in the winter months of 2012, 2016, and 2017. Tags revealed movements ranging from local-scale around the Djibouti aggregation site, regional movements along the coastline of Somaliland, movements north into the Red Sea, and a large-scale (>1,000 km) movement to the east coast of Somalia, outside of the Gulf of Aden. Vertical movement data revealed high occupation of the top ten meters of the water column, diel vertical movement patterns, and deep diving behavior. Long-distance movements recorded both here and in previous studies suggest that connectivity between the whale sharks tagged at the Djibouti aggregation and other documented aggregations in the region are likely within annual timeframes. In addition, wide-ranging movements through multiple nations, as well as the high use of surface waters recorded, likely exposes whale sharks in this region to several anthropogenic threats, including targeted and bycatch fisheries and ship-strikes. Area-based management approaches focusing on seasonal hotspots offer a way forward in the conservation of whale sharks in the Western Indian Ocean.

The use of machine learning to detect foraging behaviour in whale sharks: a new tool in conservation

In this study we present the first attempt at modelling the feeding behaviour of whale sharks using a machine learning analytical method. A total of eight sharks were monitored with tri-axial accelerometers and their foraging behaviours were visually observed. Our results highlight that the random forest model is a valid and robust approach to predict the feeding behaviour of the whale shark. In conclusion this novel approach exposes the practicality of this method to serve as a conservation tool and the capability it offers in monitoring potential disturbances of the species.

Donsol: an important reproductive habitat for the world's largest fish Rhincodon typus?

Little is known about the reproductive ecology of the whale shark Rhincodon typus. On 15 March 2020, a free-swimming neonate whale shark was found in the shallow, coastal waters of Donsol, Philippines. The total length of the male shark was 60 cm, falling within the species known size-at-birth. This is the third occurrence of neonatal whale sharks reported in the area, and coupled with the occurrence of juveniles and adults, and anecdotal mating and precopulatory behaviour reported herein, the importance of the site for this endangered species is highlighted.

The impact of injury on apparent survival of whale sharks (Rhincodon typus) in South Ari Atoll Marine Protected Area, Maldives

The whale shark (Rhincodon typus) is an endangered species with a declining global population. The South Ari Atoll Marine Protected Area (SAMPA), Maldives, is one of few locations globally where year-long residency of individuals occurs. This SAMPA aggregation appears to consist almost exclusively of immature males. Due to its year-round residency, this local aggregation is subjected to a high degree of tourism pressure. This ecotourism contributes to the high level of interest and protection offered to whale sharks by the local community. Unfortunately, if regulations are not followed or enforced, tourism can bring with it major stressors, such as accidental injuries. We used POPAN capture-mark-recapture models and lagged identification rate analysis to assess the effect of major injuries on whale shark residency within SAMPA. Injuries may be obtained outside SAMPA. We found individuals with major injuries had a higher apparent survival in the area than those without. Lagged identification rates also demonstrated that sharks with major injuries are more likely to return to the area. We suggest that major injuries result in sharks prolonging their time in the developmental habitat. These findings have implications for individual fitness and the population viability of this endangered species. We propose targeted conservation strategies be considered to protect sharks from further injury. Based on the presented spatio-temporal distributions of sharks, and current local knowledge of sighting patterns, speed limit zones and propeller-exclusion zones should be implemented and enforced. If carried out alongside tourist education, these measures will contribute to the protection of whale sharks within SAMPA and beyond. Furthermore, our results can aid research direction, alongside regulation and enforcement development, at similar sites worldwide.

In-water observations highlight the effects of provisioning on whale shark behaviour at the world's largest whale shark tourism destination

The whale shark is the world's largest fish that forms predictable aggregations across its range, many of which support tourism industries. The largest non-captive provisioned whale shark destination globally is at Oslob, Philippines, where more than 500 000 tourists visit yearly. There, the sharks are provisioned daily, year-round, allowing the human-shark interaction in nearshore waters. We used in-water behavioural observations of whale sharks between 2015 and 2017 to understand the relationship between external stimuli and shark behaviour, whether frequency of visits at the site can act as a predictor of behaviour, and the tourist compliance to the code of conduct. Mixed effects models revealed that the number of previous visits at the site was a strong predictor of whale shark behaviour, and that provisioned sharks were less likely to exhibit avoidance. Compliance was poor, with 93% of surveys having people less than 2 m from the animal, highlighting overcrowding of whale sharks at Oslob. Given the behavioural implications to whale sharks highlighted here and the local community's reliance on the tourism industry, it is imperative to improve management strategies to increase tourist compliance and strive for sustainable tourism practices.

Cell culture-based karyotyping of orectolobiform sharks for chromosome-scale genome analysis

Karyotyping, traditionally performed using cytogenetic techniques, is indispensable for validating genome assemblies whose sequence lengths can be scaled up to chromosome sizes using modern methods. Karyotype reports of chondrichthyans are scarce because of the difficulty in cell culture. Here, we focused on carpet shark species and the culture conditions for fibroblasts and lymphocytes. The utility of the cultured cells enabled the high-fidelity characterization of their karyotypes, namely 2n = 102 for the whale shark (Rhincodon typus) and zebra shark (Stegostoma fasciatum), and 2n = 106 for the brownbanded bamboo shark (Chiloscyllium punctatum) and whitespotted bamboo shark (C. plagiosum). We identified heteromorphic XX/XY sex chromosomes for the two latter species and demonstrated the first-ever fluorescence in situ hybridization of shark chromosomes prepared from cultured cells. Our protocols are applicable to diverse chondrichthyan species and will deepen the understanding of early vertebrate evolution at the molecular level.

Changes in diving behaviour and habitat use of provisioned whale sharks: implications for management

Whale shark (Rhincodon typus) tourism is increasingly popular at predictable aggregations around the world, but only a few use provisioning to ensure close interactions. Understanding the effects of provisioning on the behaviour of this endangered species is critical to manage this growing industry. We recorded the diving behaviour and habitat use of juvenile whale sharks (n = 4) for a mean of 49.5 provisioned and 33.8 non-provisioned days using temperature-depth-recorders. We found that time spent at the surface (< 2 m) between 6 am and 1 pm increased ~ sixfold, while timing of deep dives shifted from 4–10 am to 10 am–2 pm, i.e. near or at the end of the provisioning activities. The shift might be related to a need to thermoregulate following a prolonged period of time in warmer water. These changes could have fitness implications for individuals frequently visiting the provisioning site. Based on recorded amount of time spent in warm waters and published Q₁₀ values for ectotherms, we estimate a 7.2 ± 3.7% (range 1.3–17.8%) higher metabolic rate when sharks frequent the provisioning site. The observed behavioural, habitat use, and potential fitness shifts should be considered when developing guidelines for sustainable tourism, particularly in light of new provisioning sites developing elsewhere.

Whale sharks Rhincodon typus get cleaned by the blue-streak cleaner wrasse Labroides dimidiatus and the moon wrasse Thalassoma lunare in the Philippines

Cleaning interactions are essential for healthy marine ecosystem communities. This study reports the first documentation of the whale shark Rhincodon typus cleaning behaviour in the Indo-West Pacific by two wrasse species, the blue-streak cleaner wrasse Labroides dimidiatus and the moon wrasse Thalassoma lunare in Cebu, Philippines. This study documented 36 cleaning interactions with 14 individual whale sharks. The cleaning interactions appear opportunistic rather than targeted by the sharks, unlike that observed in other species of elasmobranchs. Further work should focus on understanding the drivers of these unique cleaning interactions.

The whale shark genome reveals how genomic and physiological properties scale with body size

The endangered whale shark (Rhincodon typus) is the largest fish on Earth and a long-lived member of the ancient Elasmobranchii clade. To characterize the relationship between genome features and biological traits, we sequenced and assembled the genome of the whale shark and compared its genomic and physiological features to those of 83 animals and yeast. We examined the scaling relationships between body size, temperature, metabolic rates, and genomic features and found both general correlations across the animal kingdom and features specific to the whale shark genome. Among animals, increased lifespan is positively correlated to body size and metabolic rate. Several genomic traits also significantly correlated with body size, including intron and gene length. Our large-scale comparative genomic analysis uncovered general features of metazoan genome architecture: Guanine and cytosine (GC) content and codon adaptation index are negatively correlated, and neural connectivity genes are longer than average genes in most genomes. Focusing on the whale shark genome, we identified multiple features that significantly correlate with lifespan. Among these were very long gene length, due to introns being highly enriched in repetitive elements such as CR1-like long interspersed nuclear elements, and considerably longer neural genes of several types, including connectivity, activity, and neurodegeneration genes. The whale shark genome also has the second slowest evolutionary rate observed in vertebrates to date. Our comparative genomics approach uncovered multiple genetic features associated with body size, metabolic rate, and lifespan and showed that the whale shark is a promising model for studies of neural architecture and lifespan.

First concurrent assessment of elemental- and organic-contaminant loads in skin biopsies of whale sharks from Djibouti

The Gulf of Tadjoura (Djibouti) is an important site where the whale shark (Rhincodon typus Smith, 1828) aggregates seasonally. Because of the proximity to the port of Djibouti that is located along one of the busiest shipping areas in the world, whale sharks are potentially exposed to relatively high levels of trace elements and organochlorines. To assess their contamination status, concentrations of 15 trace elements, DDTs and PCBs were concurrently assessed for the first time in skin biopsies of 20 whale sharks. Additionally, 12 zooplankton samples were collected and analyzed for trace elements content. Concentration of As, Cu, Zn and Se were higher than in previous studies on this species. Whale shark samples exceeded the maximum allowable limits for foodstuffs for Cr, Pb Se, Cd and Zn. Results from this study suggests that Cr, Ni and Mo may biomagnify in this species. With regard to PCBs, the predominant congener were Tetra-CB, accounting for 41% of total PCBs and Penta-CB for 23%, while for DDTs, the predominant metabolite was DDE contributing for 51% of total compounds. The maximum residue limit for DDTs and for the ind-PCBs in fish set by U.S. EPA and by the EU regulation, respectively, was exceeded in 62% of whale sharks. Moreover, the p,p'DDE/ΣDDT ratios were higher than the critical value of 0.6, indicating possible recent inputs of technical DDTs in the area. This study suggests some concerns for the contamination status of whale sharks in Djibouti and, since major threats for this species include bycatch and illegal fisheries, highlights some level of risk from the exposure to elemental- and organic-contaminant via shark consumption.

Whale sharks increase swimming effort while filter feeding, but appear to maintain high foraging efficiencies

Whale sharks (Rhincodon typus) - the largest extant fish species - reside in tropical environments, making them an exception to the general rule that animal size increases with latitude. How this largest fish thrives in tropical environments that promote high metabolism but support less robust zooplankton communities has not been sufficiently explained. We used open-source inertial measurement units (IMU) to log 397 h of whale shark behavior in Yucatán, Mexico, at a site of both active feeding and intense wildlife tourism. Here we show that the strategies employed by whale sharks to compensate for the increased drag of an open mouth are similar to ram feeders five orders of magnitude smaller and one order of magnitude larger. Presumed feeding constituted 20% of the total time budget of four sharks, with individual feeding bouts lasting up to 11 consecutive hours. Compared with normal, sub-surface swimming, three sharks increased their stroke rate and amplitude while surface feeding, while one shark that fed at depth did not demonstrate a greatly increased energetic cost. Additionally, based on time-depth budgets, we estimate that aerial surveys of shark populations should consider including a correction factor of 3 to account for the proportion of daylight hours that sharks are not visible at the surface. With foraging bouts generally lasting several hours, interruptions to foraging during critical feeding periods may represent substantial energetic costs to these endangered species, and this study presents baseline data from which management decisions affecting tourist interactions with whale sharks may be made.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

Photo-ID and telemetry highlight a global whale shark hotspot in Palawan, Philippines

The Philippines is home to the second largest known population of whale sharks in the world. The species is listed as endangered due to continued population declines in the Indo-Pacific. Knowledge about the connectivity within Southeast Asia remains poor, and thus international management is difficult. Here, we employed pop-up archival tags, data mining and dedicated effort to understand an aggregation of whale sharks at Honda Bay, Palawan, Philippines, and its role in the species' conservation. Between Apr and Oct 2018, we conducted 159 surveys identifying 117 individual whale sharks through their unique spot patterns (96.5% male, mean 4.5 m). A further 66 individual whale sharks were identified from local operators, and data mined on social media platforms. The satellite telemetry data showed that the whale sharks moved broadly, with one individual moving to Sabah, Malaysia, before returning to the site <1 year later. Similarly, another tagged whale shark returned to the site at a similar periodicity after reaching the Malay-Filipino border. One individual whale shark first identified in East Kalimantan, Indonesia by a citizen scientist was resighted in Honda Bay ~3.5 years later. Honda Bay is a globally important site for the endangered whale shark with connectivity to two neighbouring countries, highlighting the need for international cooperation to manage the species.

Trace Elements in Tissues of Whale Sharks (Rhincodon typus) Stranded in the Gulf of California, Mexico

Concentration of essential (Se, Zn and Cu) and non-essential (As, Cd, Hg and Pb) trace elements were measured in selected tissues of two dead whale sharks (Rhincodon typus) stranded in the Gulf of California (GC) in 2017 and 2018. Concentrations of Cd and Pb in the skeletal muscle of the whale shark from La Paz Bay, GC were higher compared to a previous study on whale shark from China. The shark from La Paz Bay also presented higher concentration of Pb in the epidermis, compared to the same tissue of the other whale shark stranded in Punta Bufeo, GC. The Hg in all analysed tissues was lower than those documented in carnivorous sharks. Molar ratio Se:Hg shows an excess of Se over Hg in all the tissues sampled in both sharks.

Seasonal trends in whale shark Rhincodon typus sightings in an established tourism site in the Gulf of California, Mexico

The number of individual whale shark Rhincodon typus sightings registered October 2015-March 2018 in Bahia de La Paz, Mexico, ranged from 73 to 129 animals per season (total 1662), with the majority of them identified as males. Density plot maps showed high correlations in sightings between sample seasons and revealed potential hotspot areas. Our study provides an essential baseline of information for the management of whale shark conservation and the associated touristic activities in the region.

Mercury and selenium in the filter-feeding whale shark (Rhincodon typus) from two areas of the Gulf of California, Mexico

Mercury and selenium were measured for first time in the endangered species whale shark (Rhyncodon typus) from two areas of the Gulf of California, Bahía Los Angeles (BLA) and Bahía La Paz (LAP) using dermal biopsies of seventy specimens. Additionally, nineteen zooplankton samples from LAP were analysed. Concentrations (ng/g, wet weight (ww)) in biopsies of BLA ranged from 1 to 40 for Hg and 100 to 680 for Se; while in LAP varied from 1 to 9 for Hg and 11 to 850 for Se. A positive correlation was found for Hg in BLA males biopsies with length. Hg and Se concentrations in the zooplankton from LAP were 1.6 ± 1.8 and 770 ± 930 ng/g, respectively. Hg biomagnification factor ranged from 0.8 to 5.3 in sharks. A molar excess of Se over Hg was found in the biopsies and the zooplankton.

Observations of stranded and swimming whale sharks Rhincodon typus in Khor Al-Zubair, NW Arabian Gulf and Shatt Al-Arab Estuary, Iraq

This study is the first record of a whale shark Rhincodon typus in Iraqi waters and the second record from the Arabian (Persian) Gulf portion of the north-west Indian Ocean. The stranding and sighting events were documented by photos of three individuals: the stranded specimen was 6.5 m total length (L_T ) and the free-swimming individuals were c. 1-2 and 2-3 m L_T respectively_. An appeal is made for Iraqi policy makers to participate in regional and international organizations for the conservation of this endangered species.

Association of the whale shark Rhincodon typus with the cannonball jellyfish Stomolophus meleagris

This study describes the first record of the whale shark Rhincodon typus association with the cannonball jellyfish Stomolophus meleagris in the Gulf of California, Mexico. Whale sharks were observed swimming and feeding among swarms of jellyfish, suggesting competition and predatory behaviour given the overlap in food preferences between both species. This finding is relevant because of the species-wide distribution and the importance of these interactions, which should be considered in conservation strategies of R. typus and management of cannonball jellyfish fisheries.

Satellite tracking of juvenile whale sharks in the Sulu and Bohol Seas, Philippines

The whale shark Rhincodon typus was uplisted to ‘Endangered’ in the 2016 IUCN Red List due to >50% population decline, largely caused by continued exploitation in the Indo-Pacific. Though the Philippines protected the whale shark in 1998, concerns remain due to continued take in regional waters. In light of this, understanding the movements of whale sharks in the Philippines, one of the most important hotspots for the species, is vital. We tagged 17 juvenile whale sharks with towed SPOT5 tags from three general areas in the Sulu and Bohol Seas: Panaon Island in Southern Leyte, northern Mindanao, and Tubbataha Reefs Natural Park (TRNP). The sharks all remained in Philippine waters for the duration of tracking (6–126 days, mean 64). Individuals travelled 86–2,580 km (mean 887 km) at a mean horizontal speed of 15.5 ± 13.0 SD km day⁻¹. Whale sharks tagged in Panaon Island and Mindanao remained close to shore but still spent significant time off the shelf (>200 m). Sharks tagged at TRNP spent most of their time offshore in the Sulu Sea. Three of twelve whale sharks tagged in the Bohol Sea moved through to the Sulu Sea, whilst two others moved east through the Surigao Strait to the eastern coast of Leyte. One individual tagged at TRNP moved to northern Palawan, and subsequently to the eastern coast of Mindanao in the Pacific Ocean. Based on inferred relationships with temperature histograms, whale sharks performed most deep dives (>200 m) during the night, in contrast to results from whale sharks elsewhere. While all sharks stayed in national waters, our results highlight the high mobility of juvenile whale sharks and demonstrate their connectivity across the Sulu and Bohol Seas, highlighting the importance of the area for this endangered species.

First results on zooplankton community composition and contamination by some persistent organic pollutants in the Gulf of Tadjoura (Djibouti)

The Gulf of Tadjoura is located in the Horn of Africa and is widely recognized as an important site where the zooplanktivorous whale sharks seasonally aggregate from October to February. The surface zooplankton community (0-3m) was weekly sampled from November 2016 to February 2017 in two sites during the whale shark aggregation period. A total of 12 phyla were identified. Copepoda represented the most abundant and diverse group with 29 different genera, and contributed with an average of 82% of the mean zooplankton density of approximately 6600indm^-3. During the sampling period, copepods were dominated numerically by Calanoida (3600indm^-3), followed by Poicilostomatatoida (1300indm^-3). Within the copepods, Paracalanidae, Calanidae, Oncaeidae and Miraciidae were the most common families. The temporal trend in zooplankton biomass at both stations revealed the highest peak in December (41.3±36.4mgm^-3), and the lowest in February (6.6±3.3mgm^-3). As no information is available on the occurrence of legacy contaminants use and release in this area, analysis revealed the consistent presence of both DDT and PCB residues in zooplankton samples in the Gulf of Tadjoura. Total PCB ranged from approximately 110 to 637ngg^-1 d.w., while total DDT from 21 to 80ngg^-1 d.w. The proportion of primary DDT in the total residue was higher than DDE and DDD, which strongly suggests that the area might actually be subjected to DDT inputs of the parent compound.

Does bathymetry drive coastal whale shark (Rhincodon typus) aggregations?

BACKGROUND

The whale shark (Rhincodon typus) is known to aggregate in a number of coastal locations globally, however what causes these aggregations to form where they do is largely unknown. This study examines whether bathymetry is an important driver of coastal aggregation locations for R. typus through bathymetry's effect on primary productivity and prey availability. This is a global study taking into account all coastal areas within R. typus' range.

METHODS

R. typus aggregation locations were identified through an extensive literature review. Global bathymetric data were compared at R. typus aggregation locations and a large random selection of non-aggregation areas. Generalised linear models were used to assess which bathymetric characteristic had the biggest influence on aggregation presence.

RESULTS

Aggregation sites were significantly shallower than non-aggregation sites and in closer proximity to deep water (the mesopelagic zone) by two orders of magnitude. Slope at aggregation sites was significantly steeper than non-aggregation sites. These three bathymetric variables were shown to have the biggest association with aggregation sites, with up to 88% of deviation explained by the GLMs.

DISCUSSION

The three key bathymetric characteristics similar at the aggregation sites are known to induce upwelling events, increase primary productivity and consequently attract numerous other filter feeding species. The location of aggregation sites in these key areas can be attributed to this increased prey availability, thought to be the main reason R. typus aggregations occur, extensively outlined in the literature. The proximity of aggregations to shallow areas such as reefs could also be an important factor why whale sharks thermoregulate after deep dives to feed. These findings increase our understanding of whale shark behaviour and may help guide the identification and conservation of further aggregation sites.

Satellite tagging highlights the importance of productive Mozambican coastal waters to the ecology and conservation of whale sharks

The whale shark Rhincodon typus is an endangered, highly migratory species with a wide, albeit patchy, distribution through tropical oceans. Ten aerial survey flights along the southern Mozambican coast, conducted between 2004–2008, documented a relatively high density of whale sharks along a 200 km stretch of the Inhambane Province, with a pronounced hotspot adjacent to Praia do Tofo. To examine the residency and movement of whale sharks in coastal areas around Praia do Tofo, where they may be more susceptible to gill net entanglement, we tagged 15 juveniles with SPOT5 satellite tags and tracked them for 2–88 days (mean = 27 days) as they dispersed from this area. Sharks travelled between 10 and 2,737 km (mean = 738 km) at a mean horizontal speed of 28 ± 17.1 SD km day⁻¹. While several individuals left shelf waters and travelled across international boundaries, most sharks stayed in Mozambican coastal waters over the tracking period. We tested for whale shark habitat preferences, using sea surface temperature, chlorophyll-a concentration and water depth as variables, by computing 100 random model tracks for each real shark based on their empirical movement characteristics. Whale sharks spent significantly more time in cooler, shallower water with higher chlorophyll-a concentrations than model sharks, suggesting that feeding in productive coastal waters is an important driver of their movements. To investigate what this coastal habitat choice means for their conservation in Mozambique, we mapped gill nets during two dedicated aerial surveys along the Inhambane coast and counted gill nets in 1,323 boat-based surveys near Praia do Tofo. Our results show that, while whale sharks are capable of long-distance oceanic movements, they can spend a disproportionate amount of time in specific areas, such as along the southern Mozambique coast. The increasing use of drifting gill nets in this coastal hotspot for whale sharks is likely to be a threat to regional populations of this iconic species.

Some like it hot: Repeat migration and residency of whale sharks within an extreme natural environment

The Arabian Gulf is the warmest sea in the world and is host to a globally significant population of the whale shark Rhincodon typus. To investigate regional whale shark behaviour and movements, 59 satellite-linked tags were deployed on whale sharks in the Al Shaheen area off Qatar from 2011–14. Four different models of tag were used throughout the study, each model able to collect differing data or quantities of data. Retention varied from one to 227 days. While all tagged sharks crossed international maritime boundaries, they typically stayed within the Arabian Gulf. Only nine sharks dispersed through the narrow Strait of Hormuz into the Gulf of Oman. Most sharks stayed close to known or suspected feeding aggregation sites over summer months, but dispersed throughout the Arabian Gulf in winter. Sharks rarely ventured into shallow areas (<40 m depth). A single, presumably pregnant female shark was the sole animal to disperse a long distance, crossing five international maritime boundaries in 37 days before the tag detached at a distance of approximately 2644 km from the tagging site, close to the Yemeni-Somali border. No clear space-use differentiation was evident between years, for sharks of different sizes, or between sexes. Whale sharks spent the most time (~66%) in temperatures of 24–30°C and in shallow waters <100 m depth (~60%). Sharks spent relatively more time in cooler (X² = 121.692; p<0.05) and deeper (X² = 46.402; p<0.05) water at night. Sharks rarely made dives deeper than 100 m, reflecting the bathymetric constraints of the Gulf environment. Kernel density analysis demonstrated that the tagging site at Al Shaheen was the regional hotspot for these sharks, and revealed a probable secondary aggregation site for whale sharks in nearby Saudi Arabian waters. Analysis of visual re-sightings data of tagged sharks revealed that 58% of tagged individuals were re-sighted back in Al Shaheen over the course of this study, with 40% recorded back at Al Shaheen in the year following their initial identification. Two sharks were confirmed to return to Al Shaheen in each of the five years of study.

Are whale sharks exposed to persistent organic pollutants and plastic pollution in the Gulf of California (Mexico)? First ecotoxicological investigation using skin biopsies

The whale shark (Rhincodon typus) is an endangered species that may be exposed to micro- and macro-plastic ingestion as a result of their filter-feeding activity, particularly on the sea surface. In this pilot project we perform the first ecotoxicological investigation on whale sharks sampled in the Gulf of California exploring the potential interaction of this species with plastic debris (macro-, micro-plastics and related sorbed contaminants). Due to the difficulty in obtaining stranded specimens of this endangered species, an indirect approach, by skin biopsies was used for the evaluation of the whale shark ecotoxicological status. The levels of organochlorine compounds (PCBs, DDTs), polybrominated diphenyl ethers (PBDEs) plastic additives, and related biomarkers responses (CYP1A) were investigated for the first time in the whale shark. Twelve whale shark skin biopsy samples were collected in January 2014 in La Paz Bay (BCS, Mexico) and a preliminary investigation on microplastic concentration and polymer composition was also carried out in seawater samples from the same area. The average abundance pattern for the target contaminants was PCBs>DDTs>PBDEs>HCB. Mean concentration values of 8.42ng/g w.w. were found for PCBs, 1.31ng/g w.w. for DDTs, 0.29ng/g w.w. for PBDEs and 0.19ng/g w.w. for HCB. CYP1A-like protein was detected, for the first time, in whale shark skin samples. First data on the average density of microplastics in the superficial zooplankton/microplastic samples showed values ranging from 0.00items/m³ to 0.14items/m³. A focused PCA analysis was performed to evaluate a possible correlation among the size of the whale sharks, contaminants and CYP1A reponses. Further ecotoxicological investigation on whale shark skin biopsies will be carried out for a worldwide ecotoxicological risk assessment of this endangerd species.

Association of whale sharks (Rhincodon typus) with thermo-biological frontal systems of the eastern tropical Pacific

Satellite tracking of 27 whale sharks in the eastern tropical Pacific, examined in relation to environmental data, indicates preferential occupancy of thermo-biological frontal systems. In these systems, thermal gradients are caused by wind-forced circulation and mixing, and biological gradients are caused by associated nutrient enrichment and enhanced primary productivity. Two of the frontal systems result from upwelling, driven by divergence in the current systems along the equator and the west coast of South America; the third results from wind jet dynamics off Central America. All whale sharks were tagged near Darwin Island, Galápagos, within the equatorial Pacific upwelling system. Occupancy of frontal habitat is pronounced in synoptic patterns of shark locations in relation to serpentine, temporally varying thermal fronts across a zonal expanse > 4000 km. 80% of shark positions in northern equatorial upwelling habitat and 100% of positions in eastern boundary upwelling habitat were located within the upwelling front. Analysis of equatorial shark locations relative to thermal gradients reveals occupancy of a transition point in environmental stability. Equatorial subsurface tag data show residence in shallow, warm (>22°C) water 94% of the time. Surface zonal current speeds for all equatorial tracking explain only 16% of the variance in shark zonal movement speeds, indicating that passive drifting is not a primary determinant of movement patterns. Movement from equatorial to eastern boundary frontal zones occurred during boreal winter, when equatorial upwelling weakens seasonally. Off Peru sharks tracked upwelling frontal positions within ~100–350 km from the coast. Off Central America, the largest tagged shark (12.8 m TL) occupied an oceanic front along the periphery of the Panama wind jet. Seasonal movement from waning equatorial upwelling to productive eastern boundary habitat is consistent with underlying trophic dynamics. Persistent shallow residence in thermo-biological frontal zones suggests the role of physical-biological interactions that concentrate food resources.

Long-term assessment of whale shark population demography and connectivity using photo-identification in the Western Atlantic Ocean

The predictable occurrence of whale sharks, Rhincodon typus, has been well documented in several areas. However, information relating to their migratory patterns, residency times and connectivity across broad spatial scales is limited. In the present study photo-identification data is used to describe whale shark population structure and connectivity among known aggregation sites within the Western Central Atlantic Ocean (WCA). From 1999 to 2015, 1,361 individuals were identified from four distinct areas: the Yucatan Peninsula, Mexico (n = 1,115); Honduras (n = 146); northern Gulf of Mexico, United States (n = 112), and Belize (n = 49). Seasonal patterns in whale shark occurrence were evident with encounters occurring in the western Caribbean Sea earlier in the year than in the GOM. There was also a significant sex bias with 2.6 times more males present than females. Seventy sharks were observed in more than one area and the highest degree of connectivity occurred among three aggregation sites along the Mesoamerican Reef. Despite this, the majority of resightings occurred in the area where the respective sharks were first identified. This was true for the WCA as a whole, with the exception of Belize. Site fidelity was highest in Mexico. Maximum likelihood modelling resulted in a population estimate of 2,167 (95% c.i. 1585.21-2909.86) sharks throughout the entire region. This study is the first attempt to provide a broad, regional population estimate using photo-identification data from multiple whale shark aggregations. Our aim is to provide population metrics, along with the description of region-scale connectivity, that will help guide conservation action in the WCA. At a global level, rapidly growing photographic databases are allowing for researchers to look beyond the description of single aggregation sites and into the ocean-scale ecology of this pelagic species.

Draft sequencing and assembly of the genome of the world's largest fish, the whale shark: Rhincodon typus Smith 1828

BACKGROUND

The whale shark (Rhincodon typus) has by far the largest body size of any elasmobranch (shark or ray) species. Therefore, it is also the largest extant species of the paraphyletic assemblage commonly referred to as fishes. As both a phenotypic extreme and a member of the group Chondrichthyes - the sister group to the remaining gnathostomes, which includes all tetrapods and therefore also humans - its genome is of substantial comparative interest. Whale sharks are also listed as an endangered species on the International Union for Conservation of Nature's Red List of threatened species and are of growing popularity as both a target of ecotourism and as a charismatic conservation ambassador for the pelagic ecosystem. A genome map for this species would aid in defining effective conservation units and understanding global population structure.

RESULTS

We characterised the nuclear genome of the whale shark using next generation sequencing (454, Illumina) and de novo assembly and annotation methods, based on material collected from the Georgia Aquarium. The data set consisted of 878,654,233 reads, which yielded a draft assembly of 1,213,200 contigs and 997,976 scaffolds. The estimated genome size was 3.44Gb. As expected, the proteome of the whale shark was most closely related to the only other complete genome of a cartilaginous fish, the holocephalan elephant shark. The whale shark contained a novel Toll-like-receptor (TLR) protein with sequence similarity to both the TLR4 and TLR13 proteins of mammals and TLR21 of teleosts. The data are publicly available on GenBank, FigShare, and from the NCBI Short Read Archive under accession number SRP044374.

CONCLUSIONS

This represents the first shotgun elasmobranch genome and will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this and other shark species, as well as providing comparative data for studies of evolutionary biology and immunology across the jawed vertebrate lineages.