Royal Society Publishing

Metabolic ‘engines’ of flight drive genome size reduction in birds

Natalie A. Wright , T. Ryan Gregory , Christopher C. Witt

Abstract

The tendency for flying organisms to possess small genomes has been interpreted as evidence of natural selection acting on the physical size of the genome. Nonetheless, the flight–genome link and its mechanistic basis have yet to be well established by comparative studies within a volant clade. Is there a particular functional aspect of flight such as brisk metabolism, lift production or maneuverability that impinges on the physical genome? We measured genome sizes, wing dimensions and heart, flight muscle and body masses from a phylogenetically diverse set of bird species. In phylogenetically controlled analyses, we found that genome size was negatively correlated with relative flight muscle size and heart index (i.e. ratio of heart to body mass), but positively correlated with body mass and wing loading. The proportional masses of the flight muscles and heart were the most important parameters explaining variation in genome size in multivariate models. Hence, the metabolic intensity of powered flight appears to have driven genome size reduction in birds.

1. Introduction

(a) Genome size evolution

The causes and consequences of variability in nuclear genome size have been a subject of research and discussion for many decades [1,2], but much remains to be learned. Several intriguing comparative patterns have been identified that provide insights into the evolutionary forces that have shaped genome size diversity. As a notable example, a proposed link between powered flight and reduced genome sizes in bats, birds and pterosaurs has been interpreted as evidence that the metabolic demands of flight exerted selective pressures for small cells with reduced DNA content [35]. These three historical instances of flight origins and genome constriction are suggestive, but comprise meagre statistical evidence for a mechanistic link. Furthermore, ancestral state estimates for the archosaur phylogeny revealed that much of the genome size reduction in the ancestors of modern birds predated the origin of flight [6].

If flight ability imposes strong evolutionary constraints on genome size, then interspecific diversity in flight ability should explain at least some of the variability in genome size within volant clades. Indeed, evidence to this effect has begun to accumulate, with some caveats. Hughes & Hughes [3] found that flightless birds have larger genomes than volant relatives, but this pattern became non-significant after accounting for phylogenetic inertia [4]. Hovering hummingbirds exhibit the highest mass-specific metabolic rates among vertebrates [7], and a large sample of hummingbird species was found to have uniformly small genomes and to include the smallest amniote genome measured to date [8]. In 74 species of temperate, migratory passerine birds, genome size was positively associated with wing-loading index, indicating that species that have evolved larger wings for flight efficiency have also evolved smaller genomes [9]. Additional components of flight performance beyond wing size and shape have yet to be investigated for potential effects on genome size.

The metabolic rate hypothesis is the leading explanation for reduced genomes in volant organisms [3]. According to this hypothesis, the size of the genome imposes a minimum size constraint on nucleated cells, and the higher proportional surface areas of smaller cells make them conducive to higher metabolic flux [10]. The energy required to produce lift and thrust requires sustained high metabolic output [11], resulting in indirect selection for smaller genomes to reduce the size of nucleated cells [10]. In keeping with this hypothesis, genome size and red blood cell size are negatively correlated with resting metabolic rate across vertebrates and within vertebrate classes, including birds [12]. Aside from the putative metabolic rate effect, a reduction in cellular DNA content and a concomitant decrease in cell size might allow flying organisms to achieve a reduction in body mass [12] and may even enhance efficiency of neural function associated with maneuverability [12,13]. Under the latter mechanisms, wing loading and wing shape would be predicted to correlate strongly with genome size.

(b) The avian flight ‘engines’

