A large body of knowledge exists on individual anthropogenic threats that have an impact on marine biodiversity in the Mediterranean Sea, although we know little about how these threats accumulate and interact to affect marine species and ecosystems. In this context, we aimed to identify the main areas where the interaction between marine biodiversity and threats is more pronounced and to assess their spatial overlap with current marine protected areas in the Mediterranean. Mediterranean Sea.\ud We first identified areas of high biodiversity of marine mammals, marine turtles, seabirds, fishes and commercial or well-documented invertebrates. We mapped potential areas of high threat where multiple threats are occurring simultaneously. Finally we quantified the areas of conservation concern for biodiversity by looking at the spatial overlap between high biodiversity and high cumulative threats, and we assessed the overlap with protected areas.\ud Our results show that areas with high marine biodiversity in the Mediterranean Sea are mainly located along the central and north shores, with lower values in the south-eastern regions. Areas of potential high cumulative threats are widespread in both the western and eastern basins, with fewer areas located in the south-eastern region. The interaction between areas of high biodiversity and threats for invertebrates, fishes and large animals in general (including large fishes, marine mammals, marine turtles and seabirds) is concentrated in the coastal areas of Spain, Gulf of Lions, north-eastern Ligurian Sea, Adriatic Sea, Aegean Sea, south-eastern Turkey and regions surrounding the Nile Delta and north-west African coasts. Areas of concern are larger for marine mammal and seabird species.\ud These areas may represent good candidates for further research, management and protection activities, since there is only a maximum 2% overlap between existing marine protected areas (which cover 5% of the Mediterranean Sea) and our predicted areas of conservation concern for biodiversity

The evolutionary dissimilarity between communities (phylogenetic beta diversity PBD) has been increasingly explored by ecologists and biogeographers to assess the relative roles of ecological and evolutionary processes in structuring natural communities. Among PBD measures, the PhyloSor and UniFrac indices have been widely used to assess the level of turnover of lineages over geographical and environmental gradients. However, these indices can be considered as ‘broad-sense’ measures of phylogenetic turnover as they incorporate different aspects of differences in evolutionary history between communities that may be attributable to phylogenetic diversity gradients. In the present study, we extend an additive partitioning framework proposed for compositional beta diversity to PBD. Specifically, we decomposed the PhyloSor and UniFrac indices into two separate components accounting for ‘true’ phylogenetic turnover and phylogenetic diversity gradients, respectively. We illustrated the relevance of this framework using simple theoretical and archetypal examples, as well as an empirical study based on coral reef fish communities. Overall, our results suggest that using PhyloSor and UniFrac may greatly over-estimate the level of spatial turnover of lineages if the two compared communities show contrasting levels of phylogenetic diversity. We therefore recommend that future studies use the ‘true’ phylogenetic turnover component of these indices when the studied communities encompass a large phylogenetic diversity gradient.

Summary1. Current global changes make it important to be able to predict which interactions will occur in the emerging ecosystems. Most of the current methods to infer the existence of interactions between two species require a good knowledge of their behaviour or a direct observation of interactions. In this paper, we overcome these limitations by developing a method, inspired from the niche model of food web structure, using the statistical relationship between predator and prey body size to infer the matrix of potential interactions among a pool of species. 2. The novelty of our approach is to infer, for any species of a given species pool, the three species-specific parameters of the niche model. The method applies to both local and metaweb scales. It allows one to evaluate the feeding interactions of a new species entering the community. 3. We find that this method gives robust predictions of the structure of food webs and that its efficiency is increased when the strength of the body-size relationship between predators and preys increases. 4. We finally illustrate the potential of the method to infer the metaweb structure of pelagic fishes of the Mediterranean sea under different global change scenarios.

Knowledge of species composition and their interactions, in the form of interaction networks, is required to understand processes shaping their distribution over time and space. As such, comparing ecological networks along environmental gradients represents a promising new research avenue to understand the organization of life. Variation in the position and intensity of links within networks along environmental gradients may be driven by turnover in species composition, by variation in species abundances and by abiotic influences on species interactions. While investigating changes in species composition has a long tradition, so far only a limited number of studies have examined changes in species interactions between networks, often with differing approaches. Here, we review studies investigating variation in network structures along environmental gradients, highlighting how methodological decisions about standardization can influence their conclusions. Due to their complexity, variation among ecological networks is frequently studied using properties that summarize the distribution or topology of interactions such as number of links, connectance, or modularity. These properties can either be compared directly or using a procedure of standardization. While measures of network structure can be directly related to changes along environmental gradients, standardization is frequently used to facilitate interpretation of variation in network properties by controlling for some co-variables, or via null models. Null models allow comparing the deviation of empirical networks from random expectations and are expected to provide a more mechanistic understanding of the factors shaping ecological networks when they are coupled with functional traits. As an illustration, we compare approaches to quantify the role of trait matching in driving the structure of plant-hummingbird mutualistic networks, i.e. a direct comparison, standardized by null models and hypothesis-based metaweb. Overall, our analysis warns against a comparison of studies that rely on distinct forms of standardization, as they are likely to highlight different signals. Fostering a better understanding of the analytical tools available and the signal they detect will help produce deeper insights into how and why ecological networks vary along environmental gradients.

