About this Journal Submit a Manuscript Table of Contents
International Journal of Evolutionary Biology
Volume 2011 (2011), Article ID 641501, 10 pages
Research Article

Modularity of the Oral Jaws Is Linked to Repeated Changes in the Craniofacial Shape of African Cichlids

Department of Biology, Syracuse University, 107 College Place, Syracuse, NY 13244, USA

Received 16 December 2010; Accepted 14 February 2011

Academic Editor: Stephan Koblmüller

Copyright © 2011 Kevin J. Parsons et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The African cichlids of the East-African rift-lakes provide one of the most dramatic examples of adaptive radiation known. It has long been thought that functional decoupling of the oral and pharyngeal jaws in cichlids has facilitated their explosive evolution. Recent research has also shown that craniofacial evolution from radiations in lakes Victoria, Malawi, and Tanganyika has occurred along a shared primary axis of shape divergence, whereby the preorbital region of the skull changes in a manner that is, relatively independent from other head regions. We predicted that the preorbital region would comprise a variational module and used an extensive dataset from each lake that allowed us to test this prediction using a model selection approach. Our findings supported the presence of a preorbital module across all lakes, within each lake, and for Malawi, within sand and rock-dwelling clades. However, while a preorbital module was consistently present, notable differences were also observed among groups. Of particular interest, a negative association between patterns of variational modularity was observed between the sand and rock-dwelling clades, a patter consistent with character displacement. These findings provide the basis for further experimental research involving the determination of the developmental and genetic bases of these patterns of modularity.

1. Introduction

Adaptive divergence is likely influenced by the coordination and integration of multiple traits. If genetic variation affecting patterns of trait covariation have fitness consequences, then a particular pattern of integration that allows for a closer match to a new local multivariate phenotypic optimum should be favoured [13]. Alternatively, ancestrally conserved patterns of integration may act to constrain the rate and direction of evolution by preventing certain functions from evolving [4, 5]. Either way, modularity may influence the pace of evolution and determine evolvability [6, 7]. It is therefore not surprising that the study of trait integration has been of interest to biologists for more than half a century [810] and has recently seen renewed attention [3, 1114].

The study of integration has more recently been extended to the closely related concept of modularity—the relative degrees of connectivity in systems. A module is a tightly integrated unit that is relatively independent from other such modules. For morphological data, modularity has been studied in a variety of contexts including those that are developmental, genetic, functional, and evolutionary in their focus [1522]. An emerging consensus is that patterns of modularity in complex phenotypes likely represents a balance between functional and developmental integration and that modularity is better viewed as a matter of degrees rather than an all-or-nothing phenomenon [12].

It has been suggested that modularity can facilitate divergence by allowing organisms to alter aspects of their phenotype without facing the developmental or fitness tradeoffs that would be present in a wholly integrated unit [12, 13]. In this way, the evolution of modularity could be tied to the idea of key innovations (see [23] for an example). The origin or evolutionary “success” of taxa is often attributed to key innovations—aspects of organismal phenotype that promote diversification [24]. Key innovations may enhance competitive ability, relax adaptive tradeoffs, or permit exploitation of a new productive resource base.

The African cichlids from lakes Victoria, Tanganyika, and Malawi in East-Africa's Great Rift Valley represent the largest extant example of vertebrate adaptive radiations known [25]. Certain anatomical features of this group have been proposed as key innovations that have facilitated the rapid evolution of these fishes [25, 26]. The best known of these represents an important example of functional modularity, wherein the highly derived cichlid pharyngeal jaw mechanism allows the processing of prey within the throat to be decoupled from prey capture by the oral jaws [12, 26]. This is thought to have allowed African cichlids to exploit a wide array of niches that would be unavailable if only one set of jaws was present [27].

The cichlid radiation of Lake Malawi is particularly interesting, because although it is intermediate to Tanganyika and Victoria in terms of age and morphological diversity [28, 34], it has produced the greatest number of endemic species (well over 700) [35, 36]. The evolutionary history of Malawi cichlids suggests that current diversity arose via three stages of diversification: (1) early divergence of the sand-dweller and rock-dweller clades, each of which has adapted to a major macrohabitat, (2) competition for trophic resources within each of these clades that caused further differentiation of trophic morphology, and (3) divergent sexual selection resulting in differentiation of male nuptial coloration [37, 38].

We recently completed an extensive analysis that explored patterns of craniofacial shape variation in African cichlids from each of the three rift lakes [28]. Our data, which represented approximately 80% of the genera across lakes, revealed that all three cichlid radiations share a common trajectory of divergence with respect to each lineage’s major axis of divergence (PC1). Our geometric morphometric analysis also showed that these changes were primarily related to changes in the relative length and size of the “preorbital region” of the skull, which encompasses the oral jaws and supporting structures, with shape posterior to the orbital region remaining relatively stable. These trends suggest that a large portion of the head diversity seen in African cichlids has been achieved by relatively simple and repeated shifts in jaw shape and that these may have happened relatively early in their evolutionary history.

Functional differences in jaw size reflect divergent foraging modes. African cichlids with longer oral jaws are either “suction feeders” and forage on zooplankton, or they are piscivorous and feed on other fishes [28]. Alternatively, species with shortened jaws are typically “biters” that possess a higher mechanical advantage to scrape algae or forage on large macrobenthic prey. In Lake Malawi, this fundamental division is reflected in the cladogenic split between rock- and sand-dwelling species. On average, rock-dwelling species have a shorter jaw in common morphospace, whereas sand-dwellers species have relatively longer jaws [28, Cooper unpublished data]. Notably, these morphological patterns seem to be a common theme in the adaptive radiation of other fish assemblages (e.g., [39]) and even in population-level divergence among ecomorphs of charr whitefish, and sunfish [4042]. Thus, a propensity for changes in the size of oral jaws seems to exist in teleosts at multiple levels of biological organization and perhaps represents a key innovation for this group as a whole. While the evolutionary origins for a preorbital module may not lie within African cichlids examining potential patterns of craniofacial modularity in cichlids may identify important targets for future developmental genetic research to understand the proximate mechanisms that have facilitated these important radiations and divergence in other groups of fishes. Cichlids may be especially useful for this research, because species with widely variable jaw morphologies can be hybridized, facilitating the creation of large populations for genetic mapping to identify the loci and genetic pathways that underlie changes in jaw shape [43, 44].

As mentioned, Liem’s [26] seminal work on the pharyngeal jaw apparatus in cichlids suggested that the functional decoupling of prey capture and processing should free the oral jaws to more readily adopt an array of niche-specific shapes for food capture, largely independent of other traits. Implicitly, this insight confers a level of modularity to the cichlid oral jaw apparatus. Recent work in our lab, as well as from others, supports this assertion by demonstrating that morphological divergence among rift lake cichlids is characterized by prodigious shifts in oral jaw shape [28, 34] and has lead to the specific hypothesis that the preorbital region of the skull represents an evolutionary module that is conserved among cichlids from each of the three East African rift lakes [28]. Here, we objectively test this hypothesis by comparing multiple combinations of models of cichlid head variational modularity. Specifically, we use an approach of model selection recently introduced by Márquez [45] to statistically assess patterns of variational modularity across a large sample of rift valley lake cichlids. To determine whether similar patterns of modularity are operating at different levels of biological organization, we also examine craniofacial modularity in each lake separately, as well as within the rock- and sand-dwelling clades of Lake Malawi.

2. Materials and Methods

2.1. Data Collection

The data used for this study has been previously published in Cooper et al. [28], where further details, including a full list of specimens sampled, can be found. Briefly, our sampling included 78.8% of the genera endemic to the three East-African rift lakes, with the following percentages from each lake: Tanganyika (74.5%), Malawi (88.5%), and Victoria (57.1%). Within Lake Malawi, 19 rock-dwelling species, representing 11 genera, were sampled, and 36 sand-dwelling species, representing 31 genera, were also sampled. Dissections were performed on cichlid heads in order to expose anatomical landmarks important for oral jaw function (Figure 1). A total of sixteen anatomical landmarks were plotted on the images of each specimen using the software program tpsDig2 [46].

Figure 1: Anatomical landmarks placed on the cichlid head. 1 = Tip of the anterior-most tooth on the premaxilla; 2 = Tip of the anterior-most tooth on the dentary; 3 = Maxillary-palatine joint (upper rotation point of the maxilla); 4 = Maxillary-articular joint (lower point of rotation of the maxilla); 5 = Articular-quadrate joint (lower jaw joint); 6 = Insertion of the interopercular ligament on the articular (point at which mouth opening forces are applied); 7 = Posterio-ventral corner of the preopercular; 8 = Most posteriorventral point of the eye socket; 9 = The most anterioventral point of the eye socket; 10 = Joint between the nasal bone and the neurocranium; 11 = Posterior tip of the ascending process of the premaxilla; 12 = Dorsal-most tip of the supraoccipital crest on the neurocranium; 13 = Most dorsal point on the origin of the A1 division of the adductor mandibulae jaw closing muscle on the preopercular; 14 = Most dorsal point on the origin of the A12 division of the adductor mandibulae jaw closing muscle on the preopercular; 15 = Insertion of the A1 division of the adductor mandibulae on the maxilla; 16 = Insertion of the A2 division of the adductor mandibulae on the articular process.
2.2. A Priori Hypotheses

Our goals were to determine first whether modularity was present in the cichlid head and second what the best-supported pattern of modularity was in our data. This required comparative testing of alternative a priori models, each of which specified a particular modular structure in the cichlid head. In this approach, each model is comprised of a series of partitions defined as anatomical regions delimited by landmarks, each representing a hypothesized module predicted to be highly integrated relative to other such partitions.

Based on knowledge of the development and biomechanical function of cichlid heads, we constructed a number of hypotheses of modularity that were intended to extensively cover potential patterns of covariance. We selected a total of five a priori models representing the spatial distribution of developmental units and functional components of the cichlid head (see Table 1). An additional “null” model representing a lack of any integration or modularity was included in our analyses. Because it is not biologically realistic to expect that patterns of modularity predicted by these developmental and functional models are mutually exclusive, all possible nonnested combinations of the modules defined by the original five hypotheses of modularity were also included in model comparisons. In total, 137 competing models were tested. It is important to note that while this list of hypotheses is far from exhaustive, it represents an extensive collection of models—likely covers a substantial proportion of the developmental and functional processes capable of affecting covariation in the cichlid head.

Table 1: A priori hypotheses of modularity in the cichlid head. Brackets denote putative modules. Note that two similar models are presented for jaw function.
2.3. Modularity Analysis

The methodology for testing a priori hypotheses of modularity was adapted from an approach proposed by Márquez [45] consisting of four basic steps implemented in the Mint software package (available at: http://www-personal.umich.edu/~emarquez/morph/).

(1) Computation of an expected covariance matrix from each model of modularity, by assuming that each module resides in its own subspace within the phenotypic space occupied by the entire structure, as described in Márquez [45].

(2) Computation of a goodness of fit statistic, , to measure the dissimilarity between observed and expected covariances for each model, as where S and S0 are the observed and modeled covariance matrices and is the transpose symbol [47]. To ensure the comparability of this statistic across models, is standardized twice: first, all values are divided by , corresponding to the null model describing complete absence of integration, so that is scaled to vary within the interval ; second, scaled is standardized via linear regression to remove the effect of the number of estimated parameters in models, which takes advantage of the linear relation observed between and the number of zeros in models.

The standardized statistic is defined as the residual m , where represents the linear function relating the values of computed from all possible models of modularity to their corresponding counts of zero elements, . Even though it would be computationally unfeasible for most studies to include all possible models, the fact that scaled values are restricted to the interval [0, 1], where 0 corresponds to the observed covariance matrix and 1 to the null model of no integration, implies that must also vary within these limits, which are sufficient to define the linear function for any given set of variables. Given a large random sample of models, with values symmetrically distributed about their mean, , and thus . Consequently, models in which correspond to comparisons where observed covariances are relatively low on average and hypothesized to be zero, and conversely, cases where occur when relatively high covariances are on average hypothesized to be zero, the best-fitting model is that with the lowest value. Note that this approach differs slightly from the one used in Márquez [45], where was estimated via regression using only the models included in a study, as opposed to all possible models. 95% confidence intervals were computed as the 2.5 and 97.5 percentiles of a distribution of 1,000 jackknife subsamples [48] formed by removing random subsets of 10% of the specimens from each sample.

(3) The statistical significance of was assessed using a parametric Monte Carlo approach. In these tests, a null distribution for the statistic is generated by comparing the original observed covariance matrix S to each of 1,000 random matrices generated from a Wishart distribution with mean vector 0 (i.e., the same mean as Procrustes residual data) and covariance matrix S [45, 49].

(4) Finally, to allow choice among the multiple models that are significantly better than chance according to the Monte Carlo approach described above, models are ranked by their goodness of fit (i.e., values, in ascending order). The relative support for each model is determined by computing the stability of its rank using a jackknife approach in which values and model ranks are recomputed after removing a random portion of the samples. In this study, we removed 10% of the data in each of 1,000 jackknife replicates.

2.4. Comparisons of Covariance among Lakes, Sand Dwellers, and Rock Dwellers

If a single model fits two of our groups (Lake Tanganyika, LT; Lake Malawi, LM; Lake Victoria, LV; Rock dwellers, RD; Sand dwellers, SD) equally well, it would not necessarily mean that they were close to each other in our model space. This is because two objects that are equally distant from a third (the best supported model) are not required to occupy the same position, especially in a high-dimensional space. In our case, the values calculated for each group represents reference points useful for determining their relative position. This vector of values can have two interpretations, the first as a set of distances between the observed covariation matrix and known patterns of modularity and the second as coordinates for the data in “model space” centered on a group covariance pattern. Because each group may be centered at a different position, only the direction of these vectors can be compared, which was achieved through the use of correlations between vectors of each group. This involved two separate analysis; first, we determined levels of correlation for across the three lakes; second, we determined levels of correlation for among LM, RD, and SD groups. However, we did not use all 137 possible gamma values in these correlations; rather, we used the ten top-ranked models in for each group. This increased the likelihood that we were testing associations between the most biologically relevant models.

3. Results

Monte Carlo tests were unable to distinguish among models, suggesting that hypotheses were too similar to distinguish amongst each other given available sample sizes. We, therefore, focus our interpretations on the basis of the relative rankings of values and their jackknife support.

Overall, there was strong support for the hypothesis that modularity is present in the heads of African cichlids. Across the three lakes the null model of no integration was ranked 57th, 100th, and 102nd, out the 137 models in LV, LM, and LT, respectively. In the RD and SD groups the null model of no integration ranked 70th and 108th, respectively. Jackknife tests provided high support for these rankings in all groups.

3.1. Top-Ranked Models

At all levels, the best supported hypothesis included one preorbital and one postorbital module. In our pooled data set across all lakes, as well as separate data sets for LV and RD, a preorbital module that defined the upper jaws and encompassed the exact same set of landmarks was identified (Figure 2). Support for these patterns of modularity was high with the top model in the pooled sample of cichlids being ranked number 1 in 96.6% of jackknife reps. LV and RD groupings had top models that were similarly highly supported with 84%, and 85% of jackknife reps, respectively.

Figure 2: The best supported hypotheses of modularity for African cichlids as a whole (a) and their respective adaptive radiations within lakes Malawi (b), Tanganyika (c), and Victoria (d). Also shown are the best-supported patterns for rock (e), and sand dwellers (f, g) within Malawi. Two models are shown for rock dwellers, because our statistical analysis was unable to discern whether one hypotheses was significantly better at describing patterns of covariance. Note that despite differences among these hypothesis, each contains a preorbital module based in the oral jaws.

For LM as a whole and the SD sample, the top ranked models displayed a preorbital module that encompassed both the upper and lower jaws (Figure 2). Statistical support for the LM model was high with 86% of jackknife reps maintaining its top ranking. In the case of the SD sample, there were two, statistically indistinguishable top models: The highest-ranked model included three modules, one encompassing the oral jaws, one defining the orbital size, and another that covers much of the posterior region of the head. The second ranked model was identical to the first with the exception that it did not possess an eye/orbital module. Support for the best SD model (i.e., three modules) was low, with only 47% of jackknife reps supporting its ranking. The second best SD model (i.e., two modules) was also ranked as the best model in 44% of jackknife reps. However, a subsequent set of analyses found that when one of these models was removed, support for the other model significantly improved to where its top ranking was supported in over 97% of jackknife reps. This analysis suggests that both models are equally valid.

The LT dataset also showed strong support for a preorbital module in its top-ranked hypothesis (supported in 98.6% of jackknife reps). However, it differed from the other groups by having a preorbital module comprised primarily of the lower jaw (Figure 2).

3.2. Relationships among Patterns of Covariance Across Lakes, and between the Lake Malawi Sand and Rock Dwellers

Across the three lakes, we observed strikingly similar patterns of covariation. We used values from a total of 23 hypotheses of modularity, reflecting the top ten ranked models for each of the three lakes, meaning that 7 of these hypotheses were shared among lakes. The -values for our tests were all extremely high, and positively correlated, indicating that in spite of differences between top-ranked models, very similar patterns of covariance underlie each of these adaptive radiations (Table 2).

Table 2: Correlation coefficients of values for the top ten ranked hypotheses of modularity between cichlids in Malawi, Tanganyika, and Victoria.

We also found that patterns of covariance may be diverging within LM. We used a total of 18 models to describe the top ten models across LM and within the RD and SD datasets. Thus, a total of 12 out of a possible 30 models were shared among these groups. The correlation between values for LM as a whole and SD dataset was particularly strong, indicating that sand dwellers may be influencing the overall pattern of modularity exhibited by Malawi cichlids. This result could be due, in part, to their larger relative sample size compared to RD cichlids. Alternatively, LM as a whole showed almost no relationship with RD species, and there was a strong negative relationship between SD and RD species (Table 3). These data suggest that patterns of trait covariance are being repelled between SD and RD.

Table 3: Correlation coefficients of values for the top ten ranked hypotheses of modularity between cichlids in Malawi as a whole, as well as the rock- and sand-dweller division within the lake.

4. Discussion

Our results demonstrate that a preorbital module is present in the oral jaws of East African rift valley cichlids and that this pattern of covariation is conserved across all lakes. This trend strongly supports the hypothesis that this pattern of modularity has influenced the rate and direction of adaptive phenotypic divergence among African cichlid radiations—an idea rooted in the proposal that the cichlid pharyngeal jaw apparatus is a key innovation that freed the oral jaws from a functional constraint [26], formalized in light of quantitative patterns of trophic divergence among cichlid lineages [28], and empirically tested here. While our results are compelling, we suggest that the comparisons of rates of evolution to other groups which lack a pharyngeal jaw apparatus (e.g., salmonids and characids), and possibly a preorbital module, may be needed to confirm whether the patterns of modularity identified in cichlids represent a key innovation.

4.1. Conserved and Divergent Patterns of Craniofacial Modularity among Cichlids

While the results of our correlation analyses indicate that general patterns of covariance are conserved across lakes, there were several notable differences in the top-ranked hypotheses of modularity among groups, suggesting that while conserved patterns exist, modularity itself is capable of evolving. The LM dataset had a pattern of modularity in which the preorbital module encompassed both the upper and lower jaws, while in the LT dataset, the preorbital module was exclusive to the mandible. For the LV dataset, integration was most prevalent for the upper jaws. The LV radiation is the youngest of the three rift lakes [50], and correspondingly, our prior analysis found relatively low levels of shape variation (disparity) in this lake compared to LM and LT cichlids [28]. Also, more than 60% of the morphological variation among species in LV can be explained by a single major axis (principal component), considerably more than was explained by this shared axis for LM and LT (i.e., Victoria cichlid head anatomies were relatively more integrated). Taken together, these results suggest that the younger divergence in LV is determined by a more limited set of strong interactions among traits. Since the upper jaw contains the anatomical linkages most responsible for highly kinetic jaw movements, such as jaw protrusion, this would imply that both the functional and morphological evolution of this lineage has been constrained. As the youngest of the three rift lake lineages, patterns observed within LV may offer insight into the proximate mechanisms that have shaped cichlid radiations in general. It is possible that the pattern of modularity we have identified in LV has played a dominant role in the early patterns of divergence of cichlids in LM and LT. Consistent with this idea, the preorbital module identified in the upper jaw for LV was very similar to one identified in the top-ranked model for our pooled data set across all cichlids (Figure 2).

The top-ranked models for the SD and RD clades within LM also exhibited notable differences. Whereas the SD group exhibited a preorbital module that included both the upper and lower jaws, RD species expressed a pattern of modularity similar to that of LV, where only the upper jaws were integrated. Moreover, the SD/RD division within Malawi was characterized by a strong negative relationship in covariance patterns, suggesting that ecological competition between these clades during the early history of the lake may have caused patterns of trait covariance to diverge. This pattern is consistent with character displacement, but at a different biological scale (groups of species or clades) than where it is usually recognized [5153]. Character displacement is often thought to occur between two closely related species; however, research suggests that character displacement can also occur between distantly related species, as well as whole communities [54, 55]; see also [56, 57] for evidence of character displacement in African cichlids. Therefore, it is appropriate to speculate that this process is contributing to divergence between SD and RD clades in LM.

Integration between the upper and lower jaws, as displayed by the SD dataset, may be especially advantageous for ram/suction-feeding predators, a predominant SD trophic niche [58], because both jaws need to work together in a highly coordinated fashion to produce kinematic force [59]. Alternatively, in RD species that most often employ a biting tactic whereby the upper jaw is relatively more stationary during foraging [26], the upper jaw is integrated, and the lower is not. This implies that the lower jaw in RD species is free to evolve a wide array of geometries, which may be advantageous for substrate feeding species, where demands on the lower jaw should be more variable relative to the upper jaw apparatus. However, this is not to say that there is a complete lack of integration between the upper and lower jaw, as modularity is a matter of degree rather than an all-or-nothing phenomenon [12] . Also, it is important to note that patterns of divergence among SD and RD are still acting within the overall context of a preorbital module (i.e., both upper and lower jaw for SD, upper jaw for RD,) suggesting that the rate and direction of phenotypic evolution is being dictated by historical constraints that are manifested in patterns of covariance and modularity. In other words, putative character displacement between SD and RD species in Malawi cichlids may be proceeding along genetic lines of least resistance [5, 53].

4.2. Origins for Adult Patterns of Modularity: Developmental Mechanisms

Although there are a number of possible functional explanations for patterns of craniofacial modularity, it is important to remember that selection must work within the context of developmental systems to improve functional performance. That is not to say development inherently constrains evolution, but rather that it can direct its outcome in concert with selection. In fact, simulations have shown that some degree of order may actually be required for evolution to proceed with ease [60]. It is, therefore, probable that the patterns of craniofacial modularity identified here, while probably causing an increased propensity for adaptations involving the oral jaws, are also dictated by underlying developmental processes. Clues to these potential processes may lie in early embryological events during the formation of craniofacial anatomy in fishes (see [61] for a similar view in mammals).

Structural progenitors of the ossified structures in the preorbital region of the skull include the trabeculae and ethmoid cartilages (i.e., anterior neurocranium), palatoquadrate (i.e., upper jaw precursor), and Meckel's cartilage (i.e., lower jaw precursor). All of these structures are derived from the same population of anterior cranial neural crest (CNC) cells that migrate away from neural tissue beginning at approximately 12 hours afterfertilization (hpf) in zebrafish [62]. Thus, the preorbital region of the skull is defined early in development, and these events may underlie the persistence of a preorbital module among African riftlake cichlids. For instance, LM cichlids show integration between the upper and lower jaws, suggesting that this developmental hypothesis may have particular merit for this adaptive radiation.

The modular divisions between the upper and lower jaws found between LV and LT may be influenced by slightly later developmental events. Fate mapping experiments in zebrafish show that at approximately 24 hpf the stomodeum forms as an invagination of the oral ectoderm, and both the pterygoid process and anterior neurocranium reside within a compact condensation of cells closely associated with dorsal edge of this structure, whereas Meckel’s cartilage forms from cells ventral to this structure [62]. Thus, while early ontogenetic events (i.e., CNC migration) regionalize the skull along the anterior-posterior axis, slightly later events (i.e., formation of the mouth) are necessary to specify the dorsal-ventral identity of the jaws within the preorbital region of the skull.

Later still in development, the sequence of ossification in bones of the craniofacial region may play a role in determining patterns of modularity. Evidence from zebrafish and Nile tilapia (Oreochromis niloticus) show that the oral jaws (premaxillae, maxillae, and dentary) are among the first structures to become mineralized in the teleost head [63, 64]. In fact, the only other structures that are ossified as early as the oral jaws include the basio-occipital and opercle. Functional reasons have been attributed to this chronological pattern in teleost development [63, 6567]. Specifically, bones involved in early basic functions such as respiration and feeding have been observed to ossify first. This suggests that the bones of the oral jaws and opercular regions of the skull are predisposed to reflect the patterns of variational modularity we have identified. Ossification sequence, and heterochronic shifts in this process, could, therefore, act as another early mechanism that sets the stage for craniofacial modularity throughout life history.

4.3. Origins for Adult Patterns of Modularity: Integrating Developmental and Functional Processes

Beyond initial ossification, bone remodeling over ontogeny could represent another means of achieving modularity of the oral jaws and a way of simultaneously integrating developmental and functional mechanisms in a straightforward way. Bone is a dynamic, metabolically active tissue that is constantly being renewed and changed. Bone cells are strain sensitive and can transduce signals from mechanical loading into cues that result in either reduced bone loss or gain [6871]. Disuse usually causes an acceleration of bone turnover, with resorption being the dominant process. Conversely, excessive strain can damage bone, which may in turn be repaired or further reinforced through remodeling. Importantly, both bone resorption and deposition involve highly conserved genetic and developmental pathways [7274].

Mechanical stimuli may be particularly important for inducing adaptive patterns of modularity through the process of bone remodelling. Bone turnover tends to be most effective in areas of high stress, thus reducing the risk of injury [69]. In teleosts, the oral jaws are used for both respiration and food acquisition, but it is likely that the oral jaws are under the highest stress during food acquisition and processing, which should in turn provide the greatest stimulus for bone remodeling [72]. Indeed, several lab-based studies on cichlids have documented that different diet treatments can induce changes in bone and head shape [75, 76], demonstrating the ability of elements in the upper and lower jaws to respond to mechanical stimuli through changes in shape. Within the RD lineage of LM, it is certainly possible that a high degree of remodeling and plasticity of the lower jaw has led to a pattern of modularity, wherein the mandible lacks a measurable degree of integration across species. The lower jaw may be more amenable to remodeling due the greater degree of movement that it is afforded in the RD lineage. Alternatively, patterns of integration within the lower jaw may differ between species, resulting in a perceived lack of integration in the combined dataset. In either case, the conclusion that must be drawn is that the lower jaw is a highly evolvable trait within the RD lineage.

Perhaps the most compelling evidence for a fundamental link between developmental and functional processes comes from work in the BMP family of signaling proteins (reviewed by [70]). Critical roles for BMP signaling during bone and cartilage development are well established (reviewed by [77]), and variation in BMP expression over ontogeny has been associated with the origin and adaptation of key vertebrate innovations including the turtle shell [78], bat wing [79], cichlid mandible [44, 80], and bird bills [8183]. All of these examples involve differential Bmp expression that is presumed to be due to mutational effects (either cis or trans), but several studies have also documented environmentally induced changes in BMP expression in skeletal tissue. Specifically, tensile stress has been shown to alter BMP expression during bone growth [84, 85], remodeling [86], and repair [87, 88]. Thus, a scenario wherein patterns of craniofacial modularity are established via early developmental mechanisms and then either reinforced or altered by functional processes might represent the true nature of variational modules within the cichlid skull. Examining how patterns of integration potentially shift over ontogeny and under different feeding regimes in different cichlid lineages would represent a fruitful line of future research.

4.4. Modularity and Evolvability of the Craniofacial Skeleton in Cichlids

Recent reviews suggest that an extended evolutionary synthesis (EES) is necessary to account for the origins of variation that is acted upon by natural selection [6, 7]. The empirical center for the EES will lie in discovering the features of organisms that determine evolvability [7]. While specific definitions of evolvability are numerous and vary according to context, modularity figures prominently in these discussions insofar as it imposes a constraint on direction or speed of evolutionary change [12, 13]. In this context, we suggest that modularity can act as a “key innovation”. While key innovations are typically defined by the appearance of an anatomical structure that precedes an adaptive radiation, as is the case for the pharyngeal jaws [26], we contend that patterns of modularity, whereby the cichlid oral jaws represent a module that allows them to change with a high degree of autonomy, have had a strong influence on the rate and direction of adaptive divergence in this group. This pattern of modularity is likely what has allowed for the rapid lengthening or shortening of the oral jaws relative to the rest of the head in cichlids and shape changes that comprise the major axes of variation in each of the three African rift lakes, and likely, it represents the template upon which additional changes in trophic morphology occur [28]. In other words, the evolution of this pattern of modularity may facilitate evolution, providing an example of the “evolution of evolvability” (see [7]). The degree to which these patterns are specific to cichlids, or represent a more generalized perciform innovation, will be an important area of future study. Several avenues may have lead to preorbital modularity; therefore, finding groups that lack this pattern of modularity and comparing rates of diversification will be important for identifying its potential role as a key innovation.

As discussed above, several avenues may have led to the consistent patterns of preorbital modularity we have discovered. In the order of their ontogenetic appearance, these include (1) migration and specification of progenitor cells, (2) dorsal-ventral division of the oral cavity, (3) sequence of ossification with early calcification of the jaws and operculum region, and (4) remodeling of bone in response to mechanical stimuli. These all represent separate hypotheses and processes that can be tested to understand the developmental and genetic basis of a preorbital module. We predict that each of these processes may play important roles in determining modularity in the cichlid head, depending on the lineage being queried. Fortunately, we have the means to assess patterns of modularity over ontogeny in cichlids and can statistically track when the patterns we have identified in adult cichlids begin to emerge. We also have the means to identify QTL associated with these anatomical modules and to track changes in gene expression during the emergence of these patterns [43, 44]. In all, cichlids represent an attractive model to reveal both the genetic basis of modularity and the evolvability of the craniofacial skeleton.


We would like to thank the curators and staff of the following museums for their invaluable assistance: the American Museum of Natural History; the Belgian Royal Museum for Central Africa; Cornell University's Museum of Vertebrates; Harvard University's Museum of Comparative Zoology; and the University of Michigan's Museum of Zoology. We would also like to thank anonymous reviewers for their commentary on manuscript drafts.


  1. E. Mayr, “Change of genetic environment and evolution,” in Evolution as a Process, J. Huxley, A. C. Hardy, and E. B. Ford, Eds., pp. 157–180, Allen and Unwin, London, UK, 1954. View at Google Scholar
  2. G. Bell, “On the function of flowers,” Proceedings of the Royal Society of London B, vol. 224, no. 1235, pp. 223–265, 1985. View at Google Scholar · View at Scopus
  3. K. J. Parsons and B. W. Robinson, “Replicated evolution of integrated plastic responses during early adaptive divergence,” Evolution, vol. 60, no. 4, pp. 801–813, 2006. View at Publisher · View at Google Scholar · View at Scopus
  4. G. Marroig and J. M. Cheverud, “Size as a line of least evolutionary resistance: diet and adaptive morphological radiation in New World monkeys,” Evolution, vol. 59, no. 5, pp. 1128–1142, 2005. View at Google Scholar · View at Scopus
  5. D. Schluter, “Adaptive radiation along genetic lines of least resistance,” Evolution, vol. 50, no. 5, pp. 1766–1774, 1996. View at Google Scholar · View at Scopus
  6. J. L. Hendrikse, T. E. Parsons, and B. Hallgrímsson, “Evolvability as the proper focus of evolutionary developmental biology,” Evolution and Development, vol. 9, no. 4, pp. 393–401, 2007. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  7. M. Pigliucci, “Is evolvability evolvable?” Nature Reviews Genetics, vol. 9, no. 1, pp. 75–82, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  8. E. C. Olson and R. L. Miller, Morphological Integration, University of Chicago Press, Chicago, Ill, USA, 1958.
  9. R. L. Berg, “The ecological significance of correlation pleiades,” Evolution, vol. 14, pp. 171–180, 1960. View at Google Scholar
  10. L. Van Valen, “The study of morphological integration,” Evolution, vol. 19, pp. 347–349, 1965. View at Google Scholar
  11. M. Pigliucci and K. Preston, Phenotypic Integration: Studying the Ecology and Evolution of Complex Phenotypes, Oxford University Press, Oxford, UK, 2004.
  12. C. P. Klingenberg, “Morphological integration and developmental modularity,” Annual Review of Ecology, Evolution, and Systematics, vol. 39, pp. 115–132, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. C. P. Klingenberg, “Evolution and development of shape: integrating quantitative approaches,” Nature Reviews Genetics, vol. 11, no. 9, pp. 623–635, 2010. View at Publisher · View at Google Scholar · View at PubMed
  14. M. Pavlicev, J. M. Cheverud, and G. P. Wagner, “Measuring morphological integration using eigenvalue variance,” Evolutionary Biology, vol. 36, no. 1, pp. 157–170, 2009. View at Publisher · View at Google Scholar · View at Scopus
  15. J. M. Cheverud, E. J. Routman, and D. J. Irschick, “Pleiotropic effects of individual gene loci on mandibular morphology,” Evolution, vol. 51, no. 6, pp. 2006–2016, 1997. View at Google Scholar · View at Scopus
  16. C. P. Klingenberg, A. V. Badyaev, S. M. Sowry, and N. J. Beckwith, “Inferring developmental modularity from morphological integration: analysis of individual variation and asymmetry in bumblebee wings,” American Naturalist, vol. 157, no. 1, pp. 11–23, 2001. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  17. C. P. Klingenberg, K. Mebus, and J. C. Auffray, “Developmental integration in a complex morphological structure: how distinct are the modules in the mouse mandible?” Evolution and Development, vol. 5, no. 5, pp. 522–531, 2003. View at Publisher · View at Google Scholar · View at Scopus
  18. C. P. Klingenberg, L. J. Leamy, and J. M. Cheverud, “Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible,” Genetics, vol. 166, no. 4, pp. 1909–1921, 2004. View at Publisher · View at Google Scholar · View at Scopus
  19. G. P. Wagner and J. G. Mezey, “The role of genetic architecture constraints in the origin of variational modularity,” in Modularity in Development and Evolution, G. Schlosser and G. P. Wagner, Eds., pp. 338–358, The University of Chicago Press, Chicago, Ill, USA, 2004. View at Google Scholar
  20. E. J. Márquez, “A statistical framework for testing modularity in multidimensional data,” Evolution, vol. 62, no. 10, pp. 2688–2708, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  21. M. L. Zelditch, A. R. Wood, R. M. Bonett, and D. L. Swiderski, “Modularity of the rodent mandible: integrating bones, muscles, and teeth,” Evolution and Development, vol. 10, no. 6, pp. 756–768, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  22. M. L. Zelditch, A. R. Wood, and D. L. Swiderski, “Building developmental integration into functional systems: function-induced integration of mandibular shape,” Evolutionary Biology, vol. 36, no. 1, pp. 71–87, 2009. View at Publisher · View at Google Scholar · View at Scopus
  23. A. S. Yang, “Modularity, evolvability, and adaptive radiations: a comparison of the hemi- and holometabolous insects,” Evolution and Development, vol. 3, no. 2, pp. 59–72, 2001. View at Publisher · View at Google Scholar · View at Scopus
  24. J. P. Hunter, “Key innovations and the ecology of macroevolution,” Trends in Ecology and Evolution, vol. 13, no. 1, pp. 31–36, 1998. View at Publisher · View at Google Scholar · View at Scopus
  25. W. Salzburger, “The interaction of sexually and naturally selected traits in the adaptive radiations of cichlid fishes,” Molecular Ecology, vol. 18, no. 2, pp. 169–185, 2009. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  26. K. F. Liem, “Evolutionary strategies and morphological innovations: cichlid pharyngeal jaws,” Systematic Zoology, vol. 22, pp. 425–441, 1973. View at Google Scholar
  27. C. D. Hulsey and J. T. Streelman, “Cichlid pharyngeal jaw fusion: is suturing the key to unparalleled trophic divergence?” Integrative and Comparative Biology, vol. 45, p. 1016, 2005. View at Google Scholar
  28. W. J. Cooper, K. Parsons, A. McIntyre, B. Kern, A. McGee-Moore, and R. C. Albertson, “Bentho-pelagic divergence of cichlid feeding architecture was prodigious and consistent during multiple adaptive radiations within African Rift-Lakes,” PLoS ONE, vol. 5, no. 3, Article ID e9551, 2010. View at Publisher · View at Google Scholar · View at PubMed
  29. G. C. Anker, “Morphology and kinematics of the stickleback, Gasterosteus aculeatus,” Transactions of the Zoological Society (London), vol. 32, pp. 311–416, 1974. View at Google Scholar
  30. M. Muller, “Optimization principles applied to the mechanism of neurocranium levation and mouth bottom depression in bony fishes (Halecostomi),” Journal of Theoretical Biology, vol. 126, no. 3, pp. 343–368, 1987. View at Google Scholar
  31. M. W. Westneat, “Kinematics of jaw protrusion in cheiline wrasses (Labridae)—testing a 4-bar linkage model,” American Zoologist, vol. 28, p. A114, 1988. View at Google Scholar
  32. M. W. Westneat, “Feeding mechanics of teleost fishes (Labridae; Perciformes): a test of four-bar linkage models,” Journal of Morphology, vol. 205, no. 3, pp. 269–295, 1990. View at Publisher · View at Google Scholar
  33. M. W. Westneat, “A biomechanical model for analysis of muscle force, power output and lower jaw motion in fishes,” Journal of Theoretical Biology, vol. 223, no. 3, pp. 269–281, 2003. View at Publisher · View at Google Scholar · View at MathSciNet
  34. K. A. Young, J. Snoeks, and O. Seehausen, “Morphological diversity and the roles of contingency, chance and determinism in African cichlid radiations,” PLoS One, vol. 4, no. 3, Article ID e4740, 2009. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  35. J. Snoeks, “Cichlid diversity, speciation and systematics: examples from the Great African Lakes,” Journal of Aquariculture and Aquatic Sciences, vol. 9, pp. 150–166, 2001. View at Google Scholar
  36. G. F. Turner, O. Seehausen, M. E. Knight, C. J. Allender, and R. L. Robinson, “How many species of cichlid fishes are there in African lakes?” Molecular Ecology, vol. 10, no. 3, pp. 793–806, 2001. View at Publisher · View at Google Scholar
  37. P. D. Danley and T. D. Kocher, “Speciation in rapidly diverging systems: lessons from Lake Malawi,” Molecular Ecology, vol. 10, no. 5, pp. 1075–1086, 2001. View at Publisher · View at Google Scholar
  38. O. Seehausen, “Hybridization and adaptive radiation,” Trends in Ecology and Evolution, vol. 19, no. 4, pp. 198–207, 2004. View at Publisher · View at Google Scholar · View at PubMed
  39. W. J. Cooper and M. W. Westneat, “Form and function of damselfish skulls: rapid and repeated evolution into a limited number of trophic niches,” BMC Evolutionary Biology, vol. 9, no. 1, Article ID 24, 2009. View at Publisher · View at Google Scholar · View at PubMed
  40. S. Skulason, D. L. G. Noakes, and S. S. Snorrason, “Ontogeny of trophic morphology in four sympatric morphs of arctic charr Salvelinus alpinus in Thingvallavatn, Iceland,” Biological Journal of the Linnean Society, vol. 38, no. 3, pp. 281–301, 1989. View at Google Scholar
  41. B. W. Robinson, D. S. Wilson, A. S. Margosian, and P. T. Lotito, “Ecological and morphological differentiation of pumpkinseed sunfish in lakes without bluegill sunfish,” Evolutionary Ecology, vol. 7, no. 5, pp. 451–464, 1993. View at Publisher · View at Google Scholar
  42. K. Østbye, T. F. Næsje, L. Bernatchez, O. T. Sandlund, and K. Hindar, “Morphological divergence and origin of sympatric populations of European whitefish (Coregonus lavaretus L.) in Lake Femund, Norway,” Journal of Evolutionary Biology, vol. 18, no. 3, pp. 683–702, 2005. View at Publisher · View at Google Scholar · View at PubMed
  43. R. C. Albertson, J. T. Streelman, and T. D. Kocher, “Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 9, pp. 5252–5257, 2003. View at Publisher · View at Google Scholar · View at PubMed
  44. R. C. Albertson, J. T. Streelman, T. D. Kocher, and P. C. Yelick, “Integration and evolution of the cichlid mandible: the molecular basis of alternate feeding strategies,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 45, pp. 16287–16292, 2005. View at Publisher · View at Google Scholar · View at PubMed
  45. E. J. Márquez, “A statistical framework for testing modularity in multidimensional data,” Evolution, vol. 62, no. 10, pp. 2688–2708, 2008. View at Publisher · View at Google Scholar · View at PubMed
  46. F. J. Rohlf, “tpsDig2,” A software used to collect xy coordinate data for use in geometric morphometrics, 2009, http://life.bio.sunysb.edu/morph/.
  47. J. T. Richtsmeier, S. R. Lele, and T. M. Cole III, “Landmark morphometrics and the analysis of variation,” in Variation: A Central Concept in Biology, B. Hallgrímsson and B. K. Hall, Eds., pp. 49–69, Elsevier, Amsterdam, The Netherlands, 2005. View at Google Scholar
  48. B. F. J. Manly, Randomization, Bootstrap and Monte Carlo Methods in Biology, Chapman and Hall/CRC Press, Boca Raton, Fla, USA, 3rd edition, 2006.
  49. W. J. Krzanowski, Principles of Multivariate Analysis: A User's Perspective, Oxford University Press, Oxford, UK, 2000.
  50. K. R. Elmer, C. Reggio, T. Wirth, E. Verheyen, W. Salzburger, and A. Meyer, “Pleistocene desiccation in East Africa bottlenecked but did not extirpate the adaptive radiation of Lake Victoria haplochromine cichlid fishes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 32, pp. 13404–13409, 2009. View at Publisher · View at Google Scholar · View at PubMed
  51. W. L. Brown and E. O. Wilson, “Character displacement,” Systematic Zoology, vol. 5, pp. 49–64, 1956. View at Google Scholar
  52. D. Schluter and J. D. McPhail, “Ecological character displacement and speciation in sticklebacks,” American Naturalist, vol. 140, no. 1, pp. 85–108, 1992. View at Publisher · View at Google Scholar · View at PubMed
  53. D. Schluter, The Ecology of Adaptive Radiation, Oxford University Press, New York, NY, USA, 2000.
  54. D. Schluter, “Character displacement between distantly related taxa? Finches and bees in the Galapagos,” American Naturalist, vol. 127, no. 1, pp. 95–102, 1986. View at Google Scholar
  55. T. Dayan and D. Simberloff, “Ecological and community-wide character displacement: the next generation,” Ecology Letters, vol. 8, no. 8, pp. 875–894, 2005. View at Publisher · View at Google Scholar
  56. L. Ruber and D. C. Adams, “Evolutionary convergence of body shape and trophic morphology in cichlids from Lake Tanganyika,” Journal of Evolutionary Biology, vol. 14, no. 2, pp. 325–332, 2001. View at Publisher · View at Google Scholar
  57. C. J. Allender, O. Seehausen, M. E. Knight, G. F. Turner, and N. Maclean, “Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 2, pp. 14074–14079, 2003. View at Publisher · View at Google Scholar · View at PubMed
  58. A. Konings, Malawi Cichlids in Their Natural Habitat, Cichlid Press, Lauenau, Germany, 2007.
  59. M. W. Westneat, “Skull biomechanics and suction feeding in fishes,” in Fish Biomechanics, G. V. Lauder and R. E. Shadwick, Eds., Academic Press, 2005. View at Google Scholar
  60. G. P. Wagner, M. Pavlicev, and J. M. Cheverud, “The road to modularity,” Nature Reviews Genetics, vol. 8, no. 12, pp. 921–931, 2007. View at Publisher · View at Google Scholar · View at PubMed
  61. W. R. Atchley and B. K. Hall, “A model for development and evolution of complex morphological structures,” Biological Reviews of the Cambridge Philosophical Society, vol. 66, no. 2, pp. 101–157, 1991. View at Google Scholar
  62. J. K. Eberhart, M. E. Swartz, J. G. Crump, and C. B. Kimmel, “Early Hedgehog signaling from neural to oral epithelium organizes anterior craniofacial development,” Development, vol. 133, no. 6, pp. 1069–1077, 2006. View at Publisher · View at Google Scholar · View at PubMed
  63. C. C. Cubbage and P. M. Mabee, “Development of the cranium and paired fins in the zebrafish Danio rerio (Ostariophysi, Cyprinidae),” Journal of Morphology, vol. 229, no. 2, pp. 121–160, 1996. View at Google Scholar
  64. K. Fujimura and N. Okada, “Shaping of the lower jaw bone during growth of Nile tilapia Oreochromis niloticus and a Lake Victoria cichlid Haplochromis chilotes: a geometric morphometric approach,” Development Growth and Differentiation, vol. 50, no. 8, pp. 653–663, 2008. View at Publisher · View at Google Scholar · View at PubMed
  65. P. Vandewalle, B. Focant, F. Huriaux, and M. Chardon, “Early development of the cephalic skeleton of Barbus barbus (Teleostei, Cyprinidae),” Journal of Fish Biology, vol. 41, pp. 43–62, 1992. View at Google Scholar
  66. P. M. Mabee and T. A. Trendler, “Development of the cranium and paired fins in Betta splendens (Teleostei: Percomorpha): intraspecific variation and interspecific comparisons,” Journal of Morphology, vol. 227, no. 3, pp. 249–287, 1996. View at Google Scholar
  67. P. Le Pabic, E. J. Stellwag, and J. L. Scemama, “Embryonic development and skeletogenesis of the pharyngeal jaw apparatus in the cichlid Nile tilapia (Oreochromis niloticus),” Anatomical Record, vol. 292, no. 11, pp. 1780–1800, 2009. View at Publisher · View at Google Scholar · View at PubMed
  68. S. W. Herring, “Epigenetic and functional influences on skull growth,” in The Skull, J. Hanken and B. K. Hall, Eds., pp. 153–206, University of Chicago Press, Chicago, Ill, USA, 1993. View at Google Scholar
  69. A. G. Robling, A. B. Castillo, and C. H. Turner, “Biomechanical and molecular regulation of bone remodeling,” Annual Review of Biomedical Engineering, vol. 8, pp. 455–498, 2006. View at Publisher · View at Google Scholar · View at PubMed
  70. R. L. Young and A. V. Badyaev, “Evolution of ontogeny: linking epigenetic remodeling and genetic adaptation in skeletal structures,” Integrative and Comparative Biology, vol. 47, no. 2, pp. 234–244, 2007. View at Publisher · View at Google Scholar
  71. K. J. Parsons and R. C. Albertson, “Roles for Bmp4 and CaM1 in shaping the jaw: evo-devo and beyond,” Annual Review of Genetics, vol. 43, pp. 369–388, 2009. View at Publisher · View at Google Scholar · View at PubMed
  72. R. C. Albertson and P. C. Yelick, “Fgf8 haploinsufficiency results in distinct craniofacial defects in adult zebrafish,” Developmental Biology, vol. 306, no. 2, pp. 505–515, 2007. View at Publisher · View at Google Scholar · View at PubMed
  73. D. M. Parichy, D. G. Ransom, B. Paw, L. I. Zon, and S. L. Johnson, “An orthologue of the kit-related gene fms is required for development of neural crest-derived xanthophores and a subpopulation of adult melanocytes in the zebrafish, Danio rerio,” Development, vol. 127, no. 14, pp. 3031–3044, 2000. View at Google Scholar
  74. P. E. Witten, A. Hansen, and B. K. Hall, “Features of mono- and multinucleated bone resorbing cells of the zebrafish Danio rerio and their contribution to skeletal development, remodeling, and growth,” Journal of Morphology, vol. 250, no. 3, pp. 197–207, 2001. View at Publisher · View at Google Scholar · View at PubMed
  75. P. H. Wimberger, “Plasticity of jaw and skull morphology in the Neotropical cichlids Geophagus brasiliensis and G. steindachneri,” Evolution, vol. 45, pp. 1545–1563, 1991. View at Google Scholar
  76. N. Bouton, F. Witte, and J. J. M. Van Alphen, “Experimental evidence for adaptive phenotypic plasticity in a rock-dwelling cichlid fish from Lake Victoria,” Biological Journal of the Linnean Society, vol. 77, no. 2, pp. 185–192, 2002. View at Publisher · View at Google Scholar
  77. X. Nie, K. Luukko, and P. Kettunen, “BMP signalling in craniofacial development,” International Journal of Developmental Biology, vol. 50, no. 6, pp. 511–521, 2006. View at Publisher · View at Google Scholar · View at PubMed
  78. J. Cebra-Thomas, F. Tan, S. Sistla et al., “How the turtle forms its shell: a paracrine hypothesis of carapace formation,” Journal of Experimental Zoology B, vol. 304, no. 6, pp. 558–569, 2005. View at Publisher · View at Google Scholar · View at PubMed
  79. K. E. Sears, R. R. Behringer, J. J. Rasweiler, and L. A. Niswander, “Development of bat flight: morphologic and molecular evolution of bat wing digits,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 17, pp. 6581–6586, 2006. View at Publisher · View at Google Scholar · View at PubMed
  80. Y. Terai, N. Morikawa, and N. Okada, “The evolution of the pro-domain of bone morphogenetic protein 4 (Bmp4) in an explosively speciated lineage of East African cichlid fishes [2],” Molecular Biology and Evolution, vol. 19, no. 9, pp. 1628–1632, 2002. View at Google Scholar
  81. A. Abzhanov, M. Protas, B. R. Grant, P. R. Grant, and C. J. Tabin, “Bmp4 and morphological variation of beaks in Darwin's finches,” Science, vol. 305, no. 5689, pp. 1462–1465, 2004. View at Publisher · View at Google Scholar · View at PubMed
  82. P. Wu, T. X. Jiang, S. Suksaweang, R. B. Widelitz, and C. M. Chuong, “Molecular shaping of the beak,” Science, vol. 305, no. 5689, pp. 1465–1466, 2004. View at Publisher · View at Google Scholar · View at PubMed
  83. P. Wu, T. X. Jiang, J. Y. Shen, R. B. Widelitz, and C. M. Chuong, “Morphoregulation of avian beaks: comparative mapping of growth zone activities and morphological evolution,” Developmental Dynamics, vol. 235, no. 5, pp. 1400–1412, 2006. View at Publisher · View at Google Scholar · View at PubMed
  84. M. Ikegame, O. Ishibashi, T. Yoshizawa et al., “Tensile stress induces bone morphogenetic protein 4 in preosteoblastic and fibroblastic cells, which later differentiate into osteoblasts leading to osteogenesis in the mouse calvariae in organ culture,” Journal of Bone and Mineral Research, vol. 16, no. 1, pp. 24–32, 2001. View at Google Scholar
  85. J. J. Mao and H. D. Nah, “Growth and development: hereditary and mechanical modulations,” American Journal of Orthodontics and Dentofacial Orthopedics, vol. 125, no. 6, pp. 676–689, 2004. View at Publisher · View at Google Scholar
  86. P. Aspenberg, N. Basic, M. Tagil, and S. Vukicevic, “Reduced expression of BMP-3 due to mechanical loading: a link between mechanical stimuli and tissue differentiation,” Acta Orthopaedica Scandinavica, vol. 71, no. 6, pp. 558–562, 2000. View at Google Scholar
  87. T. Nakase, S. Nomura, H. Yoshikawa et al., “Transient and localized expression of bone morphogenetic protein 4 messenger RNA during fracture healing,” Journal of Bone and Mineral Research, vol. 9, no. 5, pp. 651–659, 1994. View at Google Scholar
  88. M. Sato, T. Ochi, T. Nakase et al., “Mechanical tension-stress induces expression of bone morphogenetic protein (BMP)-2 and BMP-4, but not BMP-6, BMP-7, and GDF-5 mRNA, during distraction osteogenesis,” Journal of Bone and Mineral Research, vol. 14, no. 7, pp. 1084–1095, 1999. View at Publisher · View at Google Scholar · View at PubMed