Unexpected high vulnerability of functions in wilderness areas: evidence from coral reef fishes.

High species richness is thought to support the delivery of multiple ecosystem functions and services under changing environments. Yet, some species might perform unique functional roles while others are redundant. Thus, the benefits of high species richness in maintaining ecosystem functioning are uncertain if functions have little redundancy, potentially leading to high vulnerability of functions. We studied the natural propensity of assemblages to be functionally buffered against loss prior to fishing activities, using functional trait combinations, in coral reef fish assemblages across unfished wilderness areas of the Indo-Pacific: Chagos Archipelago, New Caledonia and French Polynesia. Fish functional diversity in these wilderness areas is highly vulnerable to fishing, explained by species- and abundance-based redundancy packed into a small combination of traits, leaving most other trait combinations (60%) sensitive to fishing, with no redundancy. Functional vulnerability peaks for mobile and sedentary top predators, and large species in general. Functional vulnerability decreases for certain functional entities in New Caledonia, where overall functional redundancy was higher. Uncovering these baseline patterns of functional vulnerability can offer early warning signals of the damaging effects from fishing, and may serve as baselines to guide precautionary and even proactive conservation actions.

Introduction

Human activities have already induced the collapse of many ecosystems around the world [1-3] and, in combination with climate change, have triggered major reductions in biodiversity globally [3-7]. Beyond the loss of species, there is a growing awareness that the loss of ecological functions may be the most critical consequence of human disturbances on ecosystems [8-12]. This diversity of ecological functions sustains ecosystem services on which humanity depends; such as biomass production [10]. Sustaining ecosystem functions requires both high functional diversity, i.e. a large breadth of ecological functions supported by species [13-17], and high functional redundancy, i.e. a large number of species supporting identical functions in the system [18].

In theory, species richness is thought to maintain a high level of both functional diversity and redundancy, thus ensuring the long-term functioning of ecosystems in a fluctuating environment [19]. Indeed, high species richness should increase the probability of having both species supporting different functions (functional diversity) and many species supporting the same ecological functions (functional redundancy) [17]. Many experiments confirm this theory, for example demonstrating the vulnerability of ecosystem functioning to species loss [11,20]. In natural systems, however, the benefits of high species richness to maintaining ecosystem functioning have recently been challenged by three patterns. First, at the scale of ocean basins (and using species checklists), some functions exhibit over-redundancy, i.e. are supported by a disproportionately high number of species, while others are realized by few or one species only, even in the richest regions [21]. Second, species that support specific or unique traits in ecosystems (no redundancy) tend to be rare owing to their low abundance in ecosystems [22]. Third, the distribution of species richness and abundance among trophic groups is more critical than simply the number of species to maintain ecosystem functioning and services [23]. These patterns demonstrate the importance of preserving both species and abundance within a wide range of functional groups.

Taken together, these results suggest that high levels of species richness and abundance may not insure ecosystems against functional diversity loss as we once hoped, owing to the high vulnerability of some functions that lack redundancy. Yet, this hypothesis has, to our knowledge, never been tested with empirical data in tropical ecosystems with marginal or no exposure to threats, i.e. where species density and abundances should be close to natural baselines. Assessing the vulnerability of ecological functions to threats in such scenarios would reveal the baseline of functional vulnerability, and the extent to which these ecosystems are buffered against even limited local species declines or extinctions.

Protected areas (PAs) are often used to assess ecological baselines against which biodiversity levels in exploited areas are compared [24]. However, recent studies have shown that even PAs cannot be considered as true ecological baselines because anthropogenic disturbances typically started long before these areas were established [2,25-27]. In addition, most of these areas are either too small or embedded in areas influenced by human activities and therefore cannot support the full range of ecological functions [24,28]. As an alternative, wilderness areas, i.e. large areas geographically isolated from humans by natural geographical barriers or with very limited human presence, may provide ecological baselines close to a ‘natural’ status [2,26,28]. Indeed, wilderness areas are traditionally viewed as areas featuring exceptional concentrations of biodiversity and abundance [29], albeit potentially suffering from global changes in a near future [30].

