- Research article
- Open Access
Beak and skull shapes of human commensal and non-commensal house sparrows Passer domesticus
BMC Evolutionary Biology volume 13, Article number: 200 (2013)
The granivorous house sparrow Passer domesticus is thought to have developed its commensal relationship with humans with the rise of agriculture in the Middle East some 10,000 years ago, and to have expanded with the spread of agriculture in Eurasia during the last few thousand years. One subspecies, P. d. bactrianus, residing in Central Asia, has apparently maintained the ancestral ecology, however. This subspecies is not associated with human settlements; it is migratory and lives in natural grass- and wetland habitats feeding on wild grass seeds. It is well documented that the agricultural revolution was associated with an increase in grain size and changes in seed structure in cultivated cereals, the preferred food source of commensal house sparrow. Accordingly, we hypothesize that correlated changes may have occurred in beak and skull morphology as adaptive responses to the change in diet. Here, we test this hypothesis by comparing the skull shapes of 101 house sparrows from Iran, belonging to five different subspecies, including the non-commensal P. d. bactrianus, using geometric morphometrics.
The various commensal house sparrow subspecies share subtle but consistent skeletal features that differ significantly from those of the non-commensal P. d. bactrianus. Although there is a marked overall size allometry in the data set, the shape difference between the ecologically differentiated sparrows cannot be explained by differences in size alone. Relative to the size allometry commensal house sparrows exhibit a skull shape consistent with accelerated development (heterochrony), resulting in a more robust facial cranium and a larger, more pointed beak.
The difference in skull shape and robustness of the beak between commensal and non-commensal house sparrows is consistent with adaptations to process the larger and rachis encapsulated seeds of domesticated cereals among human associated populations.
The avian beak and associated traits have long been popular targets for studies of adaptive trait evolution. The beak is the main food-processing tool of a bird. Hence, its size and shape are expected to be chiselled by natural selection to fit the feeding ecology of its carrier. One celebrated example of adaptive beak evolution is the radiation and associated diversification of feeding niches and beak shapes among Darwin’s finches [1, 2]. Long-term field studies of these finches have provided strong evidence that beak morphology evolves rapidly in response to changing ecological conditions, such as food type, food availability and interspecific competition [2–4]. Morphological evolution may be constrained by various factors, however . For instance, the size and shape of morphological traits in animals may be strongly constrained by body size and may exhibit more or less fixed allometric size relationships across species [6, 7]. Yet, the avian beak appears to be a trait that often exhibits high levels of evolvability [8–11].
The house sparrow Passer domesticus is a granivorous, widely distributed, human associated bird species [12, 13]. The origin of its commensal relationship with humans has been debated for some time [14–16]. However, in a recent paper Sætre et al.  used population genetic data, as well as previously published fossil and ecological data, to suggest that the house sparrow became associated with early human agricultural societies that rose in the Middle East about 10,000 years ago and then experienced a massive size and range expansion from about 4,000 years ago as agricultural civilizations spread through the Palearctic and Oriental regions. Although the house sparrow comprises several phenotypically distinct subspecies, these are virtually identical at neutral genetic markers, suggesting that phenotypic differentiation is of very recent origin .
One of the subspecies has apparently maintained the ancestral ecology of the species, however. The subspecies P. d. bactrianus is not associated with humans but breeds in natural or seminatural habitats, such as riverine scrub, hedgerows and trees near pastures and grassland, often far from human habitation [17, 18]. It is similar to other house sparrows at neutral genetic markers, however, consistent with the hypothesis that the ecological difference is associated with a recent ecological transition of the commensal house sparrows . In autumn P. d. bactrianus move southwards in large flocks from their breeding grounds in Central Asia to reach their wintering grounds in south-east Iran and western parts of the Indian sub-continent [19–21]. Significant for this study, P. d. bactrianus feed mainly on wild grass seeds whereas the commensal subspecies prefer seeds from cultivated crops, such as wheat and oats . The evolution of bony structures is possibly less constrained by weight in sedentary than in migratory birds because they are less dependent on endured flight. Hence, the abandonment of migratory behaviour in commensal house sparrows may have facilitated any evolutionary trend towards increased size and robustness of the beak and skull structures.
The abandonment of migratory behaviour in commensal house sparrows is likely to be an adaptation to the year-round supply of food provided by sedentary human societies through storage of cereals, spilling and feeding of domestic animals . Further, the transition from a diet of mainly wild seeds to one consisting mainly of grains from cultivated cereals  is likely to have imposed selection pressures on beak and skull morphology of the sparrows. Domestication of plants is associated with an increase in seed size [22–25], but domesticated cereals also differ from their wild relatives in several other traits, such as seed hardness and texture . Moreover, domesticated cereal seeds remain encapsulated in a tough rachis that holds the seed together, whereas in wild cereals the rachis fragmentizes at ripening . These trait differences would most certainly affect seed processing in a granivorous bird and hence the optimal size and shape of its beak and skull. In Darwin’s finches handling of larger and harder seeds is associated with larger, deeper beaks and a correlated increase in robustness of the muscle and jaw architecture of the birds .
Here, we apply geometric morphometrics to study variation in the beak and skull of commensal and non-commensal house sparrows from the same overall geographic region, namely Iran. The diversity of house sparrow subspecies is high in Iran (five recognized subspecies, including the non-commensal P. d. bactrianus) . Moreover, climatic and environmental conditions vary extensively regionally, from the moist and fertile temperate shore of the Caspian Sea in the north, through dry steppes and deserts in central parts, to tropical habitats in the south of the country. Climatic conditions are likely to affect agricultural practice, the choice of which crop species to cultivate and seed characteristics, all of which may affect optimal beak and skull morphology . However, cultivated crops differ from their wild ancestors in consistent ways , suggesting that adaptations to commensalism per se may have larger effects on morphological evolution than such local climatic variation. If adaptations to cultivated cereals have played a significant role in the recent morphological evolution of the house sparrow, we expect to find relatively larger divergence between commensal and non-commensal sparrows in beak and skull morphology than among the various commensal house sparrow subspecies, despite extensive environmental variation. We further predict that the commensal types should share (derived) skeletal features. Finally, we specifically predict that commensal house sparrows would have larger, more robust beaks and skull features associated with a more robust muscle and jaw architecture compared to the non-commensal P. d. bactrianus, as adaptations to process the larger and tougher seeds of cultivated cereal crops. At the proximate level, accelerated development (heterochrony) is a well-characterized ontological process that one could predict would yield the necessary changes in robustness.
101 completely ossified adult skulls from the house sparrow were collected from 16 populations in Iran (Figure 1). The specimens represent five of the eleven house sparrow subspecies described by Summers-Smith  that have been defined based on differences in plumage coloration, size and ecology; the non-commensal and migratory P. d. bactrianus (N = 12) and the four commensal subspecies P. d. biblicus (N = 33), P. d. persicus (N = 30), P. d. hyrcanus (N = 12) and P. d. indicus (N = 14). The specimens were obtained from different museum and university collections in Iran.
For geometric morphometrics on bird skulls and beaks, landmark analysis has been the most common approach (e.g. ). For this study, we instead applied outline analysis of the skull shapes. Outlines may be better adapted to complex shapes having important shape variation contained in smooth areas between any well-defined homologous landmarks. For cranial shapes, outline analysis can outperform landmark analysis in terms of statistical discrimination power . Outline analysis circumvents the subjective selection of shape features inherent to landmark analysis. Finally, outline analysis allows visualization of statistical results as full biological shapes rather than as deformation grids or truss networks of a small collection of landmarks.
Digital images (dorsal view) were prepared using an Olympus DP71i camera connected to an Olympus BX51 stereomicroscope (Olympus Corporation, Tokyo, Japan). The camera was fixed to a custom made rack that ensured that the skulls where photographed from the same distance and 90˚ angle. The skulls were fixed on a plasticine platform and care was taken to orient the skull consistently relative to the camera. Any residual random distortion is expected to produce unbiased noise. To quantify morphological variation the x- and y-coordinates of each skull outline were digitized using the tpsDig 2.2 software . The starting point of the outline was defined as the tip of the bill. 150 roughly equally spaced points were measured on each outline (the subsequent shape analysis does not require equal spacing).
Elliptic Fourier Analysis, EFA  is a well-established technique for outline analysis. The procedure involves spectral (Fourier) decomposition of the x and y increments along the contour. Each harmonic is represented by four coefficients, i.e. the cosine and sine amplitudes of Δx and Δy. This is equivalent to decomposition into a sum of harmonically related ellipses. EFA was carried out in the software Past, version 2.07 , with the “invariant to starting position and rotation” option. Size was recorded by the length of the semimajor axis of the first harmonic ellipse . After this calculation, size was removed prior to further EFA analysis, to ensure that the EFA coefficients themselves contained only shape variation, in order to reduce high-frequency noise and limit the number of variables. A plot of the average EFA spectrum (Figure 2) showed very small amplitudes beyond 15 harmonics and we thus chose this as the cut-off point for inclusion in the analyses.
The resulting 60 coefficients were subjected to further multivariate analysis in Past. For testing of multivariate group means we used a permutation test, PerMANOVA,  based on Euclidean distance. For regression of shape on size we used ordinary least squares multivariate regression with MANOVA testing. For tests on allometry in commensal versus non-commensal subspecies we included size, commensality and size-commensality interaction as independent variables.
The allometric trajectories of the different subspecies were also investigated using the Common Allometric Component (CAC) method of Mitteroecker et al. . The CAC is a vector in shape space computed by regression of shape variables (with centered species means) on log size. Additional residual shape components, orthogonal to the CAC, are calculated by PCA on the residuals after CAC regression. The resulting shape vectors were visualized by generating synthetic skull outlines in Past, version 3.0.
Size and shape allometry
The sizes (Figure 3) differed only slightly between the non-commensal P. d. bactrianus with mean size 12.13 mm and the four commensal subspecies with pooled mean size 12.38 mm (t = 2.00, p = 0.048). However, there was considerable size variation among the subspecies (one-way ANOVA, F = 9.85, p < 0.001) and particularly P. d. biblicus stood out as larger than P. d. bactrianus (Tukey’s Q = 5.97, p = 0.00063) and P. d. persicus (Q = 5.91, p = 0.00072).
A multivariate regression of EFA coefficients with size as independent variable was highly significant (Wilks’ λ = 0.22, F = 2.39, p = 0.0021), demonstrating clear size allometry for the data pooled over all subspecies. No obvious nonlinearity or heteroskedasticity was observed in the regression, and the data were therefore left untransformed. Figure 4 shows the shape predicted from the regression for semimajor axis sizes of 11 and 14 mm. The cranial shape predicted for a size of 11 mm (S11) show some morphological differences with regard to the one predicted for a size of 14 mm (S14): (a) S14 have relatively smaller braincases than S11, (b) the interorbital width should be more reduced proportionately in S11 sparrows than in S14 ones, and (c) the analysis also predicts a more significantly developed facial cranium (splachnocranium), and a longer bill in particular, in S14 compared to S11 individuals.
Shape differences between subspecies
The EFA coefficients of the five subspecies are significantly different (one-way PerMANOVA, F = 1.89, p = 0.023). To visualize the shape difference between P. d. bactrianus and the four commensal subspecies, we show the average (consensus) shapes from the two groups in Figure 5. The commensal subspecies are associated with a more robust facial skeleton and a larger beak. The axis of shape variation in Figure 5 is visually similar to the size allometric axis shown in Figure 4. Accordingly, there is positive correlation between the vector of differences between group means and the size regression slopes (R = 0.32, p < 0.001).
Multivariate, multiple regression of shape with size, commensality and interaction as independent variables was highly significant (Wilks’ λ = 0.020, F = 1.85, p < 0.001). The analysis showed that shape was significantly dependent on all of size (λ = 0.26, F = 1.94, p = 0.014), commensality (λ = 0.29, F = 1.68, p = 0.043) and size-commensality interaction (λ = 0.29, F = 1.65, p = 0.048). In other words, there is both an overall size allometry and a difference in allometric trajectories between commensal and non-commensal subspecies.
A plot of the scores on the Common Allometric Component axis and the second Residual Shape Component (RSC 2, explaining 18% of residual shape variation) is shown in Figure 6. The first RSC (RSC 1, 33%) was found to describe left-right asymmetry, and showed no structure with respect to taxonomy. The visualization of overall size allometry (CAC) confirms the result shown in Figure 4, obtained similarly by regression but with linear size and without within-group centering. RSC 2 describes shape variation mainly restricted to the orbital area, with higher scores associated with narrower interorbital area and more pronounced postorbital processes. Although there is considerable overlap between the subspecies, the non-commensal P. d. bactrianus occupies a region of particularly low scores on RSC 2. It may also be noted that P. d. bactrianus, and to some extent the similarly small-sized P. d. persicus, shows correlation between CAC and RSC 2, not apparent in the other subspecies.
Further details are apparent when plotting size against CAC (Figure 7). The different subspecies show some vertical displacement along the CAC axis, which could be referred to as heterochrony in a very loose sense, although the differences in allometric trajectories make this term less appropriate . In addition, the larger subspecies P. d. indicus and P. d. biblicus have smaller slopes in the size-CAC plot than the smaller three subspecies.
The morphological differences described statistically above and depicted in Figure 5 are easily visible when inspecting the actual bone structures of the different skulls. In Figure 8 we show the skull of a non-commensal P. d. bactrianus and a commensal house sparrow (P. d. persicus). P. d. persicus was chosen because it is the subspecies most similar in size to P. d. bactrianus. The chosen pictures were otherwise randomly selected. The skull of P. d. persicus appears overall more robustly constructed than the skull of P. d. bactrianus, and the splachnocranium and the beak are both elongated and more heavily built. The postorbital process (a) and the lachrymal (b) are more expanded laterally in P. d. persicus than in P. d. bactrianus (cf. Figure 6). Although not included in the statistical analysis above since the features do not affect the outline of the skull, it also appears that the oral surface of the maxillare, as one can see through the nasal openings, is wider (c), and that the nasal process of the rostrum (d) is larger in P. d. persicus compared to P. d. bactrianus. Differences in the occipital region (e) are also consistent with increased robustness in the commensal subspecies (Figure 6).
We demonstrate subtle but consistent morphological differences in the skull and beak of commensal and non-commensal house sparrows. Across subspecies residing in different environments commensal house sparrows have consistently more robust facial craniums and larger, more pointed beaks compared to the non-commensal P. d. bactrianus. We further found a clear, overall allometric relationship between skull size and shape, larger skulls having a relatively smaller brain case and a larger, more robust facial cranium than smaller ones. The commensal house sparrows were on average slightly larger than the non-commensal P. d. bactrianus, but the difference in cranial shape persisted after controlling statistically for size. The morphological changes thus appear consistent with accelerated development (heterochrony) of commensal house sparrows, assuming that the more “juvenile-like” skull of P. d. bactrianus represents the ancestral state as suggested by a previous study . Our statistical analyses indicate a significant interaction between size and ecology (commensal vs. non-commensal). An allometric shift like the one demonstrated here is often considered a hallmark of an adaptive morphological change (e.g. [35, 36]). Additionally however, we found an orthogonal shape variation of commensality only, suggesting a relatively complex ontology in the development of beak and skull shape in the case of these sparrows. When plotting size against Common Allometric Component (Figure 7) there appears to be a gradient from the smaller subspecies to the larger ones with decreasing slopes but increasing intercepts. This may be referred to as a heterochronic pre-displacement, but note that strictly speaking, heterochrony is defined with respect to age and not size.
The house sparrow is primarily a seed eater and vegetable matter makes up 85-90% of the food, complimented with more protein rich food of animal origin, particularly in the diet provided to the nestlings . Although commensal house sparrows can be quite opportunistic in what they eat they are regarded as specialized on cultivated cereal crops and will prefer oats and wheat when there is a choice . The larger and rachis encapsulated seeds of domesticated cereals [24, 26] are therefore likely to have changed the optimal feeding apparatus of the commensal house sparrows in the direction of increased robustness.
Many of the skeletal features of the commensal house sparrows as compared to the non-commensal P. d. bactrianus can be understood as adaptations to process the larger and tougher seeds of domesticated cereals. The lachrymal and postorbital processes (both larger in commensal sparrows) are connected to the posterior region of the mandible through ligaments . The attachment of these ligaments together with the musculus depressor mandibulae, play an important role in the protraction of the upper jaw . The musculus depressor mandibulae, which depresses the lower jaw, originates in the occipital region of the cranium (again more developed in commensal sparrows) and inserts in the mandible . Another important muscle, the aductor mandibulae externus has its origin in a wide area of the posterior part of the skull, passes forwards and downwards between the postorbital process, and inserts in the posterior part of the mandible . The chewing area of the beak is stronger (more bony-built) in commensal house sparrows than in P. d. bactrianus and the larger nasal process in the commensal sparrows would provide more robustness in the kinetic area between the beak and the rest of the skull.
In conclusion, our results are consistent with the hypothesis that adaptations to food sources provided by human societies, and particularly cultivated cereals, have shaped the beak and skull of commensal house sparrows to become more robust and the facial cranium and beak to become elongated and more pointy. In a future study it would be interesting to study the feeding ecologies of the non-commensal P. d. bactrianus and the various commensal subspecies in greater detail in order to link ecology and morphology more directly. Further, it would be interesting to compare the anatomy of fossil house sparrow skulls with the various contemporary ones. House sparrow fossils are relatively common in archaeological excavation sites  as are cereal seeds . Hence, in principle, it should be possible to trace evolutionary changes in the morphology of the beak and skull of the house sparrow through time, from the rise of agriculture and beyond. Combining such studies with ontological observations may further improve our understanding of the developmental changes that has taken place as the house sparrow beak and skull became adapted to process cultivated cereal grains.
Schluter D: The Ecology of Adaptive Radiation. 2000, Oxford: Oxford University Press
Grant PR, Grant BR: Unpredictable evolution in a 30-year study of Darwin’s finches. Science. 2002, 296: 707-711. 10.1126/science.1070315.
Boag PT, Grant PR: Intense natural selection in a population of Darwin's finches (Geospizinae) in the Galápagos. Science. 1981, 214: 82-85. 10.1126/science.214.4516.82.
Schluter D, Grant PR: Determinants of morphological patterns in communities of Darwin's finches. Am Nat. 1984, 123: 175-196. 10.1086/284196.
Antonovics J, van Tienderen PH: Ontoecogenophyloconstraints? the chaos of constraint terminology. Trends Ecol Evol. 1991, 6: 166-168. 10.1016/0169-5347(91)90059-7.
Reiss MJ: The Allometry of Growth and Reproduction. 1989, Cambridge: Cambridge University Press
Brakefield PM: Evo-devo and constraints on selection. Trends Ecol Evol. 2006, 2006 (21): 362-368.
Benkman CV: Divergent selection drives the adaptive radiation of crossbills. Evolution. 2003, 57: 1176-1181.
Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ: The calmodulin pathway and evolution of elongated beak morphology in Darwin's finches. Nature. 2006, 442: 563-567. 10.1038/nature04843.
Hendrikse JL, Parsons TE, Hallgrímsson B: Evolvability as the proper focus of evolutionary developmental biology. Evol Devolop. 2007, 9: 393-409. 10.1111/j.1525-142X.2007.00176.x.
Eroukhmanoff F, Hermansen JS, Bailey R, Sæther SA, Sætre G-P: Local adaptation within a hybrid species. Heredity. 2013, doi:10.1038/hdy.2013.47
Summers-Smith JD: The Sparrows: A Study of the Genus Passer. 1988, Calton: T & AD Poyser
Anderson TR: Biology of the Ubiquitous House Sparrow: From Genes to Populations. 2006, Oxford: Oxford University Press
Vaurie C: Notes on some Ploceidae from western Asia (Notes from the Walter Koelz collections 1). Am Mus Nov. 1949, 1406: 1-41.
Vaurie C: Systematic notes on Palearctic birds, No. 24. Ploceidae: the genera Passer, Petronia and Montifringilla. Am Mus Nov. 1956, 1814: 1-27.
Johnston RF, Klitz WJ: Variation and evolution in a granivorous bird: the house sparrow. Granivorous Birds in Ecosystems. Edited by: Pinowski J, Kendeigh SC. 1977, Cambridge: Cambridge University Press, 15-51.
Sætre G–P, Riyahi S, Aliabadian M, Hermansen JS, Hogner S, Olsson U, Gonzales Rojas MF, Sæther SA, Trier CN, Elgvin TO: Single origin of human commensalism in the house sparrow. J Evol Biol. 2012, 25: 788-796. 10.1111/j.1420-9101.2012.02470.x.
Korelov MN, Kuzmina MA, Gavrilov EI, Kovshar AF, Gavrin G, Borodikhin IF: 1974. The Birds of Kazakhstan, vol. 5. 1974, Akad. Nauk Kazakhskoi SSR: Alma Ata, [in Russian]
Zarudny NA, Härms M: Bemerkungen über einige Vögel Persiens II, Die Sperlinge Persiens. J Ornithol. 1913, 1913 (61): 630-661. [In German]
Paludan K: On the Birds of Afghanistan. Vid Medd Dan Naturh For. 1959, 122: 1-332.
Ali S, Ripley SD: Handbook of the Birds of India and Pakistan, Volume 10. 1974, Bombay: Oxford University Press
Lush WM, Wien HC: The importance of seed size in early growth of wild and domesticated cowpeas. J Agric Sci. 1980, 94: 177-182. 10.1017/S0021859600028033.
Smith BD: The initial domestication of Cucurbita pepo in the Americas 10,000 years ago. Science. 1997, 276: 932-934. 10.1126/science.276.5314.932.
Salamini F, Özkan H, Brandolini A, Schäfer-Pregl R, Martin W: Genetics and geography of wild cereal domestication in the near east. Nature Rev. 2002, 3: 429-441.
Peng J, Ronin Y, Fahima T, Röder MS, Li Y, Nevo E, Korol A: Domestication quantitative trait loci in Triticum dicoccoides, the progenitor of wheat. Proc Natl Acad Sci USA. 2003, 100: 2489-2494. 10.1073/pnas.252763199.
Kilian B, Ozkan H, Pozzi C, Salamini F: Domestication of the Triticea in the fertile crescent. Genetics and Genomics of the Triticeae. Edited by: Feuillet C, Muehlbauer GJ. 2009, Dordrecht: Springer, 81-119.
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. 10.1111/j.0269-8463.2005.00923.x.
Foster DJ, Podos J, Hendry AP: A geometric morphometric appraisal of beak shape in Darwin’s finches. J Evol Biol. 2008, 21: 263-275.
Baylac M, Friess M: Fourier descriptors, Procrustes superimposition, and data dimensionality: an example of cranial shape analysis in modern human populations. Modern Morphometrics in Physical Anthropology. Edited by: Slice DE. 2005, New York: Kluwer Academic/Plenum Publishers, 145-165.
Rohlf FJ: tpsDig software, version 2.12. 2008, University of New York at Stony Brook: Department of Ecology and Evolution
Kuhl FP, Giardina CR: Elliptic Fourier features of a closed contour. Comp Graph Im Process. 1982, 18: 236-258. 10.1016/0146-664X(82)90034-X.
Hammer Ø, Harper DAT, Ryan PD: PAST: Paleontological Statictics software package for education and data analysis. Palaeont Electron. 2001, 4: 9-
Anderson MJ: A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 2001 (26): 32-46.
Mitteroecker P, Gunz P, Bernhard M, Schaefer K, Bookstein FL: Comparison of cranial ontogenetic trajectories among great apes and humans. J Human Evol. 2004, 46: 679-698. 10.1016/j.jhevol.2004.03.006.
Nunn CL, Barton RA: Comparative methods for studying primate adaptation and allometry. Evol Antrhropol. 2001, 10: 81-98. 10.1002/evan.1019.
Hollander J, Adams DC, Johannesson K: Evolution of adaptation through allometric shifts in a marine snail. Evolution. 2006, 60: 2490-2497.
Wiens JA, Dyer MI: Assessing the potential impact of granivorous birds in ecosystems. Granivorous Birds in Ecosystems. Edited by: Pinowski J, Kendeigh SC. 1977, Cambridge: Cambridge University Press, 205-266.
Baumel JJ, Raikow RJ: Arthrologia. Handbook of Avian Anatomy: Nomina Anatomica Avium. Edited by: Baumel JJ, King AS, Breazile JE, Evans HE, Vanden Berge JC. 1993, Cambridge, MA: Publications of the Nuttall Ornithological Club, 133-187. 2
Zusi R: The role of the depressor mandibulae muscle in kinesis of the avian skull. Proc US Nat Mus. 1967, 123: 1-28.
Van den Berge JC, Zweers GA: Myologia. Handbook of Avian Anatomy: Nomina Anatomica Avium. Edited by: Baumel JJ, King AS, Breazile JE, Evans HE, Vanden Berge JC. 1993, Cambridge, MA: Publications of the Nuttall Ornithological Club, 189-247. 2
Sims RW: The morphology of the head of the hawfinch (Coccothraustes coccothraustes). Bull Brit Mus (Nat Hist). 1955, 2: 371-393.
Ericson PGP, Tyrberg T, Kjellberg AS, Jonsson L, Ullén I: The earliest record of house sparrows (Passer domesticus) in northern Europe. J Archaeol Sci. 1997, 24: 183-190. 10.1006/jasc.1996.0102.
We thank R. I. Bailey, T. O. Elgvin, F. Eroukhmanoff, J. S. Hermansen, J. H. Hurum, K. L. Voje and two anonymous reviewers for comments on the manuscript and useful discussions. The Norwegian Research Council provided financial support.
The authors declare that they have no competing interests.
MA conceived of the study, data was collected by SR and TA under supervision of MA, analysis and interpretation of results were done by SR, ØH and GPS, the manuscript was written by GPS, ØH and SR with contributions from AS and CSR. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.