Next Article in Journal
Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds
Next Article in Special Issue
Insights on the Evolutionary History and Genetic Patterns of Octopus vulgaris Cuvier, 1797 in the Northeastern Atlantic Using Mitochondrial DNA
Previous Article in Journal
Spatial Variability of External Egg Quality in Vertical Naturally Ventilated Caged Aviaries
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 13
Issue 4
10.3390/ani13040752
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ontogenetic Variation and Sexual Dimorphism of Beaks among Four Cephalopod Species Based on Geometric Morphometrics

by
Chao Wang
1 and
Zhou Fang
1,2,3,4,5,*
1
College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China
2
National Engineering Research Center for Oceanic Fisheries, Shanghai Ocean University, Shanghai 201306, China
3
Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China
4
Key Laboratory of Oceanic Fisheries Exploration, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China
5
Scientific Observing and Experimental Station of Oceanic Fishery Resources, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China
*
Author to whom correspondence should be addressed.
Animals 2023, 13(4), 752; https://doi.org/10.3390/ani13040752
Submission received: 4 January 2023 / Revised: 10 February 2023 / Accepted: 16 February 2023 / Published: 19 February 2023
(This article belongs to the Special Issue Assessment and Management of Cephalopod Fisheries and Ecosystems)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Octopus minor, Uroteuthis edulis, Sepia esculenta and Sthenoteuthis oualaniensis are important economic cephalopod species in coastal waters of China. They are very important role in marine ecosystems as relevant prey for large marine fish, and marine mammals are a critical component of the food chain. As the main feeding organ of cephalopods, the beak has a stable structure and is resistant to corrosion. However, the sexual dimorphism on the beak of the O. minor, U. edulis and S. esculenta and Sthenoteuthis oualaniensis at different ontogenetic stages is unknown, neither has it been determined whether variation of beak shape relates to changes in habitat environment and feeding preference at different ontogenetic stages. Using a geometric morphometrics approach, we found that the beaks of O. minor, U. edulis, S. esculenta, and Sthenoteuthis oualaniensis showed a pattern of variation, displaying sexual dimorphism and allometry between ontogenetic stages. Habitat may drive variation in beak shape. This study has furthered our understanding of the beak shape ontogenetic variation among the four species. We discuss the potential factors underlying the beak shape variation and provide a basis for the understanding of cephalopod phenotypic plasticity and its ecological significance.

Abstract

Investigating the ontogenetic variation of biological individuals helps us to fully understand the characteristics of evolution. In order to explore the ontogenetic variation and sexual dimorphism of the beak shape in Octopus minor, Uroteuthis edulis, Sepia esculenta and Sthenoteuthis oualaniensis of the China’s coastal waters, the differences between immature and mature stages and the sex-linked differences in the beak shape and size were analyzed with geometric morphometrics methods in this study. The results of Procrustes analysis of variance, principal component analysis and multivariate regression showed that the shapes of the upper beaks of O. minor, U. edulis and S. esculenta differed significantly among various ontogenetic stages (p < 0.05). The shapes of the lower beaks of U. edulis, S. esculenta and Sthenoteuthis oualaniensis were also significantly different among various ontogenetic stages (p < 0.05). The results of thin-plate spline deformation grids showed that the beaks of the four cephalopod species presented different variation patterns. This study gives us basic beak geometry morphology information for Octopus minor, Uroteuthis edulis, Sepia esculenta and Sthenoteuthis oualaniensis present in China’s coastal waters. The ontogenetic differences in beak shape might be related to extrinsic factors (diet difference and intra and interspecific competition) in habitat.

1. Introduction

Biological individuals in nature often face different habitat environments which can drive phenological variation [1] and the evolution of organisms. Shape variations of biological individuals also reflect variation in their habitats [2]. In addition, the phenomenon of sexual dimorphism is also common in biological groups due to natural selection and gender selection. Sexual selection usually affects males through female choice and may also be caused by differences in lifestyle habits [3,4,5,6].
Two methods are generally used to explore ontogenetic and sexual dimorphism. First, based on traditional morphometrics where the linear distance difference among individuals is analyzed by multivariate statistical methods. However, the shape variation of biological individuals involves geometric information, which traditional morphometrics tend to ignore [7]. Secondly, geometric morphometrics (GM) was proposed in the 1990s, which addressed the problems of traditional morphometrics [8,9,10]. This method can quantify the geometric variation of organisms and accurately measure complex biological shapes. Therefore, it is widely used in morphological studies [11,12,13]. GM is the statistical analysis of the form based on Cartesian landmark coordinates [7,10,14]. InGM, by using a superimposition method, for instance, Procrustes, which is used to remove the effect of size, position and orientation of the landmark configurations, exclusively shape variables are obtained. Then, the Procrustes shape coordinates variation and influencing factors are studied by multivariate statistical methods, such as principal component analysis (PCA) and discrimination analysis (DA). GM is an effective method to evaluate shape variations at the intraspecific level as it allows the assessment of ontogenetic changes and sexual dimorphism [15], as has been evidenced for Diplodus puntazzo [16], Munida rugosa [17], Lates niloticus [1] and Hysterocarpus traskii [18]. Assessing ontogenetic changes and sexual dimorphism of organism structures can help distinguish taxa (e.g., species, subspecies and population) and understand their possible ecological role [18,19,20].
Cephalopoda contains around 850 extant species [21], which are widely distributed in the world. Cephalopods play an important role in marine ecosystems as relevant prey for large marine fish, and marine mammals are also a critical component of the food chain [22,23,24,25]. Phenotypic plasticity in response to environment variability is one of the main characteristics of cephalopods [26,27]. Storero et al. [26] showed the variability and plasticity of Octopus tehuelchus in response to the environment. Fang et al. [27] discussed the possible phenotypic plasticity of body and beak of Ommastrephes bartramii in different aspects. Therefore, exploring the life cycle of different cephalopods is prerequisite for sustainable exploitation and utilization of this resource.
As the main feeding organ of cephalopods, the beaks are located in the buccal mass and divided into the upper beak and the lower beak [15,28,29,30]. A beak is a hard tissue structure, which has a stable structure and is resistant to corrosion [28]. Therefore, it is widely used in the study of feeding ecology [31,32] and stock discrimination [28,29,30]. Crespi-abril et al. [33] discovered that the shape of the beaks of Illex argentinus did not vary between cohorts and sexes nor between different ontogenetic phases. Jones et al. [19] confirmed that Patagonian squid (Doryteuthis gahi) populations had complex structures and high intra-species variations in body shape. Sexual dimorphism of the beak is a common phenomenon in some cephalopods [15,34,35]. However, the sexual dimorphism of the beak in some cephalopods at different ontogenetic stages is unknown, and little is known about the inter-population variation in beak shapes of Cephalopods species.
Cephalopods are diverse in China seas [23,36,37,38]. O. minor, U. edulis, S. esculenta, and Sthenoteuthis oualaniensis, representatives of Octopodiformes and Decapodiformes in the subclass Coleoidea, are important economic cephalopod species in coastal waters of China [39,40,41]. This study aims to assess variation in the shape of beaks in different ontogenetic stages of these four squid species and to analyze the sex-linked differences in beak shape and size, relating these differences to relevant factors in the life history of these species. Our analysis will expand the current knowledge about the beak development pattern of cephalopods.

2. Materials and Methods

2.1. Specimens

