Introduction
Ankylosaurian dinosaurs are one of the most peculiar groups of reptiles in the fossil record due to their extensive body covering of osteoderms creating a massive body armor that is unlike any other in the evolution of dinosaurs (e.g., Kirkland et al. Reference Kirkland, Alcalá, Loewen, Espílez, Mampel and Wiersma2013; Brown et al. Reference Brown, Henderson, Vinther, Fletcher, Sistiaga, Herrera and Summons2017). Osteoderm size and shape are diverse among ankylosaurian groups and even across distinct body regions within the same individual (e.g., tail knobs, ossicles, spines, and thoracic osteoderms; Brown et al. Reference Brown, Henderson, Vinther, Fletcher, Sistiaga, Herrera and Summons2017). However, the usage of osteoderm features to better understand ankylosaurian phylogeny and evolutionary patterns is still in its infancy. For instance, using osteoderm characters for phylogenetic reconstructions of ankylosaur phylogeny is relatively recent (Hill Reference Hill2005; Burns and Currie Reference Burns and Currie2014; Arbour and Currie Reference Arbour and Currie2016; Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021). In studies of osteoderm variation among distinct ankylosaurian groups, paleohistological analyses indicate some major differences in microstructure patterns (Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014). However, such differences have rarely been investigated in a macroevolutionary context using quantitative methods.
Additionally, the ankylosaurian fossil record is marked by major geographic sampling discrepancies, such as the larger sampling of ankylosaurians in the Northern Hemisphere compared with the southern continents (Arbour and Currie Reference Arbour and Currie2016; Arbour et al. Reference Arbour, Zanno and Gates2016). This biogeographic sampling discrepancy inherently limits the taxonomic and morphological information available on this group on a global scale, which may impact our understanding of their taxonomic and morphological diversity, as shown by recent discoveries (e.g., Maidment et al. Reference Maidment, Strachan, Ouarhache, Scheyer, Brown, Fernandez, Johanson, Raven and Barrett2021; Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; Frauenfelder et al. Reference Frauenfelder, Bell, Brougham, Bevitt, Bicknell, Kear, Wroe and Campione2022). Therefore, further sampling of ankylosaur remains from continents derived from Gondwana is critical to better understand their evolutionary trajectories.
The Antarctic continent is the most underexplored frontier of paleontology in the Southern Hemisphere (and perhaps globally). Antarctica represents a significant portion of all landmasses in the Southern Hemisphere, but only very few regions have exposed sedimentary outcrops that can provide insights into the evolution of life on Earth (e.g., Olivero Reference Olivero2012b; Burton-Johnson and Riley Reference Burton-Johnson and Riley2015; Piovesan et al. Reference Piovesan, Correia Filho, Melo, Lacerda, Dos Santos, Pinheiro, Costa, Sayão and Kellner2021). Despite such limitations, several vertebrate remains have been discovered in Antarctica over the past few decades. These discoveries include plesiosaurs (e.g., O'Gorman et al. Reference O'Gorman, Otero, Reguero and Gasparini2019; Brum et al. Reference Brum, Simões, Souza, Pinheiro, Figueiredo, Caldwell, Sayão and Kellner2022), pterosaurs (Kellner et al. Reference Kellner, Rodrigues, Costa, Weinschütz, Figueiredo, Souza, Brum, Eleutério, Mueller and Sayão2019), and dinosaurs (Lamanna et al. Reference Lamanna, Case, Roberts, Arbour, Ely, Salisbury, Clarke, Malinzak, West and O'Connor2019); for an updated survey of fossil vertebrates in the Upper Cretaceous of Antarctica, see Reguero et al. (Reference Reguero, Gasparini, Olivero, Coria, Fernández, O'Gorman, Gouiric-Cavalli, Hospitaleche, Bona, Iglesias, Gelfo, Raffi, Moly, Santillana and Cárdenas2022). Importantly, this includes one ankylosaurian species, Antarctopelta oliveroi Salgado and Gasparini, Reference Salgado and Gasparini2006— although there is some dispute on the validity of this taxon (Arbour and Currie Reference Arbour and Currie2016; Rozadilla et al. Reference Rozadilla, Rolando, Motta, Novas and Board2016).
Since 2007, the project PALEOANTAR (organized by the Museu Nacional/Universidade Federal do Rio de Janeiro, Brazil) has been working in the Antarctic Peninsula to shed light on the diversification and evolutionary history of Antarctic ecosystems in the deep past (Lima et al. Reference Lima, Sayão, Ponciano, Weinschütz, Figueiredo, Rodrigues, Bantim, Saraiva, Jasper, Uhl and Kellner2021; Kellner Reference Kellner2022; Santos et al. Reference Santos, Piovesan, Guzmán, Usman, Weinschütz, Santos, Oliveira, Figueiredo, Ricetti, Wilner, Sayão and Kellner2022), including studies on Mesozoic vertebrates (e.g., Kellner et al. Reference Kellner, Simões, Riff, Grillo, Romano, Silva, Ramos, Carvalho, Sayão, Oliveira and Rodrigues2011, Reference Kellner, Rodrigues, Costa, Weinschütz, Figueiredo, Souza, Brum, Eleutério, Mueller and Sayão2019; Brum et al. Reference Brum, Simões, Souza, Pinheiro, Figueiredo, Caldwell, Sayão and Kellner2022) and invertebrates (e.g., Pinheiro et al. Reference Pinheiro, Saraiva, Santana, Sayão, Figueiredo, Rodrigues, Weinschütz, Ponciano and Kellner2020; Videira-Santos Reference Videira-Santos, Scheffer, Ponciano, Weinschütz, Figueiredo, Rodrigues, Sayão, Riff and Kellner2020; Piovesan et al. Reference Piovesan, Correia Filho, Melo, Lacerda, Dos Santos, Pinheiro, Costa, Sayão and Kellner2021). Here, we report on new ankylosaurian materials recovered from Antarctica during the PALEOANTAR expedition in 2015. The material consists of ankylosaurian osteoderms found on Santa Marta Cove, James Ross Island (James Ross Basin, Antarctic Peninsula), which come from the same levels where Antarctopelta was previously found. Combining paleohistological, phylogenetic, and morphometric approaches, we assign these new specimens to late-diverging nodosaurids, thus expanding the phylogenetic diversity of Antarctic ankylosaurians. Further, we compiled osteoderm data from the literature that—combined with the new materials described here—provide the first detailed quantitative analysis of the evolution of the ankylosaurian body armor microstructure. We show that late-diverging nodosaurids (including our new specimens) expanded ankylosaur morphospace by decreasing the cortical thickness, which may be associated with a structural shift in their osteoderms by a major structural fiber rearrangement. We also indicate that the expansion of the core and its vascularization in late-diverging nodosaurids is similar to that observed in some titanosaurs, suggesting both an effective body armor function and a co-opted function of calcium remobilization for physiological balance. Such plastic physiology may have played a key role in the dispersal potential of this group toward higher latitudes, which present a greater disparity of climatic conditions throughout the year compared with lower latitudes.
Materials and Methods
Materials
The new specimens reported here were collected on James Ross Island and comprise three isolated osteoderms (CAV-A4, CAV-A5, and CAV-A10). The specimens and paleohistological slides are housed in the Laboratório de Paleobiologia e Microestruturas, Centro Acadêmico de Vitória, Universidade Federal de Pernambuco (Vitória de Santo Antão, State of Pernambuco, Brazil—CAV).
Anatomical Terminology
We followed the terminology of Scheyer and Sander (Reference Scheyer and Sander2004) and Burns and Currie (Reference Burns and Currie2014) to describe the anatomical osteoderm orientation, in which the core is the internal cancellous region (Burns and Currie Reference Burns and Currie2014; Fig. 1), while the cortex comprises the basal (toward the internal region of the bone) and superficial layers of compact bone (or external, toward the surface of the bone; Scheyer and Sander Reference Scheyer and Sander2004). The core could be filled by cancellous bone, with a high density of resorption cavities, or be more compacted. The cortex comprises bone that is more compact and, in some cases, exhibits a decrease in vascularization. The microstructural nomenclature followed Francillon-Vieillot et al. (Reference Francillon-Vieillot, de Buffrénil, Castenet, Geraudie, Meunier, Sire, Zylberberg, de Ricqlès and Carter1990). We used the term “structural fibers” to refer to the mineralized collagen fibers fully incorporated in the mineralized bone matrix (Scheyer and Sander Reference Scheyer and Sander2004). We adopted the term “interwoven structural fiber bundles” (ISFB; sensu Scheyer and Sander Reference Scheyer and Sander2004) to describe bundles of structural fibers within the bone matrix. In the core, the structural fibers are sparser, while in the cortex, they are more organized (e.g., the orthogonal arrangement in nodosaurids; perpendicular to the external surface in ankylosaurids; Burns and Currie Reference Burns and Currie2014).
Schematic osteoderm showing the anatomical terminology adopted and the linear measurements performed here. Abbreviations: B, basal thickness; Co, core thickness; S, superficial thickness; T, total thickness.
Paleohistology
The paleohistological protocols followed those commonly used in current literature (Chinsamy and Raath Reference Chinsamy and Raath1992; Chinsamy Reference Chinsamy2005; Lamm Reference Lamm, Padian and Lamm2013). We embedded the samples into clear epoxy resin Resapol T-208 and sectioned them with a precision router. The plane of section in CAV-A4 followed the longitudinal keel, which exhibits the widest preserved bone area to section. The sectioning plane was indeterminate in the other specimens due to their fragmentary nature. The blocks were fixed in a slide and polished in a wet metallographic polishing machine, Arotec AROPOL E until they reached a final thickness of 30–60 microns. The photomicrographs were taken under a cross-polarized microscope (Zeiss) and photographed by an Axiocam.
Morphometric Dataset
Most of the ankylosaur dataset was compiled from Burns and Currie (Reference Burns and Currie2014). We added new measurements from the literature, comprising the nodosaurid Ahshislepelta (Burns and Sullivan Reference Burns and Sullivan2011), the ankylosaurine UMNH 12675 (Loewen et al. Reference Loewen, Burns, Getty, Kirkland, Vickaryous, Titus and Loewen2013), Antarctopelta (Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019) and the titanosaurs Saltasaurus (Chinsamy et al. Reference Chinsamy, Cerda and Powell2016), and MCS-Pv 181 and MCS-Pv 182 (Cerda et al. Reference Cerda, García, Powell and Lopez2015). We included titanosaurs due to their cortex and core microstructure being quite like the one observed in some ankylosaurians, therefore having the potential of being mistakenly identified as ankylosaurian osteoderms or vice versa. Although more titanosaur specimens have been described with osteoderms (e.g., D'Emic et al. Reference D'Emic, Wilson and Chatterjee2009; Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011), we could not include them here, because (1) the photomicrograph does not figure the entire section—with cortex and core (e.g., Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011); (2) the published figures do not identify the boundaries of cortex and core—as in Rapetosaurus (see Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011)—and the figure is not in an adequate quality to delimit them; and (3) the literature data only reported the external anatomy of the osteoderms (e.g., D'Emic et al. Reference D'Emic, Wilson and Chatterjee2009).
Our dataset is available in Brum et al. (Reference Brum, Eleutério, Simões, Whitney, Souza, Sayão and Kellner2023) and comprises 8 ankylosaurid, 16 nodosaurid, 3 titanosaurian, 2 indeterminate ankylosaurian, and 2 Antarctopelta (of controversial taxonomy and systematic affinities) osteoderms, besides 2 of the new specimens reported here—we could not measure CAV-A10, as it is too fragmentary. The images obtained from the literature were measured using ImageJ 1.52a (Schneider et al. Reference Schneider, Rasband and Eliceiri2012).
Linear measurements of the osteoderms were converted into relative osteoderm ratio thickness (Fig. 1), following Burns and Currie (Reference Burns and Currie2014), for a direct shape comparison (i.e., excluding size as a variable): superficial index—superficial cortex/total osteoderm (SI = 100*S/T); basal index—basal cortex/total osteoderm (BI = 100*B/T); cortical index (CI = 100*(S + B)/T); and core index—core/total osteoderm (CoI = 100*Co/T). All values were log-transformed (log10) for subsequent statistical analyses.
Statistical Analyses
We performed ordinary least-square (OLS) linear regressions to evaluate the correlation between morphometric variables and to find a predictor model for dinosaurian osteoderm shape. Although some studies suggest the utilization of reduced major axis (RMA) linear regressions instead of OLS for inferring species allometric relationships (e.g., Warton et al. Reference Warton, Wright, Falster and Westoby2006), we followed the recommendation of Kilmer and Rodríguez (Reference Kilmer and Rodríguez2017) on using OLS instead of RMA. The RMA uses the ratio of the standard deviation between the y- and x-axes to calculate the slope, as it assumes that the error distribution should be proportional between the two axes (i.e., assumes symmetry in the error distribution between axes). However, this is rare in allometric relationships, as the x-axis in these almost always refers to body size (constant variable), while the y-axis refers to various trait values. On the other hand, OLS considers the covariation between axes relative to the variation on the x-axis to calculate the slope (not assuming symmetry in error distribution between axes). This reduces error when estimating the angular coefficient, achieving functional scaling relationships (Kilmer and Rodríguez Reference Kilmer and Rodríguez2017). Our study focuses on the trait–trait relationship (e.g., SI × CI), with most of our variables being interchangeable between axes. However, we cannot assume a symmetric distribution of error among variables, because measurements for each variable were taken by different authors, and so error distribution is likely to be asymmetric between variables (and thus between axes). Therefore, we chose to employ OLS instead of RMA.
To test whether each osteoderm ratio variable can discriminate between different groups of dinosaurs (ankylosaurids, nodosaurids, and titanosaurs) for taxonomic purposes, we performed nonparametric Kruskal-Wallis (KW) tests. Our choice for nonparametric tests is due to the small sample size for each group and the nonnormal distribution of the data (see the tests for normality in Supplementary Fig. 2). Group assignment was based on the results of the phylogenetic analysis conducted here (see “Ankylosaur Phylogeny and Assignment of CAV Specimens”). The indeterminate ankylosaurians (TMP 1998.98.1 A.2 and TMP 1987.113.4 A.1) and the osteoderms of Antarctopelta were not included within any of the three dinosaurian groups due to (1) the uncertain phylogenetic placement of Antarctopelta between previous and more recent phylogenetic analyses (e.g., Arbour and Currie Reference Arbour and Currie2016; Arbour et al. Reference Arbour, Zanno and Gates2016; Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; Frauenfelder et al. Reference Frauenfelder, Bell, Brougham, Bevitt, Bicknell, Kear, Wroe and Campione2022); and (2) recent analyses inferring Antarctopelta as external to the two main ankylosaurid groups, being the only early-diverging ankylosaurian included here (Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; Frauenfelder et al. Reference Frauenfelder, Bell, Brougham, Bevitt, Bicknell, Kear, Wroe and Campione2022; herein).
Further, we compared the decline in relative overall cortical thickness in osteoderms among groups using an exponential bivariate regression (EBR) between CI and T, as previously performed by Burns and Currie (Reference Burns and Currie2014). In the latter, specimens were assigned to three ankylosaurian groups—Ankylosauridae, Nodosauridae, and “Polacanthidae” (Burns and Currie Reference Burns and Currie2014: fig. 8). Most of our group assignments correspond to theirs (see Supplementary Table 1), except “Polacanthidae,” which are assigned to nodosaurids herein based on their more recent phylogenetic classification (e.g., Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; this study, see “Ankylosaur Phylogeny and Assignment of CAV Specimens”). Our new dataset is further distinct from Burns and Currie (Reference Burns and Currie2014), as it includes the ankylosaurid UMNH 12675, the nodosaurids Ahshislepelta, and the ankylosaur Antarctopelta (see Supplementary Table 1). To test which group best fit the EBR, we used the Akaike information criterion (AIC).
To assess the morphospace occupation and test group discrimination, we performed multivariate linear discriminant analyses (LDA). Although the principal component analysis (PCA) partially overlaps the observed data distribution depicted by OLS, we report it to provide an additional and easier visualization of morphospace occupation of ankylosaurian groups. All statistical analyses were performed in the software PAST v. 4.02 (Hammer et al. Reference Hammer, Harper and Ryan2001), with α = 0.05.
Phylogenetic Analysis
We scored the CAV specimens in the data matrix of Soto-Acuña et al. (Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021), which resulted in 68 taxa and 189 characters (all treated as unordered, thus treating all possible state transitions as equally likely). The dataset includes 23 characters related to postcranial osteoderms and 4 related to osteoderm microstructure (characters 157 to 160). The entire data matrix is available in Brum et al. (Reference Brum, Eleutério, Simões, Whitney, Souza, Sayão and Kellner2023).
We performed equal weights maximum parsimony analyses in TNT v. 1.5 (Goloboff and Catalano Reference Goloboff and Catalano2016) using traditional search (tree bisection reconnection), with 1 random seed, 6000 replicates, and 100 trees to save per replication. We retained suboptimal trees by 10 steps with a 0.01 relative difference. The first analysis was conducted without the new CAV specimens to track the microstructural osteoderm characters. The second analysis included the CAV specimens to verify that we could assign them to an ankylosaur ingroup solely by osteoderm microstructural characters.
Institutional Abbreviations
Centro Acadêmico de Vitória, Universidade Federal de Pernambuco, Pernambuco, Brazil (CAV-A); Colección Paleontológica de Coahuila, Museo del Desierto, Saltito, Mexico (CPC); Delaware Museum of Natural History, Delaware, U.S.A. (DMNH); Goldfuss-Museum, Steinmann Institute for Geology, Mineralogy and Paleontology, University of Bonn, Germany (IPB); Museo Regional Cinco Saltos, Río Negro Province, Argentina (MCS-Pv); Museo Municipal Carmen Funes, Plaza Huincul, Neuquén, Argentina (MLP); Museum of Western Colorado, Grand Junction, Colorado, U.S.A. (MWC); Natural History Museum, London, U.K. (NHMUK); Instituto Miguel Lillo, Tucumán, San Miguel Tucumán, Argentina (PVL); State Museum of Pennsylvania, Harrisburg, Pennsylvania, U.S.A. (SMP); Royal Tyrrell Museum of Palaeontology, Drumheller, Alberta, Canada (TMP); University of Alberta Laboratory for Vertebrate Paleontology, Edmonton, Alberta, Canada (UALVP); Utah Museum of Natural History, University of Utah, Utah, U.S.A. (UMNH); Zoological Institute of Paleobiology, Polish Academy of Sciences, Warsaw, Poland (ZPAL).
Results
Systematics
Dinosauria Owen, Reference Owen1842
Ornithischia Seeley, Reference Seeley1888
Thyreophora Nopsca, Reference Nopsca1915
Ankylosauria Osborn, Reference Osborn1923
Nodosauridae indet. Marsh, Reference Marsh1890
See Figures 2–4 and Table 1.
Photomicrographs of the transversal section slide of the osteoderm CAV-A4. A, Detail of the specimen and the panoramic slide. B, Detail of the alignment of osteocyte lacunae and the bone matrix organized into bundles with structural fibers, forming diffuse interwoven structural fiber bundles (ISFBs). Structural fibers assuming an orthogonal arrangement in the core under normal transmitted (C) and cross-polarized light (D). E, Detail of the ISFBs in the core region close to the basal cortex. F, Orthogonal arrangement of the structural fibers in the core, close to the margin of the osteoderm, under cross-polarized light. G, Superficial cortex in detail, with structural fibers being more perpendicular to the external surface. H, Erosion cavity with irregular surfaces, indicating an active resorption process in the core of the osteoderm. Abbreviations: RVC, reticular vascular canals; SF, structural fibers.
Photomicrographs of the transversal section slide of the keeled osteoderm CAV-A5. A, Detail of the specimen and panoramic slide. B, The interwoven structural fiber bundles (ISFBs) are more disorganized in the central region of the core. Orthogonal arrangement of the structural fibers in the basal region of the core, rich in reticular vascular canals, under normal transmitted (C, E) and cross-polarized light (D, F). G, Detail of the superficial cortex, with structural fibers more perpendicular to the external surface. H, The erosional cavity in the core of the osteoderm, with detail on its irregular inner margins and close to a micro-crack (arrow). Abbreviations: RVC, reticular vascular canals; SF, structural fibers.
Photomicrographs of the osteoderm CAV-A10. A, Detail of the specimen and panoramic slide. B, Photomicrograph of the external cortex rich in vascular canals with ISFBs. Detail of the orthogonal arrangement of the interwoven structural fiber bundles (ISFBs) in the contact between the core and external cortex under normal transmitted (C) and cross-polarized light (D). E, F, A core region rich in longitudinal vascular canals and showing an orthogonal arrangement of ISFBs under cross-polarized light. Abbreviations: LVC, longitudinal vascular canals; SF, structural fibers; SO, secondary osteons.
Thickness measurements of the specimens from James Ross Island. The values in W, S, B, S+B, and Co are in millimeters. Abbreviations: B, basal cortical bone thickness; BI, basal index; CI, cortical index; Co, core thickness; CoI, core index; S, surface cortex thickness; SI, superficial index; T, total osteoderm thickness; W, maximum width of whole osteoderm. *Partial measurement.
Locality and Horizon
The James Ross Basin outcrops are mostly represented by sequences of marine siliciclastic and volcanoclastic rocks (Barremian–Eocene; Olivero Reference Olivero2012b; Burton-Johnson and Riley Reference Burton-Johnson and Riley2015; Piovesan et al. Reference Piovesan, Correia Filho, Melo, Lacerda, Dos Santos, Pinheiro, Costa, Sayão and Kellner2021). The new specimens came from Santa Marta Cove (northeast of James Ross Island), which is marked by exposures of the upper levels of the Snow Hill Island Formation (for age and paleoenvironmental discussions, see Piovesan et al. Reference Piovesan, Correia Filho, Melo, Lacerda, Dos Santos, Pinheiro, Costa, Sayão and Kellner2021) and corresponds to the upper levels of the MG sequence (late Campanian; Olivero and Medina Reference Olivero and Medina2000; Olivero Reference Olivero2012a,Reference Oliverob). These are the same region and levels in which Antarctopelta had been originally found (Salgado and Gasparini Reference Salgado and Gasparini2006; Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019).
Specimen Description
The new specimens are fragmentary and isolated unmodified osteoderms. The superficial surface of CAV-A4 is abraded, partially exposing the core region. Specimen CAV-A5 is rounded and exhibits a longitudinal keel on the surface. Both cross sections are triangular (Figs. 2, 3). CAV-A10 has an abraded basal region (Fig. 4), and all specimens exhibit a similar microstructure, with the core filled by a compact bone interstitial to resorption cavities. These cavities are wide and accumulate close to the basal cortex. The transition from cortical bone to the core is gradual, characterized by resorption cavities within the cortex. The osteoderms also lack densely abundant secondary (Haversian) bone, with sparse and isolated secondary osteons in the core. Few resorption cavities and rich ISFBs fill the core, which is rich in compact bone (Figs. 2–4). The ISFBs in the core are more scattered (Figs. 2B–D, and 3B), whereas they are more organized close to the external/basal and marginal surfaces of the cortex (Figs. 2E–G, and 3C–G). Specimen CAV-A10 exhibits wider erosional cavities than those observed in the other CAVs (Fig. 4A). Cross-polarized light reveals that the fibers in both core and cortex are highly orthogonally organized (Figs. 2C–F, 3B–F, and 4D–F). The bone matrix is rich in flattened osteocyte lacunae, which are aligned with the main directions of the ISFBs (Fig. 4C–F). Additionally, osteocyte lacunae lack preserved canaliculi. The vascular canals are mostly reticular (Figs. 2G, and 3B–G), but in the core of CAV-A10, they are mostly longitudinal (Fig. 4E,F).
The cortical bone is compact and poorly vascularized. The cortical bone of the core is thicker than the superficial cortical bone (Table 1). Further, some erosional cavities exhibit irregular inner margins, indicating an active resorption process (Figs. 2H, and 3H). The basal cortex is mostly absent in both CAV-A5 and CAV-A4. In cross section, the basal region is slightly abraded and does not show the decrease in resorption cavities observed in the superficial cortex. The structural fibers are more numerous and more organized in the basal region (Fig. 3C,D). Despite the core potentially being reduced in thickness due to weathering of the osteoderms, the alignment of the structural fibers in the basal region indicates that it was close to the margin—such well-organized structural fibers are characteristic of regions close to the cortex. If not absent, the basal cortex was at least relatively very thin compared with other ankylosaurian osteoderms.
Bone Origin, Growth, and Ontogeny
The microstructure of the CAV specimens resembles metaplastic bone features in osteoderms, as pointed out by Main et al. (2005: p. 303), characterized by “poorly vascularized, amorphous dermal bone tissue in which fibers are very numerous and oriented in many directions.” Another feature is the lack of canaliculi in the osteocyte lacunae (Levrat-Calviac and Zylberberg Reference Levrat-Calviac and Zylberberg1986). Different from neoplastic bones, metaplastic bone lacks the development of a periosteum, being formed from dense connective tissues, such as ligaments and tendons (Haines and Mohuiddin Reference Haines and Mohuiddin1968; Organ and Adams Reference Organ and Adams2005) and—in the case of ankylosaurian and archosaur osteoderms—from stratum compactum of the dermis (Scheyer and Sander Reference Scheyer and Sander2004; Main et al. Reference Main, de Ricqlès, Horner and Padian2005; Vickaryous and Sire Reference Vickaryous and Sire2009; Cerda and Powell Reference Cerda and Powell2010; Scheyer et al. Reference Scheyer, Desojo and Cerda2014; Cerda et al. Reference Cerda, García, Powell and Lopez2015, Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019; Bellardini and Cerda Reference Bellardini and Cerda2017; Ponce et al. Reference Ponce, Cerda, Desojo and Nesbitt2017). Although ossicles and osteoderms are both types of metaplastic bone, CAV specimens can be diagnosed as osteoderms due to the higher degree of vascularization in comparison with that observed in ossicles, as well as the structural fiber bundles being less organized. In ossicles, structural fibers form a three-dimensional mesh with bundles perpendicular to each other in all directions (e.g., de Ricqlès et al. Reference Ricqlès, Pereda Suberbiola, Gasparini and Olivero2001; Cerda and Powell Reference Cerda and Powell2010; Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019).
The CAV specimens lack growth marks and dense secondary osteons (only CAV-A10 shows some isolated secondary osteons). Some ankylosaurian osteoderms are devoid of such features (Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010), whereas others have been previously described as having primary lamellar bone tissue with lines of arrested growth (LAGs) and extensive secondary reconstruction (Scheyer and Sander Reference Scheyer and Sander2004; Main et al. Reference Main, de Ricqlès, Horner and Padian2005). Such a difference was regarded as representing distinct ontogenetic stages, with juvenile/subadult stages being characterized by the rapid growth of the bone—woven or ISFB—and further substituted in the adult stage by lamellar bone with LAGs and secondary bone through erosion/reconstruction processes (Main et al. Reference Main, de Ricqlès, Horner and Padian2005; Hayashi et al. Reference Hayashi, Carpenter and Suzuki2009, Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010). Considering the histological ontogenetic stages proposed for Stegosaurus (see Hayashi et al. Reference Hayashi, Carpenter and Suzuki2009), CAV-A4 and CAV-A5 are assigned to stage 1, characterizing juvenile/subadult individuals. CAV-A10 exhibits some secondary osteons and longitudinal vascular canals, which makes it assignable to stage 2, a subadult individual. In Antarctopelta—an adult close to fully grown (Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019)—the erosion/reconstruction process is advanced, with some LAGs. This is markedly different from the tissues observed among CAVs. On the other hand, armadillos and crocodilians show an asynchronous development of the osteoderms along the anteroposterior regions of the body (Vickaryous and Hall Reference Vickaryous and Hall2006, Reference Vickaryous and Hall2008), and a similar pattern has already been observed in dinosaurs among juveniles of Pinacosaurus (Burns et al. Reference Burns, Currie, Sissons and Arbour2011; Burns and Currie Reference Burns and Currie2014). Therefore, this ontogenetic inference is only tentative, and we urge caution in inferring an ontogenetic stage to the whole organism based solely on isolated osteoderms, due to the different timing in osteoderm development in relation to the whole organism (Main et al. Reference Main, de Ricqlès, Horner and Padian2005; Hayashi et al. Reference Hayashi, Carpenter and Suzuki2009, Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010) and the lack of detailed studies in the variation along the anteroposterior axis of ankylosaurians.
Microstructural Comparisons
Osteoderms rich in ISFBs are common to dinosaurs (titanosaurs and thyreophorans; Salgado Reference Salgado2003; Scheyer and Sander Reference Scheyer and Sander2004; Cerda and Powell Reference Cerda and Powell2010; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014; Chinsamy et al. Reference Chinsamy, Cerda and Powell2016). The core of the CAV specimens is filled with a few cavities, as similarly observed in other ankylosaurians (Burns and Currie Reference Burns and Currie2014) and in the titanosaurs MCS-Pv 181 and MCS-Pv 182 (Cerda et al. Reference Cerda, García, Powell and Lopez2015). However, the bone matrix of all CAVs is organized in bundles with an alignment of the osteocyte lacunae. Such a feature more closely resembles the pattern found in ankylosaurians (see Burns and Currie Reference Burns and Currie2014) than in titanosaurs, which do not show a clear bundle pattern (osteocyte lacunae are less organized) and are richer in woven bone (see Cerda and Powell Reference Cerda and Powell2010; Cerda et al. Reference Cerda, García, Powell and Lopez2015; Chinsamy et al. Reference Chinsamy, Cerda and Powell2016; Maidment et al. Reference Maidment, Strachan, Ouarhache, Scheyer, Brown, Fernandez, Johanson, Raven and Barrett2021).
The compact bone in the core is present in both late-diverging ankylosaurids and nodosaurids (Burns and Currie Reference Burns and Currie2014). The few resorption cavities and secondary osteons differ from the ones observed in the early-diverging nodosaurids Gastonia (DMNH 49754-1) and Gargoyleosaurus (DMNH 27726) and in ankylosaurids, in which this region of the osteoderm is marked by large cavities and remodeled trabecular bone (Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Sullivan Reference Burns and Sullivan2011; Burns and Currie Reference Burns and Currie2014). Osteoderm thickness is like that of nodosaurids (Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014), just like the orthogonal arrangement of the IFSBs (Scheyer and Sander Reference Scheyer and Sander2004; Burns and Currie Reference Burns and Currie2014). The thin (or absent) basal cortex is regarded as synapomorphic to nodosaurids (Scheyer and Sander Reference Scheyer and Sander2004; Burns and Currie Reference Burns and Currie2014).
Among nodosaurids, the shape of the cross section in all CAVs differs from that observed in the thoracic osteoderm of Edmontonia (TMP 1998.98.1) and from osteoderms of Glyptodonpelta (SMP VP-1580; Burns and Currie Reference Burns and Currie2014), in which the superficial and basal cortices are aligned. However, the microstructure of CAVs is like that of Edmontonia and Glyptodonpelta, in that there are ISFBs, no clear stratification between the cortex and core, and few scattered secondary osteons. The CAV specimens differ from Sauropelta (DMNH 18206; Burns and Currie Reference Burns and Currie2014) in the stratification between core and cortex, marked by secondary osteons. Although the shape of the CAVs is quite like the dermal plate of Antarctopelta (Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019), it differs from the latter by the presence of trabeculae composed of secondary lamellar bone in Antarctopelta.
Dinosaurian Osteoderm Thickness Patterns and Linear Model
The OLS indicates a strong positive isometric relationship between Co and T among dinosaurian osteoderms (Fig. 5A)—see “Morphometric Dataset” for morphometric variable abbreviations. Among the relative thickness indices, the SI exhibits a positive isometric relationship with CI and a positive allometry relationship with CoI. The strongest correlation was observed between Co and T (R 2 = 0.962, p = 4.235 × 10−19; Fig. 5A), of which its linear function is:
$$\log {\rm Co} = 1.169\log T-0.304$$Ordinary least squares (OLS) based on core osteoderm (Co), total thickness of osteoderm (T), cortical index (CI), core index (CoI), superficial index (SI), and basal index (BI) among dinosaurian osteoderms. A, The correlation between Co and T is strong in general and between all groups. B, Although osteoderms, in general, have a clear correlation between CI and SI, only nodosaurids have a clear correlation in comparison to ankylosaurids. C, All the osteoderm groups have a moderate correlation between CoI and SI. Only the osteoderms of nodosaurids have a moderate correlation with BI and CoI (D) and BI and CI (E). F, The relationship between BI and SI has no clear correlation. The silhouette of the osteoderm cross section in the graphs indicates the osteoderm shape deformation along the main axes.
This function provides a model to estimate the mean of the core thickness or the mean of the total thickness in cases in which the cortical bones are only partially preserved.
Osteoderm Thickness Variations between Groups
We detect a linear correlation between CI and SI when considering all specimens in general and in the subsample, including nodosaurids only (R 2 = 0.613, p = 7.391 × 10−8; and R 2 = 0.748, p = 3.593 × 10−6, respectively; Fig. 5B), but not when considering Ankylosauridae only (R 2 = 0.271, p = 0.186). The CoI and SI are slightly correlated between all groups (R 2 = 0.586, p = 2.069 × 10−7 for all specimens; R 2 = 0.578, p = 0.029 for ankylosaurids; and R 2 = 0.575, p = 0.003 × 10−1 for nodosaurids; Fig. 5C). However, there is no clear correlation between BI and the other parameters (considering all specimens: R 2 = 0.271, p = 0.002 for CoI; R 2 = 0.329, p = 0.005 × 10−1 for CI; and R 2 = 0.005, p = 0.685 for SI; Fig. 5D–F). The only exceptions are BI × CoI and BI × CI in nodosaurids (R 2 = 0.625, p = 9.407 × 10−5; and R 2 = 0.604, p = 0.001 × 10−1, respectively; Fig. 5D,E). Therefore, the most informative parameters of covariation among titanosaurs and ankylosaurians comprise SI, CI, and CoI. The OLS indicates a slight shift between ankylosaurids and nodosaurids when contrasting SI with both CI and CoI—especially SI and CI, despite the lack of correlation between these two parameters among ankylosaurid osteoderms (R 2 = 0.271, p = 0.186; Fig. 5B).
The univariate KW analysis indicates no clear difference between medians among taxonomic groups under any of the four parameters tested here (Table 2). The box plot indicates a slight difference between ankylosaurids and nodosaurids under BI and CoI, but there is a large overlap in the range between these groups (Fig. 6). Therefore, osteoderm morphometric data are of very limited value to distinguish between ankylosaurian groups given the most recent systematic classification of this group (Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; see “Ankylosaur Phylogeny and Assignment of CAV Specimens”).
Box plots of dinosaurian osteoderm shape based on basal index (BI), cortical index (CI), core index (CoI), and superficial index (SI) ratios. The dinosaurian groups comprise ankylosaurids, nodosaurids (including the new CAV specimens), titanosaurs, and the parankylosaurian Antarctopelta. There is no significant difference between group medians; see Table 2 for statistical tests. The colors and symbols are the same as in the legend in Fig. 5.
Univariate Kruskal-Wallis test considering ankylosaurids, nodosaurids, and titanosaurs. No value was significant (p < 0.05). Abbreviations: BI, basal index; CI, cortical index; CoI, core index; SI, superficial index.
The relative thickness of the cortex (CI) decreases exponentially in comparison to the total thickness (T) of the osteoderm among ankylosaurians, being the thinnest among nodosaurids (Fig. 7), as initially suggested by Burns and Currie (Reference Burns and Currie2014). Considering T as a proxy for osteoderm size, such a relationship suggests an allometric decrease. Our comparisons between our EBR and that observed in Burns and Currie (Reference Burns and Currie2014) show higher AIC indices than in the latter. Ankylosaurids were the group that best fits the EBR model (AIC = 653.71; Fig. 7B). The inclusion of “polacanthids” into nodosaurids could contributed to the increase of AIC compared with the models in Burns and Currie (Reference Burns and Currie2014).
Exponential bivariate regression (EBR) between cortical index (CI) and total osteoderm (T), comparing the results of Burns and Currie (Reference Burns and Currie2014: fig. 8) (A) with our results (B). In both analyses, nodosaurids exhibit the lowest Akaike information criterion (AIC), indicating the best fit to the EBR. High CI marks early-diverging ankylosaurians, whereas lower CI marks nodosaurids and titanosaurs. Note that nodosaurids and titanosaurs exhibit higher T and lower CI. Names in bold were originally analyzed by Burns and Currie (Reference Burns and Currie2014) and in this study. The colors and symbols in B are the same as in the legend in Fig. 5.
The lowest CI is observed among most nodosaurids and titanosaurs (Fig. 7), with CAV-A4 exhibiting the lowest CI and T compared with all other osteoderms analyzed. The nodosaurids Edmontonia, Sauropelta, TMP 1967.10.29, and the titanosaur Saltasaurus exhibit the highest T values, although their CI values are nearly the same as most other nodosaurids. The CI overlaps between all ankylosaurian groups and indicates that the assessment of ankylosaurids as “thin-walled” osteoderms is subjective (Burns Reference Burns2008; Burns and Currie Reference Burns and Currie2014). The expression “thin-walled” is also vague and could also refer to low CI values rather than T, but our assessments do not indicate low CI for ankylosaurids. Instead, we find that low CI is common among nodosaurids and ankylosaurids, with later-diverging nodosaurids having an even lower CI than ankylosaurids.
Morphospace Occupation of Dinosaurian Osteoderms
We performed the PCA and LDA with the indices that have a significant linear relationship—CI, CoI, and SI (Fig. 8, Supplementary Fig. 3). SI exhibits a positive linear relationship with CI and a weak negative linear relationship with CoI (Figs. 5B,C, and 8). SI is the main predictor of group membership (LD 1; Supplementary Fig. 3); although the LDA indicates that it could not strongly discriminate between groups (58.62% of certainty). The OLS and PCA reveal that most of the data variation is explained by changes in relative cortical thickness (CI) mostly on the superficial cortex (SI) (Fig. 5; see also PC 1 in Fig. 8). Relative core thickness (CoI) is the main factor responsible for the overlap between ankylosaurids, nodosaurids, and titanosaurs.
Result of the PCA. The silhouette of the osteoderm cross section within graphs indicates the osteoderm shape deformation along the main axes. Both superficial index (SI) and cortical index (CI) contribute to most of the variation in the morphospace observed among ankylosaurian osteoderms. Abbreviations: Coa, circular off-apex osteoderm; F, flat osteoderm; Fk, flat keeled osteoderm; k, keeled osteoderm, but unknown shape; oa, off-apex osteoderm, but unknown shape; Ooa, oval off-apex osteoderm; Ok, oval keeled osteoderm; Sp, spine. The colors and symbols are the same as in the legend for Fig. 5.
The morphospace of ankylosaurids considerably overlaps with the morphospace of nodosaurids (Fig. 8). However, nodosaurids have a larger variance of osteoderm shape along the main axes of variation (PC 1), which reflects the enormous observed variation in cortical thickness in this group (e.g., CoI) and accounts for the vast majority of the total variation in the data (88.7%). The somewhat broader morphospace of ankylosaurids along PC 2 is mostly explained by the great degree of variation in CI and SI introduced by Euoplocephalus and Pinacosaurus (Fig. 8, Supplementary Fig. 3). Additionally, two of the indeterminate ankylosaurid specimens (TMP 1987.113.4A.1 and TMP 1998.98.1A.2) are outliers, located in a much more distant region of the morphospace relative to all other osteoderms that could be confidently identified to the family level.
Notably, the larger morphospace occupation of nodosaurids along the first major axis (PC 1) is highly influenced by the inclusion of the new materials described here from Antarctica (most notably, CAV-4). The latter has expanded the nodosaurid morphospace into a region previously unoccupied by any other group of dinosaurs. This suggests a unique osteoderm morphology for the new Antarctic specimens compared with other nodosaurids, ankylosaurians, and dinosaurs. The large variation in CoI among nodosaurids is responsible for the unique area of the morphospace occupied by this group, including the new CAV specimens. This suggests a reduction in cortical thickness on the osteoderms of some nodosaurids compared with other ankylosaurians, which have retained much thicker cortical bone.
Titanosaurs occupy a much smaller morphospace than all three groups of ankylosaurians, also partially overlapping with nodosaurids. However, we note that we have few data points (only three) for titanosaurs; the morphospace occupation by the group could much broader.
Ankylosaur Phylogeny and Assignment of CAV Specimens
Both phylogenetic and paleohistological data indicate a division of ankylosaurians into two main groups: Ankylosauridae and Nodosauridae (Hill et al. Reference Hill, Witmer and Norell2003; Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014; Arbour and Currie Reference Arbour and Currie2016; Arbour et al. Reference Arbour, Zanno and Gates2016; Brown et al. Reference Brown, Henderson, Vinther, Fletcher, Sistiaga, Herrera and Summons2017; Rivera-Sylva et al. Reference Rivera-Sylva, Frey, Stinnesbeck, Carbot-Chanona, Sanchez-Uribe and Guzmán-Gutiérrez2018; Fig. 9A,B). The group “Polacanthidae” (or “Polacanthinae”) was recovered as paraphyletic here, as in previous phylogenies (Hill et al. Reference Hill, Witmer and Norell2003; Arbour and Currie Reference Arbour and Currie2016; Arbour et al. Reference Arbour, Zanno and Gates2016; Brown et al. Reference Brown, Henderson, Vinther, Fletcher, Sistiaga, Herrera and Summons2017; Rivera-Sylva et al. Reference Rivera-Sylva, Frey, Stinnesbeck, Carbot-Chanona, Sanchez-Uribe and Guzmán-Gutiérrez2018; Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; Frauenfelder et al. Reference Frauenfelder, Bell, Brougham, Bevitt, Bicknell, Kear, Wroe and Campione2022). Most “polacanthids” are inferred as early-diverging nodosaurids (e.g., Polacanthus, Gastonia, and Mymooropelta; Fig. 9). Regarding microstructural differences, ankylosaurids have been described as possessing thin-walled osteoderms, with extensive Haversian bone and perpendicular organized structural fibers; nodosaurids exhibit thin/absent basal cortex and highly ordered orthogonal structural fibers; while “polachantids” have a similar thickness with ankylosaurid osteoderms, but the structural fibers are more diffuse (Scheyer and Sander Reference Scheyer and Sander2004; Burns and Currie Reference Burns and Currie2014). Our results indicate that these osteoderm features are widespread among ankylosaurians, and only the presence of Haversian bone and lamellar bone in the superficial cortex distinguish nodosaurids from other ankylosaurians. Despite the presence of few osteoderm phylogenetic characters and taxa scored in the matrix, these were sufficient to recover the family-level placement of the new (CAV) specimens from Antarctica (Fig. 9B).
Phylogenetic results from the dataset of Soto-Acuña et al. (Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021). Strict consensus of the three most parsimonious trees (MPTs) of 695 steps, without CAV specimens in (A). Strict consensus of seven MPTs of 697 steps each including CAV specimens (B). Osteoderm characters mapped on the trees: Ch. 157: external cortical histology of skeletally mature osteoderms: no osteoderms (0) lamellar bone (1), ISFB (2); Ch. 158: Haversian bone in osteoderms: no osteoderms (0) absent in core of skeletally mature osteoderms (1), maybe present in the core of skeletally mature osteoderms (2); Ch. 159: basal cortex of skeletally mature osteoderms: no osteoderms (0) present (1), absent or poorly developed (2); Ch. 160: structural fiber arrangement in osteoderms: no osteoderms (0) structural fibers absent (1), reaches orthogonal arrangement near osteoderm surfaces (2), diffuse throughout (3), highly ordered sets of orthogonally arranged fibers in the superficial cortex (4). See Brum et al. (Reference Brum, Eleutério, Simões, Whitney, Souza, Sayão and Kellner2023) for a detailed common synapomorphy list.
Discussion
Osteoderm Microstructure and Evolution
Nodosaurids and titanosaurs have much thinner osteoderm cortices compared with early-diverging ankylosaurids (Figs. 5–8). This could be associated with structural constraints, such as the limitation of nutrient circulation in osteoderms due to the lack of canaliculi (Haines and Mohuiddin Reference Haines and Mohuiddin1968). A few titanosaur osteoderms exhibit wide vascular canals and a network of cavities in their osteoderms (Cerda and Powell Reference Cerda and Powell2010; Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011; Cerda et al. Reference Cerda, García, Powell and Lopez2015), but this is not the case among the titanosaurs sampled here, whose core structure is much more like that found in ankylosaurians. Although some previous studies indicated the presence of vascular pits in ankylosaurian osteoderms, this is restricted to only a few specimens (Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014). Therefore, late-diverging nodosaurids, and the sampled titanosaurs, may have compensated for the nutritional circulatory limitations imposed by the increase in core volume (Co; Fig. 5A) and total osteoderm thickness (T) by reducing their superficial and basal cortical thickness. The convergent evolution of cortical thickness reduction in these distantly related lineages suggests that the decrease of core density is a physiological adaptation for coping with large absolute osteodermal size.
However, we caution against establishing a direct association between cortical thickness and core density. CAV materials have a thin cortex and compact cores (Figs. 2–4, 7), whereas other nodosaurid and titanosaur osteoderms have a thinner cortex and less dense core compared with the osteoderms of other ankylosaurians (e.g., Gastonia, Euoplocephalus, Nodocephalosaurus, and Saichania). Although the physiological adaptation is a strong constraint to the osteoderm size and cortical and core thicknesses, other factors could be contributing to this balance (e.g., osteoderm shape, biomechanics, variation in developmental timing along the body axis).
Additionally, we observe that “polacanthid” Gastonia, Mymoorapelta, Polacanthus, Gargoyleosaurus, and Ahshislepelta (now considered early-diverging nodosaurids; Fig. 9) mainly occupy a region of morphospace that overlaps with some ankylosaurids (Fig. 8). This partially reflects their phylogenetic instability (Thompson et al. Reference Thompson, Parish, Maidment and Barrett2012; Arbour and Currie Reference Arbour and Currie2016; Arbour et al. Reference Arbour, Zanno and Gates2016; Brown et al. Reference Brown, Henderson, Vinther, Fletcher, Sistiaga, Herrera and Summons2017; Rivera-Sylva et al. Reference Rivera-Sylva, Frey, Stinnesbeck, Carbot-Chanona, Sanchez-Uribe and Guzmán-Gutiérrez2018) and supports previous claims that there is no clear distinction between ankylosaurid and “polacanthid” osteoderms (Burns and Currie Reference Burns and Currie2014). Our results suggest that the region of morphospace occupied by both early-diverging nodosaurids (“polacanthids”) and ankylosaurids—high thickness of compact bone layers, especially the superficial cortex—represents the ancestral morphospace in the early evolution of ankylosaurians. These hypotheses would be supported even in the former classification of “polacanthid” as early-diverging ankylosaurians (outside ankylosaurids and nodosaurids; Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Burns and Currie Reference Burns and Currie2014).
Our phylogenetic results agree with recent analyses (Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021) in recovering Antarctopelta as closely related to Stegouros and Kunbarrasaurus, forming the new early-diverging ankylosaur group Parankylosauria, and also reveal the historically phylogenetic unstable “polacanthids” to be early-diverging nodosaurids (Soto-Acuña et al. Reference Soto-Acuña, Vargas, Kaluza, Leppe, Botelho, Palma-Liberona, Simon-Gutstein, Fernández, Ortiz, Milla, Aravena, Manríquez, Alarcón-Muñoz, Pino, Trevisan, Mansilla, Hinojosa, Muñoz-Walther and Rubilar-Rogers2021; Fig. 9A,B). However, our results indicate that the osteoderm characters of Antarctopelta resemble those observed in late-diverging nodosaurids and in CAV specimens (Fig. 9B).
The main difference between ankylosaurids and late-diverging nodosaurids is the basal cortex (BI; as in Burns and Currie Reference Burns and Currie2014) and the thickness of the core (CoI; not the occurrence of cavities), as discussed later. The new specimens described here (CAVs) were phylogenetically inferred to be nodosaurids, and their cortex thickness indicates they have some of the thinnest cortices among nodosaurids, even if the cortex may have been partially lost due to weathering—see “Specimen Description.” We also noticed the presence of highly ordered orthogonally arranged fibers (Fig. 9B) and the wide core thickness (Fig. 8) in CAV specimens and Antarctopelta. This pattern and the phylogenetic analysis suggest that the CAV specimens are from late-diverging nodosaurids. Considering the recent suggestion of Antarctopelta as an early-diverging ankylosaurian, our results suggest a larger phylogenetic diversity of Antarctic ankylosaurians, including both early- and late-diverging forms.
The Function of Body Osteoderms
The several differences observed in relative thickness and osteoderm microstructure among various groups of dinosaurs (Figs. 5–8) have led to several hypotheses surrounding their inferred function, including body armor, display, and thermoregulation (Scheyer and Sander Reference Scheyer and Sander2004; Cerda and Powell Reference Cerda and Powell2010; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010; Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011; Burns and Currie Reference Burns and Currie2014; Cerda et al. Reference Cerda, García, Powell and Lopez2015; Brown Reference Brown2017; Vidal et al. Reference Vidal, Ortega, Gascó, Serrano-Martínez and Sanz2017). The osteoderms analyzed here exhibit different morphologies (e.g., spikes, oval keeled, circular with offset apex; Fig. 8; see also Supplementary Table 1), indicating they may come from different body regions, excluding tail clubs and ossicles. The fragmentary preservation of the osteoderms prevents us from a detailed assessment of osteoderm microstructure and function across individual body regions, but the available data still provide informative insights into the function of ankylosaurian body osteoderms in general.
The association between microstructure and morphology in osteoderms enables the interpretation of multiple adaptive functions, as observed among extant reptile groups. For instance, in extant crocodilians, Sharpey's fibers in the cortex are perpendicular to the external margin of the osteoderm, being anchored by tendon attachment—as also inferred for stegosaurids (Scheyer and Sander Reference Scheyer and Sander2004). The presence of these fibers in crocodylomorphs contributes to the flexibility of the whole armor and improves its strength (Sun and Chen Reference Sun and Chen2013). In addition, keeled body osteoderms of crocodylomorphs are more resistant to breakage relative to non-keeled forms, whereas stress tends to concentrate near pits and ridges in ornamented osteoderms, which represent areas of crack initiation and propagation (Clarac et al. Reference Clarac, Goussard, De Buffrénil and Sansalone2019). However, the porosity of the core and the 3D vascular network are factors that enhance the bending stiffness and pressure absorption, as well as representing a trade-off between resistance and the physiological balance in crocodylomorphs—all of which may change throughout their ontogeny (Sun and Chen Reference Sun and Chen2013; Chen et al. Reference Chen, Yang and Meyers2014; Clarac et al. Reference Clarac, Goussard, De Buffrénil and Sansalone2019). In cordylid lizards, porous core osteoderms are fractured under relatively low-stress pressure and exhibit lower thermal conductivity. On the other hand, compact core osteoderms are comparatively more resistant to higher pressures and exhibit higher thermal conductivity (Broeckhoven et al. Reference Broeckhoven, du Plessis and Hui2017). Therefore, in both extant crocodiles and cordylid lizards, osteoderms have multiple functions with distinct physiological trade-offs (see Gould and Vrba Reference Gould and Vrba1982; Gould Reference Gould1991; Buss et al. Reference Buss, Haselton, Shackelford, Bleske and Wakefield1998).
Ankylosaurians and titanosaurs differ from crocodilians and stegosaurids in their extensive distribution of structural fibers, which indicates that osteoderms were fully surrounded by dermis (Scheyer and Sander Reference Scheyer and Sander2004; Cerda and Powell Reference Cerda and Powell2010; Cerda et al. Reference Cerda, García, Powell and Lopez2015). In small osteoderms of both nodosaurids and ankylosaurids, the abundance of structural fibers indicates efficient lightweight body armors (Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010), but their arrangement within late-diverging nodosaurids is distinct in comparison to the other ankylosaurians. The highly ordered orthogonal sets of structural fibers in late-diverging nodosaurids (Burns and Currie Reference Burns and Currie2014) suggest more efficiency in pressure dissipation on the dermis due to their connectivity to ossicles and adjacent osteoderms (de Ricqlès et al. Reference Ricqlès, Pereda Suberbiola, Gasparini and Olivero2001; Scheyer and Sander Reference Scheyer and Sander2004). Considering their external morphology, the osteoderms of ankylosaurians and titanosaurs are keeled and are devoid of extensive superficial ornamentation, thus suggesting high-stress resistance in these groups based on data from extant crocodiles (Clarac et al. Reference Clarac, Goussard, De Buffrénil and Sansalone2019). In ankylosaurians, the morphology of small body osteoderms (e.g., rounded keeled) resembles the plesiomorphic condition found in the early-diverging thyreophoran Scutellosaurus, which was retained (or even augmented) in ankylosaurians (Main et al. Reference Main, de Ricqlès, Horner and Padian2005). The acquisition of abundant structural fibers among small osteoderms of ankylosaurians (Scheyer and Sander Reference Scheyer and Sander2004; Hayashi et al. Reference Hayashi, Carpenter, Scheyer, Watabe and Suzuki2010) further suggests the body-armoring function of these osteoderms. Accordingly, the erosional cavities with irregular margins in CAVs (Figs. 2H, and 3H) indicate an active resorption process, which could reflect a response to microstructural fractures resulting from such a defensive role (see Robling et al. Reference Robling, Castillo and Turner2006).
In titanosaurs, the gross microstructure resembles that observed in nodosaurids—the high CoI is also accompanied by high vascularization (e.g., Cerda and Powell Reference Cerda and Powell2010; Cerda et al. Reference Cerda, García, Powell and Lopez2015). However, the vascularization pattern is distinct from that of nodosaurids, as titanosaurs have wide chambers and a main vascular canal running along the main axis of the osteoderm (Cerda and Powell Reference Cerda and Powell2010; Curry Rogers et al. Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011; Cerda et al. Reference Cerda, García, Powell and Lopez2015; Vidal et al. Reference Vidal, Ortega, Gascó, Serrano-Martínez and Sanz2017). Such vascularization in titanosaurs has led to the hypothesis that physiological balance/calcium remobilization is the primary role of osteoderms in this group (Curry-Rogers Reference Curry Rogers, D'Emic, Rogers, Vickaryous and Cagan2011; Vidal et al. Reference Vidal, Ortega, Gascó, Serrano-Martínez and Sanz2017). We argue, however, that this morphology in titanosaurs is consistent with the trade-off between protection/calcium balance found in crocodiles throughout their ontogeny (Sun and Chen Reference Sun and Chen2013; Chen et al. Reference Chen, Yang and Meyers2014; Clarac et al. Reference Clarac, Goussard, De Buffrénil and Sansalone2019). Therefore, the similar degree of high vascularization and higher SI found in nodosaurids, suggests that, besides a role in body armoring, their body osteoderms would also have an important physiological balance role in calcium remobilization. This later physiological role is clear among the new CAV specimens from Antarctica and published images of other ankylosaurian osteoderms. High vascularization, resorption cavities, and high SI support the interpretation of the role of body osteoderms in calcium balance (Cerda et al. Reference Cerda, Gasparini, Coria, Salgado, Reguero, Ponce, Gonzalez, Jannello and Moly2019: fig. 2A). This extremely important but previously overlooked physiological function of ankylosaurian osteoderms may well have played a key role as a preadaptation to extreme environmental conditions, such as those represented by the low light conditions and lower temperatures during the Cretaceous Antarctic winter.
Conclusions
Understanding the function and evolution of the unique body armor of ankylosaurian dinosaurs is still in its infancy and is highly limited by phylogenetic instability and low sample sizes from geographic regions in the Southern Hemisphere. We shed some light on those topics by describing new specimens recovered from Antarctica, combining their data with previous data, and assessing the evolution of osteoderms across morphospace and evolutionary history based on an updated phylogenetic hypothesis.
Our results indicate that it is possible to use osteoderms to differentiate between ankylosaurid and nodosaurid osteoderms based on their microstructures. However, we could not differentiate ankylosaurian osteoderms from those of titanosaurs. The newly discovered Antarctic material (CAV specimens) share a microstructural pattern with late-diverging nodosaurids and are recovered phylogenetically within this group. Additionally, all CAVs share the highly ordered arrangement of the structural fibers in the superficial cortex with the parankylosaurian Antarctopelta and the late-diverging nodosaurid Sauropelta (Fig. 9B). Finally, CAVs possess a histological microstructure pattern compatible with an early ontogenetic stage.
We provide quantitative support for previous hypotheses that ankylosaurids share the microstructure and the morphospace of early-diverging nodosaurids—previously recognized as “polacanthids” (Burns and Currie Reference Burns and Currie2014). Late-diverging nodosaurids are very distinct from all other ankylosaurians by their thin or absent basal cortex and thick core. Such osteoderm structure is similar to that observed in the osteoderms of titanosaurs that lack wide chambers and/or wide vascular networks, (e.g., Saltasaurus). As previously hypothesized, the primary function of osteoderms of both ankylosaurids and nodosaurids could be to act body armor. Additionally, we provide histological evidence combined with information from the literature that nodosaurids underwent a rapid early growth strategy combined with an increase of osteoderm vascularization, suggesting some role in calcium remobilization for physiological balance, as in titanosaur sauropods. Therefore, the function of body osteoderms in dinosaurs, especially in ankylosaurians, seems to be characterized by a more complex trade-off mechanism between biomechanical and physiological functions than simplistic explanations related to body protection or display, as previously thought, and the growth pattern in long bones of Antarctopelta indicate they were preadapted to the colonization of higher-latitude environments.
Acknowledgments
This project was supported by Programa Antártico Brasileiro (PROANTAR) through the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq: grant nos. 420687/2016-5 and 313461/2018-0 to A.W.A.K.; no. 314222/2020-0 to J.M.S.) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES: PROANTAR fellowship no. 88887.336584/2019-00 to A.S.B.). A.W.A.K. acknowledges funding from Fundação de Desenvolvimento Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ no. E-26/202.905/2018) and T.R.S. acknowledges a postdoctoral fellowship provided by the Natural Sciences and Engineering Research Council of Canada (NSERC). We thank T. M. Scheyer, T. G. Frauenfelder, an anonymous reviewer, and N. Campione for comments that helped improve the quality of the article. We acknowledge the Willi Henning Society for the free availability of the software TNT. The team of the PALEOANTAR Project thanks the entire NApOc Ary Rongel military group and the pilots of the HU-1 helicopter squadron for their logistic support during fieldwork on the Antarctic Peninsula.
Declaration of Competing Interests
The authors declare no competing interests.
Data Availability Statement
Data available from the Dryad Digital Repository: https://doi.org/10.5061/dryad.ht76hdrkx.
Open access
