The traditional gobid taxonomy is dependent on morphological and meristic characters with varying degree of accuracy. The morphological characteristic of otolith is a rigid and species specific and not affected by external factors
Geffen (1983). Sometimes identification of species based on morphological characters of body leads to ambiguity; Hence the present study is based on otolith morphometric features. The left and rights are very similar, hence only right side otoliths used for comparison.
Otolith morphology
Morphological characteristics were studied in five gobid fishes (Fig 2) and the characters are described below -
B. dussumieri
Otolith cuneiform in shape, notch present on anterior side, crista superior larger than crista inferior, ostium and cauda not separated well, the dorsal side of otolith with groove, periphery of otolith rough on dorsal and smooth on ventral side.
T. vagina
Kidney shaped otolith, the posterior part broader than anterior, periphery smooth, a small notch observed on anterior side, external outgrowth on posterior side, core point thick, structure on otolith not clear.
O. roseus
Otolith bowl shaped, two flat notch on anterior side, three small grooves on ventral side, few annuli clearly visible near periphery on ventral side, core point thicker than periphery, peripheral margin smooth, anterior side broader than posterior
.
P. polynema
Otolith oval in shape, clear notch on anterior side, rostrum larger than anti rostrum, few redial suture on ventral side, anterior side broader than posterior, otolith more or less elongated, rostrum on ventral side with few grooves.
G. giuris
Otolith squared in shape, pointed extension on anterior side, some irregular precipitation of calcium on posterior side, otolith with moderate thickness, the periphery smooth.
In the present study, mean values of nine morphometric characters for the five species were compared (Table 2). Otolith of
G.
giuris species have maximum otolith area (9.352 mm
2), length (4.276 mm) and ellipticity (3.358). Maximum elongation (1.425), form factor (0.697) and rectangularity (0.717) was observed in
P.
polynema, maximum perimeter length (14.645 mm) in
B.
dussumieri, while maximum width (2.753 mm) and circularity (25.541) in
T.
vagina (Table 2).
Interrelationship between otolith characters
The regression of otolith parameters with length, stabilised on 95% level of confidence, showed isometric growth (2.994) of the otolith area of
O.
roseus with their length with lower p-value (2.128E-17) and higher goodness of fit (R
2=0.8012). The perimeter, width, ellipticity, rectangularity and elongation also fitted well, with significant relationship. The perimeter of
P.
polynema shows isometric growth (b=3.0071) with least P value (P=7.251E-31); apart from this area, width and ellipticity also showed good relationship with the otolith length at lower P-value and higher goodness of fit.
The otolith area of
T.
vagina revealed close to isometric growth (b=2.90) with their length at lower P-value (7.637E-10), while perimeter showed positive correlation. The otolith area of
B.
dussumieri revealed higher regression coefficient ‘b’ (4.23077) at lower P- value of coefficient b (2.901E-07) and higher goodness of fit (R
2=0.8557). The otolith of
G.
giuris showed positive allometric correlation of area (4.4937) and perimeter (3.71202) with length, at lower P- value 1.816E-42 and 7.672E-28, respectively, indicating higher goodness of fit (R2) 0.7946 and 0.9151, respectively. The width and ellipticity also showed good correlation at low P- value and high goodness of fit (Table 3).
SDFA
After factor analysis seven morphometric variables were sorted and subjected to Stepwise Discriminant Function Analysis (SDFA) to find out relevant characteristics for their discrimination power. The factor analysis showed that SDFA incorporated seven out of nine fed characters in the model. Comparative importance of these mentioned variables in the SDFA model, is indicated by their loading value on the function (Roots). The factor loading for perimeter, area and length on both root emphasises on its higher separating power in the comparison of characters (Table 4).
The factor structure matrix showed higher contribution of ellipticity and rectangularity to Root 1 and rectangularity elongation to Root 2 which stress their significance in discriminating the different species (Table 5).
A classification matrix, generated by SDFA, revealed a 97.173% correct classification. Out of 213 specimens, only 3 specimens of
T.
vagina were found with overlapping characters with
O.
roseus. Rest of the specimens were correctly (100%) classified (Table 6). An instance of misclassification between
T.
vagina and
O.
roseus showed the morphometric proximity of the species that was again affirmed by the less square Mahalanobis distance (9.022). Maximum squared mahalanobis distance in the morphometric term have been observed between
G.
giuris and
B.
dussumieri followed by
O.
roseus and G.
giuris pair (Table 7).
The adequacy of the variables in species differentiation and classification power of the SDFA model is also indicated by the scatterplot, plotted based on canonical scores in different cases of same species, where the values of the species were aggregated and well-spaced from other species of the group (Fig 3).
The results of the present study show that three methods (morphology, shape indices and step wise discriminant function analysis), applied for classifying sagittal otoliths, provide potential power for identification and discrimination of five species of family Gobiidae. The results of present study also show that otolith characteristics are species specific that can be applied for discrimination of species and, otolith can play a vital role in the study of species ecology
(Nazir and Khan, 2021; Reichenbacher et al., 2007; Volpedo and Echevarría, 2003; Paxton, 2000; Aguirre and Lombarte, 1999). The factors such as sound, hearing, depth, habitat, behaviour, swimming activity of species affect the shape of otolith
(Wilson, 1985; Lombarte et al., 2003; Volpedo and Echeverría, 2003; Sadighzadeh et al., 2014; Tuset et al., 2015; Cruz and Lombarte, 2004). The data was fitted well to the regression model (R≥0.90) for area of all species, for perimeter for
G.
giuris and
P.
polynema, for width
G.
giuris,
P.
polynema and
B.
dussumieri. The outcome of regression analysis revealed more than 95% level of confidence (P<.05) for area and width of all the species; perimeter of all species except
B.
dussumieri, rectangularity of all species except
B.
dussumieri and
G.
giuris, circularity of
B.
dussumieri and
P.
polynema. ellipticity of
O.
roseus,
P.
polynema and
G.
giuris, form factor of
P.
polynema and elongation of
O.
roseus and
P.
polynema. The area and perimeter of all the species and rectangularity of
T.
vagina and
O.
roseus showed positive correlation with the otolith length and out of these. area of
O.
roseus,
T.
vagina and perimeter of
P.
polynema showed isometric growth, while area of
G.
giuris and
B.
dussumieri and perimeter of
G.
giuris and
T.
vagina showed positive allometric growth.
In the present study, differences were observed in almost all morphometric parameters (shape indices). The morphological variation in otolith may be due the depth of inhabitation, skull size and somatic growth pattern
(Aguirre and Lombarte, 1999; Tuset et al., 2006).
The SDFA performed to find correct classification of all five gobid species, based on sagittal otolith shape, revealed 97.1831% accuracy.
Lombarte and Lleonart (1993) reported around 97% accuracy in salmodis from NW Atlantic.
Bani et al., (2013) reported 94.7% accuracy in discrimination of gobid species from Anzali coast, Guilan, Iran. There was no significant difference between left and right otolith in all five species, but both otoliths were not mirror image as also reported by
Harvey et al., (2000); Hunt (1979); Waessle et al., (2003).
Results encourage the use of these relationships as a tool in predicting fish size from the otoliths and in calculating the biomass during feeding studies and also these data can be used as baseline information for taxonomic discrimination of the species.