Skip to main content

Evolution of body morphology and beak shape revealed by a morphometric analysis of 14 Paridae species



Morphological characters of birds reflect their adaptive evolution and ecological requirements and are also relevant to phylogenetic relationships within a group of related species. The tits (Paridae) are known to be outwardly homogeneous in shape, with one aberrant member, the Ground Tit (Pseudopodoces humilis), which is quite different from its relatives in both body morphology and beak shape. We combined traditional measurements and geometric morphometrics to quantify the variation in body morphology and beak shape of 14 Paridae species distributed in China. Based on these results, we sought to assess the contribution of phylogeny, altitude and species interactions to the evolution of morphological traits.


The basic features for discriminating among the 14 species studied here were overall body size, the ratio of body and tail length to culmen and tarsus length, and beak shape (long/slender/pointy vs. short/robust/blunt). These dimensions clearly separate Ps. humilis and Melanochlora sultanea from the other species in shape space. Body length and PC3 of beak shape (round outline vs. straight outline) show significant phylogenetic signals. Across 14 species, altitude is related to tarsus, culmen length and PC1 of beak shape. Within Parus major, altitude is related to body weight, body length, culmen length and PC1 of body morphology. Morphological distances and geographic distances among species are positively correlated.


The body morphology of Paridae species shows extensive evolutionary changes, while their beak has mainly evolved along the long/slender/pointy vs. short/robust/blunt dimension. Only body length and beak curvature show a phylogenetic signal. Altitude correlates with multiple traits both across and within species, suggesting that altitude is an important factor in promoting morphological divergence. The deviant appearance of Ps. humilis corresponds to its foraging and feeding adaptations to high-altitude steppe habitats. Our results also show a higher level of morphological divergence with greater difference in distribution ranges among the Paridae species involved in this study.


Diversity among a group of related organisms usually shows a pattern of descent with modification, with both phylogenetic history and natural selection influencing a species’ form [1, 2]. A bird’s body size and shape are related to its migratory habits, flight mode, habitat use, and foraging behaviour, as well as the effects of sexual selection [35]. Bird beaks are considered relevant to foraging modes, food types, parental care and song structure [68]. Morphological divergence, which enables related organisms to avoid competition and exploit available niches, plays an important role in adaptive radiation [9, 10]. Previous studies have shown that rapid diversification in adaptive radiations usually involve accelerated rates of morphological evolution [11, 12]. Therefore, morphological characters are frequently used to investigate adaptation radiation, such as beak shape adapting to different food sources in Darwin’s finches [13], extremely diverse beak morphology of the Hawaiian honeycreepers [14], and the speciation bursts of Malagasy vangas related to morphological diversification [10].

The avian Family Paridae, which includes c. 55 species that are distributed across the northern hemisphere, tropical Africa and Indonesia and are found from sea level to 5500 m [15, 16], is a well-investigated example of adaptive radiation. The hotspot of Paridae species richness is found in China, which is also considered as the centre of origin of this group and plays a role as a “cradle of diversity” [17]. Previous studies have proposed that the diversification of Paridae is primarily driven by vicariance related to orogenesis in the Himalayas and Mountains of Southwest China concurrent with the uplift of the Qinghai-Tibet Plateau (QTP) [17, 18]. Paridae species have diverse plumage within their large family but outwardly appear rather uniform in shape, except for the Sultan Tit (Melanochlora sultanea), the Yellow-browed Tit (Sylviparus modestus) and the Ground Tit (Pseudopodoces humilis) [16, 19]. Pseudopodoces humilis is endemic to the QTP and dwells exclusively on the treeless steppes, which has been identified as an extreme example in the evolution of Paridae species [16, 20]. It had long been misclassified as a small species of ground jay due to its strikingly aberrant appearance, while recent studies have revealed that this species is closely related to Parus major and P. monticolus, with a divergent time of approximately 7.7 - 9.9 million years [2123]. The extremely deviant morphology of Ps. humilis has received extensive attention, and qualitative descriptions of its morphological differences from other parids have been documented [21]. In addition, univariate analyses of linear characters, such as body size and beak length, have been performed on a few Paridae species to investigate interspecific interactions, foraging behaviour and niche shifts [24, 25]. However, few studies have quantified morphological variation across the main clades and investigated potential factors that could shape the patterns of diversity in this family.

In this study, we first aimed to quantify the variation in body morphology and beak shape of 14 Paridae species distributed in China. We applied traditional morphometric methods to analyse body morphology. For beak morphology, we first measured culmen length, but linear distances such as beak length, width and depth are insufficient to fully describe beak shape [2628]. Geometric morphometric methods can capture the explicit geometry of a morphological structure by examining associations among an entire set of landmarks and has the additional advantages of improved statistical power and fewer a priori assumptions about what should be measured [27, 29, 30]. These methods have increasingly become a powerful and useful tool to study morphological variation and adaptive radiation, especially when investigating features such as bird beaks [30] and skulls [8, 31, 32]. Therefore, we applied the improved inferential resolution of geometric morphometric methods, through which we attempted to present visible and interpretable variations in beak shape.

