INTRODUCTION
⌅Most marine organisms produce a larva that develops for a certain time period in the water column. During this planktonic period, larval size and morphology are key life-history traits that can affect all aspects of performance and survival (Pettersen et al. 2015Pettersen A.K., White C.R., Marshall D.J. 2015. Why does offspring size affect performance? Integrating metabolic scaling with life-history theory. Proc. Royal Soc. B. 282: 20151946. https://doi.org/10.1098/rspb.2015.1946 ). For many taxa, bigger larvae facilitate the access of a wider prey size range, store higher endogenous reserves and avoid predators (Pepin 1989Pepin P. 1989. Predation and starvation of larval fish: A numerical experiment of size and growth-dependent survival. Biol. Oceanogr. 6: 23-44., Park et al. 2004Park S., Epifanio C.E., Grey E.K. 2004. Behavior of larval Hemigrapsus sanguineus (de Haan) in response to gravity and pressure. J. Exp. Mar. Biol. Ecol. 307: 197-206. https://doi.org/10.1016/j.jembe.2004.02.007 , Bashevkin et al. 2020Bashevkin S.M., Christy J.H., Morgan S.G. 2020. Adaptive specialization and constraint in morphological defences of planktonic larvae. Funct. Ecol. 34: 217-228. https://doi.org/10.1111/1365-2435.13464 ).
Variation in progeny size among species, and among populations within species, was initially attributed to variation in natural selection (Fox and Czesak 2000Fox C.W., Czesak M.E. 2000. Evolutionary ecology of progeny size in arthropods. Annu. Rev. Entomol. 45: 341-369. https://doi.org/10.1146/annurev.ento.45.1.341 ). However, this variability is more and more viewed as a physiological constraint that prevents mothers from producing homogeneous progeny size (Marshall et al. 2008Marshall D.J., Bonduriansky R., Bussière L.F. 2008. Offspring size variation within broods as a bet-hedging strategy in unpredictable environments. Ecology 89: 2506-2517. https://doi.org/10.1890/07-0267.1 ). In decapod crustaceans, variations in progeny size and morphology have been associated with maternal size (or age) and/or the environmental conditions experienced by the mother during the embryonic development (Anger 2001Anger K. 2001. The Biology of Decapod Crustacean Larvae. Crustacean Issues, vol 14. A. A. Balkema, Rotterdam). For instance, lower temperature at higher-latitude populations is related to the production of larger crab larvae in the Pacific sub-Arctic/Arctic sector (Chionoecetes opilio; Landeira et al. 2017Landeira J.M., Matsuno K., Yamaguchi A., et al. 2017. Abundance, development stage, and size of decapod larvae through the Bering and Chukchi Seas during summer. Polar Biol. 40: 1805-1819. https://doi.org/10.1007/s00300-017-2103-6 ), NE Pacific (Metacarcinus magister; Shirley et al.1987Shirley S.M., Shirley T.C., Rice S.D. 1987. Latitudinal variation in the Dungeness crab, Cancer magister: zoeal morphology explained by incubation temperature. Mar. Biol. 95: 371-376. https://doi.org/10.1007/BF00409567 ), SE Pacific (Cancer setosus; Weiss et al. 2010Weiss M., Thatje S., Heilmayer O. 2010. Temperature effects on zoeal morphometric traits and intraspecific variability in the hairy crab Cancer setosus across latitude. Helgol. Mar. Res. 64: 125-133. https://doi.org/10.1007/s10152-009-0173-8 ) and NE Atlantic (Macropodia rostrata; Marco-Herrero et al. 2012Marco-Herrero E., Rodríguez A., Cuesta J.A. 2012. Morphology of the larval stages of Macropodia czernjawskii (Brandt, 1880) (Decapoda, Brachyura, Inachidae) reared in the laboratory. Zootaxa 3338: 33-48. https://doi.org/10.11646/zootaxa.3338.1.2 ). It has also been observed that seasonal changes in temperature, salinity, and food availability can influence the production of different larval phenotypes throughout the year (Urzúa and Anger 2013Urzúa Á., Anger K. 2013. Seasonal variations in larval biomass and biochemical composition of brown shrimp, Crangon crangon (Decapoda, Caridea), at hatching. Helgol. Mar. Res. 67: 267-277. https://doi.org/10.1007/s10152-012-0321-4 ; González-Ortegón and Giménez 2014González-Ortegón E., Giménez L. 2014. Environmentally mediated phenotypic links and performance in larvae of a marine invertebrate. Mar. Ecol. Prog. Ser. 502: 185-195. https://doi.org/10.3354/meps10708 ). In some caridean shrimps these environmental factors are linked to maternal size, as carapace length also varies significantly between breeding seasons (González-Ortegón et al. 2018González-Ortegón E., Vay L. Le, Walton M.E.M.K., Giménez L. 2018. Maternal trophic status and offspring phenotype in a marine invertebrate. Sci. Rep. 8: 1-9. https://doi.org/10.1038/s41598-018-27709-2 ; Oliphant and Thatje 2021Oliphant A., Thatje S. 2021. Variable shrimp in variable environments: reproductive investment within Palaemon varians. Hydrobiologia 848: 469-484. https://doi.org/10.1007/s10750-020-04455-z ). Despite the growing body of evidence supporting the maternal influence in the larval phenotype, few studies have reported this effect extracting the environmental factors. Maternal influence on larval morphology has been seen in the coconut crab, Birgus latro (Sato and Suzuki 2010Sato T., Suzuki N. 2010. Female size as a determinant of larval size, weight, and survival period in the coconut crab, Birgus latro. J. Crust. Biol. 30: 624-628. https://doi.org/10.1651/10-3279.1 ), the European lobster, Homarus gammarus (Moland et al. 2010Moland E., Olsen E.M., Stenseth N.C. 2010. Maternal influences on offspring size variation and viability in wild European lobster, Homarus gammarus. Mar. Ecol. Prog. Ser. 400: 165-173. https://doi.org/10.3354/meps08397 ), the kuruma prawn, Marsupenaeus japonicus (Sato et al. 2017Sato T., Hamano K., Sugaya T., Dan S. 2017. Effects of maternal influences and timing of spawning on intraspecific variations in larval qualities of the Kuruma prawn Marsupenaeus japonicus. Mar. Biol. 164: 70. https://doi.org/10.1007/s00227-017-3118-9 ), and the blue crab, Callinectes sapidus (Caracappa and Munroe 2018Caracappa J.C., Munroe D.M. 2018. Morphological variability among broods of first-stage blue crab (Callinectes sapidus) zoeae. Biol. Bull. 235: 123-133. https://doi.org/10.1086/699922 ). However, other species such as the red king crab, Paralithodes camtschaticus, and the shore crab, Hemigrapsus crenulatus, have shown no such relationship (Swiney et al. 2013Swiney K.M., Eckert G.L., Kruse G.H. 2013. Does maternal size affect red king crab, Paralithodes camtschaticus, embryo and larval quality? J. Crust. Biol. 33: 470-480. https://doi.org/10.1163/1937240X-00002162 , Urzúa et al. 2018Urzúa Á., Bascur M., Guzmán F., Urbina M. 2018. Carry-over effects modulated by salinity during the early ontogeny of the euryhaline crab Hemigrapsus crenulatus from the Southeastern Pacific coast: Development time and carbon and energy content of offspring. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 217: 55-62. https://doi.org/10.1016/j.cbpa.2018.01.001 ), suggesting that the consistency of the pattern across the different crustacean groups is not well supported yet.
The brush-clawed shore crab, Hemigrapsus takanoi Asakura and Watanabe, 2005Asakura A., Watanabe S. 2005. Hemigrapsus takanoi new species, a sibling species of the common Japanese intertidal crab H. penicillatus (Decapoda: Brachuyra: Grapsoidea). J. Crust. Biol. 25: 279-292. https://doi.org/10.1651/C-2514 , is a common brackish and estuarine invertebrate found under intertidal boulders in its native distribution range of the northwest Pacific (Asakura and Watanabe 2005Asakura A., Watanabe S. 2005. Hemigrapsus takanoi new species, a sibling species of the common Japanese intertidal crab H. penicillatus (Decapoda: Brachuyra: Grapsoidea). J. Crust. Biol. 25: 279-292. https://doi.org/10.1651/C-2514 ). Multiple introductions, likely by shipping lines, have facilitated the expansion of H. takanoi along the northern European coasts from the Bay of Biscay to the Baltic Sea (Makino et al. 2018Makino W., Miura O., Kaiser F., et al. 2018. Evidence of multiple introductions and genetic admixture of the Asian brush-clawed shore crab Hemigrapsus takanoi (Decapoda: Brachyura: Varunidae) along the Northern European coast. Biol. Invasions 20: 825-842. https://doi.org/10.1007/s10530-017-1604-0 ). Hemigrapsus takanoi is an active predator and its large populations are modifying the population dynamics of native mussels (Nour et al. 2020Nour O.M., Stumpp M., Lugo S.C.M., et al. 2020. Population structure of the recent invader Hemigrapsus takanoi and prey size selection on Baltic Sea mussels. Aquat. Invasions 15: 297-317. https://doi.org/10.3391/ai.2020.15.2.06 ). As a euryhaline species, H. takanoi is an excellent osmoregulator, which allows it to colonize variable-salinity environments (Shinji et al. 2009Shinji J., Strüssmann C.A., Wilder M.N., Watanabe S. 2009. Short-term responses of the adults of the common Japanese intertidal crab, Hemigrapsus takanoi (Decapoda: Brachyura: Grapsoidea) at different salinities: osmoregulation, oxygen consumption, and ammonia excretion. J. Crust. Biol. 29: 269-272. https://doi.org/10.1651/08-2998R.1 ). During the larval phase, this crab develops through five zoeal stages and one megalopal stage (Landeira et al. 2019Landeira J.M., Cuesta J.A., Tanaka Y. 2019. Larval development of the brush-clawed shore crab Hemigrapsus takanoi Asakura & Watanabe, 2005 (Decapoda, Brachyura, Varunidae). J. Mar. Biol. Assoc. U.K. 99: 1153-1164. https://doi.org/10.1017/S002531541900002X ). The zoeal stages of H. takanoi are not euryhaline and need to be exported from the brackish water (spawning grounds) into the marine environment for larval development and to recolonize brackish water in the megalopal stage (Mingkid et al. 2006Mingkid W.M., Yokota M., Watanabe S. 2006. Salinity tolerance of larvae in the Penicillate crab Hemigrapsus takanoi (Decapoda: Brachyura: Grapsidae). La mer 44: 17-21. ). The first zoeal stage is especially vulnerable since survival, swimming performance and feeding behaviour are compromised under low-salinity conditions (Landeira et al. 2020Landeira J.M., Liu B., Omura T., et al. 2020. Salinity effects on the first larval stage of the invasive crab Hemigrapsus takanoi: Survival and swimming patterns. Estuar. Coast. Shelf Sci. 245: 106976. https://doi.org/10.1016/j.ecss.2020.106976 ).
The present study aimed to evaluate the morphological variability among broods in Hemigrapsus takanoi zoea I and the effect of the size of the mother on the size and morphology of the progeny. To this end, we performed morphometrics on the carapace of the larvae and their mother. We hypothesized that there is a positive relationship between the size of the mother and the progeny’s morphometrics.
MATERIALS AND METHODS
⌅Sampling and study area
⌅Crabs were collected in Daiba Park (35°38’04”N; 139°46’26”E), located in the inner part of Tokyo Bay, where H. takanoi is the dominant intertidal crab. On 26 May 2017, ovigerous crabs (N = 120) were collected during the low-tide period from the intertidal zone (approx. 3 m2 area) by hand flipping cobbles. In the field, the specimens were identified following the key characters of pigmentation pattern on the abdominal somites and on the cephalothorax described by Asakura and Watanabe (2005)Asakura A., Watanabe S. 2005. Hemigrapsus takanoi new species, a sibling species of the common Japanese intertidal crab H. penicillatus (Decapoda: Brachuyra: Grapsoidea). J. Crust. Biol. 25: 279-292. https://doi.org/10.1651/C-2514 . Only crabs with embryos in an advanced stage of development (eyes visible) were used for the experiment. The crabs were selected to cover the widest size range possible (carapace length ranged from 8 to 15 mm). After collection, the crabs were transported to the aquarium facilities of the Tokyo University of Marine Science and Technology at Shinagawa Campus.
Culture and maintenance
⌅Inside a temperature-controlled room, the crabs were placed individually in 1 L plastic buckets containing 0.8 L of 20°C and 25 salinity seawater (field conditions at time of collection) with aeration and under natural daylight. Every morning, the water was changed, and the crabs were fed pieces of the wakame seaweed, Undaria pinnatifida (Harvey) Suringar 1873. Before that, the buckets were checked to collect newly hatched larvae. Then, both the larvae and the mother were preserved in 80% ethanol for morphometric analysis. To minimize the potential effect of long-term maintenance conditions, we used only 35 broods of larvae from crabs that released the larvae within four days of incubation. The rest of the crabs were discarded and returned to the field.
Morphometric analysis
⌅We applied morphometric measurements of 35 ovigerous crabs and 20 preserved larvae (zoea I) from their broods. Females were placed next to a scale for calibration and photographed. Then, the carapace length was measured on the digitalized images using ImageJ 1.50i software (http://rsb.info.nih.gov/ij/). The larvae were processed using a stereomicroscope with a calibrated eye-piece graticule, and the following measurements were obtained: carapace length (CL) measured from the base of the rostral spine to the posterior dorsal margin of the carapace; length of dorsal spine of carapace (DSL); length of lateral spine of carapace (LSL); rostral spine length (RSL); rostral-dorsal spine distance (RDSD), measured from the tip of the rostral spine to the tip of the dorsal spine; lateral spine distance (LSD), measured from the tip of the lateral spines of the carapace. All larval measurements were completed within one month after specimen preservation.
Statistical analysis
⌅Multivariate statistical analysis was used to identify morphological patterns. In a first step, to classify the individual larvae of all broods, hierarchical clustering analysis was performed using the log-formed morphometric data of each measurement for each larva in a Euclidean distance matrix. After this, two one-way analysis of similarity (ANOSIM) tests were performed to evaluate differences between morphogroups of larvae and between broods (Clarke and Gorley 2015Clarke K.R., Gorley R.N. 2015. PRIMER v7: User Manual/Tutorial. PRIMER-E, Plymouth.). The morphometric measurements contributing most to the differences between morphogroups were identified with the similarity percentage (SIMPER) procedure. Using the same Euclidean distance matrix, a principal component analysis (PCA) was used to examine patterns in larval morphometry among broods and to characterize which morphometric measurements were driving them. To link larval morphometric patterns among broods with the size of the mother, the PCA ordination was represented by superimposing circles of increasing size related to the CL of the female crabs. Moreover, for an easier pattern visualization, cluster morphogroups were also overlaid (Clarke and Gorley 2015Clarke K.R., Gorley R.N. 2015. PRIMER v7: User Manual/Tutorial. PRIMER-E, Plymouth.). To analyse inter-brood differences for each morphometric measurement, the nonparametric Kruskal-Wallis ANOVA test was used, because of the non-normal and heteroscedastic nature of the data set (Kruskal and Wallis 1952Kruskal W.H., Wallis W.A. 1952. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47: 583-621. https://doi.org/10.1080/01621459.1952.10483441 ). Parametric tests were not used because most of the variables did not meet the underlying conditions of normality (Kolmogorov-Smirnov) and homogeneity of variances (Levene’s test). Simple linear regressions were carried out to ascertain the extent of the effect of mother size on the morphometry of zoea I larvae. The statistical analyses were carried out using the PRIMER v7 software and the IBM SPSS Statistics package v27.
RESULTS
⌅The classification of individual larvae by dendrogram, based on morphometric variables, revealed clear grouping patterns. Indeed, two main morphogroups of larvae (Group A and Group B) were categorized by a Euclidean distance of 0.5 (Fig. 1A). The ANOSIM test supported the classification, since these morphogroups were statistically different (R=0.79; P<0.1%). The SIMPER procedure identified that LSL, RSL and DSL were the variables contributing most to the dissimilarities between morphogroups, with 47.2%, 15.0% and 14.6%, respectively. The PCA ordination allowed larval morphometric variability to be reduced to two principal components (Fig. 1B), which explained 85.6% of the cumulative variation (PC1, 75.5% and PC2, 10.1%). Broods showed a patchy distribution (with certain overlapping), suggesting high morphological similarity among the larvae of the same brood. Differences between broods were confirmed by the ANOSIM test (R=0.58; P<0.1%). Eigenvectors of morphometric measurements pointed towards the same direction in the PC1, indicating a similar effect in the ordination pattern (Fig. 1B). The superimposition of the mother carapace size on the PCA plot showed no clear pattern among broods, but the smallest crabs were mainly located on the negative side of the PC2 axis, towards the direction of the LSL eigenvector (Fig. 1C). Correlations between maternal size and each PC axis showed a moderate significant correlation with PC1 (Pearson’s correlation, r=0.457, P<0.001), but no correlation with PC2 (Pearson’s correlation, r=0.127, P<0.01). However, when the morphogroups were coded with colours in the PCA, a clear separation of the larvae was visible, showing larvae from Group A located on the negative side of PC1 and those of Group B on the positive side (Fig. 1C).
Inter-brood variability was visible by observing them independently (Fig. 2) and statistically different (KW-ANOVA, P<0.001) for most of the morphometric measurement (CL, DSL, LSL, RDSL, and RSL), but LSD showed no significant trend (KW-ANOVA, P=0.192). We also found significant inter-brood differences between morphogroups for each morphometric measurement (KW-ANOVA, P<0.001), because larvae of morphogroup A were bigger and had longer spines than those of morphogroup B (Table 1, Supplementary Table 1, 2). When the broods were arranged according to the mother CL, a pattern came out which relates the mother size and the larval morphometry by brood (Fig. 2). Thus, smaller females tended to produce larvae of morphogroup B, whereas larger females tended to produce larvae of morphogroup A.
Morphogroup A | Morphogroup B | |||||
---|---|---|---|---|---|---|
min | max | mean | min | max | mean | |
CL | 0.29 | 0.49 | 0.41 ± 0.04 | 0.29 | 0.50 | 0.39±0.04 |
LSD | 0.42 | 0.69 | 0.54 ± 0.05 | 0.39 | 0.64 | 0.51±0.05 |
RDSD | 0.69 | 1.08 | 0.89 ± 0.07 | 0.67 | 1.09 | 0.86±0.08 |
DSL | 0.23 | 0.41 | 0.31 ± 0.03 | 0.22 | 0.39 | 0.29±0.04 |
RSL | 0.16 | 0.34 | 0.26 ± 0.03 | 0.19 | 0.36 | 0.25±0.03 |
LSL | 0.08 | 0.22 | 0.13 ± 0.03 | 0.06 | 0.19 | 0.11±0.03 |
Correlations | Linear regression model | ANOVA | ||||||
---|---|---|---|---|---|---|---|---|
N | Pearson | P | r2 | intercept | slope | F | p | |
DSL | 35 | 0.475 | 0.002 | 0.226 | 0.188 | 0.01 | 9.64 | 0.004 |
RSL | 35 | 0.551 | <10-4 | 0.302 | 0.148 | 0.01 | 14.311 | 0.001 |
LSL | 35 | 0.468 | 0.002 | 0.219 | 0.039 | 0.007 | 9.247 | 0.005 |
RDSD | 35 | 0.451 | 0.003 | 0.204 | 0.661 | 0.019 | 8.449 | 0.006 |
LSD | 35 | 0.438 | 0.004 | 0.192 | 0.383 | 0.013 | 7.848 | 0.008 |
CL | 35 | 0.879 | <10-4 | 0.772 | 0.314 | 0.011 | 111.638 | <10-4 |
To determine whether maternal body size can predict the size of the larvae, we performed linear regressions (Fig. 3). We found significant positive correlations between the CL of the female crab and each morphometric variable measured in the zoea I larvae (Table 2). This correlation was particularly strong between the CL of the zoea I and the mother (Pearson, r=0.879). This regression model predicted 77% of the variance and showed a good fit for the data (F=111.638, P<10-4).
DISCUSSION
⌅Variations in the larval morphology of decapod crustaceans have been accepted as a common phenomenon, though only a few experiments have tried to explain the mechanisms behind them. Using the first zoeal stage of H. takanoi hatched in the laboratory from ovigerous crabs collected at the same location and time, we found differences in the larval morphology between broods. Interestingly, our results also showed a consistent positive relationship between the size of the larvae and the crab mother.
Our results support the initial hypothesis that bigger mothers produce bigger larvae in the crab H. takanoi. Interestingly, previous results with a sibling species, Hemigrapsus crenulatus, reported that offspring size was not related to female body size (Urzúa et al. 2018Urzúa Á., Bascur M., Guzmán F., Urbina M. 2018. Carry-over effects modulated by salinity during the early ontogeny of the euryhaline crab Hemigrapsus crenulatus from the Southeastern Pacific coast: Development time and carbon and energy content of offspring. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 217: 55-62. https://doi.org/10.1016/j.cbpa.2018.01.001 ). In this case, it is possible that the narrow size range (22-26 mm carapace width) of the female crabs used in that experiment masked the maternal influence on the progeny. In caridean shrimps, this pattern seems more consistent and it has been suggested that larger mothers may provision offspring more efficiently than smaller mothers (Oliphant and Thatje 2021Oliphant A., Thatje S. 2021. Variable shrimp in variable environments: reproductive investment within Palaemon varians. Hydrobiologia 848: 469-484. https://doi.org/10.1007/s10750-020-04455-z ). In fact, González-Ortegón et al. (2018)González-Ortegón E., Vay L. Le, Walton M.E.M.K., Giménez L. 2018. Maternal trophic status and offspring phenotype in a marine invertebrate. Sci. Rep. 8: 1-9. https://doi.org/10.1038/s41598-018-27709-2 demonstrated the importance of diet on offspring provisioning. These authors found that larger female Palaemon serratus feed at higher trophic levels, which was related to the production of bigger eggs with higher content of carbon and nitrogen. In our case, it is not easy to find this diet-mediated relationship, but it is possible that bigger H. takanoi females can benefit from the high densities of crabs occurring in the sampling location at Tokyo Bay. Thus, under these conditions there may occur a density-dependent regulation, in which the big females prey more frequently on juveniles and new recruits, as observed in other cannibalistic crab species (Eggleston and Armstrong 1995Eggleston D.B., Armstrong D.A. 1995. Pre- and post-settlement determinants of estuarine Dungeness crab recruitment. Ecol. Monogr. 65: 193-216. https://doi.org/10.2307/2937137 ).
After hatching as larva, swimming to the upper layers of the water column is a key behaviour that facilitates the offshore exportation into the marine environment (Queiroga and Blanton 2005Queiroga H., Blanton J. 2005. Interactions between behavior and physical forcing in the control of horizontal transport of decapod crustacean larvae. Adv. Mar. Biol. 47: 107-213. https://doi.org/10.1016/S0065-2881(04)47002-3 ), since the zoea stages cannot develop under low-salinity conditions (Mingkid et al. 2006Mingkid W.M., Yokota M., Watanabe S. 2006. Salinity tolerance of larvae in the Penicillate crab Hemigrapsus takanoi (Decapoda: Brachyura: Grapsidae). La mer 44: 17-21. , Landeira et al. 2020Landeira J.M., Liu B., Omura T., et al. 2020. Salinity effects on the first larval stage of the invasive crab Hemigrapsus takanoi: Survival and swimming patterns. Estuar. Coast. Shelf Sci. 245: 106976. https://doi.org/10.1016/j.ecss.2020.106976 ). Few studies deal with the morphology and swimming during the larval stages of decapod larvae. However, Caracappa and Monroe (2019)Caracappa J.C., Munroe D.M. 2019. Variability in swimming behavior among broods of blue crab (Callinectes sapidus) zoeae. J. Exp. Mar. Biol. Ecol. 518: 151176. https://doi.org/10.1016/j.jembe.2019.151176 reported interesting variability in the swimming behaviour among broods of blue crab Callinectes sapidus zoeae. For example, they found that larger zoeae also displayed faster vertical pathways, spent more time swimming upward, and spent more time in motion. Our data cannot support that bigger zoeae of H. takanoi can also perform a better transport offshore than smaller ones. However, swimming has a higher energy demand, and bigger larva can benefit from a better energy allocation due to a maternal influence. In any case, future experiments testing this size-swimming performance relationship would help to understand the maternal influence in this topic.
We also found a positive allometric relationship between the size of the larvae and the length of the carapace spines. It is known that the rostral, dorsal and lateral spines of the carapace are involved in controlling the position and orientation of the zoea while swimming (Smith and Jensen 2015Smith A.E., Jensen G.C. 2015. The role of carapace spines in the swimming behavior of porcelain crab zoeae (Crustacea: Decapoda: Porcellanidae). J. Exp. Mar. Biol. Ecol. 471: 175-179. https://doi.org/10.1016/j.jembe.2015.06.007 ). In relation with this, Caracappa and Monroe (2019)Caracappa J.C., Munroe D.M. 2019. Variability in swimming behavior among broods of blue crab (Callinectes sapidus) zoeae. J. Exp. Mar. Biol. Ecol. 518: 151176. https://doi.org/10.1016/j.jembe.2019.151176 found that zoeae with longer spines were also more likely to swim upward. In our case it is likely that longer spines could help keep zoeae oriented upward, maximizing the propulsive thrust while swimming vertically towards the upper layer of the water column. Moreover, the development of long spines may help to defend the larva from gape-limited planktivorous fishes (Bashevkin et al. 2020Bashevkin S.M., Christy J.H., Morgan S.G. 2020. Adaptive specialization and constraint in morphological defences of planktonic larvae. Funct. Ecol. 34: 217-228. https://doi.org/10.1111/1365-2435.13464 ). Predation pressure seems to be an important driver for morphological plasticity, since fish kairomones may induce spine elongation during the development of crab larvae (Charpentier et al. 2017Charpentier C.L., Wright A.J., Cohen J.H. 2017. Fish kairomones induce spine elongation and reduce predation in marine crab larvae. Ecology 98: 1989-1995. https://doi.org/10.1002/ecy.1899 ).
The present study adds new insights to understanding the morphological variability of crustacean larvae and how these patterns can be linked to maternal size. In the context of reproducibility of animal research, morphological differences between broods should be considered for future experimental designs.