Emergent community architecture despite distinct diversity in the global whale shark (Rhincodon typus) epidermal microbiome.

Microbiomes confer beneficial physiological traits to their host, but microbial diversity is inherently variable, challenging the relationship between microbes and their contribution to host health. Here, we compare the diversity and architectural complexity of the epidermal microbiome from 74 individual whale sharks (Rhincodon typus) across five aggregations globally to determine if network properties may be more indicative of the microbiome-host relationship. On the premise that microbes are expected to exhibit biogeographic patterns globally and that distantly related microbial groups can perform similar functions, we hypothesized that microbiome co-occurrence patterns would occur independently of diversity trends and that keystone microbes would vary across locations. We found that whale shark aggregation was the most important factor in discriminating taxonomic diversity patterns. Further, microbiome network architecture was similar across all aggregations, with degree distributions matching Erdos-Renyi-type networks. The microbiome-derived networks, however, display modularity indicating a definitive microbiome structure on the epidermis of whale sharks. In addition, whale sharks hosted 35 high-quality metagenome assembled genomes (MAGs) of which 25 were present from all sample locations, termed the abundant 'core'. Two main MAG groups formed, defined here as Ecogroup 1 and 2, based on the number of genes present in metabolic pathways, suggesting there are at least two important metabolic niches within the whale shark microbiome. Therefore, while variability in microbiome diversity is high, network structure and core taxa are inherent characteristics of the epidermal microbiome in whale sharks. We suggest the host-microbiome and microbe-microbe interactions that drive the self-assembly of the microbiome help support a functionally redundant abundant core and that network characteristics should be considered when linking microbiomes with host health.

Microplastic occurrence and phthalate ester levels in neuston samples and skin biopsies of filter-feeding megafauna from La Paz Bay (Mexico).

The impacts of microplastics on filter feeders megafauna have recently received increased attention. These organisms are potentially exposed to plastic ingestion and the release of added/sorbed contaminants during feeding activities. An assessment of microplastic abundance and the chemical impact of Phthalates esters (PAEs) were performed in neustonic samples and skin biopsies of Balaenoptera physalus and Rhincodon typus inhabiting the Gulf of California (Mexico). Sixty-eight percent of the net tows contained plastics with a maximum of 0.24 items/m3 mainly composed of polyethylene fragments. PAE levels were detected both in environmental and skin biopsy samples, with the highest values in the fin whale specimens (5291 ng/g d.w). Plasticizer fingerprint showed a similar distribution pattern between neustonic samples and filter-feeding species, with DEHP and MBP having the highest concentrations. The detection of PAE levels confirmed their potential role as plastic tracers and give preliminary information about the toxicological status of these species feeding in La Paz Bay.

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.

Oceanic adults, coastal juveniles: tracking the habitat use of whale sharks off the Pacific coast of Mexico.

Eight whale sharks tagged with pop-up satellite archival tags off the Gulf of California, Mexico, were tracked for periods of 14-134 days. Five of these sharks were adults, with four females visually assessed to be pregnant. At least for the periods they were tracked, juveniles remained in the Gulf of California while adults moved offshore into the eastern Pacific Ocean. We propose that parturition occurs in these offshore waters. Excluding two juveniles that remained in the shallow tagging area for the duration of tracking, all sharks spent 65 ± 20.7% (SD) of their time near the surface, even over deep water, often in association with frontal zones characterized by cool-water upwelling. While these six sharks all made dives into the meso- or bathypelagic zones, with two sharks reaching the maximum depth recordable by the tags (1285.8 m), time spent at these depths represented a small proportion of the overall tracks. Most deep dives (72.7%) took place during the day, particularly during the early morning and late afternoon. Pronounced habitat differences by ontogenetic stage suggest that adult whale sharks are less likely to frequent coastal waters after the onset of maturity.

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/m3 to 0.14items/m3. 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.

Genetic structure of populations of whale sharks among ocean basins and evidence for their historic rise and recent decline.

This study presents genetic evidence that whale sharks, Rhincodon typus, are comprised of at least two populations that rarely mix and is the first to document a population expansion. Relatively high genetic structure is found when comparing sharks from the Gulf of Mexico with sharks from the Indo-Pacific. If mixing occurs between the Indian and Atlantic Oceans, it is not sufficient to counter genetic drift. This suggests whale sharks are not all part of a single global metapopulation. The significant population expansion we found was indicated by both microsatellite and mitochondrial DNA. The expansion may have happened during the Holocene, when tropical species could expand their range due to sea-level rise, eliminating dispersal barriers and increasing plankton productivity. However, the historic trend of population increase may have reversed recently. Declines in genetic diversity are found for 6 consecutive years at Ningaloo Reef in Australia. The declines in genetic diversity being seen now in Australia may be due to commercial-scale harvesting of whale sharks and collision with boats in past decades in other countries in the Indo-Pacific. The study findings have implications for models of population connectivity for whale sharks and advocate for continued focus on effective protection of the world's largest fish at multiple spatial scales.

Isolation and characterization of microsatellite loci in the whale shark (Rhincodon typus).

In preparation for a study on population structure of the whale shark (Rhincodon typus), nine species-specific polymorphic microsatellite DNA markers were developed. An initial screening of 50 individuals from Holbox Island, Mexico found all nine loci to be polymorphic, with two to 17 alleles observed per locus. Observed and expected heterozygosity per locus ranged from 0.200 to 0.826 and from 0.213 to 0.857, respectively. Neither statistically significant deviations from Hardy-Weinberg expectations nor statistically significant linkage disequilibrium between loci were observed. These microsatellite loci appear suitable for examining population structure, kinship assessment and other applications.