The pectoral and cardiac muscles are the metabolic ‘engines’ of avian flight, and therefore might illuminate the causes of genome size evolution. The avian pectoral flight muscles, comprising the pectoralis major and supracoracoideus, provide the power for downstroke and upstroke, respectively. The pectoral flight muscles vary greatly in size across birds, ranging from approximately 5% of the body mass in rails and ground cuckoos to over 30% in some hummingbirds and pigeons [14,15]. The ratio of heart mass to body mass, the heart index [16], indicates the relative allocation of resources for oxidative metabolism. Both maximum cardiac output and aerobic power input scale with heart mass to a power of approximately 0.88 [17,18], and heart mass accurately predicts maximum metabolic rate and aerobic scope [19]. Although the flight muscles and heart generate and supply power for a substantial portion of locomotory and thermogenic activity [2023], species variation in relative flight muscle size and heart index indicates the degree of importance of burst power [24] and endurance flight performance [17], respectively. Previous studies using non-phylogenetic methods have found that genome size was negatively correlated with heart index in birds (n = 53), mammals, non-avian reptiles and amphibians [16].

(c) This study

The link between flight and reduced genome size has received substantial attention [35,12,25], but there is not yet a broad comparative analysis of variation in the flight and metabolic structures as they relate to genome size within a volant clade. In this study, we developed the largest internally consistent dataset on bird genome sizes, comprising 18 orders, 76 families and 422 species. Species were selected to represent a diversity of environments, ecological niches and life histories. We sought to elucidate the effects of flight on genome size evolution during the diversification of birds using measurements of three flight-related structures: the heart, the flight muscles and the wings.

2. Material and methods

(a) Data collection

We captured wild birds and drew blood from the brachial vein to prepare blood smear slides. Birds were euthanized, prepared as museum specimens with frozen tissues and deposited at the Museum of Southwestern Biology (University of New Mexico, USA) and/or Centro de Ornithología y Biodiversidad (Lima, Peru). Slides of blood blotted from the liver were prepared for individuals from which we were unable to draw blood. One investigator (N.A.W.) estimated genome size by Feulgen image analysis densitometry using the protocol of Hardie et al. [26], with chicken, Gallus gallus (1.25 pg), as the standard, and more than 200 nuclei measured per individual. For species represented by more than one individual, mean within-species standard deviation of genome size estimates was 0.035, whereas the among-species standard deviation in genome size was 0.138. The low intraspecific variation, relative to interspecific variation, suggests that estimates of genome size using only one individual should be sufficiently representative of the species. Genome size estimates are available in the electronic supplementary material, table S1.

We weighed each bird and extracted and weighed whole hearts (after blotting to remove blood) and the pectoralis major and supracoracoideus flight muscles. Relative flight muscle size was calculated by dividing total flight muscle mass by body mass. Heart index was calculated by dividing the heart mass by the body mass [16]. We used proxies for aspect ratio and wing loading that could be measured from museum study skins. The hand-wing index, or Kipp's index [27,28], is equal to the distance between the tip of the longest primary and the tip of the first secondary divided by the wing chord [27]; larger values indicate more pointed wings. This can be calculated as (WL − SL)/WL, where WL is the wing chord and SL is the distance from the wrist joint to the tip of the first secondary [27]. This index is correlated with dispersal ability in interspecific comparisons [27,29,30]. WL comprises the lengths of the manus and the longest primary feather, both of which are strongly correlated with overall wing length and scale isometrically with other wing components such as the ulna [3133]. Bird wings are generally shaped as one-quarter of an oval. Thus, wing area was estimated as the area of an oval divided by four (WL × SL × π/4). Wing area estimates were used to approximate wing loading (body mass/total wing area).

Species averages for the above morphological measurements were obtained from specimens at the Museum of Southwestern Biology and Florida Museum of Natural History. We were unable to obtain all morphological measurements for all species. Thus, with genome size estimates for 422 species, we analysed two overlapping datasets of species-average values: (i) relative flight muscle size, heart index and body mass for 289 species; and (ii) relative flight muscle size, heart index, body mass, hand-wing index and wing loading for 193 species.

(b) Data analysis

