Introduction

A comparative analysis of life-history traits of the Superorder Batoidea offers the opportunity for fisheries managers to take appropriate action to prevent their overfishing. The members of the Superorder Batoidea are widely distributed throughout the Northeast Atlantic Ocean (NEA) and the Mediterranean Sea (MED) (Ebert and Stehmann 2013; Last et al. 2016) and are being exploited by commercial fisheries as target and by-catch species. Generally recognized as K-selected species (MacArthur and Wilson 1967; Camhi et al. 1998), they are characterized by slow growth, late sexual maturity, low fecundity and long life (Hoenig and Gruber 1990), limiting their population’s capacity to recover from overfishing (Camhi et al. 1998; ICES 2020a, b, c). To reduce the fisheries pressure on their populations, batoids have been temporarily exempted from the European Landing Obligation policy (Uhlmann et al. 2019; ICES 2020a, b, c), as they appear to have a high chance of survival after being discarded overboard at sea (e.g. Enever et al. 2009; Ellis et al. 2012; Catchpole et al. 2017). Conditional to this temporal exemption, the European Commission initiated a research trajectory called “Roadmap for rays and skates" (Commission of the European Union 2018) to increase current scientific knowledge, and fill data gaps concerning these species and their survivability. It is therefore of scientific importance to determine the state of the art of batoids research and data availability in the NEA, MED and to analyse the available data to provide useful outputs.

Currently, due to the limited data available on commercial fisheries, most batoids in the NEA and the MED are managed under a common Total Allowable Catch (TAC) and all stocks are assessed as data-limited (ICES 2012, 2020; Lart 2019), due to a lack of robust and long-term length-age data. Although, sampling and compiling such data for each species may be very complicated and/or costly, finding methods to distinguish and/or group species based on their vulnerability to fisheries and environmental changes and capacity to recover from overfishing would be of great value to increase the capacity to understand the true fishing pressure on each “group” of species, hence the quality of stock assessments and our capacity to ensure the sustainability of the fishing stocks.

Life-history traits vary widely between taxonomic and/or functional groups, based on species-specific adaptations and strategies to respond to prevailing environmental conditions (e.g. Walsh 1994). Functionally and taxonomically closely related species and even sub-populations of the same species can have significantly different life-history traits due to regional adaptations to environmental gradients (Fromentin and Fonteneau 2001; Félix et al. 2011), and/or due to different degrees of fishing pressure (Devine et al. 2012; Claireaux et al. 2018). Such variations can be used to understand the variability between and within species and to provide preliminary evidence for the adaptability to changes in the environment, such as climate change (Dalgleish et al. 2010), habitat degradation (Öckinger et al. 2010), and overfishing (Hoenig and Gruber 1990; Cortés et al. 2010; Patrick et al. 2010; Hobday et al. 2011). These features make life-history traits useful proxies for estimating the tolerance of a species to natural and/or anthropogenic disturbances, and they could be used to identify common life strategies among groups of species (Félix et al. 2011; Pecuchet et al. 2017).

As far as batoids in the NEA and the MED are concerned, studies of life-history traits usually report values for a limited number of traits for a single or a few species based on data from samples collected on a restricted regional scale (e.g. Du Buit 1977; Nottage and Perkins 1983; Ryland and Ajayi 1984; Brander and Palmer 1985). Interspecific (e.g., Gallagher et al. 2005) and/or interregional (e.g., McCully et al. 2012) analyses of batoids are scarce and outdated, but what is known about them was included within a recent synopsis (Last et al. 2016). Reviews and meta-analyses of fish taxa such as mackerel, tuna, and bonitos (Scombridae) (Juan-Jordá et al. 2013), flatfish (Pleuronectidae) (Félix et al. 2011), and multi-species analyses (Rochet 2000; King and McFarlane 2003; Frédou et al. 2016) resulted in the delineation of clusters of species with similar life-history strategies. Although batoids are generally identified as a K-selected species (Camhi et al. 1998), this probably oversimplifies the variety of adaptations exhibited in the wild (Nichols et al. 1976; Bielby et al. 2007). Life-history traits can be used to identify similar species or groups of species with higher vulnerability (e.g., late maturity and slow growth) at an individual (e.g., discard survival) and a population level (e.g., overfishing, climate change) and establish spatial–temporal management measures accordingly (i.e., TAC, maximum allowed by-catch, gear ban, gear selectivity, spatiotemporal avoidance). For instance, based on the life-history traits values of a species and/or groups of species, these could be identified as more likely/unlikely to die after being discarded back at sea. These could ultimately provide proxies and/or be used to tune discard survival estimation, to be used and accounted for stock assessments for elasmobranchs in the future.

Considering the above, this study aimed to investigate inter-and intraspecific patterns in life-history traits of skate and other batoid species in the NEA and the MED through a literature review. Via principal component analysis (PCA) of the compiled life-history traits, species, and/or groups of species to which specific management strategies could be recommended due to their similar life-history strategies were identified. Furthermore, the compiled traits were then summarized by region of origin of the studies to analyse the plasticity of the traits across the distribution range of a species. Finally, both the compilation of available information on life-history traits and the detection of data gaps is essential for prioritizing future research efforts on batoids in the NEA and the MED.

Materials and methods

Literature review

A literature review (Khan et al., 2003; Pascual et al. 2017; Fennell et al. 2021) was conducted to identify life-history trait published values for batoids in the NEA and MED. First, the literature was screened to identify commonly reported life-history traits parameters for 14 batoids species: 11 skates (Rajidae) and 3 rays (Torpedinidae and Dasyatidae) in the NEA and the MED (Table 1). Batoid species of interest included: Starry ray (Amblyraja radiata); Common skate complex (Dipturus batis complex); Sandy ray (Leucoraja circularis), Shagreen ray (Leucoraja fullonica); Cuckoo ray (Leucoraja naevus); Blonde ray (Raja brachyura); Thornback ray (Raja clavata); Small-eyed ray (Raja microocellata); Spotted ray (Raja montagui); Undulate ray (Raja undulata); White skate (Rostroraja alba); Electric ray (Tetronarce nobiliana); Marbled electric ray (Torpedo marmorata); and Common stingray (Dasyatis pastinaca). From this first screening, a total of ten quantitative life-history parameters describing growth, reproduction and feeding habits of batoids in the NEA and MED were identified (Table 2).

Table 1 Common, scientific, and three-letter identifier codes (FAO 2021) of batoids (Batoidea) of the Northeast Atlantic Ocean and the Mediterranean Sea were analysed in the literature review
Table 2 List of life-history traits and associated search terms that were used to screen published primary literature for life-history trait values of batoid skate species. Each combination is preceded by the combination of the common and the scientific name of each species (see Table 1)

Next, the Web of Science (WoS), Scopus and Google scholar databases were screened for relevant published values, using combinations of species- and trait-specific search terms for each of the identified life-history trait parameter identified previously. The scientific and common names (e.g., Raja clavata AND thornback ray), followed by combinations of keywords and derivatives related to each of the identified life-history traits (Table 2) were one by one used to identify published life-history trait values. Grey literature sources, such as species catalogues (e.g., Ebert and Stehmann 2013; Froese and Pauly 2019), and reports from the ICES Working Group for Elasmobranch Fishes (WGEF), were also consulted. Finally, the references list of the identified publications (i.e. papers, books, reports, etc.) were screened identify relevant peer-reviewed publications not included until this moment (Khan et al. 2003; Pascual et al. 2017; Fennell et al. 2021). Publications were excluded from the study based on the type of publication, geographic location of the study and whether the study has been cited in others. Only original research and other review paper were considered, while opinion and editorial publications were discarded. Values reported for a species were only included in our data base if the study was performed in the NEA and MED. All obtained values were compiled in a unique data base, together with information on the geographical location in which the study was performed, species and sex of the individuals analysed. The location was defined based on FAO major fishing areas of the Northeast Atlantic and Mediterranean Sea (FAO 2021). Values that could not be linked to a specific region/sub-region were defined as “Global”, and were compiled in the Supplementary Table 1, but were not included in the data analysis. Wherever possible, separate values were recorded for sex. When this was impossible, the same value was assigned to both males and females.

Growth traits were described by the length at infinity (L), and growth rate (K) (Table 2) parameters from the von Bertalanffy growth model (von Bertalanffy 1938). Both parameters have been widely published in the scientific fisheries literature and are being used in age-structured models. For elasmobranchs these parameters are not used the stock assessment due to data limitations and biological characteristics, but are key for these more robust assessment methods. Additionally, the lifespan of the species corresponds to the estimated maximum expected age. Natural mortality, was defined as the natural (i.e., predation, disease, and old age) loss rate of the population.

