Next Article in Journal
Bioenergetic Pathways in the Sperm of an Under-Ice Spawning Fish, Burbot (Lota lota): The Role of Mitochondrial Respiration in a Varying Thermal Environment
Previous Article in Journal
Aldehyde Dehydrogenase Mutation Exacerbated High-Fat-Diet-Induced Nonalcoholic Fatty Liver Disease with Gut Microbiota Remodeling in Male Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Egg Incubation Mechanics of Giant Birds

1
Department of Mechanical Engineering, National Taiwan University, Taipei 10617, Taiwan
2
Department of Life Science, National Taiwan University, Taipei 10617, Taiwan
3
Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
4
Degree Program of Genome and Systems Biology, National Taiwan University, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 30 June 2021 / Revised: 28 July 2021 / Accepted: 29 July 2021 / Published: 1 August 2021
(This article belongs to the Section Biophysics)

Abstract

:

Simple Summary

Extinct giant birds have been a source of imagination, and knowledge of their incubation mechanics is crucial to our understanding of the evolution of avian reproduction. Despite the extensive studies on avian eggs, our understanding of the eggshell mechanics of giant birds, particularly the extinct ones, remains incomplete—most of these prior works were based on empirical or allometric relationships with limited quantitative analysis. In the present study, with the help of advanced three-dimensional computer simulation using data from published fossil records, we obtain more comprehensive quantitative analysis to answer important questions related to contact incubation of giant birds. Specifically, how much safety margin does the reversed sexual size dimorphism (RSSD) of moas provide? What is the theoretical upper limit of body mass for contact incubation? Is the Williams’ egg, or the putative Genyornis oological material (PGOM), really the egg of the extinct giant bird Genyornis newtoni, as commonly accepted since its discovery in 1981?

Abstract

Finite element analysis (FEA) was used to conduct mechanical analyses on eggshells of giant birds, and relate this to the evolution and reproductive behavior of avian species. We aim to (1) investigate mechanical characteristics of eggshell structures of various ratite species, enabling comparisons between species with or without reversed sexual size dimorphism (RSSD); (2) quantify the safety margin provided by RSSD; (3) determine whether the Williams’ egg can have been incubated by an extinct giant bird Genyornis newtoni; (4) determine the theoretical maximum body mass for contact incubation. We use a dimensionless number C to quantify relative shell stiffness with respect to the egg size, allowing for comparison across wide body masses. We find that RSSD in moas significantly increases the safety margin of contact incubation by the lighter males. However, their safety margins are still smaller than those of the moa species without RSSD. Two different strategies were adopted by giant birds—one is RSSD and thinner shells, represented by some moa species; the other is no RSSD and regular shells, represented by the giant elephant bird. Finally, we predicted that the upper limit of body mass for contact incubation was 2000 kg.

Graphical Abstract

1. Introduction

Birds have arguably been very successful in evolution [1]. In addition to the extraordinary capability of flying, they also have a special and efficient way of reproduction [2]. Most birds adopt bird-egg contact incubation and “egg turning” to maintain an adequate environment for embryonic development [2]. During such a process, eggs are subjected to the weight of the incubating bird and possible impact between eggs. Thus, paradoxically, an eggshell has to be robust enough to withstand the weight of its parent bird during incubation while also breakable for the chick to emerge; these, apparently, contradictory demands suggest an optimal design. In our previous work, we proposed a dimensionless metric, C number, to characterize egg stiffness with respect to egg mass. This metric facilitates the comparison of eggshell stiffness across a wide range of body mass [3]. Based on the data of 463 bird species in 36 orders across five orders of magnitude in body mass, we found that C number is nearly invariant for most species, including tiny hummingbirds and giant elephant birds [3]. However, that study only has limited results on giant birds, despite the fact that incubation mechanics of giant birds is crucial to our understanding of the evolution of avian reproduction.
Here, we use published data to study the incubation mechanics of ratites and an extinct giant bird Genyornis newtoni—a Galloansere. Ratites are large, flightless birds that are essential to understanding the early evolution of birds. The origin of the word “ratite” is the Latin word ratis (raft); these birds obtained their name because of their common character, a flat (raft-like) sternum without a keel to anchor wing muscles, which are necessary for flying. Ostrich, the world’s largest living bird, can be taller than a human and weigh up to 150 kg. In addition, based on fossil evidence, even larger ratite species once lived on earth. For example, the extinct elephant birds from Madagascar, the heaviest birds known to have lived, weighed up to 450 kg for the well-known Aepyornis maximus [4,5] and up to 650 kg for a recently recognized taxon Vorombe titan [6]. Both extant ratites (ostrich in Africa, rheas in South America, kiwi in New Zealand, and emu and cassowaries in Australia) and extinct ratites (moa in New Zealand and elephant bird in Madagascar) are separated from one another by oceans Figure 1. That ratites are flightless but geographically widespread has induced many discussions about their evolutionary relationships [7].
According to the empirical relationship between eggshell thickness and body mass [3,8], larger birds tend to produce thicker eggshells. Therefore, the dilemma of eggshells from large-size birds may be more notable. Worthy and Holdaway [9] and Worthy et al. [10] used bone material to estimate the body mass of moa and calculated a few species reached >200 kg. Next, based on sex chromosome evidence, some moa species had reversed sexual size dimorphism (RSSD), i.e., females were larger than males [11]. For extant ratites with RSSD, e.g., emu and cassowaries, lighter males take responsibility for incubating eggs. Similarly, some researchers [8] argued that moa eggshells were not capable of withstanding the weight of the heavier sex (female), so their eggs were solely incubated by the lighter sex (male). Extreme RSSD in moa species included the strategy of compensating the mechanical discrepancy of eggshells. However, Huynen et al. [12] suggested that males of larger moa still could not perform contact incubation like extant birds. In contrast to the giant elephant bird (Aepyornis maximus), which had a larger size but no RSSD [6], what might be the underlying reasons for distinct mechanisms between RSSD and non-RSSD species?
Another special case related to the incubation mechanics of giant birds is the iconic Australian Genyornis newtoni and the issue of its associated eggshell. Genyornis newtoni is a Galloansere with an estimated mass of 275 kg [13]. Fossil eggshell material, known as the “Williams’ egg,” was attributed to G. newtoni in 1981 [14] and had been commonly accepted for many years. However, Grellet-Tinner et al. recently conducted comprehensive measurement and analysis on fossil eggshell, and concluded that this “Genyornis” eggshell, or the putative Genyornis oological material (PGOM), is unexpectedly small given the size of G. newtoni. They suggest that PGOM is more likely to have been laid by the giant extinct Progura, a genus of extinct giant megapodes. Other giant birds had also lived on earth but are now extinct. Examples include the dromornithids, gastornithids, and phorusrhacids, with some giants exceeding 500 kg [15]. Those are not included in the present study due to lack of eggshell fossils.
The objective of the present study is to use a computer simulation technique called the finite element analysis (FEA) to (1) investigate mechanical characteristics of eggshell structures of various ratite species, enabling comparisons between species with or without RSSD; (2) quantify the safety margin provided by RSSD; (3) determine whether the Williams’ egg can have been incubated by an extinct giant bird Genyornis newtoni; (4) determine the theoretical maximum body mass for contact incubation. Answering these questions are important to our understanding of the evolution and reproductive behavior of birds.
FEA is conventionally used by engineers to design and analyze man-made structures or products, and has recently been applied to biomechanical problems in animals because of its powerful capability in predicting structure–function relationship [16]. For example, FEA has been used to infer the biting-induced stress distribution in skulls for moas (Dinornithiformes) [17], the sauropod taxa Camarasaurus and Diplodocus [18], a large carnivorous theropod dinosaur Allosaurus fragilis [16], the American sabercat Smilodon fataqlis [19], and Neanderthals (Homo neanderthalensis) and modern humans [20]. In our previous work, we applied three-dimensional FEA to study the eggshell mechanics of 430 avian species (36 orders, 104 families); a side-by-side comparison of eggshell stiffness of several representative species between FEA, experiments, and shell theory confirmed the accuracy of FEA [3]. This experimentally validated modelling approach was adopted here to analyze the eggshell mechanics of extinct giant birds.