Marine megafauna, the largest animals in the oceans, serve key roles in ecosystem functioning. Yet, one-third of these animals are at risk of extinction. To better understand the potential consequences of megafaunal loss, here we quantify their current functional diversity, predict future changes under different extinction scenarios, and introduce a new metric [functionally unique, specialized and endangered (FUSE)] that identifies threatened species of particular importance for functional diversity. Simulated extinction scenarios forecast marked declines in functional richness if current trajectories are maintained during the next century (11% globally; up to 24% regionally), with more marked reductions (48% globally; up to 70% at the poles) beyond random expectations if all threatened species eventually go extinct. Among the megafaunal groups, sharks will incur a disproportionate loss of functional richness. We identify top FUSE species and suggest a renewed focus on these species to preserve the ecosystem functions provided by marine megafauna.

The Mediterranean Sea is a hotspot of biodiversity, and climate warming is expected to have a significant influence on its endemic fish species. However, no previous studies have predicted whether fish species will experience geographic range extensions or contractions as a consequence of warming. Here, we projected the potential future climatic niches of 75 Mediterranean Sea endemic fish species based on a global warming scenario implemented with the Mediterranean model OPAMED8 and a multimodel inference, which included uncertainty. By 2070-2099, the average surface temperature of the Mediterranean Sea was projected to warm by 3.1 1C. Projections for 2041-2060 are that 25 species would qualify for the International Union for the Conservation of Nature and Natural Resources (IUCN) Red List, and six species would become extinct. By 2070-2099, 45 species were expected to qualify for the IUCN Red List whereas 14 were expected to become extinct. By the middle of the 21st century, the coldest areas of the Mediterranean Sea (Adriatic Sea and Gulf of Lion) would act as a refuge for cold-water species, but by the end of the century, those areas were projected to become a 'cul-de-sac' that would drive those species towards extinction. In addition, the range size of endemic species was projected to undergo extensive fragmentation, which is a potentially aggravating factor. Since a majority of endemic fishes are specialists, regarding substratum and diet, we may expect a reduced ability to track projected climatic niches. As a whole, 25% of the Mediterranean Sea continental shelf was predicted to experience a total modification of endemic species assemblages by the end of the 21st century. This expected turnover rate could be mitigated by marine protected areas or accelerated by fishing pressure or competition from exotic fishes. It remains a challenge to predict how these assemblage modifications might affect ecosystem function.

The Mediterranean Sea (0.82% of the global oceanic surface) holds 4%-18% of all known marine species (~17,000), with a high proportion of endemism [1, 2]. This exceptional biodiversity is under severe threats [1] but benefits from a system of 100 marine protected areas (MPAs). Surprisingly, the spatial congruence of fish biodiversity hot spots with this MPA system and the areas of high fishing pressure has not been assessed. Moreover, evolutionary and functional breadth of species assemblages [3] has been largely overlooked in marine systems. Here we adopted a multifaceted approach to biodiversity by considering the species richness of total, endemic, and threatened coastal fish assemblages as well as their functional and phylogenetic diversity. We show that these fish biodiversity components are spatially mismatched. The MPA system covers a small surface of the Mediterranean (0.4%) and is spatially congruent with the hot spots of all taxonomic components of fish diversity. However, it misses hot spots of functional and phylogenetic diversity. In addition, hot spots of endemic species richness and phylogenetic diversity are spatially congruent with hot spots of fishery impact. Our results highlight that future conservation strategies and assessment efficiency of current reserve systems will need to be revisited after deconstructing the different components of biodiversity.

Marine protected areas (MPAs) are major tools to protect biodiversity and sustain fisheries. For species with a sedentary adult phase and a dispersive larval phase, the effectiveness of MPA networks for population persistence depends on connectivity through larval dispersal. However, connectivity patterns between MPAs remain largely unknown at large spatial scales. Here, we used a biophysical model to evaluate connectivity between MPAs in the Mediterranean Sea, a region of extremely rich biodiversity that is currently protected by a system of approximately a hundred MPAs. The model was parameterized according to the dispersal capacity of the dusky grouper Epinephelus marginatus, an archetypal conservation-dependent species, with high economic importance and emblematic in the Mediterranean. Using various connectivity metrics and graph theory, we showed that Mediterranean MPAs are far from constituting a true, well-connected network. On average, each MPA was directly connected to four others and MPAs were clustered into several groups. Two MPAs (one in the Balearic Islands and one in Sardinia) emerged as crucial nodes for ensuring multi-generational connectivity. The high heterogeneity of MPA distribution, with low density in the South-Eastern Mediterranean, coupled with a mean dispersal distance of 120 km, leaves about 20% of the continental shelf without any larval supply. This low connectivity, here demonstrated for a major Mediterranean species, poses new challenges for the creation of a future Mediterranean network of well-connected MPAs providing recruitment to the whole continental shelf. This issue is even more critical given that the expected reduction of pelagic larval duration following sea temperature rise will likely decrease connectivity even more.