Life-history traits of batoids (Superorder Batoidea) in the Northeast Atlantic and the Mediterranean

Life-history traits provide a way to estimate the vulnerability of both individuals and populations of a species to disturbance (e.g., overexploitation, climate change). Life-history traits compilations for species of batoids in the Northeast Atlantic (NEA) and the Mediterranean Sea (MED) are scarce, outdated, and generally limited to a local or a regional scale. A literature review compiling values of 10 life-history traits describing the growth, reproductive and feeding strategies of 14 batoids in the NEA and the MED was performed. Via a principal components analysis (PCA) the main drivers of variance and (groups of) similar species were identified. Significant data gaps were revealed for natural mortality and lifespan, for most of the life-history traits of small-eyed (Raja microocellata), shagreen (Leucoraja fullonica), and sandy (Leucoraja circularis) ray, and specifically for the English Channel and Bay of Biscay. The common skate complex (Dipturus batis complex) and white (Rostroraja alba) skate were found to be different from the rest of the species due to their larger sizes, relatively slow growth, and late maturity, which may make them more vulnerable to overexploitation, while the Raja spp. and Leucoraja spp. rays tend to be smaller species growing at faster rates. Growth and reproductive life-history should be considered as vulnerability to disturbance proxies, be used to update/establish management measures (i.e. minimum size) and enhance the stock assessment predictive ability (i.e. length based-indicators). Future research should be directed towards filling important regional data gaps and providing robust estimation for unreported parameters, such as natural mortality.


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 datalimited (ICES 2012(ICES , 2020Lart 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 lifehistory 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.,  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 metaanalyses 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). Lifehistory 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 lifehistory 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.
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 (Age m ) and length at maturity (L m ) ( 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., [ 10 3 ) 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.

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 smalleyed ray; were removed from the summarized data set because of the absence of published values for age at maturity (Age m ). 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 (Age m , L m , 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 lifehistory 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. All analysis and plotting were performed in R statistical software version 4.0.2 (R Core Team 2020).

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).
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). Lifehistory 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).

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. 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; L m: 119.30-197.50 cm) and late maturation (Age m: 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  Table 2. Direct (0,1) and inverse (0,-1) correlation, are represented with blue and red tones respectively Fig. 4 Biplot of the two first principal components of the PCA of age at maturity (AgeM), length at maturity (L m ), length at infinity (L inf ), 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 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 (Age m :4.7-8.8) and smaller sizes (L m :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).

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   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 Values in bold represent the mean value for the region. Age m : Age at maturity; L m : Length at maturity; L ? : Length at infinity; K: Growth rate for males (#) and females ($).  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).
Thornback ray Both females and males thornback rays were found to in average mature at later ages in the North Sea, than in the Mediterranean Sea. However the strong deviation in age at maturity values in studies of the Mediterranean could mean that the values from both regions (i.e. NEA and MED) overlap. In the Celtic sea, female thornback rays matured at younger ages than in the North Sea, while males seemed to mature at similar ages than in the North Sea and Mediterranean. It is impossible to identify where thornback rays exhibit larger length at maturity and infinity due to the dispersion between reported values. A southward, from the NEA toward the MED, decreasing pattern can be outlined for the growth rate, while values were found to be constant within each of these regions. No studies performed in the Bay of Biscay reported values for any of the analysed life)history traits (Table 4).
Blonde ray Based on the two studies in which the age at maturity was reported, blonde rays mature earlier and at smaller sizes in the Celtic Sea than the Mediterranean Sea. Female blonde rays matured at larger length in the Mediterranean than in the NEA (i.e. North Sea and Celtic Sea), while males length seemed to not vary between the mentioned regions. In average, in NEA waters, blonde rays seem to reach larger lengths at infinity compared to the MED. The large deviation between reported values for length at infinity in the NEA could represent similar length at infinity in both regions. Similarly, in average blonde rays of the NEA (i.e. Celtic Sea) grow faster than rays of MED. Individuals of Portuguese water grow at a similar rate than in the Mediterranean sea. No values were found for the English Channel and Bay of Biscay despite being areas of its natural distribution (Table 5).
Spotted ray Spotted rays of both sexes were found to reach maturity at earlier ages and smaller lengths in the Celtic Sea than in the North Sea. No distinguishable pattern can be recognize for length at infinity between sub-regions of the NEA. Values from all sub-regions were found to be within the same range. However, in Portuguese waters, blonde rays appear to grow slightly faster, exhibiting smaller  (Table 7).
Cuckoo ray Despite distributing toward in the English Channel, Bay of Biscay, Iberian Peninsula and even in the Mediterranean Sea, cuckoo rays life history trait values were only found to have been reported for the North Sea and Celtic Sea. Although reaching maturity at a younger age, Celtic Sea cuckoo ray were found to mature at larger sizes and reached larger lengths at infinity than in the North Sea. But, the strength of this tendency is again limited by the strong data dispersion in the reported values for both parameters. In both areas males matured faster than females and males matured faster in the North Sea (Table 8).
Starry ray Northerner starry rays studied in the northern North Sea and Skagerrak sub-regions were found to mature at larger sizes than individuals from Values in bold represent the mean value for the region. Age m : Age at maturity; L m : Length at maturity; L ? : Length at infinity; K: Growth rate for males (#) and females ($). Values in bold represent the mean value for the region. Age m : Age at maturity; L m : Length at maturity; L ? : Length at infinity; K: Growth rate for males (#) and females ($).
the North Sea. No other parameter was reported for other sub-regions. Values for the length at maturity and infinity seemed to be consistent in their values between different studies as shown by the relatively low standard deviation. Contrarily, a relatively high standard deviation was found between studies reporting the growth rate of starry rays in the North Sea (Table 9).

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 (L m ), age (Age m ) 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 lifehistory 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 B TL \ 4), such as sandy, shagreen and undulate ray, to carnivores that feed on other carnivores (TL C 4), such as common stingray, electric ray and white skate (e.g. Last et al. 2016, Froese andPauly 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 lifehistory 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 speciesspecific stock assessments. Incorporating a measure of the species vulnerability through the usage of lifehistory 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 (L mat ), 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) andteleosts (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 subregions 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 missinterpretation. For instance, the length at infinity values reported for the thornback ray in the Celtic Sea by Holden (1972;L inf$ = 63.3, L inf# 49.5) are smaller than the average of the remaining reported values performed in this sub-region (L inf$ = 120.1, L inf# 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 lifehistory 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 subregions 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.
Acknowledgements Special thanks to the many scientists who indirectly contributed to this article through their reports and published studies. We would like to give special thanks to Els Torreele for her valuable advice, and to Marc Delaruelle for his time to proofread our manuscript. Finally, we very much appreciated the comments and suggestions of the anonymous reviewers and the editorial team of the journal. This study has received funding from the Interreg 2 Seas programme 2014-2020 co-funded by the European Regional Development Fund under subsidy contract No 2503-024 SUMARiS'' (Sustainable Management of rays and skates).
Funding (Information that explains whether and by whom the research was supported) This study has received funding from the Interreg 2 Seas programme 2014-2020 co-funded by the European Regional Development Fund under subsidy contract No 2503-024 SUMARiS'' (Sustainable Management of rays and skates).
Availability of data and material (data transparency) All data generated or analysed during this study are included in this published article and its supplementary information files.

Declarations
Conflict of interest (include appropriate disclosures) The authors have no conflicts of interest to declare that are relevant to the content of this article.
Human and animal participate No experiments and/or samples were taken from humans and/or animals throughout this study.

Consent to participate (include appropriate statements)
All authors agreed to participate in the elaboration of this study.

Consent for publication (include appropriate statements)
All authors edited the manuscript and agreed to its publication in its current form.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.