Genetic Demography of the Blue and Red Shrimp, Aristeus antennatus: A Female-Based Case Study Integrating Multilocus Genotyping and Morphometric Data

In this study, we quantified the three key biological processes, growth, recruitment, and dispersal pattern, which are necessary for a better understanding of the population dynamics of the blue and red shrimp Aristeus antennatus. This marine exploited crustacean shows sex-related distribution along the water column, being females predominate in the middle slope. The present study attempts to fill the existing gap in the females’ genetic demography, as scarce knowledge is available despite being the most abundant sex in catches. We analyzed morphometric data and genotyped 12 microsatellite loci in 665 A. antennatus females collected in two consecutive seasons, winter and summer 2016, at the main Mediterranean fishing ground as a model. Almost every female in summer was inseminated. Five modal groups were observed in both seasons, from 0+ to 4+ in winter and from 1+ to 5+ in summer. Commercial-sized sorting based on fishermen’s experience resulted in a moderate-to-high assertive method concerning cohort determination. Genetic data pointed out females’ horizontal movement between neighboring fishing grounds, explaining the low genetic divergence detected among western Mediterranean grounds. Our results could represent critical information for the future implementation of management measures to ensure long-time conservation of the A. antennatus populations.


Introduction
The blue and red shrimp, Aristeus antennatus (A. antennatus) (Risso, 1816) (Crustacea, Aristeidae) (En-Blue and red shrimp; Fr-Crevette rouge; Sp-Gamba roja) is the eurybathic species that achieves the broadest vertical distribution in the Mediterranean Sea with reported records from 80 down to around 3000 m depth [1,2]. However, the species shows a demographic structure along the water column. The proportion of males and juveniles tends to increase with depth, whereas females are more abundant in the middle slope [3,4], with peaks of abundance around 700 m [1]. The mating activity starts from late winter to late spring (January-May) when adults of both sexes aggregate at 400-900 m depth, and after mating, males return to deeper waters [4,5]. Although females' maturity begins in spring, they can store sperm to fertilize the eggs in the absence of males [5,6]. Hence, during the spawning season, shoals are composed mainly of mature adult females [4,7,8]. Then, in autumn and winter, there is a decrease in population density in the middle slope as reproductive shoals break up [9].
Larval dispersal by marine currents during the A. antennatus pelagic stage has been pointed out as the main responsible for the high level of genetic connectivity found across was covered by a codend cover with a mesh size of 12 mm stretched (see Abras et al. [20] for further details). As the codend cover is prohibited for commercial hauls, it was specifically authorized for these two scientific campaigns.
Immediately after being captured, the crew sorted by eye all the individuals caught in the codend into four commercial sizes (small, medium, large and extra-large) according to their expertise. A random sample of around 100 adult females from each size group was kept for subsequent analyses, but due to limited availability, all captured females were retained for the largest sizes (Table 1). Females were identified through the presence of thelycum (morphological criteria [33]). In addition, one or more spermatophores adhered to the thelycum were detected in summer individuals, as July corresponds to the peak of the spawning season [19]. In winter, females in the codend cover were identified on board by members of the scientific team according to the morphological criteria. The presence of transparent gonads suggested they were juvenile females, and a sample of around 100 was randomly kept for analysis. In summer, there were no females captured in the codend cover. All kept individuals of each campaign were immediately stored on ice on board and quickly transported to the laboratory. Once there, the carapace length (CL) was measured from the orbital margin to the mid posterior edge of the carapace using a digital vernier caliper, and a piece of approximately 10 mg of the muscle of each individual was preserved in 70% ethanol and stored at room temperature for further genetic analyses. Table 1. Commercial size groups of Aristeus antennatus (A. antennatus) females, including sampling season, commercial category, mean carapace length (CL) with standard deviation (SD), number of females (N), number of females with one or more spermatophores adhered to the thelycum, and estimates of genetic diversity (excluding locus Aa818).

DNA Extraction and Molecular Markers
Genomic DNA was extracted from the muscle tissue of each female following the phenol-chloroform protocol adjusted to A. antennatus by Fernández et al. [41]. Twelve polymorphic microsatellite loci developed for this species were amplified by means of three multiplex PCRs and one PCR singleplex: Aa123, Aa1255, Aa138, Aa1444, Aa496b, Aa667, Aa681, Aa751, Aa956, Aa1061, Aa1195, and Aa818 [19,42]. PCR products were analyzed in an ABI PRISM 3730 × l sequencer (Applied Biosystems, Foster City, CA, USA) at the Sequencing Unit of the Universidade de Santiago de Compostela (Campus Lugo, Lugo, Spain). Genotyping was carried out using GeneMapper 4.0 software with GeneScan TM 500 LIZ ®® size standard (Applied Biosystems).