2. Methods

Data required for our analyses, including (1) body mass, (2) egg mass, and (3) eggshell geometry (including length, width, and thickness), were obtained from the literature Table 1. Egg mass was estimated by Hoyt’s equation [21]. To increase the reliability of biological data, we gathered information from multiple sources. The eggshell geometry data of extant species was obtained from Dickison [22] and Juang et al. [3]. For extinct species, data for moa were obtained from Gill [23,24] and Worthy and Holdaway [9], whereas data for the elephant bird were from Mlíkovsky [25] and Grealy et al. [26]. Regarding body mass, we obtained not only average mass of the species but the range of both sexes, consulting three books [27,28,29], and materials from Hansford and Turvey [6], Bunce et al. [30], and Huynen et al. [12]. Note that populations of Dinornis (Dinornis robustus and Dinornis novaezealandiae) have markedly different size parameters for sexes [10] so our analyses are based on the correct association of egg and most relevant population data (Table 1). Population-level body mass was also used for another moa species Emeus crassus to ensure that its relatively low safety margin is not resulted from using unmatched species level data [23,24,31,32]. In total, data for nine extant and nine extinct avian species were collected. The data used for analyzing G. newtoni and PGOM were obtained from references [13,14,33,34], Table 2.

2.1. Finite Element Analysis (FEA)

We used a finite element package ANSYS to simulate the quasi-static compression on an eggshell model along the major (long) axis, which has higher geometric stiffness [35]. Recall that the stiffness, KF/δ, of a load-bearing body relates to its ability to resist deformation, where F and δ are respectively the applied load and the corresponding deformation Figure 2. First, we constructed a three-dimensional eggshell model and meshed it with 4-node structural shell elements (SHELL181). The convergent test was performed to ensure that the element size was sufficiently small. We used mesh refinement around the top near where local deformation will occur due to the 0.1-N point force applied on top of the eggshell model, causing a microscopic deformation by which we obtained stiffness.
For simplicity, we assumed a perfect ellipse profile and uniform thickness. To construct an eggshell model, we used only: (1) polar diameter, (2) equatorial diameter, and (3) thickness. To save calculation time, we only created the upper half of the model and fixed the bottom ring as the boundary condition. Thereafter, the stiffness of the entire model can be obtained by dividing that of the upper-half model by two. The effect of egg content is negligible and is not included in the model [36]. The shell model is assumed linearly elastic, homogeneous and isotropic, with elastic modulus 30 GPa [3] and Poisson’s ratio 0.307 [37].

2.2. Dimensionless Number C, Critical Thickness, and the Factor of Safety

We used the dimensionless number defined by Juang et al. [3] as
C A 2 K B W
with the length of major axis and minor axis A and B, respectively (unit: m), the weight of egg W in N, and stiffness K in N m−1 calculated from the FEA result. In an extreme scenario of the compression test, the eggshell is subjected to a point force equivalent to the body weight of its parent bird. By adjusting the shell thickness, the minimum thickness, referred to as the critical thickness tcr, that prevents the eggshell from buckling can be obtained (see Supplementary Figure S8 in reference [3]). We then created another eggshell model with the critical thickness, conducted the 0.1-N compression simulation again and obtained the critical dimensionless number. Note that C represents relative stiffness with respect to the egg size, which enables a direct comparison across a wide range of body mass. By contrast, the absolute stiffness K cannot be easily used for this purpose. For example, the stiffness of a massive elephant bird egg (Aepyornis maximus, K ≈ 6605 N mm−1) is much larger than that of a hummingbird egg (Mellisuga minima, K ≈ 11 N mm−1), but is it large enough to resist the incubating bird’s weight, as the elephant bird is also much heavier than the hummingbird? Interestingly, their C numbers are comparable (elephant bird ≈ 12,000 versus hummingbird ≈ 17,500) and similar to the average value, ~15,000, of 463 species [3]. Finally, a factor of safety, defined as the ratio of C to critical Ccr indicates the breakability of the eggshell [3].
F . S . C C c r
The higher the F.S., the less likely the eggshell is to break (F.S. < 1 indicates that the parent bird cannot contact-incubate the egg without breaking it).