Here, using coral reef ecosystems along a geographical gradient, we propose to test two hypotheses. First disparate wilderness areas tend to host a similar level of functional diversity, uncovering a consistent baseline for the breadth of functions in ecosystems, despite a high turnover in species composition. Second, this common level of functional diversity remains highly vulnerable to species declines or losses owing to a disproportional over-redundancy in some functions and a ‘natural’ lack of redundancy for some critical functions.

Coral reefs are the most diverse marine ecosystems on Earth [31] and support key services for half a billion people, such as food and income [32]. We quantified the baseline vulnerability to fishing of fish functional diversity in coral reef ecosystems across the Indo-Pacific geographical gradient, taking advantage of extensive surveys in French Polynesia, New Caledonia and Chagos. These three wilderness areas all benefit from a high level of isolation from humans [33] and a high level of enforcement owing to the presence of military forces, thereby limiting illegal fishing activities. As the ecological knowledge to assess the functions carried by individual species is limited, using species functional traits to infer functions offers a viable alternative [34]. Here, we assume that species with more diverse combinations of functional traits are more likely to support different functions (e.g. [35-37]).

Material and methods

Study regions

Remote atolls and islands in three regions (Chagos Archipelago, New Caledonia, and French Polynesia) were sampled along the Indo-Pacific biogeographic gradient, encompassing 130° of longitude (16 000 km) (electronic supplementary material, figure S1). None of these atolls and islands are inhabited: the northern Chagos Archipelago (the Great Chagos Bank, Peros Banhos and Salomon Island) is more than 650 km south of the Maldives and personnel at the Diego Garcia atoll navy base are not permitted to the northern Archipelago other than for fishery patrols; isolated atolls and islands in New Caledonia (Entrecastaux Archipelago, Astrolabe Reef, and Beautemps-Beaupré) are located between 300 and 600 km and more than 20 h by boat from the capital Nouméa [25]; two atolls at the southeast end of the Tuamotu Archipelago (Paraoa and Ahunui) are located approximately 950 km from Papeete, the capital of French Polynesia; and the Acteon Group, a cluster of atolls, is located between 200 and 500 km north of Gambier Island, French Polynesia (figure 1).

Figure 1.

Sampled coral reef ecosystems across three regions of the Indo-Pacific. Fishes were identified and counted at the outer reefs of remote atolls (blue stars) in Chagos (79 transects), New Caledonia (18 transects) and French Polynesia (37 transects).

Fish surveys

Fish data for the Pacific Islands were collected on outer reef slopes using Underwater Visual Census (UVC) along 50 m transects. Briefly, this method involved two divers each recording species identity, abundance and body length [38]. Transects in New Caledonia (18 transects) and French Polynesia (37 transects) were truncated at 10 m wide (10 × 50 m strip transects). In the Chagos Archipelago, fishes were surveyed along outer reef slopes using 50 × 5 m strip transects (79 transects). To make surveys comparable among regions, a common area of 500 m2 was obtained by randomly aggregating two Chagos transects (250 m2). Accumulation curves of species richness were performed for each region to test for potential biases owing to survey techniques (electronic supplementary material, figure S1). Species densities and abundances were estimated for 500 m2 transects and averaged for each region. Sharks and rays were removed from the main species list because of their specific traits, some poaching in Chagos, and difficulties in assessing their abundance using UVC [39]. As such, this study focused on 412 fish species belonging to 35 teleost families.

Functional traits and entities

To estimate functional diversity, we used six functional traits related to major fish attributes: (i) maximum body size, (ii) diet, (iii) home range, (iv) position over the reef, (v) activity, and (vi) gregariousness [21,40].

