Morphological differences between species of the sea bass genus Lateolabrax (Teleostei, Perciformes), with particular emphasis on growth-related changes

Abstract Morphological differences, including growth-related changes, were examined in three morphologically similar East Asian sea bass species, Lateolabraxjaponicus, L.maculatus and L.latus. In many cases, body measurements indicated specific patterns of growth-related proportional changes. Lateolabraxlatus differed from the other two species in having greater body depth, caudal peduncle depth, caudal peduncle anterior depth, snout length, and upper and lower jaw length proportions. In particular, scatter plots for caudal peduncle anterior depth relative to standard length (SL) in that species indicated complete separation from those of the other two species, being a new key character for identification. Comparisons of L.japonicus and L.maculatus revealed considerable proportional differences in many length-measured characters, including fin lengths (first and second dorsal, caudal and pelvic), snout length, post-orbital preopercular width (POPW) and post-orbital length. In particular, snout length (SNL) and POPW proportions of the former were greater and smaller for specimens >200 and ≤ 200 mm SL, respectively. Because the scatter plots of these proportions for the two species did not overlap each other in either size range, identification of the species was possible using a combination of the two characters. In addition, scatter plots of the POPW / SNL proportion (%) of L.japonicus and L.maculatus were almost completely separated throughout the entire size range examined (border level 90%), a further aid to identification. The numbers of pored lateral line scales and scales above the lateral line tended to increase and decrease with growth, respectively, in L.japonicus, whereas scales below the lateral line and gill raker numbers tended to increase with growth in L.maculatus. Because the ranges of these meristic characters may therefore vary with specimen size, they are unsuitable for use as key characters. Accordingly, a new key is proposed for the genus Lateolabrax.