3. Results and Discussion

The dimensionless number C can be interpreted as a measure of the stiffness of eggshell, with the geometry-induced stiffness eliminated and reflecting only contributions of shell thickness and material properties. In contrast, the critical C (or the further induced F.S.) is based on structural stability.
Scatter log-log plots for C number versus body mass are illustrated in Figure 3. Juang et al. [3] reported an averaged C of ~15,000 for 463 avian species and a slightly decreasing allometric trend of C with increasing body mass. Based on our results, we inferred that most ratites have a relatively lower C number, with ostrich (Struthio camelus) being an exception. Therefore, if we combined the data in Figure 3a, the slope of the blue fitting line will be slightly increased Figure 3b.
Figure 4 shows F.S., C number, and the phylogenetic relationships of 18 ratite species, nine of which are extinct, including eight moa and one elephant bird species. Our previous study [3] shows that all of the 463 species analyzed have a F.S. > 2, which we adopt here as a criterion to determine the feasibility of contact incubation for extinct species. Moreover, F.S. in the present work was analyzed separately for the two sexes of RSSD and SSR species to investigate the feasibility of contact incubation by each sex from the mechanical point of view. Two rhea species (Rhea americana and Rhea pennata) lack sexual size dimorphism (SSD) [27], as do three moa species (Anomalopteryx didiformis, Megalapteryx didinus, and Pachyornis elephantopus) and elephant bird [6]. The above species were represented as green bars in the figure. Only ostrich has SSD, whereas all 11 other species are RSSD [27,38]. Costal moa (Euryapteryx curtus) were classified into two subspecies [39], with Euryapteryx curtus curtus from the North Island smaller than Euryapteryx curtus gravis from the South Island; they were analyzed and presented separately.
An investigation of the results highlights the following observations: First, the F.S. values for the non-RSSD and non-SSD species (green bars) are mostly greater than two, indicating that they can contact-incubate their eggs. This is supported by the fact that two of them are extant rhea species and indeed perform contact incubation. Second, species with RSSD incubated by the lighter sex indeed provide a greater safety margin. For example, in the case of emu and the cassowaries the lighter males incubate their eggs with a safe F.S. > 2; the F.S. would be <2 (marginal) if the incubating sex were the heavier female. However, three moa species (Dinornis robustus, Dinornis novaezealandiae, and Emeus crassus) had F.S. < 2, even if being incubated by the lighter sex (male). Their male F.S was categorized as marginal, which means parent birds could perform contact incubation, but it was approaching unsafe. Furthermore, the presence of male-specific DNA on an outer eggshell supports the notion that for larger moas, their eggs were incubated by males [12]—care should be taken when interpreting this results as the shells have to be well etched before analysis to remove contaminants. In addition, based on the allometric relationship between eggshell thickness and body mass Figure 5b, the shells of the above-mentioned three moa species (Dinornis robustus, Dinornis novaezealandiae, and Emeus crassus) appear too thin [8]. For example, the estimated eggshell thickness of South Island giant moa (Dinornis robustus) ranged from 1.31 to 2.16 mm. The possibility that weathering thin the fossilized shells of those three moas was ruled out based on careful examinations of the shell cross-section of a South Island giant moa [43]; it was concluded that its shell was indeed thinner than the trend line predicted. As a result, those three moa species have C numbers smaller than the average of ~15,000 Figure 4c. Generally, large-size moas produced large eggs Figure 5a. If the egg was incubated by lighter males, the maternal moa could produce thinner eggshells to conserve energy. By contrast, extinct taxa with no RSSD, e.g., the elephant bird, had shell thickness 3.7 mm and C number ~12,000 that lie on the trend lines Figure 5b; these species had to produce stiffer eggs to withstand the weight of both parents, as compared to the RSSD moas. The elephant bird produced the largest egg (303 by 224 mm) with the thickest shell known, even larger than those of non-avian dinosaurs. Note that although ostriches exhibit SSD, both male and female birds contact-incubate their eggs, and the shell thickness, C number, and F.S. also lie on the trend lines. Thus, two different reproduction strategies have been adopted by giant birds—one is RSSD and thinner shells, represented by some moa species, emus and cassowaries; the other is no RSSD and shells with regular thickness, represented by the giant elephant bird and ostriches—the thickness of their shells follows the allometric trend in Figure 5b.

3.1. Kiwi and PGOM

Kiwis were excluded from the fitting line due to their extreme reproductive characteristics. The kiwi is similar in height to a domestic chicken and weighs ~2 kg on average [28,44]. However, they lay an egg with approximately one-quarter of their body mass, whereas for most other avian species, the ratio of egg mass to body mass is <5%. Based on our calculations, their F.S. is around 2 to 5, even higher than some of the other ratites Figure 4b; they should have no problem contact-incubating eggs. However, since egg weight W is in the denominator of the C number definition, the kiwi’s large egg results in a much smaller C number, leading to an excessive deviation from the allometric trend Figure 4c.
Another special case herein is an extinct Australian avian species, Genyornis newtoni, which exhibits SSD with the smaller and larger sexes having an estimated mean mass of 192.8 kg (range 180–203 kg) and 238.4 kg (range 214–262 kg), respectively [33]. Grellet-Tinner et al. [13] suggested that an extinct large megapode, which belonged to the genus Progura, was more likely to be the parent bird of the putative Genyornis oological material (PGOM, to use their abbreviation), which was assigned to Genyornis by Williams [14], a conclusion accepted for decades. Shute et al. [34] taxonomically described/revised several large extinct megapodes from Australia. The species of Progura mentioned by Grellet-Tinner et al. [13] were classified as separate species called Latagallina naracoortensis and Progura gallinacea. In the calculation of C number and F.S. for PGOM, we used the two above-mentioned species and Genyornis newtoni, then computed with Williams’ egg [14] and the egg of Grellet-Tinner et al. [13], the so-called Spooner egg (126 by 97 mm). In addition, Grellet-Tinner et al. suggested that the size of the complete shell of Williams’ egg (155 by 125 mm), estimated based on shell curvature, may have been overestimated [33]. They considered that Williams’ egg might be the same length and width as the Spooner egg. Therefore, we additionally simulated the specimen with length and width of the Spooner egg but using the shell thickness of the Williams’ egg. Figure 3d shows the C numbers of various egg-bird combinations. The combination of Williams’ egg and Genyornis newtoni lies on the allometric trend line, whereas the combinations of the Spooner egg and the two megapodes were consistent with the trend. However, the F.S. values of Genyornis newtoni are <1, suggesting that neither the Williams’ egg nor the Spooner egg could be contact incubated by Genyornis newtoni. On the other hand, the F.S. values of the two extinct megapodes are >2 and are even higher than those of some extant counterparts Figure 6, which suggests that megapodes may have evolved a less stiff eggshell to adapt to their new environment during the Pleistocene era. Accordingly, we conclude that the Williams’ egg is unlikely associated with Genyornis newtoni and is likely to be the egg of extinct megapodes with a dimension similar to the Spooner egg’s. Table 3 lists eggs of extant birds that have similar size to the Spooner egg for comparison.

