Bipartite life cycle of coral reef fishes promotes increasing shape disparity of the head skeleton during ontogeny: an example from damselfishes (Pomacentridae)
© Frédérich and Vandewalle; licensee BioMed Central Ltd. 2011
Received: 16 November 2010
Accepted: 30 March 2011
Published: 30 March 2011
Quantitative studies of the variation of disparity during ontogeny exhibited by the radiation of coral reef fishes are lacking. Such studies dealing with the variation of disparity, i.e. the diversity of organic form, over ontogeny could be a first step in detecting evolutionary mechanisms in these fishes. The damselfishes (Pomacentridae) have a bipartite life-cycle, as do the majority of demersal coral reef fishes. During their pelagic dispersion phase, all larvae feed on planktonic prey. On the other hand, juveniles and adults associated with the coral reef environment show a higher diversity of diets. Using geometric morphometrics, we study the ontogenetic dynamic of shape disparity of different head skeletal units (neurocranium, suspensorium and opercle, mandible and premaxilla) in this fish family. We expected that larvae of different species might be relatively similar in shapes. Alternatively, specialization may become notable even in the juvenile and adult phase.
The disparity levels increase significantly throughout ontogeny for each skeletal unit. At settlement, all larval shapes are already species-specific. Damselfishes show high levels of ontogenetic allometry during their post-settlement growth. The divergence of allometric patterns largely explains the changes in patterns and levels of shape disparity over ontogeny. The rate of shape change and the length of ontogenetic trajectories seem to be less variable among species. We also show that the high levels of shape disparity at the adult stage are correlated to a higher level of ecological and functional diversity in this stage.
Diversification throughout ontogeny of damselfishes results from the interaction among several developmental novelties enhancing disparity. The bipartite life-cycle of damselfishes exemplifies a case where the variation of environmental factors, i.e. the transition from the more homogeneous oceanic environment to the coral reef offering a wide range of feeding habits, promotes increasing shape disparity of the head skeleton over the ontogeny of fishes.
A primary aim of evolutionary biology is to explain the origin, structure and temporal patterns of phenotypic diversity. Many studies have proposed adaptive explanation for phenotypic variation and have focused on the role of selection in shaping patterns of diversification. Divergent selection is expected to enhance adaptive differences through time (evolutionary and ontogenetic scales). Selection may occur at any life stage and be specific to one life stage, leading to the potential for stage-specific adaptation . Despite its obvious importance, very little is known about the role of ontogeny in adaptive divergence and rigorous studies focusing on the variation of phenotypic diversity over ontogeny in various zoological groups are needed.
Coral reef fishes represent one of the most diverse assemblages of vertebrates, and moreover studies dealing with their morphological diversity are numerous [2–10]. These studies have mainly focused on the adult stages and fewer have addressed ecomorphological variation through the ontogeny of reef fishes [e.g. [11, 12]]. The majority of coral reef fishes have a complex life-cycle with two distinct phases: (1) a dispersive pelagic larval phase and (2) a sedentary demersal adult phase associated with the coral reef environment. The larval phase ends at reef settlement . In contrast to reefs, the open water environment is thought be more homogenous, especially with regard to the diversity of habitats. Conversely, intrinsic factors of the coral reefs such as high productivity, high spatial and ecological complexity, and high trophic diversity may be involved in promoting the high standing levels of fish diversity . It is thus surprising that quantitative studies on the ontogeny of the radiation of coral reef fishes are lacking. Such studies dealing with the variation of disparity, i.e. the diversity of organic form , over the ontogeny could be a first step in detection of evolutionary mechanisms of diversification in these fishes. By studying their ontogeny, via changes in shape with size (i.e. allometry), it is possible to gain a clearer understanding of the timing of selective pressures in coral reef fishes .
To our knowledge, very little quantitative data exists on adaptive allometry and on the ontogeny of shape diversification within zoological groups in general and fish clades especially [but see [1, 16, 17]]. Zelditch et al.  have focused on the dynamic of body shape disparity over ontogeny in piranhas and highlighted that the disparity decreases significantly and substantially over ontogeny within this particular clade. Adams and Nistri  investigated ontogenetic trajectories of foot morphology in eight species of European plethodontid cave salamander and showed the disparity of adult foot morphology was significantly lower than in juveniles. Two broad categories of factors have been suggested for explaining these ontogenetic patterns of disparity : the external or "ecological" constraints such as the availability of ecological space [21, 22], and the internal constraints such as developmental, genetic or functional constraints [18, 23, 24]. New studies in various taxa selected according to their life-cycle and their ecological diversity are needed. They should give further insights into the roles of external and internal constraints shaping the levels and patterns of disparity over ontogeny.
Including more than 300 species living in coral reef environments, the damselfishes (Pomacentridae) represent one of the most successful radiations of coral reef fishes [3, 14, 25, 26]. Similar to the great majority of coral reef fishes, they have a bipartite life-cycle. In the pelagic environment, all damselfishes larvae feed on planktonic copepods . At the adult stage, three trophic guilds are commonly recognized [25, 28]: the pelagic feeders sucking planktonic copepods; the benthic feeders grazing filamentous algae or biting coral polyps; and an intermediate group feeding on planktonic prey, small benthic invertebrates and algae in variable proportions. Consequently, we can consider that damselfishes have a higher trophic diversity at their adult stage, and thus we can expect that functional demands lead to an increase in morphological disparity over ontogeny. Indeed, while all species might be expected to be adapted to catch small planktonic prey during the pelagic larval phase, adult damselfishes would be expected to have specific morphological adaptations allowing optimal prey catching and processing. For example, planktivorous and algivorous species have muscles and skeletal shapes optimizing suction feeding and grazing, respectively [3, 4, 29].
The evolution of damselfishes exemplifies a case of reticulate adaptive radiation  in which morphological divergence at speciation has been associated with the repeated convergence on a limited number of ecotypes: algivory, omnivory and planktivory. Among the Pomacentridae, algivory and omnivory have both arisen seven times, planktivory four times, and feeding on scleractinian coral polyps twice [3, 25, 26].
The study of ontogenetic allometry can reveal differences in developmental patterns among species underlying their morphological differentiation. How differences in developmental patterns generate species divergence has been successfully addressed in many groups, including trilobites [e.g. ], fishes [e.g. ], newts [e.g. ], reptiles [e.g. ], rodents [e.g. ], primates [e.g. ] and humans [e.g. ]. The pattern of ontogenetic shape changes can be described by the allometric trajectory of an organism plus the rate at which it proceeds along the trajectory . Morphological divergence among taxa sharing a common allometric trajectory could be due to changes in the rate or the duration of development (i.e. rate and event heterochrony ). Divergence could result from the directional change in the allometric trajectories (i.e. allometric repatterning ), coupled or not to alterations in the rate and/or the timing of development. In the present study, we explore whether and how differences in allometric trajectories can shape the ontogenetic dynamic of morphological disparity in damselfishes.
We used geometric morphometrics to examine the patterns of morphological diversification among damselfishes throughout their ontogeny. Four functional units of the head skeleton were separately studied: the neurocranium, the unit «suspensorium and opercle», the mandible and the premaxilla. All these skeletal units are movable elements involved in prey catching . The aims of the present study were to: (1) test the hypothesis that larval morphologies are more similar (less disparate) than the morphologies of juveniles and adults, (2) compare the patterns of shape disparity, i.e. the distribution of shapes and the dimensions along which shapes are most disparate, at different ontogenetic stages, (3) compare ontogenetic trajectories to identify the evolutionary changes in developmental parameters (i.e. allometric patterns, rate of shape changes, amount of shape changes undergone over the course of ontogeny) shaping the level and the pattern of morphological disparity, and (4) explore the relationships between the developmental parameters, phylogenetic data and ecological data.
Sample and data collection
Size range (SL, mm)
We used Procrustes-based geometric morphometrics to study shape variation and dynamics of shape disparity throughout ontogeny [43–46]. For each skeletal unit, the digitized landmark configurations were subjected to a Generalized Procrustes Analysis (GPA) in order to remove non-shape variation (location, orientation and scale) [47, 48]. The "grand mean", i.e. the consensus shape of all specimens and shape variables were then generated as partial warp scores (PWs) including both uniform and non-uniform components [44, 49]. The centroid size (CS) of the structure was also computed as the square root of the sum of the squares of the distances from all LMs to their centroid . Allometry refers to the pattern of covariation between measures of size and shape . Age information is not available for our specimens, so we use size as ontogenetic scale (i.e. ontogenetic allometry).
1. Comparing models of ontogenetic allometries
The allometric patterns of shape variation were analyzed using linear multivariate regression of PWs on log-transformed size (ln-CS) [48, 50–52]. The null hypothesis that shape develops isometrically was tested in all species using TpsRegr (Version 1.34). The fit of the regression models was evaluated by the explained variance of the model and by a permutation test based on a Generalized Goodall's F-Test with 10,000 permutations.
Differences in allometric models among species were tested by a MANCOVA, testing the null hypothesis of homogeneity of linear allometric models. In these tests, shape variables (PWs) are considered as dependent variables, size (ln-CS) as covariate and species is grouping factor. As suggested by Zelditch et al.  and Webster and Zelditch , two factors can explain differences in allometric models: (a) the divergence of allometric trajectories and (b) the rate of shape changes. Consequently, we estimated and compared both factors plus two other parameters: shape at settlement (i.e. the starting point of the ontogenetic trajectories) and the length of the ontogenetic trajectory, which is a function of the rate and duration of development. As these parameters are estimated by linear multivariate regressions, they assume a linear relationship between shape and ln-CS. Plots of Procrustes distance from the mean larval shape (see below), and the variance explained by the models are used to assess the validity of log-linear models (Goodall's F-test). This is reasonable for our data because the regression of shape on ln-CS explains a large proportion (up to 88%) of shape variance for each structure in every species (see Additional File 1).
1.a. Comparing larval shapes at settlement (starting point of ontogenetic trajectories)
The samples of settling larvae are limited (≤ 10 specimens). Consequently, we used a standardized regression residual analysis to estimate and compare larval shapes at settlement [18, 54, 55]. From the multivariate regression of shape on ln-CS, the non-allometric residual fraction is standardized by Standard6 (IMP-software). «Standardized» data sets of larvae with their respective SL (see SLsettlement in Table 1), which are the predicted shapes of the entire population at these sizes, are generated. Multivariate analyses of variance (MANOVA) were performed using partial warp scores, followed by a posteriori tests of pairwise comparisons to test shape differences among species. The statistical significance of the pairwise differences was tested by a resampling-based F-test and the results of CVA assignment tests were also examined. MANOVAs were performed using Statistica 8.1 (Statsoft 2007). Pairwise F-tests were done in TwoGroup6 (IMP-software) and misclassification rates were given by CVAGen (IMP-software).
1.b. Comparing ontogenetic trajectories in shape space
The differences in trajectories of shape changes were analyzed by comparing the angle between the species-specific multivariate regression vectors using VecCompare6 (IMP-software). This test is described in detail in [48, 53] and was already exemplified in a previous study of allometry in damselfishes . Here, a within-species vector is composed of all regression coefficients of the shape variables (PWs) on the log-transformed CS. The range of angles between such vectors within each species is calculated using a bootstrapping procedure (N = 400). This range was then compared with the angle between the vectors of both species. If the between-species angle exceeds the 95% range of the bootstrapped within-species angles, the between-species angle is considered significantly different, and thus the allometric trajectories are different. Angles were computed pairwise between allometric vectors, and the resulting interspecific dissimilarity relationships (angles) between the allometric trajectories were summarized with scatter plots calculated using nonmetric multidimensional scaling (NMDS). Differences between a reference shape (i.e. larval shape at settlement) and a target shape (i.e. adult shape) can be illustrated with deformation grids (interpolating function "thin plate spline" or TPS ). Multivariate regression models of shape on size for each unit were used to show graphical illustrations of ontogenetic allometries in each species. Deformation grids illustrating these ontogenetic allometries were obtained from TpsRegr (Version 1.34).
1.c. Comparing rates of shape changes during ontogeny
The dynamics of shape change (developmental rate), defined as the rate of shape changes per unit of size in this study, was estimated for each species using the Procrustes distance (PD), the metric defining shape dissimilarity in the Kendall shape space . PD between each specimen and the average larvae were regressed on ln-CS in all species separately. The rate of divergence away from the average larval shape was compared among species using the slope of the regressions with Regress6 (IMP-software) ([for detailed explanations on this methodology, see [48, 56]). Because the relationships between PD and ln-CS are close to linear (see Additional File 2), these can be statistically compared by ANCOVA.
1.d. Comparing lengths of ontogenetic trajectories
The length of the ontogenetic trajectories is used to compare the net amount of shape change undergone over post-settlement ontogeny in the eight species. These lengths were calculated as the PD between the average larval shape at settlement and at the maximum adult body size. Confidence intervals are placed on these lengths by a bootstrapping procedure detailed in . Because our analyses are based on size-standardized data, so the bootstrapping procedure also takes the uncertainties of the regression into account. The calculations were done using DisparityBox (IMP-software).
We focused on the relationships between the studied developmental parameters, phylogenetic data and the diet of each species. Moreover, correlation analyses (Pearson r) were used to test the relationships between some developmental parameters, i.e. the dynamic of shape changes or the length of ontogenetic trajectories, and other species-specific characters, i.e. the pelagic larval duration (PLD) and the size variation undergone during post-settlement growth.
2. Measuring and analyzing shape disparity
We used the same methodology as  to measure morphological disparity (MD) at three ontogenetic stages: (1) at the mean size observed at settlement (Table 1), (2) at a common size of 60 mm SL for every species and (3) at the species-specific maximum adult body size (MAX SL). The intermediate size (60 mm SL; at such a size, some species are already adults and others are always at the juvenile stage) was arbitrarily chosen to estimate the MD when fishes are already settled on coral reef and when size difference among species are eliminated. For these analyses, we used "standardized" data sets, which are the predicted shapes of the entire population (see Table 1) at these stages.
where d j = PD between the mean shape of species j and the grand mean shape (i.e. consensus shape). N = total number of species. The level of disparity among the eight studied species was calculated for each skeletal unit (i.e. disparity at the family level). Then, we calculated the disparity level at the three studied ontogenetic stages within each trophic group. As only two species of the intermediate group were studied (i.e. D. aruanus and P. pavo), we could not calculate the disparity level of this trophic group. Consequently, both species were related to the two other trophic groups. As P. pavo commonly grazes algae, it was grouped with algivorous species. Conversely, D. aruanus was groupd with the zooplanktivrous species. All the calculations were done by DisparityBox (IMP-software), which also uses a bootstrapping procedure to place confidence intervals on this measure. We refer to [18, 48] for a detailed explanation of this methodology.
In addition to measuring the level of shape disparity, we also examined its pattern of variation, i.e. the distribution of shapes and the dimensions along which shapes are most disparate. This exploration is informative about the dynamic nature of disparity. For example, larvae and adults can have the same level of disparity but they can occupy two different sub-spaces (hyperplanes, i.e. "flat" surfaces of more than two dimensions embedded in higher dimensional space) of the shape space. Biologically, this would mean that the patterns of variation have been re-organized or re-structured to lie along very different pathways. The structures of disparity in different ontogenetic stages can be compared by comparing variance-covariance matrices. First, visual exploration of shape variations in the sub-space defined by the first two relative warps (Relative Warps analysis is equivalent to a Principal Components Analysis [PCA] of shape variables when the scaling factor α = 0, ) allow to check if the three ontogenetic stages (i.e. settling larvae, 60 mm SL, MAX SL) share the same hyperplane. Then, we used the program SpaceAngle (IMP-software) to compute the angle between hyperplanes and determine if the between-ontogenetic stages angle is no larger than the within-ontogenetic stages range. This program allows to specifically test if every ontogenetic stage occupied the same subspaces of the morphospace. SpaceAngle (IMP-software) used a method developed by  based on PCA. See  for detailed explanations on this approach but the angle between two subspaces embedded in a common higher dimensional space can be defined as the angle through which one subspace must be rotated to match the other. If the angle between two subspaces (i.e. between-ontogenetic stages angle) does not significantly differ from zero, we may not reject the null hypothesis (H0) in which specimens of two ontogenetic stages occupy the same subspace. The significance of the angle was determined by a bootstrapping procedure (N = 400) similar to that applied for comparisons of allometric vectors (see above). Relative warps analyses of standardized shape data at the three ontogenetic stages were performed using PCAGen (IMP-software) and the angles between subspaces defined by the first two PCs were calculated by SpaceAngle (IMP-software).
Geometric morphometric analyses were performed using computer programs from the TPS series (TpsDig and TpsRegr), written by F.J. Rohlf (freely available at: http://life.bio.sunysb.edu/morph/) and the IMP series (CVAGen, DisparityBox, PCAGen, Regress6, Standard6, TwoGroup6, VecCompare, SpaceAngle), created by H.D. Sheets (freely available at: http://www2.canisius.edu/~sheets/morphsoft.html). TPS deformation grids were generated in the program MORPHEUS (Slice, 1999; http://www.morphometrics.org/morpheus.html. STATISTICA, version 8.1 (Statsoft 2005) was used for other statistical analyses (i.e., MANCOVA and MANOVA) and NMDS plots were generated using Matlab (The MathWorks 2007).
Comparing ontogenetic allometries
The post-settlement ontogeny of the cephalic region is highly allometric in damselfishes. Shape variation in each skeletal unit is significantly correlated with log-transformed size (ln-CS) (all p levels of the Generalized Goodall's F test < 0.05, see Additional File 1) and a large proportion (up to 88%) of variation in total shape change during damselfish post-settlement development is explained by ontogenetic allometries.
Tests for common linear allometric models in the eight damselfish species
Suspensorium and opercle
Comparisons between ontogenetic trajectories
Suspensorium & opercle
Using deformation grids, differences in allometric patterns among species are often visually striking (see Additional File 3). However, at large-scale, some ontogenetic shape changes are common to all species. Adults had shorter and higher neurocranium than larvae while these differences are limited in the two Chrysiptera species and P. pavo. In all other species, the heightening is mainly explained by the growth of the supraoccipital crest (LMs 5-7). Adults showed higher opercle and suspensorium than larvae. The suspensorium of adults is always shorter in its central part. Indeed, the distance between the articulation of the palatin (LM 2) and the two articulations of the hyomandibular (LMs 4, 5) is always shorter in adults. All larvae had a maxillary process of the palatin (LMs 1, 2) rostro-dorsally directed and an articulation quadrate-mandible (LM 8) behind a hypothetical vertical bar passing through the articulation of the palatin. The mandible of all larvae is less high showing a shorter symphisis mandibulae than the adult one. Except in Chromis sp., adults showed a shorter dentigerous process of the premaxilla (LMs 1, 6) than larvae. The ascending process of the premaxilla (LMs 1, 2) is always longer in adults.
Tests for differences in developmental rates in the eight damselfish species
Suspensorium and opercle
Comparing levels and patterns of shape disparity at three ontogenetic stages
The shape disparity significantly increased over ontogeny within both main trophic groups (i.e. mainly zooplanktivorous and mainly algivorous species, Figure 6B-C). Within each trophic group, the increasing of disparity is most important for the neurocranium, the mandible and the premaxilla. Conversely to zooplanktivorous species, the group of mainly algivorous species did not show significant variation of shape disparity of the unit «suspensorium and opercle» over ontogeny.
At settlement, the first two relative warps RW1 and RW2 explain 54.7% of the total shape variance of the neurocranium (Figure 7). Having the lowest scores on the RW1 axis, the two Dascyllus species showed a higher neurocranium than the others. The second component (RW2), explaining 11% of the shape variance, highlights small differences in the shape of the supraoccipital crest and in the relative length of the orbital and the post-orbital regions. The variance explained by RW1 and RW2 is higher at the size of 60 mm SL and at the maximum adult body size (~72% of total shape variance) (Figure 7). Three groups can be easily distinguished along the RW1 axis at 60 mm, a first including the two Chrysiptera species and P. pavo, a second grouping A. sexfasciatus, Chromis sp and S. nigricans, and a third including the two Dascyllus species. Pomacentrus pavo and both Chrysiptera have a proportionally longer neurocranium. Abudefduf sexfasciatus and Chromis sp show a larger supraoccipital crest than the other species of the second group (RW2). At the maximum size, S. nigricans shares a more similar shape to the two Dascyllus species allowing the clear distinction of three groups in the shape space defined by RW1 and RW2. The two Chrysiptera species, P. pavo and S. nigricans, have a shorter and ventrally directed vomer (low RW1 and RW2 scores, see LMs 1-2).
The shape variance of the unit «suspensorium and opercle» is relatively low at the settlement stage (Figure 8). The two Dascyllus species have the lowest RW1 values and P. pavo the highest. Both Dascyllus species show the narrowest opercle and the highest unit «suspensorium and opercle». The second component (RW2) opposes S. nigricans, which shows a horizontal maxillary process of the palatin, to A. sexfasciatus and Chromis sp., in which this process is rostro-dorsally oriented. At 60 mm, the first component (RW1) directly opposes the two Dascyllus species to the two Chrysiptera species and P. pavo; A. sexfasciatus and S. nigricans occupy an intermediate position (Figure 8). The second component (RW2) allows the discrimination of Chromis sp. having a lower suspensorium and a maxillary process rostro-dorsally oriented. At the maximum size, S. nigricans has a ventrally bent maxillary process of the palatin (-RW2) (Figure 8).
At the settlement stage, the first two components explain nearly equal amounts of variance of mandible shape (RW1 = 39% and RW2 = 30%) (Figure 9). Only Chromis sp. and S. nigricans show distinct shapes from the other species. Chromis sp. has the highest RW1 scores, revealing a lower mandible with a short retroarticular. Stegastes nigricans differs from the others along RW2, showing a shorter angular in its ventral part. At 60 mm and maximum body size, the first two relative warps seem to express almost the same percentage of total shape variance and the same pattern of shape variation (Figure 9). For both stages, RW1 (~ 50% of the total shape variance) mainly differentiates species according to the height of the mandible and the length of the symphisis mandibulae, distinguishing A. sexfasciatus and Chromis sp. from the other species. The second relative warp (24% of the shape variance) distinguishes the three herbivorous species (the two Chrysisptera and S. nigricans) from P. pavo and the two Dascyllus by their massive mandible (Figure 9).
The dynamic of the pattern of premaxilla shape disparity seems similar to that of the mandible (Figure 10). At settlement stage, Chromis sp. shows a highly divergent shape along RW1 revealing a thin dentigerous process and an obtuse angle between the ascending and the dentigerous processes. The second relative warp distinguishes A. sexfasciatus and S. nigricans from the two Chrysiptera, the two Dascyllus species and P. pavo. Indeed, the first ones have a shorter dentigerous process and a longer ascending process than the others. The pattern of shape disparity is the same for the intermediate and the adult stages (Figure 10). Chromis sp. is highly divergent from all other species along RW1; and A. sexfasciatus and the two Chrysiptera species are distinguished from the two Dascyllus species, S. nigricans and P. pavo along RW2.
Comparisons between hyperplanes of skeletal shapes
Susp. & opercle
The level and the pattern of head skeleton shape disparity significantly vary throughout the post-settlement ontogeny of damselfishes. The disparity level increases over ontogeny both at the family level and within each trophic group (Figure 6). Adults are more disparate in shapes and occupy different sub-spaces of the total shape spaces than settling larvae, implying reorganization of variance with growth. Although larval shapes are more similar, they are already species-specific at the time of the coral reef settlement probably due to differences in the larval growth and/or the pelagic larval duration (Table 1). Indeed, all species do not have the same age when they settle on the coral reef .
Diversity of developmental parameters observed in the eight studied species
Rates of development
Lengths of ontogenetic trajectories
Susp. & opercle
The length of ontogenetic trajectory and the rate of shape change appear to be conservative developmental parameters. Depending on the skeletal unit studied, 1 to 4 values have been highlighted for 8 species (Table 6). Thus some species share the same length of ontogenetic trajectory or the same dynamics of shape change. However the variation of both parameters could be underestimated by our methodology. Indeed, having no age information on the studied specimens, compensatory changes in growth rate (defined as the relationships between age and size) could not be explored. These two developmental parameters are not correlated to the duration of the larval phase or the size variation of post-settlement ontogeny. The reasons for such a low variability of these parameters are not obvious. During settlement, the great majority of damselfishes directly recruit to the adult populations and thus directly use the same habitat of their congeners [61, 62] unlike other taxa which undergone one or some habitat changes before their recruitment (e.g. Labridae, Acanthuridae, Serranidae). Among the studied species, only A. sexfasciatus deviates from this general trend: the juveniles settle in micro-atolls of the fringing reef and adults live on the barrier reef . This consistent behavior pattern observed in pomacentrids could partially explain the low inter-specific variability of the rate of development. Both parameters are not related to diet and phylogeny. Intuitively, we hypothesized that planktivorous species at the adult stage, thus having no ontogenetic change of diet, undergo lower amount of shape changes than herbivorous species but that is not verified. The ontogenetic shape changes that are related to a shift of feeding strategies [38, 50] and other functional demands might also explain shape changes related to respiration  or sound production .
The present study describes qualitatively and quantitatively the diversification of the head skeleton during the post-settlement growth of damselfishes. The morphological disparity of the head skeleton is higher at the adult stage associated with the coral reef environment in comparison with the settling larval stage. At reef settlement, larval shapes are already species-specific probably due to differences in age. The process of diversification throughout the ontogeny of damselfishes is mainly explained by the divergence of allometric trajectories in shape space. The length of ontogenetic trajectory and the rate of development seem to be more constrained developmental parameters. The bipartite life-cycle of damselfishes exemplifies a case where the environment (i.e. the coral reef) promotes the increasing shape disparity of the head skeleton over ontogeny. Linked to different developmental parameters, we show adult head morphology diverged on various phenotypes across species, suggesting that functional demands for varied feeding strategy are more intense in adults than in larvae. Our study demonstrates that both selection and developmental processes have influenced phenotypic evolution in this fish group.
We would like to thank Y. Chancerelle, P. Ung (CRIOBE, Moorea, French Polynesia), J.M. Ouin (Aqua-Lab, Institut Halieutique et des Sciences Marines, Toliara, Madagascar) and C. Brié (Tropical Fish Tahiti, Rangiroa atoll, French Polynesia) for providing hospitality, laboratory facilities, and helping to collect the fishes. Special thanks to H. D. Sheets for his guidance with geometric morphometric analyses. We are pleased to acknowledge M. Chardon and E. Parmentier for valuable suggestions during the study. The manuscript benefited greatly from the critical comments on an earlier version by M. L. Zelditch. We also gratefully acknowledge Felix Breden and two anonymous reviewers for their insightful comments and helpful criticism of the original manuscript. This research received financial support from La Communauté Française de Belgique (Concours des bourses de voyage 2007) and from The Belgian National Fund for Scientific Research (FRS-FNRS) (FRFC contract no. 2.4.583.05).
- Woods PJ: Habitat-dependent geographical variation in ontogenetic allometry of the shiner perch Cymatogaster aggregata Gibbons (Teleostei: Embiotocidae). Journal of Evolutionary Biology. 2007, 20 (5): 1783-1798. 10.1111/j.1420-9101.2007.01386.x.View ArticlePubMedGoogle Scholar
- Bellwood DR, Wainwright PC, Fulton CJ, Hoey AS: Functional versatility supports coral reef biodiversity. Proceedings of the Royal Society B-Biological Sciences. 2006, 273 (1582): 101-107. 10.1098/rspb.2005.3276.View ArticlePubMed CentralGoogle Scholar
- Cooper WJ, Westneat MW: Form and function of damselfish skulls: rapid and repeated evolution into a limited number of trophic niches. Bmc Evolutionary Biology. 2009, 9: 10.1186/1471-2148-9-24.Google Scholar
- Frédérich B, Pilet A, Parmentier E, Vandewalle P: Comparative trophic morphology in eight species of damselfishes (Pomacentridae). Journal of Morphology. 2008, 269 (2): 175-188.View ArticlePubMedGoogle Scholar
- Fulton CJ, Bellwood DR, Wainwright PC: Wave energy and swimming performance shape coral reef fish assemblages. Proceedings of the Royal Society B-Biological Sciences. 2005, 272 (1565): 827-832. 10.1098/rspb.2004.3029.View ArticlePubMed CentralGoogle Scholar
- Konow N, Bellwood DR, Wainwright PC, Kerr AM: Evolution of novel jaw joints promote trophic diversity in coral reef fishes. Biological Journal of the Linnean Society. 2008, 93 (3): 545-555. 10.1111/j.1095-8312.2007.00893.x.View ArticleGoogle Scholar
- Motta PJ: Functional morphology of the feeding apparatus of ten species of Pacific butterflyfishes (Perciformes, Chaetodontidae): an ecomorphological approach. Environmental Biology of Fishes. 1988, 22 (1): 39-67. 10.1007/BF00000543.View ArticleGoogle Scholar
- Riedlecker E, Herler J: Trophic morphology of the coral-associated genus Gobiodon (Teleostei: Gobiidae) from the Red Sea. Journal of Zoological Systematics and Evolutionary Research. 2009, 47 (2): 160-170. 10.1111/j.1439-0469.2008.00497.x.View ArticleGoogle Scholar
- Wainwright PC, Bellwood DR: Ecomorphology of feeding in coral reef fishes. Coral reef fishes: dynamics and diversity in a complex ecosystem. Edited by: Sale PF. 2002, San Diego: Academic Press, 33-56.View ArticleGoogle Scholar
- Wainwright PC, Bellwood DR, Westneat MW, Grubich JR, Hoey AS: A functional morphospace for the skull of labrid fishes: patterns of diversity in a complex biomechanical system. Biological Journal of the Linnean Society. 2004, 82 (1): 1-25. 10.1111/j.1095-8312.2004.00313.x.View ArticleGoogle Scholar
- Fulton CJ, Bellwood DR: Ontogenetic habitat use in labrid fishes: an ecomorphological perspective. Marine Ecology-Progress Series. 2002, 236: 255-262. 10.3354/meps236255.View ArticleGoogle Scholar
- Wainwright PC, Richard BA: Predicting patterns of prey use from morphology of fishes. Environmental Biology of Fishes. 1995, 44 (1-3): 97-113. 10.1007/BF00005909.View ArticleGoogle Scholar
- Leis JM: The pelagic stage of reef fishes. The ecology of fishes on coral reefs. Edited by: Sale PF. 1991, San Diego: Academic Press, 183-230.View ArticleGoogle Scholar
- Bellwood DR, Wainwright PC: The history and biogeography of fishes on coral reefs. Coral Reef Fishes: Dynamics and Diversity in a Complex Ecosystem. Edited by: Sale PF. 2002, London: Academic Press, 5-32.View ArticleGoogle Scholar
- Erwin DH: Disparity: Morphological pattern and developmental context. Palaeontology. 2007, 50: 57-73. 10.1111/j.1475-4983.2006.00614.x.View ArticleGoogle Scholar
- Holtmeier CL: Heterochrony, maternal effects, and phenotypic variation among sympatric pupfishes. Evolution. 2001, 55 (2): 330-338.View ArticlePubMedGoogle Scholar
- Loy A, Bertelletti M, Costa C, Ferlin L, Cataudella S: Shape changes and growth trajectories in the early stages of three species of the genus Diplodus (Perciformes, Sparidae). Journal of Morphology. 2001, 250 (1): 24-33. 10.1002/jmor.1056.View ArticlePubMedGoogle Scholar
- Zelditch ML, Sheets HD, Fink WL: The ontogenetic dynamics of shape disparity. Paleobiology. 2003, 29 (1): 139-156. 10.1666/0094-8373(2003)029<0139:TODOSD>2.0.CO;2.View ArticleGoogle Scholar
- Adams DC, Nistri A: Ontogenetic convergence and evolution of foot morphology in European cave salamanders (Family: Plethodontidae). BMC Evolutionary Biology. 2010, 10: 216-10.1186/1471-2148-10-216.View ArticlePubMedPubMed CentralGoogle Scholar
- Ciampaglio CN: Determining the role that ecological and developmental constraints play in controlling disparity: examples from the crinoid and blastozoan fossil record. Evolution & Development. 2002, 4 (3): 170-188.View ArticleGoogle Scholar
- Foote M: Morphological disparity in ordovician-devonian crinoids and the early saturation of morphological space. Paleobiology. 1994, 20 (3): 320-344.Google Scholar
- Foote M: Morphological diversification of paleozoic crinoids. Paleobiology. 1995, 21 (3): 273-299.Google Scholar
- Eble GJ: Contrasting evolutionary flexibility in sister groups: disparity and diversity in Mesozoic atelostomate echinoids. Paleobiology. 2000, 26 (1): 56-79. 10.1666/0094-8373(2000)026<0056:CEFISG>2.0.CO;2.View ArticleGoogle Scholar
- Hall BK: Bauplane, phylotypic stages, and constraint - Why there are so few types of animals. Evolutionary Biology, Vol 29. 1996, 29: 215-261.Google Scholar
- Allen GR: Damselfishes of the world. 1991, Melle: Publication of natural history and pets book, MergusGoogle Scholar
- Cooper WJ, Smith LL, Westneat MW: Exploring the radiation of a diverse reef fish family: Phylogenetics of the damselfishes (Pomacentridae), with new classifications based on molecular analyses of all genera. Molecular Phylogenetics and Evolution. 2009, 52 (1): 1-16. 10.1016/j.ympev.2008.12.010.View ArticleGoogle Scholar
- Sampey A, McKinnon AD, Meekan MG, McCormick MI: Glimpse into guts: overview of the feeding of larvae of tropical shorefishes. Marine Ecology-Progress Series. 2007, 339: 243-257. 10.3354/meps339243.View ArticleGoogle Scholar
- Frédérich B, Fabri G, Lepoint G, Vandewalle P, Parmentier E: Trophic niches of thirteen damselfishes (Pomacentridae) at the Grand Récif of Toliara, Madagascar. Ichthyological Research. 2009, 56 (1): 10-17.View ArticleGoogle Scholar
- Emery AR: Comparative ecology and functional osteology of fourteen species of damselfish (Pisces: Pomacentridae) at Alligator Reef, Florida Keys. Bulletin of Marine Science. 1973, 23: 649-770.Google Scholar
- Webster M: Ontogeny and evolution of the early Cambrian trilobite genus Nephrolenellus (Olenelloidea). Journal of Paleontology. 2007, 81: 1168-1193. 10.1666/06-092.1.View ArticleGoogle Scholar
- Ivanovic A, Vukov TD, Dzukic G, Tomasevic N, Kalezic ML: Ontogeny of skull size and shape changes within a framework of biphasic lifestyle: a case study in six Triturus species (Amphibia, Salamandridae). Zoomorphology. 2007, 126 (3): 173-183. 10.1007/s00435-007-0037-1.View ArticleGoogle Scholar
- Monteiro LR, Cavalcanti MJ, Sommer HJS: Comparative ontogenetic shape changes in the skull of Caiman species (Crocodylia, Alligatoridae). Journal of Morphology. 1997, 231 (1): 53-62. 10.1002/(SICI)1097-4687(199701)231:1<53::AID-JMOR5>3.0.CO;2-P.View ArticleGoogle Scholar
- Cardini A, O'Higgins P: Post-natal ontogeny of the mandible and ventral cranium in Marmota species (Rodentia, Sciuridae): allometry and phylogeny. Zoomorphology. 2005, 124 (4): 189-203. 10.1007/s00435-005-0008-3.View ArticleGoogle Scholar
- Collard M, O'Higgins PO: Ontogeny and homoplasy in the papionin monkey face. Evolution & Development. 2001, 3 (5): 322-331.View ArticleGoogle Scholar
- Bastir M, O'Higgins P, Rosas A: Facial ontogeny in Neanderthals and modern humans. Proceedings of the Royal Society B-Biological Sciences. 2007, 274 (1614): 1125-1132. 10.1098/rspb.2006.0448.View ArticlePubMed CentralGoogle Scholar
- Klingenberg CP: Heterochrony and allometry: the analysis of evolutionary change in ontogeny. Biological Reviews. 1998, 73 (1): 79-123. 10.1017/S000632319800512X.View ArticlePubMedGoogle Scholar
- Webster M, Zelditch ML: Evolutionary modifications of ontogeny: heterochrony and beyond. Paleobiology. 2005, 31 (3): 354-372. 10.1666/0094-8373(2005)031[0354:EMOOHA]2.0.CO;2.View ArticleGoogle Scholar
- Liem KF: Ecomorphology of the teleostean skull. The skull: functional and evolutionary mechanisms. Edited by: Hanken J, Hall BK. 1993, Chicago: The University of Chicago Press, 422-452. vol. 3Google Scholar
- Frédérich B, Lehanse O, Vandewalle P, Lepoint G: Trophic niche width, shift, and specialization of Dascyllus aruanus in Toliara lagoon, Madagascar. Copeia. 2010, 218-226. 2
- Froukh T, Kochzius M: Species boundaries and evolutionary lineages in the blue green damselfishes Chromis viridis and Chromis atripectoralis (Pomacentridae). Journal of Fish Biology. 2008, 72 (2): 451-457. 10.1111/j.1095-8649.2007.01746.x.View ArticleGoogle Scholar
- Dufour V, Galzin R: Colonization patterns of reef fish larvae to the lagoon at Moorea Island, French Polynesia. Marine Ecology-Progress Series. 1993, 102 (1-2): 143-152.View ArticleGoogle Scholar
- Taylor WR, VanDyke GC: Revised procedure for staining and clearing small fishes and other vertebrates for bone and cartilage study. Cybium. 1985, 9: 107-121.Google Scholar
- Adams DC, Rohlf FJ, Slice DE: Geometric morphometrics: ten years of progress following the 'revolution'. Italian Journal of Zoology. 2004, 71 (1): 5-16. 10.1080/11250000409356545.View ArticleGoogle Scholar
- Bookstein F: Morphometric tools for landmark data: geometry and biology. 1991, Cambridge University PressGoogle Scholar
- Bookstein F: Combining the tools of geometric morphometrics. Advances in morphometrics. Edited by: Marcus LF, Corti M, Loy A, Naylor G, Slice D. 1996, New York: Plenum Press, 131-151.View ArticleGoogle Scholar
- Rohlf FJ, Marcus LF: A revolution in morphometrics. Trends in Ecology & Evolution. 1993, 8 (4): 129-132.View ArticleGoogle Scholar
- Rohlf FJ, Slice D: Extensions of the Procrustes method for the optimal superimposition of landmarks. Systematic Zoology. 1990, 39 (1): 40-59. 10.2307/2992207.View ArticleGoogle Scholar
- Zelditch ML, Swiderski DL, Sheets HD, Fink WL: Geometric morphometrics for biologists: A primer. 2004, San Diego: Elsevier Academic PressGoogle Scholar
- Rohlf FJ: Relative warps analysis and an example of its application to mosquito wings. Contributions to morphometrics. Edited by: Marcus LF, Bello E. 1993, A. G-V. Madrid: Monografias del Museo Nacional de Ciencias Naturales, CSIC, 131-159.Google Scholar
- Frédérich B, Adriaens D, Vandewalle P: Ontogenetic shape changes in Pomacentridae (Teleostei, Perciformes) and their relationships with feeding strategies: a geometric morphometric approach. Biological Journal of the Linnean Society. 2008, 95 (1): 92-105.View ArticleGoogle Scholar
- Mitteroecker P, Gunz P, Bookstein FL: Heterochrony and geometric morphometrics: a comparison of cranial growth in Pan paniscus versus Pan troglodytes. Evolution & Development. 2005, 7 (3): 244-258.View ArticleGoogle Scholar
- Monteiro LR: Multivariate regression models and geometric morphometrics: The search for causal factors in the analysis of shape. Systematic Biology. 1999, 48 (1): 192-199. 10.1080/106351599260526.View ArticlePubMedGoogle Scholar
- Zelditch ML, Sheets HD, Fink WL: Spatiotemporal reorganization of growth rates in the evolution of ontogeny. Evolution. 2000, 54 (4): 1363-1371.View ArticlePubMedGoogle Scholar
- Bastir M, Rosas A: Facial heights: Evolutionary relevance of postnatal ontogeny for facial orientation and skull morphology in humans and chimpanzees. Journal of Human Evolution. 2004, 47 (5): 359-381. 10.1016/j.jhevol.2004.08.009.View ArticlePubMedGoogle Scholar
- Darlington RB, Smulders TV: Problems with residual analysis. Animal Behaviour. 2001, 62: 599-602. 10.1006/anbe.2001.1806.View ArticleGoogle Scholar
- Zelditch ML, Lundrigan BL, David Sheets H, Garland T: Do precocial mammals develop at a faster rate? A comparison of rates of skull development in Sigmodon fulviventer and Mus musculus domesticus. Journal of Evolutionary Biology. 2003, 16 (4): 708-720. 10.1046/j.1420-9101.2003.00568.x.View ArticlePubMedGoogle Scholar
- Mezey JG, Houle D: Comparing G matrices: Are common principal components informative?. Genetics. 2003, 165 (1): 411-425.PubMedPubMed CentralGoogle Scholar
- Wellington GM, Victor BC: Planktonic larval duration of 100 species of Pacific and Atlantic damselfishes (Pomacentridae). Marine Biology. 1989, 101 (4): 557-567. 10.1007/BF00541659.View ArticleGoogle Scholar
- Cardini A, Thorington RW: Postnatal ontogeny of marmot (Rodentia, Sciuridae) crania: Allometric trajectories and species divergence. Journal of Mammalogy. 2006, 87 (2): 201-215. 10.1644/05-MAMM-A-242R1.1.View ArticleGoogle Scholar
- Sanfelice D, de Freitas TRO: The ontogeny of shape disparity in three species of Otariids (Pinnipedia: Mammalia). Latin American Journal of Aquatic Mammals. 2007, 6 (2): 139-154.View ArticleGoogle Scholar
- Lecchini D, Galzin R: Spatial repartition and ontogenetic shifts in habitat use by coral reef fishes (Moorea, French Polynesia). Marine Biology. 2005, 147 (1): 47-58. 10.1007/s00227-004-1543-z.View ArticleGoogle Scholar
- McCormick MI, Makey LJ: Post-settlement transition in coral reef fishes: overlooked complexity in niche shifts. Marine Ecology-Progress Series. 1997, 153: 247-257. 10.3354/meps153247.View ArticleGoogle Scholar
- Osse JWM: Form changes in fish larvae in relation to changing demands of function. Netherlands Journal of Zoology. 1990, 40 (1-2): 362-385.Google Scholar
- Parmentier E, Colleye O, Fine ML, Frédérich B, Vandewalle P, Herrel A: Sound production in the clownfish Amphiprion clarkii. Science. 2007, 316 (5827): 1006-10.1126/science.1139753.View ArticlePubMedGoogle Scholar
- Pratchett MS, Gust N, Goby G, Klanten SO: Consumption of coral propagules represents a significant trophic link between corals and reef fish. Coral Reefs. 2001, 20 (1): 13-17. 10.1007/s003380000113.View ArticleGoogle Scholar
- Frédérich B, Parmentier E, Vandewalle P: A preliminary study of development of the buccal apparatus in Pomacentridae (Teleostei, Perciformes). Animal Biology. 2006, 56 (3): 351-372.View ArticleGoogle Scholar
- Gluckmann I, Vandewalle P: Morphofunctional analysis of the feeding apparatus in four Pomacentridae species: Dascyllus aruanus, Chromis retrofasciata, Chrysiptera biocellata and C. unimaculata. Italian Journal of Zoology. 1998, 65: 421-424. 10.1080/11250009809386858.View ArticleGoogle Scholar
- Alfaro ME, Bolnick DI, Wainwright PC: Evolutionary consequences of many-to-one mapping of jaw morphology to mechanics in labrid fishes. American Naturalist. 2005, 165 (6): E140-E154. 10.1086/429564.View ArticlePubMedGoogle Scholar
- Wainwright PC, Alfaro ME, Bolnick DI, Hulsey CD: Many-to-one mapping of form to function: A general principle in organismal design?. Integrative and Comparative Biology. 2005, 45 (2): 256-262. 10.1093/icb/45.2.256.View ArticlePubMedGoogle Scholar
- Collar DC, Wainwright PC: Discordance between morphological and mechanical diversity in the feeding mechanism of centrarchid fishes. Evolution. 2006, 60 (12): 2575-2584.View ArticlePubMedGoogle Scholar
- Hulsey CD, Wainwright PC: Projecting mechanics into morphospace: disparity in the feeding system of labrid fishes. Proceedings of the Royal Society of London Series B-Biological Sciences. 2002, 269 (1488): 317-326. 10.1098/rspb.2001.1874.View ArticleGoogle Scholar
- Coughlin DJ, Strickler JR: Zooplankton capture by a coral-reef fish - an adaptive response to evasive prey. Environmental Biology of Fishes. 1990, 29 (1): 35-42. 10.1007/BF00000566.View ArticleGoogle Scholar
- Kavanagh KD, Alford RA: Sensory and skeletal development and growth in relation to the duration of the embryonic and larval stages in damselfishes (Pomacentridae). Biological Journal of the Linnean Society. 2003, 80 (2): 187-206. 10.1046/j.1095-8312.2003.00229.x.View ArticleGoogle Scholar
- Lo-Yat A: Variabilité temporelle de la colonisation par les larves de poissons de l'atoll de Rangiroa (Tuamotu, Polynésie Française) et utilisation de l'outil "otolithe" de ces larves. 2002, Tahiti: Université de Polynésie françaiseGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.