We then examined the factors that could potentially affect morphological variation. Firstly, we tested for phylogenetic signal in each morphological trait to evaluate whether species that descended from a recent common ancestor closely resembled each other. Secondly, we assessed the correlation between morphology and altitudinal distribution. Altitude is associated with predictable changes in temperature, precipitation and habitat type, which result in different selective pressures and consequently influence the fauna that are present along altitudinal gradients [33, 34]. Adaptation to different altitudes may result in the divergence of phenotypes between populations and eventually contribute to speciation. To test this hypothesis, we looked for morphological variation along the altitudinal gradient both among Paridae species and within P. major, which occurs from 400 m to as high as 4400 m above sea level. Lastly, we examined the correlation between character divergence and distribution overlap to test the influence of species interactions on morphology. Ecological differences are essential to coexistence, and closely related, sympatric species need to differ in at least one niche dimension of habitat selection, prey size and feeding method to avoid competition [3537]. Ecological variation along these dimensions is associated with predictable morphological changes [36]. Thus, we expected greater character divergence among Paridae species with greater overlap in their distribution range.


Data collection

We measured 376 museum specimens belonging to 14 species (see Additional file 1: Table S1 for a species list) from the National Zoological Museum of China, Institute of Zoology, Chinese Academy of Sciences. Data for body weight and body length were retrieved from the original records taken during field collection. Wing length, tail length, tarsus length and culmen length were measured to the nearest 0.1 mm with a digital caliper (see Additional file 1: Table S2 for detailed methods of measurement and Additional file 2 for a complete description of the data). For geometric morphometric analyses, the beak profiles of 306 specimens were captured using an image collecting system consisting of a Keyence VHX-1000C microscope and a Nikon Ni-E camera. Only adult males were included. All individuals of P. major were from minor section [19].

Body morphology

Before statistical analyses, all measurements were log-transformed to normalize their distributions [4]. We performed one-way ANOVAs to compare each character among the species. In addition, we conducted a canonical variate analysis (CVA) on the data of all individuals to extract axes that best discriminated among the groups and to generate a matrix of pairwise Mahalanobis distances [38].

A principal component analysis (PCA) was applied to species mean values of each character to reduce the number of variables and visualize the variation. Before the PCA, all data were standardized to a mean of zero and a variance of one to minimize the effect of different initial units. Body size is important for evolutionary and ecological studies [39]; therefore, the data were not size-corrected before the PCA to preserve the information of body size.

Beak shape variation

To characterize the shape of the beaks, we analysed the profile of part of the upper mandible from the nares to the tip because the whole lateral view of the beak is difficult to capture due to the varied rictal bristles and feather coverage. TpsDig [40] was used to place 3 landmarks and 18 semi-landmarks [29] along the outline of the beak (see Fig. 1). The semi-landmarks, which were important for quantifying the shape of the beak, which lacks clear homologous points, were placed equidistantly. We used tpsRelw [41] to slide the semi-landmarks and then output the aligned specimens (see Additional file 3 for the txt file of aligned specimens), which was later analysed for shape variation using MorphoJ [42]. The CVA was performed on all configurations to extract the axes with the greatest interspecific differences and generate a matrix of pairwise Mahalanobis distances [38]. A Procrustes ANOVA [43] was conducted as extra effects. The PCA was performed on species mean shapes to visualize the shape changes that accounted for most of the variation.

Fig. 1
figure 1

Landmarks and semi-landmarks used for the geometric morphometric analysis. A line perpendicular to the suture was drawn across the rostral edge of the nares. Two landmarks were placed where this line intersects the outline of the upper mandible, whereas the third was placed at the tip of the beak. Nine semi-landmarks were placed equidistantly between LM1 and LM3, and the other nine between LM2 and LM3

Geometric morphometrics use superimposition methods to eliminate non-shape variation in the configurations of the landmarks by overlaying them according to some optimization criterion [26]. In this study, we applied Procrustes generalized least squares. This method first translates the centroid of each configuration to the origin, then divides the configurations by centroid size to scale them to a common unit of size and finally optimally rotates each configuration to minimize the squared differences between corresponding landmarks [44, 45]. This procedure eliminates the effects of position, scale and orientation so that the Procrustes coordinates of landmarks of a configuration can represent size-free shape, while the centroid size represents overall beak size.

Phylogenetic and comparative approach

To serve as the backbone for a comparative study, a phylogenetic tree needs to show a branching pattern and estimates of the branch lengths, which are assumed to be proportional to the expected variance in the amount of evolutionary change along each branch under a Brownian-motion model [46]. The phylogeny of Paridae has been extensively investigated in recent decades, and some taxonomic revisions have been proposed [15, 47, 48]. In this study, we followed the complete phylogeny of Paridae presented by Johansson et al. [47] and constrained the topology of the 14 species involved. The branch lengths were then calculated using two nuclear genes (Myoglobin and Ornithine Decarboxylase) and two mitochondrial genes (NADH dehydrogenase subunit 2 and Cytochrome b). We used MrBayes 3.2 [49] for the phylogenetic construction with gene partitions, using best-fit models based on BIC model selection criteria (GTR + I + G, TIM + I + G, K80 + G and HKY + G for Cytb, ND2, Myo, and ODC, respectively). The sequences were downloaded from GenBank (see Additional file 1: Table S3 for accession numbers and references).