3.2. Maximum Body Mass for Contact Incubation

When considering the maximum body mass of avian species, the egg is assuredly an important factor [5]. Deeming and Birchard [45] proposed that contact incubation prevents avian species from attaining great sizes. In addition, they also explained that extinct “gigantic” birds were not truly gigantic, compared to mammals and non-avian dinosaurs. Finally, they proposed a maximum mass of 500 kg.
In this study, we estimate the maximum avian body mass by considering the body mass at which the allometric line of F.S. becomes smaller than one—contact incubation is no longer possible. Recall that F.S. ≡ C/Ccr and from Figure 3b we observe that the two solid fitting lines for C and Ccr intersect at ~2000 kg, which may be regarded as an estimate of the theoretical upper limit of body mass for contact incubation. This estimate, of course, may vary within a certain range according to the dash lines of 95% confidence interval. Note that the critical thickness was defined by considering the failure mode of buckling instead of fracture. The former is easier to predict, whereas the latter is sensitive to micro-cracks and thus not easy to determine [3,36]. Since the buckling force is generally greater than the fracture force, i.e., an egg under compression fractures before buckling, our upper limit was assuredly high enough to be a valid estimate. It is noteworthy that 2000 kg is well above the maximum estimated mass for the largest known birds, e.g., Vorombe titan, ~650 kg. Such a discrepancy might be explained by the fact that our analysis only considers whether an eggshell is strong enough to sustain the weight of the incubating bird. Another key factor that constrains the maximum body size—the hatching process—is not considered in the present study. Larger eggs have thicker shells, and result in a more challenging and prolonged hatching process. For example, it can take several hours for an ostrich chick to hatch naturally. In other words, the allowable body mass (and thus the maximum egg size) for contact incubation, constrained by hatching mechanics, is likely lower than 2000 kg. As an interesting analogy, gigantic size of 2000 kg or heavier is not rare and has been achieved by extinct and extant mammals, as well as many theropod dinosaurs—close relatives of birds. One reason they can evolve such gigantic sizes, compared to birds, is that they are not constrained by mechanics of contact incubation.

4. Conclusions

Size is a critical factor in evolution. Although extensive research on allometry related to avian reproduction has been performed, our research provides mechanical analyses on eggshells of both extant and extinct ratites and galloanseres. Based on our data, RSSD in moas significantly increases the safety margin of contact incubation if the incubation is indeed conducted solely by the lighter males. However, their safety margins are still, in general, smaller than those of the moa species without RSSD. Thus, two different reproduction strategies have been adopted by giant birds—one is RSSD and thinner shells, represented by some moa species, emus, and cassowaries; the other is no RSSD and shells with regular thickness, represented by the giant elephant bird and ostriches. In addition, our FEA shows that the Williams’ egg could not be contact incubated by Genyornis newtoni and is more likely associated with an extinct species of megapodes. Finally, we proposed a theoretical upper limit of body mass, ~2000 kg, for contact incubation based on numerical simulation (F.S. = 1). One of the most important obstacles of investigating extinct species is that materials are insufficient since fossil formation and discovery are both incidental. Regardless, our model can be used as a tool to analyze additional fossil evidence, obtained either by using DNA sequences to identify those eggshell fragments previously unidentified or by discovering new eggshell fragments accompanied by their respective adult skeleton, and thereby strengthen the allometric relationship. The mechanisms of evolution are inherently complex. In that regard, specific species might deviate from the allometric trend to some extent due to particular adaptations (e.g., kiwis, as mentioned above). However, if we assess a large range across body mass of many orders of magnitude, the interspecific allometric trend indeed provides some evolutionary insights regarding reproductive strategies of oviparous animals laying eggs protected by hard shells.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biology10080738/s1, Table S1: copyright information about bird images.

Author Contributions