Size Frequency Distributions, Cohort Identification, and Comparison with Commercial Categories
The CL-frequency distributions in 1 mm class intervals of the females in winter and summer were studied through a modal progression analysis (MPA) implemented in the FISAT II software [43]. The Bhattacharya method [44] was used for the visual identification of the different normal distributed size groups in the polymodal CL-frequency distributions. For each season, results were displayed as the computed mean of each size group with its standard deviation (SD) and a separation index (SI). The latter should be ≥2 between consecutive size groups to assume each one as a single age group; otherwise, it means substantial overlap between cohorts. Specific cohorts in our samples (0+, 1+, 2+, 3+, 4+, and 5+) were assigned, combining MPA results with CL-cohort relationships previously reported for western Mediterranean A. antennatus populations.
A one-way ANOVA and Scheffe post-hoc tests included in the IBM SPSS Statistics version 25 package (Armonk, NY, USA) were performed to compare the mean CL in winter and summer of each size group present in both campaigns (small, medium, large, and extra-large). The same analysis was performed for the MPA-identified groups shared in the two sampling seasons (1+, 2+, 3+, and 4+). Cohen's kappa coefficient (K) [45], which described the level of concordance among methods relating to the observed agreement and the agreement expected by chance, was calculated between the groups of females identified by the MPA and the categories derived from the commercial criteria, including juveniles, in each season. The interpretation of K values was as follows: 0 to 0.20 indicates slight agreement, 0.21 to 0.40 fair agreement, 0.41 to 0.60 moderate agreement or good agreement beyond chance, 0.61 to 0.80 substantial agreement, and 0.81 to 1 almost perfect agreement [45]. Calculations were performed with VassarStats, which is available online at http://www.vassarstats.net/kappa.html (accessed on 23 May 2022).

Genetic Analysis
Genetic diversity at each female commercial size group (including winter juveniles caught in the codend cover) as well as at each age group was estimated as the number of alleles per locus (N A ) and allelic richness (A R ) calculated using FSTAT 2.9.3 [46], and observed (H O ) and expected (H E ) heterozygosities using GENEPOP 4.7.0 [47]. Deviations from Hardy-Weinberg equilibrium (HWE) were summarized using the inbreeding coefficient F IS and tested using the exact probability test of Guo & Thompson [48], implemented in GENEPOP 4.7.0. Bonferroni correction was used to adjust significance levels. Null allele frequency for each locus was assessed using Micro-Checker 2.2.3 [49]. Genetic distinctions between female commercial size groups and also between age groups in winter and summer were assessed by pairwise F ST using the program FSTAT 2.9.3 [46].
As previously done with A. antennatus males in Abras et al. [20], the parental contribution of 2015 local spawners to the 2016 juvenile and small female groups were analyzed using assignment tests. Putative parental individuals from 2015 included adult females (F; N = 52; mean CL = 38.27 ± 3.76 mm), sympatric males captured with females (M; N = 61; mean CL = 22.78 ± 2.59 mm), and mating males inferred from spermatophores attached to females' thelycum (S; N = 59) reported by Planella et al. [19]. HYBRIDLAB 1.0 software [50] was used to generate independent F1 offsprings from these 2015 parent genotypes involving mating between females and males (F × M F1) and between females and spermatophore genotypes (F × S F1). Ten sets of paired baselines of 100 simulated offspring individuals from the resulting females and spermatophores (F × S F1 baseline) and 100 individuals from the admixture of females and males (F × M F1 baseline) were used in subsequent assignations. Thus, 25 more independent sets of 100 individuals for each combination (F × S F1 and F × M F1) were also simulated and independently assigned to each of the 10 baseline sets, producing a total of 250 replicates to evaluate the accuracy of the assignments to simulated baselines. Thereafter, the 2016 female groups (winter juveniles, winter small, and summer small) were assigned to each of the 10 baselines. Assignments were carried out using the Bayesian Rannala and Mountain method [51], as implemented in GENECLASS 2 [52]. Individuals with an assignment probability < 0.01 to both baselines were considered as having another source. Comparisons of the results of these latter assignations with those from simulated data were performed using one-way ANOVA and Scheffe post-hoc tests included in IBM SPSS Statistics.
In addition, the geographical origin of juvenile and small females was analyzed through assignment tests. Winter juveniles and small females were used as a baseline together with two sets of small females collected during the same winter of 2016 at two nearby fishing grounds: one from Roses (Cap de Creus Canyon, n = 55) situated to the north of Palamós, and the second from Blanes (Blanes Canyon, n = 54), located to the south (Figure 1) [20]. Chi-square tests using VassarStats were computed to compare the assignment pattern of the summer females to winter baselines. Finally, assignment tests were also used to obtain the distribution to winter baselines of female summer groups up to 3+ derived from the MPA. Fisher exact probability tests were used to compare summer assignment patterns to winter baselines. in GENECLASS 2 [52]. Individuals with an assignment probability < 0.01 to both b were considered as having another source. Comparisons of the results of these l signations with those from simulated data were performed using one-way ANO Scheffe post-hoc tests included in IBM SPSS Statistics.
In addition, the geographical origin of juvenile and small females was a through assignment tests. Winter juveniles and small females were used as a bas gether with two sets of small females collected during the same winter of 2016 nearby fishing grounds: one from Roses (Cap de Creus Canyon, n = 55) situate north of Palamós, and the second from Blanes (Blanes Canyon, n = 54), located to t ( Figure 1) [20]. Chi-square tests using VassarStats were computed to compare the ment pattern of the summer females to winter baselines. Finally, assignment te also used to obtain the distribution to winter baselines of female summer groups derived from the MPA. Fisher exact probability tests were used to compare sum signment patterns to winter baselines.