We used a time-calibrated phylogeny based on published DNA sequences for over 6000 species of birds [34] that we trimmed to include only the species in our genome size dataset. We placed the approximately 5% of the species in our dataset that were not represented in this tree following the methods of Sibly et al. [35]. We used this tree to calculate phylogenetic independent contrasts, test for phylogenetic signal in the data and perform phylogenetic generalized least-squares (PGLS) models under a Brownian motion model in R v. 2.15.1 [36] using packages ape [37], caper [38] and nlme [39]. We log-transformed body mass prior to analyses. We used Akaike's Information Criterion (AIC) to select the best model(s) for non-phylogenetic linear regression and PGLS. Multi-parameter models were discarded from consideration if a nested model, containing a subset of the same parameters, had a better AIC score [40,41].

We tested for constraints on maximum genome size that might be associated with high flight ability (i.e. high values of relative heart mass, relative flight muscle size and hand-wing index or low values of body mass and wing loading). For each measure of flight ability, we binned values into 15 equally spaced bins, following Derryberry et al. [42]. We selected the species with the maximum genome size value in each bin to include in the upper-bound set. We then tested whether these upper-bound values of genome size decline with increasing flight ability more dramatically than do all other values. To do this, we used ANCOVA to test for the significance of an interaction term between the measure of flight ability and whether a sample was included in the upper-bound set.

3. Results

There was significant phylogenetic inertia in genome size across the entire 422-species dataset. Pagel's lambda was 0.88 (95% CI: 0.788, 0.943), indicating that phylogenetic non-independence should be accounted for in comparative analyses. Blomberg's K, which estimates the amount of phylogenetic signal in the data relative to the amount expected for a trait evolving along the same tree by Brownian motion [43,44], was 0.451 (significantly greater than zero, p = 0.001).

Body mass was negatively correlated with heart index (PGLS: d.f. = 286, p < 0.001, R2 = 0.043), but not with relative flight muscle size (PGLS: d.f. = 286, p = 0.13, R2 < 0.01). Heart mass scaled with body mass to a power of 0.78 in our dataset (95% CI: 0.75, 0.81).

In single variable phylogenetic models, genome size was negatively correlated with heart index and relative flight muscle size and positively correlated with body mass and wing loading (figures 1 and 2; the electronic supplementary material, table S2). In the larger dataset (289 species), the top ranking PGLS models (i.e. those within 95% cumulative AIC weight after excluding models with uninformative parameters) included heart index and relative flight muscle size as significant predictors of genome size (table 1). In the reduced dataset (193 species), the top ranking PGLS models included relative flight muscle size, wing loading, heart index and body mass as important predictors of genome size; wing loading was unique among these variables in that it was not statistically significant (p < 0.05) in any of the top models (tables 2 and 3).

View this table:
Table 1.

Models predicting genome size within 95% cumulative AIC weight after excluding models with uninformative parameters for the 289-species dataset (f, relative flight muscle size; h, heart index; m, body mass). PGLS, phylogenetic generalized linear model; regression, conventional (non-phylogenetically corrected) linear regression.

View this table:
Table 2.

Models predicting genome size within 95% cumulative AIC weight after excluding models with uninformative parameters for the reduced dataset (N = 193 species; f, relative flight muscle size; h, heart index; wl, wing loading; hw, hand-wing index; m, mass). PGLS, phylogenetic generalized linear model; regression, conventional (non-phylogenetically corrected) linear regression.

View this table:
Table 3.

Cumulative model weights for each predictor variable in the five-parameter (N = 193 species) and three-parameter (N = 289 species) datasets, calculated by summing Akaike model weights for all models that included the variable of interest, following Arnold [40]. PGLS, phylogenetic generalized linear model; regression, conventional (non-phylogenetically corrected) linear regression.

Figure 1.

Phylogenetic independent contrasts (PICs) of species mean values for genome size plotted as a function of heart index, relative flight muscle size, body mass, hand-wing index and wing loading.

Figure 2.

Species mean values for genome size plotted as a function of heart index, relative flight muscle size, wing loading, hand-wing index and body mass, colour-coded by order. Solid lines are all-points regression lines; dashed lines are upper bound regression lines.

We found evidence of constraints on genome size with respect to some measures of flight ability. The slopes of regressions of genome size by heart index and hand-wing index, respectively, were more steeply negative for upper-bound points than for other points (table 4 and figure 2), indicating that large genome sizes were particularly unlikely to occur in species with large hearts or pointed wings.

