International Journal of Evolutionary BiologyVolume 2012 (2012), Article ID 538146, 9 pagesdoi:10.1155/2012/538146

Research Article

More than Meets the Eye: Functionally Salient Changes in Internal Bone Architecture Accompany Divergence in Cichlid Feeding Mode

R. Craig Albertson,¹ W. James Cooper,² and Kenneth A. Mann³

¹Department of Biology, University of Massachusetts Amherst, 611 North Pleasant Street, Amherst, MA 01003, USA²School of Biological Sciences, Washington State University Tri-Cities, 2710 Crimson Way, Richland, WA 99354, USA³Institute for Human Performance, Department of Orthopedic Surgery, SUNY Upstate Medical University, 750 East Adams Street, Syracuse, NY 13210, USA

Received 17 January 2012; Accepted 14 March 2012

Academic Editor: Tetsumi Takahashi

Copyright © 2012 R. Craig Albertson 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.

Abstract.

African cichlids have undergone extensive and repeated adaptive radiations in foraging habitat. While the external morphology of the cichlid craniofacial skeleton has been studied extensively, biomechanically relevant changes to internal bone architecture have been largely overlooked. Here we explore two fundamental questions: (1) Do changes in the internal architecture of bone accompany shifts in foraging mode? (2) What is the genetic basis for this trait? We focus on the maxilla, which is an integral part of the feeding apparatus and an element that should be subjected to significant bending forces during biting. Analyses of CT scans revealed clear differences between the maxilla of two species that employ alternative foraging strategies (i.e., biting versus suction feeding). Hybrids between the two species exhibit maxillary geometries that closely resemble those of the suction feeding species, consistent with a dominant mode of inheritance. This was supported by the results of a genetic mapping experiment, where suction feeding alleles were dominant to biting alleles at two loci that affect bone architecture. Overall, these data suggest that the internal structure of the cichlid maxilla has a tractable genetic basis and that discrete shifts in this trait have accompanied the evolution of alternate feeding modes.

1. Introduction

Adaptive radiations involve the concomitant evolution of ecological and phenotypic diversity within a rapidly multiplying lineage [1], and many of the most notable adaptive radiations are characterized by divergence in functional morphology. Hawaiian silverswords, for example, have evolved a suite of morphological traits associated with adaptations to an extreme range of environmental moisture (mesic to xeric) [2, 3]; Anolis lizards have diversified in regard to traits involved in clinging and climbing abilities [4–8]; both Galápagos finches and African cichlids are renowned for their extensive (and in the case of cichlids, repeated) adaptive radiations in trophic morphology that parallel, and presumably contribute to, microhabitat divergence in foraging niches [9–16]. Not surprisingly, considerable attention has been given to characterizing the phenotypic diversity associated with these extraordinary radiations [3, 4, 9, 14, 17–19].

In the case of the multiple adaptive radiations of East-African cichlids, extensive analyses of their anatomical diversity have only recently been undertaken [14, 20, 21]. Among the notable findings from this body of work is that patterns of diversification within each of the three large lakes in the region (Victoria, Tanganyika, and Malawi) are statistically similar to one another [14, 21]. In particular, previous work from our group has found that trophic variation among cichlid radiations is characterized by divergence along a conserved ecomorphological axis [14]. One end of this axis is defined by species that forage in the water column and possess elongated jaws, while the opposite end is characterized by species that feed on benthic prey items using significantly shorter jaws. Thus, the primary axis of craniofacial variation defined by East-African rift-lake cichlids distinguishes benthic from pelagic ecotypes. The concordance between morphology and foraging mode observed in this study makes sense within the more general context of teleost functional morphology. Fish with short jaws have the potential, all other factors being equal, to produce bites that are proportionally more powerful due to an increased mechanical advantage employed by the jaw adductor muscles during biting, which is advantageous for herbivores that scrape tough, filamentous plant material from the substrate and for benthic predators that generate larger bite forces in order to crush, detach, or uncover their prey [22–25]. Longer jaws, on the other hand, facilitate the capture of more elusive prey by increasing bite speed and promoting greater jaw protrusion [23, 26–31].

While the functional implications of variation in external craniofacial geometry have been extensively studied in fishes [24, 28–30, 32–34] and many other vertebrates [18, 35–39], the examination of internal bone architecture has been less prominent with respect to adaptive radiations in fishes (but see [40, 41]). This paucity of data likely reflects the effort and expense associated with obtaining descriptions of these phenotypes. Specifically, while μCT scanning is becoming increasingly accessible to more research labs and is therefore being applied to the study of a steadily increasing number of taxa, the collection and processing of this type of data remains time-consuming, computationally intensive, and expensive in comparison to imaging external bone shape (which may only require light photography). Since adaptive radiations, by definition, result in species-rich and/or highly diverse lineages, the scanning of large numbers of skeletons is simply not feasible for most labs. Here we mitigate these limitations by focusing on one of the more stalwart, and functionally relevant, bones in the face (the maxilla) and by taking advantage of our current knowledge of cichlid adaptive radiations. In particular, we focus our analyses on species that define opposite ends of the functional continuum that characterizes the primary axis of craniofacial variation among Lake Malawi cichlid species [14]. In this way we can identify and describe trends that are associated with the primary axis of diversification of this lineage as well as generate a predictive framework for more global patterns of functional divergence among cichlids and other fish species. We find that discrete changes in the internal architecture of the maxilla have accompanied shifts in foraging mode within this group. These anatomical changes are biomechanically relevant and predict that biting species possess bone that is more resistant to force transmission compared to pelagic suction feeders. Finally, we show that variation in this trait has a relatively simple genetic basis, which suggests that it can respond quickly to natural selection. We submit that a more comprehensive understanding of the genetic architecture and phenotypic variation of this functionally important trait should be a priority of future research in this and other adaptive radiations defined by divergence in feeding morphology.

2. Methods

2.1. Focal Species

Two closely related Lake Malawi cichlid species that employ alternate modes of feeding (biting versus suction) were analyzed for this study. Labeotropheus fuelleborni (LF) is a member of the rock-dwelling clade of Malawi cichlids that is specialized to scrape tough, filamentous algae from the substrate [42, 43], and it has one of the most extreme craniofacial architectures of any lake-dwelling cichlid species within this region [14]. It possesses a short, stout head, steeply rounded craniofacial profile and wide jaws that are configured to employ high mechanical advantage during biting. Maylandia zebra (MZ; the genus name Metriaclima, which the authors have used elsewhere, is a junior synonym of Maylandia; [44]) is a closely related, but more generalized rock-dwelling species that collects plankton from the water column and brushes loose algae and detritus from rocky substrates. To accommodate this alternate mode of feeding, MZ has a relatively long head, shallow skull profile, and elongated jaws that are configured to produce faster but weaker bites (i.e., lower MAs) relative to LF [42, 43].

We have shown previously that the forces generated during biting will be transmitted from the lower jaw, through the maxillae, and to the anterior portions of the neurocranium and palatine ([34], Figure 1). Bending force load should also be high in the maxilla, since it acts as a lever that pivots around the pterygoid process of the palatine (Figure 1), and which is moved by the A1 division of the adductor mandibulae muscle during biting and by its connection to the lower jaw during mouth opening [32]. The shape of the maxillae is conspicuously different in LF and MZ (Figures 1 and 2(a)), with LF possessing an element that is much wider and more conspicuously bent along the medial-lateral axis compared to MZ, where the maxilla is thin and straight. See Albertson and Kocher [42] for a more comprehensive discussion of the craniofacial anatomy of these two species and Otten [45] and Cooper et al. [34] for a description of the functional anatomy of the cichlid oral jaws.

(a)

(b)

(c)

Figure 1: (a) Illustration of cichlid craniofacial anatomy in the lateral view. (b) Micro-CT scan of the oral jaws and associated elements. (c) Anatomy of a cichlid maxilla (left side, lateral view) showing the region imaged using μCT scanning. In panels (a) and (b) the maxilla (mx) is highlighted orange. Drawing by Kristen Ann Tietjen. AM₁: first division of the adductor mandibulae muscle; AM₂: second division of the adductor mandibulae; dnt: dentary; iop: interopercle; ncm: neurocranium; op: opercle; pal: pterygoid process of the palatine; pmx: premaxilla; soc: supraoccipital crest of the neurocranium; sop: subopercle.

(a)

(b)

(c)

(d)

Figure 2: (a) LF and MZ exhibit clear differences in internal bone architecture. Both frontal and transverse sections are shown. Frontal sections were taken approximately halfway through the bone. Lines through the elements show the level at which transverse sections were taken. Differences in bone architecture were quantified as bone bending stiffness (i.e., mm⁴, (b)) and bone cross-sectional area (mm², (c)). For both measures, the F₁ and F₂ hybrid generations are intermediate, with a statistical bias toward MZ values. For both (b) and (c), the “a, b, and c” indicate statistical groupings according to a two-tail t-test, and bars indicate standard errors. (d) Linear regression of bone area on bending stiffness. The relationship between bone bending stiffness and area is approximately the same for MZ and both hybrid generations, but different for LF, which are more efficient in terms of generating greater bending stiffness via the distribution of bone.

2.2. Microcomputed Tomography and Quantification of Internal Bone Architecture

Maxillae were scanned at 12-micron resolution with a microcomputed tomography (μCT) scanner (μCT 40, SCANCO Medical, Wayne, PA). Maxillae were oriented for scanning such that cross-sectional images were perpendicular to the long axis of the articular head of the bone (Figure 1(c)). These cross-sections were taken through the thinnest portion of the caudal (i.e., “neck”) region of the articular head of the maxilla. This region is roughly semicircular in cross-section (Figure 2), and is caudal to the maxilla’s premaxillary and palatinad wings and rostral to its dorsal wing (anatomy after [46]). It lies between the articulation of the palatinad wing of the maxilla with the maxillad process of the palatine (the fulcrum for maxillary rotation) and the two regions where closing and opening forces are applied to the maxilla: the insertion of the A1 division of the adductor mandibulae muscle on the medial surface of the dorsal wing (closing) and the ligamentous attachments between the shank of the maxilla and the lower jaw (opening). The maxillary shank is connected to the lower jaw by connective tissue that attaches to both the primordial process of the articular and the coronoid process of the dentary in the fishes we examined (the maxillary connection to the dentary is sometimes less extensive in other fish species). The cross-sectional areas imaged lie almost immediately between the maxillary fulcrum and the point where biting (i.e., closing) forces are directly applied to this bone (Figure 1(c)), and an ability to resist bending should therefore be a particularly important aspect of the functional morphology of this region. Image sets were exported to ImageJ (http://rsbweb.nih.gov/), and a lower threshold was defined as 400 mg/cc hydroxyapatite equivalent to isolate bone. The BoneJ plug-in to ImageJ was used to quantify bone cross-sectional area (CSA, mm²) and principal area moment of inertia (, mm⁴) of a single 2D slice within the articular neck of the maxilla (bracketed region, Figure 1(c)). CSA is a measure of the quantity of bone while is a measure of the ability of bone to resist bending loads. In all, 7 LF, 7 MZ, 7 F₁, and 49 F₂ were scanned and analyzed in this study.

2.3. Pedigree and Linkage Analysis

Details concerning the mapping population, construction of the linkage map, and quantitative trait locus (QTL) analysis have been described elsewhere [13, 34, 47, 48]. In brief, we used a pedigree derived from crossing a single LF male to a single MZ female to generate an F₂ mapping population () and a linkage map that assigned 165 markers (both microsatellites and SNPs) to 25 linkage groups using JoinMap 3.0 [49]. A linkage analysis was performed using MapQTL 4.0 [50] with as the mapping variable. Because of the time and expense required to μCT scan cichlid maxillae, we chose 49 F₂ animals with a wide range of external maxillary thicknesses for our QTL analysis.

It is important to note that, while our experimental design (i.e., bulk segregants) captured much of the variance in maxillary width among our F₂, the relatively small number of F₂ used in this experiment makes the results susceptible to the Beavis effect, in which the number of QTL tends to be underestimated and QTL effects tend to be overestimated [51]. These specific variables should therefore be interpreted with caution as they likely represent a simplified view of the genetic architecture of these traits. However, both modeling and empirical data indicate that the accuracy of QTL localization is less affected by small sample sizes [52, 53]. Nevertheless, we consider this analysis to be largely a proof of concept and the results to be preliminary.

3. Results

3.1. Distinct Internal Bone Architectures Are Associated with Divergent Feeding Modes

Micro-CT scanning revealed clear qualitative differences in internal bone morphology between LF and MZ (Figure 2(a)). We found that the maxilla in LF is hollow, with an internal bone architecture that closely resembles that of trabecular bone in mammals. The maxilla of MZ, on the other hand, is comparatively thin and solid. We reported previously that skull bone hydroxyapatite (HA) densities are roughly similar between these two species (LF:  mg/cc HA; MZ:  mg/cc HA; [34]), which suggests that any difference in biomechanical performance should be due to geometry rather than substance. To this end, we quantified differences in bone area moment of inertia (i.e., estimated bending stiffness, mm⁴, Figure 2(b)) and bone cross-sectional area (CSA, mm², Figure 2(c)). We found that the maxilla in LF contains significantly more bone than MZ and is also significantly more resistant to bending forces. Moreover, the relationship between bending stiffness and CSA suggests that the internal architecture of the maxilla in LF is more structurally stalwart per unit of bone compared to MZ. This assertion is supported by a steeper slope describing the relationship between bone stiffness and area in LF compared to that in MZ (Figure 2(d)). For both measures, the F₁ and F₂ hybrid generations were statistically biased toward MZ, suggesting a role for dominance in the inheritance of these traits. Moreover, the relationship between bending stiffness and CSA was approximately the same for the MZ, F₁ and F₂ populations.

3.2. Genetic Architecture of a Biomechanical Trait

Two significant QTL were detected for bone bending stiffness (Table 1). The first ( 1) localized to a narrow region on linkage group 7. The second QTL ( 2) localized to the distal end of linkage group 11. Both loci showed evidence for dominance of the MZ allele, which is consistent with the mean values for each population reported in Figure 2. The LF/LF genotype increased mean bending stiffness at both loci, although the mean phenotypic values of all genotypic classes were lower than what would be expected based on parental averages. This is likely due, at least in part, to our low F₂ sample size, which has also likely acted to inflate the percent variance explained (PVE) by each QTL. We cannot, however, rule out the possibly that other factors are leading to a downward bias in our F₂ values of stiffness, including environmental effects, or allometry. While we made every attempt to maintain constant rearing conditions across populations in terms of tank densities, substrate type, and diet, the F₂ were raised a couple of years after the parental and F₁ populations making it possible that there were unaccounted for differences in environment. Allometry could also be biasing our F₂ values. The average size of our F₂ population was smaller than that for either parental species or the F₁ (average standard length of 8.0 cm (F₂) versus 9.2 cm (LF), 9.3 cm (MZ), and 8.8 cm (F₁)). However, when using residuals from a regression of stiffness on size, the QTL results did not change. Clearly, this observation warrants further investigation.

Table 1: Two distinct QTL on two linkage groups (LGs) were detected for bone bending stiffness. Both loci show evidence for dominance of the MZ allele, which is consistent with the mean values for each population reported in Figure 2. The LF/LF genotype increases mean stiffness at both loci, although the mean phenotypic values of each genotypic class were lower than what would be expected based on parental averages. This is likely due to our low F₂ sample size, which has also likely inflated the percent variance explained (PVE) by each QTL.

Trait
Parental means (SE)
95% range* peak
Mean phenotype/F₂ genotype

MZ
LF
QTL
LG
cM
cM
LOD
MZ/MZ
MZ/LF
LF/LF
PVE

Stiffness [ (mm⁴)]
0.055 (0.017)
0.207 (0.024)
1
7
51–57
54
3.80
0.0326
0.0361
0.0606
23.4

2
11
49–50
50
3.10
0.0298
0.0378
0.0630
38.5

*Significance at the genomewide level.

We chose F₂ individuals for this analysis that exhibited a wide range of maxillary widths, with the intention of maximizing variance and thus the power to detect QTL. However, once these elements were scanned and stiffness was estimated, we found that width was only a weak predictor of bone stiffness (, ). In other words, these traits are segregating largely independent of one another, which suggests that they are under separate genetic control and that external skeletal anatomy cannot predict internal bone architecture. This assertion is supported by the observation that neither of the bending stiffness QTL fell within intervals that were previously implicated in maxillary shape [13]. In fact, QTL 2 localized to a region that is distinct from all other cichlid craniofacial QTL identified to date [13, 34, 48, 54, 55]. QTL 1, on the other hand, did localize to a region that overlaps with a QTL for jaw width [13] and exhibits a nearly identical LOD distribution with a QTL for the length of the retroarticular (RA) process of the lower jaw [48]. Similar to QTL 1, LF alleles at this locus act to increase the trait value for RA length and MZ alleles appear to be dominant.

4. Discussion

4.1. Divergence in Bone Strength and Weight among Vertebrates

Bone strength and stiffness are critical for optimizing the function of skeletal elements associated with feeding and locomotion, and natural selection will favor animals that perform these functions with greater efficiency [56–58]. While both bone density and shape contribute to stiffness and strength, dense bone is heavier than less dense bone. Vertebrate bone therefore tends to be designed such that strength and stiffness are maximized and weight is minimized [58, 59]. This trade-off is especially important in flighted vertebrates, where skeletons must be lightweight to minimize the metabolic cost of flight but strong enough to withstand the torsion and shearing forces associated with powered flight. As a result, birds and bats have evolved bones that are hollow but more dense compared to those of terrestrial vertebrates [60–62].

There is also a dynamic relationship between bone strength and weight among aquatic vertebrates. Specifically, across a spectrum of vertebrate classes the modulation of bone density appears to be a mechanism for buoyancy control [63–65]. This trend is beautifully illustrated by the evolutionary history of whales, which is marked by discrete shifts in habitat from terrestrial, to semiaquatic, and finally to fully aquatic life histories. The mechanical constraints associated with locomotion in each of these habitats are very different, and as a result these evolutionary transitions were accompanied by dramatic changes in bone architecture. For example, the shift from terrestrial to semiaquatic habitats in ancient whales (i.e., archaeocetes) was accompanied by a dramatic increase in bone density. Like other large semiaquatic mammals, this adaptation was for increased mass, which is associated with benthic foraging [63]. Modern whales, on the other hand, are fully aquatic and possess a number of adaptations for life in the open water, including a largely osteoporotic skeleton [63]. While functional parameters including bone stiffness have not been examined in modern cetaceans, it is notable that osteoporotic bone in cetaceans is not observed in elements associated with feeding or locomotion (i.e., skull and vertebrae), where functional demands remain high [64]. Thus, a balance has been struck between increased buoyancy and efficient foraging and locomotion in the skeletons of modern whales.

The evolutionary history of Antarctic notothenioid fishes represents another striking example of how bone development has been modified to affect buoyancy. Antarctic notothenioids represent one of the best described adaptive radiations among marine fishes [66], and the hallmark of their evolution is the development of secondary pelagicism via alteration of buoyancy [67]. This lineage is thought to have evolved from a robustly mineralized bottom-dwelling perciform species beginning 40–60 mya when the waters of the Antarctic continental shelf were still temperate [67]. The grounding of the ice sheet on the continental shelf and changing trophic conditions led to the local extinction of the diverse late Eocene fish fauna, thus freeing pelagic niches into which the notothenioids radiated [68]. About 50% of notothenioid species now either live or forage in the pelagic habitat [69]. In many instances, secondary pelagicism has been achieved through pedomorphism, including the complete or partial retention of the notochord, delayed ossification of the skeleton, and replacement of bone by connective tissue [65, 70–72]. Similar to cetaceans, osteoporotic bone in pelagic notothenioids is most pronounced in areas of the skeleton that are not intimately associated with foraging (e.g., oral jaws) or locomotion (e.g., pectoral fins) [70].

While the examples above represent changes in bone structure at the macroevolutionary level, it is reasonable to assume that similar trends underlie microevolutionary divergence. As mentioned above, cichlids have diverged along a benthic-pelagic ecomorphological axis, and extensive modifications to the skeletal system have accompanied this divergence [14]. LF and MZ are closely related species that lie on opposing ends of this continuum, and while bone density does not appear to be different between these two species [34], LF has a more extensively mineralized skeleton (i.e., more bone in more places) [42], which is commensurate with other adaptations for a benthic mode of feeding. These findings suggest that levels and patterns of bone deposition are more evolvable in this group than are the material properties of bone (although a more rigorous survey of HA density in a greater number of elements and across more taxa is needed). We also show here that internal bone architecture appears to be surprisingly malleable among cichlids, as strikingly different cross-sectional bone shapes exist between species that employ alternate modes of feeding. This sets up clear predictions that can be tested in a larger number of species. For example, if species were arrayed along a benthic-pelagic ecomorphological axis, one might expect that this would establish a continuum of internal bone architectures. Alternatively, since LF represents a highly derived species, it is also possible that the internal bone architecture described here (i.e., high stiffness) is unique to this species. Clearly, this would be a fruitful area of future research.

4.2. Roles for the Environment versus Genetics in Determining Internal Bone Architecture

Bone geometry influences stiffness such that bone with a solid cross-section is less rigid whereas hollow bone with the same cross-sectional area is more rigid. Natural selection should therefore favor one configuration over the other depending on the task to be performed (e.g., biting versus sucking). Alternatively, given the varying functional demands imposed on the vertebrate skeletal system over ontogeny, or from season to season, selection might favor a plastic skeletal system that can respond to different mechanical stimuli. Distinguishing between these alternatives represents an important, but muddled area of research. In other words, the degree to which internal bone architecture is genetically preprogrammed or mechanically regulated remains unclear.

On one hand, both computational modeling and empirical studies offer strong support for the assertion that internal bone architecture responds to mechanical stimuli [59, 73, 74]. Alternatively, disparate vertebrate taxa have modified internal bone geometry due to novel functional demands (e.g., powered flight in birds and bats) [35, 61, 62], which suggests a genetic component for this trait. Unfortunately, compared to the relatively large body of literature dedicated to the study of the genetics of bone material properties (focused mainly on mouse mutants, reviewed by [75]), less is known about the genetic basis of bone geometry. Moreover, mutations that lead to aberrant bone architectures usually also affect material properties. “Wolff’s Law” suggests that bone adapts to mechanical stimuli to maintain a narrow range of strain (reviewed by [76]). It is therefore thought that for many/most mouse mutants where both bone material and geometry are affected, deficient material properties are the primary defect and altered geometry represents a secondary response to compensate for abnormal bone strains [76].

Cichlids offer a genetic system where internal bone architecture varies independently from material properties, thus mitigating the confounding issues associated with Wolff’s Law. Whereas HA density appears relatively conserved between the species examined here, internal geometry differs dramatically. This could be due to a fundamental constraint associated with changing material properties in fishes (e.g., higher-density bone is more brittle) or because altering bone architecture is a more efficient way to affect stiffness. For example, bending stiffness of a round bone is equal to , where the elastic modulus () is proportional to HA density of bone and area moment of inertia () is proportional to radius⁴. Doubling HA would lead to a doubling of stiffness, whereas doubling the radius would lead to a 16-fold increase in stiffness. Changing bone architecture is therefore a more efficient way to change bone function due to increased demands on bending load. Either way, the decoupling of these two properties of bone, as well as the ability to perform genetic mapping studies, offers an excellent opportunity to examine the genetic basis of internal bone geometry. Moreover, the ability to rear cichlids on a range of diets (e.g., hard versus soft), thereby altering the mechanical environment in which the jaws develop, would enable an assessment of the degree to which this trait responds to the environment. Thus, cichlids represent an ideal system in which to characterize the genetic and environmental factors that influence this functionally salient trait.

4.3. Conclusions

We demonstrate that cichlids with divergent feeding morphologies and behaviors exhibit different internal bone architectures that translate to different estimates of load-bearing function. We show further that this functional trait has a tractable genetic basis. Since bone geometry has a profound effect on skeletal performance and since performance determines resource use in nature, examining the genetic basis for this trait has the potential to yield important new insights into the mechanisms that have contributed to several notable adaptive radiations (e.g., those that involve divergence in feeding systems or locomotion). Moreover, continued work in the cichlid system may also contribute to an understanding of many unresolved issues in the biomedical literature, especially those focused on decoupling the genetic from epigenetic influences on internal bone geometry. In conclusion, external bone shape and size, while important in determining skeletal function, do not tell the entire story [76]. Future research should therefore be aimed at elucidating a better understanding of (1) the material properties and (2) internal geometry of skeletal elements associated with feeding and locomotion in this and other adaptive radiations.

Acknowledgments

The authors thank James Wernle for technical assistance with μCT scanning. R. C. Albertson’s cichlid research program is supported in part by a grant from the NSF (IOS-1054909).

References

D. Schluter, The Ecology of Adaptive Radiations, Oxford University Press, Oxford, UK, 2000.

B. G. Baldwin and M. J. Sanderson, “Age and rate of diversification of the Hawaiian silversword alliance (Compositae),” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 16, pp. 9402–9406, 1998.

G. D. Carr, Monograph of the Hawaiian Madiinae (Asteraceae): Argyroxiphium, Dubautia, and Wilkesia. Allertonia, vol. 4, Pacific Tropical Botanical Garden, 1985.

R. Calsbeek and D. J. Irschick, “The quick and the dead: Correlational selection on morphology, performance, and habitat use in island lizards,” Evolution, vol. 61, no. 11, pp. 2493–2503, 2007.

K. Huyghe, A. Herrel, B. Vanhooydonck, J. J. Meyers, and D. J. Irschick, “Microhabitat use, diet, and performance data on the Hispaniolan twig anole, Anolis sheplani: pushing the boundaries of morphospace,” Zoology, vol. 110, no. 1, pp. 2–8, 2007.

R. B. Langerhans, J. H. Knouft, and J. B. Losos, “Shared and unique features of diversification in greater Antillean Anolis ecomorphs,” Evolution, vol. 60, no. 2, pp. 362–369, 2006.

J. B. Losos, T. W. Schoener, K. I. Warheit, and D. Creer, “Experimental studies of adaptive differentiation in Bahamian Anolis lizards,” Genetica, vol. 112-113, pp. 399–415, 2001.

J. B. Losos, “Integrative approaches to evolutionary ecology: anolis lizards as model systems,” Annual Review of Ecology and Systematics, vol. 25, pp. 467–493, 1994.

P. R. Grant and B. R. Grant, “Evolution of character displacement in Darwin's finches,” Science, vol. 313, no. 5784, pp. 224–226, 2006.

A. Sato, C. O'huigin, F. Figueroa et al., “Phylogeny of Darwin's finches as revealed by mtDNA sequences,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 9, pp. 5101–5106, 1999.

A. Sato, H. Tichy, C. O'Huigin, P. R. Grant B, R. Grant, and J. Klein, “On the origin of Darwin's finches,” Molecular Biology and Evolution, vol. 18, no. 3, pp. 299–311, 2001.

R. C. Albertson and T. D. Kocher, “Genetic and developmental basis of cichlid trophic diversity,” Heredity, vol. 97, no. 3, pp. 211–221, 2006.

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.

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.

I. Kornfield and P. F. Smith, “African cichlid fishes: model systems for evolutionary biology,” Annual Review of Ecology and Systematics, vol. 31, pp. 163–196, 2000.

T. D. Kocher, “Adaptive evolution and explosive speciation: the cichlid fish model,” Nature Reviews Genetics, vol. 5, no. 4, pp. 288–298, 2004.

K. F. Liem, “Functional morphology,” in Cichlid Fishes: Behavior, Ecology and Evolution, M. H. A. Keenleyside, Ed., pp. 129–150, Chapman and Hall, London, UK, 1991.

A. Herrel, J. Podos, B. Vanhooydonck, and A. P. Hendry, “Force-velocity trade-off in Darwin's finch jaw function: a biomechanical basis for ecological speciation?” Functional Ecology, vol. 23, no. 1, pp. 119–125, 2009.

J. Podos and S. Nowicki, “Beaks, adaptation, and vocal evolution in Darwin's finches,” BioScience, vol. 54, no. 6, pp. 501–510, 2004.

P. Chakrabarty, “Testing conjectures about morphological diversity in cichlids of Lakes Malawi and Tanganyika,” Copeia, no. 2, pp. 359–373, 2005.

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.

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 24, 2009.

P. C. Wainwright and B. A. Richard, “Predicting patterns of prey use from morphology of fishes,” Environmental Biology of Fishes, vol. 44, no. 1–3, pp. 97–113, 1995.

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.

M. W. Westneat, M. E. Alfaro, P. C. Wainwright et al., “Local phylogenetic divergence and global evolutionary convergence of skull function in reef fishes of the family Labridae,” Proceedings of the Royal Society B, vol. 272, no. 1567, pp. 993–1000, 2005.

D. C. Collar, B. C. O'Meara, P. C. Wainwright, and T. J. Near, “Piscivory limits diversification of feeding morphology in centrarchid fishes,” Evolution, vol. 63, no. 6, pp. 1557–1573, 2009.

A. N. Rice, W. J. Cooper, and M. W. Westneat, “Diversification of coordination patterns during feeding behaviour in cheiline wrasses,” Biological Journal of the Linnean Society, vol. 93, no. 2, pp. 289–308, 2008.

M. W. Westneat, “Transmission of force and velocity in the feeding mechanisms of labrid fishes (Teleostei, Perciformes),” Zoomorphology, vol. 114, no. 2, pp. 103–118, 1994.

C. D. Hulsey and F. J. García De León, “Cichlid jaw mechanics: linking morphology to feeding specialization,” Functional Ecology, vol. 19, no. 3, pp. 487–494, 2005.

T. B. Waltzek and P. C. Wainwright, “Functional morphology of extreme jaw protrusion in neotropical cichlids,” Journal of Morphology, vol. 257, no. 1, pp. 96–106, 2003.

L. A. Ferry-Graham, P. C. Wainwright, and D. R. Bellwood, “Prey capture in long-jawed butterflyfishes (Chaetodontidae): the functional basis of novel feeding habits,” Journal of Experimental Marine Biology and Ecology, vol. 256, no. 2, pp. 167–184, 2001.

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.

M. W. Westneat, “Skull biomechanics and suction feeding in fishes,” in Fish Biomechanics, R. E. Shadwick and G. V. Lauder, Eds., Elsevier Academic Press, San Diego, Calif, USA, 2006.

W. J. Cooper, J. Wernle, K. Mann, and R. C. Albertson, “Functional and Genetic Integration in the Skulls of Lake Malawi Cichlids,” Evolutionary Biology, vol. 38, no. 3, pp. 316–334, 2011.

E. R. Dumont, “Feeding mechanisms in bats: variation within the constraints of flight,” Integrative and Comparative Biology, vol. 47, no. 1, pp. 137–146, 2007.

D. S. Strait, G. W. Weber, S. Neubauer et al., “The feeding biomechanics and dietary ecology of Australopithecus africanus,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 7, pp. 2124–2129, 2009.

A. Herrel, S. E. Vincent, M. E. Alfaro, S. Van Wassenbergh, B. Vanhooydonck, and D. J. Irschick, “Morphological convergence as a consequence of extreme functional demands: examples from the feeding system of natricine snakes,” Journal of Evolutionary Biology, vol. 21, no. 5, pp. 1438–1448, 2008.

A. Herrel, B. Vanhooydonck, K. Huyghe, R. Van Damme, and D. J. Irschick, “Rapid evolutionary divergence in feeding mechanics after the colonization of new environments in lizards,” Comparative Biochemistry and Physiology, vol. 148, p. S52, 2007.

A. Herrel, P. Aerts, J. Fret, and F. de Vree, “Morphology of the feeding system in agamid lizards: ecological correlates,” Anatomical Record, vol. 254, no. 4, pp. 496–507, 1999.

C. D. Hulsey, R. J. Roberts, A. S. P. Lin, R. Guldberg, and J. T. Streelman, “Convergence in a mechanically complex phenotype: detecting structural adaptations for crushing in cichlid fish,” Evolution, vol. 62, no. 7, pp. 1587–1599, 2008.

A. P. Summers, R. A. Ketcham, and T. Rowe, “Structure and function of the horn shark (Heterodontus francisci) cranium through ontogeny: development of a hard prey specialist,” Journal of Morphology, vol. 260, no. 1, pp. 1–12, 2004.

R. C. Albertson and T. D. Kocher, “Assessing morphological differences in an adaptive trait: a landmark-based morphometric approach,” Journal of Experimental Zoology, vol. 289, no. 6, pp. 385–403, 2001.

A. J. Ribbink, B. A. Marsh, A. C. Marsh, A. C. Ribbink, and B. J. Sharp, “A preliminary survey of the cichlid fishes of rocky habitats in Lake Malawi,” South African Journal of Zoology, vol. 18, no. 3, pp. 149–310, 1983.

Catalog of Fishes electronic version, http://research.calacademy.org/ichthyology/catalog/fishcatsearch.html

E. Otten, “The jaw mechanism during growth of a generalized Haplochromis species: H. Elegans trewavas 1933 (Pisces, Cichlidae),” Netherlands Journal of Zoology, vol. 33, pp. 55–98, 1983.

C. D. N. Barel, F. Witte, and M. J. W. Van Ooijen, “The shape of the skeletal elements in the head of a generalized Haplochromis species: H. elegans Trewavas 1933 (Pisces, Cichlidae),” Netherlands Journal of Zoology, vol. 26, no. 2, pp. 163–265, 1976.

R. C. Albertson, J. T. Streelman, and T. D. Kocher, “Genetic basis of adaptive shape differences in the cichlid head,” Journal of Heredity, vol. 94, no. 4, pp. 291–301, 2003.

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.

J. W. Van Ooijen and R. E. Voorrips, “JoinMap 3.0: software for the calculation of genetic linkage maps,” in Plant Research International, Wageningen UR, Wageningen, The Netherlands, 2001.

J. W. Van Ooijen, M. P. Boer, R. C. Jansen, and C. Maliepaard, “QTLMap 4.0: Software for the calculation of QTL positions on genetic maps,” in Plant Research International, Wageningen UR, Wageningen, The Netherlands, 2002.

W. D. Beavis, “The power and deceit of QTL experiments: lessons from comparative QTL studies,” in Proceedings of the 49th Annual Corn and Sorghum Research Conference, pp. 252–268, American Seed Trade Association, Washington, DC, USA, 1994.

W. D. Beavis, “QTL analyses: power, precision, and accuracy,” in Molecular Dissection of Complex Traits, A. H. Paterson, Ed., pp. 145–161, CRC Press, Boca Raton, Fla, USA, 1998.

A. H. Price, “Believe it or not, QTLs are accurate!,” Trends in Plant Science, vol. 11, no. 5, pp. 213–216, 2006.

T. A. Stewart and R. C. Albertson, “Evolution of a unique predatory feeding apparatus: functional anatomy, development and a genetic locus for jaw laterality in Lake Tanganyika scale-eating cichlids,” BMC Biology, vol. 8, article 8, 2010.

K. Parsons, E. Marquez, and R. C. Albertson, “Constraint and opportunity: the genetic basis and evolution of modularity in the cichlid mandible,” American Naturalist, vol. 179, no. 1, pp. 64–78, 2012.

R. A. Anderson, L. D. Mcbrayer, and A. Herrel, “Bite force in vertebrates: opportunities and caveats for use of a nonpareil whole-animal performance measure,” Biological Journal of the Linnean Society, vol. 93, no. 4, pp. 709–720, 2008.

A. A. Biewener, “Biomechanical consequences of scaling,” Journal of Experimental Biology, vol. 208, no. 9, pp. 1665–1676, 2005.

J. D. Currey, Bones: Structure and Mechanics, Princeton University Press, Princeton, NJ, USA, 2002.

R. Huiskes, R. Rulmerman, G. H. Van Lenthe, and J. D. Janssen, “Effects of mechanical forces on maintenance and adaptation of form in trabecular bone,” Nature, vol. 405, no. 6787, pp. 704–706, 2000.

J. D. Currey and R. M. Alexander, “The thickness of the walls of tubular bones,” Journal of Zoology, vol. 206, pp. 453–468, 1985.

E. R. Dumont, “Bone density and the lightweight skeletons of birds,” Proceedings of the Royal Society B, vol. 277, no. 1691, pp. 2193–2198, 2010.

S. M. Swartz, M. B. Bennett, and D. R. Carrier, “Wing bone stresses in free flying bats and the evolution of skeletal design for flight,” Nature, vol. 359, no. 6397, pp. 726–729, 1992.

N. M. Gray, K. Kainec, S. Madar, L. Tomko, and S. Wolfe, “Sink or swim? Bone density as a mechanism for buoyancy control in early cetaceans,” Anatomical Record, vol. 290, no. 6, pp. 638–653, 2007.

A. Ricqles and V. Buffrenil, “Bone histology, heterochronies and the return of Tetrapods to life in water: where are we?” in Secondary Adaptation of Tetrapods to Life in Water, J.-M. Mazin and V. Buffrenil, Eds., pp. 289–306, Verlag Dr. Friedrich Pfeil, Munchen, Germany, 2001.

O. S. Voskoboynikova, “Evolution of the visceral skeleton and phylogeny of Nototheniidae,” Journal of Ichthyology, vol. 33, pp. 23–47, 1993.

J. T. Eastman and A. R. McCune, “Fishes on the Antarctic continental shelf: Evolution of a marine species flock?” Journal of Fish Biology, vol. 57, pp. 84–102, 2000.

J. T. Eastman, “The nature of the diversity of Antarctic fishes,” Polar Biology, vol. 28, no. 2, pp. 93–107, 2005.

J. T. Eastman, Antarctic Fish Biology: Evolution in a Unique Environment, Academic Press, San Diego, Calif, USA, 1993.

J. T. Eastman, “Phyletic devergence and specialization for pelagiclife in the Antarctic nototheniid fish Pleuragrammaan tarcticum,” Comparative Biochemistry and Physiology, vol. 118, no. 4, pp. 1095–1101, 1997.

R. C. Albertson, Y. L. Yan, T. A. Titus et al., “Molecular pedomorphism underlies craniofacial skeletal evolution in Antarctic notothenioid fishes,” BMC Evolutionary Biology, vol. 10, no. 1, article 4, 2010.

O. S. Voskoboynikova, “Rates of individual development of the bony skeleton of eleven species of the family Nototheniidae,” Journal of Ichthyology, vol. 34, pp. 108–120, 1994.

O. S. Voskoboynikova, O. Y. Tereshchuk, and A. Kellermann, “Osteological development of the Antarctic silverfish Pleuragramma antarcticum (Nototheniidae),” Cybium, vol. 18, pp. 251–271, 1994.

J. N. Dowthwaite, R. M. Hickman, J. A. Kanaley, R. J. Ploutz-Snyder, J. A. Spadaro, and T. A. Scerpella, “Distal radius strength: a comparison of DXA-derived vs pQCT-measured parameters in adolescent females,” Journal of Clinical Densitometry, vol. 12, no. 1, pp. 42–53, 2009.

J. N. Dowthwaite and T. A. Scerpella, “Skeletal geometry and indices of bone strength in artistic gymnasts,” Journal of Musculoskeletal Neuronal Interactions, vol. 9, no. 4, pp. 198–214, 2009.

P. Chavassieux, E. Seeman, and P. D. Delmas, “Insights into material and structural basis of bone fragility from diseases associated with fractures: how determinants of the biomechanical properties of bone are compromised by disease,” Endocrine Reviews, vol. 28, no. 2, pp. 151–164, 2007.

M. C. H. Van der Meulen, K. J. Jepsen, and B. Miki, “Understanding bone strength: size isn't everything,” Bone, vol. 29, no. 2, pp. 101–104, 2001.