Non-technical Summary
Theropods are bipedal dinosaurs that appeared 230 million years ago and are still extant as birds. Their history is characterized by extreme variations in body mass, with gigantism evolving independently and on several occasions between many theropod groups. However, no study has shown whether all theropods evolved the same limb adaptations to high body mass or whether they had different morphologies. Here we studied the shape variation across 68 femora from 41 species of theropods using a 3D comparative approach, multivariate statistics, and phylogenetically informed analyses. We demonstrated that all the heaviest theropods evolved similar adaptations regardless of their phylogenetic affinities by enlarging muscular attachments and articular surfaces. We also highlighted that the lightest theropods evolved femoral adaptations to miniaturization, which occurred close to the bird lineage (Avialae). In addition, our results support a gradual evolution of known “avian” features, independent from body mass variations, which may relate to a more “avian” type of locomotion, where the knee drives hindlimb movement instead of the hip, like in earlier theropod relatives. The distinction between body mass variations and a more “avian” locomotion is represented by a decoupling in the mediodistal crest morphology, whose biomechanical nature should be studied to better understand the importance of its functional role in gigantism, miniaturization, and the evolution of a more “avian” type of locomotion.
Introduction
Theropod dinosaurs first appeared in the Late Triassic and radiated in the Early Jurassic, before giving rise to birds (Ostrom Reference Ostrom1976; Gatesy and Middleton Reference Gatesy and Middleton1997; Padian and Chiappe Reference Padian and Chiappe1998; Holtz and Osmólska Reference Holtz, Osmólska, Weishampel, Dodson and Osmolska2004; Brusatte et al. Reference Brusatte, Nesbitt, Irmis, Butler, Benton and Norell2010; Irmis Reference Irmis2011; Benson Reference Benson2018), one of the most diverse clades of extant vertebrates. Non-avian theropods are characterized by a strictly bipedal locomotor habit and extreme body-size variations that ranged from a few grams to several tons, making them the largest terrestrial obligate bipeds ever (Sereno Reference Sereno1999; Carrano Reference Carrano2000, Reference Carrano, Carrano, Blob, Gaudin and Wible2006; Farlow et al. Reference Farlow, Gatesy, Holtz, Hutchinson and Robinson2000; Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014; Hutchinson Reference Hutchinson2021; D'Emic et al. Reference D'Emic, O'Connor, Sombathy, Cerda, Pascucci, Varricchio, Pol, Dave, Coria and Curry Rogers2023). Extreme theropod body-size increases appeared after the Triassic/Jurassic transition and followed a general trend toward gigantism that was common in several dinosaurian lineages (e.g., Theropoda, Sauropodomorpha, Thyreophora, Iguanodontia, Ceratopsidae) from the Early Jurassic to the Late Cretaceous, but also within Theropoda itself (Sereno Reference Sereno1999; Holtz and Osmólska Reference Holtz, Osmólska, Weishampel, Dodson and Osmolska2004; Carrano Reference Carrano, Carrano, Blob, Gaudin and Wible2006; Irmis Reference Irmis2011; Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014; Hutchinson Reference Hutchinson2021). Indeed, body mass increased convergently among several theropod lineages, reaching numerous occurrences of gigantism (Zanno and Makovicky Reference Zanno and Makovicky2013; Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014; Lee et al. Reference Lee, Cau, Naish and Dyke2014; Tsai et al. Reference Tsai, Middleton, Hutchinson and Holliday2018, Reference Tsai, Middleton, Hutchinson and Holliday2020; Cullen et al. Reference Cullen, Canale, Apesteguía, Smith, Hu and Makovicky2020); usually defined as masses greater than 1000 kg (Sereno Reference Sereno1999; Benson et al. Reference Benson, Hunt, Carrano and Campione2018; Cullen et al. Reference Cullen, Canale, Apesteguía, Smith, Hu and Makovicky2020; Hutchinson Reference Hutchinson2021). More specifically, Benson et al. (Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014) demonstrated that, out of 16 exceptional shifts (i.e., defined as shifts in evolutionary rates at nodes significantly differing from the overall pattern) in theropod body mass, 12 represented an increase. The remaining four occurrences of extreme shifts in body mass involved size decreases, with some linked to the miniaturization trend on the line to birds (Carrano Reference Carrano, Carrano, Blob, Gaudin and Wible2006; Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014, Reference Benson, Hunt, Carrano and Campione2018; Brusatte et al. Reference Brusatte, Lloyd, Wang and Norell2014; Lee et al. Reference Lee, Cau, Naish and Dyke2014). This miniaturization seems to have paralleled the evolution of a more knee-driven locomotion closer to that of crown birds (Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Hutchinson and Allen Reference Hutchinson and Allen2009; Allen et al. Reference Allen, Kilbourne and Hutchinson2021).
Inevitably, these changes of body mass are associated with changes in the morphology of hindlimb bones (Gregory Reference Gregory1912; Biewener Reference Biewener1983, Reference Biewener1989; Carrano Reference Carrano1998, Reference Carrano2001; Campione and Evans Reference Campione and Evans2012; Mallet et al. Reference Mallet, Cornette, Billet and Houssaye2019; Etienne et al. Reference Etienne, Filippo, Cornette and Houssaye2021). Archosaurian femoral morphology was highly constrained by several biological factors such as variations of body size/mass, shifts in locomotor habits, variation from a more sprawling to erect posture, and a shift from a more hip- to knee-driven gait, as demonstrated through morphometrics, biomechanics, microanatomy, and embryology (Charig Reference Charig1972; Bakker and Galton Reference Bakker and Galton1974; Parrish Reference Parrish1986; Gatesy Reference Gatesy1990; Christiansen Reference Christiansen1998; Carrano Reference Carrano1999; Farlow et al. Reference Farlow, Gatesy, Holtz, Hutchinson and Robinson2000; Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Maidment and Barrett Reference Maidment and Barrett2012; Bishop et al. Reference Bishop, Hocknull, Clemente, Hutchinson, Farke, Barrett and Lloyd2018a,Reference Bishop, Hocknull, Clemente, Hutchinson, Farke, Beck, Barrett and Lloydb; Allen et al. Reference Allen, Kilbourne and Hutchinson2021; Egawa et al. Reference Egawa, Griffin, Bishop, Pintore, Tsai, Botelho and Smith-Paredes2022; Lefebvre et al. Reference Lefebvre, Houssaye, Mallison, Cornette and Allain2022). The resulting shifts in the functional morphology and biomechanics of the femur have been investigated at the scale of Archosauria as a whole, but also between several dinosaurian clades and in the context of major faunal turnover across biological crises (Kubo and Kubo Reference Kubo and Kubo2012; Sookias et al. Reference Sookias, Butler and Benson2012; Bishop et al. Reference Bishop, Hocknull, Clemente, Hutchinson, Farke, Barrett and Lloyd2018a,Reference Bishop, Hocknull, Clemente, Hutchinson, Farke, Beck, Barrett and Lloydb, Reference Bishop, Bates, Allen, Henderson, Randau and Hutchinson2020; Cuff et al. Reference Cuff, Demuth, Michel, Otero, Pintore, Polet, Wiseman and Hutchinson2022; Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). However, how the femoral anatomy evolved along with repeated body mass increases among theropods remains poorly studied despite the high number of convergences. Furthermore, theropods provide a singular opportunity to study how femora of obligate bipeds evolved specialization to gigantism without shifting to quadrupedal habits, as in all other clades, which evolved gigantism primarily as quadrupeds.
Here, we investigate femoral changes that occurred within Theropoda during variations of body mass and the evolution of a more knee-driven locomotor habit toward the avian line from the Late Triassic to the Late Cretaceous. We quantify the morphological variation of theropod femora using 3D geometric morphometrics (3D GMM), which is well-suited to infer aspects of the functional morphology of hindlimb bones from their shape in extant and extinct vertebrates (Martín-Serra et al. Reference Martín-Serra, Figueirido and Palmqvist2014; Hedrick et al. Reference Hedrick, Dickson, Dumont and Pierce2020; Lefebvre et al. Reference Lefebvre, Allain, Houssaye and Cornette2020; Mallet et al. Reference Mallet, Billet, Cornette and Houssaye2022; Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). Following a similar approach among early archosauriform femora, Pintore et al. (Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b) demonstrated that the width of the epiphysis and the fourth trochanter position and shape were correlated to variations in body mass and independent from shifts in locomotor habits, as indicated by a decoupling in the fourth trochanter morphology between variations of body mass (i.e., more distally located in the heaviest taxa) and locomotor habits (i.e., proximo-distally symmetric in quadrupedal taxa and asymmetric with a steeper distal slope in bipedal taxa). In addition, features associated with body mass variations were demonstrated as evolving convergently between avemetatarsalians and pseudosuchians, whereas a high phylogenetic signal was recovered for features correlated to locomotor shifts. This finding sheds light on homoplastic features, which cladistic analyses desire to avoid. Therefore, we presume that features correlating with body mass variations in earlier archosauriforms would also vary convergently between the major theropod lineages included in our study, which independently evolved gigantism while remaining exclusively bipedal. Furthermore, the 3D GMM approach allows us to highlight several variations of commonly discussed femoral features relevant to theropod locomotion (e.g., proximo-distal shift of the fourth trochanter, offset of the femoral head, proximal extension of the lesser trochanter) and to test for their relationships with evolutionary allometry, phylogeny, and the evolutionary continuum from of more hip- to knee-driven type of locomotion within the lineage leading to extant birds. Investigating how femoral morphology accommodates convergent changes of body mass along with hindlimb orientation and more hip-/knee-driven locomotion in strictly bipedal animals is of major interest for understanding the extent of decoupling between morphological variations that are size related (homoplasy) versus those related to hindlimb orientation.
Materials and Methods
Sample
Our sample integrates 68 femora from 41 species of theropods, including 20 femora of specimens with an estimated body mass exceeding 1000 kg (which are hereby referred as “giants”) following the body mass estimation function for bipedal dinosaurs based on the minimal diaphyseal circumference (MDC) from Campione et al. (Reference Campione, Evans, Brown and Carrano2014; Table 1). The error range of body mass estimation was investigated by plotting a regression between log-transformed MDC and estimated body masses across all specimens from the sample (Supplementary Fig. S1). Specimens were selected to best document the multiple occurrences of gigantism across theropod phylogeny, within the constraint of availability and suitability for 3D surface digitization. We define our outgroup as all non-averostran theropods represented in our sample (Coelophysis bauri and Liliensternus liliensterni from the Late Triassic; Dilophosaurus wetherilli from the Early Jurassic; Fig. 1; Table 1). Ceratosauria is represented by several individuals from both Ceratosauridae and Abelisauroidea. Ceratosauridae includes two species of Ceratosaurus, the relatively small C. nasicornis and the more robust species C. dentisulcatus (body mass indicators are indicated in Table 1), which could be synonyms from different ontogenetic stages (Carrano and Sampson Reference Carrano and Sampson2008; Carrano et al. Reference Carrano, Benson and Sampson2012; Fig. 1). Our assemblage of Abelisauroidea is representative of the body-size range estimated within this clade by Grillo and Delcourt (Reference Grillo and Delcourt2017), because it includes several specimens of one of the smallest noasaurids, Masiakasaurus knopfleri, the large abelisaurid Carnotaurus sastrei, and a specimen from Argentina, currently under description, similar to Eoabelisaurus mefi and Skorpiovenator bustingorryi, which are also included in our study (Fig. 1; Table 1). Early Tetanurae are represented by a juvenile individual of the bizarre, herbivorous Chilesaurus diegosuarezi (Novas et al. Reference Novas, Salgado, Suárez, Agnolín, Ezcurra, Chimento, de la Cruz, Isasi, Vargas and Rubilar-Rogers2015), also hypothesized to represent the earliest-diverging ornithischian (Baron and Barrett Reference Baron and Barrett2017; Fig. 1; Table 1). Other tetanurans includes several specimens of Megalosauroidea, Allosauroidea, Tyrannosauroidea, Ornithomimidae, Therizinosauridae, Oviraptorosauridae, and paravians (Fig. 1; Table 1). Megalosauroidea are represented by Piatnitzkysaurus floresi from the Early Jurassic and Torvosaurus tanneri from the Late Jurassic (Fig. 1; Table 1). Megalosauroidea also includes two large Spinosauridae that each preserve a complete femur: Suchomimus tenerensis and Spinosaurus aegyptiacus, which may or may not have shared the same ecological (e.g., semi-aquatic) habits (Amiot et al. Reference Amiot, Buffetaut, Lécuyer, Wang, Boudad, Ding and Fourel2010; Ibrahim et al. Reference Ibrahim, Sereno, Dal Sasso, Maganuco, Fabbri, Martill, Zouhri, Myhrvold and Iurino2014, Reference Ibrahim, Maganuco, Dal Sasso, Fabbri, Auditore, Bindellini and Martill2020; Hone and Holtz Reference Hone and Holtz2021; Fabbri et al. Reference Fabbri, Navalón, Benson, Pol, O'Connor, Bhullar, Erickson and Norell2022; Sereno et al. Reference Sereno, Myhrvold, Henderson, Fish, Vidal, Baumgart, Keillor, Formoso and Conroy2022; Fig. 1; Table 1). Early Allosauroidea comprises Erectopus superbus, which may be an early metriacanthosaurid or an allosaurid (Allain Reference Allain and Carpenter2005; Carrano et al. Reference Carrano, Benson and Sampson2012; Fig. 1; Table 1). Later allosaurids comprise Allosaurus fragilis, A. jimmadseni, and an indeterminate specimen from the Late Jurassic of the Morrison Formation informally referred to as “Wyomingraptor,” which is relatively robust for the clade (Fig. 1; Table 1). Allosauridae's sister group Carcharodontosauridae includes the giant taxa Tyrannotitan chubutensis and Giganotosaurus carolinii (Fig. 1; Table 1). We include the proceratosaurid Guanlong wucaii as a small early-diverging tyrannosauroid (Fig. 1; Table 1). Tyrannosauridae comprises the albertosaurine Gorgosaurus libratus, the relatively small tyrannosaurine Alioramus altai, and its larger relatives Teratophoneus curriei, Daspletosaurus torosus, a juvenile Tarbosaurus bataar, and six specimens of Tyrannosaurus rex (Brusatte and Carr Reference Brusatte and Carr2016; Fig. 1; Table 1). Ornithomimosauria are represented by an assemblage of an undescribed early Ornithomimidae from the Early Cretaceous of Angeac-Charente (France; Allain et al. Reference Allain, Vullo, Le Lœuff and Tournepiche2014) and Gallimimus bullatus, a later relative of approximately the same size (Fig. 1; Table 1). Ornithomimosauria is the only major theropod clade that evolved gigantism (>1000 kg) for which we could not obtain representatives, such as deinocheirids. Therizinosauria comprises several specimens of Falcarius utahensis, an early representative of the clade, and one specimen of the larger therizinosaurid Nothronychus graffami (Fig. 1; Table 1). Oviraptosaurians are represented by both Caenagnathidae and Oviraptoridae (Fig. 1; Table 1). Oviraptorids include the heyuannian Conchoraptor gracilis and the relatively larger Citipatiinae Citipati osmolskae, represented here by the largest individual known for this species (Clark et al. Reference Clark, Norell and Barsbold2001). Caenagnathidae comprises an indeterminate species of Caenagnathus and its relatively larger relative Anzu wyliei. Paravians are represented by the small troodontids Sinovenator changii and an indeterminate species (Fig. 1; Table 1). Paravians also include Dromaeosauridae with a range of body sizes spanning from the small Bambiraptor feinbergi and Velociraptor mongoliensis to the middle-sized Deinonychus antirrhopus and several specimens of the relatively large Utahraptor ostrommaysi (Fig. 1; Table 1). When available, well-preserved left and right femora from the same specimens were used (3 out of 65 individuals).
Taphonomic alterations have been demonstrated to affect the results of geometric morphometric and statistical analyses (Hedrick et al. Reference Hedrick, Schachner, Rivera, Dodson and Pierce2019; Lefebvre et al. Reference Lefebvre, Allain, Houssaye and Cornette2020; Wynd et al. Reference Wynd, Uyeda and Nesbitt2021, Pintore et al. Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a). Following results from Pintore et al. (Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a), we included only complete femora without visible distortion of the anteroposterior and lateromedial curvature of the diaphysis, as these types of deformations were demonstrated to have the greatest impact on 3D GMM analyses conducted on asymmetrical objects such as femora. The whole femur of Nothronychus (giant Therizinosauridae) is taphonomically flattened anteroposteriorly but preserves the relative proportions of condyles and trochanters. This type of affine deformation was documented as only having a slight effect on spline deformations in 3D GMM (Pintore et al. Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a). Therefore, we chose to integrate it into our dataset to ensure representativeness of giant Therizinosauridae and to test for biological variations linked to an anteroposterior flattening of the femur in our results. Known ontogenetic stages are shown in Table 1.
Institutional Abbreviations
AMNH: American National History Museum, New York, NY, USA; ANG: Angeac-Charente Collection, Musée d'Angoulême, Angoulême, France; BYU: Brigham Young University Museum of Paleontology, Provo, UT, USA; CM: Carnegie Museum of Natural History, Pittsburgh, PA, USA; FIP: Florida Institute of Paleontology, Wellington, FL, USA; FMNH: Field Museum of Natural History, Chicago, IL, USA; FSAC: Faculté des Sciences Aïn Chock, Casablanca, Morocco; GPM: Glenrock Paleontological Museum, Glenrock, WY, USA; HMN: Museum für Naturkunde, Berlin, Germany; IGM: Mongolian Geological Institute, Ulaanbaatar, Mongolia; IVPP: Institute of Vertebrate Paleontology and Paleoanthropology, Beijing, China; MACN: Museo Argentino de Ciencias Naturales Bernardino Rivadavia, Buenos Aires, Argentina; MMCH: Museo Paleontológico Municipal Ernesto Bachmann, Villa El Chocón, Argentina; MNBH: Musée National Boubou-Hama, Niamey, Niger; MNHN: Muséum National d'Histoire Naturelle, Paris, France; MCZ VPRA: Vertebrate Paleontology (Reptiles and Amphibians) Collection, Museum of Comparative Zoology, Harvard University, Cambridge, MA, USA; MOR: Museum of the Rockies, Bozeman, MT, USA; MPEF: Museo Paleontológico Egidio Feruglio, Trelew, Argentina; MUCP: Museo de la Universidad Nacional del Comahue, Neuquén, Argentina; NCSM: North Carolina Museum of Natural Sciences, Raleigh, NC, USA; RAM: Raymond M. Alf Museum of Paleontology, Claremont, CA, USA; SNGM: Servicio Nacional de Geología y Minería, Santiago, Chile; TATE: Tate Museum, Casper College, Casper, WY, USA; TMM: Jackson School of Geosciences Vertebrate Paleontology Laboratory, University of Texas, Austin, TX, USA; TMP: Royal Tyrrell Museum of Palaeontology, Drumheller, Canada; TPI: North American Museum of Ancient Life, Lehi, UT, USA; UA: Université d'Antananarivo, Antananarivo, Madagascar; UCMP: University of California Museum of Paleontology, Berkeley, CA, USA; UMNH: Utah Museum of Natural History, Salt Lake City, UT, USA; ZPAL: Institute of Paleobiology of the Polish Academy of Sciences, Warsaw, Poland.
3D Digitization
Various digitization approaches were used to create the 3D models (Table 1). The surface scanners Next Engine (NextEngine Inc., Santa Monica, CA, USA), Artec EVA and Space Spider (Artec 3D, Luxembourg), and Polhemus FastSCAN (Polhemus, Inc., Colchester, VT, USA) were used with their respective software for surface reconstructions (19 out of 68 femora): ScanStudio Pro (NextEngine Inc., Santa Monica, CA, USA) and Artec Studio Professional (Artec 3D, Luxembourg). Photogrammetry was performed using a Sony DSC-F828 (29 out of 68). Surface reconstructions were conducted using the package Bundler,Footnote * PMVS (Patch-based Multi-view Stereo Software),Footnote ** Metashape Professional (Agisoft LLC, Russia), and Meshlab (Cignoni et al. Reference Cignoni, Callieri, Corsini, Dellepiane, Ganovelli and Ranzuglia2008) to compute dense point clouds, create scaled 3D meshes, and mirror all the left femora. CT scans were segmented using Mimics (Materialise NV, Belgium) and Avizo (Thermo Fisher Scientific Inc., MA, USA) software to create 3D meshes (20 out of 68; see scan parameters in Supplementary Table S1). Several studies have demonstrated that the difference in the geometry of the resulting 3D morphologies was low when comparing photogrammetry, CT scans, and surfaces scans, especially when analyzing specimens at a large comparative scale (Falkingham Reference Falkingham2012; Fau et al. Reference Fau, Cornette and Houssaye2016; Soodmand et al. Reference Soodmand, Kluess, Varady, Cichon, Schwarze, Gehweiler, Niemeyer, Pahr and Woiczinski2018; Díez Díaz et al. Reference Díez Díaz, Mallison, Asbach, Schwarz and Blanco2021). Furthermore, Waltenberger et al. (Reference Waltenberger, Rebay-Salisbury and Mitteroecker2021) demonstrated that 3D GMM could properly analyze combined 3D models digitized from various approaches.
Geometric Morphometrics
We used 3D GMM to investigate the morphological variation within our dataset. This method relies on the digitization of anatomically (or geometrically) homologous landmarks (Zelditch et al. Reference Zelditch, Swiderski and Sheets2012). We digitized single anatomical landmarks and sliding semilandmarks along curves and surfaces following the protocol described in Gunz et al. (Reference Gunz, Mitteroecker, Bookstein and Slice2005), Gunz and Mitteroecker (Reference Gunz and Mitteroecker2013), and Botton-Divet et al. (Reference Botton-Divet, Cornette, Fabre, Herrel and Houssaye2016). High-density 3D GMM is more effective than using anatomical landmarks alone to capture the entire morphology of biological objects and precisely describe shape variations, which is particularly useful, for example, to circumvent the lack of anatomical landmarks along the diaphysis of long bones (Gunz et al. Reference Gunz, Mitteroecker, Neubauer, Weber and Bookstein2009; Zelditch et al. Reference Zelditch, Swiderski and Sheets2012; Gunz and Mitteroecker Reference Gunz and Mitteroecker2013; Botton-Divet et al. Reference Botton-Divet, Houssaye, Herrel, Fabre and Cornette2015; Goswami et al. Reference Goswami, Watanabe, Felice, Bardua, Fabre and Polly2019).
In total, 662 landmarks, including 27 anatomical single landmarks, 163 sliding semilandmarks along curves, and 472 along surfaces, were used to measure the femoral variation across our dataset. We manually digitized anatomical landmarks and sliding semilandmarks along curves only on every specimen using the Landmark software v. 3.0.0.6 (Wiley et al. Reference Wiley, Amenta, Alcantara, Ghosh, Kil, Delson, Harcourt-Smith, Rohlf, St John and Hamann2005). To minimize the impact of taphonomically eroded structures, we relied on concavities rather than convexities (e.g., borders rather than the tip of tuberosities; Supplementary Table S2). Unfortunately, we could not digitize landmarks along all key structures relevant to theropod locomotor functional morphology and biomechanics, because their borders were not clearly visible throughout the entirety of the specimens included in this study. For example, the trochanteric crest, on which the iliofemoralis externus (IFE) is thought to have inserted, was strongly reduced in Tetanurae (Hutchinson Reference Hutchinson2001) and therefore could not be included in our 3D GMM protocol. Additionally, we digitized the entirety of landmarks and semilandmarks (surfaces included) on one specimen hereby referred as “the template.” We chose the femur of Tyrannosaurus rex FMNH PR2081, as it represents one of the best preserved specimens (i.e., with the most prominent structures of our dataset), ensuring that sliding semilandmarks will be correctly projected along the surfaces of every other specimen (Fig. 2). The projection was performed semi-automatically using the function placePatch of the R package Morpho v. 2.8 (Schlager Reference Schlager, Zheng, Li and Szekely2017). Then, we performed two sets of spline relaxations minimizing the bending energy of a thin plate spline: (1) five iterations between the configurations of the template and every other specimen using the function relaxLM of Morpho; and (2) three iterations between a mean configuration and every specimen using the function slideLM. A more detailed workflow is presented in Pintore et al. (Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). These steps ensured that all corresponding semilandmark coordinates were geometrically homologous between each specimen (Gunz et al. Reference Gunz, Mitteroecker, Bookstein and Slice2005). We performed a generalized Procrustes analysis (GPA) using the gpagen function from geomorph v. 3.3.1 (Adams and Otárola-Castillo Reference Adams and Otárola-Castillo2013). GPA enables superimposition of every specimen in order to homogenize their positions in the Cartesian coordinate system and isolate the shape from the size component (Gower Reference Gower1975; Rohlf and Slice Reference Rohlf and Slice1990; Zelditch et al. Reference Zelditch, Swiderski and Sheets2012). The remaining differences between superimposed configurations (i.e., Procrustes residuals) could then be interpreted as the shape variation between all specimens.
We computed a principal component analysis (PCA) to reduce dimensionalities of the variation and isolate several morphological components (Gunz and Mitteroecker Reference Gunz and Mitteroecker2013). We also computed neighbor-joining trees using the total PCA scores to represent global morphological distances between each bone. Trees were computed in unrooted type based on Euclidean distances using the distance, nj, and plot.phylo functions of the ape R package (Paradis et al. Reference Paradis, Claude and Strimmer2004). Then, we computed regressions between PC axes and log-transformed estimated body masses to investigate whether one of these morphological components would highlight shape variations linked with body mass increase (Table 1). We performed repeatability testing by digitizing the same landmark configuration 10 times on three Tyrannosaurus femora (i.e., FMNH PR2081, TMM 40811, MOR 1125), which was the taxon with some of the least intraspecific variability between specimens from our dataset (Table 1). After computing a GPA and a PCA, we highlighted three distinct clusters corresponding to the three different specimens, each composed of the 10 iterations of the same landmark configurations (Supplementary Fig. S2). This demonstrated that our landmark configuration was repeatable, because the biological variation was greater than the possible operator bias in reproducing the same landmark configuration multiple times on the same specimen. Once we identified the PC axis that polarized variations related to body mass increases, we computed 3D visualizations in order to highlight the associated femoral features. To do so, we calculated a mean shape between all the femora of our dataset by performing a spline relaxation between the template landmark configuration and a mean landmark configuration extracted from the GPA. The associated TPS deformation was used to warp the template mesh onto the mean landmark configuration in order to create the mean shape of all the specimens. From this step, we could interpolate the mean landmark configuration with the configurations at the most negative and positive extremes of the chosen PC axis in order to create minimal and maximal theoretical shapes. To further highlight which features varied the most, we computed vectors of displacement between every landmark and semilandmark of the two extreme morphologies using the function segments3d of the rgl R package v. 0.100.54 (Adler and Murdoch Reference Adler and Murdoch2020) that we colored with a heat-map gradient using the blue2red functions of the ColorRamps R package (Keitt Reference Keitt2008; Botton-Divet Reference Botton-Divet2017).
The allometric effect (i.e., evolutionary and ontogenetic allometry; Klingenberg Reference Klingenberg2016) was analyzed using different sets of metrics (Table 1). We measured the femoral lengths (FL; i.e., maximal height of the proximal end relative to the distal end) across all our specimens using Meshlab v. 2020.06 (Cignoni et al. Reference Cignoni, Callieri, Corsini, Dellepiane, Ganovelli and Ranzuglia2008). Campione and Evans (Reference Campione and Evans2012, Reference Campione and Evans2020) demonstrated that the MDC of the femur was a reliable indicator of body mass (BM) in dinosaurs. We thus extracted MDC using the Cross Section and Extract Contours tools of the software CloudCompare v. 2.12.Footnote *** We computed more precise estimations of BM following the methodology provided by (Campione et al. Reference Campione, Evans, Brown and Carrano2014), which was developed specifically for the femur of bipedal non-avian dinosaurs, using the function bipeds from the R package MASSTIMATE v. 2.0 (Campione Reference Campione2016). The equation relies on MDC, and the output is corrected by an index (“cQE.cor” in the bipeds function) based on the eccentricity at the MDC of the femur, which we measured using Meshlab v. 2020.06 (Cignoni et al. Reference Cignoni, Callieri, Corsini, Dellepiane, Ganovelli and Ranzuglia2008). We also investigated the size effect within our dataset by computing Pearson's correlation tests between the log-transformed centroid size and estimated body masses of each specimen and their distribution along the chosen PC axis using the R function cor.test, for which a significant result would indicate that shape variation along that axis had an allometric component (Mitteroecker and Gunz Reference Mitteroecker and Gunz2009).
We constructed a phylogeny following Carrano and Sampson (Reference Carrano and Sampson2008), Carrano et al. (Reference Carrano, Benson and Sampson2012), Brusatte and Carr (Reference Brusatte and Carr2016), Zanno and Makovicky (Reference Zanno and Makovicky2013), Turner et al. (Reference Turner, Makovicky and Norell2012), Novas et al. (Reference Novas, Salgado, Suárez, Agnolín, Ezcurra, Chimento, de la Cruz, Isasi, Vargas and Rubilar-Rogers2015), and Funston (Reference Funston2020) using Mesquite software v. 3.6.1 (Maddison and Maddison Reference Maddison and Maddison2019) with all branch lengths set to one (Fig. 1). We mapped this phylogeny onto the PCA in order to compute a phylomorphospace using the function plot.gm.prcomp of the geomorph package. We used the K mult statistics, an extension of Blomberg's K (Blomberg et al. Reference Blomberg, Garland and Ives2003), to quantify how much the constructed phylogeny compares with its expectation under a Brownian motion model based on the distribution of specimens in the complete morphospace and along chosen PC axes (Adams Reference Adams2014a). When significant, the K mult statistics indicates a relationship between phylogeny and the variation of femoral morphology. A value below one would indicate that the phylogeny did not predict the variation of femoral morphology under a Brownian model of evolution well enough, likely suggesting convergence, whereas a value equal to or above one would suggest variation between clades and according to the phylogeny (i.e., phylogenetic signal). In addition, we conducted phylogenetic generalized least squares regressions (PGLS) to investigate the covariation between our different size metrics and the femoral morphology while accounting for phylogenetic relationships under Brownian motion across our theropod dataset at the unidimensional (i.e., along one isolated PC axis) and multidimensional (i.e., across all PC axes) levels (Adams Reference Adams2014b).
To assess branch length sensitivity over our phylogenetically based statistical analyses (i.e., phylogenetic signal and PGLS), we compared results obtained from a time-calibrated tree (with real and transformed branch lengths) with our results obtained from our phylogeny with equal branch lengths (Supplementary Material). The time-calibrated tree was created based on the minimal and maximal known ages from each taxon in our sample. We used the R package paleotree (Bapst Reference Bapst2014) to optimize date ranges per taxon in branch lengths and to estimate node dates from these time ranges using the dateNodes function, which provided the time-calibrated tree with “real” branch lengths. We then reran the analyses on the time-calibrated tree following the same protocol as described earlier. The sensitivity analyses yielded similar results when investigating phylogenetic signal and evolutionary allometry using equal and real branch lengths (Supplementary Material). Therefore, we favored using equal branch lengths over “real” branch lengths, as the latter represented a greater assumption, given that our sample is phylogenetically large and composed of fossils with uncertainty regarding their ages (Supplementary Material).
Results
Morphological Variation
First, we investigated the global morphological variation using the neighbor-joining tree (Supplementary Fig. S3). This visualization highlights that specimens from the same species are consistently clustered together but that clades seem spread across the neighbor-joining tree (Supplementary Fig. S3). This observation is consistent with the distribution of specimens along the first axis of the PCA (PC 1; Fig. 3A). The first two axes of the PCA (PC 1 and PC 2) represent more than 50% of the global morphological variation (40.0% and 11.3%, respectively; Fig. 3A). PC 1 represents essentially femoral robusticity (i.e., epiphyseal width relative to the diaphyseal length; Figs. 3, 4A), whereas PC 2 broadly represents the degree of fusion of trochanters and curvature of various femoral features (Figs. 3, 4B).
The distribution of the specimens along the first axis follows an increase of body mass, as shown by the increase of the circle size toward the positive side (Fig. 3A) and its significant association with log-transformed centroid size (see details in “Evolutionary Allometry”). Indeed, every specimen with a body mass estimation > 1000 kg is located on the positive side of PC 1, except the large abelisaurid Carnotaurus (Fig. 3A; Car). All early non-averostran neotheropods and maniraptoriforms (i.e., pennaraptoran and non-pennaraptoran maniraptoriforms) are located on the negative side of PC 1 (except the giant therizinosaurid Nothronychus, the relatively large dromaeosaurid Utahraptor, and, surprisingly, the smallest of our oviraptosaurids, Conchoraptor; Figs. 3A, 5; Not, Con). All troodontids and the ~2 kg dromeosaurid Bambiraptor form a cluster together on the most negative extreme of PC 1 and are clearly isolated from the rest of the sample (Figs. 3A, 5; Bam, Sin, Tro). All specimens of Tetanurae, other than maniraptoriforms (i.e., Chilesaurus; allosauroids; megalosauroids, except the early Erectopus; and tyrannosauroids, except the smallest taxa Guanlong and Alioramus), are located on the positive side of PC 1 (Fig. 3A; Chi, Ere, Gua, Ali). Ceratosaurs are located close to the center of PC 1, with the smallest (Masiakasaurus) and the largest (Carnotaurus) taxa having negative values and the other, medium-sized taxa, having positive values (Fig. 3A; Mas, Car).
Theoretical shapes and vectors displaying the amount of landmark displacement highlight that the most important morphological variation is located on the distal ridge of the fourth trochanter (Fig. 4A). The ridge is also more prominent and proximo-distally wider on the maximal theoretical shape, whereas it is flatter and proximal-distally narrower on the minimal theoretical shape (Figs. 3B, 4A). The distribution along PC 1 also varies according to the epiphyseal width relative to the diaphyseal length (i.e., femoral robusticity; Fig. 3B,D,E). In the proximal epiphysis, the longest vectors are on the medial side of the femoral head, which is upturned on robust femora and downturned on gracile ones (Fig. 4A). In the distal epiphysis, the longest vectors are on the anterior part of the mediodistal crest, which is greatly proximo-distally enlarged toward the middle of the shaft in the most robust femora (Fig. 4A). Finally, the anterior condyles are more prominent on the maximal theoretical shape than on the minimal one (Figs. 3E, 4A).
The distribution along the second axis does not follow variation in body mass (see “Evolutionary Allometry” for details), and specimens with a body mass estimated to more than 1000 kg are located on both sides of the axis, as shown by the circle size (Fig. 3A). However, the distribution along PC 2 highlights a relatively clear separation between pennaraptorans (i.e., the clade closest to birds in our dataset) and other major neotheropod clades along this axis (Fig. 5). Early non-averostran neotheropods from the Late Triassic and Early Jurassic (except one Coelophysis specimen out of five) and all ceratosaurians are located on the negative part of the morphospace, whereas all pennaraptorans are located on the positive part of the morphospace (Figs. 3A, 5; Coe). The rest of the dataset (i.e., non-pennaraptoran Tetanurae) is spread across the positive and negative parts but mostly closer to the middle part of the morphospace (Figs. 3A, 5). Among the non-pennaraptoran tetanurans, the earliest representatives of each clade are located closer to the negative side of the axis than the more derived relatives from the same clade, as seen among megalosauroids (i.e., Piatnitzkysaurus located more negatively than other megalosauroids or more negatively than Torvosaurus and Spinosauridae), tyrannosauroids (i.e., Guanlong) and ornithomimosaurians (i.e., ornithomimid indet.; Figs. 3A, 5; Supplementary Fig. S3; Pia, Gua, Ali, Orn, Gal). We made the opposite observations for allosauroids (i.e., Erectopus) and therizinosaurs (i.e., Falcarius; Fig. 3A; Ere, Fal). Among ceratosaurians, all abelisauroids are closer to the negative side than ceratosaurids (Figs. 3A, 5).
Theoretical shapes along PC 2 highlight that the longest displacement vectors are localized along the proximal part of the mediodistal crest, which is more convex and prominent toward the medial side in the minimal theoretical shape, making it look more sigmoidal (i.e., more curved mediolaterally) than in the maximal theoretical shape (Figs. 3B,F, and 4B). In addition, some of the longest vectors are localized on the proximal area of the lesser trochanter, which is more prominent anteriorly and has a proximal extent equal to that of the greater trochanter on the maximal theoretical shape than on the minimal one (Figs. 3G, 4B). Similarly, the anteroposterior width of the greater trochanter is wider on the maximal shape than on the minimal one, which results in a narrower gap between the greater and lesser trochanter, indicating the fusion of these two tubercles in the taxa on the most positive side of PC 2 (Figs. 3G,H, and 4B). This mostly represents a fusion between the lesser and the greater trochanters toward the pennaraptoran clade (Fig. 5). Finally, there is a greater offset between the epiphyses in the minimal theoretical shape than in the maximal one (Figs. 3H,I, and 4B).
Phylogenetic Signal
Results from the multivariate K statistics highlighted a significant association between femoral morphology and phylogeny when performed on the global variation (K mult) and the two first PC axes individually (K). The phylogenetic signal ratio was clearly below one for the global morphological variation (K mult = 0.39, p-value < 0.01), and for the first two PC axes (PC 1: K = 0.77, p-value < 0.01; PC 2: K = 0.49, p-value < 0.01). These results demonstrated that the morphological variation had a lower phylogenetic signal than predicted by the Brownian motion, indicating that it was also structured within clades (i.e., homoplastic) (Fig. 5).
The phylomorphospace highlights this within-clade variation along PC 1, except for five (out of 15) pennaraptorans (all troodontids and the smallest dromaeosaurid Bambiraptor) that are clearly isolated on the negative side from the rest of the specimens (Fig. 5; Bam, Sin, Tro). All other pennaraptorans share similar levels of femoral robusticity with the smallest representatives of the other tetanurans (Chilesaurus, Piatnitzkysaurus, Erectopus, Guanlong, Alioramus, Falcarius, ornithomimid indet.), but also with all non-averostran neotheropods and all ceratosaurians, despite early theropods, ceratosaurs, and maniraptorans being the clades most phylogenetically distant from each other among our dataset (Fig. 5; Chi, Pia, Ere, Gua, Ali, Fal, Orn). The largest representatives of several clades of tetanurans all group together on the most positive side of PC 1, regardless of phylogenetic distinction, as shown by the distance between the two carcharodontosaurids Tyrannotitan and Giganotosaurus (Fig. 5).
Along PC 2, the earliest representatives of each theropod lineage are located on the opposite side (whether negative or positive) of their later relatives (i.e., Erectopus on the most positive part of the allosauroids’ morphospace; Guanlong on the most negative part of the tyrannosauroids’ morphospace; Piatnitzkysaurus on the most negative side of the megalosauroids’ morphospace; Fig. 5). However, despite the low phylogenetic signal indicated by the K statistics, we observed that the earliest-diverging taxa (non-averostran neotheropods) are mostly grouped on the most negative side of PC 2, whereas later taxa such as non-avian pennaraptorans that are the closest relatives to the avian lineage group on the most positive side of PC 2 (Fig. 5). All other tetanurans group in the middle part of the phylomorphospace (Fig. 5).
Evolutionary Allometry
We found a strong and positive association when performing Pearson's correlation test between log-transformed centroid size and the distribution along PC 1 (Fig. 6), which mostly represented femoral robusticity (r 2 = 0.6, p-value < 0.01), but no significant correlation along PC 2 (p-value > 0.05). We demonstrated an even stronger association when testing for the association between log-transformed estimated body mass (based on MDC) and PC 1 (r 2 = 0.7, p-value < 0.01). We obtained similar results when factoring out the phylogenetic resemblance between investigated taxa by computing PGLS at the unidimensional (i.e., isolated PC axis) level on log-transformed centroid size (PC 1: r 2 = 0.2, p-value <0.01; PC 2: p-value > 0.05) and log-transformed estimated body mass (PC 1: r 2 = 0.4, p-value < 0.01; PC 2: p-value > 0.05). In addition, the PGLS at the multidimensional level (i.e., all PC axes; the complete morphological variation) revealed a significant but low association of the global morphological variation with both log-transformed centroid size (r 2 = 0.07, p-value < 0.01) and estimated body mass (r 2 = 0.09, p-value < 0.01).
Discussion
Convergent Specialization to Gigantism in Theropods and Other Archosauriforms Relative to the Evolution of Bipedalism
Femoral Specializations to Gigantism in Theropods
A positive relationship between femoral robusticity (i.e., increase of width relative to length) and body mass is well known to occur among quadrupedal mammals and archosaurs (Gregory Reference Gregory1912; Biewener Reference Biewener1983, Reference Biewener1989; Carrano Reference Carrano1998, Reference Carrano2001; Christiansen Reference Christiansen1999; Campione and Evans Reference Campione and Evans2012; McPhee et al. Reference McPhee, Benson, Botha-Brink, Bordy and Choiniere2018; Mallet et al. Reference Mallet, Cornette, Billet and Houssaye2019; Etienne et al. Reference Etienne, Filippo, Cornette and Houssaye2021, Pintore et al. Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a). Furthermore, the increase of femoral robusticity in mammals is often coupled with enlarged muscle insertions located closer to the middle of the shaft, which enables application of greater moments around the hip joint of heavy animals (Hildebrand Reference Hildebrand1974; Polly Reference Polly and Hall2007; Mallet et al. Reference Mallet, Cornette, Billet and Houssaye2019). Accordingly, the placement of the fourth trochanter (i.e., attachment site for muscles between the femur and the tail) correlated similarly with body mass among non-avian dinosaurs and other archosauriforms, regardless of their locomotor habit (Coombs Reference Coombs1978; Parrish Reference Parrish1986; Carrano Reference Carrano1999; Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). In addition, the integration of both femoral robusticity and displacement of the fourth trochanter proved to be a reliable indicator of a more cursorial or graviportal morphology in bipedal saurischians (i.e., non-sauropod sauropodomorphs and theropods; Coombs Reference Coombs1978; Christiansen Reference Christiansen1998; Carrano Reference Carrano1999; Lefebvre et al. Reference Lefebvre, Houssaye, Mallison, Cornette and Allain2022; Pintore et al. Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a). Therefore, this well-established link between femoral robusticity, fourth trochanter displacement, and increase of body mass gives a broader comparative context for studying the evolution of gigantism across various taxa of obligate bipedal dinosaurs.
We highlighted that the same femoral specialization to gigantism occurred convergently within lineages of theropods, as shown by the significant K statistics below one and the positive allometric relationship between femoral robusticity and the increase of body mass (Figs. 3, 4A, 5, 6). Furthermore, we demonstrated that this convergent femoral robusticity evolved different degrees of specializations depending on the maximum body mass within each theropod lineage. Indeed, the same degree of femoral robusticity and fourth trochanter displacement was shared between several theropods with an estimated body mass of approximately 500 kg despite being distantly related within the theropod clade, such as the largest dromaeosaurid Utahraptor, the ceratosaurids Ceratosaurus nasicornis and Skorpiovenator, the early neotheropod Dilophosaurus, and the relatively small tyrannosaurid Alioramus (Table 1; Figs. 1, 3, 4A, 5, 6; Uta, Dil, Cer, Sko, Ali). Moreover, the level of femoral robusticity was globally lower in pennaraptorans and ceratosaurids (Figs. 3A, 5, 6). Therefore, the convergent increase of body mass led to similar morphological specializations in the femora of several theropod lineages, whether it led to gigantism (e.g., Ceratosauridae, Megalosauridae, Tyrannosauridae, Allosauridae, non-pennaraptoran maniraptoriforms) or not (pennaraptorans) and regardless of their phylogenetic relationships (Fig. 6). Further, theropods evolved various growth strategies that were previously demonstrated to link with the evolution of gigantism and miniaturizations in theropods (Padian et al. Reference Padian, de Ricqlès and Horner2001). For example, prolonged growth appears to have been favored by some carcharodontosaurids, whereas accelerated growth was conversely noted in other allosauroids as well as in some tyrannosauroids (Cullen et al. Reference Cullen, Canale, Apesteguía, Smith, Hu and Makovicky2020; Canale et al. Reference Canale, Apesteguía, Gallina, Mitchell, Smith, Cullen, Shinya, Haluza, Gianechini and Makovicky2022; D'Emic et al. Reference D'Emic, O'Connor, Sombathy, Cerda, Pascucci, Varricchio, Pol, Dave, Coria and Curry Rogers2023). Hence, femoral shape specializations to different body masses evolved convergently through different developmental strategies among theropods.
Femoral Specializations to Increasing Body Mass in Earlier Archosaurs
Early archosaur femora evolved similar femoral specializations to body mass increase regardless of their locomotor habits, as shown by large bipedal dinosaurs and quadrupedal pseudosuchians similarly evolving mediolaterally enlarged epiphyses and a more distally located fourth trochanter, much as in later theropods (Pintore et al. Reference Pintore, Delapré, Lefebvre, Botton-Divet, Houssaye and Cornette2022a; Figs. 3B,C, and 4A). Moreover, there appears to be a decoupling in fourth trochanter morphology that was more asymmetric in bipedal early archosauriforms than in quadrupeds and more rounded in the most robust femora than in the most gracile ones (Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). However, we did not highlight any variation of symmetry in the fourth trochanter of theropods, which were exclusively bipedal. Conversely, the fourth trochanter was prominent and rounded in larger theropods and flatter in the smaller ones (Figs. 3C, 4A). Even if the fourth trochanter evolved a similar prominence with increasing body mass in both bipedal non-sauropod sauropodomorphs and theropods, theropods did not reach the enormous body masses estimated for sauropods; hence, they never reduced their fourth trochanter to the rather flat morphology seen in sauropods (Lefebvre et al. Reference Lefebvre, Houssaye, Mallison, Cornette and Allain2022). This difference could indicate either that (1) beyond a certain body mass (i.e., greater than the heaviest theropods), fourth trochanter specializations to increasing body mass evolved from distal migration along the shaft to a more flat morphology in large secondary quadrupedal sauropodomorphs, as the fourth trochanter was already located close to the middle of the diaphysis in earlier bipedal sauropodomorphs, at least in Mussaurus and Plateosaurus compared with smaller, earlier-diverging taxa such as Eoraptor and Anchisaurus (Sereno et al. Reference Sereno, Martínez and Alcober2012; Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b), and was morphofunctionally constrained from shifting more distally; or (2) a common evolutionary pattern existed between the two clades of saurischians, which diverged after a certain amount of time and/or differences in traits.
However, there were a few differences between femoral specializations to body mass observed between early archosauriforms and theropods. The medial border of the femoral head was proximo-dorsally deflected (i.e., “upturned”) among the largest theropods, as seen in the large carcharodontosaurid Giganotosaurus, and ventrodistally deflected (i.e. “downturned”) among the smallest ones (Figs. 3A,B, and 4A; Gig), whereas it did not vary across early archosauriforms (Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). An upturned femoral head was discussed as an unambiguous synapomorphy of Carcharodontosauridae, which includes many giants (Harris Reference Harris1998; Brusatte and Sereno Reference Brusatte and Sereno2008; D'Emic et al. Reference D'Emic, Melstrom and Eddy2012; Canale et al. Reference Canale, Novas and Pol2015, Reference Canale, Apesteguía, Gallina, Mitchell, Smith, Cullen, Shinya, Haluza, Gianechini and Makovicky2022) but also some smaller, earlier relatives such as Concavenator and Neovenator (the earliest carcharodontosaurids; Hocknull et al. Reference Hocknull, White, Tischler, Cook, Calleja, Sloan and Elliott2009; Cuesta et al. Reference Cuesta, Ortega and Sanz2018). This feature was also present in Australovenator (the sister taxa of all carcharodontosaurids) but not in smaller, earlier relatives from other tetanuran clades. Harris (Reference Harris1998) noted that early non-averostran theropods and ceratosaurids had a downturned femoral head, whereas abelisaurids and non-carcharodontosaurid tetanurans had a rather horizontal femoral head. However, we observed that the femoral head of all early non-averostran neotheropods and ceratosaurians appeared relatively similar to the rather horizontal femoral head orientation of small-bodied taxa in non-carcharodontosaurid tetanurans (e.g., ornithomimids, therizinosaurids, and dromaeosaurids) instead of downturned (Figs. 3A,B, 4A, 5, and 6). We explain these observations contrasting with those of Harris (Reference Harris1998) by the inclusion of “miniaturized” theropods (i.e., less than 20 kg) in our dataset. Some of these taxa had even more downturned femoral heads than those of early non-averostran neotheropods and ceratosaurians (Table 1; Figs. 3A,B, 4A, 5, and 6). In addition, according to our analysis, large-bodied taxa (<1000 kg) within Allosauridae, Tyrannosauroidea, Megalosauroidea, and Therizinosauridae had an upturned femoral head similar to the condition described in carcharodontosaurids (Figs. 3A,B, 4A, 5, and 6). Bates et al. (Reference Bates, Benson and Falkingham2012) suggested that an upturned femoral head could have evolved as an adaptation to greater body mass in carcharodontosaurids, reducing the amount of stress that might be experienced at the femoral head and shaft by decreasing the bending moment arm of hip joint forces applied to the femoral head. Our results also suggest that this feature could result from a progressive specialization to large body mass that may have appeared convergently in several lineages of giant theropods and not exclusively in carcharodontosaurids and that it should be regarded as a continuum rather than a discrete character. Furthermore, the ancestral presence of an upturned femoral head in the carcharodontosaurid lineage but not in smaller, earlier relatives from other tetanuran clades suggests that the original morphology of these animals could have facilitated the evolution of gigantism within this clade, shedding light on why carcharodontosaurids evolved such large body size.
Another difference in the femoral specializations between theropods and early archosauriforms was the enlargement and proximal shift of the mediodistal crest (also termed “medial epicondyle” [Carrano and Sampson Reference Carrano and Sampson2008] and “craniomedial crest” [Hutchinson Reference Hutchinson2001]; Figs. 3A,C, and 4A). A well-developed mediodistal crest is a feature exclusive to neotheropods and most likely acted as part of the origin of the femorotibialis musculature (Hutchinson Reference Hutchinson2001; Carrano and Hutchinson Reference Carrano and Hutchinson2002). However, a well-developed crest can be characterized by different morphotypes, with a morphology varying from a smooth, subdued crested to a hypertrophied flange, and with varied extent along the shaft, traits described to distinguish various neotheropod clades (Ezcurra Reference Ezcurra2006; Carrano and Sampson Reference Carrano and Sampson2008). First, its enlargement toward the middle of the shaft by 25% at least of the femoral length was described as a carcharodontosaurid synapomorphy (Holtz et al. Reference Holtz, Molnar, Currie, Weishampel, Dodson and Osmolska2004; Brusatte and Sereno Reference Brusatte and Sereno2008). Second, a hypertrophied flange-like mediodistal crest (i.e., more prominent toward the medial side) is a diagnostic feature of abelisauroids (Carrano et al. Reference Carrano, Sampson and Forster2002; Tykoski and Rowe Reference Tykoski, Rowe, Weishampel, Dodson and Osmolska2004; Ezcurra Reference Ezcurra2006; Martinelli et al. Reference Martinelli, Marinho, Egli, Hechenleitner, Iori, Veiga, Basilici, Soares, Marconato and Ribeiro2019), and a long and crest-like crest morphology is present in ceratosauroids (Carrano and Sampson Reference Carrano and Sampson2008) and other non-abelisauroid neotheropods (Ezcurra Reference Ezcurra2006), whereas the crest is absent in all coelurosaurs except tyrannosauroids (Novas et al. Reference Novas, Ezcurra, Agnolin, Pol and Ortíz2012). Because we demonstrated a decoupling in the mediodistal crest extension between PC 1 (i.e., proximo-distal extension) and PC 2 (i.e., mediolateral extension), this observation is therefore consistent with our results showing that the two abelisauroids Masiakasaurus (a small noasaurid) and Carnotaurus (a large abelisaurid) clustered together in the PC 1–PC 2 morphospace (Figs. 3A,B, and 4A; Cer, Mas). Accordingly, these specimens should have displayed the greatest development of the mediodistal crest of our dataset (i.e., prominent, hypertrophied, and flange-like), but this was not the case. Indeed, we demonstrated that its proximal extension was similar to those of early neotheropods and other ceratosaurians, and to a lesser degree than in any other neotheropods with a more robust femoral morphology (Figs. 3A,B, and 4A). In addition, and even if its medial extension was greater than that of Tetanurae, it was similar to those of other ceratosaurians and early neotheropods too (Figs. 3A,B, and 4A). This inconsistency indicates that a “well-developed mediodistal crest” should be interpreted as a continuum from either more or less expanded toward the medial side (i.e., prominent) or more or less elongated toward the middle of the shaft (i.e., proximally extended). We demonstrated that the mediodistal crest extended toward the middle of the shaft in large theropods and that it was less prominent in maniraptoran theropods that were the closest relatives to birds (Figs. 3A,B,F, 4A,B, and 5). Moreover, a proximally elongated mediodistal crest in large theropods could indicate biomechanical modifications of femorotibial muscles and/or knee extensor tendons (Hutchinson Reference Hutchinson2001). Stronger muscular forces along the medial axis of the legs could improve mediolateral stability, vital for the bipedal locomotion of large theropods (Bishop et al. Reference Bishop, Clemente, Weems, Graham, Lamas, Hutchinson and Rubenson2017, Reference Bishop, Falisse, De Groote and Hutchinson2021). Therefore, the morphology of the mediodistal crest was decoupled between the evolution of gigantism (i.e., variation of proximal extension) and the phylogenetic proximity to the bird lineage (i.e., prominence/medial extension), as in the case of the fourth trochanter between small versus large and bipedal versus quadrupedal avemetatarsalians, pseudosuchians, and early archosauriforms (Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b).
In conclusion, the femur of several lineages of theropods convergently specialized to the increase of body mass and similarly to other saurischians and archosauriforms, suggesting that a similar femoral Bauplan was replicated during the evolution of archosauriform gigantism while integrating specific innovations for large bipeds.
The Influence of Ecological Habits on Femoral Robusticity
Cursoriality in Carnotaurus
Our results highlighted that not all giant theropods had a robust femoral morphology. Indeed, the largest abelisauroid, Carnotaurus (~1500 kg), had a femoral robusticity similar to that of smaller theropods from ~20 to ~500 kg with cursorial morphology, such as Velociraptor and Gallimimus (Figs. 3A, 6; Car, Vel, Gal). The ceratosaurian with the closest femoral robusticity was the smallest of all abelisauroids, Masiakasaurus (~20 kg), while medium-sized abelisauroids such as Skorpiovenator and Eoabelisaurus (~800 kg) had greater femoral robusticity than Carnotaurus (Figs. 3A, 6; Mas, Sko, Eoa). It was previously suggested that Carnotaurus had an unusually athletic morphology, on the basis of its tibial length being relatively longer than its femoral length and on the prominent, dorsally inclined transverse processes on the caudal vertebrae; the latter suggests an enlarged caudofemoralis muscle, which attaches on the fourth trochanter at the posterior side of the femur (Bonaparte et al. Reference Bonaparte, Novas and Coria1990; Mazzetta et al. Reference Mazzetta, Fariña and Vizcaíno1998; Persons and Currie Reference Persons and Currie2011). Persons and Currie (Reference Persons and Currie2011) speculated that the possible heightened athleticism of Carnotaurus might have evolved because of the competition they may have encountered with carcharodontosaurids, which lived in the same geographic area and were much larger than any abelisauroid. Whether this hypothesis is valid remains untested, but even if the carcharodontosaurids Giganotosaurus and Tyrannotitan appeared earlier than the first known occurrence of Carnotaurus, we show a definite morphofunctional distinction between those taxa. Indeed, the femoral morphology of carcharodontosaurids was among the most robust (e.g., Tyrannotitan and Giganotosaurus), contrasting with Carnotaurus having the most gracile morphology of all giant theropods (Table 1; Figs. 3A, 6; Car, Gig, Tyr).
Semi-aquatic Lifestyle in Spinosaurus
The femur of the megalosauroid Spinosaurus led to a surprisingly small body mass estimation (~600 kg; Table 1) given its extreme estimated body length of approximately 14 m (Sereno et al. Reference Sereno, Myhrvold, Henderson, Fish, Vidal, Baumgart, Keillor, Formoso and Conroy2022) and its high level of femoral robusticity (Figs. 3A, 6). Indeed, Spinosaurus's femoral robusticity was similar to that of the other large spinosaurid genera Suchomimus and megalosauroid Torvosaurus. Ibrahim et al. (Reference Ibrahim, Sereno, Dal Sasso, Maganuco, Fabbri, Martill, Zouhri, Myhrvold and Iurino2014) described the Spinosaurus femur as short and robust, resulting in an unusually low ratio of hindlimb length to total body length relative to other theropods (27% in Spinosaurus vs. 40% in Suchomimus and even greater in other theropods). This reduction of the hindlimb length while preserving a high femoral robusticity was suggested as indicating an aquatic lifestyle for Spinosaurus and potentially a quadrupedal habit (Ibrahim et al. Reference Ibrahim, Sereno, Dal Sasso, Maganuco, Fabbri, Martill, Zouhri, Myhrvold and Iurino2014, Reference Ibrahim, Maganuco, Dal Sasso, Fabbri, Auditore, Bindellini and Martill2020; Fabbri et al. Reference Fabbri, Navalón, Benson, Pol, O'Connor, Bhullar, Erickson and Norell2022). However, this hypothesis was extensively discussed by other authors, who favored a semi-aquatic to terrestrial lifestyle (Amiot et al. Reference Amiot, Buffetaut, Lécuyer, Wang, Boudad, Ding and Fourel2010; Hone and Holtz Reference Hone and Holtz2021; Sereno et al. Reference Sereno, Myhrvold, Henderson, Fish, Vidal, Baumgart, Keillor, Formoso and Conroy2022). Furthermore, this contrast between high femoral robusticity and low estimated body mass might be explained by the MDC-based body mass estimation method we used not being appropriate for a theropod with ecological and locomotor habits different from fully terrestrial bipedality. Another explanation could be that FSAC KK11888 was actually estimated to be only 76% of the maximum known body length for Spinosaurus (Sereno et al. Reference Sereno, Myhrvold, Henderson, Fish, Vidal, Baumgart, Keillor, Formoso and Conroy2022). Ibrahim et al. (Reference Ibrahim, Sereno, Dal Sasso, Maganuco, Fabbri, Martill, Zouhri, Myhrvold and Iurino2014) demonstrated that this individual was not a juvenile by estimating its age to be around 17 years old through histological analysis, and suggested that FSAC KK11888 was most likely a subadult, as Sereno et al. (Reference Sereno, Myhrvold, Henderson, Fish, Vidal, Baumgart, Keillor, Formoso and Conroy2022) also did. An absence of variation of femoral robusticity along ontogeny was also observed among juvenile and adult semi-aquatic Nile crocodiles (Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). Indeed, the femoral robusticity of these taxa was as high as in the Late Triassic quadrupedal semi-aquatic phytosaurs, quadrupedal pseudosuchians such as armored aetosaurs (e.g., Typothorax and Paratypothorax), and large bipedal sauropodomorph dinosaurs (e.g., Plateosaurus and Mussaurus; Pintore et al. Reference Pintore, Houssaye, Nesbitt and Hutchinson2022b). In light of the presumed at least semi-aquatic lifestyle of Spinosaurus, the high robusticity of its femur relative to its estimated age, body mass, and length appears consistent with the condition observed across different ontogenetic stages of semi-aquatic quadrupedal archosaurs. This observation reinforces the hypothesis that both the evolution of some degree of amphibious/semi-aquatic lifestyle and gigantism convergently led to a high femoral robusticity in archosauriforms, regardless of locomotor habit. Although a similar convergent phenomenon was demonstrated in the internal bone structure of both semi-aquatic and graviportal amniotes (Houssaye et al. Reference Houssaye, Waskow, Hayashi, Cornette, Lee and Hutchinson2016, Reference Houssaye, Martin, Boisserie and Lihoreau2021), this parallel has not been highlighted in the external morphology of their bones.
The Evolution of Miniaturization and Its Associated Femoral Specialization
As discussed by Benson et al. (Reference Benson, Hunt, Carrano and Campione2018), the multiple originations of gigantism represent only one facet of the evolutionary history of dinosaur body mass. Whereas it seems that there was a single allometric trajectory to evolve femoral robusticity along increasing body mass within theropods, it does not seem to be the case with femoral gracility along decreasing body mass (Figs. 3A, 6). Instead, the rather gracile femora were subdivided into two clusters, one containing most of the medium-sized cursorial taxa in our dataset (e.g., Gallimimus, Deinonychus, Masiakasaurus, Falcarius; Fig. 3A; Gal, Dei, Mas, Fal), and one containing only “miniaturized” theropods (i.e., with body masses between 1 to 20 kg, lighter than the lightest earlier taxa; e.g., the dromaeosaurid Bambiraptor and the troodontid Sinovenator; Figs. 3A, 6; Bam, Sin). Whether theropods experienced a sustained size decrease during their evolutionary history (Turner et al. Reference Turner, Pol, Clarke, Erickson and Norell2007; Lee et al. Reference Lee, Cau, Naish and Dyke2014) or a more punctuated mode of evolution (Carrano Reference Carrano, Carrano, Blob, Gaudin and Wible2006; Benson et al. Reference Benson, Hunt, Carrano and Campione2018), most studies agree that a drastic reduction in body mass occurred close to the bird lineage of Avialae (Brusatte et al. Reference Brusatte, Lloyd, Wang and Norell2014), specifically around the paravian node (Turner et al. Reference Turner, Pol, Clarke, Erickson and Norell2007; Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014). This theropod miniaturization was even previously documented as the plesiomorphic condition for all earlier bird relatives, which later facilitated the evolution of key innovations such as powered flight (Sereno Reference Sereno1999; Xu and Norell Reference Xu and Norell2004; Carrano Reference Carrano, Carrano, Blob, Gaudin and Wible2006; Turner et al. Reference Turner, Pol, Clarke, Erickson and Norell2007). We suggest that this accelerated miniaturization close to the bird lineage was associated with certain femoral specializations (i.e., long diaphysis relative to epiphyseal width, proximo-distally reduced fourth trochanter located close to the femoral head, downturned femoral head, mediodistal crest proximo-distally reduced), as represented by the cluster integrating the most gracile femora. Furthermore, this specialization is not observed among small early-diverging non-tetanuran theropods with an estimated body mass less than ~20 kg, because they are not included in the cluster (e.g., the Late Triassic Coelophysis and Masiakasaurus, the smallest ceratosaurian from the Late Cretaceous; Figs. 3A, 5; Coe, Mas). Thus, this specialization was not linked only to small body size but rather seemed restricted only to a few paravian clades, consistent with the main miniaturization event on line to birds (Figs. 3A, 6, 7).
However, some “miniaturized” paravians and non-paravian pennaraptorans had an unusually high femoral robusticity relatively to their estimated body mass (e.g., the smallest oviraptorid Conchoraptor and the dromaeosaurid Velociraptor; Figs. 3A, 6; Con, Vel) and did not plot along with other “miniaturized” theropods such as Bambiraptor and Sinovenator (Figs. 3A, 6; Bam, Sin). For example, Conchoraptor was the only “miniaturized” non-paravian pennaraptoran from our sample, as we estimated its body mass to be around 13 kg (Table 1; Figs. 3A, 7), which was consistent with other studies (Kundrát Reference Kundrát2007; Balanoff et al. Reference Balanoff, Bever, Rowe and Norell2013, Reference Balanoff, Bever and Norell2014). Therefore, the timing of the main miniaturization event on line to crown birds could have occurred earlier than the paravian node, at the pennaraptoran node (as displayed in Fig. 7), although the “miniaturization” of Conchoraptor could also be an independent convergent occurrence. Nevertheless, the femur of Conchoraptor did not specialize to “miniaturization” as in other dromaeosaurids and troodontosaurids from our dataset. Hence, specialization to miniaturization did not evolve in small oviraptorosaurids, whereas it did in the smallest dromaeosaurids and troodontids (Fig. 7). Conversely, some troodontids had an unusually low femoral robusticity given their estimated body mass (Table 1; Fig. 6; Tro), which may be related to the fairly cursorial morphology, such as relatively elongated legs and metatarsals for their body size observed for this clade (Holtz Reference Holtz1994). This observation further highlights that, while there was a main miniaturization event that was ancestral for several pennaraptoran clades, its associated femoral specialization occurred only later within paravians and may have been preserved during subsequent distinct events of body mass increase. This could indicate either that: (1) femoral specialization to miniaturization evolved once around the paravian node and was preserved along the later increase of body mass in troodontids (Sinovenator ~ 3 kg, Troodontidae indet. ~ 40 kg) but not in the later increase of body mass in dromaeosaurids (Velociraptor ~ 20 kg; Deinonychus ~ 100 kg; Utahraptor ~ 550 kg; Fig. 7); or (2) femoral specialization to miniaturization convergently evolved only among some paravians (Bambiraptor and troodontids) from a rather medium-sized cursorial morphology (Fig. 7). However, these hypotheses are complicated by uncertainties in dromaeosaurid phylogenetic relationships (e.g., Turner et al. Reference Turner, Makovicky and Norell2012). We caution that these evolutionary scenarios could vary depending on the phylogenetic positions of Bambiraptor and Utahraptor (the smallest and largest dromaeosaurids within our sample), which could have diverged earlier within dromaeosaurids than in our composite phylogeny (Fig. 1). Nevertheless, our results improve understanding of how miniaturization evolved close to the bird lineage. The femoral specialization to miniaturization evolved only from the paravian node and onward, potentially convergently between some dromaeosaurids and troodontids, but was not shared with non-paravian theropods of similar body mass or with paravians closely related to the relatively large dromaeosaurid Utahraptor (Figs. 3A, 6, 7). This raises the question of whether the origins of a rather “avian” femoral morphology (i.e., fused greater and lesser trochanters, no epiphyseal offset, reduced and concave mediodistal crest) across the pennaraptoran node were independent from the femoral specialization to miniaturization.
Was the Evolution of the Avian Femoral Morphology Independent from the Specialization to Miniaturization?
Taxa that were the most distantly related to the avian lineage had the femora with the most sigmoidal shape, with a slight medial offset of the femoral head in relation to the distal epiphysis, a convex and prominent mediodistal crest, and well-separated greater and lesser trochanters. Moreover, the lesser trochanter had a lower proximal extension than the greater trochanter (Figs. 3A,F–I, 4B, and 5). Conversely, taxa that were closer to the bird lineage had straighter femora with little to no epiphyseal offset, a reduced and concave mediodistal crest, along with extension of the lesser trochanter to the same proximal level as that of the greater trochanter (Figs. 3A,F–I, 4B, and 5). These latter morphological variations mostly occurred convergently, but to a greater extent in clades that were phylogenetically closer to the bird lineage (Fig. 5). Therefore, we interpret the combination of these femoral features to represent a more “avian”-like morphology.
Lesser Trochanter
The lesser trochanter in dinosaurs is thought to have been the attachment site for the iliotrochantericus caudalis (ITC) muscle (Hutchinson Reference Hutchinson2001), which was key in the transition from hip-based to knee-based locomotion toward crown birds and is thus characteristic of avian locomotion (Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Allen et al. Reference Allen, Kilbourne and Hutchinson2021). The proximal extension of the lesser trochanter, and its subsequent fusion with the greater trochanter, gradually evolved from earlier theropods toward birds, perhaps to accommodate an anteroproximal shift of the ITC insertion (Carrano Reference Carrano2000; Hutchinson Reference Hutchinson2001). This shift may have increased the capacity of the ITC to cause internal (medial) long-axis rotation of the femur, providing mediolateral forces in order to keep legs close to the sagittal midline and maintain balance during the stance phase in bipedal and erect archosaurs, especially with a more or less knee-based locomotion (Hutchinson and Gatesy Reference Hutchinson and Gatesy2000). It was previously demonstrated that an actual fusion between greater and lesser trochanters only evolved within maniraptoran theropods, except in therizinosaurs, which preserved a narrow intertrochanteric groove (Chiappe et al. Reference Chiappe, Norell and Clark1997; Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Holtz and Osmólska Reference Holtz, Osmólska, Weishampel, Dodson and Osmolska2004). Yet the proximal extension of the lesser trochanter (i.e., without evolving an actual fusion with the greater trochanter) evolved earlier than Maniraptora in Coelurosauria (Gauthier Reference Gauthier1986; Hutchinson Reference Hutchinson2001). Accordingly, we have highlighted that the greatest level of proximal extension of the lesser trochanter evolved earlier than Maniraptora along the theropod phylogeny in tyrannosaurids (Figs. 3A,C, and 8C). The proximal extension of the lesser trochanter was low in the Late Jurassic early-branching tyrannosaurid Guanlong and gradually increased toward the Late Cretaceous giant Tyrannosaurus and its close relatives. The ITC is an important component of weight support in single-limb stance (regardless of mechanism; i.e., abduction or long-axis rotation), and changes in its insertion reflect increases in the weight this muscle may have been supporting—either via absolutely larger masses (e.g., as in tyrannosaurids; Figs. 3A,C, and 8C) or via doing a greater fraction of the load-sharing as part of more knee-driven locomotion (Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Allen et al. Reference Allen, Kilbourne and Hutchinson2021). This suggests that a shift of the ITC insertion correlated with variation in body mass within at least one lineage (e.g., tyrannosauroids) not closely related to the avian lineage (Figs. 3A, 8, Supplementary Fig. S4; Gua, Trx). Therefore, this observation suggests a more complex link between body mass support, limb kinematics, and ITC insertions. However, we did not observe any variation of proximal extension in the earliest-branching neotheropods (non-averostran neotheropods and ceratosaurians) or in non-coelurosaurian tetanurans (megalosauroids and allosauroids; Figs. 3A, 8B). This observation suggests a mosaic evolution of the proximal extension of the lesser trochanter, which independently evolved to a different degree exclusively among clades that were closer to the bird lineage (Fig. 8B).
Epiphyses Offset
While no variation of lesser trochanter extension and mediodistal crest prominence was visible along the megalosauroid lineage, there was a greater epiphyseal offset in the early Lower Jurassic megalosauroid Piatnitzkysaurus than in the Late Jurassic megalosaurid Torvosaurus and spinosaurids from the Lower Cretaceous, Spinosaurus and Suchomimus, the latter having the most parallel alignment between the two epiphyses (Fig. 3A,H,I; Pia, Tor, Suc, Spi; Fig. 8, OFF). A lower epiphyseal offset reflects improved protraction/retraction motions of the femur, which is associated with more parasagittal limb posture (Carrano Reference Carrano2000; Hutchinson and Gatesy Reference Hutchinson and Gatesy2000), commonly associated with the archosaurian erect posture and, more relevant here, the evolution of more knee-based locomotion toward birds (Farlow et al. Reference Farlow, Gatesy, Holtz, Hutchinson and Robinson2000; Hutchinson and Allen Reference Hutchinson and Allen2009). Yet although a shift in epiphyseal offset evolved within theropods toward crown-group birds (Egawa et al. Reference Egawa, Griffin, Bishop, Pintore, Tsai, Botelho and Smith-Paredes2022), it also convergently appeared several times within other dinosaur clades (e.g., Ornithischia and Sauropodomorpha; Carrano Reference Carrano2000; Hutchinson Reference Hutchinson2001). Therefore, our observations further characterize the convergent nature of femoral epiphyseal offset by demonstrating its occurrence along non-coelurosaurian lineages, which diverged earlier than the avian lineage within Theropoda.
Mediodistal Crest
Finally, the variation of the mediodistal crest prominence had a wide phylogenetic span, with no apparent allometric relationship with body mass. Indeed, only the clades most distantly related to the bird lineage had a prominent convex mediodistal crest (Figs. 3A,F, 4B, and 8A). Among tetanurans, non-pennaraptorans had a somewhat intermediate “smooth” condition, with little to no medial extension along each lineage (megalosauroids, allosauroids, tyrannosauroids, therizinosauroids, ornithomimosaurids), except for some Allosaurus specimens having a surprisingly prominent crest. In contrast, all pennaraptorans had a clearly concave mediodistal crest (Fig. 8D). Therefore, modification from a rather convex to concave mediodistal crest may reflect a functional shift from earlier theropods to birds, which is consistent with its variation along the phylogeny of theropods of our dataset. However, and because we demonstrated that two (out of three) “avian” femoral features (i.e., lesser trochanter extension and epiphyseal offset) evolved convergently within Theropoda, we suspect that the phylogenetic signal of the mediodistal crest variation was overlooked by the K statistics (K < 1). Nevertheless, we discussed earlier that the proximal extension, instead of the prominence, of the mediodistal crest was associated with an increase of body mass. Thus, our result highlighted a decoupling between variations of body mass and the evolution of avian/knee-based locomotor mechanisms in the mediodistal crest morphology (Figs. 3, 4, 8). This may suggest a mosaic evolution in femorotibialis muscle anatomy along with body mass and locomotor variations.
Mosaic Evolution of a More “Avian” Femoral Morphology
Carrano (Reference Carrano2000) stated that dinosaurs had three fundamental transformations during the evolution of parasagittal hindlimb function: (1) a reorientation of the femoral head toward a fully medial direction, (2) the proximal extension of the lesser trochanter, and (3) changes of the ilium, the last of which is not directly relevant to our study. Our results highlighted that the two major femoral transformations were independent from variations of body mass (Figs. 3A,G–I, and 4B; Supplementary Fig. S4). Similar changes (i.e., lesser trochanter extension and epiphyseal offset) evolved convergently between three major clades of dinosaurs (theropods, sauropodomorphs, and ornithischians; Coombs Reference Coombs1978; Gauthier Reference Gauthier1986; Novas Reference Novas1996; Carrano Reference Carrano2000; Hutchinson Reference Hutchinson2001), and we demonstrated additional convergent evolutions of the same features within neotheropods as well (Fig. 8). We interpret this as further emphasizing that transformations experienced by Mesozoic theropods along the lineage to birds did not follow a gradual increase from the taxa most distantly related to the most closely related to birds, as noted in Carrano (Reference Carrano2000; e.g., megalosauroids from deeply nested lineages had a more similar femoral “avian” morphology to pennaraptorans and tyrannosauroids than other theropods more closely related to birds). Instead, we suggest that repeated iterations of femoral transformations within several lineages may have led to more derived bipedal functions (i.e., more “avian” style of locomotion) in dinosaurs, ultimately shifting from more hip- to more knee-driven locomotor habits in paravians. Consequently, the “avian” femoral morphology may have evolved in a mosaic fashion, which may have started to evolve deep within the evolutionary history of dinosaurs and followed, just like miniaturization, either a long sustained trend or a somewhat stepwise pattern, but without any clear, binary distinction between “non-avian” and “avian” femoral morphology (Hutchinson and Allen Reference Hutchinson and Allen2009; Brusatte et al. Reference Brusatte, Lloyd, Wang and Norell2014; Benson et al. Reference Benson, Hunt, Carrano and Campione2018).
Specialization to Miniaturization Evolved Only in the Most “Avian” Femoral Morphology
Although it is evident that some of the most “avian” femoral morphologies (i.e., straight femur with no epiphyseal offset, reduced and concave mediodistal crest, lesser trochanter proximally extending to the level of greater trochanter) appeared only in “miniaturized” taxa and close relatives, we demonstrated no apparent link with variation in body mass (Figs. 3A, 5–7). However, while the first occurrence of miniaturization may have evolved in maniraptorans, we demonstrated that the actual femoral specializations to miniaturization (i.e., long diaphysis relative to epiphyseal width, proximo-distally reduced fourth trochanter located close to the femoral head, downturned femoral head, mediodistal crest proximo-distally reduced) appeared only in some paravians, which already integrated fundamental modifications toward a rather “avian” femoral morphology (Figs. 3, 5, 7, 8; Supplementary Fig. S4; i.e., a concave mediodistal crest appeared only in pennaraptorans; a complete fusion between greater and lesser trochanter evolved in maniraptorans; a full medial offset of the femoral head first evolved in megalosauroids). Thus, only theropods with the most “avian”-like femoral morphology evolved a highly derived specialization to miniaturization (Figs. 3A, 7, 8). This pattern is consistent with findings from Benson et al. (Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014) that maintaining small body size was key to evolving essential avian features. Furthermore, an increase of body mass evolved within paravians, as seen with Utahraptor having one of the most “avian” femoral morphologies of the dataset while having an estimated body mass of ~550 kg. This striking example highlights that theropods closely related to birds could evolve “avian”-like femoral features while independently evolving the same femoral specializations to gigantism as any other earlier giant relatives, but without reverting to a plesiomorphic condition of overall femoral morphology. Therefore, even if a more “avian” femoral morphology is a specialization to miniaturization, it did not limit later evolutionary increases of body mass. Further, large body mass also evolved in Cenozoic crown-group birds to an estimated maximum of 500 to 700 kg (e.g., Dromornis stirtoni; see Chinsamy et al. Reference Chinsamy, Handley and Worthy2023). Although much larger than the inferred mass of the ancestral bird (Benson et al. Reference Benson, Campione, Carrano, Mannion, Sullivan, Upchurch and Evans2014), their estimated maximum body mass is still well below the upper limit reached among non-avian theropods. This has led to different hypotheses that body mass increases in crown birds were constrained by various factors, such as the resistance of their eggshell presumably limiting the body mass of the incubating individual (Deeming and Birchard Reference Deeming and Birchard2009). Another possible constraint would have been the more anteriorly located center of mass associated with a more “avian” knee-based locomotion, which would have constrained the femoral morphology to withstand specific bending stresses, hence limiting the evolution of specialization to extreme body mass increases, which may have been more possible in the typically more vertically oriented femora of many non-avian theropods (Gatesy Reference Gatesy1991; Chan Reference Chan2017). Indeed, a more knee-based locomotion (i.e., more flexed hindlimb, more horizontal femur) is associated with increased torsional loads on the femur and relies less on femoral abductors than hip-driven bipedal locomotion does (Carrano Reference Carrano2000; Hutchinson and Gatesy Reference Hutchinson and Gatesy2000; Hutchinson and Allen Reference Hutchinson and Allen2009), hence differing from the condition in some earlier-diverging theropods. We suggest that this relationship between specialization of body mass increase and a somewhat “avian” femoral morphology may have been facilitated by the morphological decoupling of mediodistal crest evolution. Further investigations of femoral morphology and how it related to biomechanical constraints across miniaturized to giant Mesozoic and Cenozoic paravians could yield new evidence on how repeated variations of body mass in theropods could have evolved along the line to birds, independent of or related to the evolution of the avian Bauplan.
Conclusion
Our study demonstrates that the femoral morphology of obligately bipedal theropods consistently evolved similar specializations to gigantism regardless of their phylogenetic affinities, maximal body mass, and (inferred) degree of hip-/knee-driven locomotor mechanism. Femoral robusticity increased convergently between several theropod lineages through a consistent shift of the fourth trochanter toward the middle of the shaft, a mediolateral enlargement of the epiphyses, a proximodorsal reorientation of the medial side of the femoral head, and a proximal enlargement of the mediodistal crest. Some of these features, such as the distal shift of the fourth trochanter and the enlargement of the epiphyses, are common specializations to gigantism known in quadrupedal dinosaurs and other vertebrates. However, variations relative to the proximodorsal orientation of the femoral head and the proximal extension of the mediodistal crest may be exclusive to gigantism in theropod dinosaurs. These features could reflect the functional constraints experienced by femora of giant bipeds with increasingly vertically oriented hindlimbs. In contrast, some paravians demonstrated an extreme degree of femoral gracility, which we interpreted as a specialization to miniaturization. This specialization was not common to all small, potentially “miniaturized” (less than ~20 kg), derived theropods (i.e., not all pennaraptorans). Only “miniaturized” pennaraptorans that acquired at least three key femoral “avian” features (i.e., fusion between the greater and lesser trochanters, fully medial offset of the femoral head, concave mediodistal crest) displayed a strong specialization to miniaturization. Accordingly, a rather “avian” femur did not prevent the later evolution of specialization to gigantism in Utahraptor, which is large relative to other dromaeosaurids (~550 kg). Furthermore, we hypothesized a decoupling between the evolution of gigantism and more knee-based locomotion in the mediodistal crest morphology. While the crest's proximal extension was linked to an increase in body mass across several lineages, its shape varied from more convex and prominent in the earliest-branching theropods to concave and reduced in the closest relatives to the bird lineage. Therefore, the evolution of a more avian femoral morphology was independent from the repeated occurrence of gigantism, but the two were not mutually exclusive. A thorough investigation of the convergent femoral specializations to large body mass within theropods, bipedal early sauropodomorphs, and more or less quadrupedal/bipedal habits in ornithischians could yield new insights into how dinosaurian hindlimbs accommodated the evolution of gigantism without gradually evolving more knee-driven “avian” locomotion.
Acknowledgments
The data acquisition step of this study was largely impeded by the COVID-19 pandemic because of traveling restrictions. Creating a dataset of this quality would have never been possible without the help of many people around the world. We are indebted to them for providing high-quality data (3D models and photographs) of theropod femora, which were essential to our study: J. Canale and D. Palombi (MMCH); D. Pol (MPEF); T. Szczygielski and J. Słowiak (ZPAL); C. Kammerer, A. Canoville and L. Zanno (NCSM); C. Levitt-Bussian, R. B. Irmis, C. Webb and P. Policelli (UMNH); J. B. Scannella (MOR); B. Strilisky and D. M. Henderson (TMP); A. A. Farke (Raymond M. Alf Museum, USA); X. Yao (Yunnan University, China); J. Molnar (New York Institute of Technology, USA); C. Yu (AMNH); H. Mallison (Palaeo 3D, Germany); P. J. Makovicky (University of Minnesota, USA); P. Sereno; S. Baumgart (FMNH); D. Vidal (Universidad Autónoma de Madrid, Spain); N. Myhrvold (Intellectual Ventures, USA); T. Souter (Independent designer, France); and X. Xing (IVPP). We also warmly thank C. Colin-Fromont, V. Pernègre, J. Barbier, and R. Allain (MNHN), J.-F Tournepiche and D. Augier (Musée d'Angoulême, France). We thank R. Allain, R. Cornette, E. Guilbert, C. Bader, C. Etienne, and R. Lefebvre (MNHN), S. C. R. Maidment (Natural History Museum, UK), S. E. Pierce (MCZ), S. L. Brusatte (University of Edinburgh, UK), and C. Mallet (Institut Royal des Sciences Naturelles, Belgium) for helpful discussions and recommendations about analyses and interpretation of data. Finally, we thank three reviewers (two anonymous and C. T. Griffin) for their constructive feedback and stimulating discussions, which improved the quality of our article.
Funding
This research was supported by the European Research Council under Horizon 2020 Starting Grant GRAVIBONE (715300 to A.H.) and Advanced Grant DAWNDINOS (695517 to J.R.H.); an Australian Government Research Training Program Scholarship (to P.J.B.), the Paleontological Society (Robert J. Stanton and James R. Dodd Award, to P.J.B.), and the International Society of Biomechanics (Matching Dissertation Grant, to P.J.B.).
Competing interests
The authors declare no competing interest.
Data Availability Statement
Data available from the Dryad Digital Repository https://doi.org/10.5061/dryad.tx95x6b52.