View this table:
Table 4.

Comparison of upper bound regressions to standard linear regressions for genome size by each flight ability predictor variable. We used an ANCOVA with an interaction term between the predictor variable and whether a sample was included in the upper-bound set.

Visualization of the evolution of genome size across the avian tree revealed lability among families and orders (figure 3). Orders differed significantly from one another in average genome size (ANOVA: p < 0.001, d.f. = 404; figure 3; the electronic supplementary material, table S3), with Apodiformes exhibiting the smallest average genome sizes. Of the orders represented by at least five species, Piciformes had the largest and most variable genomes (figures 2 and 3; the electronic supplementary material, table S3). Passerine genomes were significantly smaller than those of non-passerine orders (t-test: p = 0.006; d.f. = 180; 95% CI of mean difference: (−0.08, −0.013)). Deep phylogenetic structure in genome size was evident among major clades of passerines. Nine-primaried oscines had larger genomes, on average, than other passerines (t-test: p < 0.001; d.f. = 128; 95% CI of mean difference: (0.019, 0.071)). However, extreme values were not necessarily consistent with these broader patterns. For example, out of all passerine families, the smallest average genome sizes were found in two oscine families, the vangas (Dicruridae; C-value 1.13 pg) and the indigobirds and whydahs (Viduidae; C-value 1.14 pg); and the largest passerine genome was not found in a nine-primaried oscine, but rather the phainopepla (Phainopepla nitens; C-value 1.83 pg).

Figure 3.

Family-level phylogeny with terminal branches colour-coded by the average genome size for each family. Internal branches are colour-coded by the estimated ancestral genome size, using maximum-likelihood ancestral character state estimation in the R package ape [37]. The tree is from Jetz et al. [45] and was downloaded from www.birdtree.org with data preference set to ‘Hackett all species’.

4. Discussion

(a) Which aspects of flight are related to genome size?

Flight muscle size, heart size, body size, wing aspect ratio and wing loading reveal different aspects of flight ability in birds. Each of these parameters, with the exception of wing aspect ratio, is individually correlated with genome size across a wide swath of the avian tree (figures 1 and 2; the electronic supplementary material, table S2), strongly supporting the idea that small genomes and overall flight ability are linked. In multivariate models, flight muscle size, heart size, body size and wing loading explain unique fractions of genome size variation, with flight muscle size and heart size being the most consistently important predictors of genome size. There was no indication that wing aspect ratio (measured by hand-wing index) explains any additional variation. These findings are consistent with the hypothesis that the cellular metabolic rate during flight, rather than constraints on aerial maneuverability or lift, was the main cause of genome size reduction in birds.

Relative flight muscle size was the single best predictor of genome size: it was significant in almost all top models in both the 289- and 193-species datasets (tables 13). Flight muscles provide the power for flapping flight, and as such are most important in strong take-offs and bursts of speed in flight, where power requirements are highest [24]. Thus, large flight muscles are expected to indicate an evolutionary response to selection for rapid take-off ability or powered acceleration during flapping flight. Additionally, the size of the flight muscles is a primary determinant of exercise-induced maximum metabolic rate or thermogenic capacity in intraspecific comparisons for at least five species of birds that have been tested, representing three orders and a wide range of body sizes [2023,46].

Heart index was an important and statistically significant predictor of genome size in top PGLS models (tables 13). The size of the heart constrains its stroke volume and, as a result, limits maximum cardiac output, aerobic power, exercise-induced maximum metabolic rate and aerobic scope [1719]. Routine powered flight incurs a 10- to 20-fold sustained increase in metabolic rate [11]. Therefore, heart index in birds should closely reflect aerobic power requirements for sustained flight.

