Geometric Morphological Differences Distinguish Populations of Scup in the Northwestern Atlantic Ocean

Abstract Scup Stenotomus chrysops, a commercially important marine species, are distributed from Nova Scotia to south Florida, and may represent morphologically distinct populations across their range. It was determined whether there were morphological differences between a North Atlantic Bight (NAB) population (41°N) and two South Atlantic Bight (SAB) populations (30°N and 32°N) of this species from April 2005 to July 2005, when the populations had formed spawning groups. Morphology was compared among populations by means of a geometric, landmark-based analysis of morphological and meristic traits for 180 individuals that were sexed and staged to maturity. A backward, stepwise discriminant functions analysis (DFA) produced a model that generated DFA scores that differed significantly between the NAB and SAB populations. Forehead and body depth dimensions differed among populations but not sexes. The DFA model predicted similar scores for 32 new NAB individuals and 31 new SAB individuals; however, scores were consistently underpredicted for the NAB individuals and overpredicted for the SAB individuals. We concluded that distinct northern and southern populations support a two-stock or two-species hypothesis; however, it is unclear whether this stock structure is related to reproductive or other life history traits. Either phenotypic plasticity or divergent selection may explain the morphological dissimilarities among these populations, but their influence on individual fitness remains unknown.

Scup Stenotomus chrysops, a commercially important marine species, are distributed from Nova Scotia to south Florida. In 2005, the U.S. commercial and recreational catch of scup was 5.3 million kg. This widely distributed species has historically been recognized as two or three distinct species (S. chrysops, S. aculeatus, and occasionally S. versicolor; Bigelow and Schroeder 1953). Early studies suggested that S. chrysops occurred north of Cape Hatteras, North Carolina, and S. aculeatus occurred south of Cape Hatteras (Johnson 1978), which may be a formidable barrier to gene flow between the North Atlantic Bight (NAB) and South Atlantic Bight (SAB) (Jones and Quattro 1999;Adams and Rosel 2006). Stenotomus aculeatus is not currently recognized as a distinct species (Robins et al. 1991;Carpenter 2002). Dispersal of larval fish between the NAB and SAB (Hare and Cowen 1996;Grothues et al. 2002) may prevent genetic isolation of these populations (Jones and Quattro 1999). However, no published studies have examined the level of morphological divergence between NAB and SAB populations. Our objective was to compare the size and shape of scup from the northern and southern extents of its range.
The science of morphometrics has a long history in quantifying phenotypes and has been radically changed by recent, widespread use of the digital processing of specimens (Cadrin and Friedland 1999). Phenotypic differences in size and shape within a species may exist because of sexual dimorphism (Love 2002) and ecological specialization (Sage and Selander 1975;Schluter and McPhail 1992). Size differences among fish within the same age-cohort may also reflect differences in reproductive capacity and growth, which are locally adaptive (Conover et al. 2005). Quantifying phenotypic differences among fish populations and identifying their causes and consequences may be informative of differences in natural history across a species' geographic range, which would have implications for both theoretical and applied work in ecological and fishery science.
This study sought to determine the extent to which morphological and meristic traits discriminate scup populations from the northern and southern extents of the species' range in the northwestern Atlantic Ocean. We tested the hypothesis that the size and shape of Subject editor: Richard Brill, Virginia Institute of Marine Science, USA similarly sized individuals differed between the populations from the NAB and the SAB.