The samples were collected from 2018 to 2021 by Chinese commercial jigging and trawl vessels in the East China Sea and South China Sea and included 198 O. minor, 205 U. edulis, 198 S. esculenta, and 218 Sthenoteuthis oualaniensis (Table 1). All samples were immediately frozen to −18 °C and shipped back to the laboratory for biological experiments.
Dorsal mantle length (ML) of each specimen was measured by measuring tape and the sex and gonadal maturity of each species were then determined [42,43,44]. In this study, two ontogenetic stages were considered, catalogued as mature and immature according to the maturity stages proposed for each taxon (O. minor and S. esculenta: [42,44]; U. edulis and Sthenoteuthis oualaniensis: [43]). After completing biological experiments as above, the beaks were removed from the buccal mass with tweezers, cleaned, placed in a glass bottle containing 75% ethyl alcohol, and numbered for subsequent analysis.

2.2. Data Acquisition

We focused on lateral profiles of the upper and lower beaks (Figure 1), which were commonly used in previous studies [15,29,45]. The beak photos were taken in a small photography shed. After debugging the small photography shed and placing the scale ruler, the NiKonD750 camera was fixed with a tripod (the parameters of shooting lens are Micro 105 mm f/2.8). To prevent shooting errors, the positions of beaks, camera and focusing point were kept unchanged during shooting. Beak samples were taken out one by one from the glass bottle and positioned, photos were taken immediately. In order to facilitate reading image data with the software and ensure the efficiency of subsequent shape analysis, the images were preprocessed in Photoshop CS6. The obtained images were then used to establish landmarks.
The images of beaks were each named according to species, sex, ontogenetic stage, and the specimen number. We added semilandmarks to better represent beak shape [46]. According to the previous studies [15], 8 landmarks and 12 semilandmarks were marked on the upper beak (Figure 1). Additionally, 10 landmarks and 10 semilandmarks were marked on the lower beak. Landmarks and semilandmarks of the beaks were digitized in 2D with the software TpsDig2 [47]. The coordinate data were read with the “readland.tps” function in the R package “geomorph” [48]. The landmarks and semilandmarks were established with the “define.sliders” function in the R package “geomorph” [48]. The process of establishing landmarks was repeated twice to minimize measurement errors [49]. The two coordinate replicates differed by less than 5%; the data were averaged for the subsequent shape analysis.

2.3. Geometric Morphometrics Analysis

The ontogenetic variation and sexual dimorphism of beak shape of O. minor, U. edulis, S. esculenta, and Sthenoteuthis oualaniensis were analyzed with GM methods.
All morphological measurements and statistical analyses were performed in the R 4.0.5 “geomorph” package [48]. Firstly, the landmarks of all samples were transformed and rotated in order to eliminate the effects of non-shape variation by generalized Procrustes analysis (GPA) with the “gpagen” function [8,48]. The squared sum of the distances from all landmarks to the centroid was defined as the centroid size, which can be used to represent the size of the organism’s structure [50]. Therefore, the centroid size obtained from the landmarks data of beaks was used as an index of the beak size in order to compare the variation in the beak size between immature and mature specimens. We conducted a Procrustes analysis of variance (ANOVA) with the “procD.lm” function to assess ontogenetic differences and sexual dimorphism in the upper and lower beak shape [48,51]. In order to reduce the spatial dimensions of shape data, principal component analysis (PCA) of landmarks data was performed with the “gm.prcomp” function to determine the main components of the shape variation in the upper and lower beak samples [48], and the beak shape variation associated with principal component 1, 2 (PC1, PC2) and principal component 3, 4 (PC3, PC4) was graphically depicted. Then, we also performed the multivariate regression analysis to compare the allometric patterns of beaks in ontogenetic stages based on the relationship between beak shape and centroid size [2,10]. Finally, the shape variation of beaks by ontogenetic stage was visualized with thin-plate spline (TPS) deformation grids of the “plotRefToTarget” function [48,52].

3. Results

3.1. Ontogenetic Variation

There are significant differences in the size of the upper beaks among the four cephalopod species (p < 0.05; Table 2), S. esculenta had the largest beaks, followed by U. edulis and Sthenoteuthis oualaniensis. The beaks of O. minor were the smallest. The upper beaks are slightly larger than the lower beaks in all four species (Figure 2). In addition, during ontogenesis, beak size also increased accordingly. The beak size of O. minor and S. esculenta showed the more significant variation (Figure 2a,c) compared to those of U. edulis and Sthenoteuthis oualaniensis between ontogenetic stages, respectively (Figure 2b,d).
The results of the Procrustes analysis of variance showed that the shapes of the upper beaks of O. minor, U. edulis, and S. esculenta differed significantly between ontogenetic stages (p < 0.05; Table 2), but the shape of the upper beak of Sthenoteuthis oualaniensis showed no significant difference between immature and mature stages (p > 0.05; Table 2). Results of principal component analysis also suggest more overlap between the immature and mature stages of the upper beak of the Sthenoteuthis oualaniensis (Figure 3d). The results of principal component analysis on the shape of the upper beak indicated that the first four principal components (PC1, PC2, PC3, and PC4) altogether accounted for 60.8%, 65.8%, 55.9% and 53.4% of the variation of the shapes of the upper beaks of O. minor, U. edulis, S. esculenta and Sthenoteuthis oualaniensis and the upper beaks between immature and mature stages could be better distinguished in the scatter plot of principal components (Figure 3a–d). The results of the deformation mesh analysis of thin-plate splines deformation grids showed that the main shape variations of the upper beaks of O. minor, U. edulis and S. esculenta between immature and mature stages occurred in the rostrum (Landmarks 1, 2, and 3), hood (Landmarks 16, 18, and 19), wing (Landmarks 4 and 20) and lateral wall (Landmarks 7, 10, 11, and 13) (Figure 4). Although the morphological variations in upper beaks of each species occurred at similar sections between immature and mature stages, four cephalopod species still had different shape variation directions in various parts between immature and mature stages. The rostra of the O. minor and S. esculenta gradually became shorter, more blunt and thicker (Figure 4a,c), but U. edulis became longer, sharper and thinner (Figure 4b). The hoods of O. minor and U. edulis became more curved and wider (Figure 4a,b), whereas the hoods of S. esculenta became curved and compressed (Figure 4c). The lateral wall of the upper beak of each species became wider in different directions. In addition, differences were apparent in changes of upper beak shape with increasing upper beak size (Figure 5a–d) and the predicted shape of the upper beak was positively correlated with the upper beak size (Figure 5a–d).
The results of the Procrustes analysis of variance showed that the size of the lower beak of four species showed significant variation (p < 0.05; Table 3). The shapes of the lower beaks of U. edulis, S. esculenta, and Sthenoteuthis oualaniensis were significantly different between immature and mature stages (p < 0.05; Table 3), but the shape of the lower beaks of O. minor showed no significant variation (p > 0.05; Table 3). The first four principal components (PC1, PC2, PC3, and PC4) altogether accounted for 57.2%, 49.9%, 54.5% and 58.0% of the variation in lower beak shape of the lower beaks of O. minor, U. edulis, S. esculenta, and Sthenoteuthis oualaniensis (Figure 6a–d). The explanatory rate of the first 4 principal components for the variation of the lower beak was slightly lower than that for the variation of the upper beak. The scatter plots of principal components for the beak shape in immature stages partially coincided with those in mature stages, but they could still explain the variation in the beak shape among various ontogenetic stages (Figure 6b–d). The results of the deformation mesh analysis of thin-plate splines deformation grids showed that the main variations of the shapes of the lower beaks of U. edulis, S. esculenta, and Sthenoteuthis oualaniensis in the immature and mature stages occurred at rostrum (Landmarks 1 and 18), hood (Landmark 17), wing (Landmarks 4, 5, 6, and 7) and lateral wall (Landmarks 9, 11, and 14) (Figure 7). In mature stages, U. edulis had sharper rostrum, wider wing, longer hood, and wider lateral wall (Figure 7b) and similar variations were observed in the lower beaks of S. esculenta and Sthenoteuthis oualaniensis (Figure 7c,d). In particular, the shapes of the lower beaks of different species showed different variations. For instance, the hoods of U. edulis and S. esculenta became more curved (Figure 7b,c), whereas the hood of Sthenoteuthis oualaniensis flattened (Figure 7d). Moreover, the variation of the shape of the lower beak with size was not significant (Figure 8b–d). Compared with the upper beak, the lower beak showed significant difference in variation rate of shape among four species in both immature and mature stages. The beak shape variation pattern of S. esculenta was different from other species. Sthenoteuthis oualaniensis had similar variation pattern of the lower beak (Figure 8d) and the predicted shape of the lower beak was positively correlated with beak size. The shape of the lower beak of S. esculenta showed no significant variation with the increase of beak size in immature stages, but the shape of the lower beak of S. esculenta was negatively correlated with beak size in mature stages (Figure 8c, right).