Conceptualization, A.Y., P.-Y.C. and J.-Y.J.; methodology, A.Y., P.-Y.C. and J.-Y.J.; software, A.Y., H.-J.W., P.-Y.C. and J.-Y.J.; validation, A.Y. and H.-J.W.; formal analysis, A.Y. and H.-J.W.; investigation, A.Y. and H.-J.W.; resources, J.-Y.J.; data curation, A.Y. and H.-J.W.; writing—original draft preparation, A.Y.; writing—review and editing, H.-J.W., H.-T.Y. and J.-Y.J.; visualization, A.Y. and H.-J.W.; supervision, H.-J.W. and J.-Y.J.; project administration, J.-Y.J.; funding acquisition, J.-Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology (MOST) of Taiwan and National Taiwan University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank John Kastelic for English editing this manuscript, Yi-Jung Lin (Museum of Zoology, National Taiwan University) and Chi-Ya Lee (National Museum of Natural Science) for help on sample collections, and Trevor Worthy and two other anonymous reviewers for offering their comments and suggestions, which greatly improve the scientific quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Brusatte, S.L.; O’Connor, J.K.; Jarvis, E.D. The Origin and Diversification of Birds. Curr. Biol. 2015, 25, R888–R898. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Deeming, D.C. Avian Incubation: Behaviour, Environment, and Evolution; Oxford University Press: Oxford, UK, 2002; p. xiv. 421p. [Google Scholar]
  3. Juang, J.Y.; Chen, P.Y.; Yang, D.C.; Wu, S.P.; Yen, A.; Hsieh, H.I. The avian egg exhibits general allometric invariances in mechanical design. Sci. Rep. 2017, 7, 14205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Dawkins, R.; Wong, Y. The Ancestor’s Tale: A Pilgrimage to the Dawn of Evolution, 2nd ed.; Houghton Mifflin Harcourt: Boston, MA, USA, 2016; p. xxi. 771p. [Google Scholar]
  5. Murray, P.F.; Vickers-Rich, P. Magnificent Mihirungs: The Colossal Flightless Birds of the Australian Dreamtime; Indiana University Press: Bloomington, IN, USA, 2004. [Google Scholar]
  6. Hansford, J.P.; Turvey, S.T. Unexpected diversity within the extinct elephant birds (Aves: Aepyornithidae) and a new identity for the world’s largest bird. R. Soc. Open Sci. 2018, 5, 181295. [Google Scholar] [CrossRef] [Green Version]
  7. Mitchell, K.J.; Llamas, B.; Soubrier, J.; Rawlence, N.J.; Worthy, T.H.; Wood, J.; Lee, M.S.; Cooper, A. Ancient DNA reveals elephant birds and kiwi are sister taxa and clarifies ratite bird evolution. Science 2014, 344, 898–900. [Google Scholar] [CrossRef] [Green Version]
  8. Birchard, G.F.; Deeming, D.C. Avian eggshell thickness: Scaling and maximum body mass in birds. J. Zool. 2009, 279, 95–101. [Google Scholar] [CrossRef]
  9. Worthy, T.H.; Holdaway, R.N. The Lost World of the Moa: Prehistoric Life of New Zealand, 1st ed.; Indiana University Press: Bloomington, IN, USA, 2002. [Google Scholar]
  10. Worthy, T.; Bunce, M.; Cooper, A.; Scofield, P. Dinornis—An insular oddity, a taxonomic conundrum reviewed. Monogr. Soc. Hist. Nat. Balear. 2005, 12, 377–390. [Google Scholar]
  11. Bunce, M.; Worthy, T.H.T.H.; Ford, T.; Hoppitt, W.; Willerslev, E.; Drummond, A.; Cooper, A. Extreme reversed sexual size dimorphism in the extinct New Zealand moa Dinornis. Nature 2003, 425, 172–175. [Google Scholar] [CrossRef] [PubMed]
  12. Huynen, L.; Gill, B.J.; Millar, C.D.; Lambert, D.M. Ancient DNA reveals extreme egg morphology and nesting behavior in New Zealand’s extinct moa. Proc. Natl. Acad. Sci. USA 2010, 107, 16201–16206. [Google Scholar] [CrossRef] [Green Version]
  13. Grellet-Tinner, G.; Spooner, N.A.; Worthy, T.H. Is the “Genyornis” egg of a mihirung or another extinct bird from the Australian dreamtime? Quat. Sci. Rev. 2016, 133, 147–164. [Google Scholar] [CrossRef] [Green Version]
  14. Williams, D.L.G. Genyornis Eggshell (Dromornithidae, Aves) from the Late Pleistocene of South-Australia. Alcheringa 1981, 5, 133–140. [Google Scholar] [CrossRef]
  15. Worthy, T.H.; Degrange, F.J.; Handley, W.D.; Lee, M.S.Y. The evolution of giant flightless birds and novel phylogenetic relationships for extinct fowl (Aves, Galloanseres). R. Soc. Open Sci. 2017, 4, 170975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Rayfield, E.J.; Norman, D.B.; Horner, C.C.; Horner, J.R.; Smith, P.M.; Thomason, J.J.; Upchurch, P. Cranial design and function in a large theropod dinosaur. Nature 2001, 409, 1033–1037. [Google Scholar] [CrossRef] [PubMed]
  17. Attard, M.R.G.; Wilson, L.A.B.; Worthy, T.H.; Scofield, P.; Johnston, P.; Parr, W.C.H.; Wroe, S. Moa diet fits the bill: Virtual reconstruction incorporating mummified remains and prediction of biomechanical performance in avian giants. Proc. R. Soc. B 2016, 283. [Google Scholar] [CrossRef] [Green Version]
  18. Button, D.J.; Rayfield, E.J.; Barrett, P.M. Cranial biomechanics underpins high sauropod diversity in resource-poor environments. Proc. R. Soc. B 2014, 281, 1–9. [Google Scholar] [CrossRef] [PubMed]
  19. McHenry, C.R.; Wroe, S.; Clausen, P.D.; Moreno, K.; Cunningham, E. Supermodeled sabercat, predatory behavior in Smilodon fatalis revealed by high-resolution 3D computer simulation. Proc. Natl. Acad. Sci. USA 2007, 104, 16010–16015. [Google Scholar] [CrossRef] [Green Version]
  20. Wroe, S.; Parr, W.C.H.; Ledogar, J.A.; Bourke, J.; Evans, S.P.; Fiorenza, L.; Benazzi, S.; Hublin, J.-J.; Stringer, C.; Kullmer, O.; et al. Computer simulations show that Neanderthal facial morphology represents adaptation to cold and high energy demands, but not heavy biting. Proc. R. Soc. B 2018, 285. [Google Scholar] [CrossRef] [PubMed]
  21. Hoyt, D.F. Practical Methods of Estimating Volume and Fresh Weight of Bird Eggs. Auk 1979, 96, 73–77. [Google Scholar] [CrossRef]
  22. Dickison, M.R. The Allometry of Giant Flightless Birds. Ph.D. Thesis, Duke University, Durham, NC, USA, 2007. [Google Scholar]
  23. Gill, B.J. Eggshell characteristics of moa eggs (Aves: Dinornithiformes). J. R. Soc. N. Z. 2007, 37, 139–150. [Google Scholar] [CrossRef]
  24. Gill, B.J. A catalogue of moa eggs (Aves: Dinornithiformes). Rec. Auckl. Mus. 2006, 43, 55–80. [Google Scholar]
  25. Mlíkovsky, J. Eggs of extinct aepyornithids (Aves Aepyornithidae) of Madagascar: Size and taxonomic identity. Sylvia 2003, 39, 133–138. [Google Scholar]
  26. Grealy, A.; Phillips, M.; Miller, G.; Gilbert, M.T.P.; Rouillard, J.-M.; Lambert, D.; Bunce, M.; Haile, J. Eggshell palaeogenomics: Palaeognath evolutionary history revealed through ancient nuclear and mitochondrial DNA from Madagascan elephant bird (Aepyornis sp.) eggshell. Mol. Phylogenetics Evol. 2017, 109, 151–163. [Google Scholar] [CrossRef]
  27. Davies, S.J.J.F.; Bamford, M. Ratites and Tinamous: Tinamidae, Rheidae, Dromaiidae, Casuariidae, Apterygidae, Struthionidae; Oxford University Press: Oxford, UK, 2002; p. xxiii. 310p. [Google Scholar]
  28. Grzimek, B. Grzimek’s Animal Life Encyclopedia 8 Birds I Tinamous and Ratites to Hoatzins, 2nd ed.; Gale: Farmington Hills, MI, USA, 2003; Volume 8. [Google Scholar]
  29. Hume, J.P.; Walters, M. Extinct Birds; Bloomsbury Publishing: London, UK, 2012. [Google Scholar]
  30. Bunce, M.; Worthy, T.H.; Phillips, M.J.; Holdaway, R.N.; Willerslev, E.; Haile, J.; Shapiro, B.; Scofield, R.P.; Drummond, A.; Kamp, P.J.J.; et al. The evolutionary history of the extinct ratite moa and New Zealand Neogene paleogeography. Proc. Natl. Acad. Sci. USA 2009, 106, 20646–20651. [Google Scholar] [CrossRef] [Green Version]
  31. Holdaway, R.N.; Worthy, T.H. A reappraisal of the late quaternary fossil vertebrates of Pyramid Valley Swamp, North Canterbury, New Zealand. N. Z. J. Zool. 1997, 24, 69–121. [Google Scholar] [CrossRef]
  32. Huynen, L.; Millar, C.D.; Scofield, R.P.; Lambert, D.M. Nuclear DNA sequences detect species limits in ancient moa. Nature 2003, 425, 175–178. [Google Scholar] [CrossRef] [PubMed]
  33. Grellet-Tinner, G.; Spooner, N.A.; Handley, W.D.; Worthy, T.H. The Genyornis Egg: Response to Miller et al.’s commentary on Grellet-Tinner et al., 2016. Quat. Sci. Rev. 2017, 161, 128–133. [Google Scholar] [CrossRef]
  34. Shute, E.; Prideaux, G.J.; Worthy, T.H. Taxonomic review of the late Cenozoic megapodes (Galliformes: Megapodiidae) of Australia. R. Soc. Open Sci. 2017, 4. [Google Scholar] [CrossRef] [Green Version]
  35. Lazarus, A.; Florijn, H.C.B.; Reis, P.M. Geometry-Induced Rigidity in Nonspherical Pressurized Elastic Shells. Phys. Rev. Lett. 2012, 109, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hahn, E.N.; Sherman, V.R.; Pissarenko, A.; Rohrbach, S.D.; Fernandes, D.J.; Meyers, M.A. Nature’s technical ceramic: The avian eggshell. J. R. Soc. Interface 2017, 14, 20160804. [Google Scholar] [CrossRef] [Green Version]
  37. Kemps, B.; De Ketelaere, B.; Bamelis, F.; Govaerts, T.; Mertens, K.; Kamers, B.; Tona, K.; Decuypere, E.; De Baerdemaeker, J. Development of a Methodology for the Calculation of Young’s Modulus of Eggshell using Vibration Measurements. Biosyst. Eng. 2004, 89, 215–221. [Google Scholar] [CrossRef]
  38. Angst, D.; Buffetaut, E. Paleobiology of Giant Flightless Birds; ISTE Press—Elsevier: London, UK, 2018. [Google Scholar]
  39. Worthy, T.; Scofield, R. Twenty-first century advances in knowledge of the biology of moa (Aves: Dinornithiformes): A new morphological analysis and moa diagnoses revised. N. Z. J. Zool. 2012, 39, 87–153. [Google Scholar] [CrossRef]
  40. Yu, G.C.; Lam, T.T.Y.; Zhu, H.C.; Guan, Y. Two Methods for Mapping and Visualizing Associated Data on Phylogeny Using Ggtree. Mol. Biol. Evol. 2018, 35, 3041–3043. [Google Scholar] [CrossRef]
  41. Yu, G.C.; Smith, D.K.; Zhu, H.C.; Guan, Y.; Lam, T.T.Y. GGTREE: An R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  42. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  43. Oskam, C.L.; Haile, J.; McLay, E.; Rigby, P.; Allentoft, M.E.; Olsen, M.E.; Bengtsson, C.; Miller, G.H.; Schwenninger, J.-L.; Jacomb, C.; et al. Fossil avian eggshell preserves ancient DNA. Proc. R. Soc. B 2010, 277, 1991–2000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Calder, W.A. The Kiwi. Sci. Am. 1978, 239, 132–142. [Google Scholar] [CrossRef]
  45. Deeming, D.C.; Birchard, G.F. Why were extinct gigantic birds so small? J. Avian Biol. 2008, 1, 187–194. [Google Scholar] [CrossRef]
