Organic contaminants as an ecological tool to explore niche partitioning: a case study using three pelagic shark species

Chemical contaminant profiles are linked to an animal’s niche, providing a potential tool by which to assess resource partitioning in pelagic species. As proof of concept, we examined contaminant signatures in three species of sharks (Isurus oxyrinchus, Prionace glauca, and Alopias vulpinus) known to overlap in both space and time. Since these sharks comprise a predatory guild within the Southern California Bight (SCB), we predicted that species may partition spatial and dietary resources to limit the extent of competitive exclusion. Indeed, species were distinguishable by both total contaminant loads and their contaminant fingerprint, as random forest analysis found that species could be correctly classified 96% of the time. Our results demonstrate the utility of chemical analyses for ecological studies, and how contaminant tracers can be used in combination with traditional methods to elucidate how species may undergo niche partitioning to reduce competition for overlapping resources within predatory guilds.

species that have a high propensity to accumulate these contaminants [25][26][27] . Similar to stable isotopes, organic contaminants are acquired mainly through diet and are incorporated into tissues based on each contaminant's physio-chemical properties 28 . In addition, organic contaminants can offer some degree of spatial resolution as different areas have unique chemical signatures based on the history of direct contaminant release or atmospheric deposition 29 . These contaminant signatures are then incorporated into local biota, which will generally reflect the environment, assuming a majority of contaminants from the diet are absorbed and metabolism by the predator is limited so as to not alter these contaminant signatures. For instance, organisms utilizing the Southern California Bight (SCB) have a strong signature of a dichlorodiphenyltrichloroethane (DDT) metabolite (e.g. 4,4′-DDE) due to the high degree of DDT release into the marine environment in the 1900s 30 . Since species' contaminant loads vary based on their use of spatial and trophic resources, chemical analysis is a potential tool to examine resource partitioning in sympatric species while taking aspects of their physiology into consideration 22 .
In the California Current Large Marine Ecosystem (CCMLE), Isurus oxyrinchus (Shortfin Mako), Alopias vulpinus (Common Thresher Shark), and Prionace glauca (Blue Shark) comprise a sympatric predator guild 12,31 , making them excellent candidates to determine whether toxicological tools can aid in our understanding of resource partitioning. Within the CCMLE, the SCB represents an important nursery area for these sharks, as juveniles of all three species are relatively abundant in the region 31,32 . Previous tagging and dietary analyses indicate that some degree of niche partitioning occurs among these three species. Spatially, A. vulpinus are the most coastal 32 , whereas P. glauca and I. oxyrinchus make more regular movements offshore 12 . With respect to diet, all three species demonstrate some degree of prey overlap, but species can be distinguished based on stomach content composition 16 . As such, we would predict these differences to be reflected in their contaminant signatures. In an accompanying study 33 , we previously demonstrate how aspects of species' physiology and ecology can result in different lifetime contaminant accumulation trajectories by measuring a suite of legacy contaminants in liver tissue across a range of animal sizes. P. glauca, the species with the lowest estimated metabolic rate and least time spent nearshore, demonstrated low potential for lifetime accumulation, although adult sampling was limited. While A. vulpinus and I. oxyrinchus both demonstrated patterns of growth dilution during the juvenile stage, only I. oxyrinchus exhibited significant increases in liver contaminant concentration in adults, likely due to their high metabolic rate and relatively high trophic position. Species characteristics that may influence contaminant input (i.e. where animals feed, what they feed on or how often they feed) may also have the potential to influence species' unique contaminant signatures. Therefore, we aimed to demonstrate the applicability of organic contaminants as a tool to examine the degree of niche partitioning among three species of sympatric sharks within the same predatory guild using contaminant information collected previously 33 .

Methods
Sample collection. Liver samples were opportunistically obtained from frozen samples archived at the Southwest Fisheries Science Center's (SWFSC) from the juvenile Shortfin Mako/Blue Shark survey 31 and Common Thresher Shark survey, the National Marine Fisheries Service (NMFS) West Coast Region Fishery Observer Program, participating drift gillnet fishermen, and recreational fishermen. A subset of these archived samples was then analyzed for organic contaminants (see below) in individuals selected to represent a range of sizes and sexes for I. oxyrinchus, P. glauca, and A. vulpinus sampled from the years 2011-2013. Archived samples were also supplemented with contaminant data from previously published reports [34][35][36] from animals sampled in the same time frame (i.e., 2011-2013). organochlorine contaminant analysis. Organic contaminants were analyzed at California State University Long Beach's IIRMES facility following previously published methods 34 . Liver subsamples (0.5-1.0 g wet weight) from each animal were extracted for 14-16 hrs via a Soxhlet apparatus in 100% methylene chloride solution, followed by subsequent evaporation and sample purification by elution through an Alumina-B/ Silica gel. Extracts were spiked with internal standards (4,4′-Dibromobiphenyl and 2,2′,5,5′-Tetrabromobiphenyl, Accustandard, Inc, New Haven, USA) and injected onto an Agilent gas chromatograph (GC; 6890N series) equipped with a mass selective detector (MSD; Agilent 5973 inert series, Santa Clara, USA). Ion peaks were then identified using gas chromatography mass spectrometry (GCMS) software for 54 PCB congeners (sum = tPCBs), DDT and its metabolites (4,4′-DDT, 4,4′-DDE, 4,4′-DDD, 2,4′-DDT, 2,4′-DDE, 2,4′-DDD; sum = tDDXs), and non-DDT chlorinated pesticides (24 compounds screened; sum = tPEST); total concentrations were expressed as the sum of all contaminant groups (tOCs). The limit of detection for all compounds was 1 ng/g. Part way through sample processing, the method at the laboratory where chemical analysis took place to include 4,4′-DDMU, a downstream break down product of 4,4′-DDE; thus, only a subset of samples (n = 44) included this metabolite and no samples from previous studies had this metabolite analyzed. Concentrations of this metabolite are reported for samples where it was measured, but were not included in overall analyses due to incongruency in measured congeners among samples. Lipid content was determined gravimetrically from split aliquots of the extracts; however, since neither wet nor lipid normalized concentrations were correlated with liver lipid content for any species (p ≥ 0.16) we reported values on a wet weight concentration basis.
To ensure quality control/quality assurance, all samples were spiked with recovery surrogate compounds (TCMX, PCB30, PCB112, PCB198) prior to extraction to measure extraction and recovery efficiency of procedures (recovery ranges for all four surrogates: 75 ± 10%). In addition, a certified reference material (Lake Michigan Trout tissue 1947, National Institute of Standards and Technology, Gaithersburg, USA), a pair of blank spikes (i.e. 400 ng of each measured congener), a blank (spiked only with recovery surrogates), and a liver sample replicate were run in tandem with each batch (i.e., ~every 24 samples). All compounds in the certified reference material were within ±30% of the reference value 37 . Recovery of measured compounds (n = 83) in blank spikes (n = 4) were 108 ± 7% for PCBs and 89 ± 20% for pesticides, with 4.2 ± 3.3% relative standard deviation between blank spike replicates. Blanks had high recovery of recovery surrogates (90 ± 28%) and showed little to www.nature.com/scientificreports www.nature.com/scientificreports/ no contamination during procedures (e.g. <2 ng/g of PCB153/PCB138); therefore, liver sample values were not blank-corrected. For liver samples where a replicate was run, the median relative significant difference between concentrations of quantifiable compounds within replicate pairs was 9.9 ± 2.8%, and mean concentrations were used for samples where replicates were available.

Data analysis.
Organic contaminants were compared among species in two ways to determine if species could be distinguished based on their contaminant concentrations or signatures. For our comparison of tOCs, individual contaminant concentrations were summed across compounds for each shark sample and expressed on a wet weight (ww) basis. tOCs were then compared among I. oxyrinchus, P. glauca, and A. vulpinus using a Kruskal-Wallis (KW) test followed by a Wilcoxon rank-sum test, with all sizes and sexes combined within species. Contribution of each contaminant group (i.e. tPCBs, tDDXs, and tPEST) was calculated as the proportion that each contaminant group contributed to the total contaminant load (tOCs). Differences in contaminant proportion was assessed using a regression modeling approach for dependent data with a beta distribution such as proportional data implemented in the betareg package in R (v 3.2.0).
To determine if species could be distinguished based on their contaminant signatures (i.e. relative differences in standardized contaminant proportion), measured contaminants for each individual were standardized by dividing each individual contaminant concentration by the sample's own PCB153 concentration (a contaminant that is consistently detected in all species 38 ) to obtain a ratio for each contaminant per individual. Since some contaminants had much larger concentrations than others (e.g. 4,4′-DDE), ratios were standardized via a z-score transformation by species. A classification random forest analysis was then used to determine if individual shark samples could be correctly assigned by species based on their z-score transformed contaminant profiles alone. Random forest analysis was conducted using the randomForest package in R. To determine species assignment, we set the number of trees (ntree) equal to 10,000, and to avoid biases caused by unequal sample sizes among the different species, we set the number of individuals used in the training set (sampsize) equal to eleven, which is half the number of individuals in the smallest group (P. glauca). To visualize the separation of species based on contaminant signatures we performed a non-metric multidimensional scaling (nMDS; vegan package in R) using contaminant ratios where each contaminant was standardized to the samples own PCB153 concentration and k was set to two.

Results
Species showed significant differences in tOCs (KW, p < 0.001; Fig. 1) despite the inclusion of both sexes and a broad range of sizes, particularly for I. oxyrinchus and A. vulpinus that encompassed young-of-the-year through mature animals. I. oxyrinchus (n = 47, 53 to 337 cm fork length [FL]) had the highest tOCs (mean ± SD: 55,186 ± 66,968 ng −1 g ww), which were approximately 15 and 90 times greater than mean tOCs in A. vulpinus (n = 51, 63 to 283 cm FL, 3,781 ± 2,226 ng −1 g ww) and P. glauca (n = 22, 50 to 178 cm FL, 611 ± 382 ng −1 g ww), respectively. Besides having the lowest tOCs, P. glauca also had the fewest number of detectable contaminants measured (approximately 33 consistently measured contaminants out of a possible 84 screened), whereas A. vulpinus and I. oxyrinchus had much higher numbers (48 and 57, respectively).
While no information is available on the health of individuals sampled for this study, species risk to the potential negative effects of contaminant exposure may not be similar due to the magnitude differences in PCB concentrations across species. With respect to tPCB concentrations eliciting negative effects in mammals [39][40][41][42] , mean I. oxyrinchus consistently exceeded these thresholds regardless of whether wet or lipid normalized values were used, while mean P. glauca concentrations fell under these thresholds. Depending on the study, tPCBs concentrations in A. vulpinus either exceeded or fell below mammalian effect thresholds. Although studies are limited, environmental PCB exposure has the capability of eliciting negative physiological responses in elasmobranchs 43,44 . Compared to these studies, mean tPCBs in all three shark species exceed the lowest wet weight tPCB value in U. halleri and two of the three species exceed lipid normalized concentrations 43 . Considering the disparity in tPCB www.nature.com/scientificreports www.nature.com/scientificreports/ concentrations, it is likely that physiological risk to contaminant exposure is not equitable across these shark species.
PCB compounds were the next highest contributor to tOCs, and contributions were not significantly different among the three species. I. oxyrinchus had the most variability in PCB contribution (33 ± 12.6%), followed by P. glauca (25 ± 8.9%) then A. vulpinus (27 ± 7.8%). Among PCBs that contributed the most to total PCB loads, PCB153 and PCB138 consistently ranked first and second among species; however, the relative contribution of these congeners to tPCBs differed (Table 1). Of the top 10 contributing PCBs, species differed in both PCB congener rank order and relative contributions of these congeners to tPCB loads (Table 1). Detection frequencies also varied among species. I. oxyrinchus had the most congeners detected in one or more samples (91% of 54 congeners measured) as well as consistently higher detection frequencies of congeners across samples (71 ± 30% of samples), followed by A. vulpinus (% congeners detected = 87%; average detection frequency = 60 ± 32%) with the least number (83%) and lowest frequency of detections (36 ± 36%) in P. glauca.
Categorical random forest analysis showed that species were distinguishable based on their contaminant signature and were correctly assigned 96% of the time ( Table 2). Contaminants that contributed to species separation can be found in the Supplemental Table 1. nMDS showed clear separation among the three species in their contaminant signatures (Fig. 2; stress = 0 www.nature.com/scientificreports www.nature.com/scientificreports/

Discussion
Several ecological tools are available for studying habitat use and sympatric species resource partitioning. The most appropriate tool depends on the question of interest, the nature of the ecological community, and the available resources. Exploiting the fact that an animal accrues different chemicals signatures based on its distribution and foraging ecology allows us to use chemical tools to determine resource partitioning in seemingly overlapping species and subsequently make inferences on species' ecological niches. Accumulation of contaminants by predators has consequences for the transport of these chemicals in the environment and the health of their populations. In turn, contaminant accumulation is a product of the interaction of species' niche, life history, and physiological characteristics. Here, we highlight the ability to use contaminants as an ecological tool to assess niche partitioning and in an accompanying study 33 , we further examine how ecological and physiological factors synergistically influence contaminant accumulation trajectories in three species.
Despite spatial and temporal overlap of I. oxyrinchus, P. glauca, and A. vulpinus 16,45 in the CCMLE, we were able to discriminate among these three species with high accuracy using both total contaminant loads and profiles. This study corroborates prior findings showing significant differences in ecological niches among these three species in the same area 16 . Preti et al. 16 found differences in stomach contents with A. vulpinus relying mainly on small epipelagic fish, P. glauca on squid and other organisms associated with the deep scattering layer and I. oxyrinchus being intermediate to A. vulpinus and P. glauca. Therefore, contaminants can be used as an ecological tool to study niche partitioning in species with varying degrees of spatial, temporal or dietary overlap 18,19,46 .
The strong southern California contaminant signal in all three species (i.e. DDX proportion) indicates the relative importance of the CCLME to these species' forage bases and the persistence of DDX in liver tissue. A. vulpinus are considered to be more coastally oriented 32,47 than either P. glauca or I. oxyrinchus. Therefore, it was not surprising that A. vulpinus had greater contaminant contributions from DDXs than both P. glauca and I. oxyrinchus. Note that details on ontogenetic-related changes in DDX signatures can be found in Lyons et al. 33 . A similar patterns relating to DDXs with respect to California coastal residency was found in pinnipeds sampled from southern California 48 , with more coastally associated species having higher DDX proportions and more pelagic species having lower DDX proportions. In comparison, P. glauca and I. oxyrinchus are known to spend protracted periods offshore in open waters 12 . However, despite the time spent offshore, DDXs still contributed to a large proportion of their total contaminant burden (>60% on average), suggesting that prey items with strong southern California influences are prominent in the diets of I. oxyrinchus, P. glauca, and A. vulpinus. By contrast, DDX signatures in Lamna ditropis (Salmon Sharks), another eastern North Pacific shark species that also utilizes offshore waters and occasionally is found in the CCLME, were much lower than that in P. glauca or I. oxyrinchus, with DDXs contributing to only 49% of the total load 34 .
Not only does geographic distribution and diet influence contaminant signatures, but these factors will also influence total contaminant concentrations, which also were distinguishable among species. There was a clear hierarchy in contaminant concentrations among the three species, with I. oxyrinchus having the highest, followed by A. vulpinus and then P. glauca. It is well known that the trophic positions of prey plays a significant role in predator concentrations, with higher trophic level predators having higher accumulated concentrations 49,50 . Of the three species, I. oxyrinchus has the highest trophic level 16 . Considering the relationship between contaminant concentrations and trophic positioning 49,51 , we confirm previous findings that contaminant concentration can be used to make inferences about ecology.
In addition to foraging ecology and distribution, physiological factors are also expected to influence contaminant accumulation 52,53 ; these can, at times, be less emphasized compared to trophic position or trophic linkages in influencing contaminant concentrations among organisms 54,55 . The physiological influence on contaminant concentration is particularly exemplified in the separation between I. oxyrinchus and P. glauca, despite overlap in their diets. Although stable isotope analysis has confirmed that juvenile P. glauca forage at a lower trophic level than juvenile I. oxyrinchus 56 , adults of both species have been documented to feed on marine mammals 16 . Therefore, while diets overlap, the magnitudes greater contaminant concentrations in I. oxyrinchus, suggests that other, physiological factors also have an influence on contaminant accumulation. For example, metabolic rates are higher in regionally endothermic I. oxyrinchus than in ectothermic P. glauca 57 , which will likely influence contaminant accumulation. In contrast, while I. oxyrinchus and A. vulpinus are more physiologically similar (both