Chinese Journal of Oceanology and Limnology   2018, Vol. 36 issue(3): 981-989     PDF       
http://dx.doi.org/10.1007/s00343-018-7082-0
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

ZHAO Bo(赵博), LIU Jinhu(刘金虎), SONG Junjie(宋骏杰), CAO Liang(曹亮), DOU Shuozeng(窦硕增)
Otolith shape analysis for stock discrimination of two Collichthys genus croaker (Pieces: Sciaenidae, ) from the northern Chinese coast
Chinese Journal of Oceanology and Limnology, 36(3): 981-989
http://dx.doi.org/10.1007/s00343-018-7082-0

Article History

Received Mar. 22, 2017
accepted in principle Jan. 19, 2017
accepted for publication Mar. 13, 2017
Otolith shape analysis for stock discrimination of two Collichthys genus croaker (Pieces: Sciaenidae, ) from the northern Chinese coast
ZHAO Bo(赵博)1,3, LIU Jinhu(刘金虎)1,2, SONG Junjie(宋骏杰)1,3, CAO Liang(曹亮)1,2, DOU Shuozeng(窦硕增)1,2,3     
1 CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2 Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China;
3 University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: The otolith morphology of two croaker species (Collichthys lucidus and Collichthys niveatus) from three areas (Liaodong Bay, LD; Huanghe (Yellow) River estuary, HRE; Jiaozhou Bay, JZ) along the northern Chinese coast were investigated for species identification and stock discrimination. The otolith contour shape described by elliptic Fourier coefficients (EFC) were analysed using principal components analysis (PCA) and stepwise canonical discriminant analysis (CDA) to identify species and stocks. The two species were well differentiated, with an overall classification success rate of 97.8%. And variations in the otolith shapes were significant enough to discriminate among the three geographical samples of C. lucidus (67.7%) or C. niveatus (65.2%). Relatively high mis-assignment occurred between the geographically adjacent LD and HRE samples, which implied that individual mixing may exist between the two samples. This study yielded information complementary to that derived from genetic studies and provided information for assessing the stock structure of C. lucidus and C. niveatus in the Bohai Sea and the Yellow Sea.
Keywords: otolith     size descriptors     elliptic Fourier coefficients (EFC)     stock structure     Collichthys lucidus     Collichthys niveatus    
1 INTRODUCTION

The otolith is a calcium carbonate structure located in the inner ear of teleosts, and it serves as a mechanoreceptor that stimulates sensory cells in the inner labyrinth system (Fay and Popper, 2000). The core of the otolith forms during embryogenesis, and the calcium carbonate, in addition to trace elements in the environment, is deposited on the core as fish grow. Once deposited, the otolith materials cannot be degraded or altered (Thresher, 1999). Therefore, the otolith is considered to be a book that records the fish's physiological, morphological and environmental information on daily and annual scales. Otolith analysis has been widely used for age determination, species or stock identification, and growth, life history and migration analyses of fish (Campana and Casselman, 1993; Campana, 1999; Ferguson et al., 2011; Avigliano et al., 2015; Zhang et al., 2016; Avigliano et al., 2017).

Otolith shape analysis is a cheap and time-efficient method that is used for fish stock identification (Castonguay et al., 1991; Zhang et al., 2014; Izzo et al., 2017). With the development of digital photographic technology, landmark-based morphometrics (e.g., Procrustes analysis and thinplate splines) and boundary outline methods (e.g., Fourier transform and wavelet transform) were introduced as part of the otolith shape analysis (Monteiro et al., 2005; Parisi-Baradad et al., 2005; Ponton, 2006; Tracey et al., 2006; Lord et al., 2012; Tuset et al., 2015). Among these methods, elliptical Fourier transform (EFT) characterizes outline shapes with a series of harmonics, each of which is represented by four coefficients. The accuracy of describing the shape outline is higher as the number of harmonics used increases (Lord et al., 2012). Therefore, EFT is frequently used in otolith analysis for species or stock identification due to its relatively high efficiency to describe outline information (Stransky and MacLellan, 2005; Tracey et al., 2006; Yu et al., 2014).

The croakers Collichthys lucidus and C. niveatus are small, demersal fish species that are widely distributed southwards from the Bohai Sea, the Yellow Sea and the East China Sea. They usually spawn in coastal waters from April to July, feed around the spawning grounds and then move to deeper water for over-wintering (He and Li, 1988). With the decline of some commercial fish populations, the croakers have become important fishery species in China in recent years (Hu et al., 2015). Over the past decade, the annual landings of the two croakers from the Chinese coast were about 250 000 t, peaking at 300 000 t in 2014 (The People's Republic of China Ministry of Agriculture, Fisheries Bureau, 2006–2016). Despite their economic importance, the stock structures of the two croakers have been relatively little investigated.

For C. lucidus, two recent studies employing mitochondrial DNA sequencing suggested that the croaker in the Chinese coastal water can be classified into two phylogenetic clades (i.e., the northern clade in the Bohai Sea and the Yellow Sea and the southern clade in the East China Sea), with the geographical boundary in Zhoushan (Yin, 2013) or Wenzhou (Zhao et al., 2015), which are two adjacent spawning grounds for the croakers in the northern East China Sea. However, homogeneous mtDNA analysis has not excluded the existence of local populations within either phylogenetic clade. Take the northern clade for example, the Bohai Sea is a shallow semi-closed sea that has limited water exchange with the Yellow Sea (Bian et al., 2016). Moreover, the discharges of large rivers into the seas, the water temperatures, climatology elements, sediment types and biotic factors of the Bohai Sea and Yellow Sea are different, which may lead to unique geographical niches for the croakers (Sun, 2005; Tang, 2005). Due to the significant ecological differences and potential geographical segregation in the northern clade of the croaker, geographical stocks that have limited individual mixture may exist within the northern phylogenetic clade. To test this hypothesis, C. lucidus of the northern clade from three most concentrated spawning grounds (Liaodong Bay and Huanghe (Yellow) River estuary of the Bohai Sea; Jiaozhou Bay of the Yellow Sea) were used for stock discrimination by otolith shape analysis.

Figure 1 Sampling sites of the two croakers LD: Liaodong Bay; HRE: Huanghe River estuary; JZ: Jiaozhou Bay.

C.niveatus is a sibling specie of C. lucidus that usually co-occurs in these regions. But information about its stock structure is not yet available. For interspecific comparison and testing the hypothesis that geographical stocks may exist in the north China seas, stock discrimination of the C. niveatus from the three sampling sites was also conducted by otolith shape analysis in this study.

2 MATERIAL AND METHOD 2.1 Sample collection

The otoliths used for this study were from an otolith collection at the Institute of Oceanology, Chinese Academy of Sciences. These otoliths were collected during fishery research surveys in the Bohai Sea (Liaodong Bay samples, LD; Huanghe River estuary samples, HRE) and the Yellow Sea (Jiaozhou Bay samples, JZ), as described in Dou et al. (2012). The numbers of occipital spines were used to differentiate between the two croakers, C. lucidus has 4–6 occipital spines while C. niveatus has only 2 occipital spines (Chen, 1991). In short, the basic biological information of the fish samples (e.g., standard length SL, weight, sex and maturity stage) was measured immediately after the fish were captured. The samples were then frozen until further analysis. In the laboratory, both sides of the sagittal otoliths were removed from each fish, washed and stored dry in glass vials.

A total of 493 specimens (260 C. lucidus and 233 C. niveatus; 110–130 mm in total length for both species) were used for the otolith shape analysis (Table 1). No significant differences in fish length were detected among the sampling sties (ANOVA, P > 0.05 for both species), assuming that the otolith shape variation caused by fish length could be reduced to some extent (Cardinale et al., 2004).

Table 1 Basic information about the two croakers from the three sampling sites*
2.2 Otolith image capture and morphological analysis

Prior to imaging, the otoliths were cleaned in distilled water using an ultrasonic cleaner to eliminate the tissues adhered to the surfaces. The otoliths were then dried for 12 hours at 30℃ in an oven and weighed to the nearest 0.01 mg. The digital image of the left otolith of each individual fish was captured using a video camera linked to a binocular microscope (ACT- 2, Nikon SMZ1000; Japan). During the imaging, the otoliths were placed on a black background, with the sulcus acusticus facing down (Fig. 2).

Figure 2 The left otoliths of the two croakers

Image-Pro Plus 6.0 was used to measure the size descriptors (SDs; area, length, width and perimeter) of the otoliths. SHAPE 1.2 was used to extract the elliptic Fourier descriptors (EFDs) from the images of the otoliths. The coefficients were normalized with the program package to be invariant to otolith size and orientation. Additionally, the Fourier power (FP) spectrum was calculated according to Crampton (1995) to determine the amplitude number of harmonics required for the best description of the otolith outline.

The cumulative variance (FPc) was calculated to define the adequate number of harmonics to describe more than 99% of the cumulative power. The first seven harmonics and the first eight harmonics were adequate to reach the threshold of 99% for C. lucidus and C. niveatus, respectively. To balance the number of variables for interspecific comparison, the first seven harmonics were taken into account for both species (i.e., 7×4=28 Fourier coefficients; Fig. 3). Since the normalization of the EFDs (NEFDs) translated the first three coefficients into fixed values (i.e., c1=1, and c2=c3=0), the number of Fourier coefficients was reduced to 25 (c4c28) in the present study (Tracey et al., 2006).

Figure 3 Otolith shape outlines used in elliptical Fourier analysis
2.3 Data analysis

Prior to the shape analysis, the normality and homogeneity of each otolith morphological variable were first examined using Kolmogoroy-Smirnov (KS) and Levene's tests, respectively. The variables that did not meet the normality or homogeneity tests were rank-transformed prior to statistical analyses. NEFDs of female and male individuals were compared via MANCOVA analysis.

Analysis of covariance (ANCOVA) was conducted on the otolith morphological variables to check whether interactions existed between fish length and sampling sites. If a significant interaction was detected on a variable (P < 0.05, groups with unequal slopes), then the variable was excluded from subsequent analyses because it could not be scaled for fish length in a consistent manner (Begg and Brown, 2000; Cardinale et al., 2004). The remaining variables were scaled using allometric scaling method if they were significantly correlated with fish length.

The scaling was conducted based on an allometric growth model, Y=aXb, where both a and b are constants. After a series of transformations and substitutions (for details refer to Lleonart et al., 2000), the data standardization equation was described as follows (Lleonart et al., 2000):

which is independent of a, and where X0 is the average length of the samples.

A principal component analysis (PCA) was conducted on the NEFDs based on the variancecovariance matrix to decrease the dimensions of the variables. Principal components with an eigenvalue equal to, or greater than, the mean eigenvalue were used for subsequent analyses. A canonical discriminant analysis (CDA) was conducted on the selected principal components scores for interspecific and intraspecific identifications. As the most commonly used analysis for discriminating among samples with multiple variables, CDA can produce linear combinations that maximize the ratio of 'betweengroup' to 'pooled with-group' sample variances. A total of three CDA trials were conducted in this study. Homogeneity of the within-group covariance matrices was tested to select either a linear or a quadratic model. Since the matrices were not homogeneous (Box's M test, P < 0.05) in all cases, quadratic discriminant functions were computed in the three CDA trials. The first trial was conducted to compare the otolith shape variation between C. lucidus and C. niveatus. The other two trials were performed independently to discriminate among the stocks of each species. In all of the trials, a stepwise method was used to identify the effective variables for inclusion in the discriminant functions. Jack-knife cross-validation was performed to generate an unbiased estimation of classification success rates. To provide an intuitive understanding of the morphological differences among sample sites, the otolith SDs and otolith weight were analysed with one-way ANOVA.

All of the statistical tests were conducted using SAS 9.4 for Windows, and the differences were considered significant at P < 0.05.

3 RESULT

No significant differences were observed between the females and the males (MANCOVA, P > 0.05). The left otoliths of both sexes were pooled and used in the statistical analyses.

The main components and coefficient used for the CDA on interspecific identification and stock discrimination are listed in Appendix Table A1.

Appendix Table A1 Standardized canonical discriminant functions and main components of each principle scores used in each CDA trial

The results of one-way ANOVA conducted on the SDs and otolith weight showed that all the variables except for otolith weight differed significantly among sample sites (P < 0.05; Table 2).

Table 2 Student-Newman-Keuls mean comparisons of the SDs and otolith weight among sample sites of the two croakers
3.1 Interspecific identification

CDA conducted on polling samples of C. lucidus and C. niveatus reached a high classification success rate of 97.8%. Principle component 1 (PC1) which derived from the NEFDs contributed most to the discriminant function 1 (F1; Appendix Table A1). 2.3% C. lucidus and 2.2% C. niveatus were misclassified into each other.

3.2 Stock discrimination

In C. lucidus, F1 and F2 explained 84.96% and 15.04% of the between-group variance, respectively, which was significant for the distinction of the three sample sites (P < 0.05 in both cases). Among the discriminant variables, PC1 derived from the NEFDs contributed the most to F1, whereas PC6 contributed the most to F2. The CDA function coefficients and the main components of each principle score used are presented in Appendix Table A1. The CDA produced an overall classification success rate of 67.7%, with the highest rate for the LD sample (88.5%, prior probability 0.535), followed by the JZ sample (50.8%, prior probability 0.258) and the HRE sample (35.2%, prior probability 0.208). The misclassification rates between the LD and YER samples (i.e., the geographically adjacent samples in the Bohai Sea) were very high. For example, of the 16 misclassified LD individuals (11.5%, n=139), 11 individuals were assigned to the HRE sample, which is in contrast to only 5 assigned to the JZ sample. Similarly, of the 35 misclassified HRE individuals (64.8%, n=54), 28 individuals were assigned to the LD sample, and only 7 was assigned to the JZ sample (Table 3; Fig. 4).

Table 3 Results of CDA running on the stock identification of the two croakers
Figure 4 Scatter plots showing the CDA scores determined from the intraspecific identification a. C. lucidus; b. C. niveatus.

In C. niveatus, F1 and F2 explained 63.2% and 36.8% of the between-group variance, respectively, which was significant for the distinction of the three sample sites (P < 0.05 in both cases). Among the discriminant variables, PC3 contributed the most to F1, whereas PC4 contributed the most to F2. The CDA function coefficients and main components of each principle score used are presented in Appendix Table A1. The CDA produced an overall classification success rate of 65.2%, with the highest rate for the JZ sample (70.6%, prior probability 0.292), followed by the HRE sample (69.1%, prior probability 0.416) and the LD sample (54.4%, prior probability 0.292). Similar to the results in C. lucidus, most of the misclassification occurred between the LD and HRE samples. For example, of the 30 misclassified HRE individuals (30.9%, n=97), 17 individuals were assigned to the LD sample (Table 3; Fig. 4).

4 DISCUSSION

Generally, interspecific differences in otolith shapes are primarily determined by the genetics of the fish and can be used for identifying fish species (Lombarte and Lleonart, 1993; Cardinale et al., 2004). These differences in shape provide a useful method for species identification, particularly for identifying fossil samples (Reichenbacher et al., 2007) or dietary items in a stomach content analysis (Škeljo and Ferri, 2012). However, interspecific differences in otolith shapes are not always significant enough to allow for identifying species by visually observing their shape contours, particularly among species of the same genus. In such cases, an otolith shape analysis is required to detect the fine, interspecific differences needed for species identification. In the present study, the two croaker species were well identified using otolith shape descriptors.

Many studies have proven the validity of otolith shape analysis for fish stock identification (Campana and Casselman, 1993; Campana, 1999; Ferguson et al., 2011; Avigliano et al., 2017). In otolith shape analyses, the NEFDs usually show great ability to discriminate fish stocks, while the SDs (e.g., area, length, width and perimeter) can provide some understanding of the morphological differences among stocks (Cadrin and Friedland, 1999). The classification success of the intraspecies samples of the two croakers was not as high as that for interspecific identification. However, the results probably reflected the real scenarios of variations in otolith shapes among the geographical samples. For instance, most of the misclassification cases occurred between the two geographically adjacent samples of the Bohai Sea. Of the 260 C. lucidus samples, 84 individuals were misclassified, with 60.7% of the cases (51 of 84) occurring between the LD and HRE samples. In C. niveatus, 75.3% (61 of 81) of the misclassified LD and HRE samples were assigned to each other. Given that LD and HRE were geographically close in the Bohai Sea, individual mixing between the two samples may occur more frequently than with the JZ sample. In addition, the geo-hydrological and trophic conditions of the LD and HRE areas were similar, which could reduce the differences in the otolith shapes of the LD and HRE samples, resulting in relatively high misclassification between the two sample sites.

It is commonly believed that the otolith shape is regulated by the genetics of the fish and the environmental history they experience during their life history. The genetics involved in the biomineralization of otoliths may remain relatively stable at the intraspecific level due to the biological constraint that the form of an otolith should be in accordance with its function as a sound transducer (Lychakov and Rebane, 2000). The species-specific eternal otolith shape is genetically regulated, while the environmental effects on the otolith shape are mainly expressed at an intraspecific level (Vignon and Morat, 2010; Vignon, 2012). Currently, no information is available in the literature for the stock discrimination of C. niveatus on the coast of China. For C. lucidus, genetic analysis indicated that the croaker in the coastal waters of China could be classified into two distinct phylogenetic clades (the northern and the southern clades), and samples within the southern or northern groups retained a high degree of genetic homogeneity (Yin, 2013; Zhao et al., 2015). In the present study, the three sample sites were all located in the areas far north of the geographical boundaries. According to the genetic results in both studies mentioned above, the three geographical samples belonged to the northern clade, which should retain a high degree of genetic homogeneity. Therefore, it was reasonable to speculate that environmental factors might play a dominant role in causing the inter-group variations of the otolith shapes. The results of the stock identification supported this speculation. For example, most misclassification cases happened between the two geographically adjacent LD and HRE samples. The inter-group variations of otolith shapes among the samples could be primarily determined by the different environmental histories the croakers experienced rather than genetic differences. Previous studies had demonstrated that environmental factors (e.g., temperature and feeding conditions) influenced the growth of otoliths and, consequently, caused intergroup variations in the otolith shapes among the geographical samples (Campana and Casselman, 1993; Capoccioni et al., 2011; Vignon, 2012), which made the otolith shape analysis a useful technique for stock discrimination. The finding that significant variations of otolith shapes occurred among samples with high degrees of genetic homogeneity indicated that an otolith shape analysis could yield information complementary to that derived from genetic studies. Therefore, a holistic approach involving a broad spectrum of complementary techniques (e.g., otolith shape analysis, otolith microchemistry and genetics) should be applied for accurate stock discrimination and structure analysis, which is essential to effective fishery management (Begg and Waldman, 1999).

5 CONCLUSION

The two croakers C. lucidus and C. niveatus could be discriminated with high classification success rate of 97.8% using an otolith shape analysis. Modest classification success rates (67.7% and 65.2% for C. lucidus and C. niveatus, respectively) for discriminating among the three geographical samples of the two species were obtained. The results provided useful information for assessing the stock delimitation of the two croakers in the Bohai and the Yellow Seas. The intraspecific variations of the otolith shapes among the geographical samples of C. lucidus were probably mainly induced by environmental factors rather than genetics.

6 DATA AVAILABILITY

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

References
Avigliano E, Domanico A, Sánchez S, Volpedo A V. 2017. Otolith elemental fingerprint and scale and otolith morphometry in Prochilodus lineatus provide identification of natal nurseries. Fish. Res., 186: 1-10. DOI:10.1016/j.fishres.2016.07.026
Avigliano E, Velasco G, Volpedo A V. 2015. Use of lapillus otolith microchemistry as an indicator of the habitat of Genidens barbus from different estuarine environments in the southwestern Atlantic Ocean. Environ. Biol. Fish., 98(6): 1623-1632. DOI:10.1007/s10641-015-0387-3
Begg G A, Brown R W. 2000. Stock identification of haddock Melanogrammus aeglefinus on Georges bank based on otolith shape analysis. Trans. Am. Fish. Soc., 129(4): 935-945. DOI:10.1577/1548-8659(2000)129<0935:SIOHMA>2.3.CO;2
Begg G A, Waldman J R. 1999. An holistic approach to fish stock identification. Fish. Res., 43(1-3): 35-44. DOI:10.1016/S0165-7836(99)00065-X
Bian C W, Jiang W S, Pohlmann T, Sündermann J. 2016. Hydrography-physical description of the Bohai Sea. In: Harff J, Zhang H eds. Environmental Processes and the Natural and Anthropogenic Forcing in the Bohai Sea, Eastern Asia. J. Coast. Res., 74: 1-12.
Cadrin S X, Friedland K D. 1999. The utility of image processing techniques for morphometric analysis and stock identification. Fish. Res., 43(1-3): 129-139. DOI:10.1016/S0165-7836(99)00070-3
Campana S E, Casselman J M. 1993. Stock discrimination using otolith shape analysis. Can. J. Fish. Aquat. Sci., 50(5): 1062-1083. DOI:10.1139/f93-123
Campana S E. 1999. Chemistry and composition of fish otoliths:pathways, mechanisms and applications. Mar.Ecol. Prog. Ser., 188: 263-297. DOI:10.3354/meps188263
Capoccioni F, Costa C, Aguzzi J, Menesatti P, Lombarte A, Ciccotti E. 2011. Ontogenetic and environmental effects on otolith shape variability in three Mediterranean European eel (Anguilla anguilla, L.) local stocks. J. Exp.Mar. Biol. Ecol, 397(1): 1-7. DOI:10.1016/j.jembe.2010.11.011
Cardinale M, Doering-Arjes P, Kastowsky M, Mosegaard H. 2004. Effects of sex, stock, and environment on the shape of known-age Atlantic cod (Gadus morhua) otoliths. Can.J. Fish. Aquat. Sci., 61(2): 158-167. DOI:10.1139/f03-151
Castonguay M, Simard P, Gagnon P. 1991. Usefulness of Fourier analysis of otolith shape for atlantic mackerel(Scomber scombrus) stock discrimination. Can. J. Fish.Aquat. Sci., 48(2): 296-302. DOI:10.1139/f91-041
Chen D G. 1991. Fisheries Ecology in the Yellow Sea and the Bohai Sea. China Ocean Press, Beijing, China: p.273-276.
Crampton J S. 1995. Elliptic Fourier shape analysis of fossil bivalves:some practical considerations. Lethaia, 28(2): 179-186. DOI:10.1111/let.1995.28.issue-2
Dou S Z, Yu X, Cao L. 2012. Otolith shape analysis and its application in fish stock discrimination:a case study. Oceanol. Limnol. Sin., 43(4): 701-712.
Fay R R, Popper A N. 2000. Evolution of hearing in vertebrates:the inner ears and processing. Hear. Res., 149(1-2): 1-10. DOI:10.1016/S0378-5955(00)00168-4
Ferguson G J, Ward T M, Gillanders B M. 2011. Otolith shape and elemental composition:complementary tools for stock discrimination of mulloway (Argyrosomus japonicus) in southern Australia. Fish. Res., 110(1): 75-83. DOI:10.1016/j.fishres.2011.03.014
He B Q, Li H Q. 1988. Stock assessment of Collichthys lucidus in pearl river esturary. J. Fish. China, 12(2): 125-134.
Hu Y, Zhang T, Yang G, Zhao F, Hou J L, Zhang L Z, Zhuang P. 2015. Assessment of resource situation of Collichthys lucidus in coastal waters of the Yangtze estuary. Chin. J.Appl. Ecol., 26(9): 2867-2873.
Izzo C, Ward T M, Ivey A R, Suthers I M, Stewart J, Sexton S C, Gillanders B M. 2017. Integrated approach to determining stock structure:implications for fisheries management of sardine, Sardinops sagax, in Australian waters. Rev. Fish. Biol. Fish., 27(1): 267-284. DOI:10.1007/s11160-017-9468-z
Lleonart J, Salat J, Torres G J. 2000. Removing allometric effects of body size in morphological analysis. J. Theor.Biol., 205(1): 85-93. DOI:10.1006/jtbi.2000.2043
Lombarte A, Lleonart J. 1993. Otolith size changes related with body growth, habitat depth and temperature. Environ.Biol. Fish., 37(3): 297-306. DOI:10.1007/BF00004637
Lord C, Morat F, Lecomte-Finiger R, Keith P. 2012. Otolith shape analysis for three Sicyopterus (Teleostei:Gobioidei:Sicydiinae) species from New Caledonia and Vanuatu. Environ. Biol. Fish., 93(2): 209-222. DOI:10.1007/s10641-011-9907-y
Lychakov D V, Rebane Y T. 2000. Otolith regularities. Hear.Res., 143(1-2): 83-102. DOI:10.1016/S0378-5955(00)00026-5
Monteiro L R, Di Beneditto A P M, Guillermo L H, Rivera L A. 2005. Allometric changes and shape differentiation of sagitta otoliths in sciaenid fishes. Fish. Res., 74(1-3): 288-299. DOI:10.1016/j.fishres.2005.03.002
Parisi-Baradad V, Lombarte A, Garcia-Ladona E, Cabestany J, Piera J, Chic O. 2005. Otolith shape contour analysis using affine transformation invariant wavelet transforms and curvature scale space representation. Mar. Freshw.Res., 56(5): 795-804. DOI:10.1071/MF04162
Ponton D. 2006. Is geometric morphometrics efficient for comparing otolith shape of different fish species? J. Morphol., 267(6): 750-757. DOI:10.1002/(ISSN)1097-4687
Reichenbacher B, Sienknecht U, Küchenhoff H, Fenske N. 2007. Combined otolith morphology and morphometry for assessing taxonomy and diversity in fossil and extant killifish (Aphanius, †Prolebias). J. Morphol., 268(10): 898-915. DOI:10.1002/(ISSN)1097-4687
Škeljo F, Ferri J. 2012. The use of otolith shape and morphometry for identification and size-estimation of five wrasse species in predator-prey studies. J. Appl. Ichthyol., 28(4): 524-530. DOI:10.1111/jai.2012.28.issue-4
Stransky C, MacLellan S E. 2005. Species separation and zoogeography of redfish and rockfish (genus Sebastes) by otolith shape analysis. Can. J. Fish. Aquat. Sci., 62(10): 2265-2276. DOI:10.1139/f05-143
Sun X P. 2005. The Bohai Sea. In: Su J L ed. China Offshore Hydrology. China Ocean Press, Beijing, China. p. 182-190. (in Chinese)
Tang Y X. 2005. The Yellow Sea. In: Su J L ed. China Offshore Hydrology. China Ocean Press, Beijing, China. p. 193-202. (in Chinese)The People's Republic of China Ministry of Agriculture, Fisheries Bureau. 2006-2016. China Fisheries Yearbook. China Agriculture Press, Beijing, China. (in Chinese)
Thresher R E. 1999. Elemental composition of otoliths as a stock delineator in fishes. Fish. Res., 43(1-3): 165-204. DOI:10.1016/S0165-7836(99)00072-7
Tracey S R, Lyle J M, Duhamel G. 2006. Application of elliptical Fourier analysis of otolith form as a tool for stock identification. Fish. Res., 77(2): 138-147. DOI:10.1016/j.fishres.2005.10.013
Tuset V M, Imondi R, Aguado G, Otero-Ferrer J L, Santschi L, Lombarte A, Love M. 2015. Otolith patterns of rockfishes from the Northeastern Pacific. J. Morphol., 276(4): 458-469. DOI:10.1002/jmor.v276.4
Vignon M, Morat F. 2010. Environmental and genetic determinant of otolith shape revealed by a non-indigenous tropical fish. Mar. Ecol. Prog. Ser., 411: 231-241. DOI:10.3354/meps08651
Vignon M. 2012. Ontogenetic trajectories of otolith shape during shift in habitat use:interaction between otolith growth and environment. J. Exp. Mar. Biol. Ecol., 420-421: 26-32. DOI:10.1016/j.jembe.2012.03.021
Wu Z X, Chen X L. 1991. Collichthy lucidus richardson:a preliminary study on the age and phasic growth of Collichthy lucidus. J. Zhejiang College Fish., 10(2): 140-143.
Yin L N. 2013. Studies on genetic diversity and population structure of Collichthys lucidus. Ocean University of China, Qingdao. (in Chinese with English abstract)
Yu X, Cao L, Liu J H, Zhao B, Shan X J, Dou S Z. 2014. Application of otolith shape analysis for stock discrimination and species identification of five goby species (Perciformes:Gobiidae) in the northern Chinese coastal waters. Chin. J. Oceanol. Limnol., 32(5): 1060-1073. DOI:10.1007/s00343-015-4022-0
Zhang C, Ye Z J, Li Z G, Wan R, Ren Y P, Dou S Z. 2016. Population structure of Japanese Spanish mackerel Scomberomorus niphonius in the Bohai Sea, the Yellow Sea and the East China Sea:evidence from random forests based on otolith features. Fish. Sci., 82(2): 251-256. DOI:10.1007/s12562-016-0968-x
Zhang C, Ye Z J, Wan R, Ma Q Y, Li Z G. 2014. Investigating the population structure of small yellow croaker(Larimichthys polyactis) using internal and external features of otoliths. Fish. Res., 153: 41-47. DOI:10.1016/j.fishres.2013.12.012
Zhao M, Song W, Ma C Y, Zhang F Y, Jiang K J, Song Z M, Ma L B. 2015. Population genetic structure of Collichthys lucidus based on the mitochondrial cytochrome oxidase subunit I sequence. J. Fish. Sci. China, 22(2): 233-242.
References
Avigliano E, Domanico A, Sánchez S, Volpedo A V, 2017. Otolith elemental fingerprint and scale and otolith morphometry in Prochilodus lineatus provide identification of natal nurseries. Fish. Res., 186: 1–10. Doi: 10.1016/j.fishres.2016.07.026
Avigliano E, Velasco G, Volpedo A V, 2015. Use of lapillus otolith microchemistry as an indicator of the habitat of Genidens barbus from different estuarine environments in the southwestern Atlantic Ocean. Environ. Biol. Fish., 98(6): 1623–1632. Doi: 10.1007/s10641-015-0387-3
Begg G A, Brown R W, 2000. Stock identification of haddock Melanogrammus aeglefinus on Georges bank based on otolith shape analysis. Trans. Am. Fish. Soc., 129(4): 935–945. Doi: 10.1577/1548-8659(2000)129<0935:SIOHMA>2.3.CO;2
Begg G A, Waldman J R, 1999. An holistic approach to fish stock identification. Fish. Res., 43(1-3): 35–44. Doi: 10.1016/S0165-7836(99)00065-X
Bian C W, Jiang W S, Pohlmann T, Sündermann J. 2016. Hydrography-physical description of the Bohai Sea. In: Harff J, Zhang H eds. Environmental Processes and the Natural and Anthropogenic Forcing in the Bohai Sea, Eastern Asia. J. Coast. Res., 74: 1-12.
Cadrin S X, Friedland K D, 1999. The utility of image processing techniques for morphometric analysis and stock identification. Fish. Res., 43(1-3): 129–139. Doi: 10.1016/S0165-7836(99)00070-3
Campana S E, Casselman J M, 1993. Stock discrimination using otolith shape analysis. Can. J. Fish. Aquat. Sci., 50(5): 1062–1083. Doi: 10.1139/f93-123
Campana S E, 1999. Chemistry and composition of fish otoliths:pathways, mechanisms and applications. Mar.Ecol. Prog. Ser., 188: 263–297. Doi: 10.3354/meps188263
Capoccioni F, Costa C, Aguzzi J, Menesatti P, Lombarte A, Ciccotti E, 2011. Ontogenetic and environmental effects on otolith shape variability in three Mediterranean European eel (Anguilla anguilla, L.) local stocks. J. Exp.Mar. Biol. Ecol, 397(1): 1–7. Doi: 10.1016/j.jembe.2010.11.011
Cardinale M, Doering-Arjes P, Kastowsky M, Mosegaard H, 2004. Effects of sex, stock, and environment on the shape of known-age Atlantic cod (Gadus morhua) otoliths. Can.J. Fish. Aquat. Sci., 61(2): 158–167. Doi: 10.1139/f03-151
Castonguay M, Simard P, Gagnon P, 1991. Usefulness of Fourier analysis of otolith shape for atlantic mackerel(Scomber scombrus) stock discrimination. Can. J. Fish.Aquat. Sci., 48(2): 296–302. Doi: 10.1139/f91-041
Chen D G, 1991. Fisheries Ecology in the Yellow Sea and the Bohai Sea. China Ocean Press, Beijing, China: p.273–276.
Crampton J S, 1995. Elliptic Fourier shape analysis of fossil bivalves:some practical considerations. Lethaia, 28(2): 179–186. Doi: 10.1111/let.1995.28.issue-2
Dou S Z, Yu X, Cao L, 2012. Otolith shape analysis and its application in fish stock discrimination:a case study. Oceanol. Limnol. Sin., 43(4): 701–712.
Fay R R, Popper A N, 2000. Evolution of hearing in vertebrates:the inner ears and processing. Hear. Res., 149(1-2): 1–10. Doi: 10.1016/S0378-5955(00)00168-4
Ferguson G J, Ward T M, Gillanders B M, 2011. Otolith shape and elemental composition:complementary tools for stock discrimination of mulloway (Argyrosomus japonicus) in southern Australia. Fish. Res., 110(1): 75–83. Doi: 10.1016/j.fishres.2011.03.014
He B Q, Li H Q, 1988. Stock assessment of Collichthys lucidus in pearl river esturary. J. Fish. China, 12(2): 125–134.
Hu Y, Zhang T, Yang G, Zhao F, Hou J L, Zhang L Z, Zhuang P, 2015. Assessment of resource situation of Collichthys lucidus in coastal waters of the Yangtze estuary. Chin. J.Appl. Ecol., 26(9): 2867–2873.
Izzo C, Ward T M, Ivey A R, Suthers I M, Stewart J, Sexton S C, Gillanders B M, 2017. Integrated approach to determining stock structure:implications for fisheries management of sardine, Sardinops sagax, in Australian waters. Rev. Fish. Biol. Fish., 27(1): 267–284. Doi: 10.1007/s11160-017-9468-z
Lleonart J, Salat J, Torres G J, 2000. Removing allometric effects of body size in morphological analysis. J. Theor.Biol., 205(1): 85–93. Doi: 10.1006/jtbi.2000.2043
Lombarte A, Lleonart J, 1993. Otolith size changes related with body growth, habitat depth and temperature. Environ.Biol. Fish., 37(3): 297–306. Doi: 10.1007/BF00004637
Lord C, Morat F, Lecomte-Finiger R, Keith P, 2012. Otolith shape analysis for three Sicyopterus (Teleostei:Gobioidei:Sicydiinae) species from New Caledonia and Vanuatu. Environ. Biol. Fish., 93(2): 209–222. Doi: 10.1007/s10641-011-9907-y
Lychakov D V, Rebane Y T, 2000. Otolith regularities. Hear.Res., 143(1-2): 83–102. Doi: 10.1016/S0378-5955(00)00026-5
Monteiro L R, Di Beneditto A P M, Guillermo L H, Rivera L A, 2005. Allometric changes and shape differentiation of sagitta otoliths in sciaenid fishes. Fish. Res., 74(1-3): 288–299. Doi: 10.1016/j.fishres.2005.03.002
Parisi-Baradad V, Lombarte A, Garcia-Ladona E, Cabestany J, Piera J, Chic O, 2005. Otolith shape contour analysis using affine transformation invariant wavelet transforms and curvature scale space representation. Mar. Freshw.Res., 56(5): 795–804. Doi: 10.1071/MF04162
Ponton D, 2006. Is geometric morphometrics efficient for comparing otolith shape of different fish species? J. Morphol., 267(6): 750–757. Doi: 10.1002/(ISSN)1097-4687
Reichenbacher B, Sienknecht U, Küchenhoff H, Fenske N, 2007. Combined otolith morphology and morphometry for assessing taxonomy and diversity in fossil and extant killifish (Aphanius, †Prolebias). J. Morphol., 268(10): 898–915. Doi: 10.1002/(ISSN)1097-4687
Škeljo F, Ferri J, 2012. The use of otolith shape and morphometry for identification and size-estimation of five wrasse species in predator-prey studies. J. Appl. Ichthyol., 28(4): 524–530. Doi: 10.1111/jai.2012.28.issue-4
Stransky C, MacLellan S E, 2005. Species separation and zoogeography of redfish and rockfish (genus Sebastes) by otolith shape analysis. Can. J. Fish. Aquat. Sci., 62(10): 2265–2276. Doi: 10.1139/f05-143
Sun X P. 2005. The Bohai Sea. In: Su J L ed. China Offshore Hydrology. China Ocean Press, Beijing, China. p. 182-190. (in Chinese)
Tang Y X. 2005. The Yellow Sea. In: Su J L ed. China Offshore Hydrology. China Ocean Press, Beijing, China. p. 193-202. (in Chinese)The People's Republic of China Ministry of Agriculture, Fisheries Bureau. 2006-2016. China Fisheries Yearbook. China Agriculture Press, Beijing, China. (in Chinese)
Thresher R E, 1999. Elemental composition of otoliths as a stock delineator in fishes. Fish. Res., 43(1-3): 165–204. Doi: 10.1016/S0165-7836(99)00072-7
Tracey S R, Lyle J M, Duhamel G, 2006. Application of elliptical Fourier analysis of otolith form as a tool for stock identification. Fish. Res., 77(2): 138–147. Doi: 10.1016/j.fishres.2005.10.013
Tuset V M, Imondi R, Aguado G, Otero-Ferrer J L, Santschi L, Lombarte A, Love M, 2015. Otolith patterns of rockfishes from the Northeastern Pacific. J. Morphol., 276(4): 458–469. Doi: 10.1002/jmor.v276.4
Vignon M, Morat F, 2010. Environmental and genetic determinant of otolith shape revealed by a non-indigenous tropical fish. Mar. Ecol. Prog. Ser., 411: 231–241. Doi: 10.3354/meps08651
Vignon M, 2012. Ontogenetic trajectories of otolith shape during shift in habitat use:interaction between otolith growth and environment. J. Exp. Mar. Biol. Ecol., 420-421: 26–32. Doi: 10.1016/j.jembe.2012.03.021
Wu Z X, Chen X L, 1991. Collichthy lucidus richardson:a preliminary study on the age and phasic growth of Collichthy lucidus. J. Zhejiang College Fish., 10(2): 140–143.
Yin L N. 2013. Studies on genetic diversity and population structure of Collichthys lucidus. Ocean University of China, Qingdao. (in Chinese with English abstract)
Yu X, Cao L, Liu J H, Zhao B, Shan X J, Dou S Z, 2014. Application of otolith shape analysis for stock discrimination and species identification of five goby species (Perciformes:Gobiidae) in the northern Chinese coastal waters. Chin. J. Oceanol. Limnol., 32(5): 1060–1073. Doi: 10.1007/s00343-015-4022-0
Zhang C, Ye Z J, Li Z G, Wan R, Ren Y P, Dou S Z, 2016. Population structure of Japanese Spanish mackerel Scomberomorus niphonius in the Bohai Sea, the Yellow Sea and the East China Sea:evidence from random forests based on otolith features. Fish. Sci., 82(2): 251–256. Doi: 10.1007/s12562-016-0968-x
Zhang C, Ye Z J, Wan R, Ma Q Y, Li Z G, 2014. Investigating the population structure of small yellow croaker(Larimichthys polyactis) using internal and external features of otoliths. Fish. Res., 153: 41–47. Doi: 10.1016/j.fishres.2013.12.012
Zhao M, Song W, Ma C Y, Zhang F Y, Jiang K J, Song Z M, Ma L B, 2015. Population genetic structure of Collichthys lucidus based on the mitochondrial cytochrome oxidase subunit I sequence. J. Fish. Sci. China, 22(2): 233–242.