Size-Age Composition of A. antennatus Females
Mean CL was significantly lower in winter than in summer for all commercia size groups (juveniles from codend cover removed) but the medium-sized catego ure 2a). Practically all the females caught in summer were inseminated (331/334 and a small proportion of them (27/331; 8.2%) presented more than one sperma (Table 1)

Size-Age Composition of A. antennatus Females
Mean CL was significantly lower in winter than in summer for all commercial female size groups (juveniles from codend cover removed) but the medium-sized category ( Figure 2a). Practically all the females caught in summer were inseminated (331/334; 99.1%) and a small proportion of them (27/331; 8.2%) presented more than one spermatophore ( Table 1).
The MPA from CL-frequency distributions showed that females in both winter and summer were composed of five separated (SI > 2) modal groups (Table 2, Figure 3a,b). Based on CL-age keys reported in previous studies on A. antennatus western Mediterranean fisheries [23,[30][31][32]34,53,54] (Table 3), the age of each MPA group was assigned to cohorts from 0+ to 4+ in winter and from 1+ to 5+ in summer. The size of the first identified cohort (19.26 ± 1.77 mm) was slightly larger than the one reported by Carbonell [34] and looks similar to the smaller modal cohort indicated by D'Onghia et al. [23] for samples at the Ionian Sea collected in April and May. Though these latter authors did not use codend cover, but a stretched mesh of 20 mm in the codend and reduced proportions of collected females were in the range between 16 and 26 mm of CL, especially in April catches. Genes 2022, 13, x FOR PEER REVIEW 6 of 18 The MPA from CL-frequency distributions showed that females in both winter and summer were composed of five separated (SI > 2) modal groups (Table 2, Figure 3a,b). Based on CL-age keys reported in previous studies on A. antennatus western Mediterranean fisheries [23,[30][31][32]34,53,54] (Table 3), the age of each MPA group was assigned to cohorts from 0+ to 4+ in winter and from 1+ to 5+ in summer. The size of the first identified cohort (19.26 ± 1.77 mm) was slightly larger than the one reported by Carbonell [34] and looks similar to the smaller modal cohort indicated by D'Onghia et al. [23] for samples at the Ionian Sea collected in April and May. Though these latter authors did not use codend cover, but a stretched mesh of 20 mm in the codend and reduced proportions of collected females were in the range between 16 and 26 mm of CL, especially in April catches. Table 2. Description of the groups identified through the modal progression analysis of the carapace length (CL) frequency distributions of A. antennatus females, including sampling season, modal and age group, mean carapace length with standard deviation (SD), separation index (SI), number of females (N), and estimates of genetic diversity (excluding locus Aa818).   The comparison of the mean CL in winter and summer of each age group present in both campaigns (1+, 2+, 3+, and 4+) resulted in lower values in winter than in summer for all age groups except for the cohort 2+ but only significant in cohorts 1+ and 3+. In the case of the age group 2+ occurred the opposite; the mean CL was significantly higher in winter than in summer (Figure 2b).

Modal
All females in the range of mean ± SD of each modal group was clearly considered as a member of the cohort and retained to have better baselines for individuals to make comparisons between cohorts and commercial sizes. Using this strategy, 77 out of 331 females captured in winter (23.3%) and 79 out of 334 of those caught in summer (23.7%) remained doubtful and were removed from cohort groups ( Table 2). The agreement between these cohort baselines and the commercial sizes was substantial according to Cohen's kappa either in winter samples (K = 0.75; 95% CI = 0.68-0.81) or in summer ones (K = 0.73; 95% CI = 0.67-0.80) (Table S1). Thus, each cohort baseline was mainly composed of females assigned to a single commercial size by the crew vessel, and only cohorts 3+ in winter and 4+ in summer were clear admixtures between two commercial size groups (Figure 3c,d).
The results indicated that the medium commercial size in winter is composed of individuals from the 2+ and 3+ cohorts, and the large and extra-large summer females are also admixtures of two consecutive cohorts (Figure 3a,b). The presence of cohort 0+ females in the small size group retained by the codend in winter indicated that a portion of these young individuals are captured even in a codend of a mesh size of 40 mm. The distribution of CL between each commercial category and the more related MPA group in winter is fairly coincident except for the extra-large females in summer that showed a mean CL value between the two older cohorts (Figure 4). All females in the range of mean ± SD of each modal group was clearly considered as a member of the cohort and retained to have better baselines for individuals to make comparisons between cohorts and commercial sizes. Using this strategy, 77 out of 331 females captured in winter (23.3%) and 79 out of 334 of those caught in summer (23.7%) remained doubtful and were removed from cohort groups ( Table 2). The agreement between these cohort baselines and the commercial sizes was substantial according to Cohen's kappa either in winter samples (K = 0.75; 95% CI = 0.68-0.81) or in summer ones (K = 0.73; 95% CI = 0.67-0.80) (Table S1). Thus, each cohort baseline was mainly composed of females assigned to a single commercial size by the crew vessel, and only cohorts 3+ in winter and 4+ in summer were clear admixtures between two commercial size groups (Figure 3c,d).
The results indicated that the medium commercial size in winter is composed of individuals from the 2+ and 3+ cohorts, and the large and extra-large summer females are also admixtures of two consecutive cohorts (Figure 3a,b). The presence of cohort 0+ females in the small size group retained by the codend in winter indicated that a portion of these young individuals are captured even in a codend of a mesh size of 40 mm. The distribution of CL between each commercial category and the more related MPA group in winter is fairly coincident except for the extra-large females in summer that showed a mean CL value between the two older cohorts (Figure 4).

Genetic Diversity and Genetic Divergence of A. antennatus Females
A total of 665 A. antennatus females were genotyped for 12 microsatellite loci (Table S2). Only locus Aa818 was ruled out from all subsequent analyses due to a null allele frequency > 0.2 in four commercial categories (Table S2) and six age groups (Table S3).
Regarding commercial categories (Table 1), females with medium size reached the maximum mean A R (5.34, Table 1) among winter samples; meanwhile, the small size group displayed the highest H O and H E (0.492 and 0.644, respectively). The females from the large and extra-large size groups showed the lowest diversity values and were the only ones displaying significant departures from HWE (F IS = 0.284 and 0.303, respectively, Table 1). In summer, diversity indexes were fairly similar to winter ones, but small females achieved the highest value of A R (5.40) and large females showed the highest H O and H E (0.511 and 0.633, respectively). A significant departure from HWE was observed in the small, large, and extra-large females in summer (Table 1). These deviations were related to positive F IS values, indicating a deficit of heterozygosity, and affect commercial sizes composed of more than a cohort (Figure 3b).

Genetic Diversity and Genetic Divergence of A. antennatus Females
A total of 665 A. antennatus females were genotyped for 12 microsatellite loci (Table  S2). Only locus Aa818 was ruled out from all subsequent analyses due to a null allele frequency > 0.2 in four commercial categories (Table S2) and six age groups (Table S3).
Regarding commercial categories (Table 1), females with medium size reached the maximum mean AR (5.34, Table 1) among winter samples; meanwhile, the small size group displayed the highest HO and HE (0.492 and 0.644, respectively). The females from the large and extra-large size groups showed the lowest diversity values and were the only ones displaying significant departures from HWE (FIS = 0.284 and 0.303, respectively, Table 1). In summer, diversity indexes were fairly similar to winter ones, but small females achieved the highest value of AR (5.40) and large females showed the highest HO and HE (0.511 and 0.633, respectively). A significant departure from HWE was observed in the small, large, and extra-large females in summer (Table 1). These deviations were related to positive FIS values, indicating a deficit of heterozygosity, and affect commercial sizes composed of more than a cohort (Figure 3b).
In relation to age groups (Table 2), in winter, the highest mean AR and HE were achieved by group 2+ and the maximum HO by group 1+. Significant departures from HWE genotype expectations were observed in groups 1+ and 4+. In summer, group 5+ showed the maximum result of AR and HE, whereas the highest HO was achieved by group 3+. Significant departures from HWE were observed in groups 3+ and 4+.
Pairwise FST between all-female commercial categories as well as age groups were low, ranging from 0.0000 to 0.0047 (Table 4). After Bonferroni correction, FST values between paired samples resulted in no significant differences. In relation to age groups ( Pairwise F ST between all-female commercial categories as well as age groups were low, ranging from 0.0000 to 0.0047 (Table 4). After Bonferroni correction, F ST values between paired samples resulted in no significant differences.

Origin of Juvenile and Small Females
In the assignment tests for the assessment of the parental contribution of 2015 local spawners to the 2016 juvenile and small female groups, the simulated individuals indicated a limited accuracy of baselines (Table 5), probably due to the low but significant divergence between males and spermatophores [19]. However, in both cases, around 60% of the simulated F1 individuals were assigned to the correct baseline (64.5% for F × S; 58.6% for F × M). Winter juvenile females captured in the codend cover and small winter and summer females showed similar assignment patterns showing a higher percentage of assignation to the F × S F1 baseline than to the F × M F1 one. The pattern for capturing females was also distinct from simulated F1 offsprings from local spawners in the larger proportion (>20%) of individuals assigned to other sources (Table 5 and Figure 5).  Assignment tests to analyze the geographical origin of juvenile showed a limited accuracy of the baselines, especially in the case of the Palamós collected in winter ( Figure 6). Even so, more than 50% of the f were correctly assigned to their baseline. As previously observed in cruitment from other grounds was confirmed by the high proportio signed to Roses, Blanes, and other sources, reaching 35% of the total a mer small females. Small summer females showed an assignment p Assignment tests to analyze the geographical origin of juvenile and small females showed a limited accuracy of the baselines, especially in the case of the small females from Palamós collected in winter ( Figure 6). Even so, more than 50% of the females in this group were correctly assigned to their baseline. As previously observed in males [20], the recruitment from other grounds was confirmed by the high proportion of individuals assigned to Roses, Blanes, and other sources, reaching 35% of the total assignments in summer small females. Small summer females showed an assignment pattern significantly different from winter juveniles (p = 0.02) and winter small (p < 0.0001) but more related to winter juveniles (44% of individuals assigned to this baseline) ( Figure 6).
Finally, the assignment distribution of female MPA groups up to 3+ to their closest winter baselines showed that, although summer cohorts had assignment patterns significantly different from winter baselines in all cases (Fisher exact test p < 0.005), they were assigned in a higher proportion to their predecessor year cohort in winter (Figure 7). However, summer cohorts had a percentage of contribution of their same year cohort in winter (from 14.3 in summer females 3+ to 28.9% in summer females 2+) and also input from other sources (up to 22.4% in summer females 3+) (Figure 7). spring sets and the 2016 A. antennatus female groups. Continuous lines indicate the median value of each group and dashed lines represent the quartiles. Letters (a, b) indicate groups identified by the post-hoc Scheffe test. F, females; S, spermatophores; M, males. Assignment tests to analyze the geographical origin of juvenile and small females showed a limited accuracy of the baselines, especially in the case of the small females from Palamós collected in winter ( Figure 6). Even so, more than 50% of the females in this group were correctly assigned to their baseline. As previously observed in males [20], the recruitment from other grounds was confirmed by the high proportion of individuals assigned to Roses, Blanes, and other sources, reaching 35% of the total assignments in summer small females. Small summer females showed an assignment pattern significantly different from winter juveniles (p = 0.02) and winter small (p < 0.0001) but more related to winter juveniles (44% of individuals assigned to this baseline) ( Figure 6). Finally, the assignment distribution of female MPA groups up to 3+ to their closest winter baselines showed that, although summer cohorts had assignment patterns significantly different from winter baselines in all cases (Fisher exact test p < 0.005), they were assigned in a higher proportion to their predecessor year cohort in winter (Figure 7). However, summer cohorts had a percentage of contribution of their same year cohort in winter (from 14.3 in summer females 3+ to 28.9% in summer females 2+) and also input from other sources (up to 22.4% in summer females 3+) (Figure 7).

Size-Age Composition and Commercial Categories of A. antennatus Females
As reported by Demestre & Lleonart [25], females and males should be analyzed separately because they have different demographic patterns: (i) Females and males show ontogenic differential distribution [55]. (ii) Females live longer than males [25]. (iii) Females are more abundant in catches [25]. (iv) The size of the catches is larger in females, indicating the size dimorphism of the species [25,56], which is consistent with both the present study and the previous one focused on males carried out during the same sampling campaigns [20]. Small females were bigger (mean CL = 25.26 ± 3.39 mm) than males in the same sampling (mean CL = 21.92 ± 2.83 mm) [20] (t-test p < 0.001). (v) Females have greater age diversity and size composition [55]. (vi) Both sexes become sexually active within the first year of their life [5,56], but females show seasonal maturity from May to October, with a peak in July-August, whereas males are mature during the whole year [33]. Size at first maturity (CL 50 ) can be influenced by environmental factors and fluctuate between geographical areas [55]. Sardà & Demestre [33] reported a CL 50 of 27 mm in females, with some mature females from 23 mm in the Catalan Sea. Demestre & Martín [57] calculated a CL 50 of 26 mm in the same area. Carbonell et al. [55] also reported a CL 50 of 26 mm in females in the Balearic Sea whereas García-Rodríguez & Esteban [56] reported a CL 50 of 21.9 mm in the same place. In this study, the mean CL in the small size category was 22.96 mm in winter and 27.65 mm in summer, with most females inseminated (Table 1), being consistent with the aforementioned studies.
Related to the MPA analysis (Table 2), the number of estimated age groups and the size of each group was quite similar to data reported by other studies conducted close to our area of study (Table 3), i.e., Gulf of Alicante, Gulf of Vera and Balearic Islands, [34,53,54]. However, our data do not differ substantially from other studies conducted in more remote areas either, i.e., Ligurian, Tyrrhenian, and Ionian Seas [23,[30][31][32]. Indeed, the largest size-age ratio differences between groups were observed in older age groups, approximately from 4+ onwards, in which it becomes difficult to distinguish the cohorts due to age-related reduction of growth rate [34]. These differences could be explained by several reasons: (i) Oceanographic conditions are different in the areas analyzed with a broad range of variables exerting influence, such as climate, hydrodynamics, trophic resources, and predators [18,31]. The species must adapt to the different hydrological conditions and environmental factors [1,58]. (ii) The lack of an appropriate experimental method to estimate growth parameters. Although other methods have been proposed for age determination, such as the quantification of the age-pigment lipofuscin or tagging studies [32,54], further investigation is needed to prove their efficiency. Most studies make age estimations based exclusively on sizes, so results could be slightly biased by the influence of molting activity [59], fishing habits [31], and the exploitation pattern, which brings increased pressure on females [34]. (iii) The widespread recruitment throughout the Mediterranean Sea. Growth taxes could differ between individuals that recruit on the lower slope (slower growth) and the ones doing it in shallower waters (faster growth). Thus, they could not be distinguished when joining since individuals of the same age would reach slightly different sizes [23]. Data reported by Carbonell [34] estimated a monthly growth rate of 2 mm for juveniles and 1 mm for adults. Nevertheless, analysis of growth in crustaceans is challenging because it is a discontinuous process carried out by molting [60], so the comparison between studies and areas is extremely difficult and not always feasible.
A. antennatus individuals are sorted after fishing by the fisherman according to commercial categories, which are not specifically based on quantitative criteria and do not depend on cohorts. On the other hand, the MPA allows an objective way of sorting, but it is not easy to identify the different groups in the polymodal CL-frequency distributions, especially in the largest sizes corresponding to the oldest cohorts. Both methods reached a substantial agreement (K = 0.75 and 0.73 in winter and summer, respectively) being almost all age groups mainly represented by a single commercial category (contribution between 79% and 100%) with the exception of cohorts 3+ in winter and 4+ in summer, thus pointing out fishermen's sorting based on years of experience as a generally assertive method. However, the pooling of different cohorts into the same commercial category also took place, as observed in the medium size in winter and the large and extra-large categories in summer (Figure 3).

Recruitment and Geographical Origin of A. antennatus Females
Genetic data pointed out for the first time that mating males consisted of males sympatric with females and also males from other spawning grounds [19]. Furthermore, a study based exclusively on A. antennatus males detected that males sampled at fishing grounds during the spawning season have a limited contribution to the spermatophores attached to females [20]. That previous study reported through assignment tests the genetic correlation of most males collected in 2016 to the offspring of 2015 spawners with high contribution from other sources. Likewise, genetic analyses in the present work based on the same simulated F1 of individuals collected in 2015 indicate that most juvenile and small females in 2016 were also correlated to the offspring of 2015 spawners with a considerable assignment to other sources (Table 5 and Figure 5). Thus, our results confirmed a displacement of males, either vertical movement inside Palamós fishing ground or horizontal displacement from adjacent fishing grounds, as previously suggested [19,20].
Since Blanes and Roses females from winter 2016 were used in the assignment tests to analyze the geographical origin of juvenile and small females had small size (mean CL of 24.23 and 27.61 mm for Blanes and Roses, respectively) [16], only juvenile and small females in the present study could be assigned. Winter juvenile females had a higher contribution than the small ones to summer small females indicating that this latter group was mainly composed of juvenile females that have grown and females that have moved from other grounds ( Figure 6). These results were confirmed by the assignment distribution of female MPA groups up to 3+ to their closest winter baselines (Figure 7). In all cases, summer female groups were assigned by over 55% to their predecessor year cohort in winter, indicating again the growth of the winter individuals to join the next summer cohort group. However, the same year cohort in winter, as well as females from other grounds, also contributed to the group in summer, as indicated by the percentage of assignation to these sources. Data reported by Abras et al. [20] in males pointed out a limited contribution of winter male juveniles to the summer pre-adult male groups, thus differing from results on females (present study). Otherwise, in both sexes, the level of contribution from adjacent fishing grounds is high. More than 40% of males in summer groups originated from outside the Palamós ground [20], whereas the percentage of contribution from neighboring grounds to summer small females rose to 35% (present study) ( Figure 6). Our results corroborate the connectivity between nearby fishing grounds of A. antennatus in the north-western Mediterranean Sea, which likely occurs through the dispersal of different stages of the life-cycle of A. antennatus, i.e., larvae, juveniles, and adults of both sexes [16,17,20,40].
The dispersal of individuals from one ground to another involves the combination of residents and immigrants, which could lead to the low genetic differentiation (F ST ) observed with no significant variation in either commercial categories or age groups (Table 4). On the other hand, Cannas et al. [13] hypothesized that A. antennatus females could disperse more easily than males due to its main presence in intermediate currents that flow faster than deep currents. Recently, Agulló et al. [17] reported the existence of horizontal displacements of A. antennatus females and a high genetic connectivity between spawning females' grounds within the Mediterranean Subarea GSA6.

Conclusions
Combining morphometric and genetic data, this study provides new insights into A. antennatus growth, recruitment, and mobility pattern with particular emphasis on females. Commercial-sized sorting based on fishermen's experience allows a fast classification of A. antennatus captures according to market indications, which indirectly results in a moderate-to-high assertive method concerning females' cohort determination. Although, it should be kept in mind that there could be individuals belonging to more than one cohort in a commercial category. The present genetic data detail the horizontal movement of A. antennatus females between neighboring fishing grounds, which is in agreement with the high genetic connectivity pattern detected for A. antennatus populations in the north-western Mediterranean Sea. Our findings could help to improve the current protection measures and ensure the long-term sustainability of the fishery.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/genes13071186/s1, Table S1: Concordances used to calculate Cohen's kappa coefficient (K) for the comparison between commercial categories and modal progression analysis (MPA) groups in winter and summer; Table S2: Genetic diversity for twelve microsatellite loci for commercial-sized A. antennatus female groups; Table S3: Genetic diversity for twelve microsatellite loci for A. antennatus female groups identified through the modal progression analysis of the carapace length (CL) frequency distributions.