In fisheries science, the age (Agem) and length at maturity (Lm) (Table 2) are commonly estimated from logistic regression between the size and the maturity states (mature/immature), providing estimates of ages and sizes at which 50% of the individuals of a species are sexually mature (Schnute and Richards 1990; Chen et al. 1994). Although maturation is a continuous process and achieved through a series of transitional states (e.g., range of values), age and length at maturity are point estimates for describing the species’ reproductive traits. These parameters have been published for most skate and batoid species. However, because their estimation relies on dissections of a sufficiently large number (e.g., > 103) of individuals, they may be of limited value for species caught in small amounts and/or subject to misidentification, such as small-eyed (Raja microocellata), shagreen (Leucoraja fullonica) and sandy (Leucoraja circularis) ray (McCully and Walls 2015; ICES 2020a, b, c). Furthermore, in contrast to bony fishes, ageing elasmobranchs commonly requires the analysis of vertebral sections and images techniques (Campana 2014; SUMARiS 2019). It is a still evolving discipline, meaning that a large proportion of the methods used in the past and their associated products (peer reviewed papers), may now be completely obsolete (Campana 2014). Fecundity data were also compiled and defined as the average total number of offspring (eggs or pups) produced per year per female. Development was defined as the time needed for the young to fully develop inside an egg or its mother’s uterus. Regional fecundity and development data are usually limited to particular studies and/or reports and are therefore commonly generalized for a species across its distribution. Differences in reproduction types were analyzed by categorizing the species as either oviparous (producing eggs) or viviparous (producing living young). The trophic level derived from FishBase (Froese and Pauly 2019) was used to describe the feeding strategy of each species.

Exploratory analysis

Interspecific analysis

