- Research article
- Open Access
Flightless birds are not neuroanatomical analogs of non-avian dinosaurs
BMC Evolutionary Biologyvolume 18, Article number: 190 (2018)
In comparative neurobiology, major transitions in behavior are thought to be associated with proportional size changes in brain regions. Bird-line theropod dinosaurs underwent a drastic locomotory shift from terrestrial to volant forms, accompanied by a suite of well-documented postcranial adaptations. To elucidate the potential impact of this locomotor shift on neuroanatomy, we first tested for a correlation between loss of flight in extant birds and whether the brain morphology of these birds resembles that of their flightless, non-avian dinosaurian ancestors. We constructed virtual endocasts of the braincase for 80 individuals of non-avian and avian theropods, including 25 flying and 19 flightless species of crown group birds. The endocasts were analyzed using a three-dimensional (3-D) geometric morphometric approach to assess changes in brain shape along the dinosaur-bird transition and secondary losses of flight in crown-group birds (Aves).
While non-avian dinosaurs and crown-group birds are clearly distinct in endocranial shape, volant and flightless birds overlap considerably in brain morphology. Phylogenetically informed analyses show that locomotory mode does not significantly account for neuroanatomical variation in crown-group birds. Linear discriminant analysis (LDA) also indicates poor predictive power of neuroanatomical shape for inferring locomotory mode. Given current sampling, Archaeopteryx, typically considered the oldest known bird, is inferred to be terrestrial based on its endocranial morphology.
The results demonstrate that loss of flight does not correlate with an appreciable amount of neuroanatomical changes across Aves, but rather is partially constrained due to phylogenetic inertia, evident from sister taxa having similarly shaped endocasts. Although the present study does not explicitly test whether endocranial changes along the dinosaur-bird transition are due to the acquisition of powered flight, the prominent relative expansion of the cerebrum, in areas associated with flight-related cognitive capacity, suggests that the acquisition of flight may have been an important initial driver of brain shape evolution in theropods.
Major behavioral transitions often correlate with neuroanatomical changes because novel sensory inputs and motor control pathways can form new or more robust connections, increasing the volume and density of associated regions of the brain [1,2,3,4,5]. One such evolutionary transition, from non-volancy to powered flight, has been acquired independently by three different vertebrate groups—pterosaurs, bats, and birds . Among these clades, birds provide arguably the best system within which to study the evolution of flight and associated morphological adaptations because they have a long stem lineage that is represented by a rich fossil record [7, 8]. This record has elucidated the origins of many of the postcranial ‘adaptations’ traditionally considered directly related to flight and has shown that these characters, such as the furcula and hollow long bones, arose not at the origin of this locomotory mode—near the appearance of Archaeopteryx lithographica—but earlier among non-avian theropods [9,10,11].
The early history of neuroanatomical adaptions related to this transition, however, is more obscure. Previous studies of avian neuroanatomical evolution indicated that the avian brain is characterized by an enlarged cerebrum and dorsoventrally flexed lateral profile. These studies divided endocasts into distinct volumes based on major brain regions (e.g. cerebrum, cerebellum, optic lobe). The fossil record shows that these ‘avian-like’ cerebral volumes evolved at least as far back as the origin of Maniraptora, ~ 160 Mya [11, 12]. This indicates that the size of the brain follows the more general trend of other ‘avian’ characteristics, evolving well before the origin of powered flight. Purely volumetric analyses, however, are limited in their capacity to characterize morphology. Neural pathways important in creating and regulating locomotor behavior often are distributed differentially among multiple regions of the brain , thus understanding the regional shape changes has the potential to inform the evolutionary tempo and mode of fight capacity along the theropod lineage.
Here, we use the loss of flight in crown group birds as a potential reverse analogue to the acquisition of flight in theropods. To determine if there is a link between neuroanatomy and loss of flight, we examined the endocasts of multiple independent pairs of flightless birds and their closest volant relatives to test for modifications in brain shape. We then compared the endocasts to the plesiomorphic condition in flightless, non-avian theropod dinosaurs. We used high-resolution X-ray computed tomographic (CT) imaging to construct three-dimensional (3-D) digital endocasts of the cranial cavity of a broad sample of modern flying and flightless members of Aves (sensu ) and extinct non-avian theropods. Shape data from these specimens were collected using a high-dimensional geometric morphometric (GM) approach [14,15,16] and they were subsequently subjected to a suite of multivariate analyses [17, 18]. With this approach, we evaluated whether (1) the loss of flight incurs predictable changes to neuroanatomical shape, such as a reduction in cerebral areas known to function in flight ; and (2) the endocast of modern flightless birds returns to a shape more similar to that of non-avian theropods than to their flying relatives.
CT imaging and digital endocast reconstruction
Taxonomic sampling included the skulls of 80 specimens representing 51 species: 25 flying and 19 flightless avians, and 7 non-avian dinosaurs: Alioramus altai  (IGM 100/1844), Khaan mckennai  (IGM 100/973), Citipati osmolskae  (IGM 100/978), Incisivosaurus gauthieri  (IVPP V 13326), Zanabazar junior  (IGM 100/1), an unnamed troodontid (IGM 100/1126 ) and Archaeopteryx lithographica  (NMNH PV OR 37001). We CT scanned each specimen using a General Electric v|tome Phoenix Computed Tomography (CT) scanner (General Electric, Heidelberg, Germany) at the Microscopy and Imaging Facility at the AMNH or the high-resolution source at the University of Texas at Austin High-Resolution X-ray Computed Tomography Facility (Additional file 1: Tables S1 and S2). Raphus cucullatus (NHMUK PV A9040)  was scanned at the NHMUK Imaging Facility. Scan data were imported into VGStudio MAX v2.2 (Volume Graphics GmbH, Heidelberg, Germany) to construct a 3-D model of the endocranial cavity (endocast) following the protocol set forth by Balanoff et al. . Because the brain fills nearly the entire cranial cavity in birds [25, 26], its morphology is accurately reflected by use of an endocast [24, 27]. The 3-D models were smoothed in VGStudio MAX and exported as PLY files for the software Landmark Editor v3.6 to virtually place landmarks .
Time-calibrated phylogenetic tree
A time-calibrated tree was used to perform comparative phylogenetic analyses on endocranial shape. For extant birds, we constructed a topology using the data provided online by Jetz et al.  (birdtree.org), pruned to include only species sampled in the current study. The Jetz et al.  trees were chosen for their dense species sampling. Our dataset has, in some cases, multiple species that are very closely related and the increased tip sampling of that tree allowed us to retain all of our species.
TreeAnnotator v1.8.1  was used to construct a maximum clade credibility tree from 1000 posterior trees based on the Hackett et al.  tree backbone . Because some of the sampled taxa are not included in the Jetz et al.  dataset, we replaced these tip labels with the closest relatives to each of these species in our analysis. Thus, Alle alle, Chunga burmeisteri, and Goura cristata, and Rollandia rolland were replaced by Pinguinus impennis, Pelecyornis australis, Raphus cucullatus, and Podilymbus gigas, respectively. We considered the above listed extinct crown group avians to exist in the Recent for the phylogeny (i.e., 0 Ma) in the context of deep time. The dataset included Struthio molybdophanes and S. camelus, however, only the latter was present in the Jetz et al.  dataset. Both species were included for the non-phylogenetically informed analyses (e.g. PCA and LDA), but only S. camelus was used for the construction of the time-calibrated tree and phylogenetically informed analyses. This tree was used for the crown-group Aves dataset (Fig. 1). For the Coelurosaur dataset, we used Mesquite v3.02  to incorporate non-avian dinosaurs into the time-calibrated avian phylogeny based on the phylogenetic relationships proposed by Brusatte et al. . Branch lengths were calculated for non-avian coelurosaurs based on oldest fossil ages from the Paleobiology Database (http://www.paleobiodb.org) in sister subclades for each internal node (e.g., the species range for Incisivosaurus is 122.46–125.45 Ma and for Oviraptorosauria is 130.0 Ma, based on the first fossil occurrence of Caudipteryx zoui; subtracting the first occurrence age from the latter gives 4.55 Myr, which is the branch length for Incisivosaurus). This combined extant and fossil tree was used for the Coelurosaur dataset (Fig. 1).
To collect landmark data, each 3-D endocranial reconstruction was imported into Landmark Editor  for digital placement of 3-D Cartesian coordinate points. The relatively featureless surface of an endocast presents an issue for landmark placement due to the lack of clear surficial unions of tissues or other large features that help define discrete landmarks . We placed evenly spaced 3-D semilandmarks on five sections of the endocast (i.e., left and right cerebrum, left and right optic lobes, and cerebellum). To do this, we first placed discrete Type I or Type II landmarks (sensu , also see ) along major endocranial features (e.g. the triple point between the optic lobe, cerebrum, and cerebellum). These landmarks were used to define patches to place sliding semilandmarks on each brain section. With this approach, each patch comprises discrete landmarks anchoring the patch boundary, 3-D semilandmarks on lines between discrete landmarks along the surface of the endocast model, and 3-D surface semilandmarks within the patch. Patches of 6 × 6 sliding semilandmarks were used for each half of the cerebrum (Fig. 2, Additional file 1: Tables S3 and S4). 4 × 4 patches were used for each optic lobe. A single patch of 4 × 5 semilandmarks was used for the cerebellum. In places where the edges of neighboring patches overlapped, we manually deleted the duplicated landmarks in R, leaving a total of 109 landmarks. We chose the number of semilandmarks for each patch so as to capture the shape of each lobe without oversampling it, which we confirmed using the function LaSEC in the LaMBDA R package . The raw landmark files are available on Dryad.
The landmark data were imported into R v3.1.1 , and aligned using “geomorph” v3.0.3  with generalized Procrustes superimposition and with sliding semilandmarks minimizing total bending energy [15, 16]. The aligned coordinates were imported into MorphoJ v1.06a  as two datasets: 1) full sample with non-avian and avian dinosaurs (N = 80; ‘Coelurosaur dataset’) and 2) subsampled dataset with the crown-group birds (N = 73; ‘Aves dataset’). An analysis using the Coelurosaur dataset indicated that landmark placement on Alioramus may have been incorrect, probably due to deformation and lack of easily identifiable landmarks. Therefore, a third dataset was created excluding Alioramus, which is discussed below as the ‘Coelurosaur dataset.’ Due to bilateral symmetry, the symmetrical component of the shape data was used for further analyses. Species with multiple specimens were aligned with the pooled data and the mean shape was calculated for each species.
We subjected each shape dataset to a principal components analysis (PCA) in MorphoJ. Plotting the scores associated with the first three PC axes creates a morphospace of the overall shape differences among specimens and identifies the morphological changes occurring along each axis. In addition to endocast images, ‘lollipop’ diagrams were used to visualize major shape changes occurring in morphospace along PC axes, where the vectors indicate the direction and magnitude of change from the mean shape. Using the “geomorph” R package , we also created a phylomorphospace based on the first two PC axes. The plot shows the correspondence of phylogenetic relatedness to morphological resemblance and can illustrate morphological innovation through the amount of morphospace explored by taxa . This analysis visualizes unequal magnitude of change per clade or branch and unequal morphological innovation by the geometry and relative length of the branches .
A dataset containing the PC axes encompassing 95% of the total shape variation was exported from MorphoJ and imported into R to run a linear discriminant analysis (LDA) using the “MASS” R package v7.3–45 . For the Coelurosaur dataset, we used the first 17 PC axes that encompass 95% of the total shape variation. This analysis looks to maximally separate the a priori locomotor groups (e.g. terrestrial, volant, secondarily flightless). If the LDA is able to adequately separate the shape data into a priori groups, then it suggests that locomotor mode induces certain neuroanatomical shapes across taxonomic groups. If not, then those groups are not supported. Group membership was cross-validated to see if the LD axes are able to correctly predict group assignment. Two locomotor categories (‘protoflying’ and ‘swimming’) had single members and therefore those groups were not included in the dataset for the analysis. The data point for these members were subsequently projected onto the morphospace based on LDA using the ‘predict’ function in R. We also conducted cross-validation analysis to assess the ability to correctly assign locomotory mode based on endocranial shape.
These analyses were run twice. The first analysis was run with crown-group avians (N = 73 specimens representing 44 species) and their locomotory mode (e.g. volant or flightless). 'Flightlessness' among extant birds was defined by a complete loss of lift generation, i.e., a total inability to create sufficient lift to raise the body off the ground for any amount of time, such as the ostrich (Struthio camelus), dodo (Raphus cucullatus), and kakapo (Strigops habroptilus) and was labeled ‘flightless’ in these analyses. A penguin, Eudyptes chrysocome, is present in the dataset. Penguins use subaqueous flight to propel through the water  and therefore it was placed in its own locomotion category (‘swimming’). As it was the only specimen for the ‘swimming’ locomotor category, it was removed from the dataset for LDA. Therefore, the first LDA was performed with two a priori groups: ‘volant’ (N = 37 specimens, representing 25 species), and ‘flightless’ (N = 35 specimens, representing 18 species), for a total of 72 specimens and 43 species.
The second LDA contained the original set of crown-group birds and six non-avian dinosaurs (three oviraptorosaurs, two troodontids, and Archaeopteryx). The non-avian dinosaurs (excluding Archaeopteryx) were categorized as ‘terrestrial’ because their primary locomotory mode was not inherited from flying ancestors. The distinction here between ‘secondarily flightless’ and ‘terrestrial’ is important because ‘secondarily flightless’ indicates that the species evolved from flight-capable ancestors, whereas ‘terrestrial’ indicates their ancestors were never capable of volant activity. Although a penguin (Eudyptes chrysocome) is present in the dataset, it was again removed for this analysis. Similarly, because of continuing debate concerning the level of flying ability of Archaeopteryx [12, 42, 43], it was defined as a ‘proto-flyer’ and was the only member of that locomotor category. As such, Archaeopteryx was also removed from the dataset for the LDA. Therefore, this LDA was performed with three a priori groups: ‘terrestrial’ (N = 5 specimens, each representing a single species), ‘volant’ (N = 37 specimens, representing 25 species), and ‘flightless’ (N = 35 specimens, representing 18 species), for a total of 77 specimens and 48 species. Once the analysis was complete, the ‘predict’ function was used to place Archaeopteryx and Eudyptes into the LD morphospace.
To test explicitly whether locomotory mode and other factors, such as allometry and phylogenetic inertia, drive predictable changes to endocranial shape, we performed several regression analyses. First, regressions of the symmetric component of the shape data against log-transformed centroid size were run for each dataset in MorphoJ using 10,000 replicates to assess the allometric effect. Second, phylogenetically-informed least-squares analyses on locomotory mode and shape data were performed using the ‘procD.pgls’ function in the “geomorph” R package. In addition, we evaluated the effect of phylogenetic signal using the ‘physignal’ function in the “geomorph” R package and allometry based on log-transformed centroid sizes exported from MorphoJ that were used as a proxy for endocranial size [44, 45].
Principal components analyses
The PCA of the Aves dataset generated 17 PC axes accounting for 95% of the symmetric component of shape variation, with the first three axes associated with 66.6% of the variation (36.4, 16.2, and 13.8%, respectively) (Additional file 1: Figure S1; overall distribution of data points resemble Coelurosaur dataset, Fig. 3). The shape of the posterior aspects of the endocast drive neuroanatomical changes along PC1, where positive scores indicate an anteroventral shift in the posterior margin of the cerebrum, a posteroventral shift in cerebellar location relative to the cerebrum, and a posterior shift and dorsoventral expansion of the optic lobes (Additional file 1: Figure S2B). PC2 correlates with lateral and anterior expansion of the anterior cerebrum, posteroventral reduction of the cerebellum, and a poteromedial shift of the optic lobes (Additional file 1: Figure S2C). Importantly, the plot of PC2 versus PC1 results in no clear visual separation between flying and flightless birds, but rather data points tend to cluster by clade (Fig. 3, Additional file 1: Figure S1). There is variable directionality in shape change within pairs of volant-flightless sister taxa. For example, there is a negative shift along PC1 in Columbiformes (from Caloenas nicobarica to Raphus cucullatus), but a positive shift along PC1 in Psittaciformes (from Nestor meridionalis to Strigops habroptilus), and a positive shift in PC2 in Cariamiformes (from Cariama cristata to Pelecyornis australis) (Additional file 1: Figure S1). In some groups, flightless members radiate in multiple directions away from volant members (e.g. Gruiformes and Anseriformes).
Within the Coelurosaur dataset the first 17 axes account for 95% of the symmetric component of shape variation, with the first three axes associated with 66.4% of the variation (38.7, 17.5, and 10.4%, respectively) (Fig. 3). The shape changes corresponding to PC1 and PC2 axes are equivalent to those in the Aves dataset. PC1 correlates with a relative anterodorsal shift of the anterior margin of the cerebrum, anteroventral shift of the anterior margin of the cerebellum and relative decrease in dorsoventral height of the cerebrum, and a posteroventral shift of the posterior margin of the cerebellum (Fig. 3 inset endocrania, Additional file 1: Figure S3B). Along PC2, more positive numbers indicate an anteroventral shift in the anterior cerebrum, a lateral reduction of the optic lobes, and a posterior shift and expansion of the cerebellum (Fig. 3 inset endocrania, Additional file 1: Figure S3C). The cumulative effect of these shifts is an anteroposteriorly longer and mediolaterally narrower endocast towards positive PC1 and a shorter, rounder endocast towards negative PC1. The PC morphospace of the Coelurosaur dataset shows no clear visual distinction between modern flying and flightless birds, however non-avian dinosaurs, as well as Archaeopteryx, occupy a unique region of morphospace (Fig. 3). Plots of PC3 versus PC1 show some separation between flying and flightless taxa in both the Aves (Additional file 1: Figure S4) and Coelurosaur datasets (Fig. 4), with non-avian dinosaurs again occupying their own morphospace. In both analyses, changes along PC3 are driven primarily by the lateral expansion of the cerebrum, a ventral shift of the anterior margin of the cerebrum, and a dorsolateral shift in the optic lobe (Fig. 4 inset endocrania, Additional file 1: Figure S5).
In phylomorphospace, the early history of the coelurosaurian lineage begins on the positive end of PC1 and shifts negatively (Fig. 3). Along PC1, Archaeopteryx occupies a position that lies between non-avian dinosaurs and crown-group birds (Fig. 3). Many internal nodes along the avian stem, as estimated by ancestral state reconstructions, appear in a small section of morphospace, which may be due to poor sampling at the base of Aves. Radiations from this small section of morphospace contain several subclades and their hypothetical ancestors that correspond to the origin of Neoaves. The test for phylogenetic signal indicates that phylogenetic signal is very low, but significant, in both the Aves dataset (Blomberg’s K = 0.04; P < 0.0001) and the Coelurosaur dataset (K = 0.05; P < 0.0001).
Linear discriminant analysis
Based on PC axes that account for 95% of symmetric component of shape, the morphospace constructed from LD axes and the Coelurosaur dataset shows crown group birds in a distinct cluster from non-avian dinosaurs and is able to visually distinguish between locomotory modes (Fig. 5). LD1 distinguishes non-avian dinosaurs and crown-group birds. However, volant and flightless crown-group birds are separated along LD2, orthogonal to major neuroanatomical changes along the dinosaur-bird transition. Shape changes towards negative LD2 are anteroposterior contraction of the cerebrum, posterolateral expansion of the cerebrum, narrowing of the cerebellum both in terms of lateral and anteroposterior extent, and an anterior and ventral shift of the optic lobes. Interestingly, the cross-validation technique based on LDA classifies Archaeopteryx as terrestrial with a posterior probability of ~1 (neither flying nor secondarily flightless, which had posterior probabilities of 3.89 e− 30 and 0.252e− 29, respectively). The Rockhopper penguin, Eudyptes chrysocome, falls among the volant birds with a posterior probability of 0.853, versus 0.147 for flightless and 2.66e− 16 for running. A cross-validation analysis indicated that some taxa were misclassified (Additional file 1: Table S5). Rallus philippensis (two of three specimens), Gallus gallus, Nestor meridionalis, Diomedea immutabilis, Phalacrocorax penicillatus (all three specimens), Coragyps atratus, Tachyeres patachonicus, Anas platyrhynchos, and Podilymbus podiceps are all volant birds that were classified as secondarily flightless. Phalacrocorax harrisi (one of two specimens), Rhea americana (one of three specimens), Pelecyornis australis, Gallirallus rovianae (two of three specimens), Raphus cucullatus, Tachyeres pteneres (one of three specimens), Strigops habroptilus, and Gallirallus australis (both specimens) are flightless birds that were classified as volant. Lastly, IGM 100/1126, a troodontid, which is terrestrial, was classified as flying. The other 53 specimens were correctly classified by this analysis. The LDA for the Aves dataset (Additional file 1: Figure S6), shows separation along the single LD axis between volant and flightless birds with small overlap in value between the two locomotory modes.
The lack of parallel neuroanatomical changes associated with the dinosaur-bird transition and volant-flightless birds implies that shifts between these locomotory modes incur contrasting changes to brain shape or these shifts have not contributed substantially to neuroanatomical variation. To explicitly assess the impact of locomotory mode on endocranial shape, we tested for the effect of locomotory mode on endocranial shape before and after correcting for phylogenetic structure and allometric signal. The ahistorical analyses indicate that size has a statistically significant effect in both the Aves dataset (R2 = 0.169; P < 0.0001) and the Coelurosaur dataset: (R2 = 0.106; P < 0.0001). However, when examined in a phylogenetic context, both size and locomotion independently and together failed to reject the null hypothesis of no effect on the shape data for the Aves dataset (P > 0.44; Table 1). The Coelurosaur dataset returned significant results for the effect of size (R2 = 0.397; P = 0.049) and locomotion (R2 = 0.007; P = 0.004) separately, but not the combined effects of size and locomotion (R2 = 0.295; P = 0.496; Table 1).
Our study presents several notable results regarding the relationship between flight capacity and brain shape evolution. These results, visually shown by the PC morphospaces, are not congruent with the hypothesis that flightless birds revert to a plesiomorphic brain shape that is similar to their terrestrial, non-avian relatives. In fact, there is comparatively little shape difference between flightless birds and their closest flying relatives (Fig. 3), indicating that flightless birds largely retain brain shapes more similar to their closest living relatives. Moreover, non-volant birds do not collectively approach the endocranial shape of non-avian theropods. Therefore, based on visual inspection of the PC morphospace, the neuroanatomy of flightless birds is clearly not a reverse analogue of flightless non-avian dinosaurs.
The LD morphospace, based on axes that maximize separation between group means, distinguishes between non-avian and avian groups, as expected from PCA. Importantly, it shows an overall separation between volant and flightless birds along LD2 with some overlap (Fig. 5). This shows that some aspects of endocranial shape can distinguish volant and flightless birds. The shape variation associated with LD2 implies that flightless birds, relative to their volant counterparts, generally exhibit anteroposterior expansion of the cerebrum, decreased dorsal convexity of the cerebrum, narrowing of the cerebellum, and inward contraction of the optic lobes.
Meanwhile, cross-validation analysis misclassifies the locomotory mode of approximately 30% of the specimens based on both the Aves and Coelurosaur dataset, indicating that neuroanatomical shape is not as reliable a neuroanatomical correlate for locomotor categories as the LD morphospace would make it seem. Among non-avian dinosaurs, IGM 100/1126 is misclassified as flying (Additional file 1: Table S5). IGM 100/1126 falls closer to the PC morphospace area of crown group birds than the other non-avian dinosaurs and Archaeopteryx (Fig. 3), even though it falls within the non-avian morphospace in the LDA (Fig. 5). However, the small number of non-avian coelurosaurs sampled for this study (N = 5) may simply have contributed greater degree of uncertainty in estimating locomotory mode. Sampling additional non-avian coelurosaurs through additional discoveries of well-preserved skulls may enhance the capacity of endocranial shape to distinguish between terrestrial and secondarily flightless modes. Importantly, none of the flightless birds were misclassified as terrestrial and based on Aves dataset, the locomotory mode of birds is misclassified for 25 of 72 taxa. This result corroborates the observation from PC morphospace that flightless and volant birds overlap in their neuroanatomical shape, but are distinct from non-avian dinosaurs. Taken together, the results from PC morphospace and LDA imply that although some neuroanatomical differences may be associated with secondary loss of flight, these changes account for very small proportion of shape variation to be a meaningful predictor for locomotory mode.
Furthermore, whereas ahistorical regression analyses on both the Aves and Coelurosaur datasets indicate a significant effect of allometric signal on brain shape, phylogenetically-informed regression analyses indicate that locomotory mode is not a significant predictor of endocranial shape with or excluding allometric signal for the Aves dataset (Table 1). The phylogenetically-informed analysis of variance of the Coelurosaur dataset returned significant results for the independent effect of locomotory mode and size on the shape of the endocast but insignificant results when tested together (Table 1). As such, there is no evidence that locomotory mode is a strong predictor of endocranial shape when corrected for size. Taken together, major shifts in locomotory mode have not driven uniform neuroanatomical changes within crown-group birds. Another possibility explaining the lack of locomotor signal is that the divergence of more recently evolved flying and flightless avian sister taxa occurred on a much shorter time-scale than the dinosaur-bird transition. Consequently, many flightless lineages may not have yet had the time to evolve clear neuroanatomical changes that would be associated with flightlessness. The latter scenario, however, is unlikely for paleognaths, which have been flightless since the early Cenozoic .
It is worth noting that the shape changes in this study only reflect surficial differences in neuroanatomy because we employed endocasts. As such, shifts in locomotory mode may have an impact on internal neuroanatomy that cannot be detected through the use of endocasts. Nevertheless, some surficial anatomy may reflect changes occurring internally. Although it cannot be confirmed with the currently available data, the separation of non-avian dinosaurs and crown-group birds in PC morphospaces indicate a clear trend towards relative expansion of the cerebrum along the dinosaur-bird transition (Figs. 3 and 4), which may be related to increased stimulus processing ability . A recent study suggests that the Wulst (a visual and somatosensory integration and processing center) and the entopallium (a visual interpretation area) are important in fast-paced visual processing during flight . Although the landmark scheme used here did not explicitly test for shape differences in the Wulst, the expansion of the posterodorsal region of the cerebrum in avian evolution (Fig. 3, Additional file 1: Figure S2), where these flight-related nuclei lie, may correspond to an increased trend in volant behaviors . Previous studies based on proportional volumes have shown that the cerebral expansion did not appear abruptly at the origin of Avialae, but rather, occurred more gradually throughout coelurosaurian evolutionary history [11, 12]. Inflated, ‘avian-like’ brain volumes first appear at the base of Maniraptora  and significant volumetric expansion of the cerebrum does not occur until well within the crown group . However, the shape data here demonstrate that the neuroanatomical shape differences between non-avian dinosaurs and crown-group birds are primarily driven by the cerebrum. Collectively, these results suggest that cerebral shape may have changed in response to the increased use of flight behaviors between non-avian dinosaurs and the modern radiation of birds even as proportional volume remained relatively constant . Intriguingly, LD2, that separates volant and flightless birds, is associated with slight reduction in the convexity of the cerebrum, including the Wulst, suggesting that its relative expansion could reflect flight capacity across coelurosaurs. Such investigation requires explicit characterization of Wulst morphology.
Although the overall predictive power of endocranial shape for locomotory mode is weak, the LDA based on current sampling classifies Archaeopteryx as terrestrial, and not volant. This is an intriguing result given the ongoing debate about its flight capabilities [42, 43]. The endocast of Archaeopteryx seems to more closely resemble its non-avian relatives rather than crown group birds, meaning it may not have been capable of extensive volant behaviors. Increased taxonomic sampling that fills in the gap between Archaeopteryx and the modern radiation of birds may strengthen our ability to classify the locomotory mode of the earliest members of Aves. The LDA classifies Eudyptes as a volant bird. Early studies note that the penguin uses subaqueous flight for locomoting underwater , and our study shows that its neuroanatomical shape corroborates the notion that subaqueous flight is equivalent to aerial flight. This outcome may be due to the fact that the evolution of wing-propelled diving (i.e. the acquisition of a new behavior) may have played a larger role in shaping the brain than the initial loss of aerial flight (i.e. the loss of an established behavior) .
If not locomotory mode, then what are the major drivers of brain shape evolution in birds? Previous GM studies on avian endocasts have shown that phylogenetic history, allometry, and orbit shape account for statistically significant (P < 0.001), but small proportions of endocranial shape variation [17, 49]. With different taxonomic sampling and use of high-dimensional shape characterization, our results corroborate these findings, indicating that there is allometric signal and a very weak, albeit statistically significant, phylogenetic signal present in the endocranial shape datasets. While these factors have variably contributed to evolutionary changes in bird brains, they collectively account for a small proportion (~ 10%) of the total variation [e.g., 49], suggesting that additional drivers are yet to be identified .
In addition, increased taxonomic sampling is still needed between the base of Avialae and the origin of living birds. Unfortunately, most of the fossils from this interval are fragmentary  or crushed  making endocranial construction extremely difficult. This paucity of specimens contributes to the apparent gap between Archaeopteryx and modern birds in the phylomorphospace, preventing a more reliable model of neuroanatomical evolution along this lineage. Early avialans undertook a variety of locomotor strategies, from the fully volant Enantiornithines  to the swimming Hesperornithiformes . Sampling this region of the tree, perhaps from recently discovered specimens , might bridge the gap from the non-volant theropods and at least partially volant Archaeopteryx to the fully volant crown group and provide examples of transitional brain morphologies not represented by living or extinct theropods.
The osteological changes associated with the evolution of flight in theropod dinosaurs is well documented in the fossil record, but the accompanying neuroanatomical changes are less known. Here, we used 3-D GM techniques on endocasts of extinct and modern dinosaurs (birds) to quantify and evaluate neuroanatomical changes related to this major locomotory change. In particular, we analyzed closely-related sets of volant and flightless birds to use as a potential reverse analogue to the initial acquisition of flight. The results demonstrate that loss of flight capacity does not incur predictive changes to endocranial shape in extant birds. In addition, the brain morphology of flightless birds do not converge to that of non-avian dinosaurs. While additional sampling may close the gap between non-avian dinosaurs and Neornithes, brain evolution from non-avian dinosaurs to birds appears to have been a unidirectional phenomenon or “point of no return.” Among crown-group birds, returning to an ancestral locomotory mode (i.e., flightlessness) does not correlate with a reversion to the more pleisiomorphic brain shape found in terrestrial non-avian dinosaurs. Modern flightless birds, therefore, are not reverse neuroanatomical analogues of non-avian dinosaurs.
- AMNH (FARB):
American Museum of Natural History, New York, NY, USA (Fossil Amphibian, Reptile, and Bird Collection)
Field Museum of Natural History, Chicago, IL, USA
Mongolian Institute of Geology, Ulaanbaatar, Mongolia
Institute of Vertebrate Paleontology and Paleoanthropology Beijing, People’s Republic of China
University of Kansas Natural History Museum, Lawrence, KS, USA
- NHMUK (PV OR, A):
Natural History Museum, London, United Kingdom (Paleontology Vertebrates and Old Register, Aves)
Smithsonian Institution National Museum of Natural History
Texas Cooperative Wildlife Collection, College Station, Texas, USA
University of Texas Memorial Museum, Austin, Texas, USA
Andrews TJ, Halpern SD, Purves D. Correlated size variations in human visual cortex, lateral geniculate nucleus, and optic tract. J Neurosci. 1997;17(8):2859–68.
Herculano-Houzal S, Manger PM, Kaas JH. Brain scaling in mammalian evolution as a consequence of concerted and mosaic changes in numbers of neurons and average neuronal cell size. Front Neuroanat. 2014;8(77):1–28.
Robinson CD, Patton MS, Andre BM, Johnson MA. Convergent evolution of brain morphology and communication modalities in lizards. Curr Zool. 2015;61(2):281–91.
Olkowicz S, Kocourek M, Lučan RK, Porteš M, Fitch WT, Herculano-Houzel S, Němec P. Birds have primate-like numbers of neurons in the forebrain. P Natl Acad Sci. 2016;113(26):7255–60.
Gold MEL, Schulz D, Budassi M, Gignac PM, Vaska P, Norell MA. Flying starlings, PET and the evolution of volant dinosaurs. Curr Biol. 2016;26:R265–7.
Smith JM. 1952. The importance of the nervous system in the evolution of animal flight. Evolution. 1952;6(1):127–9.
Xu X, Zhou Z, Dudley R, Mackem S, Chuong C-H, Erickson GM, Varricchio DJ. An integrative approach to understanding bird origins. Science. 2014;346(6215):1253293.
Brusatte SL, Lloyd GT, Wang SC, Norell MA. Gradual assembly of avian body plan culminated in rapid rates of evolution across the dinosaur-bird transition. Curr Biol. 2014;24(20):2386–92.
Nesbitt SJ, Turner AH, Spaulding M, Conrad JL, Norell MA. The theropod furcula. J Morphol. 2009;270(7):856–79. https://doi.org/10.1002/jmor.10724.
Prum RO, Brush AH. The evolutionary origin and diversification of feathers. Q Rev Biol. 2002;77(3):261–95. https://doi.org/10.1086/341993.
Balanoff AM, Bever GS, Rowe TB, Norell MA. Evolutionary origins of the avian brain. Nature. 2013;501(7465):93–6. https://doi.org/10.1038/nature12424.
Balanoff AM, Smaers JB, Turner AH. Brain modularity across the theropod–bird transition: testing the influence of flight on neuroanatomical variation. J Anat. 2015;229:204–14.
Turner AH, Makovicky P. Norell MA. A review of dromaeosaurid and paravian phylogeny. B Am Mus Nat Hist. 2012;371:1–209.
Adams D, Rohlf FJ, Slice DE. A field comes of age: geometric morphometrics in the 21st century. Hystrix. 2013;24(1):7–13.
Gunz P, Mitteroecker P, Bookstein FL. Semilandmarks in three dimensions. In: Slice DE, editor. Modern Morphometrics in physical anthropology. Developments in primatology: Progress and prospects. Springer, Boston, MA; 2005. p. 73–98.
Gunz P, Mitteroecker P. Semilandmarks: a method for quantifying curves and surfaces. Hystrix. 2013;24(1):103–9.
Kawabe S, Shimokawa T, Miki H, Matsuda S, Endo H. Variation in avian brain shape: relationship with size and orbital shape. J Anat. 2013;223(5):495–508.
Ahrens HE. Morphometric study of phylogenetic and ecologic signals in procyonid (Mammalia: Carnivora) endocasts. Anat Record. 2014;297(12):2318–30.
Brusatte SL, Carr TD, Norell MA. The osteology of Alioramus, a gracile and long-snouted tyrannosaurid (Dinosauria: Theropoda) from the late cretaceous of Mongolia. B Am Mus Nat Hist. 2016;366:1–197.
Clark JM, Norell MA, Barsbold R. Two new oviraptorids (Theropoda: Oviraptorosauria), upper cretaceous Djadokhta formation, Ukhaa Tolgod, Mongolia. J Vertebr Paleontol. 2001;21(2):209–13.
Balanoff AM, Xu X, Kobayashi Y, Matsufune Y, Norell MA. Cranial osteology of the theropod dinosaur Incisivosaurus gauthieri (Theropoda: Oviraptorosauria). Am Mus Novit. 2009;3651:1–35. https://doi.org/10.1206/644.1.
Norell MA, Makovicky PJ, Bever GS, Balanoff AM, Clark JM, Barsbold R, Rowe T. A review of the mongolian cretaceous dinosaur Saurornithoides (Troodontidae: Theropoda). Am Mus Novit. 2009;3654:1–63.
Gold MEL, Bourdon E, Norell MA. The first endocast of the extinct dodo (Raphus cucullatus) and an anatomical comparison amongst close relatives (Aves, Columbiformes). Zool J Linnean Soc. 2016;177(4):950–63. https://doi.org/10.1111/zoj.12388.
Balanoff AM, Bever GS, Colbert MW, Clarke JA, Field DJ, Gignac PM, Ksepka DT, Ridgely RC, Smith NA, Torres CR, Walsh S, Witmer LM. Best practices for digitally constructing endocranial casts: examples from birds and their dinosaurian relatives. J Anat. 2015;229:173–90. https://doi.org/10.1111/joa.12378.
Jerison HJ. Brain evolution and dinosaur brains. Am Nat. 1969;103(934):575–88.
Iwaniuk AN, Nelson JE. Can endocranial volume be used as an estimate of brain size in birds? Can J Zool. 2002;80(1):16–23.
Witmer LM, Ridgely RC, Dufeau DL, Semones MC. Using CT to peer into the past: 3D visualization of the brain and ear regions of birds, crocodiles, and nonavian dinosaurs. In: Anatomical imaging. Japan: Springer; 2008. p. 67–87.
Wiley DF, Amenta N, Alcantara DA, Ghosh D, Kil YJ, Delson E, Harcourt-Smith W, Rohlf FJ, St John K, Hamann B. Evolutionary morphing. In Visualization, 2005. VIS 05 IEEE (p. 431–438).
Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491(7424):444–8.
Rambaut A, Drummond AJ. 2013. TreeAnnotator v1. 7.0. Available as part of the BEAST package at http://beast.bio.ed.ac.uk.
Hackett SJ, Kimball RT, Reddy S, Bowie RC, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han KL, Harshman J, Huddleston CJ. A phylogenomic study of birds reveals their evolutionary history. Science. 2008;320(5884):1763–8.
Maddison WP, Maddison DR. 2015. Mesquite: a modular system for evolutionary analysis. Version 3.04 http://mesquiteproject.org
Bookstein FL. Morphometric tools for landmark data: geometry and biology. New York: Cambridge University Press; 1991.
Adams DC, Rohlf FJ, Slice DE. Geometric morphometrics: ten years of progress following the ‘revolution’. Ital J Zool. 2004;71(1):5–16.
Watanabe A. In review. How many landmarks are enough to characterize size and shape variation? PLoS One.
Core Team R. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016. https://www.R-project.org/
Adams DC, Otárola-Castillo E. Geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4(4):393–9.
Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11(2):353–7.
Sidlauskas B. Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008;62(12):3135–56.
Venables WN, Ripley BD. Modern applied statistics with S. 4th ed. New York: Springer; 2002.
Townsend CW. The use of the wings and feet by diving birds. Auk. 1909;26(3):234–48.
Foth C, Tischlinger H, Rauhut OW. New specimen of archaeopteryx provides insights into the evolution of pennaceous feathers. Nature. 2014;511(7507):79–82.
Feo TJ, Field DJ, Prum RO. Barb geometry of asymmetrical feathers reveals a transitional morphology in the evolution of avian flight. P R Soc B. 2015;282(1803):20142864.
Webster M, Sheets HD. A practical introduction to landmark-based geometric morphometrics. Quantitative methods in Paleobiology. 2010;16:168–88.
Zelditch ML, Swiderski DL, Sheets HD. Geometric morphometrics for biologists: a primer: Academic Press; 2012.
Ksepka DT, Phillips MJ. Avian diversification patterns across the K-Pg boundary: influence of calibrations, datasets. and model misspecification Ann Mo Bot Gard. 2015;100(4):300–28.
Carril J, Tambussi CP, Degrange FJ, Benitez Saldivar MJ, Picasso MBJ. Comparative brain morphology of Neotropical parrots (Aves, Psittaciformes) inferred from virtual 3D endocasts. J Anat. 2015;229:239–51.
Proffitt JV, Clarke JA, Scofield RP. Novel insights into the early neuroanatomical evolution in penguins from the oldest described penguin brain endocast. J Anat. 2016;229(2):228–38. https://doi.org/10.1111/joa.12447.
Marugán-Lobón J, Watanabe A, Kawabe S. 2016. Studying avian encephalization with geometric morphometrics. J Anat. 2016;229(2):191–203.
O’Connor JK, Li D-Q, Lamanna MC, Wang M, Harris JD, Atterholt J, You H-L. A new early cretaceous enantiornithine (Aves, Ornithothoraces) from northwestern China with elaborate tail ornamentation. J Vertebr Paleontol. 2016;36(1):e1054035. https://doi.org/10.1080/02724634.2015.1054035.
Ji S-H, Atterholt J, O’Conner JK, Lamanna MC, Harris JD, Lo D-Q, You H-L, Dodson P. A new, three-dimensionally preserved enantiornithine bird (Aves: Ornithothoraces) from Gansu Province, North-Western China. Zool J Linnean Soc. 2011;162:201–19.
Chiappe LM. Cretaceous birds of Latin America. Cretac Res. 1991;12(1):55–63.
Marsh OC. Odontornithes: a monograph on the extinct toothed birds of North America: with thirty-four plates and forty woodcuts, vol. 18: US Government Printing Office; 1880.
Xing L, O'Connor JK, McKellar RC, Chiappe LM, Tseng K, Li G, Bai M. A mid-cretaceous enantiornithine (Aves) hatchling preserved in Burmese amber with unusual plumage. Gondwana Res. 2017;49:264–77.
We thank A. Balanoff for sharing endocast data. We thank A. Balanoff and M. Norell for their comments on early drafts. We thank the Microscopy and Imaging Facility at the American Museum of Natural History, especially M. Hill and H. Towbin, and the University of Texas High Resolution CT Laboratory, especially J. Maisano and M. Colbert for imaging the specimens. E. Stanley and M. Colbert taught MELG the process of making endocasts. A. Perrard, J. Foox, M. Siddall, A. Laing, A. West, and K.J. Soda assisted with phylogenetic and statistical analyses. E. Wilberg provided helpful feedback on figures. We thank four anonymous reviewers for their constructive feedback.
MELG was supported by the Richard Gilder Graduate School, American Museum of Natural History, NSF DDIG (1311790). MELG and AW were supported by NSF DEB (1801224). The funding agencies had no role in the design of the study, in the collection, analysis, or interpretation of the data, nor in writing the manuscript.
Availability of data and materials
The landmark data sets and tree file supporting the results of this article are available in the Dryad repository, https://doi.org/10.5061/dryad.n1k111m.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. CT scanning parameters used for each specimen scanned for this project. Table S2. CT scanning parameters used for additional specimens used for this project from Balanoff et al. . Table S3. Discrete landmarks and anatomical descriptions. Table S4. Sliding semilandmarks used and anatomical descriptions. Table S5. Posterior probabilities from the cross-validation analysis of the LDA of the Coelurosaur dataset. Figure S1. Endocranial shape variation in the Aves dataset, illustrated in a morphospace constructed from PC1 and PC2. Figure S2. Shape changes along positive PC1 and PC2 for the Aves dataset shown as lollipop diagrams in dorsal and right lateral views. Figure S3. Shape changes along positive PC1 and PC2 for the Coelurosaur dataset shown as lollipop diagrams in dorsal and right lateral views. Figure S4. Endocranial shape variation in the Aves dataset, illustrated in a morphospace constructed from PC1 and PC3. Figure S5. Shape changes along positive PC3 for the Aves and Coelurosaur datasets shown as lollipop diagrams in dorsal and right lateral views. Figure S6. Bar plot of the single linear discriminant axis constructed from the LDA of endocranial shape variation in the Aves dataset. (DOCX 1520 kb)