Introduction
The sea basses of the genus Lateolabrax (Lateolabracidae) are common East Asian coastal marine fishes (occasionally also occurring in fresh water). Bleeker (1854-57) established the genus for a single species, Lateolabrax japonicus (Cuvier, 1828), Katayama (1957) later describing a second species, Lateolabrax latus, from Japan. More recently, Yokogawa and Seki (1995) concluded that differences between the Japanese and Chinese forms of "L. japonicus" were sufficient for the Chinese form to be recognized as a distinct species, being referred to as "spotted sea bass" by Yokogawa and Tajima (1996). Finally, it was formally redescribed as Lateolabrax maculatus (McClelland, 1844) in Yokogawa's (2013b) revision, where Lateolabrax lyiuy (Basilewsky, 1855), which is incorrectly treated as valid and applied to the Chinese form (Kottelat 2013;Eschmeyer 2019), was regarded as a junior synonym of L. maculatus. At this point, three valid species of Lateolabrax are recognized (Fig. 1).
Lateolabrax latus has been distinguished from L. japonicus by having greater proportions of body and caudal peduncle depth (BD and CPD), more dorsal and anal fin rays (≥15 and ≥9, respectively), fewer scales below the lateral line (≤16) and possessing ventromandibular scale rows (VSRs) (Katayama 1957). Furthermore, the range of dorsal fin ray (DFR) counts in L. latus, which had been considered to not overlap that of L. japonicus, had become established as a key identification character (e.g., Katayama 1960aKatayama , 1965Katayama , 1984Hatooka 1993). However, subsequent finding of L. latus individuals with 14 DFRs [overlapping the range in L. japonicus (12-14 DFRs)] (Hatooka 2000(Hatooka , 2013Murase et al. 2012) made this character an incomplete key for identification. In addition, VSRs have not been adopted in recent keys proposed for Lateolabrax (Hatooka 2000(Hatooka , 2013, because they have been found in some specimens of the other two Lateolabrax species (Paxton and Hoese 1985;Hirota et al. 1999;Kang 2000;Murase et al. 2012). On the other hand, recent keys have included "caudal fin notch depth," being shallower in L. latus than in the other two species (Hatooka 2000(Hatooka , 2013, although the lack of any proportional information means that verification following examination of possible growth-related changes is necessary. Furthermore, proportional differences in BD and CPD appear to be based on the premise that their proportions are stable (isometric growth), although this has not been verified to date.
Lateolabrax maculatus has been characterized by many clear black spots on the body, but this character is also problematic as a few individuals of the species lack such spots (Yokogawa and Seki 1995), whereas some individuals of the other two Lateolabrax species have dots (Fig. 2). Although Yokogawa and Seki (1995) revealed differences between L. japonicus and L. maculatus in some meristic characters, including counts of lateral line scales, gill rakers and vertebrae, overlapping ranges between the two species result in no single character separating them completely. Proportional snout length (SNL), also recently used to separate the two species [SNL of L. maculatus relatively shorter than in L. japonicus (Hatooka 2000(Hatooka , 2013Yamada et al. 2007)], may also be problematic due to lack of proof of isometric growth. Furthermore, morphology of the first anal pterygiophore (arched and straight in L. japonicus and L. maculatus, respectively), proposed by Kang (2000), still needs to be validated due to possible growth-related changes.

Morphological measurements
Methods of measurements and counts followed Hubbs and Lagler (1970). Dimensions were taken with calipers (minimum scale 0.1 mm), with particular care for smaller specimens due to the effect of even a small error on the calculated proportion. The characters examined are listed with abbreviations in Table 1 and illustrated in Figure  3. New or uncommon length-measured characters included: post-orbital preopercular width (horizontal distance from orbit posterior edge to preopercle posterior margin), post-orbital length (distance from orbit posterior edge to opercle posterior angle), caudal peduncle anterior depth (distance between posterior ends of dorsal and anal fin bases), caudal fin notch depth (horizontal distance from bottom of notch to margin of naturally spread fin) and pectoral scaly area length (defined by Yokogawa and Seki 1995) (see Figure 3).
Scale row and paired fin ray counts were made on the left side of the body, whereas gill rakers were counted on the first gill arch on the right side by separating the upper and lower limbs of the gill arch. Because counts of pelvic fin-spine (P 2 FS) and soft rays (P 2 FRs) showed no variation (P 2 FS: 1, P 2 FRs: 5 in all specimens), these counts were omitted from the statistical analyses. Abdominal and caudal vertebrae were counted, and first anal fin pterygiophore morphology observed from radiographs.
Total numbers of recognizable black or faint spots / dots on the left side of the body and mid-dorsal aspect of the caudal peduncle ( Fig. 2) were counted. Dorsal head squamation [reported as differing between L. japonicus and L. maculatus (Yokogawa and Seki 1995)], was examined in all three species. Ventromandibular scale rows were also examined on the left side by separating the anterior and posterior parts following Murase et al. (2012), and their status recorded as present, vestigial or absent.

Statistical computations
For a length-measured dimension (LD), a growth-related proportional change pattern is given by the relationship between base dimension [e.g., standard length (SL) or head length (HL)] and the LD proportion (LD / SL or LD / HL). Because the relationship between SL (or HL) and LD is generally expressed by a power regression formula (LD = a SL b ) (allometric growth), the following formula was used (LD / SL = a SL b-1 ). Accordingly, power regressions were applied for the relationships between SL (or HL) and the LD proportions (Table 2). Characteristics that changed with growth were evaluated so as to determine if the changes were isometric or allometric, i.e., regressions between SL (or HL) and LD were transformed into natural logarithms (ln) (lnLD = a lnSL + b), and a t test was used to examine slope significance for the null hypothesis (a = 1), according to Zar (2010). When a differed significantly from 1, the character was considered to have changed allometrically, i.e., its proportion had increased or decreased with growth. Meristic counts (MC) were regressed using SL (MC = a SL + b), and a t test used to examine slope significance for the null hypothesis (a = 0) (Zar 2010). When a differed significantly from 0, the character was considered to have changed with growth. In addition, standard errors, which indicated data variation from the regression lines, were calculated during the above analyses (Zar 2010).
To examine inter-specific differences in length-measured characters, regressions between SL (or HL) and LD were also logarithm-transformed (lnLD = a lnSL + b), since most characters showed allometric growth (Table 3). Parameters of the regressions (a and b) were compared by analysis of covariance (ANCOVA) (t test), following the methods of Yamada and Kitada (2004).
Because some meristic counts tended to increase significantly with growth (Table  4), they were compared using the Mann-Whitney U test (Iwasaki 2006). Example numbers for the U test being >20 for all species, z values (instead of U values) for the normal distribution were calculated after correction for distribution continuity, following Iwasaki (2006).
In the above statistical inferences, due to multiple tests being applied simultaneously in each case, multiple comparisons were introduced for the t test results, risk percentages for the t values being corrected according to total test counts, using the Holm-Bonferroni method (Holm 1979).

Growth-related proportional changes Body characters
In the three Lateolabrax species, slopes of the logarithm-transformed regressions were significantly different from 1 (allometric growth) for most characters (Table 3), indicating that most of the body proportions changed with growth. Relationships between standard length (SL) and length-measured body proportions are shown graphically by species in Figure 4, those with prominent plot separation between species being shown with multiple specific plots in Figure 5.  Table 1.
Similar patterns of growth-related proportional changes common to the three species were observed for some characters, viz., significant positive allometric growth (proportions increased with growth) in body width and significant negative allometric growth (proportions decreased with growth) in head (HL) and pre-dorsal length (PDL), and second dorsal, anal and pelvic fin (longest ray) lengths (SDFL, AFL and P 2 FL), although patterns of the regression curves or plot distributions for the three spe-  Table 3). Differing specific growthrelated proportional changes were evident for some other characters, e.g., pre-anus length (PAL), isometric growth in L. japonicus, negative and positive allometric growth in L. maculatus and L. latus, respectively ( Fig. 4A-C, Table 3); and caudal fin notch depth (CFND), modestly and highly negative allometric growth in L. japonicus and L. latus, respectively, and isometric growth in L. maculatus ( Fig. 4G-I, Table 3). In the latter, however, despite specific growth-related patterns, ranges of the CFND / SL proportions taken over the entire range of SLs were similar to one another, viz., 2.0-8.4%, 1.9-7.4% and 2.9-7.9%, in L. japonicus, L. maculatus and L. latus, respectively (    . Relationships between standard length and proportions of some length-measured body characters of three Lateolabrax species. For character abbreviations, see Figure 3 and Table 1. Solid lines indicate power regression curves (parameters given in Table 2).

Cephalic characters
For length-measured dimensions (LD) of cephalic characters, SL-based (SL-LD / SL) and HL-based relationships (HL-LD / HL) are illustrated in pairs with multiple specific plots in Figure 6. In each species, significant allometric growth was recognized in most length-measured cephalic characters as well as length-measured body characters (Table 3). In particular, negative allometric growth was so significant for orbital diameter (OD) (very high t values, see Table 3) that the plots for each all formed typical arched curves (Fig. 6C, D), indicating rapid decrement of OD proportions with growth. Such acute relative OD decrement in the three species was clearly apparent from photographs ( Fig. 1). Growth-related proportional change patterns based on SL and HL were inconsistent with each other for some characters in L. japonicus and L. latus, e.g., snout length (SNL) of L. japonicus was isometric and positively allometric for SL and HL,  Table 1. Solid lines indicate power regression curves (parameters given in Table 2).

Figure 5.
Relationships between standard length and proportions of some length-measured body characters which exhibited prominent plot separation among three Lateolabrax species. For character abbreviations, see Figure 3 and Table 1. Solid lines indicate power regression curves (parameters given in Table 2) for each species. respectively; that of L. latus was negatively and positively allometric for SL and HL, respectively ( Fig. 6A, B, Table 3). While the patterns were consistent between the SL-and HL-based relationships in L. maculatus for all cephalic characters (Fig. 6A-P,  Table 3), allometric increment / decrement rates varied in the two-way relationships e.g., proportions of post-orbital preopercular width (POPW) decreased with growth acutely and slightly for SL and HL, respectively ( Fig. 6I, J, Table 3).
As well as some body characters, specific proportional change patterns were recognized for some characters, e.g., SL-based relationships of POPW, exhibiting isometric growth in L. japonicus, and positive and negative allometric growth in L. maculatus and L. latus, respectively (Fig. 6K, Table 3); and SNL, exhibiting isometric growth in L. japonicus, and high and modest negative allometric growth in L. maculatus and L. latus, respectively (Fig. 6A, Table 3).

Pectoral scaly area length
The relationship between SL and pectoral scaly area length (PSAL) in L. latus was well fitted to a power regression (like many other body and cephalic length-measured characters), the PSAL / SL proportion gradually decreasing with growth ( Fig. 4X, Table  2). In the other two species, however, proportional PSAL rapidly increased from the smallest specimens to a peak and thereafter gradually decreased (Fig. 4V, W), therefore being unsuitable for simple patterned regressions. Synchronous plotting for the two species showed the proportional PSAL of L. maculatus to be distinctly less than that of L. japonicus during the initial stage (< ca. 150 mm SL), although plots of the two species largely overlapped during the subsequent decreasing stage (Fig. 7). The proportional PSAL of L. latus during the former stage was much greater than in the other two species ( Fig. 4V-X).

Inter-specific differences Length-measured body and cephalic characters
Plot separation of L. latus from the other two species was prominent for vertical body dimensions of body depth (BD), caudal peduncle depth (CPD) and caudal peduncle anterior depth (CPAD), L. japonicus and L. maculatus both showing significant negative allometric growth, the degree of relative decrease being especially acute in the former. Although BD of L. latus showed slight negative allometric growth, CPD and CPAD were regarded as isometric ( Fig. 5A-C, Table 3). However, despite considerable plot separation of BD and CPD between L. latus and the other species, plots of the three species overlapped for the smaller size class (< ca. 200 mm SL) (Fig. 5A, B). CPAD plots for L. latus were entirely separated from those of the other two species (border level 15%) (Fig.  5C). Although similar plot separation for caudal peduncle length (CPL) in L. latus was Figure 6. Relationships between standard length or head length and proportions of length-measured cephalic characters of three Lateolabrax species. For character abbreviations, see Figure 3 and Table 1. Solid lines indicate power regression curves (parameters given in Table 2) for each species.  Table 1. Solid lines indicate power regression curves (parameters given in Table 2) for each species. also apparent, ranges of proportional CPL of the three species almost overlapped due to considerable variation in plot distribution in L. japonicus and L. maculatus ( Fig. 4G-I).
Plot separation of first and second dorsal (FDFL and SDFL), caudal (CFL) and pectoral (P 1 FL) fin lengths was also apparent between L. japonicus and L. maculatus ( Fig. 5E-H), the former showing significant negative allometric growth of each feature, whereas the latter exhibited isometric growth for all, except SDFL (Table 3). Proportions in the former were distinctively greater than in the latter in the smaller size class (< ca. 200 mm SL), although plots of the two species overlapped in the larger size class (> ca. 200 mm SL), since fin length proportions decreased and did not change with growth, respectively ( Fig. 5E, G, H). Such proportional differences in fin length in the smaller size class between the two species were clearly apparent from photographs ( Fig. 1A, C).
Upward plot separation of L. latus from the other two species was prominent for SNL and upper and lower jaw lengths (UJL and LJL), there being almost no overlap with L. maculatus and only modest overlap with L. japonicus (Fig. 6A, B, M-P). Plots of OD for L. latus were similarly upwardly separated from those of the other two species (Fig. 6C, D), especially in the HL-based graph (Fig. 6D). Post-orbital length (POL) plots for L. latus were shifted downward from those of the other two species (Fig. 6K, L), plot separation being more prominent in the HL-based graph (Fig. 6L).
On the other hand, plot separation between L. japonicus and L. maculatus was prominent for SNL, POPW and POL (Fig. 6A, B, I-L). SNL plots for the two species overlapped in the smaller size class (< ca. 200 mm SL), subsequently progressively separating with growth due to the proportional SNL of L. maculatus decreasing with growth (negative allometry), to a border level of ca. 7.7% (Fig. 6A) in the larger size class (> ca. 200 mm SL). This phenomenon was more apparent in the HL-based relationship because proportional SNL in L. japonicus increased with growth (positive allometry) (Fig. 6B, Table 3), unlike that for the SL-based relationship (isometric growth) (Fig. 6A, Table 3). Similar patterns were observed for POL, plots of the two species overlapping in the smaller size class (< ca. 200 mm SL), but subsequently separating to a certain extent with growth due to a proportional POL increase in L. maculatus (positive allometry) (Fig. 6K, L). In contrast, POPW plots of the two species were completely separated from each other in the smaller size class (< ca. 200 mm SL), having a border level of ca. 7.5%, but progressively overlapped with growth due to the proportional POPW of L. japonicus and L. maculatus increasing and decreasing with growth, respectively (Fig. 6I, J). POPW proportional to SNL is shown graphically in Figure 8. The SL-POPW / SNL regressions were positively allometric for L. japonicus and L. latus, and isometric for L. maculatus (Table 3). Plots for L. japonicus and L. maculatus were separated from each other almost entirely throughout all size ranges (border level 90%), following a slight plot overlap at ca. 100 mm SL (Fig. 8). In addition, plots for L. latus were displaced well downward from the other two species, despite some overlap with L. japonicus (Fig. 8).

Meristic characters
The t tests of regressions between SL and meristic counts (null hypothesis, slope = 0) proved significant for scales on (LLS) and above the lateral line (SAL) in L. japonicus, and scales below the lateral line (SBL) and gill raker counts [lower limb and total (LGR and TGR, respectively)] in L. maculatus (Table 4). Whereas SAL counts in L. japonicus tended to decrease with growth ( Fig. 9), having negative slope values (Table 4), the remaining characters tended to increase (Fig. 9, Table 4). No significant differences in any meristic characters were found in L. latus (Table 4), indicating that none changed with growth in that species. Figure 10 shows multiple specific frequency histograms for all meristic characters, L. latus clearly differing from the other two species in dorsal (DFR) and anal fin ray (AFR) counts (there being only slight range overlaps), as well as in pectoral fin ray (P 1 FR) and SBL counts, again with some range overlaps. Notably, DFRs (14) in L. latus had only a 7.4% overlap of the ranges of the other two species, the latter differing significantly in vertebral counts [caudal and total (CV and TV, respectively)] and ranges of LLS, LGR and TGR. However, no species had a meristic character count range that was entirely separated from those of the other species.

Spots / dots on lateral body region
Some examples of L. japonicus and L. latus had small and fine dots, respectively, on the lateral body region ( Fig. 2A, C), whereas L. maculatus usually had many clear black spots (Fig. 2B). In both of the former, dots appeared to be limited to some smaller specimens (Fig. 11A, C), the maximum sizes of specimens with dots being 260.6 mm SL (BSKU 100765) and 254.8 mm SL (KAUM-I. 66393), respectively. The t tests indicated significant regressions between SL and dot counts for the two species (null hypothesis, slope = 0 rejected), both indicating negative correlations (minus slope values) (Fig. 11A, C, Table 4). The proportions of dotted specimens of the total material examined were 35.6% and 46.3% (51.9 and 60.0% for specimens <250 mm SL) in L. japonicus and L. latus, respectively. In L. maculatus, spot counts were typically abundant (ca. 40 on average), but variable (absent in 4.9% of specimens) (Fig. 11B) and not related to body size, a t test (null hypothesis, slope = 0) indicating no significant regression between SL and spot counts (Table 4).

Squamation on dorsal head region
Post-juvenile specimens (> ca. 70 mm SL) of the three Lateolabrax species had a pair of scale rows (dorsocephalic scale rows, DSRs) extending forward from the inter-orbital area, which was densely covered with fine scales (Fig. 12). DSRs in L. japonicus and L. latus were well developed distally, with anterior edges always beyond the anterior nostril position (ANP) (Fig. 12A, B, E, F), and almost reaching the upper lip in large specimens of L. latus (Fig. 12F). On the other hand, DSRs in small specimens of L.  Table 2) for each species. . Relationships between standard length and some meristic characters which exhibited growthrelated changes in some Lateolabrax species. For character abbreviations, see Figure 3 and Table 1. Solid lines indicate linear regressions (parameters given in Table 4). maculatus were almost entirely restricted to the inter-orbital region, not extending beyond ANPs (Fig. 12C), although gradual development with growth resulted in DSRs extending beyond the ANP in specimens > ca. 150 mm SL (Fig. 12D).  Table 4).

Squamation on ventral head region
Some individuals of the three Lateolabrax species had a pair of ventromandibular scale rows (VSRs), VSR status by body size being summarized in Table 5. In L. japonicus, although VSRs were entirely absent in specimens ≤ 100 mm SL, a few ca. 150 mm SL had vestigial VSRs. Subsequently, the proportion of specimens with VSRs gradually increased with growth, those lacking anterior and posterior VSRs comprising 25.0% and 0%, respectively, of the largest size class (> 400 mm SL). VSRs were entirely absent in L. maculatus specimens < 200 mm SL, appearing in a few just over 200 mm SL. Subsequently, the proportion of specimens with VSRs gradually increased with growth, those without anterior and posterior VSRs comprising 36.4% and 0%, respectively, of the largest size class (> 400 mm SL). Although VSRs were absent in most L. latus specimens ≤ 100 mm SL, a few over 90 mm SL had incipient or established VSRs. Subsequently, the proportion of specimens with VSRs rapidly increased with growth, including most up to 300 mm SL and all > 300 mm SL. Notably, 100-200 mm SL specimens with VSRs showed greater development of the anterior portion, contrary to the developmental pattern displayed by the other two species. The prominence of VSR appearance was ranked: 1 L. latus, 2 L. japonicus, 3 L. maculatus.

Morphology of first anal pterygiophore
All three Lateolabrax species had a well-developed first anal pterygiophore (FAP), which comprised a short thin plate-like anterior part and a long thick spiny posterior part (Fig. 13). In L. japonicus, although the FAPs were straight in small specimens (< ca. 90 mm SL) (Fig. 13A), they became modestly arched in larger specimens (Fig.  13B-D), suggesting a growth-related morphological change. In contrast, the FAPs in L. maculatus remained straight (morphologically stable) regardless of body size (Fig.  13E-H). In L. latus, on the other hand, although the FAPs were straight in some specimens (Fig. 13I, K), they were slightly arched distally in others (Fig. 13J, L), thus showing neither growth-related morphological change nor morphological stability. As such, relationships between body size and FAP morphology were specifically unique.

Statistical differences
Analyses of covariance (ANCOVA) for regressions of logarithm-transformed lengthmeasured characters by pairwise comparisons for the three Lateolabrax species indi-  cated significant differences in the slopes or intercepts of all such characters (Table 6). In general, significance (t values) between L. japonicus and L. latus, and L. maculatus and L. latus were greater than those between L. japonicus and L. maculatus, suggesting that L. latus exhibited greater morphological differences from the other two species (Table 6). High significance levels between the species were apparent for the SNL-POPW relationship (t values for intercepts ca. 28-44), in which the scatter plots were almost entirely separated from one another (Fig. 8). The next highest significance levels between the species were for vertical body dimensions (BD, CPD and CPAD), which also exhibited considerable plot separation from one another (Fig. 5) (t values of ca. 10 for slopes between L. japonicus and L. maculatus and between L. japonicus and L. latus and ca. 27-37 for intercepts between L. maculatus and L. latus) ( Table 6).
Although the Mann-Whitney U tests for pairwise comparisons of meristic characters of the three species found significant differences in many, significance was not apparent for others, including counts of vertical fin rays [dorsal fin spines (DFSs), DFRs and AFRs] between L. japonicus and L. maculatus, and vertebrae [abdominal vertebrae (AVe), CVe and TVe] between L. japonicus and L. latus (Table 7).   Standard errors (SEs) for regression lines between logarithm-transformed SL and length-measured characters, and between SL and meristic characters are summarized in Table 8. For many characters, L. latus had the lowest SE values among the three species, followed by L. japonicus (Table 8). In general, degrees of SE could be ranked: 1 L. maculatus, 2 L. japonicus, 3 L. latus.
On the other hand, taxonomic and related literature on Lateolabrax have commonly noted the diagnostic importance of ranges and / or averages of body proportions (e.g., Katayama 1960a, b;Yokogawa and Seki 1995;Kim and Jun 1997;Yamada et al. 2007;Murase et al. 2012), although such, being commonly subject to allometric growth, are largely biased by the body sizes of specimens examined. For example, Figure 14 summarizes proportional body depth (BD) and orbital diameter (OD) ranges previously reported for L. japonicus and L. maculatus, respectively, compared with the present study. The smaller proportional ranges previously reported were all less than those presented here, representing many variously-sized specimens, suggesting that the former were based on relatively few specimens. Also, the variations in published proportional ranges, in some cases showing no range overlap (e.g., Fig. 14J vs K; L vs M), suggested differing body size ranges of the material studied. Although such proportional data has often been included in taxonomic diagnoses, the inherent inconsistencies have made specimen comparisons and specific identifications problematic. In fact, the use of proportions subject to isometric growth in species diagnoses is a legitimate procedure, although such proportions are rare in both Lateolabrax species (Table 3) and the other species listed above. However, the use of non-isometric proportional data, traditionally under the premise of (presumed) isometric growth, in species diagnoses is inappropriate.
Differing growth-related proportional change patterns in the three Lateolabrax species include pre-anus length (PAL) (Fig. 4A-C, Table 3) and post-orbital preopercular width (POPW) (Fig. 6K, Table 3). Similarly, the very similar East Asian frog flounders Pleuronichthys lighti and P. cornutus have the caudal fin, and dorsal and anal fins shortened with growth in the former and latter, respectively (Yokogawa 2015), indicating the potential for differing specific patterns, even between closely related species. Comparisons of black bass congeners (genus Micropterus) have shown the upper jaw length proportion to increase with growth in M. salmoides (Yokogawa 2014), while remaining stable in M. dolomieu (Senou 2002). Although the three Lateolabrax species share a similar "bass shape" with M. salmoides, the upper and lower jaw length (UJL and LJL) / standard length (SL) proportions decreased with growth in the former (Fig. 6M, O, Table 3), unlike the latter (Yokogawa 2014). Also, it is notable that BD and head length (HL) proportions of the three Lateolabrax species decreased with growth ( Fig. 5A, Table  3), in contrast to the centrarchids M. salmoides and L. macrochirus (Yokogawa 2013a), in which BD and HL increased with growth (Yokogawa 2014). This suggests that some phylogenetic factors may be responsible for growth-related proportional change patterns.
As in many other fishes (Okiyama 1988), BD of L. japonicus increased relatively with growth during the larval stage (from 13-16 to 26-30% of SL) until ca. 25 mm SL, thereafter being "stable," according to Tanaka and Matsumiya (1982) and Tamura et al. (2013), although subsequently decreasing from ca. 30 to ca. 21% of SL (Fig. 5A). Similarly, HL of L. japonicus and L. latus increased relatively with growth during the larval stage (Kinoshita 1988), in contrast to the growth-related acute decrement of HL during the juvenile and adult stages (Fig. 4S, U). During the larval stage of L. japonicus and L. latus (11-19 mm SL), the greater HL / SL proportion of the latter compared with the former in same-sized larvae, enabled ready distinction of the two species from each other (Kinoshita 1988). Although a similar distinction was observed in juvenile fishes (ca. 40-100 mm SL), very similar growth-related HL decreasing patterns between the two species in the adult stage (Fig. 4S, U) made it clear that Kinoshita's (1988) criterion for separation was applicable only for larvae of the two species.
Growth-related proportional change patterns of length-measured cephalic characters (based on SL and HL) were sometimes inconsistent in L. japonicus and L. latus (Fig. 6, Table 3), possibly due to HL being negatively allometric with SL (decreasing with growth) (Fig. 4S, U, Table 3) and paralleling or exceeding the change rate of some cephalic characters, resulting in negative allometry and isometry in SL-based relationships appearing as isometry and positive allometry in the HL-based ones, respectively. However, OD was negatively allometric relative to both SL and HL (Fig. 6C, D, Table  3), due to their degree of allometry relative to SL exceeding that of HL to SL. On the other hand, the consistency of the growth patterns between the two-way relationships in L. maculatus ( Fig. 6A-P, Table 3) may be due to the growth-related decreasing rate of proportional HL being less apparent in this species (Fig. 4T) than in the others (Fig.  4S, U) and therefore less influential on the relative growth of the cephalic characters. Although HL-based proportions of cephalic characters have been frequently used for cephalic characters in taxonomic studies on Lateolabrax (e.g., literature cited in Fig.  14), it should be recognized that the base dimension (HL) is not a stable character.
The proportional values (percentages) of proportions subject to allometric growth are correlated with the base dimension (e.g., SL and HL). In Figure 14, because both BD and OD were negatively allometric in both L. japonicus and L. maculatus (Figs 5A, 6D, Table 3), high and low proportional values are regarded as representing small and large size specimens, respectively. McClelland (1844) noted in the original description of L. maculatus (as Holocentrum maculatum) that the eyes were large, indicating that his description was based on a small specimen(s). The OD / SL proportion taken from his specimen illustration (pl. 21, fig. 1) was 6.4%, whereas the SL calculated by the inverse function of the SL-OD / SL regression (Fig. 6C, Table 2) was ca. 184 mm, agreeing with the above suggestion. This suggests that length-measured characters (including OD) subject to allometric growth can be utilized for estimation of body size. Hirota et al. (1999) compared their morphometric data for L. maculatus (as Lateolabrax sp.) from Kanto region, Japan [n = 6, 151-451 (average 298.3) mm SL] with those examined by Yokogawa and Seki (1995) (Hirota et al. 1999, table 1). Such inconsistency was clearly due to body size differences of the specimens examined in the two studies, i.e., the larger specimens in the former study provided lower OD proportions (Fig.  6D). Nevertheless, Hirota et al. (1999) suggested that the different OD proportions resulted from Yokogawa and Seki (1995) having measured eye diameter rather than OD, which was incorrect. Kim and Jun (1997) examined the morphology of Korean  OD (19.7 and 19.8) for the respective lots (Kim and Jun 1997, table 1). However, those degrees of proportional similarity between such different-sized specimens is extremely unlikely due to the highly negatively allometric proportions of those characters in this species (Figs 5A, B, 4S, 6C, Table 3).
Because most of the length-measured characters of the three Lateolabrax species were subject to allometric growth (Table 3), raw dimension measurement data were logarithm-transformed in order to transform the data distribution to be symmetric for statistical analysis, including canonical discriminant analysis (Bae et al. 2016) and analysis of covariance (ANCOVA), performed in the present study. Although Wang et al. (2016) provided multiple-regression analyses between body weight (BW) and some body dimensions for L. maculatus using raw data, the approach was problematic, because the raw dimension data (including BW) needed to have been logarithm-transformed before analysis, as done for M. salmoides by Yokogawa (2014).
Counts of pored scales on the lateral line (LLSs) and scales above the lateral line (SALs) tended to increase and decrease with growth, respectively, in L. japonicus (Fig.  9A, D, Table 4), those of scales below the lateral line (SBLs) and lower-limb and total gill rakers (LGRs and TGRs) tending to increase with growth in L. maculatus (Fig. 9H, K, N, Table 4). By contrast, overall meristic counts (except dots) did not change with growth in L. latus (Table 4), implying some phylogenetic determination of growthrelated meristic characters, as in the case of PSAL change patterns. Although the mechanism by which such counts increase or decrease with growth is uncertain, an SBL count increase with growth has been reported for L. macrochirus (Yokogawa 2013a), M. salmoides (Yokogawa 2014) and P. cornutus, in which gill raker numbers also increased with growth (Yokogawa 2015), suggesting that such phenomena are not so rare in fishes. Although meristic characters have been frequently used as important keys in taxonomic studies on the premise that they are stable at any body size, the potential for growth-related changes should be considered and actively assessed in taxonomic studies. Nozaka (1995) examined the morphology of L. japonicus fingerlings from the eastern Seto Inland Sea (n = 112, average 141.1 mm SL), comparing his data with Yokogawa and Seki (1995) [n = 65, 122.8-417.0 (average 301.4) mm SL] and noting differences in LLS and gill raker numbers (average LLSs = 76.4 and 83.1, average TGRs = 24.9 and 27.2, in the former and latter, respectively). Inconsistency in LLS counts may have resulted from body size differences in specimens examined, larger specimens resulting in higher LLS counts (Fig. 9A, Table 4). On the other hand, the difference in gill raker counts, which do not change with growth in L. japonicus (Fig.  9J, M, Table 4), may have resulted from the non-inclusion of rudiments located on the gill arch edges, since low gill raker counts as reported by Nozaka (1995) have not been found in the many other L. japonicus samples examined from around Japan (Yokogawa, unpublished data).
The growth-related status of dots / spots on the lateral body region also varied among the three Lateolabrax species. In L. japonicus and L. latus, although dots appeared in some smaller specimens (up to 260.6 and 254.8 mm SL, respectively), they disappeared with growth (Fig. 11A, C), a well-known phenomenon in the former species (e.g., Katayama 1960aKatayama , 1960bYokogawa 1995;Kim and Jun 1997;Kim et al. 2004;Ishikawa and Senou 2010), but barely noted in taxonomic descriptions of the latter species, other than Katayama (1957Katayama ( , 1960a and Murase et al. (2012). This may have been due to such dots being so fine or faint (Fig. 2C) that they were overlooked, or because descriptions were based only on large specimens. However, spot counts were not related to body size in L. maculatus, which typically had many clear spots in both large and small specimens (Fig. 2B, Table 4). Although many taxonomic descriptions of this species have incorrectly noted that spot counts decreased gradually with growth (Tchang et al. 1955;Chu et al. 1962Chu et al. , 1963Chu 1985;Chen 1987;Chen et al. 1990;Li and Zhang 1991;Feng and Jiang 1998;Chen and Fang 1999;Wang et al. 2001;Zhao and Zhong 2006), such may have been based only on subjective observations without statistical analysis, unlike the present study. On the other hand, large individuals of this species tend to have smaller and more rounded (non-jagged) spots than in small individuals [e.g., Katayama 1984, plate 108-I, 52 cm, as a variation of L. japonicus;   fig. 1, 600 mm in total length (TL), as L. sp.], which may have provided some grounds for the above views. Descriptions of L. maculatus (as L. japonicus) from Hong Kong noted that in young specimens, spots were larger and fewer in number, whereas with advancing fish age the large spots become smaller and more numerous (Chan and Tang 1968;Sadovy and Cornish 2000). However, although growth-related spot size decrement is correct, growth-related spot number increment is not.
The proportional growth-related change pattern of pectoral scaly area length (PSAL) in L. latus closely fitted a power regression (Fig. 4X, Table 2). However, simple patterned regressions could not be applied to L. japonicus and L. maculatus since they exhibited inverted V-shaped changes (Fig. 4V, W). This may reflect the phylogenetic status of the three species, L. latus being genetically further from the other two species (Yokogawa 1998;Song et al. 2008;Shan et al. 2016). A similar growth-related change pattern was also observed for the maximum blotch diameter on the dorsal fin (% of SL) in Banjos banjos banjos (Matsunuma and Motomura 2017, fig. 8d), inferring that such non-linear patterns arise in some characters in which dimensions are not determined by internal bony structure, rather than in normal body portions. Although PSAL, as defined by Yokogawa and Seki (1995) (see above), was examined in L. japonicus and L. maculatus, overall growth-related change patterns were limitedly revealed for both at that time due to size-biased samples. Accordingly, Nozaka's (1995) examination of L. japonicus fingerlings (see above) resulted in a much smaller proportional PSAL range and average than those given by Yokogawa and Seki (1995) for larger examples of that species. Such disagreement was regarded as arising from body size differences in the material specimens between the two studies.

Inter-specific differences and taxonomy
Lateolabrax latus is typically characterized by a deeper body, represented by BD and CPD. However, neither character provides unequivocal identification due to the range overlap for proportional BD and CPD between L. latus and L. japonicus (Katayama 1957(Katayama , 1965. In the present study, although the scatter plots for proportional BD and CPD of L. latus were well separated from those of the other two species, some overlap occurred in the smaller size class (< ca. 200 mm SL) (Fig. 5A, B). However, the newly defined dimension caudal peduncle anterior depth (CPAD), located between BD and CPD (Fig. 3), is suitable for distinguishing L. latus from the other two species, there being no plot overlap with the latter (border level 15%) (Fig. 5C).
The CPAD proportion may be a useful feature for specific identification, since it can also be determined from illustrations and photographs of Lateolabrax species. For instance, an illustration of "L. japonicus (as Perca-labrax japonicus)" in Fauna Japonica (Temminck and Schlegel 1846, pl. II, fig. 1, drawn by Keiga Kawahara) may, in fact, be L. latus, because the proportional CPAD (% of SL) measured from the illustration was 15.4%, falling within the range of the latter (Fig. 5C). Because the SL of the illustrated specimen estimated by the earlier-described procedure (use of an inverse function of SL-OD / SL regression for L. latus) was ca. 336 mm, proportional BD and CPD (% of SL), which had no plot overlap with the larger size classes (>200 mm SL) of the other two Lateolabrax species, may also be used for specific identification. The proportional BD and CPD of the illustrated specimen were 29.3 and 12.1%, respectively, corresponding with the ranges of L. latus (Fig. 5A, B). Although Katayama (1960b) also recognized the greater BD and CPD proportions of Temminck and Schlegel's (1846) specimen, he identified it as L. japonicus because the dorsal and anal fin ray counts (xiv, 13 and iii, 8, respectively) corresponded to the ranges for L. japonicus. In fact, he may have counted 12 spines in the first dorsal fin, and 2 spines plus 13 rays in the second (SDF). However, the SDF should be regarded as comprising 1 spine and 14 rays, the ray next to the first SDF spine having a distal branch. Specimens examined in the present study included 6 L. latus with 14 dorsal fin rays (DFRs) and 8 anal fin rays (AFRs), supporting the opinion that Temminck and Schlegel's (1846) illustration was of L. latus. Similarly, illustrations of L. japonicus and L. latus in Katayama (1965, figs 520 and 521) should actually be reversed, since their proportional CPAD (% of SL) values were 15.1 and 13.5%, respectively, falling within the respective ranges of L. latus and L. japonicus.
In addition to caudal peduncle stoutness in L. latus, Hatooka (2000Hatooka ( , 2013 proposed peduncle shortness as a diagnostic character of the species. Similarly, Murase et al. (2012) recorded proportions of caudal peduncle length (CPL) (% of SL) for L. latus (n = 27, 18.7-20.9), L. japonicus (n = 25, 20.0-23.4), and L. maculatus (n = 7, 20.7-22.3), indicating a clear difference between L. latus and the other two species. However, despite the distinctly downward shift in plot distribution in L. latus from the other two species found here, the CPL proportion range (n = 136, 18.3-22.7) largely overlapped those of L. japonicus (n = 229, 18.5-24.6) and L. maculatus (n = 170, 18.6-25.3), owing to considerable variation in plot distribution in the latter two species (Fig.  4G-I). The disagreement between the above two studies and the present one is likely to have resulted from differing numbers of specimens examined. In conclusion, although the proportional CPL of L. latus tended to be lower than in the other species, adoption of the feature as a diagnostic key for L. latus is problematic.
Caudal fin notch depth (CFND) has been recently proposed as a new character for distinguishing L. latus from the other two species, the former having a shallower CFND than the others (Hatooka 2000(Hatooka , 2013. However, although growth-related patterns of proportional CFND (% of SL) differed from one another among the three species (Fig.  4J-L) and ANCOVA for the logarithm-transformed regressions indicated significant differences of the slopes between any two species (Table 6), the ranges relative to overall SL (2.9-7.9, 2.0-8.4 and 1.9-7.4% for L. latus, L. japonicus and L. maculatus, respectively) were similar (Fig. 4J-L) and unable to distinguish between species. In fact, the proportional CFND of L. latus decreased acutely with growth, with relatively little variation owing to high correlation with SL (Fig. 4L, Table 3), being almost stable at low values (around 4-5%) in specimens > ca. 200 mm SL (Fig. 4L). In contrast, the other two species had highly variable proportional CFND, up to ca. 8% at any body size (Fig.  4J, K). Therefore, individual specimens of L. japonicus and L. maculatus with greater CFND may give the impression that L. latus has a shallower CFND than the others, as emphasized by some photographs of L. latus in which the caudal fins are so well opened that CFND decreases considerably (nearly truncate) (e.g., Masuda et al. 1975, pl. 42E;Ishikawa and Senou 2010;Murase et al. 2012, fig. 2C). It is possible that the caudal fin of L. latus may spread more than that of the other two species owing to broader membrane between the fin rays (Fig. 1E, F), particularly when fresh (when specimens were photographed). Notwithstanding, the results herein clearly indicate that CFND is problematic as a key character. Although Shimose et al. (2011) made underwater observations of and photographed a single Lateolabrax fish at Ishigaki Island, Okinawa, Japan, suggesting it to likely be L. latus based on some visually-recognized features, including CFND, the influence of such a key in the popular media is unfortunate.
Among the length-measured cephalic characters of L. latus, plot separation of that species from the others was marked for snout length (SNL) (Fig. 6A, B), post-orbital length (POL) (Fig. 6K, L), and upper and lower jaw lengths (UJL and LJL) ( Fig.  6M-P). In particular, SNL may be a practical means of distinguishing L. latus from the others because plots were vertically separated for both in the SL-and HL-based relationships (border levels ca. 9 and 28%, respectively) (Fig. 6A, B), which were similar to Murase et al.'s (2012) results. However, POL may not be practical for identification because the plots and vertical axis ranges overlapped considerably with those of L. japonicus (Fig. 6K, L). Although Murase et al. (2012) showed an unequivocal difference in POL (% of SL) between L. latus (n = 27, 14.1-15.8) and the other species [L. japonicus (n = 25, 16.1-18.5), L. maculatus (n = 7, 16.4-20.2)], such may have been due to the low numbers specimens examined, as in the case of CPL. The fact that SNL and POL of L. latus are greater and shorter, respectively, than in the other species infers that the eyes of L. latus are located more posteriorly than in the latter.
The UJL and LJL plots for all three species (SL-based relationships) were well clustered around their regression curves (high negative allometry), but could not be distinguished from one another vertically (Fig. 6M, O). On the other hand, since the UJL and LJL plots of L. latus in the HL-based relationships formed almost horizontal clusters, they could be vertically distinguished from those of the other two species (border levels of ca. 45 and 49%, respectively) ( Fig. 6N, P). Despite Murase et al.'s (2012) proposal of some diagnostic characters for L. latus including greater SNL and shorter POL, they excluded UJL, despite having measured that dimension. Although Hirota et al.'s (1999) (see above) examination of L. maculatus recorded SNL and UJL proportions (% of HL) [23.2-30.0 (average 26.3) and 39.4-46.4 (average 42.5), respectively], the maximum values of both fell within the ranges peculiar to L. latus (Fig. 6B, N). Assuming correct calculations, their catalogued "L. maculatus" specimens (whereabouts unknown) may have included L. latus. This possibility is also suggested by their higher counts of DFRs [13-14 (average 13.3)] and AFRs [8-9 (average 8.2)], including a small proportion of specimens (n = 6) with minor counts in L. maculatus [14 DFRs (16.6%) and 9 AFRs (5.3%)] (Fig. 10B, D).
The original description of L. latus included several diagnostic meristic characters, including counts of DFRs, AFRs and SBLs (Katayama 1957). In particular, DFR numbers =15, considered peculiar to the species, have subsequently been noted as an important diagnostic key (Katayama 1960a(Katayama , 1965(Katayama , 1984Masuda et al. 1975;Araga 1981;Hatooka 1993). However, because some L. latus specimens with 14 DFRs (overlapping the ranges of the other two Lateolabrax species) have been recognized (Sakai et al. 1998;Hatooka 2000Hatooka , 2013Murase et al. 2012), including 7.4% of L. latus specimens in the present study, DFR counts alone cannot absolutely distinguish L. latus from the others, although higher DFR counts may be useful (Fig. 10B). In contrast, AFR and SBL counts have rarely been adopted as diagnostic for L. latus, inferring that the count range overlaps between L. latus and the other two species are problematic for specific identification. In the present study, L. latus was well separated from the other species by AFRs (Fig. 10D) and DFRs, whereas SBL counts broadly overlapped (Fig. 10H). On the other hand, pectoral fin ray (P 1 FR) counts, which have not been emphasized as having taxonomic significance for L. latus, showed a strong modal shift between L. latus and L. japonicus (16 and 17, respectively) (Fig. 10E). Although the large range overlap of P 1 FR counts in L. japonicus and L. maculatus preclude their diagnostic use, they may be useful in the case of L. latus. For example, the two Lateolabrax specimens collected from Tanegashima Island both having 16 P 1 FRs (Sakai et al. 1998) are likely referable to L. latus.
In addition to length-measured and meristic characters in the original description of L. latus a further diagnostic feature proposed was the possession of ventromandibular scale rows (VSRs) (Katayama 1957). Although frequently noted as diagnostic for L. latus until recent years (e.g., Katayama 1960aKatayama , 1965Katayama , 1984Masuda et al. 1975;Araga 1981;Hatooka 1993), the possession of such scales has subsequently been omitted from keys to the genus Lateolabrax (Hatooka 2000(Hatooka , 2013 owing to the presence of VSRs in some specimens of L. japonicus and L. maculatus (Table 5) (Paxton and Hoese 1985;Hirota et al. 1999;Kang 2000;Murase et al. 2012). Furthermore, the lack of VSRs in some small L. latus (mainly ≤100 mm SL) (Table 5) underlines the unsuitability of this feature as a diagnostic character for L. latus. It was clear in the present study that VSRs did not exist in larvae and juveniles of all Lateolabrax species, but first appeared in L. latus at ca. 90 mm SL, thereafter rapidly developing with growth until present in almost all large individuals. In L. japonicus and L. maculatus, the appearance of VSRs was delayed, beginning from around 150 and 200 mm SL, respectively, and thereafter gradually developing with growth, although still absent in some large individuals. Such specific differences in squamation development may be common for PSAL (Fig. 4X) and dorsocephalic scale rows (DSRs) (Fig. 12), development being greatest in L. latus and least in L. maculatus, as indicated by Murase et al. (2012).
The diagnosis accompanying the original description of L. latus included ventral (pelvic fins) generally dusky, unlike in L. japonicus (Katayama 1957), followed by Katayama (1965) and Araga (1981). Although such coloring was infrequent in preserved L. latus specimens examined here, it has been noted in some large fresh adult specimens [e.g., photographs in Araga (1981) and Ishikawa and Senou (2010)]. However, non-dusky (pale) pelvic fins have been commonly observed in small L. latus (to fingerling size) (Fig. 1E, Murase et al. 2012, fig. 2A, B) and some large fresh condition specimens (Fig. 1F, Murase et al. 2012, fig. 2C). Possibly based on this supposed feature, the English name "blackfin sea bass" has been employed for L. latus (e.g., Matsuyama et al. 2002;Arakaki et al. 2014;FishBase 2018), however, such naming is not suitable, because it suggests that all fins were black, and many L. latus specimens including the large individual (915 mm TL) figured in FishBase (2018) do not have dusky ("black") pelvic fins. Instead, "flat sea bass," which describes the deeper body, a common feature of the species, should be applied for L. latus, following Yokogawa and Kishimoto (2012).
Recent keys for identification of L. japonicus and L. maculatus have adopted SNL, that of L. maculatus supposedly being relatively shorter than that of the former (Yamada et al. 2007;Hatooka 2000Hatooka , 2013. However, plots of proportional SNL largely overlapped in smaller size classes (< ca. 200 mm SL) of the two species, although plots for L. maculatus shifted downward (highly negative allometry) and were clearly separated from those of L. japonicus in specimens > ca. 200 mm SL (border levels ca. 7.7% and 24% for SL-and HL-based relationships, respectively) (Fig. 6A, B). Accordingly, SNL proportions enable separation only of large specimens (> ca. 200 mm SL) of the two species; e.g., Wakabayashi and Nakamura's (2003) L. maculatus specimen (as L. sp.) from Shima Peninsula, Japan (381 mm SL) was identifiable by its SNL proportions (7.1 and 22.8% of SL and HL, respectively).
On the other hand, post-orbital preopercular width (POPW) is a notable dimension, showing a contrasting pattern to SNL, i.e., plots of proportional POPW in small (< 200 mm SL) L. maculatus shifted upward and separated completely from those of similar sized L. japonicus (border levels ca. 7.5% and 23% for SL-and HLbased relationships, respectively), although larger specimens (> 200 mm SL) of the two species had some overlap due to the relative decrease of POPW with growth (highly negative allometry) in the former (Fig. 6I, J). Thus, a combination of SNL and POPW proportions [for small (< ca. 200 mm SL) and large (> ca. 200 mm SL) specimens, respectively] enables the two species to be separated unequivocally for their entire size range. Furthermore, the POPW / SNL proportion, which largely separates the two species throughout their entire size range (border level 90%) (Fig.  8), can also be adopted.
Proportional differences between L. japonicus and L. maculatus were also apparent in many of the fin lengths (first and second dorsal, caudal and pectoral), proportions of the former being distinctly greater than those of the latter in smaller specimens (< ca. 200 mm SL), although plots of the two species overlapped in the larger size class (> ca. 200 mm SL), due to the relative fin lengths decreasing and not changing with growth in the former and latter species, respectively (Fig. 5E-H). That this means of distinguishing between small specimens of L. japonicus and L. maculatus has largely gone unrecognized is probably due to a lack of morphological examination of small Lateolabrax specimens. The benchmark size of 200 mm SL being common to SNL, POPW and fin lengths of the two species suggests some synchronization of specific growth-related morphological changes.
Although Yokogawa and Seki (1995, figs 6, 7) proposed that considerable differences in LLS and gill raker numbers were sufficient for unequivocal differentiation of L. japonicus and L. maculatus when used in combination, the present study has demonstrated greater count range overlaps between the two species (70-84 LLSs and 24-26 TGRs, vs 76-82 LLSs and 24 TGRs) (Fig. 10F, K), due to LLS and gill raker counts increasing with growth in L. japonicus and L. maculatus, respectively (Fig. 9A, M, Table 4). Similarly, Kang's (2000) comparable frequency distributions of LLS and gill raker counts between the two species from Korean waters may have resulted from a size bias in specimens examined, his L. maculatus material including only very large specimens (ca. 500-750 mm SL). Accordingly, counts of LLSs and gill rakers, which can be biased by specimen size, are now likely to be unsuitable for distinguishing between the two species. In fact, Lou et al. (2002), who compared morphology between L. japonicus (1 sample lot from Tokyo, Japan) and L. maculatus (5 sample lots from Beihai, Xiamen, Fuzhou, Zhoushan and Weihai, China), showed considerable range overlaps for LLS and TGR counts, although the average values of those counts for L. maculatus were unequivocally lower than those for L. japonicus. Although Iseki et al. (2010) identified 263 Lateolabrax specimens from western Japan as L. maculatus (as L. sp.) based on LLS and gill raker counts proposed by Yokogawa and Seki (1995), some difficulties in identification may have been encountered due to some of their specimens being very large (up to 1130 mm SL), with gill raker counts that approached or overlapped the range for L. japonicus.
On the other hand, caudal and total vertebral counts (CV and TV, respectively), in which dominant counts were almost completely replaced between L. japonicus and L. maculatus (20 and 19 CVe, 36 and 35 TVe, for the former and latter, respectively) (Fig. 10M, N), may be useful for specific identification because they do not change with growth (Table 4). A modal count of 35 TVe in L. maculatus was indicated by Lou et al. (2002) (see above), who recorded average TV counts for 5 sample lots from China, viz., 34.75 (Beihai, n = 40), 34.64 (Xiamen, n = 19), 34.90 (Fuzhou, n = 10), 34.98 (Zhoushan, n = 27) and 35.07 (Weihai, n = 50), in spite of a geographic cline that suggested a trend towards lower and higher TVe in sample lots from southern and northern regions, respectively. Notwithstanding, Chen et al. (2001) recorded an average TV count of 35.31 (n = 98) for a sample lot from Laizhou, China, inferring that approximately 30% of their specimens had 36 TVe, which largely contradicts the present results (Fig. 10N). However, the former average count is suspect, differing considerably from the sample lot from Weihai (Lou et al. 2002), located close to Laizhou. In fact, such a high average TV value has not been recorded elsewhere at any time for L. maculatus (Yokogawa and Seki 1995;Lou et al. 2002). Although vertebral counts [abdominal (AV), CV and TV, respectively] of L. japonicus and L. latus are similar to each other, those of L. maculatus stand apart (Fig. 10L-N, Table 7), in contrast to their phylogenetic relationship (Yokogawa 1998;Song et al. 2008;Shan et al. 2016). In this case, since the difference in L. maculatus was primarily due to a difference in CV counts, which generally reflect inter-specific differences or lower, unlike AV counts which may reflect differences at a higher taxonomic level (Takahashi 1962), the vertebral count peculiarity in L. maculatus may not have phylogenetic significance.
Although L. maculatus typically possessed many black spots on the body, individual spot counts and patterns varied considerably (Yokogawa 2013b, fig. 2), a few specimens (4.9% of total) entirely lacking spots. In addition, the proportion of dotted L. japonicus specimens (35.6% of total) made visual separation of the two species difficult, the use of color pattern for specific identification being of value only as an accessory character. Youn's (2002) key, however, distinguished between the two species on the presence or absence of black spots, may causing mis-identification. Yokogawa and Seki (1995) demonstrated differences between L. maculatus and L. japonicus in some newly-demonstrated characters, including PSAL and DSRs (scale development in these characters being poorer in L. maculatus). However, because their examined material was size-biased (see above), overall growth-related change patterns were still unclear. Examination of PSAL and DSR in the present study have overcome that problem. Although differences between the two species were apparent in specimens < ca. 150 mm SL, squamation developed thereafter with growth in L. maculatus, the two species consequently having similar degrees of squamation in large specimens (Figs 7, 12). Notwithstanding, specific differences in specimens < ca. 150 mm SL can be used to identify Lateolabrax individuals up to fingerling size. Growth-related squamation development has been examined in laboratory-reared larval and juvenile L. japonicus (Fukuhara and Fushimi 1982) and L. maculatus (Kang 2000). Although squamation initially occurred on the caudal peduncle at ca. 19 mm SL in both species, body squamation was completed earlier in the former (ca. 35 mm SL vs 47 mm SL) (Fukuhara and Fushimi 1982;Kang 2000), indicating delayed development in L. maculatus. The slower development in PSAL and DSRs in L. maculatus might be an extension of such squamation delay, which is a characteristic peculiar to that species.
A morphological difference in the first anal pterygiophore (FAP) between L. japonicus and L. maculatus was initially noted by Kang (2000) during his detailed osteological observations of the three Lateolabrax species, and included in one of his keys (for adults) to the genus Lateolabrax; FAPs were arched and straight in L. japonicus and L. maculatus, respectively (Kang 2000). However, FAPs of small L. japonicus specimens (< ca. 90 mm SL) were found here to be straight (Fig. 13A), a condition not found by Kang (2000) due to his examining only larger specimens (minimum size 185.5 mm TL). Although Kang (2000) also described FAP in L. latus as straight, some examples of that species examined here had the FAP slightly arched distally (Fig. 13J, L). Because Kang (2000) examined only three L. latus specimens, ontogenetic morphological variations were not considered at that time. However, despite the growth-related morphological changes now apparent in L. japonicus, morphological differentiation of FAP is stable in specimens of L. japonicus and L. maculatus > 90 mm SL (Fig. 13B-D, F-H), enabling separation of the two species. Yokogawa and Kishimoto's (2012) identification of a long-finned Lateolabrax specimen from Japan (SPMN-h 40001, 331 mm SL) as L. japonicus was based on its genetic characteristics, although morphological identification of the specimen was equivocal, the TV count of 35 being suggestive of L. maculatus (Fig. 10N). However, identification of the specimen as L. japonicus was settled by the FAP being arched (Yokogawa and Kishimoto 2012, fig. 2a).
Standard errors (SEs) for the length-measured and meristic character regressions, which indicated degrees of morphological variation, were generally lowest in L. latus (Table 8), suggesting less morphological variation in that species. This may be due to less genetic variation, average observed heterozygosity for 28 isozymic loci in L. latus being 0.033, much lower than that of L. japonicus (0.095) and L. maculatus (0.103) (Yokogawa 1998). Usually, lower genetic diversity occurs in a small or reduced population, but the L. latus specimens examined in the present study were from a broad area around southern Japan. Possibly, in spite of the species' broad distribution, L. latus resources may not be so abundant, since the species is much less popular than L. japonicus in Japanese commercial markets. In contrast, SEs were generally highest in L. maculatus (Table 8), inferring considerable morphological variation. The significant geographical differences in otolith morphology among some L. maculatus samples from China (Ye et al. 2007) may have also resulted from its genetic diversity. This is supported by L. maculatus being broadly distributed along the east Asian coast, with some local populations being so genetically divergent from one another as to form a genetic / geographic cline, unlike L. japonicus, which is genetically stable (Yokogawa 2004;Liu et al. 2006;Han et al. 2015). In this regard, Zhao et al. (2018) reasonably considered that the Leizhou Peninsula, Hainan Island and Shandong Peninsula were major physical barriers, substantially blocking gene flow and genetic admixture among local L. maculatus populations.
The present study demonstrated a number of growth-related morphological changes in the three Lateolabrax species, including some new key characters for identification. Despite the number of taxonomic descriptions and studies of Lateolabrax, such features have remained obscure due to the limited numbers of specimens examined and an inherent belief that fish morphology is stable regardless of growth, notwithstanding some Endo, Mr. Hirokazu Kishimoto (Shizuoka City, Japan), Mr. Taiji Kurozumi (Natural History Museum and Institute, Chiba), Dr. Brian L. Sidlauskas (Oregon State University) and Mr. Ikuo Wakabayashi (Wildlife Research Society of Shima Peninsula) helped with provision of literature. Finally, I wish to thank Dr. Graham S. Hardy (Ngunguru, New Zealand) for checking the manuscript.