Methods
Scup were collected from May to July 2005, when spawning aggregations were formed for the NAB (418N) and SAB (308N-348N) populations ( Figure 1). We compared populations from two SAB latitudes to evaluate the level of morphological difference across a relatively small spatial area. Adults were collected by means of a combination of trawling and trapping methods. Nine NAB specimens were collected during the Massachusetts Division of Marine Fisheries 2005 Spring Inshore Bottom Trawl Survey (1.8-cm-mesh otter trawl; 20-min tow at 2.5 knots during the day). Additional NAB specimens were collected by means of pot traps (60 cm wide 3 60 cm deep 3 60 cm high) that were soaked for 3 h. All SAB specimens were collected by means of wire fish traps (Collins 1990) set by the South Carolina Department of Natural Resources during surveys for their Marine Resources Monitoring Assessment and Prediction program. Collection depths ranged from 6 to 36 m. All individuals were frozen before processing. After thawing, each individual was measured, weighed, photographed for digital analysis (see below), and dissected to determine sex. Lateral line scales and other meristic traits were counted following Hubbs and Lager (2004). We counted the number of spines and rays for the dorsal, anal, and pelvic fins. Measurements for morphometric characters were made with a ruler (to the nearest 0.10 cm). These measurements included the length of the longest unbranched dorsal ray, the length of the pectoral fin, and the total length. Additional measurements of shape were recorded from digital images. Digital images were taken of thawed fish with a Nikon Cool Pix 5.0 Mega Pixels camera and analyzed with TPSdig (version 2.04; Rohlf 2005). Before digital imaging, fish were pinned to a dissecting pan and 10% solutions of formalin were applied to their fins for stiffening. For each digital image, 17 landmarks that could be clearly represented on the body of the fish were chosen out of 55 originally examined ( Figure 2; Table 1). A box truss of 28 lines connecting these landmarks was generated for each fish to represent the basic shape of the fish (Strauss and Bookstein 1982;Cadrin and Friedland 1999). Not all landmarks were used to generate a box truss in order to reduce the number of linear dimensions subjected to analysis and ease interpretation of the results. The effect of length on the linear dimensions was removed by means of linear regression. The residual variance was saved from individual regression models of length (independent variable) and the linear dimension (dependent variable). Kolmogorov-Smirnov tests for normality were performed for all variables. As there were significant departures from normality for all variables, the linear dimensions were log 10 transformed to help normalize the variances, but in no case was normality achieved. However, histograms indicated approximate normality for each variable, probably as a result of the large sample sizes.
Thirty-two response variables subjected to a backward, stepwise discriminant functions analysis (DFA), the P-value to leave being 0.15 or more. The variables included 28 linear measurements from the box truss network, mass, the length of the longest unbranched dorsal ray, the length of the pectoral fin, and the lateral line scale count. The multivariate DFA reduced the dimensionality of the data set by systematically removing response variables that contributed little to discriminating among latitudinal populations. The relative importance of each response variable in the classification of the independent variable was indicated by its discriminant coefficient standardized by a constant. The coefficients were used to generate a DFA model for assigning additional individuals to particular latitudes (see below).
The performance of the DFA was evaluated in three ways. First, scores from the DFA were saved after the analysis and plotted to interpret the level of separation among populations. If morphology differs between the NAB and SAB populations, the scores for individuals from these two regions should not overlap. The scores for individuals from both SAB latitudes were predicted to overlap. Second, Wilk's k was used to evaluate the significance of the discriminant functions model for discriminating among populations. Significant k values would indicate that discriminant function scores differ between groups. Third, the error rate for classifying individuals to each latitudinal population was evaluated by means of a matrix of classification errors computed from a DFA of a subset of the data (via jackknifing).
The variables important for discriminating among populations were determined by two methods. First, we used the magnitude of the standardized discriminant coefficient as an indication of the importance of each Distal tip of premaxillary bone, closed mouth 2 Nasal bone, bump over eye 3 Deepest part of frontal bone 4 Tip of supraoccipital bone 5 Origin of dorsal fin 6 Origin of unbranched rays 7 Termination of dorsal fin 8 Dorsal side of caudal peduncle, at the nadir 9 End of lateral line 10 Ventral side of caudal peduncle, at the nadir 11 Termination of anal fin 12 Origin of anal fin 13 Origin of pelvic fin 14 Dorsal origin of pectoral fin 15 Dorso-posterior tip of interopercle 16 Articulation of maxillary and dentary bones, closed mouth 17 Center of eye  338N and 11 from 348N). Following the procedures mentioned above, morphometric data were generated for new individuals collected from May to July 2005. Because new scores were calculated from the existing DFA model, the performance of the model could be assessed by (1) plotting the new predicted scores against the scores assigned originally (see above), (2) comparing the predicted DFA scores for the NAB and SAB by means of analysis of variance (ANOVA), and (3) conducting a nested ANOVA that tested for effects due to position (NAB versus SAB) and prediction error (original scores versus predicted scores) within positions. If the DFA model was effective, the plots of the new predicted scores should overlap the plot of the original scores. In addition, we hypothesized that there would be a significant difference in scores between the NAB and SAB populations but no significant nested effect (i.e., no difference between original and predicted scores for either the NAB or SAB). All analyses were performed in SYSTAT (version 10).