Body mass generally shows strong phylogenetic signal [43] and correspondingly was present and significant in all top ranking non-phylogenetic models, but not all PGLS ones (tables 13). Either body mass or heart index was significant (or nearly so) in every top ranking model, but in no model were both parameters significant (tables 1 and 2). Heart mass scales with body mass to a power of 0.78 in our dataset, very close to the three-quarter power scaling of metabolic rate [47]. Heart size may be a more precise predictor of metabolic rate than body mass because it is more directly, mechanistically connected to metabolism than is body mass, in part because it scales isometrically with cardiac stroke volume [18,48]. Indeed, heart mass has been shown to explain variation in exercise-induced maximum metabolic rate after accounting for the effects of body mass [19]. Consistent with this, for our largest dataset, AIC decisively indicated that heart index was the best single predictor of genome size and that body mass did not explain residual variation after heart index was included in the model (table 1).

Hand-wing index did not appear in top ranking PGLS models (tables 2 and 3). The link between genome size and hand-wing index in conventional regression (figure 2 and table 3; the electronic supplementary material, table S2) appears to have been driven by phylogenetic inertia, specifically the confluence of extreme values in the hummingbirds and swifts (Apodiformes). Accordingly, the effect of hand-wing index disappears when covariance due to phylogeny is taken into account. Long, pointed wings improve high-speed flight performance by reducing drag, especially at higher speeds [49]. High aspect ratio is typical of bird species that are long-distance migrants or aerial foragers [50,51]. However, some birds with high aspect ratios engage in flight styles with low energetic demand (e.g. dynamic soarers) and accordingly would not achieve the metabolic intensity associated with flapping flight.

Wing loading was included as an important predictor variable in several top PGLS models, but was never a significant variable in these models (p > 0.25; tables 2 and 3). Birds with low wing loading are able to produce lift with low metabolic energy input [52]. The tendency for aerial specialists to have low wing loading is thought to be the basis for the previous finding that passerines with low wing loading tend to have small genomes [13]. However, wing loading reflects a balance of ecological and biomechanical pressures and many strong, energetic flyers have high wing loading (e.g. pigeons, doves and many ducks [52]). Thus, wing loading, like aspect ratio, is at least partly decoupled from flight metabolism and is unlikely to explain genome size variation if the maximum sustained rate of cellular metabolism is the key driver of genome size reduction.

(b) The importance of accounting for phylogenetic inertia

Genome size exhibited phylogenetic signal across the avian tree, as indicated by Pagel's lambda and Blomberg's K-values that were significantly greater than zero. Related taxa were less similar in genome size than expected under a model of evolution by Brownian motion, as indicated by the Blomberg's K-value being less than one [43,44]. This apparent lability is more typical of behavioural traits than constrained life-history traits such as body mass [43] and it suggests that divergent selection has acted on genome size, although measurement error is also expected to contribute to this effect. Nonetheless, as a result of the pervasive phylogenetic inertia in genome size, the PGLS analyses are expected to be more informative regarding potential mechanisms in cases where they differ from non-phylogenetically controlled analyses [53]. Hence, in the 289-species analysis, the inclusion of body mass and exclusion of heart index in the top ranking non-phylogenetic model are likely examples of type I and type II errors, respectively (table 1).

(c) Constraints on genome size associated with extreme flight ability

Genome size appears to be constrained with respect to heart index and hand-wing index, as indicated by the different slope of the regression for upper bound points (figure 2 and table 4). Species with small hearts or rounded wings may have either large or small genomes, but strong flyers with large hearts or pointed wings appear to be constrained to have small genomes. This pattern suggests that genome size evolves neutrally in poor flyers, but that large genomes impose a fitness cost on birds that experience high metabolic intensity during flight.

(d) Genome size variation across the tree of birds

The hummingbirds have exceptionally small genomes with low variation in genome size [8]. This study is the first to show that the nearest relatives of hummingbirds, the swifts (Apodidae), also have small genomes (figure 3 and the electronic supplementary material, table S3). Like many other traits of hummingbirds that are thought to be adaptive for metabolically intensive flight (e.g. high aspect ratio, reduced hindlimbs, large pectoral flight muscles), reduced genome size probably evolved in the common ancestor of hummingbirds and swifts (figure 3).

