A new ankylosaurid from the Upper Cretaceous Nemegt Formation of Mongolia and implications for paleoecology of armoured dinosaurs

A new ankylosaurid dinosaur, Tarchia tumanovae sp. nov., has been recovered from the Upper Cretaceous Nemegt Formation of Mongolia. It includes a well-preserved skull, dorsal, sacral, caudal vertebrae, sixteen dorsal ribs, ilia, a partial ischium, free osteoderms, and a tail club. The squamosal horns of T. tumanovae are divided into two layers, the external dermal layer and the underlying squamosal horn proper. The irregular ventral margin of the base of the upper dermal layer may represent a resorption surface, suggesting that the squamosal horns of some ankylosaurids underwent extreme ontogenetic remodeling. Localized pathologies on the dorsosacral ribs and the tail provide evidence of agonistic behaviour. The tail club knob asymmetry of T. tumanovae resulted from restricted bone growth due to tail club strikes. Furthermore, T. tumanovae had an anteriorly protruded shovel-shaped beak, which is a morphological character of selective feeders. Ankylosaurid diets shifted from low-level bulk feeding to selective feeding during the Baruungoyot and the Nemegt “age” (middle Campanian-lower Maastrichtian). This ankylosaurid niche shifting might have been a response to habitat change and competition with other bulk-feeding herbivores.

(e-f) Skeletal diagram of the specimen in dorsal (e) and left lateral (f) views (white bones represent recovered elements). (g) Skeletal reconstruction with dermal armour. Abbreviations: 7os, Type 7 osteoderm; csr, caudosacral vertebra; dr, dorsal rib; dsr, dorsosacral vertebra; il, ilium; maj os, major osteoderm; ot, ossified tendon; sk, skull. Adobe Illustrator CC (version 24.0.1, https:// www. adobe. com/ kr/ produ cts/ illus trator. html) was employed to produce (a-g).   www.nature.com/scientificreports/ Description. The shape of the skull is trapezoidal and broader than long in dorsal view (Fig. 3, see Supplementary Information S1 for measurements). All caputegulae are pitted externally. The neuroanatomy of MPC-D 100/1353 was fully described before by Paulina-Carabajal et al. 9 . The premaxillae are fused dorsally, but the palatal surfaces are separate (Figs. 2, 3). The rostral tip is protruded anteroventrally, and a premaxillary notch is present. The narrow internarial bar is oriented posterodorsally. A thin, rugose ossification is present on each premaxilla below the external naris. The anterior margin of the premaxillary ornamentation is convex, whereas the posterior margin is concave. The palatal surface of the premaxillae is shovel-like with a round anterior boundary. No premaxillary teeth are present. The subcircular external nares face anteriorly. The entrance to the airway (aperture A, sensu 14 ) is large and subcircular. The airway is filled with a matrix. A supranarial notch is present on the medial wall of the maxillary region, lateral to the entrance to the airway. A short, medioventrally sharp intranasal process is present beneath the entrance of this aperture. A larger oval dorsolaterally-facing paranasal aperture (aperture C) is situated on the ventral wall of the external nares behind the internasal bar. The external nares are rimmed dorsomedially, dorsally, and laterally by supranarial caputegulae (sensu 15 ). The dorsomedial and dorsal portion of the supranarial caputegulum is thin, whereas the wide ventrolateral part gives the caputegulum a boot-like appearance in lateral view. The anterior end of the nasals contacts the internasal bar anteromedially. A single medium-sized bulbous internarial caputegulum is situated above the contact between the internarial bar of the premaxillae and the anterior nasals and between the thin medial portion of the two supranarial caputegulae. Behind the internarial caputegulum, eleven pyramidal nasal caputegulae are present. Most of these caputegulae are large and surrounded by a broad The maxillae are anteroposteriorly elongate, extending to below the orbits. In lateral view, the anterior portion of the maxillae is covered by the loreal caputegulum, whereas the posterior portion is exposed. In palatal view, the convex maxillary tooth row is situated medial to the buccal emargination. Nineteen alveoli are present in each maxilla. A single loreal caputegulum is present on each side. The loreal caputegulum is large, rhomboid, and has a posterior keel that posterolaterally juts out. The lacrimal caputegulum is large and flat laterally and has a keeled edge on the dorsal margin. A single medium-sized possible prefrontal caputegulum is present on each side, forming the lateral margin of the skull. The left possible prefrontal caputegulum is damaged, but the right is well preserved. This caputegulum is keeled dorsally and has a lateroposteriorly pointing apex in close contact with the anterior supraorbital caputegulum.
Paired, transversely oriented frontal caputegulae are pyramidal and rectangular (Fig. 3a,c). The nasofrontal sagittal furrow has a weak Z-shaped offset as in ZPAL MgD I/111 (holotype of T. kielanae). Although the parietal caputegulae are poorly defined, three shallow dorsal furrows diverge anteroposteriorly from each other at an angle of about 35°. Two medial supraorbital caputegulae are present on the right, whereas only the anteromedial one is preserved on the left. The longitudinally oriented anteromedial caputegulum is sigmoidal. The transversally oriented posterolateral caputegulum has a rounded anterior surface, a flat, anteriorly-inclined posterior surface. The left supraorbital caputegulae have damaged apices, whereas the right ones are well preserved with a distinct apex. The anterior supraorbital caputegulum is narrow and keeled, with an anterolaterally directed apex. The posterior supraorbital caputegulum is laterally pointed and keeled dorsally and is about five times larger than the anterior supraorbital caputegulum. The orbits have a posteriorly thick orbital rim and face anterolaterally. The quadratojugal horns are triangular and project ventrolaterally (Figs. 2, 3). A quadratojugal "neck" (sensu 16 ) is present at the base of the horns. Postocular caputegulae are absent. The squamosal horns are pyramidal and posteriorly recurved (Figs. 3a,c). These horns are divided into the upper external layer of the squamosal horn and the underlying squamosal horn proper (sensu 17 ). The former is dorsolaterally keeled with a longitudinal furrow present on the posterior half of the keel. The narrow, sharp, and medially curved anterior portion of the caputegulae lies in a broad, deep postorbital fossa posterior to the supraorbital. The anterior tip of the right horn was slightly broken sometime after CT scanning of the braincase by Paulina-Carabajal et al. 9 . A narrow, deep sulcus separates the irregular ventral margin of the base of the external layer of the squamosal horn and the underlying squamosal horn proper. The surface of the external dermal layer is pitted, whereas the squamosal horn proper has a granular texture. The nuchal shelf is dorsally uplifted and does not overhang the outer rim of the skull. Two nuchal caputegulae are present on each side. Both anterior and lateral nuchal caputegulae are elongate and transversely positioned. However, the lateral nuchal caputegulae are about four times larger than www.nature.com/scientificreports/ the anterior nuchal caputegulae, forming the posterior margin of the cranium. In occipital view, the nuchal shelf is not fused with the supraoccipital and paroccipital processes (Fig. 2d,h). The rostral extension of the vomer is dorsoventrally thin, splayed, and fused with the posteromedial region of the premaxillae (Fig. 3b,d). The osseous nasal septum (sensu 18 ) extends dorsally but does not meet the skull roof. The vomerine keel ends ventral to the alveolar ridge. The palatine extends posteromedially from the maxilla and gently projects dorsally, forming a posteroventral secondary palate (sensu 18 ). The ectopterygoid is small and wedge-like. The pterygoid has a vertical anterior surface with a foramen pierced through the central body. The pterygoid flange projects anterolaterally and contacts the dorsally positioned ectopterygoid. The quadrate ramus contacts the posterolaterally positioned quadrate. The posterolateral edge of the quadrate ramus is damaged on both sides. The posteromedial margin of the main pterygoid body is not fused with the basipterygoid processes of the basisphenoid. The basipterygoid processes are divided from each other. An interpterygoid vacuity is present between the paired pterygoids.
In palatal view, the contact between the basisphenoid and the basioccipital forms a rugose transverse ridge (Fig. 3b,d). A basioccipital foramen is present on the convex ventral surface of the basioccipital. The basioccipital and the exoccipitals are entirely fused and form a reniform occipital condyle oriented posteroventrally. The ovoid foramen magnum is taller than wide (Fig. 2d,h). A small hill-like process lies between the foramen magnum and the nuchal shelf. A horizontal groove is present below the paired exoccipital protuberances. The lateral terminus of the paroccipital process is long but not fused to the quadrate reaching laterally to the squamosal horns. The transversely broad quadrates are inclined anteroventrally toward the distal articular condyles in lateral view. The medial condyle of the quadrate is larger than the lateral condyle.
Only the third and eighth right maxillary teeth and the seventh left maxillary tooth are preserved (Fig. 3e-g). Although these teeth are partially embedded within the sockets, up to eleven marginal denticles can be observed. Shallow vertical grooves are present between the denticles, and a shelf-like labial cingulum is also present.
Two isolated dorsal vertebrae with fused ribs (Fig. 4), probably presenting the fourth and eleventh dorsal vertebrae based on the length of the ribs, are preserved. Only the right rib is preserved on the fourth dorsal www.nature.com/scientificreports/ vertebra (Fig. 4e,f), whereas both ribs are present on the eleventh (Fig. 4g,h,j).  The synsacrum is well preserved and includes three dorsosacrals, one parasacral, three sacrals, and two caudosacral vertebrae (Fig. 6). These vertebrae are fused along the centra and the neural spines. The centra of these vertebrae are spool-shaped and laterally constricted. A single medial ridge is present on the ventral surface of the centra of the dorsosacrals and the parasacral. Two medial ridges are present on the first and second sacrals and both caudosacrals. The neural arches of the dorsosacrals are centrally located dorsal to the centra, whereas the neural arches of the parasacral and sacrals are shifted anteriorly, projecting beyond the anterior margin of the centra. The neural arches of the caudosacrals are also shifted anteriorly but not projected as much as in the parasacral and sacrals. The fused, laterally compressed neural spines are directed posterodorsally. All dorsosacral, parasacral, sacral, and caudosacral ribs are fused to the vertebrae, although the dorsosacral ribs were separated during preparation. The first pair of dorsosacral ribs are arced similar to the posterior dorsal ribs, whereas the second and third pairs are rod-like and much shorter in length. Bone healing is observed on both sides of the distal portion of the first dorsosacral ribs (Fig. 5i,j). In dorsal view, the ribs become shorter in mediolateral length from the parasacral to the second caudosacral. The parasacral ribs are slightly projected posterolaterally and contact the ilium. The sacral ribs are projected laterally, whereas the caudosacral ribs are projected anterolaterally. All sacral and caudosacral ribs are fused to the ilium. In lateral view, the ribs from the parasacral to caudosacrals project at a lower angle. The parasacral and the first sacral ribs have anteroposteriorly expanded dorsal and ventral www.nature.com/scientificreports/ surfaces, which give the ribs an I-shaped proximal cross-section (Fig. 6a,b). The cross-section of the proximal ribs becomes trapezoid from the second sacral to the second caudosacral (Fig. 6a,b). Fourteen caudal vertebrae are fused into the handle (sensu 19 ) of the tail club ( Fig. 7h-j). The handle is slightly curved in a posterodorsal fashion. The first to second handle vertebrae are spool-shaped, similar in length and height, and laterally compressed. The centra are elongated from the third back, except the fourteenth handle vertebra is short and knob-like. The low neural spines are blade-like and decrease in height posteriorly. The paddle-like prezygapophyses project beyond the anterior margin of the centra and diverge at an angle of about 20º in dorsal view. These interlock with the distally elongate, wedge-like postzygapophyses of the adjacent vertebra. The small nub-like transverse processes are only present on the first to third handle vertebrae. The chevrons are elongate, blade-like, and fused onto the ventral surface of the centra. Well-preserved ossified tendons were uncovered along the lateral surfaces of the tail handle (Fig. 8a-d). These tendons are arranged in an imbricating pattern. Most of these tendons are flattened and periodically bifurcate into as many as four branches. Some of the tendons are elliptic in cross-section and have longitudinal striae on the external surfaces. A poorly healed pathology is preserved on one tendon. This tendon was fractured in the middle portion and fused at an angle of about 125º.
The ilia diverge at an angle of 27º from the body midline (Fig. 6). The long and blade-like preacetabular processes extend anterolaterally (Figs. 6g, 7a-d). The closed acetabulum is level with the medially positioned first and second sacral ribs. The short buttress-like postacetabular processes extend posterolaterally from the acetabulum. Only the proximal portion of the left ischium is preserved (Fig. 8e-g). The proximal margin of the ischium is convex.
No cervical half-rings were preserved. A single isolated Type 2 osteoderm (sensu 6 ) is preserved (Fig. 8e-g). The external surfaces are damaged, whereas the medial surface is well preserved and smooth. It is polygonal in dorsal view, sharply keeled dorsally, and thin-walled. The surface texture of this osteoderm is pitted. Two isolated Type 7 osteoderms are sub-circular, dorsoventrally flattened, with a keel near one edge that tapers in height anteriorly (Fig. 8h-m). They are thin-walled, rugose externally, and have a slightly concave ventral surface. The distal end of the tail is enveloped with two major and one minor osteoderms, forming a tail club knob (Fig. 7j). The length and width of the tail club knob are nearly equal. The major osteoderms are hemispherical in dorsal view and dorsolaterally keeled; the right is slightly larger than the left. Pathological grooves are present along the lateral surfaces of the two major osteoderms (Fig. 7k). The minor osteoderm is a rhomboidal shape in dorsal view. The surface texture of the tail club knob is pitted, rugose, and spongy. Phylogenetic analysis. The phylogenetic analysis resulted in a single most parsimonious tree (tree length = 32 steps, consistency index = 0.813, and retention index = 0.769) (Fig. 9a). T. tumanovae is a sister to the clade that includes T. kielanae and T. teresae. T. kielanae and T. teresae share only one synapomorphy: moderate-sized basioccipital foramen (character 10: state 1). The three Tarchia species shares five synapomorphies: a "neck" present at the base of the quadratojugal horn (7:1) (ambiguous in T. kielanae); tall foramen magnum (9:1); tall braincase (15:1); posteroventrally oriented occipital condyle (17:1); no postocular caputegulae (21:0) (ambiguous in T. kielanae). Saichania, Zaraapelta, and Minotaurasaurus were recovered as successive outgroups to the clade containing three Tarchia species.

Discussion
The squamosal horns of MPC-D 100/1353 are divided into the external layer of the squamosal horn and the underlying squamosal horn propers (Fig. 2a,b,e,f). These two layers are tightly fused in the holotype of Zaraapelta (MPC-D 100/1338) 17  Evidence of fracture healing can be observed on both sides of the first dorsosacral ribs of MPC-D 100/1353, in the anterolateral part of the pelvic area (Fig. 5i,j). Arbour et al. 21 suggested that the localized injuries on the pelvic area in ankylosaurines are likely caused by intraspecific combat inflicted by the tail club knob. This idea was supported by the fact that there was no relationship between tail club knob size and predator body mass, and concentrated pathologies observed in free caudal vertebrae, tail club knobs, and pelvic osteoderms confined to mature individuals 21,22 . A poorly healed ossified tendon on the tail knob handle is present in the holotype of T. tumanovae (Fig. 8d), which is a possible injury due to active tail use during combat.
The tail club knob of the holotype also has pathologies. Grooves are present along each lateral surface of the two major osteoderms (Fig. 7k). Moreover, the tail club knob is asymmetric in dorsal view, the left major osteoderm being shorter in mediolateral width than the right (Fig. 7j). Similar asymmetric bone growth has been observed on the postorbital horn of a male Dall sheep (Ovis dalli dalli), which impact each other with their heads at intraspecific combat 23  www.nature.com/scientificreports/ and UALVP 16247 (Ankylosauridae indet.) 22,24,25 . UALVP 16247 is similar to MPC-D 100/1353 by having the left major osteoderm smaller in volume than the right 24 . In ROM 788 and ROM 75860, however, the left major osteoderm is larger than the right 22,25 . Modern African elephants (Loxodonta) show tusk asymmetry due to side preferences in tusk use in stripping bark, digging root, and during agonist interactions [26][27][28] . Comparable to the tusk asymmetry in African elephants, tail asymmetry may relate to side preferences in tail use among ankylosaurine taxa or individuals. Pathologies found on the pelvic area and tail of MPC-D 100/1353 provide additional evidence of agonistic behaviour in ankylosaurines. Laterally wide trunks of ankylosaurines could have protected vital organs from being ruptured during conspecific tail club strikes 29 . Sub-rectangular broad muzzles are a morphological character of low-level bulk feeders, whereas anteriorly protruded shovel-shaped muzzles are selective feeders in ankylosaurines [30][31][32] . Similar dietary adaptations based on rostral shape are also known in mammalian herbivores, such as ungulates and ground sloths [33][34][35][36] . Based on these examples, both Tarchia species from the Nemegt Formation were probably selective feeders (Fig. 9b).
Bulk feeding ankylosaurines were present before the Baruungoyot and the Nemegt "age" (middle Campanian-lower Maastrichtian) based on known skull specimens (Fig. 9b). On the other hand, all ankylosaurines from the Nemegt and the Baruungoyot formations were adapted for selective feeding. These dietary shifts in ankylosaurines probably relate to habitat change, the shift from semi-arid (Bayanshiree Formation) and arid (Djadokhta and Baruungoyot formations) to more humid climates (Nemegt Formation) 31,37 . Climate-driven habitat change alters the plant communities in the environment 38 , and niche shifting in ankylosaurids might have responded to this. Recently, Jerzykiewicz et al. 39 proposed that Djadokhta, Baruungoyot, and Nemegt Formations are coeval, and the Nemegt Gobi Basin can be visualized as an ephemeral lake surrounded by semi-arid alluvial plains and arid dune fields. If this is the case, the dietary difference between these Mongolian taxa may result from habitat differentiation within the same basin.

Methods
Fossil preparation of the studied specimen (MPC-D 100/1353) was done at a laboratory in Hwaseong City of South Korea in 2012. The specimen was returned to Mongolia in 2016 and is now permanently held in the Institute of Paleontology in Ulaanbaatar, Mongolia. All measurements were taken using a measuring tape and a digital caliper. Comparisons to other ankylosaurid taxa were made by examining some specimens in the Institute of Paleontology, Mongolian Academy of Sciences, Mongolia, or were extracted from published literature. The term 'caputegulum' (sensu 15 ) was used to refer to the cranial ornamentation of ankylosaurs. Osteoderm types correspond to the terms used by Arbour et al. 6 .
For the phylogeny of the new T. tumanovae, the character list and data matrix used in this study was modified from that of Penkalski and Tumanova 8 (Supplementary Data S2). The modifications include the following: modified one character (5); added two new characters (22 and 23); revised a few character states of Pinacosaurus grangeri (character state 5:0 to 15:1) and PIN 3142/250 (holotype of Tarchia teresae) (5:? to 5:1, and 6:0 to 6:?). Including T. tumanovae, seven taxa with 23 characters (Supplementary Data S3) were analyzed in TNT (Tree Analysis Using New Technology) version 1.1 51 . A traditional search was performed, using implied character weighting and a k-value of 3.