Adaptive phenotypic plasticity in the Midas cichlid fish pharyngeal jaw and its relevance in adaptive radiation
© Muschick et al; licensee BioMed Central Ltd. 2011
Received: 18 January 2011
Accepted: 30 April 2011
Published: 30 April 2011
Phenotypic evolution and its role in the diversification of organisms is a central topic in evolutionary biology. A neglected factor during the modern evolutionary synthesis, adaptive phenotypic plasticity, more recently attracted the attention of many evolutionary biologists and is now recognized as an important ingredient in both population persistence and diversification. The traits and directions in which an ancestral source population displays phenotypic plasticity might partly determine the trajectories in morphospace, which are accessible for an adaptive radiation, starting from the colonization of a novel environment. In the case of repeated colonizations of similar environments from the same source population this "flexible stem" hypothesis predicts similar phenotypes to arise in repeated subsequent radiations. The Midas Cichlid (Amphilophus spp.) in Nicaragua has radiated in parallel in several crater-lakes seeded by populations originating from the Nicaraguan Great Lakes. Here, we tested phenotypic plasticity in the pharyngeal jaw of Midas Cichlids. The pharyngeal jaw apparatus of cichlids, a second set of jaws functionally decoupled from the oral ones, is known to mediate ecological specialization and often differs strongly between sister-species.
We performed a common garden experiment raising three groups of Midas cichlids on food differing in hardness and calcium content. Analyzing the lower pharyngeal jaw-bones we find significant differences between diet groups qualitatively resembling the differences found between specialized species. Observed differences in pharyngeal jaw expression between groups were attributable to the diet's mechanical resistance, whereas surplus calcium in the diet was not found to be of importance.
The pharyngeal jaw apparatus of Midas Cichlids can be expressed plastically if stimulated mechanically during feeding. Since this trait is commonly differentiated - among other traits - between Midas Cichlid species, its plasticity might be an important factor in Midas Cichlid speciation. The prevalence of pharyngeal jaw differentiation across the Cichlidae further suggests that adaptive phenotypic plasticity in this trait could play an important role in cichlid speciation in general. We discuss several possibilities how the adaptive radiation of Midas Cichlids might have been influenced in this respect.
Adaptive radiations arise through the rapid divergence of an ancestral species into a multitude of morphologically and ecologically differentiated taxa . This process is assumed to be driven by divergent natural selection and ecological speciation where the adaptation to different niches eventually results in the evolution of reproductive isolation . For example, specialization to certain food resources might lead to divergent habitat preferences, which in turn might isolate the populations reproductively [reviewed in ]. Specialization in diet is usually accompanied by morphological adaptations facilitating resource exploitation as has been shown in some textbook examples of adaptive radiation, e.g. the Darwin finches on the Galapagos Islands , the cichlid fishes in East African lakes [5–7], or the cosmopolitan tiger beetles .
Often, adaptive radiations are triggered by an altered adaptive landscape providing opportunity to invade previously not encountered ecological niches (e.g. after colonization of a new environment) or not accessible niches (e.g. after evolution of a 'key innovation') [9, 10]. Recent studies showed that these adaptive peak shifts might happen rapidly [reviewed in ], and raise the question of how the adaptive morphological change drives the shift from one peak to another on the adaptive surface [12, 13]. Mutation in coding and regulatory sequences and selection might not be sufficient to explain the rapidity of ecological adaptation seen in some instances . Adaptation from standing genetic variation is also not likely to apply to all cases of adaptive radiations, particularly those with only a small number of founders . Adaptive phenotypic plasticity might play a key role allowing populations to enter the 'realm of attraction' of a new adaptive peak, in which genetic assimilation occurs through directional selection favoring genotypes that produce even more extreme phenotypes than what would be possible by plastic response of the ancestral genotype alone [16, 17]. Baldwin discussed this topic already in 1896 and described it as 'a new factor in evolution' [18, 19]. Although its importance meanwhile became evident, phenotypic plasticity and genetic assimilation were dismissed as being unimportant during the modern evolutionary synthesis . There has been a recent resurgence of interest in these phenomena [21–25], but the link to diversification is still little explored and under debate [26–28]. Not many investigations of phenotypic plasticity in model systems for speciation research, such as cichlid fishes, have been attempted (but see [29–33]).
The Neotropical Midas Cichlid species complex (Amphilophus spp.), is recognized among evolutionary biologists for its rapid phenotypic diversification and speciation [6, 34]. This species complex has its center of its distribution in Nicaragua, and is comprised of an array of very young species that inhabit both the large Nicaraguan lakes, and several volcanic crater-lakes that contain small scale adaptive radiations [35, 36]. The large Nicaraguan lakes, characterized by relatively turbid and shallow waters, have repeatedly acted as source populations for the colonization of nearby crater-lakes newly formed in the calderas of extinguished volcanoes. In these lakes the Midas cichlids encountered novel environmental conditions - i.e. presence of deeper zones and clearer water - and speciated in situ [34, 35, 37–41]. Crater-lake species have separated along depth and benthic-limnetic axes [34, 35], with the open water column apparently being the first novel habitat invaded. Also, the Midas cichlid species have differentiated in their trophic adaptations. Usage of food sources like stonewort, Aufwuchs, evasive invertebrate prey, fish or snails differs species-specifically . The Midas cichlids species, as well as other Neotropical and Old World cichlids, often differ in the relative degree of hypertrophy of a second set of jaws in the throat - the pharyngeal jaw - derived from branchial arch components and important for food mastication [reviewed in ]. Specialization for feeding on hard-shelled prey like snails, mussels, or crustaceans (durophagy) through this hypertrophy of the pharyngeal jaw apparatus (PJA) has been found to be a common axis of differentiation in crater-lake Midas cichlids as well as in other cichlid groups [5, 31, 32, 34, 42–44]. Its frequency and independency of acquisition across the phylogenetic tree suggests an important role of this adaptation in cichlid speciation [, , reviewed in ]
The Midas cichlid species in the crater lakes are often well differentiated in the trophic apparatus and only a few thousand years old [34–37]. The trophic polymorphism in the Midas crater-lake species could be derived from standing genetic variation, since the polymorphism is present in the large lakes, too [31, 32, 38, 41]. However, the probably limited number of colonizing individuals would render a scenario of the evolution of trait divergence subsequent to colonization also plausible. This scenario is arguably more likely for remote crater-lakes with a monophyletic Midas cichlid assemblage, e.g. Lake Apoyo (see ). A plausible scenario could be that the divergence in the pharyngeal jaw apparatus in the crater lake Midas cichlid species might have been initiated by phenotypic plasticity in the ancestor. Reproductive isolation might then have occurred via habitat isolation through the heterogeneous distribution of snails in Nicaragua's volcanic crater-lakes, where densities appear to be dependent on depth and substrate type . During times of low food availability otherwise opportunistic individuals adapted for durophagy might confine to areas of high snail density and thereby encounter mates non-randomly in respect to their pharyngeal jaw type [31, 32, 46, 47]. If the ancestor of derived species was phenotypically plastic in ecologically relevant traits, this plasticity might have triggered the diversification. The "flexible stem" model, proposed by West-Eberhard , predicts that the directions in phenotypic space in which plasticity is expressed influence the trajectories of phenotypic evolution via genetic accommodation, similar to evolution along "genetic lines of least resistance" . Therefore, it also predicts the outcomes of adaptive radiations seeded by the same ancestor and evolving in similar environments to be similar in terms of their phenotype composition.
In several cichlid fish species (family Cichlidae), plasticity in different traits has been demonstrated: Meyer experimentally induced changes in the oral jaw morphology in the Neotropical cichlid Parachromis managuensis by feeding different diets , a similar procedure was followed by Bouton and coworkers using the African cichlid Neochromis greenwoodi . The Lake Victoria cichlid Haplochromis pyrrhocephalus was almost driven to extinction by the upsurge of the introduced, predatory Nile perch in the 1980s, but was able to adapt morphologically to the new environmental conditions of high predatory pressure and eutrophication in only two decades . It has been interpreted that the speed and complexity of these morphological changes relied on a joined action of phenotypic plasticity and genetic change. The molluscivorous Astatoreochromis alluaudi naturally exhibits molariform pharyngeal jaws (i.e. stout, broad and strong jaw-bones with wide and flat teeth) . However, when raised on soft artificial food under laboratory conditions , in natural conditions in lakes not inhabited by snails , or in lakes inhabited by snails but also with a molluscivorous competitor present , they develop less stout pharyngeal jaws with cuspid teeth (papilliform).
Specializations matter most during ecological "crunch times", when resource availability is low and opportunistic feeding is precluded [42, 46]. The ability to exploit resources then at all or more efficiently than other species can, matters for the individual's survival. But specializations come with a trade-off. The specialization of being able to feed on particular diets especially efficiently often comes at the cost of being much less efficient when dealing with alternative diets. Apparently, such a trade-off exists in the Neotropical Midas Cichlid (Amphilophus cf. citrinellus) between two different types of pharyngeal jaws, molariform and papilliform. Individuals with papilliform lower pharyngeal jaws are more effective when dealing with soft food items . Individuals with molariform jaws, on the other hand, can crack larger and harder snail shells and do this faster than papilliform individuals .
These cases of phenotypic plasticity, the basis of lacustrine cichlid radiations on trophic specialization [44, 55, 56] and the possible causal linkage of plasticity and diversification [23, 30, 31, 57] call for examination of adaptive phenotypic plasticity in trophic traits in an adaptive radiation of cichlids comprising species differentiated in these traits. The lower pharyngeal jaw (LPJ) might constitute 'an ideal component of cichlid trophic morphology' to be investigated in this respect . Preferably, the case in study should have a known and young history, involve colonization of new habitats and tests for plasticity in the ancestral or similar to the ancestral source population.
Here, we tested in a common garden experiment the developmental plasticity of the lower pharyngeal jaw of Amphilophus citrinellus (Günther, 1864) exposed to diets differing in hardness. Earlier work  had suggested that the species in this species complex are phenotypically plastic and that the abundance of molariform fish correlates with the abundance of their major prey item, hard-shelled snails.
The experiment was performed on a laboratory stock derived from the crater Lake Masaya, which was bred in captivity for several decades. Although Lake Masaya is a volcanic crater-lake, its A. citrinellus population is very close to the populations of the Lake Nicaragua - which is probably the ancestral source population of most crater-lake radiations - in terms of body shape  and phylogenetic relationships . Furthermore, it has been suggested that Lake Masaya might have been colonized as recently as 450 years ago .
We investigated whether the development of pharyngeal jaws differed between three types of diets: (1) intact snails with shell, (2) peeled snails without shell, and (3) finely ground up whole snails frozen in pellets, from which fish could nibble off the thawed, soft outer layer when those were given into the water. We aimed to verify whether a hard diet could induce changes in the pharyngeal jaw of the fish, and whether the generation of robust pharyngeal jaws with stout teeth (molariform jaws) was determined by higher calcium content in the diet, or by mechanical stimulation of the jaws when crushing hard food items.
Our study finds that diet can induce changes on the trophic apparatus of the Midas cichlids, and that this changes are related to the mechanical stimulation of the jaws.
Geometric morphometric analyses
Distances in LPJ shape
diet group comparison
'with shell' vs. 'no shell'
'with shell' vs. 'ground'
'no shell' vs. 'ground'
Analyses of weights and lengths
Group comparisons for morphometric data (non-geometric)
WS vs. G
NS vs. G
NS vs. WS
LPJ centroid size
body weight (Ln)
weight × group
LPJ weight (Ln)
body weight (Ln)
weight × group
LPJ height (Ln)
body weight (Ln)
weight × group
Otolith weight (Ln)
body weight (Ln)
weight × group
LPJ height (Ln)
LPJ weight (Ln)
weight × group
Otolith weight (Ln)
LPJ weight (Ln)
weight × group
Phenotypic plasticity has been hypothesized to be able to promote divergence only if it is not complete, i.e. sufficient to achieve the same fitness as if the trait was expressed constitutively . A plastic response would be adaptive if it shifts the phenotype in the direction of a new peak on the adaptive surface, and non-adaptive or maladaptive responses to stressful environments would place the phenotype away from any optimum . Here, we were able to induce an adaptive plastic response in the LPJ of A. citrinellus by feeding different diets. It qualitatively resembles interspecies differences found in nature, although less pronounced.
In our common garden experiment, the changes induced on the fish exposed to a hard shell diet - i.e. horns of the LPJ pointing more outwards, posterior outline less concave, LPJ relatively heavier and possibly increased bone density - mirror those identified as adaptations for mollusk crushing in several other cichlid sister-species pairs [42, 43],very closely related species in the Midas cichlid complex in several crater lakes [34, 37] and in constitutively expressed  or induced  phenotypes in other species. The expression of a relatively hypertrophied pharyngeal jaw due to durophagy resembling adaptations found in specialized molluscivorous fish, and the result that hypertrophication is much weaker when fish are fed with high-calcium, low-impact diet leads to the conclusion that the observed phenotypic plasticity is indeed adaptive. The trade-off in feeding performance between different phenotypes further evidences the adaptive nature of plasticity in this trait .
A surprising finding is that LPJ height did not differ between the experimental groups, since along this dimension divergence is commonly found in non-molluscivorous/molluscivorous species pairs . A possible explanation would be that this trait behaves allometrically with larger and older molariform fishes expressing more re-growing molars thickening the LPJ. A longer common garden experiment might reveal plasticity in this trait as well. An alternative is that LPJ height is simply not plastic, and its evolution is solely governed by mutation and selection that might bring about developmental constraints. Structural constraints and the lack of phenotypic accommodation would be a possible explanation as well. Under this scenario, an increase in LPJ height would not be possible due to prohibitive spatial demands.
Several findings suggest that no specific and adaptive shape difference was induced by a high-calcium diet alone. Only small differences in shape were observed between 'no shell' and 'ground' groups, and those differences did not resemble known adaptations for durophagy. Furthermore, the comparisons including otolith weight show that calcium allocation is strongly biased towards the LPJ in the 'with shell' group but not in the 'ground' group. There, it appears to affect the skeleton evenly as indicated by the group comparison for sagittae weight when correcting for LPJ weight. This corroborates the finding that the mechanical impact on the LPJ during feeding triggers increased calcium allocation towards the jaw and suggests that a high-calcium diet leads to an unspecific increase in calcium deposition.
The sagittae, as well as the other otoliths, grow in small increments throughout the fish's life  and their weight is considered to reflect weight of the individual and availability of calcium during its life. However, Ichii and Mugiya  showed that fish raised on a calcium depleted diet did not show different bone densities after a period of 58 days, but were able to substitute the lacking dietary input of calcium by increasing uptake through the gills from the water. Farrell and Campana  observed that environmental availability of calcium does not affect its deposition on the otolith. These studies have background levels of calcium in both, supplied diet and water, which might differ from levels in our experiment, involve different species and their experiments were conducted significantly shorter. These differences in experimental setup might explain why in our study an effect of calcium availability on bone and otolith growth was observed as opposed to the other studies.
The effects of the mechanical impact were strong enough to exceed anticipated effects of a higher availability of calcium in the 'ground' diet due to facilitated uptake of minerals from the readily processed shells. 'With shell' fish regularly spat out shell fragments during mastication, and Hoogerhoud  reports snail shell pieces to pass the digestive tract of cichlids apparently unharmed. Such observations might explain the slight and non-significant shift towards relatively heavier otoliths in 'ground' fish when accounted for body weight. Several studies on phenotypic plasticity express concerns about the influence of diet quality on developmental differences between treatment groups, so that detrimental effects of a low-quality diet might be mistaken for (adaptive) phenotypic plasticity [33, 47, 65, 66]. Here, we addressed these concerns with our feeding regime. Specifically, we are able to show that induced differences were not due to a lack of calcium in the diet. Even though the studied individuals descended from an inbred line, which has not been subject to artificial selection favoring plasticity in the pharyngeal jaw apparatus, ability to express this trait plastically persisted. This suggests that the plasticity of the LPJ in A. citrinellus might not be a trait under selection itself, but more likely an instance of a hidden reaction norm .
Similarly to the Midas Cichlid, other cichlid species show PJA adaptable or adapted to durophagy: in Neotropical cichlids non-molluscivorous and molluscivorous species, having papilliform and molariform LPJs respectively, often represent closely related sister species pairs . The same trajectory of divergence has been found between trophic morphs of the same species, Herichthys minckleyi, occurring in the Cuatro Ciénegas basin, Mexico. Along the same axis allometric changes happen during the ontogeny of the Mayan Cichlid Cichlasoma urophthalmus, introduced in Southern Florida . The presence of hypertrophied pharyngeal jaws is not restricted to cichlids, or even to freshwater fishes: members of the marine families Sciaenidae, Haemulidae and Carangidae express a similar type of PJA, allowing them to feed on hard-shelled prey. The phylogenetic relationship to species with non-hypertrophied pharyngeal jaws can be close, e.g. congeneric, in these cases as well .
The number of cases of closely related species or trophic morphs of a single species exhibiting such divergent morphologies, as well as their phylogenetic dispersal, is astonishing. This trajectory in morphospace might be similarly important as the well-known deep-bodied vs. elongated body trajectory found in many benthic-limnetic fish species pairs (e.g. [69–71], and those reviewed in ). Both phenotypic contrasts are usually accompanied by extensive diet and/or habitat preference differences, respectively. Such ecological diversification has been shown to be a major factor in empirically studied speciation events and its importance in speciation is well supported by theoretical models [34, 73–75]. In the Midas Cichlid species complex, ecological diversification has been shown to occur along both axes, even in correlation , and probably led to speciation in several cases .
Phenotypic plasticity and rates of diversification
The importance of phenotypic plasticity in population divergence and speciation gained increasing attention in the last years [22, 23, 26, 33, 47, 57, 76–81]. Both studies focusing on single species and studies within a larger comparative framework investigated this link: Nylin & Wahlberg found support for a 'plasticity scenario' for the diversification of nymphaline butterflies during the Tertiary and argued that herbivorous taxa able to occupy several niches were more likely to diversify along with the angiosperm radiation . In coastal San Diego a population of montane dark-eyed juncos (Junco hyemalis, Aves) was able to establish itself due to an adaptive plastic response in reproductive effort . A recent review by Pfennig et al.  summarizes theoretical and empirical studies and diagnoses an important, but largely underappreciated, role of phenotypic plasticity in speciation and adaptive radiation. Comparing sister clade pairs - with one clade being known to include cases of resource polyphenism, while the other does not - Pfennig and McGee found evidence that resource polyphenism is associated with greater species richness in fishes and amphibians .
The role of phenotypic plasticity in population divergence appears to be at least twofold: (1) plasticity increases the probability of population persistence after colonization of a new environment, thus making its split from the ancestral population more likely [83, 84], and (2) plasticity provides means of conquering other peaks on the adaptive landscape, possibly leading to assortative mating and speciation with parallel outcomes in repeated cases [12, 14, 17, 23, 33].
Theoretical investigations support these predictions. Probability of population persistence increases with plasticity while being dependent on the amount of environmental change and the costliness of plasticity . At a moderate rate of environmental change and if plasticity is costly, high levels of plasticity are expected to lead to an increased probability of extinction while an intermediate level improves the ability of persistence . Access to novel ecological niches is improved because an increase in epigenetic variability does facilitate the circumvention of adaptive valleys and smoothes the fitness landscape [13, 86, 87]. Using numerical simulations Thibert-Plante and Hendry  find plasticity to commence reduction in gene flow between populations in contrasting environments. To do so, plasticity must occur before dispersal but could then lead to reproductive isolation even prior to any adaptive genetic divergence.
Our demonstration of adaptive phenotypic plasticity in the LPJ of A. cf. citrinellus suggests that this could be a crucial factor in ecological speciation and adaptive radiation in the repeated Amphilophus crater-lake radiations and possibly in other cichlid clades as well [30, 31]. The results of our experiment support the "flexible stem" hypothesis, in that the induced differences between treatment groups - more robust LPJs in the 'with-shell' group, less robust LPJs in groups fed soft food - resemble between-species differences in crater-lake radiations. However, we did not test for plasticity in the ancestor itself, nor in fish derived from the large Nicaraguan Lakes, but in a stock derived from Lake Masaya. The different history might have caused an alteration of the plastic response in experimental groups compared to the real ancestor. But since there is a considerable chance of the Lake Masaya A. citrinellus population being very young and since plasticity here seems not be lost easily (at least not over several generations), we suggest that our results endorse the "flexible stem" hypothesis for the Midas Cichlid assemblage. Because the induced plasticity does not reach the extent of morphological divergence found between species in nature we conclude, that the expectations from the "adaptive surface model" are fulfilled as well.
In which way exactly phenotypic plasticity and genetic accommodation in the pharyngeal jaw might abet diversification in the Amphilophus species complex remains speculative. A direct influence on the formation of reproductive isolation might be given through enhancement of habitat preference. If individuals expressing the same type of pharyngeal jaw have a higher chance of mating with each other, and gene flow between groups is hampered strongly enough, population subdivision might be initiated. The heterogeneous distribution of snails, if it is stable over time and patches are sufficiently large, might be the basis for habitat preference by jaw type. Alternatively, the hypothesized function of the pharyngeal jaw apparatus in sound production, e.g. during courtship, might bring about assortative mating according to jaw type if female sound preference is divergent as well .
However, even if phenotypic plasticity is less important in sympatric speciation scenarios it might still influence diversification in allopatry [reviewed in ]. By augmenting the probability of population persistence after colonization of a new environment, e.g. a crater-lake, and the possibility of genetic accommodation of plastic trait changes the likelihood of allopatric speciation between ancestral source population and the new colonizing population is increased. It remains unclear, whether or not the repeated endemic radiations of Midas cichlids in Nicaraguan crater-lakes are facilitated by phenotypic plasticity in the pharyngeal jaw or if the constitutively expressed differences in jaw shape between species are a secondary result of speciation driven by other factors. The best documented case of an in-crater-lake diversification, the origination of the Arrow cichlid Amphilophus zaliosus in Lake Apoyo, seems to have been driven by diverging habitat preferences with differences in pharyngeal jaw shape being probably secondary . However, in other, less-well documented cases the hypothesis that adaptations in the pharyngeal jaw apparatus triggered divergence remains valid, but would need to be further investigated.
We demonstrated phenotypic plasticity in the pharyngeal jaw of the cichlid fish Amphilophus citrinellus that is due not to differences in nutritional composition of the diet, but brought about largely by the mode of feeding. This finding might suggest that plasticity plays an important role in diversification.
Future research on how a plastic reaction in one trait could impact the expression of other traits through correlated plastic responses might contribute to the understanding of parallelisms so often encountered in nature. For example, it seems the papilliform pharyngeal jaw type is correlated with fusiform limnetic body shape whereas the molariform jaw type is correlated with deeper, benthic body shape . The extent to which this 'integration of plastic responses'  is determined, and by which factors, still remains to be elucidated. Also, what role a stage of fixed polymorphism plays in the process of diversification, whether it is an intermediate step  or a 'dead-end', remains to be investigated.
How adaptive phenotypic plasticity is mediated genetically is another important issue. In cichlids, the family of bone morphogenetic proteins (BMPs) is known to be involved in shaping bones of the oral and pharyngeal jaws  and might constitute good candidates, along with respective transcription factors and ligands, for the elucidation of the genetics of phenotypic plasticity in the PJA.
Cichlids, are a prime system for speciation research and have an important trophic trait expressed plastically, and therefore constitute a cogent group for investigating the role of adaptive phenotypic plasticity in diversification. Research combining experimental and field studies with modern tools of analysis, such as sensitive group assignment methods or gene expression quantification, will be most rewarding avenues of research to elucidate the link between plasticity and speciation
Common garden experiment
We divided fry of a single Amphilophus citrinellus brood from an inbred line into three similarly sized groups and fed them on diets differing in mechanical durability and calcium content. The three study groups of 30 A. citrinellus individuals each were kept under standardized laboratory conditions with 12 h daylight for a period of six month. The fish stock used (AM-stock at the University of Konstanz) derives from Lake Masaya, a volcanic crater-lake in Nicaragua. Originally, these fish came from the Berkeley stocks of George Barlow who gave some of these fish to the Steinhard Aquarium in San Francisco. In 2001 fish from there were brought to Konstanz and are the stock of A. cf. citrinellus that were used in these experiments. This fish stock has been bred in captivity on soft artificial food for several decades. Moreover, in Lake Masaya no snails occur and neither are cichlids with molariform pharyngeal jaws reported [54, 91].
The fish groups were raised on different diets: (1) Melanoides tuberculata snails, laboratory grown, with intact bodies and intact or slightly damaged shells (in case the snail was deemed too large), (2) snail bodies, where the shells were manually removed, and (3) M. tuberculata with shell but ground to fine paste using mortar and pestle, which was given frozen in pieces to large to be swallowed as a whole. Food amount was adjusted to match group's estimated size gain. Fish were kept in one large tank (1.8 × 0.5 × 0.5 meter, 450 l) and perforated walls allowed water exchange between the compartments containing the three experimental groups. To counteract position bias, we swapped groups between compartments several times throughout the experiment.
Measurements & analyses
Fishes were sacrificed and weighed, and standard and total length were recorded. We excised LPJs and sagittae, and cleaned and dried them. LPJs and otoliths were weighed to the nearest milligram. LPJs were scanned on a standard desktop scanner. Coordinates of 8 landmarks were recorded for each LPJ using tpsDig 2.11 (, for landmark positions see Figure 1). Landmarks represented homologous, defined locations on the jaws outline. Their positioning followed Klingenberg et al.  with the exception of their landmarks 5 and 6 - instead the anterior tip was covered by our landmark 8. Otherwise landmark position were the same, though differently numbered. Landmark arrangements were procrustes aligned, i.e. their positional, rotational and size information was removed from the dataset. However, size information was recorded in centroid size and was used for joint analysis with other data. Since the LPJ is a symmetrical structure we extracted the symmetric component of shape variation using MorphoJ . We conducted discriminant function analyses (DFA) for each pair of groups to produce Figure 1. A canonical variates analyses (CVA) using residuals of a pooled-within-diet-groups regression on body weight (Ln) yielded mean shape distances and their significance levels were assessed by permutation testing (10.000 permutations).
Fish body weight, LPJ weight, height, and centroid size, and otolith weight were evaluated via analysis of variance (ANOVA) and group-pairwise differences of residuals means were assessed for significance using Tukey's honest significant difference-test. All these measures were Ln transformed prior to analysis. For otoliths the weight of the heavier sagitta was used, to minimize influence of preparation damage.
All statistical tests on length and weight data were performed using the R statistical environment .
MM was a Master's student in AM's laboratory and is now a Ph.D. student with WS and interested in the phenomenon of convergent evolution and its implications for speciation and adaptive radiations.
MB was a postdoc in AM's laboratory when this study was conducted. She is an evolutionary ecologist interested in speciation and the origin of adaptive radiations.
WS was a postdoc in AM's laboratory when this study was conducted. He is an evolutionary biologist interested in the evolution of adaptive radiations of cichlid fishes.
AM is an evolutionary biologist interested in speciation and the origin of evolutionary diversity.
We would like to thank Nathalie Feiner for assistance with the experiments, Christian Klingenberg for advice on geometric morphometrics and Daniel Berner for support in statistics, as well as the staff of University of Konstanz' animal care facility for caretaking of the fish. Dean Adams, Hannes Svardal, Michael Matschiner and four anonymous reviewers provided helpful comments on earlier versions of this manuscript. This work was funded by the Deutsche Forschungsgemeinschaft priority program 1127 ('Radiations') (to MB and AM) and by the University of Konstanz.
- Schluter D: The Ecology of Adaptive Radiation. 2000, New York: Oxford University PressGoogle Scholar
- Skulason S, Smith TB: Resource Polymorphisms in Vertebrates. Trends Ecol Evol. 1995, 10: 366-370. 10.1016/S0169-5347(00)89135-1.View ArticlePubMedGoogle Scholar
- Coyne JA, Orr HA: Speciation. 2004, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- Grant PR, Grant BR: How and Why Species Multiply: The Radiation of Darwin's Finches. 2007, Princeton, New Jersey: Princeton University PressGoogle Scholar
- Meyer A: Phylogenetic relationships and evolutionary processes in East African cichlid fishes. Trends Ecol Evol. 1993, 8: 279-284. 10.1016/0169-5347(93)90255-N.View ArticlePubMedGoogle Scholar
- Barlow GW: The Cichlid Fishes: Nature's Grand Experiment in Evolution. 2000, New York: Perseus PublishingGoogle Scholar
- Fryer G, Iles TD: The Cichlid Fishes of the Great Lakes of Africa: Their Biology and Evolution. 1972, Edinburgh: Oliver & BoydGoogle Scholar
- Pearson DL, Vogler AP: Tiger Beetles: The Evolution, Ecology, and Diversity of the Cicindelids. 2001, Ithaca, New York: Cornell University PressGoogle Scholar
- Yoder JB, Clancey E, Des Roches S, Eastman JM, Gentry L, Godsoe W, Hagey TJ, Jochimsen D, Oswald BP, Robertson J, et al: Ecological opportunity and the origin of adaptive radiations. J Evol Biol. 2010, 23: 1581-1596. 10.1111/j.1420-9101.2010.02029.x.View ArticlePubMedGoogle Scholar
- Gavrilets S, Vose A: Dynamic patterns of adaptive radiation. Proc Natl Acad Sci USA. 2005, 102: 18040-18045. 10.1073/pnas.0506330102.View ArticlePubMedPubMed CentralGoogle Scholar
- Hendry AP, Nosil P, Rieseberg LH: The speed of ecological speciation. Funct Ecol. 2007, 21: 455-464. 10.1111/j.1365-2435.2007.01240.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Fear KK, Price T: The adaptive surface in ecology. Oikos. 1998, 82: 440-448. 10.2307/3546365.View ArticleGoogle Scholar
- Gavrilets S: Fitness Landscapes and the Origin of Species. 2004, Princeton, New Jersey: Princeton University PressGoogle Scholar
- Schlichting CD, Pigliucci M: Phenotypic Evolution: A Reaction Norm Perspective. 1998, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Barrett RDH, Schluter D: Adaptation from standing genetic variation. Trends Ecol Evol. 2008, 23: 38-44. 10.1016/j.tree.2007.09.008.View ArticlePubMedGoogle Scholar
- Waddington CH: Genetic assimilation. Adv Genet. 1961, 10: 257-293.PubMedGoogle Scholar
- Price TD, Qvarnstrom A, Irwin DE: The role of phenotypic plasticity in driving genetic evolution. Proc R Soc Lond, Ser B: Biol Sci. 2003, 270: 1433-1440. 10.1098/rspb.2003.2372.View ArticleGoogle Scholar
- Baldwin JM: A new factor in evolution. Am Nat. 1896, 30: 441-451. 10.1086/276408.View ArticleGoogle Scholar
- Baldwin JM: A new factor in evolution (Continued). Am Nat. 1896, 30: 536-553. 10.1086/276428.View ArticleGoogle Scholar
- Schlichting CD: The role of phenotypic plasticity in diversification. Phenotypic Plasticity: Functional and Conceptual Approaches. Edited by: DeWitt TJ, Scheiner SM. 2004, New York: Oxford University Press, 191-200.Google Scholar
- Crispo E, DiBattista JD, Correa C, Thibert-Plante X, McKellar AE, Schwartz AK, Berner D, De Leon LF, Hendry AP: The evolution of phenotypic plasticity in response to anthropogenic disturbance. Evol Ecol Res. 2010, 12: 47-66.Google Scholar
- West-Eberhard M: Phenotypic plasticity and the origins of diversity. Annu Rev Ecol Syst. 1989, 20: 249-278.View ArticleGoogle Scholar
- West-Eberhard MJ: Developmental Plasticity and Evolution. 2003, New York: Oxford University PressGoogle Scholar
- West-Eberhard MJ: Developmental plasticity and the origin of species differences. P Natl Acad Sci USA. 2005, 102: 6543-6549. 10.1073/pnas.0501844102.View ArticleGoogle Scholar
- Lande R: Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation. J Evol Biol. 2009, 22: 1435-1446. 10.1111/j.1420-9101.2009.01754.x.View ArticlePubMedGoogle Scholar
- Thibert-Plante X, Hendry AP: The consequences of phenotypic plasticity for ecological speciation. J Evol Biol. 2010, 24: 326-342.View ArticlePubMedGoogle Scholar
- Crispo E: The Baldwin effect and genetic assimilation: Revisiting two mechanisms of evolutionary change mediated by phenotypic plasticity. Evolution. 2007, 61: 2469-2479. 10.1111/j.1558-5646.2007.00203.x.View ArticlePubMedGoogle Scholar
- Pfennig DW, McGee M: Resource polyphenism increases species richness: a test of the hypothesis. Philos T R Soc B. 2010, 365: 577-591. 10.1098/rstb.2009.0244.View ArticleGoogle Scholar
- Crispo E, Chapman LJ: Geographic variation in phenotypic plasticity in response to dissolved oxygen in an African cichlid fish. J Evol Biol. 2010, 23: 2091-2103. 10.1111/j.1420-9101.2010.02069.x.View ArticlePubMedGoogle Scholar
- Meyer A: Phenotypic plasticity and heterochrony in Cichlasoma managuense (Pisces, Cichlidae) and their implications for speciation in cichlid fishes. Evolution. 1987, 41: 1357-1369. 10.2307/2409100.View ArticleGoogle Scholar
- Meyer A: Ecological and evolutionary consequences of the trophic polymorphism in Cichlasoma citrinellum (Pisces: Cichlidae). Biol J Linn Soc. 1990, 39: 279-299. 10.1111/j.1095-8312.1990.tb00517.x.View ArticleGoogle Scholar
- Meyer A: Morphometrics and allometry in the trophically polymorphic cichlid fish, Cichlasoma citrinellum: Alternative adaptations and ontogenetic changes in shape. J Zool. 1990, 221: 237-260. 10.1111/j.1469-7998.1990.tb03994.x.View ArticleGoogle Scholar
- Wund MA, Baker JA, Clancy B, Golub JL, Foster SA: A test of the "flexible stem" model of evolution: ancestral plasticity, genetic accommodation, and morphological divergence in the threespine stickleback radiation. Am Nat. 2008, 172: 449-462. 10.1086/590966.View ArticlePubMedGoogle Scholar
- Barluenga M, Stölting KN, Salzburger W, Muschick M, Meyer A: Sympatric speciation in Nicaraguan crater lake cichlid fish. Nature. 2006, 439: 719-723. 10.1038/nature04325.View ArticlePubMedGoogle Scholar
- Elmer KR, Kusche H, Lehtonen TK, Meyer A: Local variation and parallel evolution: morphological and genetic diversity across a species complex of neotropical crater lake cichlid fishes. Philos T R Soc B. 2010, 365: 1763-1782. 10.1098/rstb.2009.0271.View ArticleGoogle Scholar
- Barluenga M, Meyer A: Phylogeography, colonization and population history of the Midas cichlid species complex (Amphilophus spp.) in the Nicaraguan crater lakes. BMC Evol Biol. 2010, 10: 326-10.1186/1471-2148-10-326.View ArticlePubMedPubMed CentralGoogle Scholar
- Elmer KR, Lehtonen TK, Kautt AF, Harrod C, Meyer A: Rapid sympatric ecological differentiation of crater lake cichlid fishes within historic times. BMC Biol. 2010, 8: --View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson AB, Noack-Kunnmann K, Meyer A: Incipient speciation in sympatric Nicaraguan crater lake cichlid fishes: sexual selection versus ecological diversification. Proc R Soc Lond B Biol Sci. 2000, 267: 2133-2141. 10.1098/rspb.2000.1260.View ArticleGoogle Scholar
- Barlow GW: The Midas Cichlid in Nicaragua. Investigations of the Ichthyofauna of Nicaraguan lakes. Edited by: Thorson TB. 1976, Lincoln, NB: University of Nebraska Press, 333-358.Google Scholar
- Barlow GW, Munsey JW: The red devil-Midas-arrow cichlid species complex in Nicaragua. Investigations of the Ichthyofauna of Nicaraguan lakes. Edited by: Thorson TB. 1976, Lincoln, NB: University of Nebraska Press, 359-369.Google Scholar
- Barluenga M, Meyer A: The Midas cichlid species complex: incipient sympatric speciation in Nicaraguan cichlid fishes?. Mol Ecol. 2004, 13: 2061-2076. 10.1111/j.1365-294X.2004.02211.x.View ArticlePubMedGoogle Scholar
- Meyer A: Trophic polymorphisms in cichlid fish: Do they represent intermediate steps during sympatric speciation and explain their rapid adaptive radiation?. Trends in Ichthyology. Edited by: Schroder J-H, Bauer J, Schartl M. 1993, Oxford, UK: Blackwell Science, Ltd, 257-266.Google Scholar
- Hulsey CD, Roberts RJ, Lin ASP, Guldberg R, Streelman JT: Convergence in a mechanically complex phenotype: Detecting structural adaptations for crushing in cichlid fish. Evolution. 2008, 62: 1587-1599. 10.1111/j.1558-5646.2008.00384.x.View ArticlePubMedGoogle Scholar
- Liem KF, Osse JWM: Biological Versatility, Evolution, and Food Resource Exploitation in African Cichlid Fishes. Am Zool. 1975, 15: 427-454.View ArticleGoogle Scholar
- McCrary JK, Madsen H, Gonzalez L, Luna I, Lopez LJ: Comparison of gastropod mollusc (Apogastropoda: Hydrobiidae) habitats in two crater lakes in Nicaragua. Rev Biol Trop. 2008, 56: 113-120.PubMedGoogle Scholar
- Liem KF, Kaufman LS: Intraspecific macroevolution: functional biology of the polymorphic cichlid species Cichlasoma minckleyi. Evolution of Fish Species Flocks. Edited by: Echelle AA, Kornfield I. 1984, Orono, Maine: University of Maine Press, 203-215.Google Scholar
- Wimberger PH: Trophic polymorphisms, plasticity, and speciation in vertebrates. Theory and Application of Fish Feeding Ecology. Edited by: Stouder DJ, Fresh KL, Feller RJ. 1994, Columbia, South Carolina: University of South Carolina Press, 19-43.Google Scholar
- Schluter D: Adaptive radiation along genetic lines of least resistance. Evolution. 1996, 50: 1766-1774. 10.2307/2410734.View ArticleGoogle Scholar
- Bouton N, Witte F, van Alphen JJ: Experimental evidence for adaptive phenotypic plasticity in a rock-dwelling cichlid fish from Lake Victoria. Biol J Linn Soc. 2002, 77: 185-192. 10.1046/j.1095-8312.2002.00093.x.View ArticleGoogle Scholar
- Witte F, Welten M, Heemskerk M, van der Stap I, Ham L, Rutjes H, Wanink J: Major morphological changes in a Lake Victoria cichlid fish within two decades. Biol J Linn Soc. 2008, 94: 41-52. 10.1111/j.1095-8312.2008.00971.x.View ArticleGoogle Scholar
- Greenwood PH: Environmental effects on the pharyngeal mill of a cichlid fish, Astatoreochromis alluaudi, and their taxonomic implications. Proc Linnean Soc Lond. 1965, 176: 1-10.View ArticleGoogle Scholar
- Huysseune A: Phenotypic plasticity in the lower pharyngeal jaw dentition of Astatoreochromis alluaudi (Teleostei: Cichlidae). Arch Oral Biol. 1995, 40: 1005-1014. 10.1016/0003-9969(95)00074-Y.View ArticlePubMedGoogle Scholar
- Hoogerhoud RJC: Ecological Morphology of Some Cichlid Fishes. unpublished PhD thesis. 1986, University of LeidenGoogle Scholar
- Meyer A: Cost of morphological specialization: feeding performance of the two morphs in the trophically polymorphic cichlid fish, Cichlasoma citrinellum. Oecologia. 1989, 80: 431-436. 10.1007/BF00379047.View ArticleGoogle Scholar
- Greenwood PH: Explosive speciation in African lakes. Proc Roy Inst Gr Brit. 1964, 40: 256-269.Google Scholar
- Greenwood PH: Morphology, endemism and speciation in African cichlid fishes. Verh Dtsch Zool Ges. 1973, 66: 115-124.Google Scholar
- Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP: Phenotypic plasticity's impacts on diversification and speciation. Trends Ecol Evol. 2010, 25: 459-467. 10.1016/j.tree.2010.05.006.View ArticlePubMedGoogle Scholar
- Villa J: Ichthyology of the lakes of Nicaragua: historical perspective. Investigations of the Ichthyofauna of Nicaraguan lakes. Edited by: Thorson TB. 1976, Lincoln, Nebraska: University of Nebraska Press, 101-113.Google Scholar
- Ghalambor CK, McKay JK, Carroll SP, Reznick DN: Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007, 21: 394-407. 10.1111/j.1365-2435.2007.01283.x.View ArticleGoogle Scholar
- Trapani J: A morphometric analysis of polymorphism in the pharyngeal dentition of Cichlasoma minckleyi (Teleostei: Cichlidae). Arch Oral Biol. 2004, 49: 825-835. 10.1016/j.archoralbio.2004.03.003.View ArticlePubMedGoogle Scholar
- Witte F, Barel CDN, Hoogerhoud RJC: Phenotypic plasticity of anatomical structures and its ecomorphological significance. Neth J Zool. 1990, 40: 278-298.View ArticleGoogle Scholar
- Pannella G: Fish otoliths - daily growth layers and periodical patterns. Science. 1971, 173: 1124-1127. 10.1126/science.173.4002.1124.View ArticleGoogle Scholar
- Ichii T, Mugiya Y: Effects of a dietary deficiency in calcium on growth and calcium-uptake from the aquatic environment in the goldfish, Carassius auratus. Comparative Biochemistry and Physiology a-Physiology. 1983, 74: 259-262. 10.1016/0300-9629(83)90597-2.View ArticleGoogle Scholar
- Farrell J, Campana SE: Regulation of calcium and strontium deposition on the otoliths of juvenile tilapia, Oreochromis niloticus. Comparative Biochemistry and Physiology a-Physiology. 1996, 115: 103-109. 10.1016/0300-9629(96)00015-1.View ArticleGoogle Scholar
- Myers P, Lundrigan BL, Gillespie BW, Zelditch ML: Phenotypic plasticity in skull and dental morphology in the prairie deer mouse (Peromyscus maniculatus bairdii). J Morphol. 1996, 229: 229-237. 10.1002/(SICI)1097-4687(199608)229:2<229::AID-JMOR7>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Wimberger PH: Effects of Vitamin-C-Deficiency on Body Shape and Skull Osteology in Geophagus-Brasiliensis - Implications for Interpretations of Morphological Plasticity. Copeia. , 1993: 343-351.Google Scholar
- Bergmann GT, Motta PJ: Diet and morphology through ontogeny of the nonindigenous Mayan cichlid 'Cichlasoma (Nandopsis)' urophthalmus (Günther 1862) in southern Florida. Environ Biol Fishes. 2005, 72: 205-211. 10.1007/s10641-004-1480-1.View ArticleGoogle Scholar
- Grubich J: Morphological convergence of pharyngeal jaw structure in durophagous perciform fish. Biol J Linn Soc. 2003, 80: 147-165. 10.1046/j.1095-8312.2003.00231.x.View ArticleGoogle Scholar
- Olafsdottir GA, Snorrason SS, Ritchie MG: Postglacial intra-lacustrine divergence of Icelandic threespine stickleback morphs in three neovolcanic lakes. J Evol Biol. 2007, 20: 1870-1881. 10.1111/j.1420-9101.2007.01375.x.View ArticlePubMedGoogle Scholar
- Mcphail JD: Ecology and Evolution of Sympatric Sticklebacks (Gasterosteus) - Morphological and Genetic-Evidence for a Species Pair in Enos Lake, British-Columbia. Can J Zool. 1984, 62: 1402-1408. 10.1139/z84-201.View ArticleGoogle Scholar
- Fenderson OC: Evidence of subpopulations of lake whitefish, Coregonus clupeaformis, involving a dwarfed form. Trans Am Fish Soc. 1964, 93: 77-94. 10.1577/1548-8659(1964)93[77:EOSOLW]2.0.CO;2.View ArticleGoogle Scholar
- Taylor EB: Species pairs of North temperate freshwater fishes: Evolution, taxonomy, and conservation. Rev Fish Biol Fish. 1999, 9: 299-324. 10.1023/A:1008955229420.View ArticleGoogle Scholar
- McKinnon JS, Mori S, Blackman BK, David L, Kingsley DM, Jamieson L, Chou J, Schluter D: Evidence for ecology's role in speciation. Nature. 2004, 429: 294-298. 10.1038/nature02556.View ArticlePubMedGoogle Scholar
- Gavrilets S, Vose A, Barluenga M, Salzburger W, Meyer A: Case studies and mathematical models of ecological speciation. 1. Cichlids in a crater lake. Mol Ecol. 2007, 16: 2893-2909. 10.1111/j.1365-294X.2007.03305.x.View ArticlePubMedGoogle Scholar
- Dieckmann U, Doebeli M: On the origin of species by sympatric speciation. Nature. 1999, 400: 354-357. 10.1038/22521.View ArticlePubMedGoogle Scholar
- Scoville AG, Pfrender ME: Phenotypic plasticity facilitates recurrent rapid adaptation to introduced predators. Proc Natl Acad Sci USA. 2010, 107: 4260-4263. 10.1073/pnas.0912748107.View ArticlePubMedPubMed CentralGoogle Scholar
- Gomez-Mestre I, Buchholz DR: Developmental plasticity mirrors differences among taxa in spadefoot toads linking plasticity and diversity. Proc Natl Acad Sci USA. 2006, 103: 19021-19026. 10.1073/pnas.0603562103.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams CE, Huntingford FA: Incipient speciation driven by phenotypic plasticity? Evidence form sympatric populations of Arctic charr. Biol J Linn Soc. 2004, 81: 611-618. 10.1111/j.1095-8312.2004.00314.x.View ArticleGoogle Scholar
- Stauffer JR, Van Snik Gray E: Phenotypic plasticity: its role in trophic radiation and explosive speciation in cichlids (Teleostei: Cichlidae). Anim Biol. 2004, 54: 137-158. 10.1163/1570756041445191.View ArticleGoogle Scholar
- Pfennig DW, Murphy PJ: How fluctuating competition and phenotypic plasticity mediate species divergence. Evolution. 2002, 56: 1217-1228.View ArticlePubMedGoogle Scholar
- Parsons KJ, Robinson BW: Replicated evolution of integrated plastic responses during early adaptive divergence. Evolution. 2006, 60: 801-813.View ArticlePubMedGoogle Scholar
- Nylin S, Wahlberg N: Does plasticity drive speciation? Host-plant shifts and diversification in nymphaline butterflies (Lepidoptera: Nymphalidae) during the tertiary. Biol J Linn Soc. 2008, 94: 115-130. 10.1111/j.1095-8312.2008.00964.x.View ArticleGoogle Scholar
- Yeh PJ, Price TD: Adaptive phenotypic plasticity and the successful colonization of a novel environment. Am Nat. 2004, 164: 531-542. 10.1086/423825.View ArticlePubMedGoogle Scholar
- Williams DG, Mack RN, Black RA: Ecophysiology of introduced Pennisetum setaceum on Hawaii: the role of phenotypic plasticity. Ecology. 1995, 76: 1569-1580. 10.2307/1938158.View ArticleGoogle Scholar
- Chevin LM, Lande R, Mace GM: Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory. PLoS Biol. 2010, 8: e1000357-10.1371/journal.pbio.1000357.View ArticlePubMedPubMed CentralGoogle Scholar
- Pal C, Miklos I: Epigenetic inheritance, genetic assimilation and speciation. J Theor Biol. 1999, 200: 19-37. 10.1006/jtbi.1999.0974.View ArticlePubMedGoogle Scholar
- Borenstein E, Meilijson I, Ruppin E: The effect of phenotypic plasticity on evolution in multipeaked fitness landscapes. J Evol Biol. 2006, 19: 1555-1570. 10.1111/j.1420-9101.2006.01125.x.View ArticlePubMedGoogle Scholar
- Rice AN, Lobel PS: The pharyngeal jaw apparatus of the Cichlidae and Pomacentridae: function in feeding and sound production. Rev Fish Biol Fish. 2003, 13: 433-444. 10.1007/s11160-004-8794-0.View ArticleGoogle Scholar
- Agrawal AA: Ecology - Phenotypic plasticity in the interactions and evolution of species. Science. 2001, 294: 321-326. 10.1126/science.1060701.View ArticlePubMedGoogle Scholar
- Albertson RC, Kocher TD: Genetic and developmental basis of cichlid trophic diversity. Heredity. 2006, 97: 211-221. 10.1038/sj.hdy.6800864.View ArticlePubMedGoogle Scholar
- Klingenberg CP, Barluenga M, Meyer A: Body shape variation in cichlid fishes of the Amphilophus citrinellus species complex. Biol J Linn Soc. 2003, 80: 397-408. 10.1046/j.1095-8312.2003.00246.x.View ArticleGoogle Scholar
- Rohlf FJ: tpsDig Version 2.11. 2008, Ecology and Evolution, SUNY, Stony Brook, New York, [http://life.bio.sunysb.edu/morph/]Google Scholar
- Klingenberg CP, Barluenga M, Meyer A: Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry. Evolution Int J Org Evolution. 2002, 56: 1909-1920.View ArticleGoogle Scholar
- Klingenberg CP: MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011, 11: 353-357. [http://www.flywings.org.uk/MorphoJ_page.htm]View ArticlePubMedGoogle Scholar
- R Development Core Team: R: a language and environment for statistical computing. 2007, R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]Google 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.