Figure 1. Extant distribution of living ratites and the Holocene distribution of moa and elephant birds. Each color relates the ratite groups to their area of origin. Kiwi birds, which range from 10–100 kg, are printed in gray color. Silhouettes of birds and eggs indicate a relative proportion of actual size. RSSD or SSD is implied by marking sex symbols on two different silhouettes within a species. For species without significant RSSD or SSD, there is only one silhouette. See Supplementary Materials Dataset Table S1 for details on the source of bird images.
Figure 1. Extant distribution of living ratites and the Holocene distribution of moa and elephant birds. Each color relates the ratite groups to their area of origin. Kiwi birds, which range from 10–100 kg, are printed in gray color. Silhouettes of birds and eggs indicate a relative proportion of actual size. RSSD or SSD is implied by marking sex symbols on two different silhouettes within a species. For species without significant RSSD or SSD, there is only one silhouette. See Supplementary Materials Dataset Table S1 for details on the source of bird images.
Biology 10 00738 g001
Figure 2. Eggshell models in finite element package ANSYS (Mechanical APDL 15.0). (a) a representative mesh of the upper-half eggshell with mesh-refinement near the top. The bottom ring (equator) is fixed. (b) a cross-sectional view of the displacement field in the y-direction (upward) subjected to a point load, F, on top. Negative value, marked in blue, represents downward displacement. δ is the displacement at the load application point. This particular case demonstrates the scenario of critical thickness.
Figure 2. Eggshell models in finite element package ANSYS (Mechanical APDL 15.0). (a) a representative mesh of the upper-half eggshell with mesh-refinement near the top. The bottom ring (equator) is fixed. (b) a cross-sectional view of the displacement field in the y-direction (upward) subjected to a point load, F, on top. Negative value, marked in blue, represents downward displacement. δ is the displacement at the load application point. This particular case demonstrates the scenario of critical thickness.
Biology 10 00738 g002
Figure 3. Allometries between body mass (BM) and the dimensionless number C (blue) and critical C (red) on a double logarithmic scale. Dash lines: 95% confidence interval. (a) two fitting lines of data from Juang et al. [3] with data points of 15 ratite species (excluding kiwi, see discussion). (b) new fitting lines combined with data points of 15 ratite species. (c) the new fitting line of C with three data points of kiwi. (d) the new fitting line of C with six data points of PGOM. Williams’ egg is marked as a light blue star, whereas the Spooner egg is marked as a red star. Purple hollow stars indicate our additional specimens. Points of larger body mass represent Genyornis newtoni, whereas smaller points represent two extinct megapodes, respectively.
Figure 3. Allometries between body mass (BM) and the dimensionless number C (blue) and critical C (red) on a double logarithmic scale. Dash lines: 95% confidence interval. (a) two fitting lines of data from Juang et al. [3] with data points of 15 ratite species (excluding kiwi, see discussion). (b) new fitting lines combined with data points of 15 ratite species. (c) the new fitting line of C with three data points of kiwi. (d) the new fitting line of C with six data points of PGOM. Williams’ egg is marked as a light blue star, whereas the Spooner egg is marked as a red star. Purple hollow stars indicate our additional specimens. Points of larger body mass represent Genyornis newtoni, whereas smaller points represent two extinct megapodes, respectively.
Biology 10 00738 g003
Figure 4. Factor of safety F.S., C number and the phylogenetic relationship of 18 ratite species. The phylogenetic tree was plotted with cladogram using ggtree [40,41] and the nearby graph was plotted by package ggplot2 [42]. (a) phylogeny of ratites is based on a tree of Mitchell et al. [7], and phylogeny of moa is supplemented by Bunce et al. [30]. (b) the pink bar indicates the case of incubation by females, whereas blue indicates incubation by males. The purple bar in some species indicates F.S. overlap of two sexes, whereas green indicates species with no sexual size dimorphism. (c) the C numbers for ratites were predicted by FEA simulations. See Supplementary Dataset for details on the source of bird images (Table S1).
Figure 4. Factor of safety F.S., C number and the phylogenetic relationship of 18 ratite species. The phylogenetic tree was plotted with cladogram using ggtree [40,41] and the nearby graph was plotted by package ggplot2 [42]. (a) phylogeny of ratites is based on a tree of Mitchell et al. [7], and phylogeny of moa is supplemented by Bunce et al. [30]. (b) the pink bar indicates the case of incubation by females, whereas blue indicates incubation by males. The purple bar in some species indicates F.S. overlap of two sexes, whereas green indicates species with no sexual size dimorphism. (c) the C numbers for ratites were predicted by FEA simulations. See Supplementary Dataset for details on the source of bird images (Table S1).
Biology 10 00738 g004
Figure 5. The regression line of (a) initial egg mass and (b) eggshell thickness on body mass. Data used to plot regression are from Juang et al. [3]. Green dots indicate species with no RSSD. Red dots (female) and blue dots (male) connected by dashed lines indicate species with RSSD. The numbers 1 to 18 represent different species, as denoted in Table 1.
Figure 5. The regression line of (a) initial egg mass and (b) eggshell thickness on body mass. Data used to plot regression are from Juang et al. [3]. Green dots indicate species with no RSSD. Red dots (female) and blue dots (male) connected by dashed lines indicate species with RSSD. The numbers 1 to 18 represent different species, as denoted in Table 1.
Biology 10 00738 g005
Figure 6. Factor of safety F.S., C number and the phylogenetic relationship of megapodes and Genyornis newtoni. (a) phylogenetic tree of megapodes is based on Shute et al. [34]. (b) F.S. of extinct species were calculated with Williams’ egg (red), Spooner egg (blue), and additional specimen (yellow), respectively. Extant megapodes are illustrated as black dots. (c) different colors of C number indicate identical implication of F.S.
Figure 6. Factor of safety F.S., C number and the phylogenetic relationship of megapodes and Genyornis newtoni. (a) phylogenetic tree of megapodes is based on Shute et al. [34]. (b) F.S. of extinct species were calculated with Williams’ egg (red), Spooner egg (blue), and additional specimen (yellow), respectively. Extant megapodes are illustrated as black dots. (c) different colors of C number indicate identical implication of F.S.
Biology 10 00738 g006
Table 1. Absent of MBM indicates that species with no SSD or RSSD. FBM: female body mass (kg), MBM: male body mass (kg), and ESR: egg size ratio (egg mass divided by average body mass). a cassowary, b emu, c kiwi bird, d ostrich, e rhea, f elephant bird, g moa.
Table 1. Absent of MBM indicates that species with no SSD or RSSD. FBM: female body mass (kg), MBM: male body mass (kg), and ESR: egg size ratio (egg mass divided by average body mass). a cassowary, b emu, c kiwi bird, d ostrich, e rhea, f elephant bird, g moa.
SpeciesEgg Length, B (mm)Egg Width, A (mm)Shell Thickness, t (mm)Egg Mass, W (g)Body Mass, M (kg)ESRRef.
FBM
(Max.)
FBM
(Min.)
MBM
(Max.)
MBM
(Min.)
1 Struthio camelus d1581311.9216008563130801.79%[3,12]
2 Rhea americana e128860.95254020 1.75%[3,28]
3 Rhea pennata e126920.865842515 2.92%[22,28]
4 Dromaius novaehollandiae b136890.94610453530201.88%[3,12]
5 Casuarius casuarius a135920.97644503530251.84%[3,12]
6 Casuarius unappendiculatus a136900.9760464.3552.6540.733.31.26%[22,27]
7 Apteryx owenii c110700.42951.911.30.923.17%[22,28]
8 Apteryx haastii c123770.534003.31.52.61.218.59%[22,28]
9 Apteryx mantelli c125780.54173.272.092.591.8217.06%[22,27]
10 Aepyornis maximus f3032243.79120541334 2.08%[6,25,26]
11 Megalapteryx didinus g1601081.1710238028 1.89%[23,30]
12 Dinornis novaezealandiae g1901501.37523431609169342.65%[10,23,24]
13 Dinornis robustus g2401781.234167275102113613.03%[10,23,24]
14 Pachyornis elephantopus g2211501.79272510649 3.52%[23,29]
15 Anomalopteryx didiformis g1651201.3913026426 2.89%[23,30]
16 Euryapteryx curtus curtus g121970.9624302020152.94%[12,23]
17 Euryapteryx curtus gravis g2051581.528041056780553.65%[12,23]
18 Emeus crassus g17913411761805250363.23%[23,24,31,32]
Note: Data used in the simulation.
Table 2. Comparison of two known eggs (with measured parameters) with possible source taxa (with estimated body mass) to generate ESR. These values were used in the FEA simulations. ESR: egg size ratio (W/M).
Table 2. Comparison of two known eggs (with measured parameters) with possible source taxa (with estimated body mass) to generate ESR. These values were used in the FEA simulations. ESR: egg size ratio (W/M).
Egg SpecimenSpeciesEgg Length, B (mm)Egg Breath, A (mm)Shell Thickness, t (mm)Egg Mass, W (g)Body Mass, M (kg)ESRRef.
WilliamsGenyornis newtoni1551251.1513271920.69%[14,33]
Latagallina naracoortensis1551251.1513276.121.75%[13,14,34]
Progura gallinacea1551251.1513277.717.23%[13,14,34]
Spooner EggGenyornis newtoni126971.36501920.34%[13,33]
Latagallina naracoortensis126971.36506.110.66%[13,34]
Progura gallinacea126971.36507.78.44%[13,34]
Table 3. Similar eggshell geometries with PGOM (putative Genyornis oological material). Data were published by Juang et al. [3]. The bold text indicates similarity to the Spooner egg <15%. (Length: 126 mm, width: 97 mm).
Table 3. Similar eggshell geometries with PGOM (putative Genyornis oological material). Data were published by Juang et al. [3]. The bold text indicates similarity to the Spooner egg <15%. (Length: 126 mm, width: 97 mm).
SpeciesCommon NameEgg Length, B (mm)Egg Width, A (mm)Shell Thickness, t (mm)Egg Mass, W (g)Body Mass, M (g)Dimensionless Number, CFactor of Safety, F.S.
Diomedea ExulansWandering Albatross129.579.70.58455.0819031852.16
Phoebastria NigripesBlack-footed Albatross108.269.20.50286.0319563284.81
Phoebetria PalpebrataLight-mantled Albatross104.064.50.48243.0315054283.95
Macronectes GiganteusSouthern Giant Petrel104.465.90.58237.0439510,8655.61
Aptenodytes PatagonicusKing Penguin104.575.80.80306.011,75113,7414.00
Cygnus ColumbianusTundra Swan106.968.20.76280.0675012,6552.41
Pinguinus ImpennisGreat auk124.075.80.74372.0500082997.56
Rhea AmericanaGreater Rhea128.086.00.90525.023,00096372.80
Casuarius CasuariusSouthern Cassowary135.092.10.97644.044,00091082.12
Dromaius NovaehollandiaeEmu136.089.00.94610.034,20080432.02
Apteryx AustralisSouthern Brown Kiwi125.878.50.50434.0233039242.79
Macrocephalon MaleoMaleo105.661.70.38222.0156434883.20
Gymnogyps CalifornianusCalifornia Condor110.266.70.92280.0845017,9116.08
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yen, A.; Wu, H.-J.; Chen, P.-Y.; Yu, H.-T.; Juang, J.-Y. Egg Incubation Mechanics of Giant Birds. Biology 2021, 10, 738. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10080738

AMA Style

Yen A, Wu H-J, Chen P-Y, Yu H-T, Juang J-Y. Egg Incubation Mechanics of Giant Birds. Biology. 2021; 10(8):738. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10080738

Chicago/Turabian Style

Yen, An, Hsiao-Jou Wu, Pin-Yi Chen, Hon-Tsen Yu, and Jia-Yang Juang. 2021. "Egg Incubation Mechanics of Giant Birds" Biology 10, no. 8: 738. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10080738

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop