Scaling of Avian Primary Feather Length

The evolution of the avian wing has long fascinated biologists, yet almost no work includes the length of primary feathers in consideration of overall wing length variation. Here we show that the length of the longest primary feather () contributing to overall wing length scales with negative allometry against total arm (ta = humerus+ulna+manus). The scaling exponent varied slightly, although not significantly so, depending on whether a species level analysis was used or phylogeny was controlled for using independent contrasts: . The scaling exponent was not significantly different from that predicted (0.86) by earlier work. It appears that there is a general trend for the primary feathers of birds to contribute proportionally less, and ta proportionally more, to overall wingspan as this dimension increases. Wingspan in birds is constrained close to mass (M 1/3) because of optimisation for lift production, which limits opportunities for exterior morphological change. Within the wing, variations in underlying bone and feather lengths nevertheless may, in altering the joint positions, permit a range of different flight styles by facilitating variation in upstroke kinematics.


Introduction
The total length of the avian wing derives from the underlying wing bones (humerus, radius/ulna and manus) and the functional primary feathers (Fig. 1). Although scaling exponents vary slightly depending upon whether the effects of common ancestry are controlled for using independent contrasts or not (M 0.35 and M 0.39 respectively, table 1 in [1]), it is well established that wingspan (b) in birds scales with slightly positive allometry with respect to body mass (M .1/3 ) [1][2][3][4]. This positive allometry, however, appears related to size dependent variation in flight behaviour [1]. Specifically, the line of best fit is depressed at lower body masses and elevated at high body masses, because slow speed flapping flight styles seen in smaller birds are associated with short-wings, while the soaring flight styles of larger birds favour longer wings [4]. Surprisingly, and in spite of variations in flight behavior, the relative contribution of the primary feathers to overall wing length has received little attention from ornithologists.
Curiously, total-arm (ta = humerus+ulna+manus) length does not scale with unity against wing semi-span b semi and instead appears to scale with positive allometry (ta!b semi 1:08 ), indicating that larger birds have longer wings relative to their M, but also have longer ta relative to their b [1]. An explanation forwarded to explain this disproportionate increase in ta with b [1] was that primary feather length (f prim ) is relatively shorter in longer winged birds (i.e., !b v1 ). The only data available, however, suggested f prim scaled as M 0.32 , which was not significantly different from the exponent predicted for isometry (M 1/3 ) [5]. Worcester's [5] study was, however, limited in taxa (n = 13) and, because the relationship between b and ta was not determined, whether increasing M in the sample correlated with a relatively longer ta was not known. Therefore, a trend towards shorter primaries in birds with longer ta remains a possibility [1] and required further investigation.
Nudds [1] also acknowledged that if elbow angle varied with b it would influence how close to parallel the leading edge of the humerus and ulna was and hence the relationship between ta and b ( Fig. 1). Elbow angle is extremely difficult to measure in live birds, however, because bones are not visible from the wing-surface. Plucking of feathers is unethical and undesirable, and x-ray not necessarily practical, but if negative allometry was found between f prim and ta then the effect of elbow angle could be ruled out.
'Stretched' or 'flat' wing preservations are rare in museum collections and those including the humerus intact within the skin are even more so (personal observations). However, a small collection of suitably stretched wing specimens was located in the Royal British Columbia Museum, Victoria, BC, Canada (RBCM). Even these had the humerus removed from the wing, but fortunately kept separately to permit all wing-bone measurements to be recorded from a homogenous specimen. These skins allowed us to test the hypothesis that f prim scales with negative allometry against ta (i.e., f prim !ta v1 ) as proposed by Nudds [1]. More specifically it was predicted that, over the range of wing semi-spans (b semi = 0.075 to 1.622 m) used in Nudds [1], the predicted scaling exponent between f prim !ta should approximate to 0.86, because ta!b semi 1:08 and, f prim zta!b semi , so f prim !b 0:93 and therefore f prim !ta 0:93=1:08~0:86 . A predicted exponent of 0.86 assumes that size dependent variation in f prim is entirely responsible for the positive allometry seen in ta (i.e., elbow angle is constant across all wingspans).

Methods
Humerus, ulna and manus lengths were measured using Vernier calipers to the nearest mm from the 'spread wing' bird skin collection at the RBCM. 'Total-arm' is the sum of humerus, ulna and manus length [3,6]. Primary feather length (f prim ) was measured from the distal end of digit 2 of the manus to the feather tip, parallel to the feather shaft ( Fig. 1). The primary feather chosen was that contributing the most distal point of the wing representing maximum b.
Because the data set comprises interspecific measures (Table 1), the effects of common ancestry must be considered to prevent spurious correlations resulting from common descent rather than from independent evolution. Here a comparative analysis using standardized independent contrasts, conducted in CAIC version 2.6.9 [7], was used. The analyses were implemented in three ways. Initially the scaling relationships were calculated using species as independent data points. The analysis was then repeated using CAIC and the phylogenetic hypotheses of Sibley and Ahlquist [8], and finally CAIC was implemented using the phylogenetic hypotheses of Livezey and Zusi [9]. A punctuated model of evolution was used in both cases: the branch length estimates of Sibley and Ahlquist [8] are disputed and none are available for the phylogeny of Livezey and Zusi [9]. The topological disagreement between these two hypotheses [8,9] is useful, because if phylogeny is going to affect the results, then using two different phylogenies is likely to have a greater effect than changes to branch lengths within a single phylogeny. Use of two different phylogenies should therefore indicate whether the scaling relationships determined are likely to be affected by future refinements of phylogenetic topology.
The relationship between f prim and ta was investigated using the empirical scaling formula y~kx a where a is the allometric exponent (slope) and k is the allometric coefficient (intercept), which was in turn determined using a Model II reduced major axis (RMA) regression [10][11][12]. Regression analyses using independent contrasts were performed through the origin [7]. The RMA slope was calculated as the ordinary least squares (OLS) Model I slope (regression coefficient) divided by the OLS correlation coefficient, and 95% confidence limits were calculated following Sokal and Rohlf [13]. The standard error (s.e.) of the RMA slope was taken as equal to that of the s.e. of the OLS slope. Two-tailed t-tests were used to test for differences between calculated slopes and the slopes predicted for geometric similarity (a = 1) or predicted from Nudds

Results
In all three analyses the relationship between f prim and ta was statistically significant, with the scaling exponent dependent upon the analysis used (Fig. 2). The scaling exponent determined using species as independent data points was significantly below (t = 26.50, p,0.001) that predicted for geometric similarity (a = 1). Similarly, both CAIC using the phylogeny of Livezey and Zusi [9] and the phylogeny of Sibley and Ahlquist [8] produced slopes significantly below 1 (t = 22.72, p,0.05 and t = 23.12, p,0.05 respectively). In all three cases the scaling exponents were below, yet not significantly different from, that predicted (i.e., f prim !ta 0:86 ) by Nudds [1] as demonstrated by the 95% confidence intervals (Fig. 2).

Discussion
As predicted by Nudds [1] and contrary to that suggested by the data of Worcester [5], there is a general trend for the primary feathers of birds to contribute proportionally less to overall winglength as b semi increases. The sample size here was relatively small (n = 34) compared to the sample sizes (n = 306) used to investigate the scaling of ta [1], which precludes any analysis of flight style or ecologically driven variation in f prim /ta ratio. The wingspan of birds is constrained close to M 1/3 , because of optimisation for lift, limiting the opportunities for exterior morphological change. Within the wing, however, variations in underlying bone ratios may permit a range of different flight styles, by possibly facilitating variation in upstroke kinematics [6]. It is not unreasonable to expect the relationship between f prim and ta to also vary depending upon the ecology or flight style of the bird.
The scaling relationship determined here between f prim and ta does not entirely exclude the possibility of size dependent variation in elbow angle. Although there were no statistical differences between the calculated scaling exponents (Fig. 2) and the 0.86 predicted [1], they were lower (0.78-0.82) and the 95% confidence intervals broad. Of course, the angle at the elbow in a stretched out wing when a bird is having its wingspan measured [14] is not Table 1. The primary feather and total arm data set.

Number in figure 2a
Species n Total-arm (m) Primary feather (m) necessarily functional. Instead, it could just be an artefact of how the bird is held by the researcher. Indeed, the elbow angle is likely varied in flight and during a wing-stroke [15,16]. This, of course, begs the question of exactly what we measure when we measure b in a bird: it may be maximum extended wingspan, but is this used during flight? In hummingbirds, span in flight is effectively the width of the body plus the distances from the wrists to wing tips [17], but in other birds the portions of the wings between the wrists and the body need to be considered [18]. Measurements of functional wingspan from birds in flight are long overdue.
In conclusion, there is a general trend for f prim to contribute relatively less to overall wingspan in larger birds. Conversely, ta contributes more to the overall length with increasing b. Why this trend exists is not immediately obvious. Although tentative at this stage, the scaling of ta and f prim may be the product of an as yet unidentified optimum ratio for feathers to wing-skeleton length within the biomechanical and aerodynamic constraints acting upon the scaling of b (M 1/3 ) [1]. Similarly, whether the scaling is driven by aerodynamics, feather biomechanical properties or a combination of both requires further investigation. Figure 2. Scatter plots of log 10 primary feather length (m) against log 10 total-arm length (sum of humerus, ulna and manus length in m). The regression lines (dashed lines are 95% C.I.s) describing the relationship were A) species treated as independent data points: y = 0.57 x 0.78 (0.71-0.85) , t = 23.35, n = 34, r 2 = 0.94, p,0.001, B) phylogenetic independent contrasts (IC) using the phylogeny of Livezey and Zusi [9]: y = x 0.82 (0.68-0.96) , t = 12.26, n = 24, r 2 = 0.84, p,0.001 and C) phylogenetic independent contrasts (IC) using the phylogeny of Sibley and Ahlquist [8]: y = x 0.79 (0.64-0.93) , t = 11.46, n = 21, r 2 = 0.84, p,0.001. See table 1 for the species corresponding to the numbers in panel A. doi:10.1371/journal.pone.0015665.g002