|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: NT MA. Performed the experiments: NT MA. Analyzed the data: NT MA. Contributed reagents/materials/analysis tools: NT MA. Wrote the paper: NT MA.
Large-scale patterns or trends in species diversity have long interested ecologists. The classic pattern is for diversity (e.g., species richness) to decrease with increasing latitude. Taxonomic distinctness is a diversity measure based on the relatedness of the species within a sample. Here we examined patterns of taxonomic distinctness in relation to latitude (ca. 32–48 °N) and depth (ca. 50–1220 m) for demersal fishes on the continental shelf and slope of the US Pacific coast.
Both average taxonomic distinctness (AvTD) and variation in taxonomic distinctness (VarTD) changed with latitude and depth. AvTD was highest at approximately 500 m and lowest at around 200 m bottom depth. Latitudinal trends in AvTD were somewhat weaker and were depth-specific. AvTD increased with latitude on the shelf (50–150 m) but tended to decrease with latitude at deeper depths. Variation in taxonomic distinctness (VarTD) was highest around 300 m. As with AvTD, latitudinal trends in VarTD were depth-specific. On the shelf (50–150 m), VarTD increased with latitude, while in deeper areas the patterns were more complex. Closer inspection of the data showed that the number and distribution of species within the class Chondrichthyes were the primary drivers of the overall patterns seen in AvTD and VarTD, while the relatedness and distribution of species in the order Scorpaeniformes appeared to cause the relatively low observed values of AvTD at around 200 m.
These trends contrast to some extent the patterns seen in earlier studies for species richness and evenness in demersal fishes along this coast and add to our understanding of diversity of the demersal fishes of the California Current.
Species are distributed neither uniformly nor randomly across the globe. Large-scale patterns in biodiversity have interested ecologists since at least the time of Wallace and Darwin , . Understanding these large-scale patterns helps us to develop hypotheses regarding how communities and ecosystems are organized on both ecological and evolutionary time scales , . For example, documenting patterns in biodiversity might also be important due to its potential relationship with ecosystem function –. Biodiversity is also a key concept for conservation and management, and is fundamental for ecosystem-based approaches , , . Conserving biodiversity through the protection of species richness is often an implicit or explicit goal of many management and conservation strategies , –.
One of the best known large-scale patterns of biodiversity is the decline in species richness from the equator towards the poles , –. This pattern has been observed in a range of taxa in many different environments—from terrestrial plants  to ants  to deep-water invertebrates  and marine fishes , , and numerous hypotheses have been proposed to explain the latitudinal gradient including: geographic area, evolutionary speed, geometric constraints and productivity . In marine systems species richness is also related to depth, typically declining with depth, although it may be highest at intermediate depths –. Additionally, the relationship between biodiversity and large-scale geological features is important for the understanding of biogeography, which can help in designating management areas.
Biodiversity is, however, a complex concept and is much more intricate than just the total number of species in a given area , . Workers have also examined patterns in species evenness , –, but standard measures of either richness or evenness treat all species as equivalent in value in their contribution to diversity. That is, two species of fish are considered as diverse as one fish and one flatworm. Increasingly, the taxonomic relationships among species have been used to describe another dimension of biodiversity, allowing diversity to be considered within the context of deeper potential functional or evolutionary lineages –.
Taxonomic distinctness quantifies diversity as the relatedness of the species within a sample, based on the distances between species in a classification tree . Average taxonomic distinctness (Δ+ or AvTD) is the mean of all species-to-species distances through the tree for all pairs of species within a sample, and represents the taxonomic breadth of the sample. The variation in taxonomic distinctness (Λ+ or VarTD) is the variation in branch lengths among all pairs of species (it is not the variance of AvTD among samples), and is a measure of the irregularities and divergences in the distribution of branch lengths within a sample. Both indices are appealing because they are based on presence/absence data, and unlike many biodiversity measures, neither is affected by the number of species or the sampling effort , .
In the marine environment, taxonomic distinctness has been used as a tool to examine environmental degradation like the effects of trawling , fishing in general (through marine reserve status, pollution , and other anthropogenic impacts , ). In addition, general trends in taxonomic distinctness with latitude and depth have been examined for some marine invertebrate taxa – and demersal fishes , . Previous studies of taxonomic distinctness in demersal fishes have been quite focused, limited to a depth range of <570 m and done within a single region spanning less than 1.5° of latitude , ). Even for studies of taxonomic distinctness for invertebrates spanning much larger latitudinal ranges, there have been limitations to inferences, due to incomplete data structures in the samples available for analysis (i.e., not all depth strata were sampled at all latitudes of interest, e.g., , ).
Here, we describe quantitatively how AvTD and VarTD varied with depth and latitude for demersal fishes on the continental shelf and slope of the western U.S.A. We also identified individual taxa that were primarily responsible for the observed overall patterns. This work is unprecedented for marine demersal fishes, covering a relatively broad range across temperate latitudes (ca. 32–48° N) and depth (ca. 50–1220 m), as well as having sufficient replication to allow investigation of potential interactions between these two gradients in their effects on taxonomic distinctness in demersal fishes.
We used data from the 1999–2002 Pacific West Coast Upper Continental Slope Trawl Survey  and the 2003 US West Coast Bottom Trawl Survey  to analyze trends in taxonomic distinctness with depth and latitude. The 1999–2002 survey was limited to the continental slope (184–1280 m) while the 2003 survey was expanded to included portions of the shelf (55–183 m). The trawl survey extends from 48°10′N to 32°30′N (Fig. 1). The trawls were carried out using Aberdeen style nets with a small mesh (5-cm stretched measure or less) liner in the cod-end. Trawl duration was approximately 15 minutes at 2.2 knots. Bottom contact and acoustic instruments were attached to the nets to record aspects of mechanical performance as well as gear depth. Catches were sorted to species level or closest taxonomic level. See Keller et al.  for details. We analyzed data consisting of only those taxa identified to species.
Although the trawl survey targets demersal species (fish typically associated with features of the ocean bottom), pelagic species are caught in the trawls as well. Preliminary analyses showed little difference between the results obtained with the full data set (all species) and those obtained with a reduced data set containing only demersal fishes. Therefore, we chose to analyze the full data set. However, the trawl survey does not sample complex, rocky habitat , , and the results of our analyses and our conclusions are limited to ‘trawlable habitat’.
We examined two measures of diversity: average taxonomic distinctness (AvTD or Δ+) and variation in taxonomic distinctness (VarTD or Λ+) . AvTD is a measure of the taxonomic breadth of a sample. It utilizes presence/absence data and is calculated based on the taxonomic distance through a classification tree between every pair of species within a sample (here a trawl):
where is the branch length between species pairs and S is the number of observed species in the sample. Here we used a standard Linnaean classification with 19 taxonomic levels [Table 1, 46]. We used the simple linear scaling  where the maximum distance through the tree is set at ω=100. When branch lengths are un-weighted, the step between each taxonomic level in the tree is considered to be equal. However, all taxa were not defined to the same level of detail. For instance, some groups had defined tribes or subfamilies while others did not. Therefore, we weighted the branch lengths by the number of species having a definition at that level—essentially the amount of information in that taxonomic level. For example, 27% of species were defined at the level of tribe and 63% at the level of subfamily, but all species had definitions at the level of family (Table 1). AvTD is independent of sample size and the number of species in a sample , .
VarTD is the variation in taxonomic distances between each pair of species about the AvTD value for that sample :
where ω is the branch length between pairs of species, S is the number of species observed in that sample and Δ+ is the average taxonomic distinctness for the sample as defined above. VarTD is also independent of sample size, number of species and the value of AvTD within a sample unit.
It is useful to examine both AvTD and VarTD because they capture independent aspects of the taxonomic diversity in a data set. Specifically, two taxonomic hierarchies can produce the same AvTD with different VarTD . For example, a species list containing several different orders each with one genus and species, but also with some genera having many species, would give a high VarTD compared to a list (of equivalent AvTD) in which most species were from different families but within the same order . Finally, the Chondrichthyes exhibited a strong influence on both AvTD and VarTD (see Results). Therefore, we also calculated the taxonomic distinctness measures for the Actinopterygii alone to examine taxonomic distinctness of this major taxon. Values of AvTD and VarTD were calculated for each trawl using PRIMER v6 .
We used generalized additive models [GAMs], [ 49,50] to examine how AvTD and VarTD varied with latitude and depth because preliminary diagnostic plots suggested that relationships were not linear. For each measure we fit a two factor GAM (with identity link and Gaussian error distribution):
where x1 was depth, x2 was latitude, β0 was the intercept parameter, and εi were random normal errors with zero mean and a common variance The response variable y was either AvTD or VarTD. The smoothing functions f1 and f2 were thin plate regression splines . Because depth and latitude were measured on different scales, we used a tensor product smooth (f12) of thin plate regression splines for the interaction term . The optimal level of smoothing was chosen with general cross validation. Data were fit in R v2.10.0 using the package ‘mgcv’ , . While the data set spans 55–1280 m, data were sparse in some of the deepest areas at some latitudes, so we limited the analyses to data collected from depths <1220 m to assure relatively even coverage across depths and latitudes. Diagnostic plots of residuals showed reasonable symmetry, but some deviation from normality in terms of kurtosis (Figure S1 & Figure S2). Although GAMs will be fairly robust to this kind of deviation, conclusions should nevertheless be taken with some caution.
To more specifically examine how AvTD and VarTD were related to latitude at particular depths (and vice versa), and to more fully investigate significant depth x latitude interactions (see Results), we conducted a second round of GAM (identity link, Gaussian distribution) analyses in which we binned data either by depth or by latitude zones. Here,
where β0 was the intercept parameter, x was depth within a specific latitude bin (or latitude within a specific depth bin), f was a thin plate regression spline and εi were random normal errors with zero mean and a common variance.
We established the following five depth zones after examining the first analyses for AvTD (see Results): 50–150 m, 151–300 m, 301–600 m, 600–900 m, 900–1220 m. We binned data into four latitude zones based on approximate location of large-scale geographic features along the coastline: South of Point Conception, Point Conception to Cape Mendocino, Cape Mendocino to Cape Blanco, and North of Cape Blanco (Fig. 1).
When the degrees of freedom for the smoothed term were near their minimum value, we examined whether the trend was described sufficiently by a linear regression using Akaike's information criterion (AIC) and analysis of deviance between the two models.
To identify which taxa influenced patterns of taxonomic distinctness the most, we used two complimentary steps: analysis of taxonomic trees and exclusion of specific taxa from the calculation of both AvTD and VarTD. First, we compared taxonomic trees across depth zones for the area between Point Conception and Cape Mendocino to examine how the branching patterns differed. We chose this region because here the r2 for depth was strongest for AvTD and high for VarTD (see results). We used the same depth strata as in previous analyses. We calculated the frequency of occurrence for each species in each depth x latitude bin and constructed trees for each bin that comprised those species which occurred in at least 50% of the hauls for that bin. These ‘50% trees’ do not completely describe the classification tree for the five depth zones, but they do provide a qualitative guide to the type of transitions in relatedness among taxa that occurs between depth strata.
Second, we examined the effect of excluding specific taxa from the calculation of AvTD and VarTD . After visually examining the taxonomic trees, we identified several taxa of interest. We then re-calculated AvTD (or VarTD) excluding the individual taxon of interest and subtracted these values from the AvTD for the full data set (e.g., AvTD - AvTDno taxon 1=AvTDdiff). While one could examine only AvTDno taxon 1, AvTDdiff givers a better visual interpretation of the effect of that taxon on the overall AvTD. For each taxon of interest, we analysed the relationship between AvTDdiff vs. depth and latitude using GAMs as described above. In these analyses, negative values of AvTDdiff indicate that the presence of that taxon decreased AvTD in the original analysis of all data because AvTD increased when it was removed, while positive values indicate the opposite.
We examined 1948 trawls from surveys conducted between 1999 and 2003. These trawls contained 243 species from 75 families, 30 orders and four classes. The most speciose families were Scorpaenidae (rockfishes, 48 spp.), Pleuronectidae (righteye flounders, 16 spp.), Zoarcidae (eelpouts, 14 spp.), Cottidae (sculpins, 15 spp.), Liparidae (snailfishes, 11 spp.) and Rajidae (skates, 8 spp.) The most speciose orders were Scorpaeniformes (mail cheeked fishes including rockfishes and scorpion fishes, 90 spp.), Perciformes (perches, 34 spp.), Pleuronectiformes (flatfishes, 21 spp.), Gadiiformes (cods, 13 spp.), Argentiniformes (marine smelts, 11 spp.), and Stomiiformes (dragonfishes, 11 spp.). By class the majority of fishes were Actinopterygiians (ray-finned fishes, 218 spp.) or Chondrichthyians (cartilaginous fishes, 23 spp.) with one Myxini (hagfishes) and one Petromyzontida (lampreys).
Average taxonomic distinctness (AvTD) varied significantly with depth (F7.61=20.73, p<0.001) and latitude (F7.84=8.05, p<0.001), and there was also a significant depth x latitude interaction (F13.38=11.26, p<0.001)(Fig. 2a,b). The GAM explained 24% of the variation in AvTD (r2=0.24, n=1948). AvTD was highest around 500 m, especially in the region around 35°N. AvTD was lowest around 200 m and intermediate at depths over approximately 800 m. AvTD in the shallowest regions was similar to that found in the deepest depths. The estimated intercept () was 49.34 (±0.09 s.e.) indicating that, on average, species were related between the level of Superorder and Subdivision (Table 1).
Variation in taxonomic distinctness (VarTD) was also significantly affected by depth (F5.25=65.24, p<0.001), latitude (F7.56=3.14, p=0.002), and their interaction (F12.16=10.78, p<0.001)(Fig. 2c,d). The GAM explained 63% of the variation in VarTD (r2=0.63, n=1948). VarTD generally increased from the shallower areas, reaching a peak around 200–300 m and then decreasing sharply with increasing depth. Overall, there was a weak correlation between AvTD and VarTD (n=1946, r=0.27, p<0.001).
To further investigate depth x latitude interactions, we established five depth zones based primarily on the AvTD patterns in Figure 2: 50–150 m, 151–300 m, 301–600 m, 600–900 m, 900–1220 m. The shallowest depth bin was limited to the shelf and was set to encompass depths at which AvTD appeared to increase with increasing latitude, in contrast to the rest of the data. The 151–300 m zone brackets an area where AvTD is lowest (around 200 m). It includes the deepest portions of the shelf and the shallow areas of the slope. AvTD was highest at depths between 301–600 m and decreased between 601–900 m. It then appeared to remain fairly stable at depth. For consistency, we utilized these same depth bins for VarTD.
When binned by depth, latitude explained significant but small amounts of variation in AvTD at all depths (Fig. 3). In the shallowest depth zone (50–150 m), a linear model was sufficient (analysis of deviance, p=0.12) to describe the relation between AvTD and latitude. AvTD increased with latitude but the explained variance was just over three percent (linear regression, F1, 169=6.26, p=0.017, r2=0.033). Between 151–300 m the relationship between AvTD and latitude was more complex with peaks in AvTD around 35°N and 40°N, and with lows around 32°N, 37°N and 45–46°N, but with no overall trend. The effect of latitude on AvTD was strongest in the 301–600 m depth zone with latitude explaining 25% of the variance in AvTD. AvTD peaked just north of 35°N and then declined overall as latitude increased, with a number of peaks and troughs. At 601–900 m a linear trend was sufficient (analysis of deviance, p=0.13) to describe a decrease in AvTD with latitude but explained only 7% of the variation (linear regression, F1, 381=31.91, p<0.001, r2=0.07). In the deepest area (901–1220 m), AvTD decreased rapidly between 32°N and 35°N. It then increased until approximately 37–38°N after which it decreased with increasing latitude.
VarTD was also only weakly related to latitude within the five depth strata (Fig. 4). At 50–150 m, VarTD increased with latitude, explaining 19% of the variance, with peaks at approximately 34°N, 39°N and 45°N. In the 151–300 m depth zone, latitude explained only 8% of the variance in VarTD. VarTD increased from 32°N to 40°N and then decreased to the north. At 301–600 m VarTD increased between 32°N and approximately 37°N and then decreased to a low at 45°N. At 601–900 there was a similar peak in VarTD around 36–37°N. At 900–1220 m, VarTD was lowest between 34–35°N and peaked around 38°N, but there was no clear overall trend with latitude.
The relationship between either AvTD or VarTD and depth was stronger (higher r2) than their respective relationships with latitude (Fig. 5). The explained variance for AvTD was highest in the Point Conception to Mendocino latitude bin and fairly low elsewhere. For VarTD, the variance explained by depth was much higher, being greater than 62% in three of the four latitude bins. There was some variation, as expected from the significant latitude x depth interaction terms in the main GAMs, with peak values for both AvTD and VarTD shifting among latitude bins.
The taxonomic trees provided here, consisting of those species found in at least 50% of the trawls within a given depth bin, comprised 12–14 species per depth bin for the region between Point Conception and Cape Mendocino. Thirty-four species were found at least 50% of the time in at least one depth bin (Fig. 6, Table 2). Several taxa (Pleuronectiformes, Scorpaeniformes and Chondrichthyes, in particular) showed differences in their branching patterns, which may help explain the variation in AvTD among depths.
At 50–150 m, the high diversity of the Pleuronectidae (six species in six genera) appears to be counter-balanced by the relative paucity of species in other families or orders. The result is intermediate AvTD and high VarTD. The low AvTD around 200 m (here the 151–300 m depth stratum) appears to be due to a combination of the Pleuronectidae (four species in four genera) and the Scorpaeniformes (four species in two families), and, in particular, the genus Sebastes, which had three species. AvTD was highest between approximately 301–600 m. Some of the increase appears due to the addition of an entire order (Carcharhiniformes) to the Chondrichthyes at these depths. Additionally, there were only two Pleuronectids, and the Scorpaeniformes added families and genera while losing species within the Sebastes.
The primary change between the shallower areas and the deeper areas is the reduction of taxa within the Chondrichthyes in the zones deeper than 600 m. This likely explains much of the decrease in AvTD with depth as the number of species in different classes dropped. Additionally, the Scorpaeniformes gained a species (Sebastolobus altivelis). Among the deeper depth zones (those greater than 600 m), the Scorpaeniformes, Pluronectiformes, and Argentiniformes showed no change in branching structure. The continued decline in AvTD with depth appears to be due to the addition of a Morid, Antimora microlepis to the Gadiformes, and the addition of Bothrocara brunneum to the Zoarcidae.
The relationship between tree structure and VarTD is easy to discern in the current example. VarTD was highest in the shallower depth strata around 200 m in particular (Fig. 2c,d & 5e,f,g,h). At these depths, the trees contained a combination of short branches and long branches. There were high numbers of related species (four Pleuronectids and three Sebastes), but three of the four other orders contained only one species. Moreover four of the thirteen species were in the Chrondrichthyes and thus related to the remaining species at the level of class. At deeper depths, species were more evenly distributed among taxa with only one to two species per family (so fewer short path-lengths) and fewer Chondrichthyes (so relatively fewer long path-lengths).
Based on the analysis of the taxonomic trees, we chose to examine the effect on AvTD of removing five taxa: Scorpaeniformes, Pleuronectiformes, Gadiformes, Perciformes and the Chondrichthyes. AvTDdiff and VarTDdiff for all five taxa showed significant depth x latitude interactions (Fig. 7 & 8, GAM, p<0.05 for all). For brevity, we do not present full GAM results for each taxon but give only the r2 for each model.
AvTDdiff for Scorpaeniformes (Fig. 7a, r2=0.37) was negative across all depths and latitudes indicating that this order generally lowered the AvTD for the full dataset. This was primarily due to there being many species in the genus Sebastes, as well as Scorpaeniformes being the most speciose order. AvTDdiff was lowest (ca. −10) at approximately 200 m across all latitudes, which coincides with the relatively high diversity of Sebastes at these depths seen in the tree analysis and the low overall AvTD at the same depth.
In regions deeper than approximately 200 m, removal of Pleuronectiformes (Fig. 7b, r2=0.50) did not cause large changes in AvTD, which was negative in value but close to zero. However at depths shallower than 200 m, AvTDdiff was close to −10 indicating that this order lowered AvTD in the shallower areas. This pattern corresponds directly with the diversity of the flatfishes seen in the tree structure, with four to seven species found to occur in the shallower two zones but only two species in the deeper areas. AvTDdiff for the Pleuronectiformes also appeared to decrease with latitude in the shallows.
AvTDdiff for both Gadiformes (Fig. 7c, r2=0.08) and Perciformes (Fig. 7d, r2=0.35) also showed significant relationships with latitude and depth. The Gadiformes raised AvTD by 15 points or more overall, but showed only small differences in relation to depth and latitude. The Perciformes also showed only minor variation with depth and latitude, generally having only a small effect on AvTD (with AvTDdiff values near zero).
Removal of the class Chondrichthyes from the calculation of AvTD resulted in the most dramatic effects (Fig. 7e, r2=0.39). AvTDdiff for Chondrichthyes was large and positive across all depths ranging from 20 to just under 40 and indicating that this taxon raised AvTD substantially in all areas. More importantly, AvTDdiff for the Chondrichthyes matched the overall pattern of AvTD from the full data set indicating that this taxon drove the general pattern of taxonomic distinctness. It was highest between about 300–600 m and declined with latitude in this depth region. It also increased with increasing latitude in the shallower areas and was lower and relatively flat at depth. These results match the tree branching patterns with the highest Chondrichthyan diversity seen in the 300–600 m depth areas and the lowest seen at depth.
AvTD calculated for only the Actinopterygii shows a substantially different pattern from the overall pattern obtained using the full data set (Fig. 7f). AvTDActin. is still lowest around 200 m depth but is higher at depths greater than 600 m than in shallower regions. The depth x latitude interaction was significant (GAM, F13.1=4.40, p<0.001), as were both main effects of depth (F6.9=57.50, p<0.001) and latitude (F4.5=2.82, p=0.02). The GAM showed that depth and latitude explained 64% of the variation in AvTDActin. (n=1948, r2=0.64).
For VarTD, removal of Scorpaeniformes (Fig. 8a, r2=0.20), Pleuronectiformes (Fig. 8b, r2=0.50), Gadiformes (Fig. 8c, r2=0.27), and Perciformes (Fig. 8d, r2=0.25) resulted in negative VarTDdiff values at most depths and latitudes, indicating that these taxa generally lowered VarTD when included in the analysis. The Scorpaeniformes did increase VarTD (positive VarTDdiff values) at approximately 200 m in the more southern areas, as did the Pleuronectiformes in the shallowest areas at higher latitudes. VarTDdiff for the exclusion of Chondrichthyes (r2=0.60), however, was strongly positive and generally matched the overall pattern of VarTD.
VarTD for the Actinopterygii was related to depth (F7.12=39.0, p<0.01), latitude (F7.74=4.0, p<0.01) and there was a significant depth x latitude interaction (F11.3=9.2, p<0.01) (Fig. 8f). The model explained 53% of the variation in VarTD for the Actinopterygii (r2=0.53, n=1948). As for VarTD calculated using all taxa, VarTD for the Actinopterygii was highest around 200 m and decreased with depth.
Large-scale gradients in taxonomic distinctness have not been extensively studied in the marine environment. Here, taxonomic distinctness (AvTD and VarTD) varied with both latitude and depth, but the latitude patterns were weak overall. This result is similar to other work in marine systems where trends of AvTD within phyla tend to be weak in relation to latitude (when controlling for depth) but variation with depth is stronger. For both latitude and depth, the two taxonomic distinctness measures showed patterns that differed from more traditional diversity measures, specifically species richness and evenness, for the same fish assemblage. Thus, AvTD and VarTD provide a perspective on diversity that differs somewhat from that obtained when analyzing richness or evenness alone. These distinctness measures provide important, complimentary information that should be of interest to biologists and managers. An important caveat to our analyses is that they are limited to trawl-sampled demersal fishes on primarily soft bottoms where trawl surveys can be conducted. Patterns on hard, complex substrata (e.g., rocky reefs) or with fishes not susceptible to capture by trawls may differ.
The relationship between taxonomic distinctness and either latitude or depth for various marine taxa has not been extensively studied. Where it has been investigated, the relationship with latitude is generally weak and relationships with depth are stronger , , , , although this weak latitude relationship is not always obvious unless one controls for depth . For annelids and crustaceans on the Norwegian continental shelf, and for macrobenthos on the European continental shelf, AvTD increases with increasing latitude but the explained variation is small , . However north-east Atlantic demersal fishes do show geographic variation in taxonomic distinctness between the western waters of the United Kingdom and the southern North Sea, with Elasmobranchs contributing substantially to the patterns . At very large spatial scales, AvTD for pelagic copepods shows low variability between 40°N and 40°S but then becomes highly variable as it declines towards the poles .
Latitudinal trends in both AvTD and VarTD were present in west coast demersal fishes but were generally very weak. On the shelf (here the 50–150 m depth stratum) both AvTD and VarTD increased with latitude, as seen with macrobenthos on the Norwegian and European shelves. While the relationship was very weak for AvTD and only moderately strong for VarTD, in both cases these trends are the reverse of the typical trend of diversity (in the sense of species richness), which decreases with increasing latitude. The only moderately strong latitudinal trend was for AvTD between 301–600 m, where latitude explained 25% of the variation. In the deeper areas, AvTD and Var TD tended to decrease with increasing latitude following the more common diversity-latitude pattern.
Latitudinal trends in taxonomic distinctness for demersal fishes on the continental slope (here deeper than 151 m, although this does include a small portion of the shelf) differed from those reported for species richness and evenness , especially in the more shallow areas of the slope (Tolimieri  did not calculate richness or evenness values for the shelf, so we cannot make comparisons with results obtained from the 50–150 m depth bin). Richness and evenness were both positively correlated with latitude between 150–349 m and richness between 350–549 m . However, AvTD varied but showed no overall trend with latitude between 151–300 m, and it was negatively correlated with latitude between 301–600 m. In deeper areas, all metrics tended to show a decline with latitude.
Relationships between each of the two distinctness measures and depth were stronger (higher r2) than they were with latitude — especially for VarTD. Perhaps the most interesting pattern relates to the 151–300 m depth stratum. Richness and evenness for slope demersal fishes  at similar depths were both relatively high, as was VarTD in the present study. However AvTD was relatively low in this depth stratum, giving a different perspective on the nature of diversity. Combining the different diversity measures, we can conclude that there were many, closely related species at these depths, in particular within the family Scorpaenidae.
In deeper areas, richness, evenness and AvTD showed a general decrease with depth from high values around 301–600 m (400–500 m stratum in Tolimieri ), while species density and VarTD peaked and declined at shallower depths, generally around 200 m. Intermediate to low diversity at depth was due to several factors: fewer species, a less even representation of those species in terms of their relative abundances, higher relatedness among those species and a more consistent branching pattern in the classification tree — the latter two features being due to there being fewer Chondrichthyes. However, within the Actinopterygii, AvTD was actually higher at depth than in more shallow areas. High Actinopterygiian AvTD at depth was due to these species being more evenly distributed among orders. This strong effect of cartilaginous fishes on taxonomic distinctness measures was also seen by Rogers et al.  for demersal fishes in the waters from the western United Kingdom to the southern North Sea. Mechanistically or perhaps mathematically, the long branch lengths that the Chondrichthyes add (due to being related to most other taxa only at the level of class) and the variability in their frequency of occurrence, especially with depth, is the root of their strong effect on taxonomic distinctness patterns.
While patterns with depth were stronger than with latitude, in all cases the depth x latitude interaction was significant for AvTD and VarTD. This result is similar to that seen for both assemblage structure  and species richness and evenness trends  in west coast demersal fishes. For example, northerly trawls had deeper/colder-water assemblages at a given depth than did sites farther to the south , although this pattern disappeared at deeper areas. Bottom temperatures decrease with increasing latitude, especially in the shallower areas of the slope, and are a potential cause of the depth x latitude interaction.
The low explained variation for AvTD (~25%) may result from a general lack of any strong pattern at this spatial scale or from missing explanatory variables, with substratum type being an obvious potential factor. Diversity can differ among habitat types and more complex combinations of multiple habitat types often leads to higher diversity in terms of species richness . While taxonomic distinctness may differ naturally among habitats, the effect does not seem to be strong, at least among soft bottom communities (e.g., mud versus sand). For example, sediment grain size explained little variation in AvTD for annelids, crustaceans and molluscs on the Norwegian Shelf. The greatest variance explained by sediment grain size was around 2.1% for crustaceans. Similarly, unless sites are degraded somewhat, AvTD for free-living nematodes  and for megafauna, macroinfauna and nematodes  tends to be similar to expectations based on the regional species pool.
The data analyzed here did not include samples from complex, rocky habitat where towing nets is difficult or impossible. As such, the largest potential difference in habitat type (rocky versus soft sediment) did not influence the analyses. The importance of rockfish to AvTD, especially around 200 m, suggests that habitat might have a strong effect on AvTD if complex, rocky habitats were included in the analyses. Rockfish are very diverse, with 72 species in the northeast pacific, and they are most common in rocky habitats . One might expect, therefore, that rocky habitats dominated by many closely related species (with lots of rockfish) would have low AvTD compared to other areas with an assemblage more diverse at higher taxonomic levels.
There is variation in oceanic habitat (upwelling intensity, chlorophyll-a, temperature, salinity, etc.) along the west coast that correlates with the assemblage structure of demersal fishes . Similarly, assemblage structure and diversity patterns of oceanic copepods appear to be related to productivity regimes . These oceanic habitats might help to explain some of the large-scale patterns seen in our data. For example the peak in AvTD in the 50–150 m depth bin between approximately 43–46 °N corresponds with oceanic habitat characterized by consistent upwelling . However large-scale changes in oceanic habitats alone would not be enough to explain the rather large spread of residuals (unexplained variation) around the trend lines in our analyses.
On the west coast of the US there are a number of geological features that may be important as biogeographic boundaries or transition zones. Point Conception (~34.5°N) has long been thought to be a boundary for fishes . However, more recent phylogeographic  and range end-point analyses  have suggested that the Los Angeles Region and Monterey Bay Region, both of which have submerged canyons, are the more likely boundaries for fishes , molluscs ,  and marine algae . Point Conception is better interpreted as a transition zone as opposed to a barrier , and there is some evidence that this transition zone may have shifted over time . For slope demersal fishes, the assemblage structure , species richness and evenness  all showed some relation to geographic features, although there is variability among depth zones. In terms of assemblage structure, Cape Mendocino was the most consistent region of change in the assemblage structure of the 26 most abundant slope demersal fishes. Richness tended to show highs around Point Conception and lows near Monterey Bay, while evenness tended to be low around Cape Mendocino.
While the patterns vary somewhat among the depth zones, AvTD and VarTD tended to show peaks around Point Conception and Cape Mendocino and showed both lows and highs around the Monterey Bay area (ca. 37–38°N). There was some variation with depth, as Point Conception tended only to be important in the more shallow depth strata. The high diversity in most metrics around Point Conception makes sense if it is a transition zone between northern and southern fauna –. Given the general weakness of latitudinal trends in taxonomic distinctness, these observations should not be overly stressed.
Conservation and management efforts often focus on maintaining biodiversity, specifically species richness , , . At one level this focus is an attempt to protect as many species as possible with limited resources . At another level, the potential relationship between ecosystem function and biodiversity also makes the maintenance of richness a priority , –.
Diversity is, however, a complex idea, and it can be measured in a number of complimentary ways . While important and essential to our understanding of communities, measures like species richness and evenness do not recognize the taxonomic relationships among species, and treat all species as equivalent. More recently a number of metrics have been developed that do take into consideration the evolutionary relationships, including taxonomic diversity , phylogenetic diversity ,  and the taxonomic distinctness measures (AvTD and VarTD) used here. The fact that these taxonomic distinctness measures are statistically independent of either the sampling intensity or richness of samples ,  makes them especially attractive tools for investigating structural biodiversity.
Taxonomic distinctness often provides complimentary or contradictory information to the more typical richness analyses. For example, species richness was lower in the Columbretes Island Marine Reserve than in reference sites. However, taxonomic distinctness was higher in the reserve . For demersal fishes on the slope in the present study, richness and evenness were highest in the shallow portion of the continental slope, but this area had low AvTD due to the high relatedness of species (mostly rockfish or flatfish) in this depth stratum. In fact, for the Actinopterygii, diversity measured as AvTD was the opposite of overall richness patterns, being highest at depth.
Focusing on species richness alone does not recognize diversity among deeper evolutionary lineages. The presence of several distinct lineages at higher taxonomic levels may be important to provide resources for future evolutionary innovation , , especially since extinction and speciation rates differ among taxa. Measures of taxonomic distinctness may also represent functional diversity to some extent . Not evaluating diversity at higher taxonomic levels fails to recognize potentially important patterns of diversity. For example, there are fewer plant species in Malesia—a region including Malaysia to Papua New Guinea and the Solomon Islands—but there are more plant families in this region than in the rest of the neotropics. This is important since the number of families may provide a better measure of both functional diversity and the evolutionary potential of the assemblage . This pattern would likely be reflected by higher AvTD values for Malesia than for other areas, as more species are related at higher taxonomic levels. We do not mean to suggest that metrics like species richness should be ignored, but a more holistic analysis of assemblage structure, including measures which incorporate the degree of relatedness among species, is necessary to fully understand patterns of diversity.
Residual plots for Generalized Additive Model (GAM) examining variation in AvTD versus depth and latitude. In the GAM, yi=β0+f1(x1i)+f2(x2i)+f12(x1i×x2i)+εi where yi was AvTD, x1 was depth, x2 was latitude, β0 was the intercept parameter and εi were random normal errors with zero mean and a common variance. The smoothing functions f1 and f2 were thin plate regression splines . Because depth and latitude were measured on different scales, we used a tensor product smooth (f12) of thin plate regression splines for the interaction term . The optimal level of smoothing was chosen with general cross validation. Data were fit in R v2.10.0 using the package ‘mgcv’ , .
(0.20 MB PDF)
Residual plots for Generalized Additive Model (GAM) examining variation in VarTD versus depth and latitude. In the GAM, yi=β0+f1(x1i)+f2(x2i)+f12(x1i×x2i)+εi where yi was VarTD, x1 was depth, x2 was latitude, β0 was the intercept parameter and ε1 were random normal errors with zero mean and a common variance. The smoothing functions f1 and f2 were thin plate regression splines . Because depth and latitude were measured on different scales, we used a tensor product smooth (f12) of thin plate regression splines for the interaction term . The optimal level of smoothing was chosen with general cross validation. Data were fit in R v2.10.0 using the package ‘mgcv’ , .
(0.27 MB PDF)
We thank one anonymous reviewer, P. Taylor, W. Wakefield, M. Ford, J. Cope, and P. Levin, for comments on various drafts. J. Samhouri, T. Francis, K. Andrews, C. Harvey, and G. Williams all provided useful discussion. We especially thank D. Jeter, A. Pettitte, J. Posada, M. Rivera, et al. for all that they do to make life enjoyable, and wish H. Matsui the best of luck. This analysis would not have been possible without the hard work of the trawl survey team.
Competing Interests: The authors have declared that no competing interests exist.
Funding: This research was supported in part by NOAA Fisheries (Northwest Fisheries Science Center). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. M. J. Anderson was partially supported by a Royal Society of New Zealand Marsden Grant, no. MAU0713 and also by a Research Associateship from the US National Research Council, through C. Jordan at NOAA.