Fish sizes were coded using six ordered categories: 0–7 cm, 7.1–15 cm, 15.1–30 cm, 30.1–50 cm, 50.1–80 cm and more than 80 cm. Diet was characterized based on the main items consumed by each species, which led to seven trophic categories: herbivorous–detritivorous (i.e. fishes feeding on turf or filamentous algae and/or undefined organic material), macroalgal herbivorous (i.e. fishes eating large fleshy algae and/or seagrass), invertivorous targeting sessile invertebrates (i.e. corals, sponges and ascidians), invertivorous targeting mobile invertebrates (i.e. benthic species such as crustaceans, echinoderms), planktivores (i.e. fishes eating small organisms in the water column), piscivorous (including fishes and cephalopods), and omnivorous (i.e. fishes for which both vegetal and animal material are important in their diet) [21,41]. Home range was coded using three ordered categories: sedentary (including territorial species), mobile within a reef, and mobile between reefs. Position in the water column was coded using three ordered categories: benthic, bentho-pelagic, and pelagic. Activity period was coded using three ordered categories: diurnal, both diurnal and nocturnal, and nocturnal. Schooling was coded using five ordered categories: solitary, pairing or living in small (3–20 individuals), medium (20–50 individuals) or large (more than 50 individuals) groups. This functional trait's database was built using information about the ecology of adult life-stages available in the literature and according to observations made in the Indo-Pacific by the survey team [40,42,43].

More detailed descriptions linking these traits to ecological processes can be found in the electronic supplementary material of published articles [21,28]. Because all traits were coded using categories, we defined functional entities (FEs) as groups of species sharing the same trait categories. In total, 412 fish species were clustered into 157 different FEs [21,28].

Fish functional diversity

In order to compare the level of functional richness (FRic) among the wilderness areas, we measured the FRic of fishes for each region defined as the volume inside the convex hull occupied by species within a functional space [44,45]. To build a functional space, we calculated pairwise functional distances between species pairs based on the six functional traits using the Gower metric, which allows mixing different types of variables while giving them equal weight [46]. A principal coordinates analysis (PCoA) was performed on this distance matrix to build a multidimensional functional space. We retained the first four principal axes (PCs), which faithfully represent the Gower distance between species [44,45,47].

Taxonomic and functional β-diversity

In order to determine if wilderness areas host different species and functional compositions across the Indo-Pacific, we assessed taxonomic and functional β-diversity among regions. We used the β-diversity partitioning framework based on the Jaccard dissimilarity index [48,49]. Taxonomic β-diversity (βjac) equals (equation (2.1))where, a is the number of species in both regions A and B, b is the number of species present in region A but not in region B and c is the number of species present in region B but not in region A.

To distinguish between the species replacement versus nestedness components of β-diversity, we decomposed the pairwise Jaccard dissimilarity index (equation (2.1)) into two additive components [48,49]. The replacement component of the Jaccard dissimilarity index (βjtu, equation (2.2)) describes species replacement without the influence of richness difference between regions. This index is formulated as follows:where a, b and c are the same as in Equation 1.

The nestedness component of the Jaccard dissimilarity index (βjne, equation (2.3)) is the difference between βjac and βjtu. This index accounts for the fraction of dissimilarity owing to richness difference and is formulated as follows:where a, b and c are the same as in equations (2.1) and (2.2). The first term in equation (2.3) expresses a measure of richness difference, whereas the second term corresponds to the dissimilarity version of βjtu that is independent of richness difference (1 − βjtu, [49]).

The functional β-diversity among regions was decomposed into functional turnover and functional nestedness-resultant components following the same framework [50,51].

Taxonomic and functional β-diversity and their respective components were computed using the R functions from the ‘betapart’ R package (R v. 2.15.2, R development Core Team, 2012).

Functional vulnerability to fishing

Despite the variety of conceptual approaches, there is growing agreement in defining vulnerability as the combination of three components: (i) sensitivity, or the susceptibility of a system to threats, (ii) exposure, or the level of those threats on a system, and (iii) adaptive capacity, or the capacity of the system to prepare for and respond to those threats [41,52,53]. By analogy, the level of ‘functional vulnerability’ in a given ecosystem relies on (i) functional sensitivity, i.e. the extent to which particular traits are more prone to decline in the face of certain threats [54], (ii) exposure or the level of threats, and (iii) functional redundancy, i.e. the degree to which the same functional traits are supported by many and/or abundant species.

In wilderness areas, exposure to fishing is absent or extremely low [29]. Therefore, the vulnerability to fishing in wilderness areas, termed here as ‘baseline vulnerability’, is solely defined as the combination of species sensitivity to fishing, driven by their biological traits (e.g. size, growth and reproductive capacity), and the level of functional redundancy which is determined by the natural distribution of species density and abundances among FEs (electronic supplementary material, figure S2).

The sensitivity of each fish species to fishing was estimated using a fuzzy logic expert system to take into account eight life-history characteristics that are linked to species productivity and other factors that make fish species more or less sensitive to fishing [55]. This indicator has accurate predictive capacity [55] and has been widely recognized as a comprehensive and suitable indicator of fish sensitivity to fishing [56]. The sensitivity to fishing was aggregated at the level of FEs by averaging the sensitivity of all species belonging to the given FE. The scale ranged from 0 to 100.

Functional redundancy is defined as the level of functional equivalence among species in an ecosystem, such that one function may be performed by one or many species, and one species may substitute for another in the latter case [57]. Here, redundancy was assessed using three complementary indices. First, we used the number of species composing each FE in each region (FRS). However, the number of species per FE is only one aspect of functional redundancy. The distribution of species abundances within FEs represents a complementary aspect of redundancy [58-60]. Therefore, secondly we took into account the number of individuals in each FE (FR).

The mean abundance per 500 m2 of a species i in each region was calculated according to the formulawhere N is the number of transects per region and n is the number of individuals of species i. Then, functional redundancy of a given FE for a given region FR was obtained as follows:with representing the mean number of individuals of species i per 500 m2 and S is the number of species in the given FE for that given region.

Thirdly, we computed a redundancy index accounting for both the number of species and their abundances in each FE. We used the Shannon entropy index FRShannon with the rationale that a given FE will have more redundancy if represented by many abundant species. Conversely, FEs will have low redundancy if represented by few and rare species. Accordingly, in each region, FRShannon for each FE was computed aswithwhere, p is the individuals' proportion of species i in the FE, and S is the number of species per FE.

We applied the correction derived from the equivalent number of species [61]:The equivalent number of species is an unbiased measure of Shannon entropy, following Hill's ‘doubling’ property which ensures that the diversity index doubles with the level of diversity, as opposed to nonlinear diversity indices that behave counterintuitively [62].

Quantitatively, the vulnerability of fish FEs to fishing for each region was assessed using a framework based on multi-criteria decision-making (MCDM) and the TOPSIS method (Technique for Order Preference by Similarity to an Ideal Solution). Applied to our specific case, this technique ranks FEs according to their relative distance to the positive and negative ideal solutions, which represent the conditions obtained when the criteria have extreme values [41,63]. The positive ideal solution (A+) corresponds to the conditions where sensitivity to threats is minimum while redundancy is maximum (electronic supplementary material, figure S3). Conversely, the negative ideal solution (A−) corresponds to the conditions where sensitivity to threats is maximum while redundancy is minimum.

Functional vulnerability was then expressed as the relative distance to these positive and negative ideal solutions according to equation (2.9):where V is the vulnerability of functional entity i, d is the distance to A+ and d the distance to A− in Euclidean space. The vulnerability index ranges from 0 if the criteria scores correspond to A+ and to 1 when the criteria scores correspond to A− [41,63].

The vulnerability of an FE is high when both the fishing sensitivity of that FE is high and when redundancy is low [41,64].

Mapping functional sensitivity, redundancy and vulnerability in functional space

The density distribution of functional redundancy, sensitivity and vulnerability within the functional space was estimated using the kernel method with a Gaussian estimation [65]. The smoothing parameter h was estimated using the ad hoc method, which is the optimum h value obtained for the normal distribution [65]:where, n is the number of FEs, and σ2 being the estimated variance for x and y coordinates:The kernel density estimation was computed using the ‘adehabitatHR’ R package (R v. 2.15.2, R development Core Team, 2014).

Functional redundancy, sensitivity and vulnerability were mapped onto the functional space for each region and their match was estimated with the Pearson product-moment correlation.

Results

Similar level of functional diversity across wilderness areas

The greatest number of species and FEs was found in New Caledonia, encompassing 83% of FEs recorded in the three regions (figure 2). By contrast, French Polynesia showed the lowest number of species (42%) and FE (69%) (figure 2; additional information in the electronic supplementary material). Overall, 45% of the total FEs (71) were common to the three regions, whereas only 13% (56) of the species were shared (figure 2; additional information in the electronic supplementary material).

Figure 2.

Sharing of species and functional entities between regions. Venn diagram (top) for the number of species (a) and the number of functional entities (b) for Chagos (orange), New Caledonia (green), and French Polynesia (blue). Percentages indicate the proportion of species and functional entities in each region compared with the total pool. Histograms (bottom) show the proportion (left-axis) and the number of species (a) or functional entities (b) (right-axis) in each region (colour bars), the proportion of unique species (a) or functional entities (b) in each region (black segments on coloured bars), the number of species (a) or functional entities (b) for each pair of regions, and the number of species (a) or functional entities (b) common to the three regions.

The FRic in each region (FRic) ranged from 82% of the global functional volume (French Polynesia) to 94.4% (New Caledonia) (electronic supplementary material, figure S4). The low variability of functional diversity between biogeographic regions is consistent with the weak level of β-functional diversity (functions turnover) between regions ranging from 0.14 to 0.18 (maximum is 1). Conversely, the β—taxonomic diversity was higher, ranging between 0.66 and 0.79 (maximum at 1) (electronic supplementary material, figure S5, table S1 and additional information), implying that most FEs are present independent of species identities.

Heterogeneous distribution of functional redundancy across wilderness areas and functional entities

Despite the overall stability of functional diversity across the biogeographic gradient, the highest functional redundancy of both species per FE (15) and Shannon entropy (17) was found in New Caledonia, whereas Chagos and French Polynesia only reached 8–9 species (figure 3b; additional information in the electronic supplementary material). In each region, more than 60% of FEs were represented by only one species, and 12–20% of the FEs were represented by only two species (figure 3a; additional information in the electronic supplementary material), showing that the majority of FEs had low species redundancy.

Figure 3.

The levels of functional redundancy in terms of species richness and Shannon entropy across functional entities for the three regions. Distribution of functional entities (in percentage) along a gradient of functional redundancy in terms of the number of species by functional entity in each region (a). Relationships between the number of species per functional entity and the Shannon entropy (expressed as equivalent number of species, Material and methods) for each functional entity are shown for each region (b). The bottom panel demonstrates the variation of Shannon entropy for a fixed number of species per functional entity. no. S/FE is the number of species per functional entity.

Mapping variation of Shannon entropy across the functional space showed a highly heterogeneous distribution of functional redundancy for each region, with the highest values disproportionally packed into some parts of the functional space (figure 4a; electronic supplementary material, figure S7), leaving most of the functional space with low or no redundancy. The extreme concentration of functional redundancy in the top left of the figure (figure 4a) was represented by sedentary small to medium size detritivorous, invertivorous, planktivores and omnivorous feeders (e.g. surgeonfishes, damselfishes, butterflyfishes) (figure 4a; electronic supplementary material, figure S6). The top right concentration of redundancy (figure 4a) was characterized by small to medium size invertivorous feeders (mobile prey) (see also additional information in the electronic supplementary material).

Figure 4.

Mapping redundancy, sensitivity and vulnerability in the functional space for fish faunas of the three regions. The top row (a) shows for each region the distribution of functional redundancy within functional entities, measured using the Shannon entropy, across the functional space. The middle row (b) shows the distribution of sensitivity to fishing across the functional space while the bottom row (c) shows the distribution of vulnerability to fishing in this space, and (d) position of vulnerable functional entities in the functional space. Colours indicate the level of vulnerability of functional entities and fish shapes were chosen to illustrate the main genus of each functional entity (see also the electronic supplementary material, figure S6).

High vulnerability of functional entities in wilderness areas

As most of the functional space had low functional redundancy, the variation in vulnerability to fishing was partly driven by sensitivity to fishing (mean Pearson's coefficient across regions = 0.93, p < 0.001). Variations in functional sensitivity and redundancy across the functional space were independently distributed (mean Pearson's coefficient across locations = −0.09, p > 0.1) albeit the most sensitive species systematically showed very low redundancy (electronic supplementary material, figure S8).

Medium to large top predators, both solitary and sedentary, such as groupers (top right) and moray eels (top right), as well as mobile and medium sized schooling species such as the bluefin trevally (Caranx melampygus) and the great barracuda (Sphyraena barracuda) (bottom right) were highly vulnerable to fishing in the three regions (figure 4c,d; electronic supplementary material, figure S6) owing to their high sensitivity to fishing (figure 4b) and low redundancy (1 species per FE (figure 4a)). For the same reasons, large invertivores species feeding on mobile prey such as the humphead wrasse (Cheilinus undulatus), large invertivores (sedentary and mobile), detritivores and planktivores, located in the upper right side of the functional space (figure 4b; electronic supplementary material, figure S6) showed a high vulnerability to fishing (figure 4b; addition information in the electronic supplementary material).

At the centre of the functional space, we observed a high fishing vulnerability peak owing to the presence of two large, mobile and medium sized schooling species; the bumphead parrotfish (Bolbometopon muricatum), an invertebrate sessile feeder and the humpback unicornfish (Naso brachycentron), a macroalgae feeder (figure 4b; electronic supplementary material, figure S6). At the right side of the functional space, we observed a second peak of vulnerability to fishing (hence high sensitivity coupled with low redundancy) characterized by the bicolour parrotfish (Cetoscarus ocellatus), a large bioeroder [66] and the bluespine unicornfish (Naso unicornis), a macroalgae feeder, both being the only species in their respective FEs and highly sensitive to fishing, explaining their high vulnerability.

Conversely, the high diversity of small to medium size herbivorous–detritivorous, omnivorous, invertivorous and planktivores species in some FEs (top left of the functional space), contributed to their high functional redundancy, counterbalancing (mean Pearson's coefficient across locations of −0.21, p < 0.001) their high sensitivity to fishing and making them less vulnerable (figure 4). The higher functional redundancy in New Caledonia decreased the vulnerability of certain FEs. For example, the functional entity composed of medium size invertebrate feeders (mainly from the triggerfishes family) was composed of nine species in New Caledonia and five in French Polynesia, thereby, in concert with Shannon entropy, reducing its vulnerability of 12% compared with French Polynesia. Variation in redundancy and vulnerability across the functional space depicted by the third and fourth axes of PCoA is showed in electronic supplementary material, figure S9.

Discussion

Homogeneous level of functional diversity and functional vulnerability across the Indo-Pacific biogeographic gradient

Surprisingly, despite the large differences in species richness and identity among the three regions, we found very consistent patterns of functional diversity along the species richness biogeographic gradient [67], suggesting that the set of ecological functions for the growth and persistence of coral reefs are supported independently of species identity and richness, thus serving as a ‘natural’ baseline indicator. As such, the presence/absence of FEs compared with wilderness areas may be a useful indicator to assess ecosystem conditions under human pressure. This result is in line with a previous study on the scale of biogeographic regions demonstrating that functions are maintained along a gradient of species richness [21]. Similar consistency was found for functional vulnerability, suggesting that combinations of fish traits have quasi-similar levels of redundancy across coral reefs and the biogeographic gradient, albeit New Caledonia had lower levels of vulnerability for several FEs owing to higher functional redundancy.

High vulnerability of ecosystem functioning

Studies of vulnerability have typically focused on species rather than ecological functions [41], and often overlook abundance distributions among species [21,41]. In addition, the assumption that functional diversity in wilderness areas, seen as a reference, is not buffered against extinction owing to relatively high functional vulnerability has never been tested. Here, we filled this gap by assessing the baseline functional vulnerability of fish assemblages to fishing in some of the most isolated coral reef ecosystems across the Indo-Pacific. We found that even in the near-absence of fishing and considering species abundance, most FEs remain highly vulnerable in each wilderness ecosystem because species diversity is overwhelmingly packed into a small set of FEs, leaving most FEs with no redundancy or functional insurance. Importantly, functional sensitivity and functional redundancy are independently distributed among the combinations of traits. As such, redundancy is seldom compensating for the high sensitivity of some traits, the most sensitive ones often having no redundancy.

Our study is complementary to previous macro-ecological investigations showing the heterogeneous distribution of redundancy across FEs [21,41], and attests that the same pattern translates to the local scale when considering species abundances. This result is not trivial because even FEs with low species redundancy may be preserved if those species are widespread, abundant and weakly sensitive to threats. This is not the case for fishes on coral reefs. Consequently, high biodiversity and abundance of fishes is unlikely to buffer coral reef ecosystems against functional diversity loss. This result is consistent with the observations that marine ecosystem functions and services scale positively and do not saturate with the level of species and FE richness within local faunas [68,69].

Vulnerability can provide early warning signals of fishing pressure, because the processes performed by vulnerable traits are most likely to be the first to decline, even under low pressure [70]. Indeed, from a sustainable management perspective, this study highlights key vulnerable FEs that should be under particular conservation measures. For example, the bluespine unicornfish (N. unicornis) has been identified as one of the most important macroalgal consumers on coral reefs [71] yet is targeted by commercial, recreational and artisanal fisheries [72]. Macroalgal feeders play a critical role in preventing and reversing coral-algal shifts and the finding of their high baseline vulnerability to fishing is worrying [8,9]. Despite recent studies demonstrating signs of overexploitation in the Pacific [72], this species is listed as ‘Least concern’ in the Red list of Threatened Species by the International Union for Conservation of Nature (IUCN). This is also the case for top predators such as the bluefin trevally (C. melampygus) and the great barracuda (S. barracuda) which are highly sensitive to fishing [73,74], while the bluefin trevally has also been found in very low abundance in marine PAs compared with wilderness [28]. The IUCN Red list criteria could benefit from taking a trait based approach to uncover the species vulnerability from a functional perspective.

Conclusion

With biodiversity being lost in the Anthropocene, the need to carefully manage functional diversity has never been greater. The slow decrease in abundance of functionally important species is insidious, despite strong implications for coral ecosystem functioning and resilience [9]. High redundancy can highlight the core processes that may persist even under high levels of fishing pressure because they are highly buffered against local species extinction. While fishing pressure should decrease on the most vulnerable FEs, it could potentially be directed toward the least vulnerable FEs highlighted in the functional space. Local adaptive management measures toward those species could focus on gear restrictions, minimum size limits [75], quotas or fishing bans during spawning times [76,77] or modified angling techniques (e.g. [78-80]). However, those species often do not garner the same fishing interest owing to their low commercial values [81]. Therefore, the burden of reducing fishing pressure on the most vulnerable FEs should be shared along market value chains [82], as fisherman targeting certain species operate under incentives determined by trade and the agents involved at both local and international scales [83,84]. Ultimately, a multilayer management approach considering the complexity of the socio-ecological system is necessary to maintain ecosystem functioning [82].