The Pici clade, comprising the woodpeckers, honeyguides, barbets and toucans, exhibits exceptionally large genomes, whereas members of its sister taxon, the Galbulae (Galbulidae and Bucconidae), have fairly small genomes (figure 3 and the electronic supplementary material, table S3). Members of the clade including Piciformes and Coraciiformes are anomalously variable, with sister lineages sometimes differing substantially in genome size (figure 3 and the electronic supplementary material, tables S1 and S2). For example, the two honeyguide species included in this study, both in the genus Indicator, exhibit strikingly different genome sizes (1.38 versus 2.05 pg; the electronic supplementary material, table S1). Han et al. [54] found large numbers of transposable elements in the Piciformes–Coraciiformes clade, which may contribute to the large average genome sizes and the high degree of variation observed in the group [4,55]. The macromutations that may have caused fluctuations in genome size among Piciformes deserve investigation by comparison of genome sequences. To this end, frozen tissues are archived for all of our voucher specimens (see electronic supplementary material, table S1).

5. Conclusion

The long-recognized link between flight and small genome size has been decisively confirmed using interspecific comparisons within the most diverse volant group. Four flight-related characteristics show that the evolution of enhanced flight ability tends to be accompanied by genome size reduction. In comparisons among phylogenetically corrected models, the best predictors of genome size proved to be the sizes of the flight muscles and heart, rather than the size or shape of the wings. The avian flight muscles generate cellular metabolic energy in proportion to their enormous bulk (5–30% of total body mass), providing the majority of power output for aerial locomotion and thermogenesis. Cardiac output constrains the power input for cellular metabolism and is in turn limited by heart size. Consistent with the metabolic rate hypothesis, these two components of the metabolic flight ‘engine’ implicate the rate of energy use as a key driver of repeated evolutionary reductions in avian genome size.

Data accessibility

Genome size estimates are available in Table S1 and on the Animal Genome Size Database (http://www.genomesize.com/). Table S1 includes museum accession numbers for each specimen used in this study. Specimen data are archived in museum databases (http://arctos.database.museum/SpecimenSearch.cfm).

Funding statement

We acknowledge NSF DEB-1146491, the Center for Evolutionary and Theoretical Immunology and R. W. Dickerman for funding support. N.A.W. is supported by the Program in Interdisciplinary Biomedical and Biological Sciences through the University of New Mexico award number T32EB009414 from the National Institute of Biomedical Imaging and Bioengineering. T.R.G. is supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. The content is the sole responsibility of the authors and does not necessarily represent the official views of the National Institute of Biomedical Imaging and Bioengineering or the National Institutes of Health.

Acknowledgements

We thank N. Jeffery for assistance in the laboratory, C. E. Gunning for help with R code and J. Hammond for statistics advice. This study was made possible by data, specimens and curatorial staff from the Museum of Southwestern Biology, the Field Museum of Natural History, and the Florida Museum of Natural History. We thank the individuals who collected key samples: J. M. Bates, E. Bautista, E. J. Beckman, P. M. Benham, S. V. Brant, D. DeSwardt, R. W. Dickerman, S. G. DuBay, S. C. Galen, E. Gendron, S. J. Hackett, A. B. Johnson, M. R. Jones, B. D. Marks, R. Nuttall, C. J. Schmitt, D. C. Schmitt, C. G. Schmitt, A. Smiley, E. W. Valdez, T. Valquí, W. Vargas, J. D. Weckstein and D. E. Willard. B. O. Wolf, J. H. Brown and two anonymous reviewers provided valuable feedback on the manuscript. This work was conducted under UNM IACUC protocol number 11-100742-mcc. We thank the management authorities of Peru (135-2009-AG-DGFFS-DGEFFS, 0377-2010-AG-DGFFS-DGEFFS, 0199-2012-AG-DGFFS-DGEFFS) and Free State, South Africa (01/8233, 01/12248).

  • Received October 24, 2013.
  • Accepted January 6, 2014.

References

View Abstract