Results
We analyzed data for 180 specimens of scup (Table  2). There were 78 males, 87 females, and 15 of unknown sex. The number of males did not differ from the number of females between the NAB and SAB (Pearson v 2 ¼ 0.81; P ¼ 0.67) or at any latitude ( Table  2). The total length of specimens also did not differ among latitudes.
Of the 32 morphological variables subjected to the backward, stepwise DFA, 20 variables had relatively high discriminant coefficients. These variables were retained in the final DFA model (Table 3), which performed well in discriminating the NAB and SAB populations (Wilks' k ¼ 0.06; F 40, 312 ¼ 24.3; P , 0.0001; Figure 3). There was fairly high accuracy in predicting individual assignment to the northern populations (90%) and southern populations (83-84%; Table 4). Error rates were lower for classifying individuals to the NAB and SAB populations than for classifying them within the SAB, supporting a twospecies or two-stock hypothesis of stock structure.
Of the 20 variables in the final DFA model, 8 variables differed significantly among populations (MANOVA: Wilk's k ¼ 0.059, F 40, 308 ¼ 23.9, P , 0.0001; Table 2). Sex did not significantly explain the variation in shape (Wilk's k ¼ 0.812, F 40, 308 ¼ 0.884, P ¼ 0.737). Most significant variables were associated with body shape (n ¼ 2) and forehead shape (n ¼ 3). Forehead shape differences were related to the distance from a bump above the eye (point 2; see Figure 2) to other landmarks on the body. Foreheads for SAB individuals were typically less convex, the distance between the eye or the pelvic fin and the forehead being smaller than for NAB individuals. In addition to shape differences, the lateral line scale count was lower for an SAB population (308N) than for the NAB population (see Table 2).

Discussion
Scup populations differed in morphology between the NAB and SAB, supporting a two-species or twostock structure for this species. However, most morphological differences were related to forehead shape, which may not be related to important growth or reproductive characteristics of interest in stock assessments. For other species of the family Sparidae (e.g., squirefish Chrysophrys auratus), forehead shape naturally differs among populations and sexes (Moran et al. 1999). Interestingly, there was no evidence of sexual dimorphism in forehead shape for our scup specimens. Other morphological differences between the NAB and SAB populations, such as mass and scale count, were slight or nonsignificant.
Morphological divergence in shape among conspe- cific populations may be caused by a phenotypic plastic response to local environmental differences or be an adaption resulting from differences in resource use (Schluter and McPhail 1992;Snorrason et al. 1994;Langerhans et al. 2003). Larval fish may develop different morphologies because of local environmental conditions that vary by latitude, leading to differences in adult morphology as well. According to Cadrin and Silva (2005), the geographic variation in adult morphology for yellowtail flounder Limanda ferruginea may be explained by differences in ontogenetic rates among local populations if morphology is a product of ontogenetic history. Different prey resources for larvae and juveniles may lead to different adult body shapes, suggesting that there is a phenotypic plastic response to resource availability (Wimberger 1990). It is not clear whether there are differences in ontogenetic rates or resource availability for larvae and juveniles between the NAB and SAB populations. Cadrin and Silva (2005) ultimately determined that all TABLE 3.-Canonical coefficients standardized for within-unit variance, as derived from a discriminant functions analysis of morphological data (see Table 1) from three populations (308N, 328N, and 408N)    yellowtail flounder constitute a single stock despite the geographic variation in their morphology, largely because there is substantial gene flow among local populations. Adult morphology may also be determined by diversifying selection and ecological adaptation. Costa and Cataudella (2007) found that shape differences were related to trophic ecology for several species of the family Sparidae, thus indicating local adaptation (Schluter and McPhail 1992;Langerhans et al. 2003) and possibly ecological radiations. However, there are no a priori reasons to expect that forehead shape is adaptive for scup. To discriminate between phenotypic plasticity and ecological adaptation as mechanisms explaining morphological divergence, a combination of laboratory studies and genetic tests using population markers (e.g., the D-loop and cytochrome b) could be useful. Laboratory studies designed to test the plastic response of larvae to diet or temperature conditions may be informative, particularly if shape differences are detectable at earlier ages than those examined here. Genetic tests may distinguish discrete population units, especially if rapidly mutating (e.g., D-loop) and relatively slowly mutating (e.g., cytochrome b) markers are used in concert to distinguish between populationlevel and species-level differences, respectively. Morphological differences, whether related to ontogenetic growth history or ecological adaptation, may be preserved under instances of low gene flow, such as that observed for metapopulations (Cadrin and Silva 2005), or high gene flow (Anderson et al. 2008); no work has addressed the gene flow between the NAB and SAB scup populations.
Our research represents an important advance in understanding scup biology and highlights the morphological divergence in this widely distributed marine species. Such divergence is commonly observed for freshwater and lacustrine populations (Schluter and McPhail 1992;Snorrason et al. 1994). We have identified significant differences in the shape of scup from populations at extremes of its range in the northwestern Atlantic Ocean. If shape is related to either environmental influences on larval development (Cadrin and Silva 2005) or diversifying selection and ecological adaptation at a trophic level (Costa and Cataudella 2007), then spatially or latitudinally different environmental factors (e.g., temperature and resource availability) may explain the variations in body shape among scup. Langerhans et al. (2003) found that the distance between habitats correlated positively with the level of divergence in body shape among conspecific populations of two neotropical fish species. Although our data support a two-stock or twospecies model for scup, the morphological variation that we observed may be the result of an environmental gradient rather than genetic isolation or differences in the age structure and reproductive fitness of adults. Additional data on life history differences and gene flow between the NAB and SAB populations will better elucidate the management units required for this species and possibly those for other wide-ranging species in the northwestern Atlantic Ocean.