We analysed the possibility of grouping and dividing the 14 batoid species based on their life-history traits values. A summarized data set was obtained by taking the mean value for each life-history trait parameter for each of the 14 batoid species, and for each sex, from all studies from the NEA and the MED. Correlation between continuous life-history parameters was first explored using the Spearman correlation index, using the corr {stats}, and plotted using the corrplot {corrplot} (Wie and Simko 2017) functions. The summarized data set was used to analyze the interspecific pattern of life-history traits through a principal component analysis (PCA), using the function “prcomp” {stats To perform the PCA, five species: common stingray (Dasyatis pastinaca), electric (Tetronarce nobiliana), sandy, shagreen and small-eyed ray; were removed from the summarized data set because of the absence of published values for age at maturity (Agem). Similarly, the natural mortality and lifespan were not included in the analysis as values for these parameters were available for only three and five of the analyzed species respectively. The remaining seven continuous life-history traits (Agem, Lm, L, K, Development, Fecundity, and Trophic level) of nine batoid species were then used to perform the PCA. Additionally, clustering patterns between batoid species were explored using the Euclidean distance, and three different linkage clustering methods (complete, single and average), using the “dist” {stats}, “hclust” {stats}, “cutree” {stats}, and “ggdendrogram” {ggdendro} (De Vries and Ripley 2020) functions. The three clustering methods differ in the distance between clusters is calculated. Single linkage merges the two clusters whose two closest members have the smallest distance, while complete linkage merges the two clusters whose merger has the smallest distance. Finally, the average linkage represents an intermediate method which takes the average distance between the members of two clusters.

Intraspecific analysis

We aimed to analyse the intraspecific variation of the life-history traits of the 14 batoid species by summarizing the compiled data per geographical region of origin of each study using PCA. Two geographical regions were defined. The North East Atlantic (NEA) region, including the: Skagerrak, North Sea, English Channel, Celtic Sea, Bay of Biscay and the Portuguese waters sub-regions, and the Mediterranean and Black Sea region (MED), including the Mediterranean Sea and the Black Sea sub-regions (Fig. 1). The lack of regional data for most species and life-history traits, impeded to perform the analysis as planned. For instance, values for at least two different sub-regions were obtained for only six (starry, Amblyraja radiata; cuckoo, Leucoraja naevus; blonde, Raja brachyura; thornback, Raja clavata; spotted, Raja montagui; and undulate, Raja undulata) of the 14 analysed species and for only four (age at maturity, length at maturity, length at infinity, and growth rate) of the 10 life-history traits. Therefore values from different regions and sub-regions for each of the six species were compared by comparing the average and standard deviation of the reported values for each parameter in each sub-region/region.

Fig. 1
figure 1

© FAO 1990–2021

FAO major fishing areas of the Northeast Atlantic. Skagerrak and Kattegat (3.a), North Sea (4.a-c), Celtic Sea (7.g-h), English Channel (7.d-7.e), Bay of Biscay (8.a-c), Portuguese Waters (9.a-b) and Mediterranean Sea.

All analysis and plotting were performed in R statistical software version 4.0.2 (R Core Team 2020).

Results

Systematic literature review

The full data set of compiled life-history traits and all associated references can be accessed in Supplementary Table 1 and Supplementary Table 2, respectively. The origin of the studies varied with the species analysed following their natural distribution (Fig. 2). Life-history traits of batoids from the genera Dipturus, Raja, and Leucoraja were mainly reported for NEA waters (e.g., North Sea, British Isles, Celtic/Irish Sea, amongst others), whereas life-history traits for more southernly distributed species, such as the white skate, common stingray, electric and marbled electric (Torpedo marmorata) ray (e.g. Dulvy et al. 2006; Notarbartolo Di Sciara et al. 2009; Serena et al. 2009), were mainly studied in the MED (e.g., Balearic Island, Tunisian coast).

Fig. 2
figure 2

Adapted from Ebert and Stehmann 2013

Distribution of the starry ray (Amblyraja radiata), common/blue skate (Dipturus batis complex), sandy (Leucoraja circularis), shagreen (Leucoraja fullonica), cuckoo (Leucoraja naevus), blonde (Raja brachyura), thornback (Raja clavata), small-eyed (Raja microocellata), spotted (Raja montagui), undulate ray (Raja undulata), white skate (Rostroraja alba), electric (Totronarce nobiliana), marbled electric ray (Torpedo marmorata) and common stingray (Dasyatis pastinaca).

For certain life-history traits, the overview (Table 3) clearly shows where essential information is still lacking (e.g., absent values for natural mortality, lifespan and, to a lesser degree, age at maturity). Life-history traits for the thornback, blonde, undulate and spotted ray have been commonly reported in the scientific literature, with studies encompassing different geographical sub-regions. Contrarily, values for rare and less abundant species, such as the sandy, shagreen, and small-eyed ray, are equally rare and limited, leaving important data gaps (Table 3).

Table 3 Averaged (Standard deviation) life-history traits values skate of the NEA and the MED Sources: See Supplementary Table 1

Strong (> 0.5) Spearman correlation values were obtained between growth and maturity parameters (Fig. 3), because species with larger lengths at infinity (L), take longer to mature (Agem), reaching larger sizes (Lm) due to a slower growth rate (K). The trophic level was found to be positively correlated to the development time, with higher trophic level species requiring more time to fully develop. Fecundity did not register a strong correlation with any of the other life-history traits analysed.

Fig. 3
figure 3

Spearman correlation matrix (corrplot {corrplot}) for life-history traits of batoids and related species in the NEA and MED. Acronyms are detailed in Table 2. Direct (0,1) and inverse (0,-1) correlation, are represented with blue and red tones respectively

Exploratory analyses

Interspecific analysis

The first two principal components of the PCA explained 78.56% of the total variance contained within the mean of the nine life-history traits analysed of nine species of batoids (starry, thornback, blonde, spotted, undulate cuckoo and electric ray, white skate and common skate complex) (Fig. 4). Growth and reproductive life-history traits (i.e., age and length at maturity, length at infinity, and growth rate), were mainly related with the first PCA axis (52.63% of total variation), and the trophic level and development by the second axis (25.93% of total variation). Fecundity was largely related with the third PC only.

Fig. 4
figure 4

Biplot of the two first principal components of the PCA of age at maturity (AgeM), length at maturity (Lm), length at infinity (Linf), growth rate, (K), Fecundity (Fec), development (Develop) and trophic level (Trophic), of male (_m) and female (_f) blonde (RJH), cuckoo (RJN), spotted (RJM), thornback (RJC) and undulate (RJU) and marbled electric (TTR) rays, common skate complex (RJB), starry (RJR) and white skate (RJA) in the NEA and MED

From this analysis two “groups” were distinguished. On one hand, the white skate and the common skate complex are considerably different from all other analyzed species, because of their very large sizes (L∞:189.90–240.90 cm; Lm: 119.30–197.50 cm) and late maturation (Agem: 11.00–23.4 years) due to a low growth rate (K:0.05–0.06 /year) (Table 3). However, these species could also be divided based on the second PCA axis, because they differ in their trophic level and time needed for the development of the young. White skate has the longest development time (15 months) of all analyzed species and one of the highest trophic levels (4.4), while the common skate complex’s eggs develop relatively quickly (4 months) and adults are lower in the trophic web (4.0). All remaining species appear to be clustered together by their smaller sizes (L:57.3–126.9 cm) and faster growth (K: 0.12–0.26 /year), reaching maturity at younger ages (Agem:4.7–8.8) and smaller sizes (Lm:27.0–83.7 cm). Nevertheless, based on the second PCA axis, two species stood out from this group: the marbled electric ray, for its long development time (10 months) and high trophic level (4.5), and at the opposite side of this axis, the undulate ray, with one of the shortest development times (3 months) and trophic levels (3.5). Similarly, the cluster analysis (Fig. 5) clearly shows the presence of at least two distinct groups of species, with, on one hand, the white skate and the common skate complex and on the other hand, the seven remaining species. In the latter group, the marbled electric ray was identified as the most different species. The described grouping pattern did not change significantly regarding the three different linkage clustering methods used (single, complete, and average linkage).

Fig. 5
figure 5

Cluster (complete linkage) plot of the Euclidean distance of age at maturity (Agem), length at maturity (Lm), length at infinity (Linf), growth rate, (K), Fecundity (Fec), development (Develop) and trophic level (Trophic) of blonde (RJH), cuckoo (RJN), spotted (RJM), thornback (RJC) and undulate (RJU) and marbled electric (TTR) rays, common skate complex (RJB), thorny (RJR) and white skate (RJA)

Intraspecific analysis

The intraspecific analysis of life-history traits between different geographical regions was severely limited by the lack of available data for most species and most life-history traits (Supplementary Table 1 and Fig. 1). A summary of the available life-history traits studies per sub-region and species is presented in Supplementary Table 3. Although, their distribution generally encompasses more than one sub-region, no comparison between sub-regions and/or regions was possible for the common skate, sandy, shagreen, small-eyed ray, white skate, electric and marbled electric ray and common stingray, due to the absence of studies in more than two regions. For example, the electric and common stingray distribute beyond the Mediterranean Sea (i.e. North Sea, Celtic Sea, English Channel and Bay of Biscay), but no studies reporting life-history trait values were found from other region/sub-regions. Similarly, common skate complex life-history traits’ studies were only found for the Celtic Sea, while the species can also be found in the North Sea (Supplementary Table 3).

Sufficient values to compare life-history traits from studies performed in the NEA and MED region were only available for the thornback and blonde ray (Tables 4 and 5). Reported values for the spotted and undulated ray covered several regions of the NEA, but did not included studies performed in the MED (Tables 6 and 7). Similarly, despite occurring from the North Sea to the Atlantic Iberian waters (ICES 2020b; 2020c), life-history trait values for the cuckoo ray (Leucoraja naevus) had only been reported for the North and Celtic Sea (Table 8). Although, the so far mentioned species, could be identified as the most data abundant species, the studies from which life-history traits were obtained did not covered homogeneously their entire distribution. For instance, although all these species occur in the Bay of Biscay and English Channel, values for these sub-regions were reported only reported for two (thornback and blonde ray) and one (blonde ray) species respectively. Values for all life-history traits of the starry rays (Amblyraja radiata) were obtained for the North Sea and Skagerrak, while no reports for the English Channel and Celtic Sea, where this species also occurs were found (Table 9).

Table 4 Cross-study average (standard deviation) life-history trait values for thornback ray (Raja clavata) in different geographic regions (Northeast Atlantic and Mediterranean and Black Sea) and sub-regions (North Sea, Celtic Sea, English Channel, Bay of Biscay, Portuguese waters, Mediterranean, and Black Sea) Sources: See Supplementary Table 1
Table 5 Cross-study average (standard deviation) life-history trait values for blonde ray (Raja brachyura) in different geographic regions (Northeast Atlantic and Mediterranean and Black Sea) and sub-regions (North Sea, Celtic Sea, English Channel, Bay of Biscay, Portuguese waters and Mediterranean Sea) Sources: See Supplementary Table 1
Table 6 Cross-study average (standard deviation) life-history trait values for spotted ray (Raja montagui) in different geographic regions (Northeast Atlantic and Mediterranean and Black Sea) and sub-regions (North Sea, Celtic Sea, English Channel, Bay of Biscay and Portuguese waters) Sources: See Supplementary Table 1
Table 7 Cross-study average (standard deviation) life-history trait values for undulate ray (Raja undulata) in different geographic regions (Northeast Atlantic and Mediterranean and Black Sea) and sub-regions (North Sea, Celtic Sea, English Channel, Bay of Biscay and Portuguese waters) Sources: See Supplementary Table 1
Table 8 Cross-study average (standard deviation) life-history trait values for cuckoo ray (Leucoraja naevus) in the NEA and sub-regions (North Sea, Celtic Sea, English Channel, Bay of Biscay, Portuguese waters) Sources: See Supplementary Table 1
Table 9 Cross-study average (standard deviation) life-history trait values for starry skate (Amblyraja radiata) in different the NEA and sub-regions (North Sea and Skagerrak) Sources: See Supplementary Table 1

Discussion

This literature review of life-history traits of 14 batoid species of the NEA and the MED identified important information gaps for specific life-history traits parameters (e.g., lifespan and natural mortality) and certain skate species (sandy, shagreen and small-eyed ray). The PCA and cluster analysis hint at the presence of at least two groups of species that differ in their length at infinity (L), sizes (Lm), age (Agem) at maturity and growth rate (K), and to a lesser degree in their development rates and trophic level. The intraspecific analysis was severely limited to particular species and life-history traits, nonetheless it provides an important overview of the available data and data gaps. The detection of intraspecific patterns of variation between region and sub-regions is not possible at this level of data availability.

Systematic literature review

The natural abundance and density batoids in the NEA and MED impacts the quantity of available data on life-history traits for each species. In the North Sea the starry and thornback ray are caught frequently and in higher number in fisheries and scientific hauls compared to the less frequent and less abundant sandy, shagreen, and small-eyed rays, consequently increasing the number of studies performed on the species and values reported for life-history traits. Although the distribution range of the species and fishing/survey efforts will also impact the probability of encountering (e.g., fishing, sampling) a species in particular regions based (Froese and Pauly 2019) relevant data gaps were found in regions and sub-regions were a species is known to naturally distribute. Herein, even for some of the most data-available species, some areas of their distribution have not be analysed so far, to provide life-history traits values. For instance, the thornback ray, distributed from the North Sea to the Mediterranean Sea and Black Sea, however no life-history traits values were found from studies coming from the Bay of Biscay. This denotes the need of more focused research on batoids in both regions.

Together with the characteristics of their populations (i.e., misidentification may explain the lack of available data for the rare and less caught sandy, shagreen and small-eyed rays. Shagreen rays appear to have been mistakenly reported as white skate in the NEA (ICES 2020a, b, c) and as sandy ray in the MED (McCully and Walls 2015). Furthermore, sandy ray landings reported in Belgium are unreliable as they are suspected to refer to small-eyed rays (ICES 2020a, b, c). However, misidentification seems to not only be limited to rare species. For instance, blonde rays have been confused with spotted and speckled rays (Raja polystigma) in the Mediterranean Sea (Bertozzi et al. 2003), while spotted rays can be confused with blonde, thornback and cuckoo rays (Ellis et al. 2007) and even thornback rays have been confounded with starry rays (Williams 2007). Although misidentification may not be the sole cause of the lack of data for batoids in the NEA and MED, it reduces the already low quantity of available data for rare species, while simultaneously impacting the accuracy of the estimated parameters for all species involved.

While most life-history traits were found to be commonly reported, natural mortality and lifespan estimates were reported for only three and five of the species analysed in our study, respectively. Because these estimations require substantial amounts of data, they are difficult to calculate (Simpfendorfer et al. 2005). Furthermore, lifespan has commonly been estimated from in-captivity individuals which may or not represent its true value for wild individuals. Batoids stock assessments in the NEA and MED are made within the framework of data-limited stocks, and limited to survey catch rates indicators. The estimation of the natural mortality, together with size and age at maturity, are key traits to progress toward analytical population assessment models.

Exploratory analysis

Interspecific analysis

Size is highly correlated with the onset of maturity in fish (Rochet 2000; Félix et al. 2011). Our review showed that larger batoids species (e.g., white skate and common skate complex) grow slower and reach maturity later and at larger sizes than smaller-sized, fast-growing species (e.g. cuckoo, spotted ray). Additionally, fish size and its growth rate are intrinsically linked to its metabolic rate. Slow metabolism species usually reach larger sizes as they do it a slower growth rate, while high metabolic species grow faster, but will be limited to smaller sizes due to the inherent limitations linked to higher metabolic activities (Clarke and Johnston 1999; Urbina and Glover 2013; Rosenfeld et al. 2015). In conclusion the two identified groups of batoids are very likely to be impacted in different ways by fishing (i.e. landings and discards), other sources of anthropogenic disturbance (i.e. habitat destruction, pollution) and environmental change (i.e. climate change).

Faster growth may be intrinsic to particular species, but is also usually related to higher metabolic rates when living in a warmer environment (Jobling 1993). Lower temperatures tend to limit the energy available for reproductive growth, as most energy is directed to ensure somatic growth (Ware 1980; Roff 1983), hence slowing down growth and delaying maturity. Water temperatures can therefore lead to an increase/reduction of the vulnerability of a species to perturbance. For instance, thornback and marbled electric ray have been found to register similar life-history traits. However, the generally colder water where the thornbacks are found (i.e. North Sea, Celtic Se) more frequently are very likely to impact on a different way the species and its resilience compared to the impact of warmer Mediterranean waters where on the latter species. Nonetheless, as it is for most fish species, although important, the temperature is not the only relevant parameter shaping life-history traits in batoids, but rather the whole ecological niche and evolutionary history which will dictate particular life-history traits and adaptations (Bijlsma and Loeschcke 2005; Schulte 2014).

The type of diet is also linked to the capacity to supply nutrients and energy, hence to grow at a faster or slower pace. Species that feed predominantly on sessile organisms tend to grow faster than ichthyophagous (Felix et al. 2011). The batoids analysed in this study vary from carnivores that feed on herbivores (3.5 ≤ TL < 4), such as sandy, shagreen and undulate ray, to carnivores that feed on other carnivores (TL ≥ 4), such as common stingray, electric ray and white skate (e.g. Last et al. 2016, Froese and Pauly 2019). The trophic level value makes references to the diet of adult individuals, while juveniles are very unlikely to feed on the same prey. Herein, most batoids analysed in this study prey in fish when adults (e.g., common stingray, electric rays and white skate), but depend on benthic invertebrates during their growth and continue to do so even after reaching maturity (Ebert and Stehmann 2013). For instance, in the Adriatic Sea, small sized thornback rays (< 25 cm Lt) tended to feed primarily on small crustaceans, while larger specimens preferred to feed on larger decapods, cephalopods and teleost (Šantić et al. 2012). Such clear ontogeny dietary shifts have been detected for several batoids species during throughout their life cycle (e.g. Farias et al. 2006; Mulas et al. 2019; Brown-Vuillemin et al. 2020), complicating the usage of indicators, such as the trophic level, to compare species. Herein specific life-stage dietary preferences may represent different limitations and/or advantages when facing sources of disturbance.

Life-history strategies have been identified as the results of from trade-offs between different traits (Pecuchet et al. 2017). Batoid life-history traits tend to show an equilibrium strategy (Pecuchet et al. 2017) and are characteristic for species with relatively large body sizes, late maturation, low fecundity, long development times and high trophic levels, making them highly susceptible to be over-fished as their stocks need more time to recover. Although sharing a common life strategy, our study highlighted clear differences in the growth (length at infinity and growth rate) and reproductive (age and length at maturity) traits between batoids species of the NEA and MED making clear that all species may not be seen as equal in term of their vulnerability. These differences appear to follow the phylogeny of the studied species, as the genera Rostroraja and Dipturus were found to be outliers that distinguished themselves based on their larger sizes and slower growth rates. In contrast, species belonging to the genera Raja, Amblyraja, and Leucoraja share intermediate to smaller sizes, faster growth and attainment of maturity at smaller sizes and earlier ages. Grouping batoids species by genus has been suggested for management considerations, as they may share a similar susceptibility to individual (post-release survival) and population (overfishing) disturbances (STECF 2017). However, the identification of these groups should not be considered as strict rules on how to manage species, but rather as a species vulnerability gradient (Walker and Hislop 1998). For instance, although a member of the Raja genus, the blonde ray has been suggested to be more susceptible to overfishing due their particularly large length at infinity (Porcu et al. 2015), correlating to its position on the PCA biplot.

Slow-growing species, such as the common skate complex and white skate, reach maturity at later ages and will therefore take longer to reproduce and produce viable recruits, making the species less capable of recovering from intensive fishing and discard, and therefore easily overexploited (e.g. Ryland and Ajayi 1984; Dulvy 2002; Iglésias et al. 2010), justifying their protection status and fisheries ban in European waters (Council of the European Union 2019, 2020). Faster growing, early maturing and generally smaller size species (e.g. Raja, Leucoraja) may be less sensitive (but not resilient, Pinsky et al. 2011) to fishing, as their populations tend to recover faster from exploitation, but face higher levels of stress when exposed to the air, due to their generally higher metabolic rates (Brown et al. 1993; Nilsson and Oslund-Nilsson 2008). Considering that all batoids species have been exempted from the landing obligation in European waters (Uhlmann et al. 2019; ICES 2020a, b, c), it is important to understand how life-history traits may relate to discard survival rates, and include these estimates into the calculation of fishing mortality (catches and discards) for stock assessments. Moreover, as survival studies are logistically complex, costly and long in duration, meaning that most species, life-history traits may provide a handy cost-effective proxy to identify species more/less likely to die once thrown back at the sea.

With few exemptions, batoids caught in European fisheries in the NEA are managed under a common TAC (ICES 2020a, b, c). The differences found between species are caused by differences in their growth and reproductive strategies, which is inherently linked to different degrees of vulnerability. This confirms the need for differential measures (i.e., TAC, maximum allowed by-catch, the prohibition of certain gear, gear selectivity, spatiotemporal avoidance) to apply for these species (as they currently exist). Reproductive traits, such as the length at maturity, are already being used to provide relevant information to continuously update minimum landings sizes. However, it is the subject of further research to conceptualize a way to include them into species-specific stock assessments. Incorporating a measure of the species vulnerability through the usage of life-history traits, such as the growth rate, could increase the realism of stock estimation and the model’s predictive ability (Gislason et al. 2010; Le Quesne and Jennings 2012). Length-based indicators, such as the proportion of fish in the catch that are larger than length at maturity (Lmat), have already been proposed to help stock assessment estimations of data-limited fisheries (e.g. Shephard et al. 2018). In conclusion, for batoids in the NEA and MED, life-history traits can and should be used as indicators and/or proxies while a robust data framework is still “on the making”.

High definition data on life-history traits of batoids in the NEA and MED is very unlikely to be sampled, compiled and made available for researcher soon. Nonetheless, even when data on key parameters needed for stock assessment (age, growth, maturity) is absent, the correlation between life-history traits parameters could be used to identify proxies and indicators, hence providing a better understanding on the vulnerability of data-limited species (Dulvy and Reynolds 2002; Hutchings et al. 2012). For instance, when measures of length at maturity, used in stock assessment models and in the definition of management measures (i.e. minimum landing size), are not available for a particular species and/or region, measures of the highly correlated length at infinity could be used as proxy. This will also be highly valuable when looking to corroborate/validate historical parameter estimation based on old and now considered obsolete, ageing methods (Campana 2014).

Intraspecific analysis

Plasticity in life-history traits across distribution ranges is common among elasmobranchs (e.g., Templeman 1987; Lombardi-Carlson et al. 2003; Neer and Thompson 2005; Licandeo and Cerna 2007; Bradley et al. 2017) and teleosts (e.g., Deniel 1990; Pauly 1994; Vinagre et al. 2009; Félix et al. 2011) and is generally attributed to environmental gradients (e.g. temperature, food sources) (Félix et al 2011) and anthropogenic disturbance (e.g. fisheries) (Devine et al. 2012; Claireaux et al. 2018) and is key to manage species stocks, as it can have important implications for stock assessment and management. Not being able to detect such variation can lead to relevant deviation in the estimation of stock dynamics, projections and reference points (Lorenzen 2016). Different batoids stocks have already been identified and are being used by the ICES to provide their advices (ICES 2021). The division and/or union of a species population into several fishing stock has to be made following the population structure of the species itself. However, this information is not always available, leading to questionable delimitations. Instead of building stock delimitations from fisheries patterns and characteristics, life-history traits can, in the absence of genetic and population dynamic studies, provide key understanding on the population structure of the species. For instance, the recompilation of life-history trait data could help to identify whether the stock delimited for cuckoo ray in the West of Scotland, southern Celtic Seas, western English Channel and Bay of Biscay (ICES 2020c) makes ecological sense, or like for most other species, if it should be split into smaller management units.

However for batoids in the NEA and MED the currently available data prevented us to analyse this. In addition to the lack of data from particular species, life-history traits and sub-regions, the differences in methodologies used in several decades of studies make straight forward comparisons impossible. Even when data is widely available (i.e. thornback ray), studies performed in the same sub-region were often found to deviate from each other. Although this differences may be caused by finer spatial plasticity (i.e. sea basins, gulfs) in batoids communities, it may also be the result of a the recompilation of a multitude of studies with non-homogeneous methodologies. Especially when comparing values between sub-regions in which the availability of data can be very different, average values have to be taken with care, as they will be influenced by a single low or high value in short series of data. Different sampling design (methodology per se, number of samples), year and season of the sampling, experience of the research team and the accuracy and precision of the reported parameters (i.e. standard deviation/error associated to the sample size) will also impact the quality of the obtained data.

Considering this, trying to explain why and how the registered values vary (or not) between regions and sub-regions would be prone to error and miss-interpretation. For instance, the length at infinity values reported for the thornback ray in the Celtic Sea by Holden (1972; Linf♀ = 63.3, Linf♂ 49.5) are smaller than the average of the remaining reported values performed in this sub-region (Linf♀ = 120.1, Linf♂ 103.1). However, there is no evidence to say that the former study and value, is less valid than the rest. Eliminating this value based solely on its deviation from the rest would therefore be wrong. We believe that valuing and assessing the (un) certainty of each reported estimations is beyond the scope of this study, and that this task, should be carried by a larger research team (i.e. ICES-WG). It is therefore beyond the possibilities of our study to accurately detect and identify general patterns of life-history traits plasticity across the studied area.

The compiled data revealed the presence of important regional data gaps for most of the analysed species. It is highly probable that this is due to the spatial limitation of scientific survey and sampling campaigns. Although, sampling batoids all across their distribution may be a complicated and costly task, obtaining a full picture of the variation of life-history trait across the species distribution is key to manage specific stocks. Hence we, as first priority, recommend to target future research to fulfil significant regional data gaps (i.e. Bay of Biscay and English Channel). It is also urgent to homogenize the methodologies for collecting and estimating life history parameters to facilitate multiple cross-study comparisons and the detection of general plasticity patterns in space and time. Finally, once an homogeneous picture of life-history traits is obtained for the studied area, the links between environmental (i.e. climate change) and anthropogenic factors (i.e. fisheries pressure) and trait plasticity could be established. This would provide a tool for sustainable management across different sub-regions and to predict how future scenarios will affect batoids fisheries stocks.

Authors' contributions (optional: please review the submission guidelines from the journal whether statements are mandatory)

DV, NVB, and SSU conceived the publication, DV, screened the literature, explored, analyzed the data, and made the figures under the supervision of NVB, SSU, and BA. DV wrote the manuscript, while PW, NVB, SSU, and BH provided critical additions and suggestions to the manuscript form and content. All authors edited the manuscript and agreed to its publication in its current form.