Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Beak and skull shapes of human commensal and non-commensal house sparrows Passer domesticus

Sepand Riyahi12, Øyvind Hammer3, Tayebeh Arbabi4, Antonio Sánchez5, Cees S Roselaar6, Mansour Aliabadian1 and Glenn-Peter Sætre7*

Author affiliations

1 Department of Biology, Faculty of Sciences, Ferdowsi University of Mashhad, Mashhad, Iran

2 Current address: Evolutionary Ecology Associate Research Unit (CSIC), Natural History Museum of Barcelona, Passeig Picasso s/n, Barcelona 08003, Spain

3 Natural History Museum, University of Oslo, Blindern, P. O. Box 1172, Oslo N-0318, Norway

4 Institute of Pharmacy and Molecular Biotechnology, Department of Biology, Im Neuenheimer Feld 364, Heidelberg D-69120, Germany

5 Catalan Institute of Paleontology, Campus de la UAB, Cerdanyola del Vallès 08193, Spain

6 Naturalis Biodiversity Center, Vertebrate Department, Darwinweg 2/PO Box 9517, RA Leiden 2300, Netherlands

7 Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biological Sciences, University of Oslo, Blindern, P. O. Box 1066, Oslo N-0316, Norway

For all author emails, please log on.

Citation and License

BMC Evolutionary Biology 2013, 13:200  doi:10.1186/1471-2148-13-200


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/13/200


Received:5 July 2013
Accepted:16 September 2013
Published:17 September 2013

© 2013 Riyahi et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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.

Conclusion

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.

Keywords:
Geometric morphometrics; Beak shape; Granivorous bird; Passer domesticus; Human commensalism

Background

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 [5]. 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. [17] 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 [17].

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 [17]. 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 [12]. 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 [17]. Further, the transition from a diet of mainly wild seeds to one consisting mainly of grains from cultivated cereals [12] 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 [26]. Moreover, domesticated cereal seeds remain encapsulated in a tough rachis that holds the seed together, whereas in wild cereals the rachis fragmentizes at ripening [24]. 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 [27].

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) [12]. 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 [11]. However, cultivated crops differ from their wild ancestors in consistent ways [24], 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.

Methods

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 [12] 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.

thumbnailFigure 1. Geographic distribution of four commensal subspecies of the house sparrow (green tones) and the non-commensal P. d. bactrianus (orange) in Iran. Sampling locations are indicated with black dots.

For geometric morphometrics on bird skulls and beaks, landmark analysis has been the most common approach (e.g. [28]). 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 [29]. 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 [30]. 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 [31] 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 [32], with the “invariant to starting position and rotation” option. Size was recorded by the length of the semimajor axis of the first harmonic ellipse [31]. 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.

thumbnailFigure 2. The average Elliptic Fourier Analysis spectrum across all specimens. The amplitude of each harmonic is calculated as the Euclidean norm of its four coefficients. The amplitude of the first harmonic (not shown) is 1.11.

The resulting 60 coefficients were subjected to further multivariate analysis in Past. For testing of multivariate group means we used a permutation test, PerMANOVA, [33] 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. [34]. 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.

Results

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).

thumbnailFigure 3. Sizes of skulls from the five subspecies of P. domesticus (means and two standard deviations).P. d. bactrianus is the only non-commensal subspecies.

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.

thumbnailFigure 4. Synthetic (predicted) skull shapes for two different skull sizes, based on linear regression of the first 15 EFA harmonics upon size. The skulls are normalized to equal size.

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).

thumbnailFigure 5. Average shape of the non-commensal P. d. bactrianus versus the commensal subspecies. Sizes are normalized.

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.

thumbnailFigure 6. CAC versus RSC 2 scores for the five subspecies of P. domesticus. Reconstructed skull outlines are shown in their respective positions in CAC-RSC 2 space.

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 [34]. 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.

thumbnailFigure 7. Log semimajor axis size versus CAC scores for the five subspecies of P. domesticus.

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).

thumbnailFigure 8. Photograph (dorsal view) of a non-commensal house sparrow P. d bactrianus (upper panel) and a similar sized commensal house sparrow P. d. persicus (lower panel). Notice the difference in robustness and elongation of the facial cranium, and in particular a) the postorbital process, b) the lachrymal process, c) the oral surface of the maxillare as seen through the nostrils, d) the nasal process and e) the occipital region.

Discussion

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 [17]. 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 [37]. 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 [12]. 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 [38]. The attachment of these ligaments together with the musculus depressor mandibulae, play an important role in the protraction of the upper jaw [39]. 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 [40]. 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 [41]. 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.

Conclusions

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 [42] as are cereal seeds [24]. 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.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

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.

Acknowledgements

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.

References

  1. Schluter D: The Ecology of Adaptive Radiation. Oxford: Oxford University Press; 2000. OpenURL

  2. Grant PR, Grant BR: Unpredictable evolution in a 30-year study of Darwin’s finches.

    Science 2002, 296:707-711. PubMed Abstract | Publisher Full Text OpenURL

  3. Boag PT, Grant PR: Intense natural selection in a population of Darwin's finches (Geospizinae) in the Galápagos.

    Science 1981, 214:82-85. PubMed Abstract | Publisher Full Text OpenURL

  4. Schluter D, Grant PR: Determinants of morphological patterns in communities of Darwin's finches.

    Am Nat 1984, 123:175-196. Publisher Full Text OpenURL

  5. Antonovics J, van Tienderen PH: Ontoecogenophyloconstraints? the chaos of constraint terminology.

    Trends Ecol Evol 1991, 6:166-168. PubMed Abstract | Publisher Full Text OpenURL

  6. Reiss MJ: The Allometry of Growth and Reproduction. Cambridge: Cambridge University Press; 1989. OpenURL

  7. Brakefield PM: Evo-devo and constraints on selection.

    Trends Ecol Evol 2006, 2006(21):362-368. OpenURL

  8. Benkman CV: Divergent selection drives the adaptive radiation of crossbills.

    Evolution 2003, 57:1176-1181. PubMed Abstract OpenURL

  9. 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. PubMed Abstract | Publisher Full Text OpenURL

  10. Hendrikse JL, Parsons TE, Hallgrímsson B: Evolvability as the proper focus of evolutionary developmental biology.

    Evol Devolop 2007, 9:393-409. Publisher Full Text OpenURL

  11. 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

    OpenURL

  12. Summers-Smith JD: The Sparrows: A Study of the Genus Passer. Calton: T & AD Poyser; 1988. OpenURL

  13. Anderson TR: Biology of the Ubiquitous House Sparrow: From Genes to Populations. Oxford: Oxford University Press; 2006. OpenURL

  14. Vaurie C: Notes on some Ploceidae from western Asia (Notes from the Walter Koelz collections 1).

    Am Mus Nov 1949, 1406:1-41. OpenURL

  15. Vaurie C: Systematic notes on Palearctic birds, No. 24. Ploceidae: the genera Passer, Petronia and Montifringilla.

    Am Mus Nov 1956, 1814:1-27. OpenURL

  16. Johnston RF, Klitz WJ: Variation and evolution in a granivorous bird: the house sparrow. In Granivorous Birds in Ecosystems. Edited by Pinowski J, Kendeigh SC. Cambridge: Cambridge University Press; 1977:15-51. OpenURL

  17. 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. PubMed Abstract | Publisher Full Text OpenURL

  18. Korelov MN, Kuzmina MA, Gavrilov EI, Kovshar AF, Gavrin G, Borodikhin IF: 1974. The Birds of Kazakhstan, vol. 5. Akad. Nauk Kazakhskoi SSR: Alma Ata; 1974.

    [in Russian]

    OpenURL

  19. Zarudny NA, Härms M: Bemerkungen über einige Vögel Persiens II, Die Sperlinge Persiens.

    J Ornithol 1913, 1913(61):630-661.

    [In German]

    OpenURL

  20. Paludan K: On the Birds of Afghanistan.

    Vid Medd Dan Naturh For 1959, 122:1-332. OpenURL

  21. Ali S, Ripley SD: Handbook of the Birds of India and Pakistan, Volume 10. Bombay: Oxford University Press; 1974. OpenURL

  22. Lush WM, Wien HC: The importance of seed size in early growth of wild and domesticated cowpeas.

    J Agric Sci 1980, 94:177-182. Publisher Full Text OpenURL

  23. Smith BD: The initial domestication of Cucurbita pepo in the Americas 10,000 years ago.

    Science 1997, 276:932-934. Publisher Full Text OpenURL

  24. 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. OpenURL

  25. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Kilian B, Ozkan H, Pozzi C, Salamini F: Domestication of the Triticea in the fertile crescent. In Genetics and Genomics of the Triticeae. Edited by Feuillet C, Muehlbauer GJ. Dordrecht: Springer; 2009:81-119. OpenURL

  27. 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. Publisher Full Text OpenURL

  28. Foster DJ, Podos J, Hendry AP: A geometric morphometric appraisal of beak shape in Darwin’s finches.

    J Evol Biol 2008, 21:263-275. PubMed Abstract | Publisher Full Text OpenURL

  29. Baylac M, Friess M: Fourier descriptors, Procrustes superimposition, and data dimensionality: an example of cranial shape analysis in modern human populations. In Modern Morphometrics in Physical Anthropology. Edited by Slice DE. New York: Kluwer Academic/Plenum Publishers; 2005:145-165. OpenURL

  30. Rohlf FJ: tpsDig software, version 2.12. University of New York at Stony Brook: Department of Ecology and Evolution; 2008. OpenURL

  31. Kuhl FP, Giardina CR: Elliptic Fourier features of a closed contour.

    Comp Graph Im Process 1982, 18:236-258. Publisher Full Text OpenURL

  32. Hammer Ø, Harper DAT, Ryan PD: PAST: Paleontological Statictics software package for education and data analysis.

    Palaeont Electron 2001, 4:9. OpenURL

  33. Anderson MJ: A new method for non-parametric multivariate analysis of variance.

    Austral Ecol 2001, 2001(26):32-46. OpenURL

  34. 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. Publisher Full Text OpenURL

  35. Nunn CL, Barton RA: Comparative methods for studying primate adaptation and allometry.

    Evol Antrhropol 2001, 10:81-98. Publisher Full Text OpenURL

  36. Hollander J, Adams DC, Johannesson K: Evolution of adaptation through allometric shifts in a marine snail.

    Evolution 2006, 60:2490-2497. PubMed Abstract OpenURL

  37. Wiens JA, Dyer MI: Assessing the potential impact of granivorous birds in ecosystems. In Granivorous Birds in Ecosystems. Edited by Pinowski J, Kendeigh SC. Cambridge: Cambridge University Press; 1977:205-266. OpenURL

  38. Baumel JJ, Raikow RJ: Arthrologia. In Handbook of Avian Anatomy: Nomina Anatomica Avium. 2nd edition. Edited by Baumel JJ, King AS, Breazile JE, Evans HE, Vanden Berge JC. Cambridge, MA: Publications of the Nuttall Ornithological Club; 1993:133-187. OpenURL

  39. Zusi R: The role of the depressor mandibulae muscle in kinesis of the avian skull.

    Proc US Nat Mus 1967, 123:1-28. OpenURL

  40. Van den Berge JC, Zweers GA: Myologia. In Handbook of Avian Anatomy: Nomina Anatomica Avium. 2nd edition. Edited by Baumel JJ, King AS, Breazile JE, Evans HE, Vanden Berge JC. Cambridge, MA: Publications of the Nuttall Ornithological Club; 1993:189-247. OpenURL

  41. Sims RW: The morphology of the head of the hawfinch (Coccothraustes coccothraustes).

    Bull Brit Mus (Nat Hist) 1955, 2:371-393. OpenURL

  42. 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. Publisher Full Text OpenURL