3.2. Sexual Dimorphism

Firstly, the results of the Procrustes analysis of variance indicated that the shapes of the upper beaks of O. minor, U. edulis and Sthenoteuthis oualaniensis were significantly different between female and male (p < 0.05; Table 2), whereas the shapes of the upper beaks of S. esculenta showed no sexual dimorphism (p > 0.05; Table 2). The scatter plot of principal components indicated that the spatial proportions of the upper beaks of female and male individuals showed some differences (Figure 3a,b,d). The upper beak of male individuals of O. minor and U. edulis had a more protruding lower recess of the lateral wall, wider hood and sharper rostra than did female individuals in immature and mature stages (Figure 4a,d). Unlike O. minor and U. edulis, the upper beaks of female individuals of Sthenoteuthis oualaniensis had sharper and longer rostra than did male individuals (Figure 4b). The multiple regression analysis results showed that the upper beaks of O. minor and Sthenoteuthis oualaniensis were not significantly different between male and female individuals (Figure 5a,d, left), but the upper beak of U. edulis was significantly different between male and female individuals (Figure 5b, left). The shapes of the upper beaks of female and male individuals in different species showed the same variation of shape with increasing size (Figure 5, right).
The results of the Procrustes analysis of variance indicated that the S. esculenta showed significant sexual dimorphism in the shape of the lower beak (p < 0.05; Table 3). The shapes of the lower beaks of O. minor, U. edulis and Sthenoteuthis oualaniensis were not significantly different between females and males (p > 0.05; Table 3). The scatter plot of principal components of S. esculenta indicated that the spatial proportions of the principal components of the lower beak shapes were significantly different between female and male individuals (Figure 6c). In immature stages, female individuals of S. esculenta had sharper rostra and more curved hoods of the lower beaks, wider wings, and slightly compressed lateral wall than male individuals. In mature stages, male individuals of S. esculenta had wider hoods, wings and lateral wall of the lower beak than did female individuals (Figure 7c). Multiple regression analysis showed that the difference in the lower beak between male and female individuals was not significant (Figure 8c, left). The shape of the lower beaks of female and male individuals showed the same variation of shape with increasing size (Figure 8c, right).

4. Discussion

4.1. Ontogenetic Pattern of Beak Shape

Cephalopods are short-lived species, so their growth and development are closely related to the habitat environment and their feeding habits are different among various ontogenetic stages [53,54,55]. As a feeding organ of cephalopods, the beak is one of the hardest parts of the cephalopod body and mainly used to bite and fix prey [56,57]. In addition, genetic factors also contribute to the shape differentiation of beaks among different species and different ontogenetic stages [58,59]. This phenomenon is possibly the consequence of the adaptive evolution of different species to habitat environment changes. Therefore, the beak shape of different species during ontogeny may be affected by many factors and display different variation patterns.
In this study, the upper and lower beaks of each cephalopod species showed significant changes in the hood, wing, and rostrum between immature and mature individuals (Figure 4 and Figure 7). The beaks in the cephalopod buccal mass are controlled by musculature in the buccal mass, including the anterior mandibular muscle, lateral mandibular muscle, posterior mandibular muscle and superior mandibular muscle [60,61,62]. The parts of the upper and lower beaks are surrounded by these muscles, include the hood, lateral wall, crest and wing [60,61,62,63]. Since the diets of each cephalopod species are gradually expanded to bigger and harder foods as they grow [41,55,64,65,66], larger beaks are required to produce a stronger bite. Therefore, we assumed that the increase in the proportion of the hood, lateral wall, crest and wing of upper and lower beaks relate to muscle growth, insertion and bite cycle and help to improve the bite force of the beaks. However, the variation pattern of different kinds of beaks were slightly different. For instance, the hoods of the lower beaks of U. edulis and S. esculenta became more curved with growth (Figure 7b,c) and the hoods of the lower beaks of Sthenoteuthis oualaniensis became flat (Figure 7d). The habitats of U. edulis and S. esculenta are more similar (Table 1). Therefore, the different development patterns suggested differences in the function and growth mechanisms of beaks among different cephalopods. These differences may result from adaptation of each species to habitat and feeding preference during ontogeny.
Larger feeding organs and larger beaks specifically would allow organisms to eat larger and harder food [41,55,65,67]. In this study, O. minor inhabits the benthic environment of temperate waters and mainly preys on crustaceans [54,64]). S. esculenta is also a demersal shallow-water species that preys on microcrustaceans in its alevin stage, and bigger crustaceans in its adult stage [53,65]. Due to the differences in feeding types and the complexity of the benthic habitat at different ontogenetic stages, the rostra of O. minor and S. esculenta may be subjected to greater abrasion during feeding and gradually become dull in the maturation stage (Figure 4a,c and Figure 7c). After sexual maturity, U. edulis in the East China Sea mainly preys on juvenile fish of Scombridae and other fishes [55], which are softer than the foods preyed by O. minor and S. esculenta. Sharper upper beaks can help to quickly secure prey and facilitate eating prey [57,68]. As a result, the rostra of O. minor and S. esculenta became dull and the rostrum of U. edulis became sharp. Sthenoteuthis oualaniensis mainly preys on microcrustacean and crustaceans in paralarval and adult stages, respectively [66]. These mature individuals need to eat more food for more energy to sustain life activities such as growth and reproduction [34]. In addition, usually cephalopod beak growth rate of the immature stage is faster. When the carcass grows to a certain stage, the growth rate of the beaks slows [33,69,70]. Therefore, the results in this study suggested that feeding strategies in different ontogenetic stages might be responsible for the difference of beak shape.
In addition to the differences in the internal formation mechanism of beaks and feeding habits, the habitat environment may also affect the variation of beaks. Firstly, the habitat environments of different cephalopod species are different. O. minor are benthic cephalopod species [53,64,65], whereas U. edulis and S. esculenta inhabit the demersal, and Sthenoteuthis oualaniensis inhabits the pelagic [71,72,73]. Secondly, in the growth process of a single species, due to migration and other living habits, the habitat environment is different in various ontogenetic stages. The East China Sea and South China Sea are greatly affected by the continental coastal currents, black tides and southwest and northeast monsoons, and the climate is complex and variable. The four species had a wide swimming range and variable habitat environment in the growth process [55]. Roscian et al.’s [68] study highlights the importance of habitat as a driver of variation in beak shape. The short, blunt, thick rostra of O. minor and S. esculenta are likely more useful for breaking shells (Figure 4a,c), while the long, sharp and thin rostra of U. edulis and Sthenoteuthis oualaniensis are likely more efficient at piercing and tearing fish (Figure 7a,c). Long and extremely mobile arms of benthic cephalopods allow them to capture, maintain and manipulate prey without needing to kill it rapidly. In contrast, pelagic cephalopod hunt by projecting the tentacles and rapidly bringing back the prey toward the beak to kill it [54,57,64,68]. Therefore, beak shape changes between immature and mature may reflect an adaptation of cephalopod species to the habitats and diet resources.

4.2. Sexual Dimorphism of Beak Shape

Gender selection, ecological differences, etc. may cause the differences in size and shape between female and male individuals and the shape differences of biological individuals between female and male individuals are called sexual dimorphism [17,18]. Some studies confirmed the presence of sexual dimorphism in cephalopod beaks [15,35]. In this study, the upper beaks of female individuals of each cephalopod species were larger than those of male individuals in the mature stages (p < 0.05) (Figure 2; Table 2 and Table 3). Female squids always grow faster than males during ontogeny, which ultimately leads to the subtle sexual dimorphism [34]. We also suggested that the shapes of the upper beaks of O. minor and U. edulis and the shape of the lower beak of S. esculenta showed gender differentiation in various ontogenetic stages (Table 2). The upper beak of O. minor and Sthenoteuthis oualaniensis females had a more protruding crest and wider hood than that of males (Figure 4a,d). In the mature stages, the hood and wing of the beaks of females of S. esculenta were wider than those of males (Figure 4d and Figure 7d). Similarly, the beak shapes of neon flying squid (Ommastrephes bartramii) also showed significant difference between females and males [15]. The female beaks, with larger, wider hood and sharp rostra, also might be related to the high quality of food requirements after sexual maturation since larger beaks provide stronger feeding ability to prey on high-level food. This also may reflect that beak shape changes in different sexes are an adaptation of cephalopod species to the habitats and diet resources.

4.3. Variations in the Size and Shape of Beaks

Allometry has long been an important focus of the study on evolution and growth and GM is a flexible and powerful tool for studying shape structure evolution and development [12,74]. The multiple regression analysis in this study revealed some differences in allometric variation of beaks of each species (Figure 5 and Figure 8). The intensity of feeding in some cephalopods varied significantly with the growth stage [55,64,66]. The intensity of feeding of individuals gradually increased in the immature stages and decreased after maturity was reached; some individuals did not eat at all after spawning [75]. In this study, the upper beaks of O. minor, U. edulis, and S. esculenta show faster variation rates of shape in their immature stages (Figure 5). It was inferred that the different feeding demands in various ontogenetic stages affected the growth rate and shape on their upper beaks, so beaks grow faster in immature stages. In addition, the rostrum is the main part for cutting foods [56,57]. Therefore, the variation of rostrum of each species was bigger, as indicated in thin-plate spline deformation grids analysis (Figure 4 and Figure 7). Hence, the phenomenon of allometry in the beak shape occurred. In this study, we found that the upper beaks had a higher variation rate of shape than the lower beaks (Figure 5b,c and Figure 8b,c). It was shown that the upper beak grew faster compared to the lower beak, which could also explain why the upper beaks of the four species were slightly larger than the lower beaks (Figure 2). In some cephalopod species, the upper beak grew faster than the lower beak [69,76]. Similar differences in other cephalopods were ascribed to the earlier formation of the upper beak during embryonic development [77]. This study suggests that this difference of size on the upper and lower beaks is also related to differences in allometry.

5. Conclusions

Through geometric morphometrics analysis, we found that the beaks of O. minor, U. edulis, S. esculenta, and Sthenoteuthis oualaniensis showed different variation patterns, displaying sexual dimorphism and allometry between ontogenetic stages. Habitat difference (diet difference and intra and interspecific competition) may drive ontogenetic variation in beak shape. This study has furthered our understanding of the beak development pattern among four cephalopods species. The detection of additional sexual differences and their expression during ontogeny in other cephalopod species is also highly interesting and necessary, considering the strongly specialized life history in cephalopods, and different habitats.

Author Contributions

Conceptualization, C.W.; Methodology, Z.F.; Software, C.W.; Formal analysis, C.W.; Investigation, Z.F.; Writing—original draft, C.W.; Writing—review & editing, Z.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2019YFD0901404), National Nature Foundation of China on the project (NSFC41876141), and Funding for the Opening of Key Laboratories for Offshore Fishery Development by the Ministry of Agriculture (LOF 2021-01).

Institutional Review Board Statement

No ethical approval was required. Beaks were only obtained from dead individuals caught. No live animals were caught specifically for this project.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The support from the captain and crews in commercial jigging vessels in scientific surveys is gratefully acknowledged. This study was financially supported by the National Key R&D Program of China (2019YFD0901404), National Nature Foundation of China on the project (NSFC41876141), and Funding for the Opening of Key Laboratories for Offshore Fishery Development by the Ministry of Agriculture (LOF 2021-01).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Nyboer, E.A.; Chapman, L.J. Ontogenetic shifts in phenotype–environment associations in Nile perch, Lates niloticus (Perciformes: Latidae) from Lake Nabugabo, Uganda. Biol. J. Linn. Soc. 2013, 110, 449–465. [Google Scholar] [CrossRef] [Green Version]
  2. Adams, D.C.; Nistri, A. Ontogenetic convergence and evolution of foot morphology in European cave salamanders (Family: Plethodontidae). BMC Evol. Biol. 2010, 10, 216. [Google Scholar] [CrossRef] [Green Version]
  3. Marian, J.E.A.R.; Apostolico, L.H.; Chiao, C.C.; Hanlon, R.T.; Hirohashi, N.; Iwata, Y.; Mather, J.; Sato, N.; Shaw, P.W. Male Alternative Reproductive Tactics and Associated Evolution of Anatomical Characteristics in Loliginid Squid. Front. Physiol. 2019, 10, 1281. [Google Scholar] [CrossRef] [PubMed]
  4. Apostolico, L.H.; Marian, J.E. Cephalopod mating systems as models for the study of sexual selection and alternative reproductive tactics: A review. Bull. Mar. Sci. 2020, 96, 375–398. [Google Scholar] [CrossRef]
  5. Guo, H.; Zhang, D.; Näslund, J.; Wang, L.; Zhang, X. Effects of spawning group sex ratio and stocking density on the outcome of captive reproduction in golden cuttlefish Sepia esculenta. Aquaculture 2022, 559, 738416. [Google Scholar] [CrossRef]
  6. Guo, H.; Zhang, D.; Wang, L.; Li, W.; He, P.; Näslund, J.; Zhang, X. Sperm competition in golden cuttlefish Sepia esculenta: The impact of mating order and male size. Aquaculture 2020, 530, 735929. [Google Scholar] [CrossRef]
  7. Rohlf, F.J.; Marcus, L.F. A revolution in morphometrics. Trends Ecol. Evol. 1993, 8, 129–132. [Google Scholar] [CrossRef] [PubMed]
  8. Rohlf, F.J. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst. Biol. 1990, 39, 40–59. [Google Scholar] [CrossRef] [Green Version]
  9. Bookstein, F.L. Morphometric Tools for Landmark Data: Geometry and Biology; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  10. Adams, D.C.; Rohlf, F.J.; Slice, D.E. A field comes of age: Geometric morphometrics in the 21st century. Hystrix Ital. J. Mammal. 2013, 24, 7–14. [Google Scholar] [CrossRef]
  11. Rohlf, F.J. On applications of geometric morphometrics to studies of ontogeny and phylogeny. Syst. Biol. 1998, 47, 147–158. [Google Scholar] [CrossRef] [Green Version]
  12. Klingenberg, C.P. Evolution and development of shape: Integrating quantitative approaches. Nat. Rev. Genet. 2010, 11, 623–635. [Google Scholar] [CrossRef] [PubMed]
  13. Echeverry, A.M.; Londoño-Cruz, E.; Benítez, H.A. Quantifying the Geometric Shell Shape between Populations of True Limpets Lottia Mesoleuca (Mollusca: Lottidae) in Colombia. Animals 2020, 10, 675. [Google Scholar] [CrossRef]
  14. Slice, D.E. Landmark coordinates aligned by procrustes analysis do not lie in kendall’s shape space. Syst. Biol. 2001, 50, 141–149. [Google Scholar] [CrossRef] [PubMed]
  15. Fang, Z.; Chen, X.; Su, H.; Thompson, K.; Chen, Y. Evaluation of stock variation and sexual dimorphism of beak shape of neon flying squid, Ommastrephes bartramii, based on geometric morphometrics. Hydrobiologia 2017, 784, 367–380. [Google Scholar] [CrossRef]
  16. Kouttouki, S.; Georgakopoulou, E.; Kaspiris, P.; Divanach, P.; Koumoundouros, G. Shape ontogeny and variation in the sharpsnout seabream, Diplodus puntazzo (Cetti 1777). Aquac. Res. 2006, 37, 655–663. [Google Scholar] [CrossRef]
  17. Claverie, T.; Smith, I.P. Allometry and sexual dimorphism in the chela shape in the squat lobster Munida rugosa. Aquat. Biol. 2010, 8, 179–187. [Google Scholar] [CrossRef]
  18. Parvis, E.S.; Coleman, R.M. Sexual Dimorphism and Size–Related Changes in Body Shape in Tule Perch (Family: Embiotocidae), a Native California Live–Bearing Fish. Copeia 2020, 108, 12–18. [Google Scholar] [CrossRef]
  19. Jones, J.B.; Pierce, G.J.; Saborido-Rey, F.; Brickle, P.; Kuepper, F.C.; Shcherbich, Z.N.; Arkhipkin, A.I. Size-dependent change in body shape and its possible ecological role in the Patagonian squid (Doryteuthis gahi) in the Southwest Atlantic. Mar. Biol. 2019, 166, 54. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, B.L.; Chen, X.J.; Wang, X.H.; Du, F.Y.; Fang, Z.; Xu, L.L. Geographic, intraspecific and sexual variation in beak morphology of purpleback flying squid (Sthenoteuthis oualaniensis) throughout its distribution range. Mar. Freshw. Res. 2019, 70, 417–425. [Google Scholar] [CrossRef]
  21. Hoving, H.J.T.; Perez, J.A.A.; Bolstad, K.S.; Braid, H.E.; Evans, A.B.; Fuchs, D.; Judkins, H.; Kelly, J.T.; Marian, J.E.; Nakajima, R.; et al. The Study of Deep-Sea Cephalopods. Adv. Cephalop. Sci. Biol. Ecol. Cultiv. Fish. 2014, 67, 235–359. [Google Scholar] [CrossRef]
  22. Clarke, M.R.; Maddock, L. The Mollusca Vol. 12. Paleontology and Neontology of Cephalopods; Academic Press: Cambridge, MA, USA, 1988; pp. 123–131. [Google Scholar]
  23. Dong, Z. Fauna of China Mollusca Cephalopoda; Science Press: Beijing, China, 1988. [Google Scholar]
  24. Arkhipkin, A.I.; Rodhouse, P.G.K.; Pierce, G.J.; Sauer, W.; Sakai, M.; Allcock, J.L.; Arguelles, J.R.; Bower, G.; Castillo, L.; Ceriola, C.S.; et al. World Squid Fisheries. Rev. Fish. Sci. Aquac. 2015, 23, 92–252. [Google Scholar] [CrossRef] [Green Version]
  25. Golikov, A.V.; Ceia, F.R.; Hoving, H.J.T.; Queirós, J.P.; Sabirov, R.M.; Blicher, M.E.; Larionova, A.M.; Walkusz, W.; Zakharov, D.V.; Xavier, J.C. Life History of the Arctic Squid Gonatus fabricii (Cephalopoda: Oegopsida) Reconstructed by Analysis of Individual Ontogenetic Stable Isotopic Trajectories. Animals 2022, 12, 3548. [Google Scholar] [CrossRef] [PubMed]
  26. Storero, L.P.; Ocampo-Reinaldo, M.; González, R.A.; Narvarte, M.A. Growth and life span of the small octopus Octopus tehuelchus in San Matias Gulf (Patagonia): Three decades of study. Mar. Biol. 2010, 157, 555–564. [Google Scholar] [CrossRef]
  27. Fang, Z.; Han, P.; Wang, Y.; Chen, Y.; Chen, X. Interannual variability of body size and beak morphology of the squid Ommastrephes bartramii in the North Pacific Ocean in the context of climate change. Hydrobiologia 2021, 848, 1295–1309. [Google Scholar] [CrossRef]
  28. Clarke, M.R. The identification of cephalopod “beaks” and the relationship between beak size and total body weight. Bull. Br. Mus. Nat. Hist. Zool. 1962, 8, 419–480. [Google Scholar] [CrossRef]
  29. Fang, Z.; Fan, J.; Chen, X.; Chen, Y. Beak identification of four dominant octopus species in the East China Sea based on traditional measurements and geometric morphometrics. Fish. Sci. 2018, 84, 975–985. [Google Scholar] [CrossRef]
  30. Pacheco-Ovando, R.; Granados-Amores, J.; González-Salinas, B.; Ruiz-Villegas, J.M.; Díaz-Santana-Iturrios, M. Beak shape analysis and its potential to recognize three loliginid squid species found in the northeastern Pacific. Mar. Biodivers. 2021, 51, 82. [Google Scholar] [CrossRef]
  31. Hobson, K.; Cherel, Y. Isotopic reconstruction of marine food webs using cephalopod beaks: New insight from captively raised Sepia officinalis. Can. J. Zool. 2006, 84, 766–770. [Google Scholar] [CrossRef] [Green Version]
  32. Fang, Z.; Thompson, K.; Jin, Y.; Chen, X.J.; Chen, Y. Preliminary analysis of beak stable isotopes (δ13C and δ15N) stock variation of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean. Fish. Res. 2016, 177, 153–163. [Google Scholar] [CrossRef] [Green Version]
  33. Crespi-Abril, A.C.; Morsan, E.M.; Barón, P.J. Analysis of the ontogenetic variation in body and beak shape of the Illex argentinus inner shelf spawning groups by geometric morphometrics. J. Mar. Biol. Assoc. U. K. 2010, 90, 547–553. [Google Scholar] [CrossRef]
  34. Brunetti, N.E.; Ivanovic, M.L.; Aubone, A.; Pascual, L.N. Reproductive biology of red squid (Ommastrephes bartramii) in the Southwest Atlantic. Rev. De Investig. Ny Desarro. Pesq. 2006, 18, 5–19. [Google Scholar]
  35. Bolstad, K.S. Sexual dimorphism in the beaks of Moroteuthisingens Smith, 1881 (Cephalopoda: Oegopsida: Onychoteuthidae). N. Z. J. Zool. 2006, 33, 317–327. [Google Scholar] [CrossRef]
  36. Norman, M.D.; Nabhitabhata, J.; Lu, C.C. An updated checklist of the cephalopods of the South China Sea. Raffles Bull. Zool. 2016, 34, 566–592. [Google Scholar]
  37. Boavida-Portugal, J.; Guilhaumon, F.; Rosa, R.; Araújo, M.B. Global patterns of coastal cephalopod diversity under climate change. Front. Mar. Sci. 2022, 8, 740781. [Google Scholar] [CrossRef]
  38. Xu, R.; Lü, Y.; Tang, Y.; Chen, Z.; Xu, C.; Zhang, X.; Zheng, X. DNA Barcoding Reveals High Hidden Species Diversity of Chinese Waters in the Cephalopoda. Front. Mar. Sci. 2022, 9, 830381. [Google Scholar] [CrossRef]
  39. Jin, Y.; Chen, X.J. Review on studies of Cephalopoda in the coastal waters of China. Mar. Fish. 2017, 39, 696–712. [Google Scholar] [CrossRef]
  40. Pang, Y.M.; Tian, Y.J.; Fu, C.H.; Wang, B.; Li, J.C.; Ren, Y.P.; Wan, R. Variability of coastal cephalopods in overexploited China Seas under climate change with implications on fisheries management. Fish. Res. 2018, 208, 22–33. [Google Scholar] [CrossRef]
  41. Zhao, C.X.; Shen, C.Y.; Bakun, A.; Yan, Y.R.; Kang, B. Purpleback flying squid Sthenoteuthis oualaniensis in the South China Sea: Growth, resources and association with the environment. Water 2020, 13, 65. [Google Scholar] [CrossRef]
  42. Arkhipkin, A.I. Reproductive System Structure, Development and Function in Cephalopods with a New General Scale for Maturity Stages. J. Northwest Atl. Fish. Sci. 1992, 12, 63–74. [Google Scholar] [CrossRef]
  43. Lipinski, M.R.; Underhill, L.G. Sexual maturation in squid: Quantum or continuum? S. Afr. J. Mar. Sci. 1995, 15, 207–223. [Google Scholar] [CrossRef] [Green Version]
  44. ICES. Report of the Workshop on Sexual Maturity Staging of Cephalopods; ICES: Livorno, Italy, 2010; pp. 1–97. [Google Scholar]
  45. Neige, P.; Dommergues, J.L. Disparity of beaks and statoliths of some coleoids: A morphometric approach to depict shape differentiation. Abh. Der Geol. Bundesanst. 2002, 57, 393–399. [Google Scholar]
  46. Gunz, P.; Mitteroecker, P. Semilandmarks: A method for quantifying curves and surfaces. Hystrix Ital. J. Mammal. 2013, 24, 103–109. [Google Scholar] [CrossRef]
  47. Rohlf, F.J. The tps series of software. Hystrix Ital. J. Mammal. 2015, 26, 9–12. [Google Scholar] [CrossRef]
  48. Adams, D.C.; Otárola-Castillo, E. Geomorph: An r package for the collection and analysis of geometric morphometric shape data. Methods Ecol. Evol. 2013, 4, 393–399. [Google Scholar] [CrossRef]
  49. Viscosi, V.; Cardini, A. Leaf morphology, taxonomy and geometric morphometrics: A simplified protocol for beginners. PLoS One 2011, 6, e25630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Chen, X.J.; Fang, Z.; Chen, Y.; Su, H. Application of Geometric Morphometrics in Aquatic Animals; Science Press: Beijing, China, 2017. [Google Scholar]
  51. Adams, D.C.; Collyer, M.L. Multivariate Phylogenetic Comparative Methods: Evaluations, Comparisons, and Recommendations. Syst. Biol. 2018, 67, 14–31. [Google Scholar] [CrossRef]
  52. Claude, J. Morphometrics with R.; Springer: New York, NY, USA, 2008. [Google Scholar] [CrossRef]
  53. Chen, X.J.; Liu, B.L.; Wang, Y.G. Cephalopods of the World; Ocean Press: Beijing, China, 2009. [Google Scholar]
  54. Jonas, W.N.; Chen, W.; Wang, C.L.; Zhang, X.M.; Li, L.G.; Xu, J.C. Habit and Adaptability of Common Octopus variabilis to Certain Environmental Factors. Fish. Sci. 2010, 29, 74–78. [Google Scholar] [CrossRef]
  55. Li, N.; Fang, Z.; Chen, X.J. Fishery of swordtip squid Uroteuthis edulis: A review. J. Dalian Ocean. Univ. 2020, 35, 637–644. [Google Scholar] [CrossRef]
  56. Hanlon, R.T.; Messenger, J.B. Cephalopod Behaviour; Cambridge University Press: Cambridge, UK, 2018; pp. 74–96. [Google Scholar]
  57. Gao, X.D.; Fang, Z.; Chen, X.J.; Li, Y.K. Variation in beak morphology of Dosidicus gigas in different waters of the eastern equatorial Pacific. J/OL. Fish. China 2021. Available online: https://kns.cnki.net/kcms/detail/31.1283.s.20210528.1454.002.html. (accessed on 28 May 2021).
  58. Tanabe, K. The jaw apparatuses of Cretaceous desmoceratid ammonites. Palaeontology 1983, 26, 677–686. [Google Scholar]
  59. Kroger, B.; Vinther, J.; Fuchs, D. Cephalopod origin and evolution: A congruent picture emerging from fossils, development and molecules. Bioessays 2011, 33, 602–613. [Google Scholar] [CrossRef] [PubMed]
  60. Kear, A.J. Morphology and function of the mandibular muscles in some coleoid cephalopods. J. Mar. Biol. Assoc. U. K. 1994, 74, 801–822. [Google Scholar] [CrossRef]
  61. Uyeno, T.A.; Kier, W.M. Functional morphology of the cephalopod buccal mass: A novel joint type. J. Morphol. 2005, 264, 211–222. [Google Scholar] [CrossRef]
  62. Uyeno, T.A.; Kier, W.M. Electromyography of the buccal musculature of octopus (Octopus bimaculoides): A test of the function of the muscle articulation in support and movement. J. Exp. Biol. 2007, 210, 118–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Miserez, A.; Schneberk, T.; Sun, C.; Zok, F.W.; Waite, J.H. The transition from stiff to compliant materials in squid beaks. Science 2008, 319, 1816–1819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Qian, Y.S. Studies on Ecological Habit and Artificial Reproductive Technique of Octopus minor; Ocean University of China: Qingdao, China, 2011. [Google Scholar]
  65. Wang, L.; Zhang, X.M.; Ding, P.W.; Liu, T.Y.; Chen, S.Q. Reproductive behavior and mating strategy of Sepia esculenta. Acta Ecol. Sin. 2017, 37, 1871–1880. [Google Scholar]
  66. Xie, J.Y.; Zhang, L.Z.; Wu, W.X.; Zhou, B.H.; Chen, Q.J.; Zhao, C.X.; He, X.B.; Xu, J.; Yan, Y.R. Feeding habit and trophic niche of purpleback flying squid (Sthenoteuthis oualaniensis) in the Nansha Islands area, South China Sea. J. Fish. China 2021, 45, 1993–2002. [Google Scholar]
  67. Castilla, A.C.; Hernández-Urcera, J.; Gouraguine, A.; Guerra, A.; Cabanellas-Reboredo, M. Predation behaviour of the European squid Loligo vulgaris. J. Ethol. 2020, 38, 311–322. [Google Scholar] [CrossRef]
  68. Roscian, M.; Herrel, A.; Zaharias, P.; Cornette, R.; Fernandez, V.; Kruta, I.; Cherel, Y.; Rouget, I. Every hooked beak is maintained by a prey: Ecological signal in cephalopod beak shape. Funct. Ecol. 2022, 36, 2015–2028. [Google Scholar] [CrossRef]
  69. Fang, Z.; Xu, L.; Chen, X.; Liu, B.; Li, J.; Chen, Y. Beak growth pattern of purpleback flying squid Sthenoteuthis oualaniensis in the eastern tropical Pacific equatorial waters. Fish. Sci. 2015, 81, 443–452. [Google Scholar] [CrossRef]
  70. Jin, Y.; Liu, B.L.; Chen, X.J.; Staples, K. Morphological beak differences of loliginid squid, Uroteuthis chinensis and Uroteuthis edulis, in the northern South China Sea. J. Oceanol. Limnol. 2018, 36, 559–571. [Google Scholar] [CrossRef]
  71. Jereb, P.; Roper, C.F.E. Cephalopods of the World. An Annotated and Illustrated Catalogue of Cephalopod Species Known to Date. Volume 1. Chambered Nautiluses and Sepioids (Nautilidae, Sepiidae, Sepiolidae, Sepiadariidae, Idiosepiidae and Spirulidae). Vol. 1. FAO Species Catalogue for Fishery Purposes. 4; FAO: Rome, Italy, 2005. [Google Scholar]
  72. Jereb, P.; Roper, C.F.E. Cephalopods of the World. An Annotated and Illustrated Catalogue of Cephalopod Species Known to Date. Volume 2. Myopsid and Oegopsid Squids. Vol. 2. FAO Species Catalogue for Fishery Purposes. 4; FAO: Rome, Italy, 2005. [Google Scholar]
  73. Iglesias, J.; Fuentes, L.; Villanueva, R. Cephalopod. Culture; Springer: New York, NY, USA, 2014; pp. 415–426. [Google Scholar] [CrossRef] [Green Version]
  74. Klingenberg, C.P. Size, shape, and form: Concepts of allometry in geometric morphometrics. Dev. Genes Evol. 2016, 226, 113–137. [Google Scholar] [CrossRef] [PubMed]
  75. Liu, M.N. Study on Morphology and Feeding Ecology of Chinese Squid (Uroteuthis chinensis); Shanghai Ocean University: Shanghai, China, 2020. [Google Scholar] [CrossRef]
  76. Ikica, Z.; Vuković, V.; Durović, M.; Joksimović, A.; Krstulovićšifner, S. Analysis of beak morphometry of the horned octopus Eledone cirrhosa, Lamarck 1798 (Cephalopoda: Octopoda), in the south–eastern Adriatic Sea. Acta Adriat. 2014, 55, 43–56. [Google Scholar]
  77. Von Boletzky, S. Origin of the lower jaw in cephalopods: A biting tissue. Palaeontol. Z. 2007, 81, 328–333. [Google Scholar] [CrossRef]