The reconstructed phylogenetic tree (see Additional file 4) was then projected into the tangent shape space of the PC scores computed from the mean beak shape of each terminal taxon to show the trajectories of specific lineages and clades. The ancestral states of internal nodes were reconstructed based on squared-change parsimony [50], which is identical to assuming Brownian motion.

Blomberg’s K [51] and Pagel’s λ [52] were used to assess the phylogenetic signal in morphological traits. For both indices, a value close to 0 indicates phylogenetic independence and a value of 1 indicates that traits are evolving according to Brownian motion in the given phylogeny. Values that are higher than 1 indicate that traits are more similar amongst species than would be expected due to a Brownian-motion model [5153]. The residuals of regressions of log-transformed linear traits on log-transformed body weight values were used as size-corrected traits [35]. We tested for phylogenetic signal in log-transformed linear traits, size-corrected traits and PC values calculated from body morphology and beak shape.

Correlation between morphology and altitude

We used regression analyses to investigate covariation between morphological traits and altitudinal distribution. Analyses were carried out at two hierarchical levels: 1) variation among the 14 species and 2) variation within P. major. For determining interspecific covariation, a regression of species mean values of each morphological trait on the mid-points of species’ altitudinal ranges was conducted. And phylogenetic generalized least square (PGLS) means were applied to correct for phylogenetic relationships [54]. The altitudinal ranges of each species were determined based on values reported in the literature [19, 55]. As for the intraspecific covariation, we used data of individuals within P. major because the large sample size of this species enabled the intraspecific analyses to be conducted. The altitudinal data of each P. major specimen was collected from the original records of the sampling sites. For specimens that lacked records of the sampling sites, altitudinal data were substituted with the mid-points of the altitudinal ranges of populations. For both interspecific and intraspecific analyses, each morphological character was separately correlated with altitude. We conducted each regression twice and then used the mean value of these replicates as the final result to minimize type I errors.

Character divergence and distribution overlap

We conducted a partial Mantel test [56] (10,000 permutations) to look for possible relationships between character divergence and the degree of sympatry between species. The test was performed on the matrices of pairwise Mahalanobis distances of the morphological characters (both body morphology and beak shape) and the matrices of distributional distance, with the phylogenetic distances included as controlled variables. Both geographic distribution and altitudinal distribution were investigated.

As all specimens we measured were from China, their morphological variation may be insufficient to explain global patterns due to the considerable intraspecific differences in morphology between populations in China and those in other areas. Therefore, for each species, we used only the distribution range within China. We drew maps of the geographical distributions for all species on a Global Information System (GIS) (ArcMap 10.2) and calculated the degree overlap using the following equation: DO = OAB / (SA + SB – OAB), where SA and SB are the areas of the ranges of species A and species B, and OAB is the area in which their ranges overlap. Furthermore, DAB = 1 – DO was used to define the geographical distance. Distribution data were extracted from BirdLife International and NatureServe [57]. The degree of altitudinal overlap between species was also calculated using this method.

We used PAST 2.17 [58] to perform the PCA, ANOVA, ANCOVA, CVA of body morphology and regression analysis; MorphoJ [42] for the geometric morphometric analysis of variation in beak shape; the package ape [59], geiger [60] and phytools [61] in R for phylogenetic signal test; the package nlme [62] in R for the PGLS analysis; and the package vegan [63] in R for the Mantel test. The phylogenetic signal test and the Mantel test were performed separately on the dataset including all 14 species and on the dataset without Ps. humilis.


Variation in body morphology

One-way ANOVAs of the linear measurements show significant interspecific differences for all six characters (see Additional file 1: Table S4 for summary statistics and Additional file 1: Figure S1 for box plots). A CVA conducted on all individuals shows significant Mahalanobis distances between every pair of species (Additional file 1: Table S5 and Figure S2). The results of both tests indicate a differentiation in body morphology among the Paridae species.

The PCA performed on species mean values of linear characters shows that most of the variation is concentrated in the first two axes (Additional file 1: Table S6). PC1 explains 84.13 % of the total variance and is interpreted as the enlargement of all six characters. PC2 explains 11.06 % of the total variance and is interpreted mainly as changes in the ratio of the body and tail length to the culmen and tarsus length (Fig. 2a). Hence, PC1 can be considered an indicator of overall body size, and PC2 represents differences in the arrangement of body parts. The first two PC axes separate Ps. humilis (a large-bodied tit with a long tarsus and culmen) and M. sultanea (a large-bodied tit with a long body and tail) from the main scatter, but the rest of the species largely cluster together (Fig. 2a).

Fig. 2
figure 2

Scatters of the principal component analysis of the morphology of the 14 species examined. a: Scatters of the PC values for body morphology. b: Scatters of the PC values for beak shape. Phylogeny is projected onto the PC values in each plot. c: Phylogeny of the Paridae, including the 14 species examined in this study. The species illustrations refer to HBW alive:

Variation in beak shape

The result of the Procrustes ANOVA shows significant interspecific differences in both shape [F (494, 12096) = 55, P < 0.0001] and size [F (13, 292) = 110, P < 0.0001] of beak morphology. The CVA shows significance for most pairwise Mahalanobis distances (see Additional file 1: Table S7), suggesting differentiation in beak shape among the Paridae species.

The PCA conducted on the covariance matrix of Procrustes coordinates concentrates more than 99 % of the total variance in the first 3 axes (Fig. 3). PC1 can be interpreted as the variation between a long, slender, pointed beak shape and a short, robust, blunt shape. PC2 is associated with the relative position of the posterior-most point of the upper and lower curves of the beak, i.e., LM1 and LM2 (Fig. 1). PC3 relates to variation between beaks with straight, sharp, triangular outlines and beaks with relatively rounder and decurved outlines (Fig. 3). Because LM1 and LM2 were artificially defined rather than strictly homologous, PC2 is likely to be related to the abrasion of the horny surface on the beaks and the preserved condition of the museum specimens. In addition, this dimension of shape variation shows less interspecific difference than the other dimensions (Additional file 1: Figure S2), which makes species mean values less relevant. Therefore, the following analyses and illustrations concentrate on PC1 and PC3. These two axes separate Ps. humilis from the main group, whereas the remaining species form a dense and narrow cluster in the tangent shape space (Fig. 2b).

Fig. 3
figure 3

Patterns of shape change associated with PCs calculated from beak shape variation

Evolutionary changes and phylogenetic signals

The projections of the phylogenetic tree into the PC plots show some visible clusters that correspond to groups of related taxa (Fig. 2). Sylviparus and Melanochlora, as the two oldest lineages of Paridae, are connected to their sister node with long branches. On the other hand, the plots show extensive crossings of branches and some relatively long branches between related species (see Ps. humilis).

Significant phylogenetic signals are found in log-transformed body length and PC3 of beak shape (Additional file 1: Table S8). After correcting for overall body size, body length loses its phylogenetic signal. In addition, when Ps. humilis, a morphological outlier, is excluded from the analysis, body weight and size-corrected culmen length show significant phylogenetic signals.

Covariation between morphology and altitude

The results of the PGLS analysis show that across the 14 species analysed, altitude is significantly correlated with tarsus and culmen length and with the short/blunt/robust vs. long/slender/pointed dimension of beak shape variation (Fig. 4 and Additional file 1: Table S9). Within individual P. major specimens, altitude is correlated with multiple morphological traits (Fig. 5 and Additional file 1: Table S9). Individuals at higher elevations tend to have a greater body weight and larger overall body size, as well as a longer body and culmen (Fig. 5).

Fig. 4
figure 4

Interspecific covariation between morphological traits and altitude of the 14 species examined. Correlations of altitude with tarsus and culmen length and PC1 of beak shape are shown. Dots with different colours and shapes represent different species. Lines represent phylogenetic generalized least square (PGLS) regressions

Fig. 5
figure 5

Intraspecific covariation between morphological traits and altitude within P. major. Correlations of altitude with body weight, body length, culmen length and PC1 of body morphology are shown. Each dot represents an individual P. major specimen. Lines represent ordinary least square (OLS) regressions

Character divergence and distribution overlap

Mantel tests detected a significant positive correlation between the morphological distance and the distributional distance (both geographic and altitudinal) after controlling for phylogenetic relationships (Additional file 1: Table S10). The results indicate that Paridae species that overlap more in their distributional range tend to have a higher level of morphological similarity than is observed for species with greater distributional differences. This pattern remains when Ps. humilis is excluded from the analysis.


In this study we quantified morphological variation among 14 species from the main clades of Paridae using both traditional and geometric morphometrics. We also assessed the contribution of phylogenetic relatedness, altitudinal distribution and/or species interactions to morphological evolution. Two major aspects of variation in body morphology were identified: overall body size and body shape (the ratio of body and tail length to tarsus and culmen length). Beak shape variation was concentrated on the long/slender/pointed vs. short/robust/blunt axis, which is in agreement with previous studies on the beak morphology of Darwin’s finches [30, 64]. In addition, we extracted other dimensions of shape change, such as the curvature of the upper mandible. Both body morphology and beak shape identify Ps. humilis and M. sultanea as obvious outliers, which can be visualized in the shape space (Fig. 2). Typical “tit-like” parids are known to be outwardly uniform in shape [16, 19], but our quantitative methods separated most of the species based on body morphology (Additional file 1: Figure S2A). However, most species do overlap in the beak variation shape space (Additional file 1: Figure S2C). Beak shape appears to vary mostly along the long/slender/pointed vs. short/robust/blunt axis, which discriminates Ps. humilis from the main cluster (Fig. 2). Paridae species are relatively conservative in beak shape, which is in contrast to other adaptive radiations that harbour a greater diversity of beak morphology, such as Darwin’s finches [13], Malagasy vangas [10] and Hawaiian honeycreepers [14].

Only body length and the curvature of the beak show a significant phylogenetic signal. The signal in body length is lost after size correction, which indicates that the signal associated with body length is primarily due to overall body size. Regarding beak morphology, the curvatures of beak outline are more phylogenetically conserved than the slenderness of the beak. The exclusion of Ps. humilis increases the phylogenetic signal in the morphological characters, such as body size and beak morphology, which suggests that the phylogenetic signals in morphological traits may be interfered by the discrepancy between the striking morphological divergence and close phylogenetic relationships. The morphological difference between Ps. humilis and its closest relatives (P. major and P. monticolus) is larger than that among other genera of Paridae [21].

Altitude is significantly correlated with tarsus length and beak morphology across the 14 species in this study. In addition, there is an intraspecific body-size cline along the altitudinal gradient within P. major, a species that covers a wide altitudinal range. Both results suggest that morphological traits of Paridae covary with altitudinal distributions. Altitude relates to the range of habitat occupied [65, 66] and acts as a selective power on morphology. The long/slender/pointed and short/blunt/robust beak morphs reflect foraging behaviours of probing on ground and pecking in wooded habitats, respectively [8]. In addition, variation in tarsus length, which was also discovered within core Corvoidea and Phylloscopus warblers [65, 66], is also considered to reflect adaptations to foraging in more open habitats at higher elevations, especially hopping and walking on the ground. The intraspecific size cline may be a response to the temperature variation that occurs along altitudinal gradients [67, 68].

Both geographical separation and ecological differences between populations contribute to speciation in birds [69]. Here, we speculate that altitude is an ecological factor that provides different selective pressures and promote the divergence of Paridae. Pseudopodoces humilis, an extremely morphologically deviant species within Paridae, is an example of speciation promoted by altitude and related habitats. The uplift of the QTP created the high steppes habitat and provided an “open” niche, which triggered the morphological evolution of the ancestor of Ps. humilis [21, 35]. Rather than being similar to its forest-living parid relatives, it is instead morphologically convergent to ground jays [21] and mynas, which are also predominantly open-country ground foragers [35]. The deviant appearance of Ps. humilis reflects adaptations to a novel adaptive zone [21, 35], including a larger body size for preserving heat; increased tarsal length for terrestrial locomotion; and a long, slender, decurved beak for probing, pecking and digging vigorously in soft earth, turf, yak dung and decaying animal corpses to find small invertebrates and larvae [19, 20]. The different selective pressures of the high steppes and the forests drive the divergence between Ps. humilis and other parids.

Closely related sympatric species usually need to diverge in ecological niches to avoid resource competition [24, 70, 71]. Previous studies have discovered niche divergence in sympatric Paridae species [72, 73] and character displacement in some local communities [39, 74]. This is supported by our results which discovered significant differences for body morphology and/or beak shape among species (Additional file 1: Table S5 and S7), albeit their superficial similarity. On the other hand, the results of partial Mantel test found further morphological divergence with distributional differences at a larger scale (Additional file 1: Table S10). Paridae shows patterns of morphological divergence related to differences in distribution ranges, likely driven by adaptation to different habitats and selective pressures at different ranges.


In this study, we quantified the morphological variation of 14 Paridae species distributed in China. The basic features for discrimination are overall body size, the ratio of body and tail length to culmen and tarsus length, and beak shape (long/slender/pointy vs. short/robust/blunt). These dimensions identify Ps. humilis and M. sultanea as obvious outliers. Body morphology shows extensive variation and can separate most species, while beak shape has evolved mostly along the long/slender/pointy vs. short/robust/blunt dimension. Only body length and beak curvature show a phylogenetic signal. The low level of phylogenetic signal in the morphological data may be due to the close phylogenetic relationship and considerable morphological divergence between Ps. humilis and its closest relatives. Altitude is correlated with multiple traits both across and within species. We speculate that the altitudinal gradient is an important factor promoting morphological divergence. The deviant appearance of Ps. humilis, including its large body size, long tarsus and long, slender, decurved beak, reflects its adaptation to the high-altitude steppes created by the uplift of the Qinghai-Tibet Plateau. We also show here that Paridae species with greater distributional difference tend to have a higher level of morphological divergence, which is probably due to different selective pressures provided by habitats at different ranges.


  1. Burns KJ, Hackett SJ, Klein NK. Phylogenetic relationships and morphological diversity in Darwin's finches and their relatives. Evolution. 2002;56(6):1240–52.

    Article  PubMed  Google Scholar 

  2. Grenier JL, Greenberg R. A biogeographic pattern in sparrow bill morphology: parallel adaptation to tidal marshes. Evolution. 2005;59(7):1588–95.

    Article  PubMed  Google Scholar 

  3. Marchetti K, Price T, Richman A. Correlates of wing morphology with foraging behavior and migration distance in the Genus Phylloscopus. J Avian Biol. 1995;26(3):177–81.

    Article  Google Scholar 

  4. Kaboli M, Aliabadian M, Guillaumet A, Roselaar CS, Prodon R. Ecomorphology of the wheatears (genus Oenanthe). Ibis. 2007;149(4):792–805.

    Article  Google Scholar 

  5. Zeffer A, Johansson LC, Marmebro A. Functional correlation between habitat use and leg morphology in birds (Aves). Biol J Linn Soc. 2003;79(3):461–84.

    Article  Google Scholar 

  6. Forstmeier W, Leisler B, Kempenaers B. Bill morphology reflects female independence from male parental help. P Roy Soc B-Biol Sci. 2001;268(1476):1583–8.

    Article  CAS  Google Scholar 

  7. Badyaev AV, Young RL, Oh KP, Addison C. Evolution on a local scale: developmental, functional, and genetic bases of divergence in bill form and associated changes in song structure between adjacent habitats. Evolution. 2008;62(8):1951–64.

    Article  PubMed  Google Scholar 

  8. Kulemeyer C, Asbahr K, Gunz P, Frahnert S, Bairlein F. Functional morphology and integration of corvid skulls - a 3D geometric morphometric approach. Front Zool. 2009; doi:10.1186/1742-9994-6-2

  9. Schluter D. The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.

    Google Scholar 

  10. Jonsson KA, Fabre PH, Fritz SA, Etienne RS, Ricklefs RE, Jorgensen TB, et al. Ecological and evolutionary determinants for the adaptive radiation of the Madagascan vangas. P Natl Acad Sci USA. 2012;109(17):6620–5.

    Article  CAS  Google Scholar 

  11. Rabosky DL, Adams DC. Rates of morphological evolution are correlated with species richness in salamanders. Evolution. 2012;66(6):1807–18.

    Article  PubMed  Google Scholar 

  12. Rabosky DL, Santini F, Eastman J, Smith SA, Sidlauskas B, Chang J, et al. Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation. Nat Commun. 2013;4:1958.

    Article  PubMed  Google Scholar 

  13. Grant PR. Ecology and evolution of Darwin's finches. Princeton: Princeton University Press; 1986.

    Google Scholar 

  14. Lerner HRL, Meyer M, James HF, Hofreiter M, Fleischer RC. Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers. Curr Biol. 2011;21(21):1838–44.

    Article  CAS  PubMed  Google Scholar 

  15. Gill FB, Slikas B, Sheldon FH. Phylogeny of titmice (Paridae): II. Species relationships based on sequences of the mitochondrial cytochrome-b gene. Auk. 2005;122(3):121–43.

    Article  Google Scholar 

  16. Harrap S, Quinn D. Tits, Nuthatches and Treecreepers. London: Helm; 1990.

    Google Scholar 

  17. Tietze DT, Borthakur U. Historical biogeography of tits (Aves: Paridae, Remizidae). Org Divers Evol. 2012;12(4):433–44.

    Article  Google Scholar 

  18. Päckert M, Martens J, Sun YH, Severinghaus LL, Nazarenko AA, Ting J, et al. Horizontal and elevational phylogeographic patterns of Himalayan and Southeast Asian forest passerines (Aves: Passeriformes). J Biogeogr. 2012;39(3):556–73.

    Article  Google Scholar 

  19. Gosler AG, Clement P. Family Paridae (Tits and Chickadees). In: del Hoyo J, Elliott A, Christie DA, editors. Handbook of the Birds of the World. Barcelona: Lynx Editions; 2007.

    Google Scholar 

  20. Lei FM, Lu TC. China endemic birds. Beijing: Science Press; 2006.

    Google Scholar 

  21. James HF, Ericson PGP, Slikas B, Lei FM, Gill FB, Olson SL. Pseudopodoces humilis, a misclassified terrestrial tit (Paridae) of the Tibetan Plateau: evolutionary consequences of shifting adaptive zones. Ibis. 2003;145(2):185–202.

    Article  Google Scholar 

  22. Yang SJ, Yin ZH, Ma XM, Lei FM. Phylogeography of ground tit (Pseudopodoces humilis) based on mtDNA: evidence of past fragmentation on the Tibetan Plateau. Mol Phylogenet Evol. 2006;41(2):257–65.

    Article  CAS  PubMed  Google Scholar 

  23. Qu Y, Zhao H, Han N, Zhou G, Song G, Gao B et al. Ground tit genome reveals avian adaptation to living at high altitudes in the Tibetan plateau. Nat Commun. 2013; doi:10.1038/ncomms3071

  24. Alatalo RV, Moreno J. Body size, interspecific interactions, and use of foraging sites in tits (Paridae). Ecology. 1987;68(6):1773–7.

    Article  Google Scholar 

  25. Gosler AG. Pattern and process in the bill morphology of the great tit Parus major. Ibis. 1987;129:451–76.

    Article  Google Scholar 

  26. Adams DC, Rohlf FJ, Slice DE. Geometric morphometrics: ten years of progress following the 'revolution'. Ital J Zool. 2004;71(1):5–16.

    Article  Google Scholar 

  27. Rohlf FJ, Marcus LF. A revolution in morphometrics. Trends Ecol Evol. 1993;8(4):129–32.

    Article  Google Scholar 

  28. Slice DE. Geometric morphometrics. Annu Rev Anthropol. 2007;36:261–81.

    Article  Google Scholar 

  29. Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Geometric morphometrics for biologists: a primer. London: Elsevier Academic Press; 2004.

    Google Scholar 

  30. Foster DJ, Podos J, Hendry AP. A geometric morphometric appraisal of beak shape in Darwin’s finches. J Evolution Biol. 2008; doi:10.1111/j.1420-9101.2007.01449.x.

  31. Degrange FJ, Picasso MBJ. Geometric morphometrics of the skull of Tinamidae (Aves, Palaeognathae). Zoology. 2010;113(6):334–8.

    Article  PubMed  Google Scholar 

  32. Klingenberg CP, Marugan-Lobon J. Evolutionary covariation in geometric morphometric data: analyzing integration, modularity, and allometry in a phylogenetic context. Syst Biol. 2013;62(4):591–610.

    Article  PubMed  Google Scholar 

  33. Barry RG. Mountain weather and climate. Cambridge, UK: Cambridge University Press; 2008.

    Book  Google Scholar 

  34. McCain CM, Grytnes J-A. Elevational gradients in species richness. In: Jonsson R, editor. Encyclopedia of life sciences (ELS). Chichester: Wiley; 2010. doi:10.1002/9780470015902.a0022548.

    Google Scholar 

  35. Price TD, Hooper DM, Buchanan CD, Johansson US, Tietze DT, Alström P et al. Niche filling slows the diversification of Himalayan songbirds. Nature. 2014; doi:10.1038/nature13272.

  36. Richman AD, Price T. Evolution of ecological differences in the Old World leaf warblers. Nature. 1992; doi:10.1038/355817a0

  37. Lack D. The ecology of closely related species with special reference to cormorant (Phalacrocorax carbo) and shag (P. aristotelis). J Anim Ecol. 1945;14(1):12–6.

    Article  Google Scholar 

  38. Campbell NA, Atchley WR. The geometry of canonical variate analysis. Syst Zool. 1981;30(3):268–80.

    Article  Google Scholar 

  39. Alatalo RV, Gustafsson L, Lundberg A. Interspecific competition and niche changes in tits (Parus Spp) - evaluation of nonexperimental data. Am Nat. 1986;127(6):819–34.

    Article  Google Scholar 

  40. Rohlf FJ. tpsDig. digitize landmarks and outlines, version 2.05. Stony Brook: Department of Ecology and Evolution, State University of New York; 2005.

    Google Scholar 

  41. Rohlf FJ. tpsRelw. relative warps analysis, version 1.36. Stony Brook: Department of Ecology and Evolution, State University of New York; 2003.

    Google Scholar 

  42. Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11(2):353–7.

    Article  PubMed  Google Scholar 

  43. Klingenberg CP, Barluenga M, Meyer A. Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry. Evolution. 2002;56(10):1909–20.

    Article  PubMed  Google Scholar 

  44. Rohlf FJ. Shape statistics: Procrustes superimpositions and tangent spaces. J Classif. 1999;16(2):197–223.

    Article  Google Scholar 

  45. Rohlf FJ, Slice D. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990;39(1):40–59.

    Article  Google Scholar 

  46. Rohlf FJ. Comparative methods for the analysis of continuous variables: geometric interpretations. Evolution. 2001;55(11):2143–60.

    Article  CAS  PubMed  Google Scholar 

  47. Johansson US, Ekman J, Bowie RCK, Halvarsson P, Ohlson JI, Price TD, et al. A complete multilocus species phylogeny of the tits and chickadees (Aves: Paridae). Mol Phylogenet Evol. 2013;69(3):852–60.

    Article  PubMed  Google Scholar 

  48. Packert M, Martens J. Taxonomic pitfalls in tits - comments on the Paridae chapter of the Handbook of the Birds of the World. Ibis. 2008;150(4):829–31.

    Article  Google Scholar 

  49. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–4.

    Article  CAS  PubMed  Google Scholar 

  50. Maddison WP. Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree. Syst Zool. 1991;40(3):304–14.

    Article  Google Scholar 

  51. Blomberg SP, Garland T, Ives AR. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57(4):717–45.

    Article  PubMed  Google Scholar 

  52. Pagel M. Inferring the historical patterns of biological evolution. Nature. 1999;401(6756):877–84.

    Article  CAS  PubMed  Google Scholar 

  53. Freckleton RP, Harvey PH, Pagel M. Phylogenetic analysis and comparative data: a test and review of evidence. Am Nat. 2002;160(6):712–26.

    Article  CAS  PubMed  Google Scholar 

  54. Mundry R. Statistical issues and assumptions of phylogenetic generalized least squares. In: Garamszegi LZ, editor. Modern phylogenetic comparative methods and their application in evolutionary biology. Berlin Heidelberg: Springer-Verlag; 2014.

    Google Scholar 

  55. Li GY, Zheng BL, Liu GZ. Fauna Sinica, Aves. Beijing: Science Press; 1982.

    Google Scholar 

  56. Mantel N. Detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27(2p1):209–20.

    CAS  PubMed  Google Scholar 

  57. Bird species distribution maps of the world. BirdLife International and NatureServe, Cambridge and Arlington. 2014. Accessed 5 May 2014.

  58. Hammer Ø, Harper DAT, Ryan PD. PAST: paleontological statistics software package for education and data analysis. Palaeontologia Electronica. 2001;4(1):9.

    Google Scholar 

  59. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2):289–90.

    Article  CAS  PubMed  Google Scholar 

  60. Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24(1):129–31.

    Article  CAS  PubMed  Google Scholar 

  61. Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.

    Article  Google Scholar 

  62. Pinheiro J, Bates D, DebRoy S, Sarkar D and R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1-122.2015.

  63. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB et al. vegan: community ecology package. R package version 2.3-3. 2016.

  64. Grant P, Abbott I, Schluter D, Curry R, Abbott L. Variation in the size and shape of Darwin's finches. Biol J Linn Soc. 1985;25(1):1–39.

    Article  Google Scholar 

  65. Kennedy JD, Weir JT, Hooper DM, Tietze DT, Martens J, Price TD. Ecological limits on diversification of the Himalayan core Corvoidea. Evolution. 2012;66(8):2599–613.

    Article  PubMed  Google Scholar 

  66. Price T. Morphology and ecology of breeding warblers along an altitudinal gradient in Kashmir, India. J Anim Ecol. 1991;60(2):643–64.

    Article  Google Scholar 

  67. Mayr E. Geographical character gradients and climatic adaption. Evolution. 1956;10(1):105–8.

    Article  Google Scholar 

  68. Ashton KG. Patterns of within-species body size variation of birds: strong evidence for Bergmann's rule. Global Ecol Biogeogr. 2002;11(6):505–23.

    Article  Google Scholar 

  69. Price T. Speciation in birds. Greenwood village. Colorado: Roberts and Company; 2008.

    Google Scholar 

  70. Barbosa A, Moreno E. Hindlimb morphology and locomotor performance in waders: an evolutionary approach. Biol J Linn Soc. 1999;67(3):313–30.

    Article  Google Scholar 

  71. Moreno E, Carrascal LM. Leg morphology and feeding postures in 4 Parus-species - an experimental ecomorphological approach. Ecology. 1993;74(7):2037–44.

    Article  Google Scholar 

  72. Dhondt AA. Effects of competition on great and blue tit reproduction: intensity and importance in relation to habitat quality. J Anim Ecol. 2010;79(1):257–65.

    Article  PubMed  Google Scholar 

  73. Sorensen MFL. Niche shifts of Coal Tits Parus ater in Denmark. J Avian Biol. 1997;28(1):68–72.

    Article  Google Scholar 

  74. Herrera CM. Niche-shift in genus Parus in southern Spain. Ibis. 1978;120(2):236–40.

    Article  Google Scholar 

Download references


We thank Dr. Per Alström for valuable comments on the manuscript, Dr. Xin-hai Li for help with the statistics and Dr. Christian Klingenberg and Dr. De-yan Ge for their useful advice on the geometric morphometric analyses. We are grateful to the National Zoological Museum, Institute of Zoology, Chinese Academy of Sciences for their generous support during the skin specimen examination. We sincerely thank the editor and reviewers for their valuable comments.


The project was funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB13020300), State Key Program of NSFC (31330073), a grant from the Ministry of Science and Technology of China (2014FY210200) to F.L., and NSFC program (J1210002) to S.S.

Availability of data and material

The datasets supporting the conclusions of this article are included within the article and the additional files.

Authors’ contributions

FL and SS conceived this study. SS, YQ, GS and FL designed the experiments. SS measured the specimens and analysed the morphological data. SS, QQ and TC conducted the statistical analysis and generated the figures. SS and FL wrote the manuscript. All authors read, revised, and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Not applicable. This research was conducted on museum skin specimens.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Fumin Lei.

Additional files

Additional file 1:

Supplementary materials. Included are 10 supplementary tables (Table S1-S10) and 2 supplementary figures (Figures S1-S2). (PDF 2118 kb)

Additional file 2:

Body measurements. Data for the linear measurements of all specimens involved in this study are shown, including body weight, body length, wing length, tail length, tarsus length and culmen. The altitude for each specimen is also included. (CSV 27 kb)

Additional file 3:

Coordinates of landmarks for geometric morphometric analysis of beak shape. (TXT 261 kb)

Additional file 4:

Phylogeny of the 14 Paridae species analysed in this study. (TRE 525 bytes)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shao, S., Quan, Q., Cai, T. et al. Evolution of body morphology and beak shape revealed by a morphometric analysis of 14 Paridae species. Front Zool 13, 30 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: