Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Is Beak Morphology in Darwin’s Finches Tuned to Loading Demands?

  • Joris Soons ,

    joris.soons@ua.ac.be

    Affiliation Laboratory of Biomedical Physics, University of Antwerp, Groenenborgerlaan 171, 2020, Antwerpen, Belgium

  • Annelies Genbrugge,

    Affiliations Laboratory of Biomedical Physics, University of Antwerp, Groenenborgerlaan 171, 2020, Antwerpen, Belgium, Evolutionary Morphology of Vertebrates, Ghent University, K.L. Ledeganckstraat 35, 9000, Gent, Belgium

  • Jeffrey Podos,

    Affiliation Department of Biology, 221 Morrill Science Center, University of Massachusetts, Amherst, Massachusetts, 01003, United States of America

  • Dominique Adriaens,

    Affiliation Evolutionary Morphology of Vertebrates, Ghent University, K.L. Ledeganckstraat 35, 9000, Gent, Belgium

  • Peter Aerts,

    Affiliations Department of Biology, University of Antwerp, Universiteitsplein 1, 2610, Antwerpen, Belgium, Department of movement and sports sciences, University of Ghent, Watersportlaan 2, 9000, Gent, Belgium

  • Joris Dirckx,

    Affiliation Laboratory of Biomedical Physics, University of Antwerp, Groenenborgerlaan 171, 2020, Antwerpen, Belgium

  • Anthony Herrel

    Affiliations Department of Biology, University of Antwerp, Universiteitsplein 1, 2610, Antwerpen, Belgium, Département d’Ecologie et de Gestion de la Biodiversité, Museum National d’Histoire Naturelle, 57 rue Cuvier, Case postale 55, 75231, Paris, Cedex 5, France

Abstract

One of nature's premier illustrations of adaptive evolution concerns the tight correspondence in birds between beak morphology and feeding behavior. In seed-crushing birds, beaks have been suggested to evolve at least in part to avoid fracture. Yet, we know little about mechanical relationships between beak shape, stress dissipation, and fracture avoidance. This study tests these relationships for Darwin's finches, a clade of birds renowned for their diversity in beak form and function. We obtained anatomical data from micro-CT scans and dissections, which in turn informed the construction of finite element models of the bony beak and rhamphotheca. Our models offer two new insights. First, engineering safety factors are found to range between 1 and 2.5 under natural loading conditions, with the lowest safety factors being observed in species with the highest bite forces. Second, size-scaled finite element (FE) models reveal a correspondence between inferred beak loading profiles and observed feeding strategies (e.g. edge-crushing versus tip-biting), with safety factors decreasing for base-crushers biting at the beak tip. Additionally, we identify significant correlations between safety factors, keratin thickness at bite locations, and beak aspect ratio (depth versus length). These lines of evidence together suggest that beak shape indeed evolves to resist feeding forces.

Introduction

The often tight correspondence between bird beaks and plant morphology well-illustrates the power and precision of natural selection [1]. Within bird populations, subtle variations in beak morphology can affect foraging efficiency [2], and in some cases tip the balance between survival and starvation [3]. For seed-crushing birds, however, analyses of feeding capacity require a broader view than just beak morphology per se. This is because a bird’s ability to crush seeds is determined mainly by bite force capacity, which in turn depends primarily on the orientation and cross sectional area of the jaw closer muscles [46]. These muscles, situated at the back of the head, generate crushing forces that are transferred to food by means of the upper and lower beak through a complex kinetic beak apparatus [611]. Beak morphology, by contrast, likely evolves to facilitate successful food manipulation and song production. However, beaks also need to avoid structural failure during loading [4]. The loading regime itself presumably varies with food manipulation strategies (e.g. tip biting versus base crushing) and may thus drive the evolution of beak shape [12].

We here test, in Darwin’s finches of the Galápagos Islands, the hypothesis that beak morphology provides structural integrity, i.e., that it evolves to resist feeding forces and avoid fracture [13]. Beaks in Darwin’s finches are known to vary broadly across and within species, and to evolve via natural selection as a response to variation in food type, food availability, and interspecific competition [3,10,1418]. Our main approach is to develop and apply finite element (FE) models of the upper beak. FE models enable exploration of the effect of complex shape variation on stress magnitude and distribution [1921], and allows us to draw inferences about patterns of loading during biting. We include here 13 species and incorporate data into our models for the bones, the jaw closer muscles, and the keratinous rhamphotheca that encapsulates the beak. Including the rhamphotheca in the model allows us to provide a more realistic estimate of stresses incurred during biting [22,23] that we achieved in previous studies [24,25].

The goal of this paper is twofold. First we examine how natural loadings are reflected in safety factors for a subset of species for which physiological FE models could be created (i.e. for which data on muscles could be obtained). Second, we created models of all 13 species scaled to a common surface area to muscle ratio [26]. This allows us to test whether variation in beak shape matches differences in food handling behavior. We predict that species will have beak shapes that allow them to minimize peak stresses during those feeding behaviors most commonly employed. In particular, we predict that base-crushing species should show the lowest stresses when biting at the base of the beak, that tip-biting species should show the lowest stresses during tip-biting, and that probing beaks should show higher stresses during both tip- and base crushing compared to those species specialized for these behaviors.

Materials and Methods

Specimen collection and CT scanning

No animals were killed for the purpose of this study. Road-killed specimens were collected during February-March of 2005 and 2006 on the main road connecting the airport on Baltra with Puerto Ayora, Santa Cruz Island. Specimens were collected under a salvage permit provided by the Galápagos National Park Service. Intact specimens were collected and preserved in a 10% aqueous formaldehyde solution for 24 h, rinsed and transferred to a 70% aqueous ethanol solution. One individual of Geospizia fortis, Geospizia fuliginosa, Geospizia scandens, Platyspiza crassirostris, Certhidea olivacea, and Camarhynchus parvulus were scanned at the UGCT scanning facility (www.ugct.ugent.be), using a micro-focus directional type X-ray tube, set at a voltage of 80 kV and a spot size of 10 mm. Specimens were mounted on a controllable rotating table (MICOS, UPR160F-AIR). For each specimen, a series of 1000 projections of 940 by 748 pixels were recorded covering 360 degrees.

Museum specimens of Geospizia magnirostris, Geospizia difficilis, Pinaroloxias inornata, Cactospiza pallida, Camarhynchus psittacula, Geospiza conirostris, and Camarhynchus pauper from the collections of the California Academy of Sciences maintained in a 70% aqueous ethanol solution were scanned at the Harvard CNS facility using an X-Tek XRA-002 micro-CT imaging system set at 75 kV. Specimens were mounted on a rotating table and a series of 3142 projections of 2000 by 2000 pixels covering 360 degrees was recorded. The voxel size of the scans was dependent on the specimen size and ranged from 2.96 μm for the smallest species to 47.00 μm for the larger species. As such, the bone shell of the bony core was represented by multiple voxels, and details within the bony beak were visible.

Muscle data

Muscle data were available for five of the species collected in the field as road-killed specimens: G. fortis, G. fuliginosa, G. scandens, C. olivacea and C. parvulus; data for three additional species were derived through dissection of museum specimens from the California Academy of Sciences (G. magnirostris, G. conirostris, and G. difficilis). These species were dissected and all muscle bundles of the jaw removed individually [27]. Muscles were blotted dry and weighed on a Mettler microbalance (±0.01mg). Next, muscles were transferred individually to Petri dishes and submerged in a 30% aqueous nitric acid solution for 18 h to dissolve all connective tissue [28]. After removal of nitric acid, muscles were transferred to a 50% aqueous glycerol solution and fibers were teased apart using blunt-tipped glass needles. Thirty fibers were selected from each muscle bundle (Table 1) and drawn using a binocular scope with attached camera lucida. A background grid was also drawn in each image to provide an object for scaling. Drawings were scanned and fiber lengths determined using imageJ (freely available at http://rsb.info.nih.gov/ij/).

thumbnail
Table 1. Table summarizing muscle mass and fiber length data used to calculate the external forces acting on the upper beak of the different species of Darwin’s finch.

https://doi.org/10.1371/journal.pone.0129479.t001

Based on muscle mass and fiber length, the physiological cross-sectional area of each muscle bundle was estimated assuming a muscle density of 1036kg/m3 [29]. Since pennate muscles were separated into their individual bundles, no additional correction for pennation angle was included. Corresponding force-generation capacities for each muscle were calculated assuming a muscle stress of 30N/cm2 [30]. As the external adductor and pseudotemporalis muscles act only indirectly on the upper mandible [6,9,10,31], the component of the muscle force transferred to the upper mandible was calculated taking into account the position of the muscles and their angles relative to the jugal bone (which transfers the forces from the lower jaw and the quadrate to the upper beak by pulling the beak downward when quadrate is rotated backward during jaw closing). The pterygoid muscle bundles act directly on the upper mandible [6,9,10,31], and muscle forces are directly transmitted through the pterygoid/palatine complex. The combined muscle forces applied to the model are shown in Table 2.

thumbnail
Table 2. Von Mises stress for different loading conditions (LC) for 13 Darwin’s finches (PB: physiological base biting; PT: physiological tip biting; FB: base biting scaled to fortis; FT: tip biting scaled to fortis).

https://doi.org/10.1371/journal.pone.0129479.t002

Segmentation and FE modelling

The segmentation of the bony core and the keratin layer (Amira 4.1 64-bit version, TGS systems) was done for all species. The tetrahedral grid construction [32], the assigned material properties (homogenuous and isotropic linear elastic with Ekeratin = 1.7 GPa and Ebone = 7.3 GPa) [22] and the boundary conditions are similar to those reported previously [2224, 33]. Convergence testing on the results, namely the von Mises (vM) stresses and the bite forces, was conducted so they change less than 5% when doubling mesh size. The final FE grids of approximately 2–3 million 4-noded linear tetrahedral elements were obtained. Consequently, multiple layers of elements were available to model the thin bony walls.

For every specimen, a FE model for base and tip loading was established. In this study we were only interested in the stresses on the upper beak, and as such elements were fully constrained (in all 6 degrees of freedom) at the back of the beak at the neurocranium. Hinging of the upper beak is possible due to bending of the thin bone (BA in Fig 1). In order to avoid singularities around point constraints, we constrained elements (in x, y and z direction, no rotation) at the ventral surface of the rhamphotheca simulating the bite constraint (T and Bs in Fig 1); the calculated reaction force of these elements is the bite force which is of equal magnitude yet of opposite sign as the seed reaction force. Bite positions were estimated using in vivo observations of animals cracking or manipulating seeds in the field [13]. Elements were constrained for a central tip bite close to the tip of the bony core, as birds are never seen to bite at the very tip in the wild (T in Fig 1). For an unilateral base bite we used a point at one fourth of the total length of the rhamphotheca, close to the edge of the bony core (Bs in Fig 1), hence creating homologous bite points across species. The available muscle forces (Table 2) estimated through dissection were applied to rigid elements at the end of the jugal and pterygoid bones in the direction of these bones (J and P in Fig 1). This resulted in physiological FE models, for both base and tip loading, for eight species (16 models). In addition we developed scaled FE models for all 13 species. We therefore measured the beak surface of all beaks and scaled them in such a way that they have the same surface area as G. fortis. Next we applied the muscle force of G. fortis to these size-scaled beaks. G. fortis was chosen because we had muscle data and scans available for multiple individuals of this species [27], but we could have picked one of the other 7 finches. As such, the effects of size are eliminated and the effect of beak shape variation on the mechanical behavior of the beak can be isolated for the 13 species of Darwin’s finches included here [26].

thumbnail
Fig 1. Schematic representation of our multi-layered (bone: B, keratin: K) finite element modeling approach, for the medium ground finch Geospiza fortis.

Bending area (Ba) and bite position (base: Bs, or tip: T) were constrained in our models for translation and rotation, and muscle forces were applied in our models via the jugal (J) and palatine (P) jaw bones (black elements are constrained). Locations of vM Stress recordings are indicated with transparant ellipses (1b: on top of bone, near base bite position; 1t: on top of bone near tip bite position; 2: on top of the beak near the nasal hinge; 3: nasal bone).

https://doi.org/10.1371/journal.pone.0129479.g001

Bite force and maximum vM stress were obtained for all 40 FE models (Table 2): 16 physiological models (tip and base loading for 8 beaks with available muscle data) and 24 models (tip and base loading for all 13 beaks size-scaled to G. fortis, minus G. fortis, see Table 2). Maximum vM stress was taken according to Saint-Venant's Principle over three different locations displaying the highest stresses on the bony core, away from areas influenced by model constraints (Figs 1 and 2, Table 2). We used 98th percentile to select maximum vM stresses presented in Fig 3 and Table 2. In these locations, the gray values of the corresponding CT images were also recorded. It has been shown [3436] that the strength of bone increases with density. Scan parameters were optimized for each CT scan and no calibration was performed, making it impossible to obtain the density directly from these gray values. Consequently, gray values were transferred to a qualitative linear scale, where 0 represents the least dense bone in the sample and 1 the densest. vM stresses on the keratin layer were not taken into account and are not displayed in Figs 2 and 4. Nevertheless, a correct modeling of the rhamphotheca is important since it has a significant influence on the stress distribution and magnitude in the bony core [23].

thumbnail
Fig 2. Physiological Finite Element Model results, side view, for 4 selected Darwin’s finch species known to use their beaks in different ways during feeding [10, 49]: G. fortis (base crushing beak), G. scandens (probing and crushing beak), C. olivacea (probing beak), and C. parvulus (tip biting beak).

Results are shown for both base (1st column) and tip (2nd column) biting simulations. Arrows indicate the location and magnitude of the calculated seed reaction forces. Warmer colors represent higher von Mises stresses. Additional FE models for the eight finches for which muscle data were available are presented in S2 and S3 Figs.

https://doi.org/10.1371/journal.pone.0129479.g002

thumbnail
Fig 3. Peak vM stress for physiological FE models (tip and base loading, eight species with available muscle data) and for scaled FE models (tip and base loading, scaled to same size as G. fortis and with the same muscle forces) for thirteen species of Darwin’s finch (* = juvenile).

https://doi.org/10.1371/journal.pone.0129479.g003

thumbnail
Fig 4. Top view of scaled FE models of the upper beaks of 13 species of Darwin’s finches.

All beaks were scaled to same size and muscle force as G. fortis, with stresses calculated for both base (left) and tip (right) biting. Warmer colors represent higher vM stresses. Note how stresses are lower during the behaviors typically employed by each species, with base crushers showing the lowest stress values during base-loading and tip crushers during tip-loading. Species with probing beaks show generally high stresses under both loading conditions. Top and side views of scaled FE models are presented in S4 and S5 Figs.

https://doi.org/10.1371/journal.pone.0129479.g004

During base loading vM stresses are not recorded for the contralateral side of the nasal bone, because a bird could decrease this stress by reducing the bite force at that side. As a consequence, the resulting reaction force on the contralateral side will be lower. It is, however, very complex to introduce this bite force tuning into our models and as such is ignored here. To quantify a birds’ risk of beak fracture during biting, we calculated safety factors for each species’ beak, by dividing bone strength by peak stress. Values for bone strength were obtained using the linear relationship of Fyhrie and Vashishth, bone strength = 0.0061 EBone [37]. In this study we used 45 MPa for a experimentally obtained Young’s modulus of 7.3 GPa for finch beak bone [22,23]. Finally keratin thickness on top and bottom of the beak, and beak dimensions (depth and width relative to length) were measured and are reported for four functional groups (Table 3), and Pearson’s correlations (with Bonferroni correction) of these values compared to the SFs are shown in Table 4.

thumbnail
Table 3. Characterization of keratin thickness and beak dimensions (relative to beak length) in the different functional groups.

https://doi.org/10.1371/journal.pone.0129479.t003

thumbnail
Table 4. Pearson correlation and p-values between beak dimensions corrected for beak length and safety factors extracted from the different models.

https://doi.org/10.1371/journal.pone.0129479.t004

Results

The calculated von Mises stresses for the physiological boundary conditions in different Darwin’s finches are presented in Figs 2, 3, S2 and S3 (see also Table 2). The force vectors representing bite forces (Figs 2, S2 and S3) correspond well with field-measured in vivo bite forces (Table 5) [4, 5, 25]. Note, however, that in vivo bite forces reported here differ proportionally from those reported previously [25], as in that study the bite forces were not corrected for the lever arms of the bite force set-up. Among finches that crush seeds at the base of their beaks, safety factors varied between 0.7 and 2.5 (Table 5). Moreover, safety factors during natural loading regimes (i.e. base biting) generally decrease as bite force increases. Importantly, safety factors are lower during tip loading behaviors suggesting that beaks in these species are not optimized for tip biting. For other species, safety factors varied between 1 and 2.3 and were generally more similar between the two loading conditions.

thumbnail
Table 5. In vivo measured bite force (6 species, means ± standard deviations, N = number of specimens) compared with the model bite force and the model safety factors (SF).

https://doi.org/10.1371/journal.pone.0129479.t005

Density variation in bone can be expected and might result in a locally higher strength at high stress locations as suggested by our preliminary data. Indeed, the gray values (Table 2) at high vM stress locations were all higher than 0.5 (0 indicates the least dense bone in the sample, 1 indicates the densest bone). For these natural loading conditions, higher relative bone densities are thus observed at locations with higher stresses which can lead to increased safety factors.

The von Mises’ (vM) stresses for the 13 size-scaled models are presented in Figs 3, 4, S4 and S5 (also see Table 2). The maximum vM stress in the size-scaled models differs broadly across different species (Table 2). For example, the beaks of C. olivacea and P. inornata in our scaled models show highest peak vM stresses of up to 100 MPa, indicating that these shapes are not suited to withstand loading regimes similar to those observed for G. fortis. The deep and wide seed crushing beaks of G. magnirostris, G. fortis and G. fulignosa, show lower stresses under base loading (below 30 MPa) than under tip loading, consistent with these birds’ emphasis on base biting strategies in nature and again in accordance with our predictions. One can also observe moderate peak stresses for the tip biting beaks (P. crassirostris, C. pauper, C. psittacula, except C. parvulus), ranging from 31 MPa to 39 MPa. The probe and crush beaks (G. scandens, G. difficilis) show higher stresses, however, ranging from 43 to 51 MPa. Finally, comparatively low maximum stresses are found during tip loading for species that use the tips of their beak during foraging and have beaks with high curvature, as predicted (e.g. Camarhynchus pauper and C. psittacula, 35 MPa). The crushing beaks (G. fortis, G. fuliginosa, G. magnirostris) on the other hand have higher stresses (42 to 47 MPa, SF = ±1) under these loading conditions. The main exceptions to the otherwise rather precise tuning of beak strength and feeding mode concerns relatively high stresses observed for G. conirostris (67 MPa), C. parvulus (55 MPa), and Cactospiza pallida (80 MPa) under different loading regimes.

The rhamphotheca was thinnest in the probers (Table 3). Base crushers have a thicker keratin layer on the bottom of the beak than tip biters, and tip biters have thicker keratin on the top of the beak. The size scaled models show a strong correlation between the beak safety factor and the thickness of the keratin at the base of the beak (Table 4). It is also interesting to note that base and tip biters seem to have similar aspect ratios of the beak (depth and width relative to length), and that these ratios are lower in the probers (Table 3). Moreover, we observe a strong and positive correlation between aspect ratio and beak safety factor (Table 4).

Discussion

For all species, the estimated maximum vM stresses, particularly near the nasal hinge, were 4–44% lower than maximum vM stresses calculated in prior models that did not take into account the keratinous rhamphotheca [24]. An earlier study on Padda oryzivora [23] showed that changing the bone modulus value also had a linear effect, albeit small, on safety factor values (SFs changed from 2.5 to 3.0 if bone modulus changed from 6.7 to 7.9 GPa). Also, variation in the keratin modulus within its measured interval (Eker = 1.3 till 2.1 GPa) had a little effect on SF. However, if we decreased the keratin modulus below this, a profound effect was observed (for Eker of 1.7 GPa the SF = 2.8, while for Eker = 0.5 GPa the SF = 1.8), indicating the importance of incorporating the keratin layer in the models. Models that integrate information about multiple layers ideally consider well-defined material properties for each constituent material. These parameters were obtained on a different species of finch (Padda oryzivora) since experiments cannot be conducted with Darwin’s finches due to the protected status of these birds. Stress regimes for multi-layered beaks using these experimentally-obtained elastic moduli of keratin and bone had been modeled with success previously for Padda oryzivora and then validated using digital speckle pattern interferometry [22,23].

Keratin itself is also a multi-layered structure [38,39] that typically shows an anisotropic mechanical behavior, with preferential directions of failure depending on cell orientation [38]. The mechanical behavior of keratin also depends on its hydration state [22,40]. However, in living animals the hydration state of the keratin is controlled to the underlying living tissues (epidermis and dermis). In birds specifically, the beak consists of four distinct layers: a layer of dead epidermal cells (stratum corneum), a series of living epidermal cells, the dermis, and finally the bony beak [41]. Keratin cells typically orient such that they are aligned with the principal deformation [41]. In contrast to bone, which is a stiff mineralized tissue prone to fracture, keratin often shows local failure which does not penetrate the entire structure due to the organization of the cells in the different layers of the structure. Moreover, keratin is abrasion resistant and worn layers are shed and renewed by the addition of new layers of the stratum corneum through deposition from the living epidermal cells. As such the keratinous rhamphotheca is a continuously growing structure where damaged cells are shed and replaced. Consequently, we here decided to focus on the mechanical behavior of the underlying bony beak while taking into account the mechanical behavior of the keratinous rhamphoteca.

We found that the beaks of different finch species operate within a range of safety factors between 0.7 and 2.5 (Table 5). From an engineering perspective this range cannot be regarded as narrow, although from a biological perspective it is rather narrow if one takes into account the large variation finches express in bite forces (1N to 65N), beak sizes (1cm to 2cm), shapes, and biological uses. The in vivo realistic range of applicable safety factors is likely even narrower than this, given that the 0.7 value was calculated for a behavior rarely observed in nature: tip-biting in the largest species, G. magnirostris. A safety factor lower than 1 implies that the beak would break if the bird would use its full muscle strength while biting at the tip. This is, however, a biting strategy rarely observed in this base-biting bird (note that this species displays a safety factor of 1.1 for base loading). Moreover, bone remodeling could further increase safety factors by decreasing vM stresses in highly loaded areas, a statement supported by the higher grayscale values in these regions. In future modeling, density dependent elastic moduli and strength could be used to provide further insights in this material tuning. Nevertheless, calculated safety factors imply that beaks can withstand loading under normal conditions, although the exceptional jaw muscle hypertrophy of the ground finches [27], particularly G. magnirostris (S1 Fig), introduces risks to the beak’s structural integrity during tip biting.

As a further caveat, the magnitude of the safety factors given here should be interpreted as relative only when comparing different species, since our measures for the Young’s modulus were based on data for a different species (P. oryzivora). Moreover, we used the linear approximation of Fyhrie and Vashishth [37]. Indeed, a wide range of strengths for denser bone, ranging from 106 to 224 MPa [4244] is available in literature and Darwin’s finches could potentially have evolved denser bone with a higher strength. Moreover, stress magnitudes are sensitive to changes in material properties [23,4548]. Another assumption we make in calculating the safety factors is that the beaks of Darwin’s finches have similar material properties as those measured for P. oryzivora. We are, however, unable to test this directly since the preservation of the tissues can have a major impact on the mechanical properties of bone and keratin. In addition, freshly killed animals cannot be obtained.

A unique advantage of FE models, applied in our study, is that they can be size-scaled, thus allowing a size-independent view on the structural merits of different beak shapes. In our set of FE models, for which we scaled the beaks of all Darwin’s finches to the same size and applied identical muscle forces [26], our primary result is that maximum vM stresses differ broadly across different species, in ways that align with the different species feeding strategies (Table 2). For example, the beaks of C. olivacea and P. inornata show highest peak vM stresses in our models (up to 100 MPa). This suggests that their beak shape is not suited to crack very hard seeds. Indeed, these birds use their beaks almost exclusively to capture insects, rather than to crush or manipulate hard objects [49].

The high stresses for G. conirostris observed in our models might be explained by the fact that the geometry included in our study was that from a juvenile, while the applied muscle forces on this juvenile geometry were those from an adult. Likely, ontogenetic changes in shape and ossification take place to optimize the beak for adult loading conditions [50] as previously demonstrated for G. fortis. Results for C. parvulus and C. pallida remain puzzling, but suggest that these species may be functionally constrained in the use of their beak when applying large muscle forces.

We observed a thicker rhamphotheca in base crushing and tip-biting birds compared to birds that probe and bite, as well as dedicated probers (Table 3). Moreover, the keratin thickness at the bottom of the beak was strongly correlated to the beak safety factor indicating that thicker keratin helps protect the beak and increase its safety factor (Table 4). Interestingly, the thicker top keratin for tip biters matches the patterns of stress generated during tip loading, suggesting that in both ecotypes keratin is an important part of the stress mitigation strategy. In addition, positive and high correlations between beak aspect ratio and beak safety factor indicate that birds with relatively deeper and wider beaks have higher safety factors (Table 4) and are thus better equipped to withstand loading due to biting. This provides evidence that beak shape affects the risk of failure.

In summary, our FE models demonstrate that beak shapes are generally well-suited for mitigating risk of fracture in accordance with a species’ predominant feeding habitat. Understanding how beaks evolve to reduce risk of fracture may help to explain patterns of selection on beak size and shape in natural populations, and ultimately should be considered as one of the axes of adaptation and specialization in the Darwin’s finch radiation.

Data availability

Data is available from biomesh.org.

Supporting Information

S1 Fig. Muscle mass and tip bite force compared to the body mass of finches.

Note how G. magnirostris and G. fortis have an exceptionally high muscle mass and bite force for their size.

https://doi.org/10.1371/journal.pone.0129479.s001

(DOCX)

S2 Fig. Top and side view for physiological FE models of upper beak during base biting for 8 Darwin finches.

https://doi.org/10.1371/journal.pone.0129479.s002

(DOCX)

S3 Fig. Top and side view for physiological FE models of upper beak during tip biting for 8 Darwin finches.

https://doi.org/10.1371/journal.pone.0129479.s003

(DOCX)

S4 Fig. Top and side view for scaled FE models of upper beak during base biting for 13 Darwin finches.

https://doi.org/10.1371/journal.pone.0129479.s004

(DOCX)

S5 Fig. Top and side view for scaled FE models of upper beak during tip biting for 13 Darwin finches.

https://doi.org/10.1371/journal.pone.0129479.s005

(DOCX)

Acknowledgments

We thank S. Maas and J. Weis (University of Utah) for allowing us to use the FEBio software package. We also thank three anonymous reviewers for constructive comments on an earlier version of the manuscript. Field work was coordinated through the Charles Darwin Research Station and the Galápagos National Park Service (GNPS). We are particularly grateful for the generosity of the GNPS in granting a salvage permit, and to the California of Sciences for the loan of material. The authors thank Luis De Leon, Ana Gabela, Andrew Hendry, Mike Hendry, Eric Hilton, Sarah Huber, Katleen Huyghe, and Bieke Vanhooydonck for their assistance in the field. This work was performed in part at the Center for Nanoscale Systems (CNS), a member of the National Nanotechnology Infrastructure Network (NNIN).

Author Contributions

Conceived and designed the experiments: JS, PA, JD, AH, DA. Performed the experiments: JS, AG, PA, AH, JP. Analyzed the data: JS, AG, PA, AH. Wrote the paper: JS, AG, JP, PA, JD, AH, DA.

References

  1. 1. Temeles EJ, Kress WJ. Adaptation in a Plant-Hummingbird Association. Science. 2003; 300:630–3. pmid:12714743
  2. 2. Temeles EJ, Koulouris CR, Sander SE, Kress WJ. Effect of flower shape and size on foraging performance and trade-offs in a tropical hummingbird. Ecology. 2009; 90:1147–61. pmid:19537537
  3. 3. Boag PT, Grant PR. Intense Natural Selection in a Population of Darwin’s Finches (Geospizinae) in the Galapagos. Science. 1981;214:82–5. pmid:17802577
  4. 4. Herrel A, Podos J, Huber SK, Hendry AP. Evolution of bite force in Darwin’s finches: a key role for head width. J. Evol. Biol. 2005; 18:669–75. pmid:15842496
  5. 5. Herrel A, Podos J, Huber SK, Hendry AP. Bite performance and morphology in a population of Darwin’s finches: implications for the evolution of beak shape. Funct Ecol. 2005; 19:43–48.
  6. 6. van der Meij M, Bout R. Scaling of jaw muscle size and maximal bite force in finches. J Exp Biol. 2004; 207:2745–53. pmid:15235003
  7. 7. Bock WJ. Kinetics of the avian skull. J Morphol. 1964; 114:1–42.
  8. 8. Bout RG Biomechanics of the avian skull. In: Gasc JP, Casinos A, Bels VL, editors. Vertebrate Biomechanics and Evolution II. Oxford: BIOS Scientific Publishers Ltd; 2002. pp. 49–62.
  9. 9. van der Meij MAA, Bout RG. The relationship between shape of the skull and bite force in finches. J Exp Biol. 2008; 211:1668–80. pmid:18456894
  10. 10. Bowman RI. Morphological Differentation and Adaptation in the Galapagos Finches. University of california press: Berkeley and Los Angeles; 1961.
  11. 11. Nuijens F, Zweers F. Characters discriminating two seed husking mechanisms in finches (Fringillidae: Carduelinae) and Estrildids (Passeridae: Estrildinae). J Morphol.1997; 232:1–33.
  12. 12. Campàs O, Mallarino R, Herrel A, Abzhanov A, Brenner MP. Scaling and shear transformations capture beak shape variation in Darwin’s finches. PNAS. 2010; 107(8):3356–60. pmid:20160106
  13. 13. Herrel A, Podos J, Vanhooydonck B, Hendry AP. Force-velocity trade-off in Darwin's finch jaw function: a biomechanical basis for ecological speciation? Funct Ecol. 2009; 23:119–125.
  14. 14. Grant PR. Inheritance of size and shape in a population of Darwin’s finches. Proc R Soc Lond B. 1983; 212:403–432.
  15. 15. Grant PR. Ecology and Evolution of Darwin’s Finches. Princeton Univ: Princeton; 1999.
  16. 16. Grant PR, Grant RB. Evolution of character displacement in Darwin’s finches. Science. 2006; 313:224–6. pmid:16840700
  17. 17. Lack D. Darwin’s Finches. Cambridge Univ Press: Cambridge; 1947.
  18. 18. Foster D, Podos J, Hendry AP. A geometric morphometric appraisal of beak shape in Darwin’s finches. J Evol Biol. 2008; 21:263–75. pmid:18021202
  19. 19. Rayfield EJ. Finite element analysis and understanding the biomechanics and evolution of living and fossil organisms. Annu. Rev. Earth Planet. Sci. 2007; 35: 541–576.
  20. 20. Ross CF. Finite element analysis in vertebrate biomechanics. Anat Rec. 2005; 283A:253–8.
  21. 21. Richmond BG, Wright BW, Grosse I, Dechow PC, Ross CF, Spencer MA, et al. Finite element analysis in functional morphology. Anat Rec Part A. 2005; The Anatomical Record Part A: 283(2), 259–74. pmid:15747355
  22. 22. Soons J, Herrel A, Aerts P, Dirckx JJJ. Determination and validation of the elastic moduli of small and complex biological samples: bone and keratin in bird beaks. J R Soc Interface. 2012; https://doi.org/10.1098/rsif.2011.0667
  23. 23. Soons J, Herrel A, Genbrugge A, Adriaens D, Aerts P, Dirckx J. Multi-layered bird beaks: a finite-element approach towards the role of keratin in stress dissipation. J R Soc Interface. 2012; 9(73), 1787–96. pmid:22337628
  24. 24. Soons J, Herrel A, Genbrugge A, Aerts P, Podos J, Adriaens D, et al. Mechanical stress, fracture risk and beak evolution in Darwin’s ground finches (Geospiza). Philos T Roy Soc B. 2010; 365:1093–8. pmid:20194171
  25. 25. Herrel A, Soons J, Aerts P, Dirckx J, Boone M, Jacobs P, et al. Adaptation and function of the bills of Darwin's finches: divergence by feeding type and sex. Emu. 2010; 110(1), 39–47.
  26. 26. Dumont E R, Grosse I R, Slater G J. Requierements for comparing the performance of finite element models of biological structures. J Theor Biol. 2009; 256:96–103. pmid:18834892
  27. 27. Genbrugge A, Herrel A, Boone M, Van Hoorebeke L, Podos J, Dirckx J, et al. The head of the finch: the anatomy of the feeding system in two species of finches (Geospiza fortis and Padda oryzivora). J Anat. 2011; 219,6: 676–95. pmid:21999913
  28. 28. Loeb GE, Gans C. Electromyography for experimentalists. University of Chicago Press, Chicago; 1986.
  29. 29. Mendez J, Keys A. Density and composition of mammalian muscle. Metabolism. 1960; 9, 184–8
  30. 30. Nigg B, Herzog W. Biomechanics of the Musculo-skeletal System. John Wiley & Sons Ltd, University of Calgary, 3th edition; 2007. pmid:18001971
  31. 31. Nuijens FW, Zweers GA. Characters discriminating two seed husking mechanisms in finches (Fringillidae: Carduelinae) and estrildids (Passeridae: Estrildinae). J. Morp. 1997; 232(1):1–33
  32. 32. Si H. TetGen: A Quality Tetrahedral Mesh Generator and three-dimensional Delaunay Triangulator v1.4. Control at <http://tetgen.berlios.de>; 2006.
  33. 33. Maas S, Weiss JA. FEBio: Finite Elements for Biomechanics. User’s Manual, Version 1.0. Online publication; 2008.
  34. 34. Carter DR, Hayes WC. Bone compressive strength: the influence of density and strain rate. Science. 1976; 194:1174–6 pmid:996549
  35. 35. Rice JC, Cowin SC, Bowman JA. On the dependence of the elasticity and strength of cancellous bone on apparent density. J biomech. 1988; 21.2 155–68. pmid:3350829
  36. 36. Ford CM, Keaveny TM. The dependence of shear failure properties of trabecular bone on apparent density and trabecular orientation. J biomech. 1996; 29.10: 1309–17. pmid:8884476
  37. 37. Fyhrie DP, Vashishth D. Bone stiffness predicts strength similarly for human vertebral cancellous bone in compression and for cortical bone in tension. Bone. 2000; 26:169–73. pmid:10678412
  38. 38. Farren L, Shayler S, Ennos AR. The fracture properties and mechanical design of human fingernails. J Exp Biol. 2004; 207.5: 735–41.
  39. 39. Chou SF, Overfelt RA, Miller ME. Anisotropic mechanical behavior of keratin tissue from quill shells of North American porcupine (Erethizon dorsatum). Mater Sci Eng A-Struct. 2012; 557:36–44.
  40. 40. Li BW, Zhao HP, Feng XQ, Guo WW, Shan SC. Experimental study on the mechanical properties of the horn sheaths from cattle. J exp biol. 2010; 213(3):479–86. pmid:20086133
  41. 41. Genbrugge A, Adriaens D, Kegel B, Brabant L, Hoorebeke L, Podos J, et al. Structural tissue organization in the beak of Java and Darwin's finches. J anat. 2012; 221, 5: 383–93. pmid:22938039
  42. 42. Yamada H. Strength of Biological Materials. Williams & Wilkins: Baltimore; 1970.
  43. 43. Vogel S. Comparative Biomechanics: Life’s Physical World. Princeton University Press: Princeton; 2003.
  44. 44. Currey JD. Bones: Structure and Mechanics. Princeton University Press: Princeton; 2006.
  45. 45. Panagiotopoulou O, Kupczik K, Cobb SN. The mechanical function of the periodontal ligament in the macaque mandible: a validation and sensitivity study using finite element analysis. J. Anat. 2011; 218,75–86. pmid:20584094
  46. 46. Berthaume MA, Dechow PC, Iriarte-Diaz J, Ross CF, Strait DS, Wang Q, et al. Probabilistic finite element analysis of a craniofacial finite element model. j theor boil. 2012; 300: 242–53. pmid:22306513
  47. 47. Kupczik K, Dobson CA, Fagan MJ, Crompton RH, Oxnard CE, O’Higgins P. Assessing mechanical function of the zygomatic region in macaques: validation and sensitivity testing of finite element models. J Anat. 2007; 210,1: 41–53. pmid:17229282
  48. 48. Fitton LC, Proa M, Rowland C, Toro-ibacache V, O’Higgins P. The impact of simplifications on the performance of a finite element model of a Macaca fascicularis cranium. The Anat Rec. 2015; 298,1: 107–21. pmid:25339306
  49. 49. Grant P R, Grant B R. How and why Species Multiply: The radiation of Darwin’s Finches. Princeton University Press: Princeton & Oxford; 2008.
  50. 50. Genbrugge A, Heyde AS, Adriaens D, Boone M, Van Hoorebeke L, Dirckx J, et al. Ontogeny of the cranial skeleton in a Darwin’s finch (Geospiza fortis). J anat. 2011; 219, 2: 115–31. pmid:21599660