Animals 13 00752 g001 550
Figure 1. Shape and digitized landmarks (hollow dots) and semi-landmarks (solid dots) of the Cephalopods beak. (a): Landmark positions. (b): Shape description—A: Rostral tip, B: Jaw angle, C: Wing, D: Lateral wall, E: Crest, F: Hood, G: Rostrum.
Figure 1. Shape and digitized landmarks (hollow dots) and semi-landmarks (solid dots) of the Cephalopods beak. (a): Landmark positions. (b): Shape description—A: Rostral tip, B: Jaw angle, C: Wing, D: Lateral wall, E: Crest, F: Hood, G: Rostrum.
Animals 13 00752 g001
Animals 13 00752 g002 550
Figure 2. Beak centroid size variation of four cephalopods in different ontogenetic stages. (a): Upper beak; (b): Lower beak. I: Immature stages; M: Mature stages; A: Octopus minor; B: Uroteuthis edulis; C: Sepia esculenta; D: Sthenoteuthis oualaniensis.
Figure 2. Beak centroid size variation of four cephalopods in different ontogenetic stages. (a): Upper beak; (b): Lower beak. I: Immature stages; M: Mature stages; A: Octopus minor; B: Uroteuthis edulis; C: Sepia esculenta; D: Sthenoteuthis oualaniensis.
Animals 13 00752 g002
Animals 13 00752 g003a 550Animals 13 00752 g003b 550
Figure 3. Results of principal component analysis (PCA) of upper beak of the four cephalopods, showing the first principal component (PC1) versus (PC2) and (PC3) versus (PC4) shape variation with 95% ellipse confidence intervals of immature (solid line) and mature (dashed line) beaks. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 3. Results of principal component analysis (PCA) of upper beak of the four cephalopods, showing the first principal component (PC1) versus (PC2) and (PC3) versus (PC4) shape variation with 95% ellipse confidence intervals of immature (solid line) and mature (dashed line) beaks. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g003aAnimals 13 00752 g003b
Animals 13 00752 g004a 550Animals 13 00752 g004b 550
Figure 4. Thin-plate spline deformation grids of upper beak of the four cephalopods considered in different ontogenetic stages. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 4. Thin-plate spline deformation grids of upper beak of the four cephalopods considered in different ontogenetic stages. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g004aAnimals 13 00752 g004b
Animals 13 00752 g005a 550Animals 13 00752 g005b 550
Figure 5. Correlativity of the regression scores of upper beak shape and predicted values versus log(Centroid Size). F–A: Immature, female; F–B: Mature, female; M–A: Immature, male; M–B: Mature, male. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 5. Correlativity of the regression scores of upper beak shape and predicted values versus log(Centroid Size). F–A: Immature, female; F–B: Mature, female; M–A: Immature, male; M–B: Mature, male. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g005aAnimals 13 00752 g005b
Animals 13 00752 g006a 550Animals 13 00752 g006b 550Animals 13 00752 g006c 550
Figure 6. Results of principal component analysis (PCA) of lower beak of the four cephalopods, showing the first principal component (PC1) versus (PC2) and (PC3) versus (PC4) shape variation with 95% ellipse confidence intervals of immature (solid line) and mature (dashed line). (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 6. Results of principal component analysis (PCA) of lower beak of the four cephalopods, showing the first principal component (PC1) versus (PC2) and (PC3) versus (PC4) shape variation with 95% ellipse confidence intervals of immature (solid line) and mature (dashed line). (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g006aAnimals 13 00752 g006bAnimals 13 00752 g006c
Animals 13 00752 g007a 550Animals 13 00752 g007b 550
Figure 7. Thin-plate spline deformation grids of lower beak of the four cephalopods considered in different ontogenetic stages. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 7. Thin-plate spline deformation grids of lower beak of the four cephalopods considered in different ontogenetic stages. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g007aAnimals 13 00752 g007b
Animals 13 00752 g008 550
Figure 8. Correlativity of the regression scores of lower beak shape and predicted values versus log(Centroid Size). F–A: Immature, female; F–B: Mature, female; M–A: Immature, male; M–B: Mature, male. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Figure 8. Correlativity of the regression scores of lower beak shape and predicted values versus log(Centroid Size). F–A: Immature, female; F–B: Mature, female; M–A: Immature, male; M–B: Mature, male. (a) Octopus minor. (b) Uroteuthis edulis. (c) Sepia esculenta. (d) Sthenoteuthis oualaniensis.
Animals 13 00752 g008
Table
Table 1. Sample information of four cephalopods in this study.
Table 1. Sample information of four cephalopods in this study.
SpeciesOntogenetic StagesQuantity of Samples (F, M)Mantle Length (mm) (F, M)
Octopus minorImmature97 (37, 60)16–69, 12–74
Mature102 (39, 63)20–72, 20–79
Uroteuthis edulisImmature110 (55, 55)92–190, 86–182
Mature95 (53, 42)103–240, 100–232
Sepia esculentaImmature120 (60, 60)51–129, 56–139
Mature78 (49, 29)62–160, 63–154
Sthenoteuthis oualaniensisImmature118 (60, 58)101–171, 100–189
Mature100 (60, 40)104–296, 109–187
Table
Table 2. Procrustes ANOVA of the shapes of upper beaks of the four cephalopods considered in different ontogenetic stages.
Table 2. Procrustes ANOVA of the shapes of upper beaks of the four cephalopods considered in different ontogenetic stages.
Octopus minordfSSMSFZp
Size10.05560.055610.91165.02870.001
Ontogenetic stages10.01090.01092.13951.81420.035
Sex10.01600.01603.14752.46760.009
Size × ontogenetic stages10.01170.01172.29592.03830.024
Size × sex10.00430.00430.8494−0.14840.549
Ontogenetic stages × sex10.00790.00791.54911.16510.125
Size × ontogenetic stages × sex10.00520.00521.02550.26760.385
Residuals1910.97260.0051
Total1981.0843
Uroteuthis edulisdfSSMSFZp
Size10.03350.03356.31053.40490.001
Ontogenetic stages10.01340.01342.52572.48510.007
Sex10.01030.01031.93451.92920.029
Size × ontogenetic stages10.01080.01082.04151.55200.066
Size × sex10.01590.01593.00312.10990.021
Ontogenetic stages × sex10.01260.01262.38242.32260.012
Size × ontogenetic stages × sex10.01270.01272.38791.85850.030
Residuals1971.04550.0053
Total2041.1548
Sepia esculentadfSSMSFZp
Size10.02860.02868.36904.96280.001
Ontogenetic stages10.02520.02517.35395.11750.001
Sex10.00470.00471.37570.98430.172
Size × ontogenetic stages10.01080.01083.16032.81670.003
Size × sex10.00690.00692.01571.85080.038
Ontogenetic stages × sex10.00190.00190.5642−1.10750.874
Size × ontogenetic stages × sex10.00260.00260.7575−0.41600.661
Residuals1900.64970.0034
Total1970.7304
Sthenoteuthis oualaniensisdfSSMSFZp
Size10.02950.02958.93395.43940.001
Ontogenetic stages10.00370.00371.10510.48860.318
Sex10.00660.00661.99382.01650.021
Size × ontogenetic stages10.00250.00240.7398−0.55310.698
Size × sex10.00440.00441.33450.85170.190
Ontogenetic stages × sex10.00480.00481.43961.10140.134
Size × ontogenetic stages × sex10.00550.00551.66371.47380.068
Residuals2100.69410.0033
Total2170.7510
df: degrees of freedom; SS: sum of squares; MS: mean squares’ F: test statistics; Z: effect size. Bold p values indicate significant at α = 0.05.
Table
Table 3. Procrustes ANOVA of the shapes of lower beaks of the four cephalopods considered in different ontogenetic stages.
Table 3. Procrustes ANOVA of the shapes of lower beaks of the four cephalopods considered in different ontogenetic stages.
Octopus minordfSSMSFZp
Size10.03990.03997.59015.98240.001 *
Ontogenetic stages10.00690.00691.31620.84050.193 ns
Sex10.00910.00911.72321.52170.064 ns
Size × ontogenetic stages10.00660.00661.24740.73340.235 ns
Size × sex10.00640.00641.21920.70830.244 ns
Ontogenetic stages × sex10.00820.00821.55491.24350.114 ns
Size × ontogenetic stages × sex10.00660.00661.26230.76950.219 ns
Residuals1921.00840.0053
Total1991.0919
Uroteuthis edulisdfSSMSFZp
Size10.03500.035010.14276.40090.001 *
Ontogenetic stages10.00700.00702.02071.97040.025 *
Sex10.00510.00511.47651.19450.133 ns
Size × ontogenetic stages10.00430.00431.25190.79940.210 ns
Size × sex10.00320.00320.9257−0.06210.523 ns
Ontogenetic stages × sex10.00840.00842.42452.41160.008 *
Size × ontogenetic stages × sex10.00630.00621.81121.73960.039 *
Residuals1970.67930.0034
Total2040.7484
Sepia esculentadfSSMSFZp
Size10.01960.01965.29654.21390.001 *
Ontogenetic stages10.01310.01313.53093.34740.001 *
Sex10.00940.00932.52192.27540.011 *
Size × ontogenetic stages10.01460.01463.92793.32260.001 *
Size × sex10.00430.00421.14610.53420.305 ns
Ontogenetic stages × sex10.00420.00421.12290.48100.320 ns
Size × ontogenetic stages × sex10.00750.00752.03091.78680.030 *
Residuals1900.70440.0037
Total1970.7770
Sthenoteuthis oualaniensisdfSSMSFZp
Size10.03150.03156.70914.50430.001 *
Ontogenetic stages10.01110.01112.35712.05590.017 *
Sex10.00430.00430.9091−0.01410.502 ns
Size × ontogenetic stages10.00470.00470.99140.18590.435 ns
Size × sex10.01510.01513.21572.75760.002 *
Ontogenetic stages × sex10.01030.01022.18532.04850.024 *
Size × ontogenetic stages × sex10.00600.00601.28250.78390.226 ns
Residuals2100.98500.0047
Total2171.0677
df: degrees of freedom, SS: sum of squares, MS: mean squares, F: test statistics, Z: effect size, ns: not significant. * p: significant at α = 0.05.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.; Fang, Z. Ontogenetic Variation and Sexual Dimorphism of Beaks among Four Cephalopod Species Based on Geometric Morphometrics. Animals 2023, 13, 752. https://doi.org/10.3390/ani13040752

AMA Style

Wang C, Fang Z. Ontogenetic Variation and Sexual Dimorphism of Beaks among Four Cephalopod Species Based on Geometric Morphometrics. Animals. 2023; 13(4):752. https://doi.org/10.3390/ani13040752

Chicago/Turabian Style

Wang, Chao, and Zhou Fang. 2023. "Ontogenetic Variation and Sexual Dimorphism of Beaks among Four Cephalopod Species Based on Geometric Morphometrics" Animals 13, no. 4: 752. https://doi.org/10.3390/ani13040752

APA Style

Wang, C., & Fang, Z. (2023). Ontogenetic Variation and Sexual Dimorphism of Beaks among Four Cephalopod Species Based on Geometric Morphometrics. Animals, 13(4), 752. https://doi.org/10.3390/ani13040752

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop