A new cynodont from the Upper Triassic Los Colorados Formation (Argentina, South America) reveals a novel paleobiogeographic context for mammalian ancestors

Probainognathia is a derived lineage of cynodonts which encompass Mammalia as their crown-group. The rich record of probainognathians from the Carnian of Argentina contrasts with their Norian representation, with only one named species. Here we describe a new probainognathian, Tessellatia bonapartei gen. et sp. nov., from the Norian Los Colorados Formation of the Ischigualasto-Villa Unión Basin of Argentina. The new taxon, represented by a partial cranium with associated lower jaws, was analyzed through neutron and X-rays micro-tomography (μCT). The high-resolution neutron μCT data allowed the identification of a unique character combination, including features inaccessible through traditional techniques. We constructed the largest phylogenetic data matrix of non-mammalian cynodonts. The new species and its sister taxon, the Brazilian Therioherpeton cargnini, are recovered as probainognathians, closely related to Mammaliamorpha. We conducted the first quantitative paleobiogeographic analysis of non-mammalian cynodonts, focusing in probainognathians. The results indicate that Probainognathia and Mammaliamorpha originated in southwestern Gondwana (in the Brazilian Paraná Basin), which was an important center of diversification during the Triassic. Finally, the Chinese Lufeng Basin is identified as the ancestral area of Mammaliaformes. These new findings, besides adding to the knowledge of the poorly represented Norian cynodonts from the Los Colorados Formation, are significant to improve our understanding of probainognathian diversity, evolution, and paleobiogeographic history.

few fragmentary specimens. This is in concordance with the general scarcity of cynodonts from Norian strata globally 3,4 (see Supplementary Information).
The Norian Los Colorados Formation of the Ischigualasto-Villa Unión Basin of Argentina is renowned for its vertebrate fossil record 5,6 . One of the youngest dicynodonts of Gondwana, Jachaleria colorata, is represented in the lower levels of the unit [7][8][9] . The upper levels of the Los Colorados Formation represent a different faunal assemblage 10 that has provided one of the oldest turtles 11,12 ; basal and derived representatives of the crocodylian lineage 5,[13][14][15][16] ; and a number of relatively well-represented dinosaurs 5,[17][18][19][20][21][22] . The latter is key to the understanding of the rise of this group as it is the earliest assemblage where dinosaurs are dominant components of the ecosystem 23 . Non-mammaliaform cynodonts (NMC) are scarcely represented in this unit 14,24,25 . Only two partial skulls of Chaliminia musteloides 24,25 and a few fragmentary postcranial elements of an unnamed taxon 5,14,26 have been reported. They are among the oldest tritheledontids, a lineage of mainly Gondwanan cynodonts that have been proposed to be closely related to basal mammals 25,[27][28][29] . Since 2014, renewed exploration efforts in the Los Colorados Formation led to new findings at the Parque Nacional Talampaya (La Rioja Province, north-western Argentina), including five cynodont specimens whose preliminary study suggests that at least three of them might represent previously unrecognized taxa 30 . Recent geochronological studies of the unit suggest that the richest and most renowned faunal assemblage ('La Esquina local fauna'), found in the upper levels of Los Colorados Formation, might be mid-Norian (~ 220-211 Ma) in age 31,32 .
Herein we describe a new species of derived NMC found in the upper levels of the Los Colorados Formation ( Fig. 1) represented by a partial cranium with associated lower jaws. The specimen was analyzed through neutron and X-rays micro-tomography (μCT). Although only rarely used to analyze fossil specimens, it has become clear that neutron μCT complements X-ray μCT, including circumventing some problems that might arise from the latter methodology depending on the characteristics of the fossil sample studied 33,34 . In addition, we constructed the most comprehensive data-matrix for non-mammalian cynodonts published to date, allowing us to analyze the phylogenetic relationships of the new species and the main hypotheses regarding the ancestry of mammaliaforms and their close relatives. We also present the first quantitative paleobiogeographic analysis of non-mammalian cynodonts, focusing in probainognathians. Our results provide new insights regarding the evolutionary and paleobiogeographic history of Probainognathia. In particular, the new finding improves our knowledge on the diversity of the poorly represented Norian forms in the rich fauna of the Los Colorados Formation, northwestern Argentina.

Type species. Tessellatia bonapartei.
Etymology. From the Latin tessella (each one of the tiles composing a mosaic), in reference to the combination of basal and derived features recognized in this taxon.
Diagnosis. Same as for species.
Etymology. In honor to the late Dr. José F. Bonaparte, who worked unrelentingly to broaden our knowledge of Mesozoic ecosystems and described the first cynodont remains from the Los Colorados Formation.
Holotype. PULR-V121, partial cranium, represented by the snout and orbital region, with associated lower jaws.
Locality and horizon. Upper third of the Los Colorados Formation, Talampaya National Park, La Rioja Province, Argentina. The specimen was found in a massive-to parallel-laminated sandy mudstone interval locally interbedded with parallel-to rippled-laminated sandstone, representing deposition in a floodplain setting sporadically affected by sandy splays from the fluvial channels. See Supplementary Information online for further details on the geological framework.
Diagnosis. Small probainognathian with the antorbital region notably longer than the height of the skull at the level of the anterior margin of the orbit. Autapomorphies of Tessellatia bonapartei are: masseteric fossa and low coronoid process posterior to the lower tooth row; angular process ventrally projected, well-developed, and semicircular in outline; broad groove between the lateral wall of the dentary and the postcanine alveoli; and upper postcanine count larger (14) than the lower postcanines (7), with the last six upper postcanines lacking an opposing lower tooth. Differs from other non-mammaliaform cynodonts except Riograndia guaibensis and the early mammaliaform Morganucodon in the presence of a large upper canine and a small lower canine. Differs from other non-mammaliaform probainognathians (except Aleodon cromptoni, Prozostrodon brasiliensis, Pseudotherium argentinus, and Tritylodon longaevus) in the subequal contribution to the secondary bony palate of the maxilla and palatine. Shares with many derived non-tritylodontid probainognathians and the mammaliaform Morganucodon an anteriorly elevated alveolar margin of the dentary. Shares with some derived probainognathians traditionally grouped in Tritheledontidae (Chaliminia musteloides, Elliotherium kersteni, Irajatherium hernandezi, and Pachygenelus monus), the presence of a ventrally bowed bony secondary palate with deep, narrow, lateral groove for the lower postcanines; posterior upper postcanines with convex labial and concave www.nature.com/scientificreports/ lingual surfaces, and bearing a centrally placed, symmetrical main cusp with convex mesial and distal margins flanked by smaller, lingually placed accessory cusps; and lower postcanines, comparatively larger than uppers, with a large, asymmetrical, mesial main cusp followed by smaller distal accessory cusps. Differs from Chaliminia musteloides, the only other named cynodont known from the Los Colorados Formation, in having the alveolar margin of the maxilla sigmoidal, with its dorsal-most point at the level of the canine; the horizontal ramus of the dentary comparatively low and lacking a well-developed platform lateral to the posterior postcanines; and the posteriormost lower incisor smaller than the canines and postcanines. www.nature.com/scientificreports/ brasiliensis, Pseudotherium argentinus, and Brasilodon tetragonus (but not Therioherpeton cargnini, Elliotherium kersteni, and Chaliminia musteloides) among derived non-tritylodontid probainognathians (NTP). The alveolar margin of the maxilla has a sigmoidal outline in lateral view ( Fig. 2b-c). The maxillary canal in Tessellatia bonapartei is completely ossified (Fig. 3) and separated from the maxillary sinus, as in derived probainognathians including mammaliaforms. The short infraorbital canal and posteriorly placed infraorbital foramen of Tessellatia bonapartei ( Fig. 2c; infraorbital foramen 2 in 39 ; see 40 ), is shared with derived probainognathians, suggesting the presence of a mobile rhinarium with sensitive vibrissae (see 40,41 ).
A remarkably large facial exposure of the lacrimal (Fig. 2b, c) is a distinctive feature of Tessellatia bonapartei that, among probainognathians, is only recognized in some tritylodontids. Unlike in other probainognathians, Tessellatia bonapartei shares with its sister-group Therioherpeton kersteni, some tritylodontids, and Morganucodon, a short lacrimal-frontal contact and frontals with short anterolateral projections.
Tessellatia bonapartei is more comparable to early-diverging probainognathians in the presence of an orbital process of the palatine that is excluded from the orbital wall and not contacting the frontal. Unlike in derived NTP, in which the secondary palate extends to the end of the tooth row or posterior to it, the osseous secondary palate of Tessellatia bonapartei ends anteriorly to the antepenultimate upper postcanine (PC13; Fig. 2e, f). On the other hand, Tessellatia bonapartei shares with some derived probainognathians (e.g., Chaliminia musteloides, Elliotherium kersteni, Irajatherium hernandezi, and Pachygenelus monus), traditionally grouped in Tritheledontidae (see 25 ), the presence of a ventrally bowed osseous secondary palate reaching the level of the crown of the upper postcanines (Fig. 2b, c). In these forms, deep, narrow grooves are present medial to the upper tooth rows, posterior to the maxilla-palatine suture (Fig. 2e, f).
Lower jaw. The dentary is relatively robust, high, and lateromedially wide (Fig. 2b, c, h). Unique to Tessellatia bonapartei is the presence of a broad groove between the lateral wall of the dentary and the postcanine alveoli (Fig. 2g). This groove starts at the posterior half of pc2 and ends laterally to the penultimate lower postcanine (pc6). The anterior region of the dentary is upwardly bent regarding the postcanine line in Tessellatia bonapartei (Fig. 2b, c, h) as well as in many derived NTP and the mammaliaform Morganucodon. The shape and development of the angular process of the dentary of Tessellatia bonapartei distinguishes it from other Triassic cynodonts.
Dentition. The upper dentition includes at least one incisor, one canine, and 14 postcanines (1 + ?I/1C/14PC) ( Fig. 2b-f). The lower dentition includes one incisor preserved in situ in the dentary and an isolated fragment of another, one canine, and seven postcanines (2 + ?i/1c/7pc) (Fig. 2b, c, g, h). Incisors are small, non-procumbent teeth (Fig. 2b-h), unlike those of some tritheledontids (see 25 ) such as Chaliminia musteloides, Pachygenelus www.nature.com/scientificreports/ monus, and Riograndia guaibensis. The presence of a large upper canine (1.5 or more the mesiodistal length of the first postcanine) in conjunction with a reduced lower one (the same mesiodistal length or less of the first postcanine) (Fig. 2b- (Fig. 2i, j) is similar to those of tritheledontids (see 25 ). A reduced number of lower postcanines (7) compared to the number of upper ones (14) is a distinctive feature of Tessellatia bonapartei among NMC (Fig. 2b-h). It is noteworthy that the lower tooth row is very short when compared to the upper one. The last lower postcanine (pc7) would have occluded approximately between PC8 and PC9. In this scenario, the last five upper postcanines lacked a lower counterpart.
Comparisons with Chaliminia musteloides. Chaliminia musteloides, represented by cranial remains, is the only previously named cynodont from the Los Colorados Formation. The known specimens are comparable in size to Tessellatia bonapartei. Hence, we deem it important to highlight some additional differences of the published specimens of Chaliminia musteloides (the holotype PVL 3857 and the referred specimen PULR 088) with the holotype of Tessellatia bonapartei. Tessellatia bonapartei has a proportionally longer snout when compared to the height of the skull at the level of the anterior margin of the orbit (Supplementary Table S3). When compared to the holotype of Chaliminia musteloides, Tessellatia bonapartei has a comparatively slenderer dentary. On the other hand, the proportions of the horizontal ramus of the dentary are similar between Tessellatia bonapartei and the referred specimen of Chaliminia musteloides (PULR 081) (Supplementary Table S4). However, it is important to take into account that the anterior end of the skull and lower jaws are missing in Tessellatia bonapartei; a longer skull and dentary would mean a longer snout and a more slender horizontal ramus of the dentary, accentuating the differences in these proportions with the specimens of Chaliminia.
In addition, Tessellatia bonapartei ( Fig. 2) differs from Chaliminia musteloides in the presence of a constriction in the snout posterior to the canines that results in an expanded rostrum anteriorly; a strongly concave posterior margin of the secondary palate; a notably short lower tooth row with the ascending coronoid process of the dentary well-posterior to the last lower postcanine, slightly sloping dorsally at the level of the orbit; the masseteric fossa does not reach the level of the last lower postcanine; a broad groove between the lateral wall of the dentary and the postcanine alveoli; and small non-procumbent lower incisors (Fig. 2). Tessellatia bonapartei lacks the strong osseous platform lateral to the last lower postcanines observed in the holotype of Chaliminia musteloides. The lower canine of Tessellatia bonapartei has approximately the same diameter than the first lower postcanine, as measured at the upper portion of the root, whereas the lower canine of Chaliminia musteloides is markedly larger than the postcanines. Unlike the condition of Chaliminia musteloides, there is no lower postcanine diastema in Tessellatia bonapartei.
Despite sharing the general morphology, there are some differences in the postcanine dentition between Tessellatia bonapartei and Chaliminia musteloides. Unlike the convex labial and concave lingual surface of the upper postcanines of Tessellatia bonapartei, in Chaliminia musteloides there is a central blunt crest separating mesial and distal depressions labially whereas the lingual surface is almost flat. Contrary to what is observed in Tessellatia bonapartei, the lower postcanines of Chaliminia musteloides are mesiodistally shorter than the upper ones. In Chaliminia musteloides, the lower postcanines show more conspicuous, higher distal accessory cusps, more closely placed to the main cusp than in Tessellatia bonapartei.
Phylogenetic results. The parsimony analysis of 151 characters and 73 operational taxonomic units (OTUs) resulted in 29,304 most parsimonious trees (MPT) of 786 steps. The strict consensus is poorly resolved. There is a large polytomy including a wide array of taxa. A basal group of Cynodontia (represented by late Permian forms), Epicynodontia, and Eucynodontia, as usually conceived 2,28,42 , are not recovered. Only Cynognathia and a clade including some derived probainognathians are recovered by our analysis. Tritylodontidae, the only clade with strong support, is nested among the latter group.
A second analysis excluding three wild-card taxa (i.e., Diegocanis elegans, Charruodon tetracuspidatus, and Microconodon tenuirostris) retrieved 200 MPT of 782 steps, whose strict consensus is better resolved than with the complete dataset (Fig. 4). A basal group of late Permian cynodonts is followed by Cynosaurus suppostus, a node that perhaps represents the Epicynodontia level. The next node is a large polytomy including several terminals. Traditional taxa of Probainognathia, such as Chiniquodon spp., Trucidocynodon riograndensis, Ecteninion lunensis, Candelariodon barberenai, Probainognathus jenseni + Bonacynodon schultzi, and Aleodon cromptoni + Aleodon brachyrhamphus, are recovered as part of the mentioned polytomy that also includes basal members of Epicynodontia. In any case, resolving this huge polytomy that includes Eucynodontia or obtaining a clade including all taxa traditionally considered Probainognathia is not our goal. Our aim is to evaluate the position of Tessellatia bonapartei among other probainognathian, which our matrix and methodology allows us to do. Cynognathia and Tritylodondidae form monophyletic groups as in the complete analysis. This is not the case for Tritheledontidae, whose members are found to belong to the Mammaliamorpha clade (as defined by 43 ), as stem-taxa to Tritylodontidae or to Mammaliaformes. According to this, Mammaliamorpha includes tritylodontids, mammaliaforms, and a number of derived probainognathians. Tessellatia bonapartei is recovered as sister-group of the Brazilian Therioherpeton cargnini based on two unambiguous synapomorphies: axis of the posterior region of the maxillary tooth row directed towards the center of the subtemporal fossa (character 98:1) and transverse axis of crown strongly oblique to midline axis of the skull (character 114:1). The clade formed by Tessellatia bonapartei + Therioherpeton cargnini is in a basal position among the reduced-probainognathian A traditional Probainognathia is portrayed in the majority rule consensus tree of the reduced, second analysis ( Supplementary Fig. S5). It shows a basal polytomy integrated by Candelariodon barberenai, a chiniquodontid clade (Chiniquodon + Aleodon), and a clade with remaining probainognathians, which includes a polytomy formed by three clades: (Bonacynodon schultzi + Probainognathus jenseni), (Trucidocynodon riograndensis + Ecteninion lunensis), and the other probainognathians. In this topology, the placement of Lumkuia fuzzi, the oldest and only basal Probainognathia from South Africa, is counterintuitive, and different to previous phylogenies in which it was recovered as the basal-most Probainognathia (e.g., 2,43 ) or as the basal-most Eucynodontia, outside of Cynognathia and Probainognathia (e.g. 42 ). In our majority rule tree, basal-most Probainognathia are represented by a group of largely South American Late Triassic taxa.

Probainognathian paleobiogeography. A Bayesian Binary Markov Chain Monte Carlo (BBM) analy-
sis of 100,000 cycles and 100 chains was performed, including nine a priori determined areas of provenance (Supplementary Table S1) and allowing the maximum number of areas per node. The study was based on the probainognathian clade as shown in the majority rule consensus tree ( Supplementary Fig. S5) obtained from the analysis of the pruned dataset (70 OTUs, 151 characters). Additionally, we performed a Hausdorf 's 44 manual optimization of Probainognathia, Mammaliamorpha, and Tritylodontidae nodes (Supplementary Tables S5, S6, S7, S8).
The results (Fig. 5) indicate that the paleobiogeographic history of probainognathians was dominated by dispersions (22 events) and sympatric speciation events (18 events), followed by vicariances (13 events) and with a single extinction event. Probainognathia originated in southwestern Gondwana (Paraná Basin, Brazil; 83.53%) and this was followed by a large and relatively quick increase in the diversity. This diversification event included    Tables S5−S6). According to our results, mammaliamorphs thrived in the Paraná Basin and experimented a few dispersions to the Ischigualasto-Villa Unión Basin (Argentina; 2 events), the Lufeng Basin (China; 1 event), and to North American basins (1 event). It is in the Lufeng Basin where the origin of Mammaliaformes is expected (61.54%). However, it must be noted that we did not include Gondwanadon tapani and Tikitherium copei, from the upper Carnian-lower Norian Tiki Formation (Madhya Pradesh, India), in our analysis as they are both represented by an isolated tooth which is even incompletely preserved in the case of Gondwanadon [45][46][47] . Although never tested phylogenetically, Gondwanadon has been considered to be closely related to Morganucodon 48  www.nature.com/scientificreports/ originally interpreted to be related to Morganucodon, is now regarded as a docodontan, more closely related to the Rhaetian Delsatia rhupotopi and Woutersia spp. from France 50 . The Tiki Formation records the only faunal assemblage from the Triassic of Gondwana featuring remains of mammaliaforms (in fact, the oldest record of the group) and the early recognition of two mammaliaforms with clearly different tooth patterns suggests an undocumented worldwide evolution of this group. By the end of the Norian, dispersion events from the Paraná Basin to the Karoo Basin (South Africa) are identified. In this latter region, they diversified and dispersed to the Lufeng Basin (1 event) and to North American basins (1 event). The origin of Tritylodontidae is inferred to have occurred in the Karroo Basin (33.53%), whilst Hausdorf 's optimization indicates that the origin of this group occurred in the Lufeng Basin (Supplementary Table S7).

Discussion
There is a rich Triassic record of cynodonts in the Ischigualasto-Villa Unión Basin of Argentina, particularly remarkable in the Carnian Chañares and Ischigualasto formations 51,52 . The representation of this group becomes extraordinarily scarce in the fossil-rich Los Colorados Formation ('La Esquina local fauna'), the youngest Triassic unit of the basin, where only one species represented by two small skulls and a second putative small taxon known by a few fragmentary postcranial remains are known to date 14,[24][25][26] . The finding of Tessellatia bonapartei, thus reinforces the knowledge of a subdued fauna of small animals in the Los Colorados Formation, only represented by cynodonts so far. Tessellatia bonapartei displays a mosaic of basal and derived features. Despite this, it is unambiguously recovered as a derived probainognathian. Tessellatia bonapartei bears a postcanine morphology very similar to that of taxa traditionally grouped in Tritheledontidae. Although this might point to a similar diet, the distinctive features in the dentition of Tessellatia bonapartei (e.g., a notably reduced number of lower postcanines regarding upper ones, large upper canine in conjunction with a small lower one, and small-sized incisors) suggest a unique feeding-habit. The cynodonts from the Los Colorados Formation fill a temporal gap in the South American record of Probainognathia, between the early Norian forms from the Brazilian Caturrita Formation and the late Norian-Rhaetian, still undescribed, cynodonts from Quebrada del Barro Formation, San Juan, Argentina 8 . The faunal association from the Los Colorados Formation also contrasts with these latter South American assemblages because it is numerically dominated by large taxa (e.g., sauropodomorphs), with only a few known specimens of the rare small-sized cynodonts present (see Supplementary Information). On the other hand, small cynodonts are diverse and/or abundant in the Caturrita and Quebrada del Barro formations 8,53 .
Contrasting with the limited record of Probainognathia (and cynodonts) of the Los Colorados Formation, the upper Carnian-Norian record of the Santa Maria and Caturrita formations (Brazil) shows a burst of diversity of small derived probainognathians with at least eight species 4 . In fact, the Brazilian record of the group provides most of the information regarding Probainognathia during that time in the world. A second, very important pulse of diversification of this lineage is clearly headed by mammaliamorphs mostly from the Lower Jurassic. In this case, tritylodontids make the great difference because extreme modifications of the skull, dentition and, together with early mammaliaforms (i.e., Morganucodon and allies), of the postcranium. Tritylodontids were also abundant and diverse, with a cosmopolitan distribution, and the only cynodont lineage of this time reaching comparatively large body sizes.
Basal and derived probainognathians are recognized during the Carnian, including some relatively earlydiverging forms that attained the largest body-masses of the Triassic representatives of the clade (e.g., Chiniquodon, Trucidocynodon riograndensis; 51,54 ). During the Norian, only small-sized probainognathians are represented. This scenario suggests that the Carnian-Norian transition represented a pivotal moment for the diversity and ecology of cynodonts in general and probainognathians in particular, probably related to the Carnian Pluvial Episode 55 .
According to our phylogenetic results (Fig. 4), the two taxa from Los Colorados are not closely related. Tessellatia bonapartei is an early-diverging taxon whereas the comparatively more crownward positioned Chaliminia musteloides is more closely related to tritylodontids than to mammaliaforms. The calibrated phylogeny together with the quantitative paleobiogeographic analysis of probainognathians (Fig. 5) presented here point to two major diversification events. One of them took place during the Middle Triassic and the earliest Late Triassic and included the origin of Mammaliamorpha. The second occurred mainly in the Norian and the Rhaetian and is represented by the origin of tritylodontids and closely allied forms. It is by mid-Norian times that the origin of Mammaliaformes would have occurred; however, mammaliaforms are not recognized until the latest Triassic and Jurassic. The results of our analysis imply long ghost-lineages for many taxa, including Tessellatia bonapartei, suggesting a still unrecorded history of Probainognathia.
Dispersal events from the Paraná Basin during the Carnian-earliest Norian were restricted to the Ischigualasto-Villa Unión Basin whereas during the late-Norian-Rhaetian dispersions are from the Paraná Basin to the Karoo Basin. This suggests that, during these two time-lapses, barriers or favorable conditions to dispersal were alternatively in place in the Ischigualasto-Villa Unión and Karoo basins preventing or encouraging dispersions from the Paraná Basin. It is also interesting to note that the three dispersion events from the Paraná Basin to the Ischigualasto-Villa Unión Basin represented by Tessellatia bonapartei, Chaliminia musteloides, and Pseudotherium argentinus were not followed by diversifications. These results highlight the importance of improving the knowledge on Norian probainognathians.
Our study reveals a new species of a probainognathian cynodont, adding to the poorly known cynodont diversity from the rich Norian faunal assemblage of the Los Colorados Formation. The phylogenetic evaluation of Tessellatia bonapartei through the most comprehensive non-mammalian cynodonts data-matrix assembled to date together with the first quantitative paleobiogeographic analysis of probainognathians highlight a previously unrecognized evolutionary and biogeographic history.

Methods
Micro-tomography. The small size and delicate bones make it impossible to remove the hard rock matrix without damaging the specimens and losing important information. In order to overcome this issue, PULR-V121 was analyzed through X-ray micro-tomography in YPF Tecnología S.A. (Y-TEC, Ensenada, Buenos Aires, Argentina) using the Bruker SkyScan 1173 instrument. The equipment was set up to 100 kV and 80µA. A total of 900 images of the specimen were captured through a 360° tomography (rotation step 0.4°) with an exposure time of 250 ms with two frames averaged. The experimental design resulted in a 40.01 µm pixel size. The tomographic reconstruction was then produced with the software NRecon 1.6.9.8 (https:// www. micro photo nics. com/ microct-syste ms/ nrecon-recon struc tion-softw are). Although obtaining acceptable results, the resolution was not ideal. Regrettably, the reconstructed images proved to be difficult to interpret, as bone and matrix were in some regions indistinguishable from each other. We acknowledged that this issue was probably technique-related due to the presence of ferruginous material in the sediment. Hence, we decided to perform a preliminary neutrontomography at the RA-6 facility (Comisión Nacional de Energía Atómica, San Carlos de Bariloche, Argentina). Considering the promissory results obtained and the hypothesis that this specimen may represent a previously unregistered taxon of the very poorly represented mammalian ancestors in the late Late Triassic of Argentina, we performed a neutron tomography with the highest possible spatial resolution at the ANTARES instrument 56,57 in the Forschungs-Neutronenquelle Heinz Maier-Leibnitz Zentrum (FRM II, Garching, Germany). The specimen did not show previous radioactivity before introducing it directly to the reactor hall at ANTARES instrument. For the neutron tomography, PULR-V121 was wrapped in aluminum foil together with two additional specimens from the same stratigraphic levels (PULR-V222 and PULR-V223) to maximize the available beam-time. The package was placed in a 5 cm long slot of an aluminum cylinder. A small aluminum plate was fixed to the cylinder using aluminum tape to act as a floor. This setup feature stabilized the specimens during the tomography and allowed to place them as close as possible to the detector ( Supplementary Fig. S1). At ANTARES, a collimation ratio of L/D = 500 was used. The Andor Neo sCMOS detector was equipped with a 100 mm Zeiss Milvus f2.0 lens which allowed us to obtain high resolution images, with a 19.74 µm pixel size. We performed a standard (white-beam), 360° tomography employing a Gd 2 O 2 S based neutron scintillator of 6 cm × 6 cm of 20 µm thickness. The exposure time was 17 s and each angular position (every ~ 0.192°) was acquired three times for improving statistics. The neutron tomography took circa 17 h and 42 min. In order to normalize the images obtained, 19 open beam (open shutter, no sample in the beam) and 5 dark field (closed shutter) images were taken before and after the tomographic acquisition of the fossil remains, respectively.
The images were normalized and filtered using Image-J v. 1.52p 58 software and then reconstructed with Octopus Reconstruction 8.9.3.4 software (https:// octop usima ging. eu) at the Heinz Maier-Leibnitz Zentrum facility. Posteriorly, the reconstructed images were subjected to a new filtering process with the Inverse Scale Space Filter (ISS) module implemented in KipTool [59][60][61] . The ISS, an edge preserving de-noising filter based on the equation formulated by Burger et al. 62 , notably increased the sharpness of the images without sacrificing morphological information. The segmentation of the specimen was performed by Dr. Gaetano who also generated the 3D surface model of PULR-V121 using Avizo 7.1 (FEI; https:// www. therm ofish er. com/ ar/ es/ home/ elect ron-micro scopy/ produ cts/ softw are-em-3d-vis/ avizo-softw are. html). Approximately two weeks after the neutron tomography, PULR-V121 showed a very low decay ratio and was possible to remove it from the reactor hall.
Phylogenetic analysis. We put together a data matrix combining those of Liu and Olsen 42 and Ruta et al. 28 considering the subsequent modifications to both of them as well as modifying or deleting some of the characters and character states. We also added new characters and included the new specimen PULR-V121 as well as other relevant taxa. As a result, a comprehensive data matrix including 73 taxa and 151 characters was produced. Previous scorings were revised, and corrections implemented (Supplementary Appendices S1-S4). The aim of our phylogenetic analysis is to test the phylogenetic relationships of the new Norian taxon, Tessellatia bonapartei.
A first analysis was produced after the complete data matrix. After the recognition of three wild-card taxa (i.e., Diegocanis elegans, Charruodon tetracuspidatus, and Microconodon tenuirostris), we pruned them from the matrix and performed a second analysis considering only the 70 remaining taxa. TNT 1.5 software 63,64 was used for searching of most parsimonious trees. Routine used was the command xmult = level 10, that produce 14 autoconstrained replications; each replication with random sectorial searches, drifting (36 iters) and fusing (10 round), finding best score 1 time; followed by bb (bbreak) command that perform branch-swapping (tree bisection reconnection) using pre-existing trees. Characters were unordered and equally weighted. The analysis of the complete matrix resulted in 29,304 most parsimonious trees (mpt) of 786 steps; whereas the second analysis resulted in 200 mpt of 782 steps.
Most scholars agree on the putative probainognathian nature of dromatheriids, a group of small cynodonts represented by isolated teeth or fragments of mandible, mostly discovered in the Laurasia subcontinent and India 47 . However, the placement of dromatheriids has yet to be tested cladistically. Unfortunately, removal of Microconodon tenuirostris (one of the wildcard taxa in the complete analysis) from our second analysis, does not allow us to evaluate the relationships of this lineage.

Paleobiogeographic analysis.
A biogeographic analysis of the probainognathian clade as represented in the majority rule consensus tree of the pruned (70-taxa) phylogenetic analysis was performed in RASP 4.2 65 . The cladogram branch-length was temporally calibrated with RStudio 66 using PaleoTree package 67 with EQUAL methodology. The First Appearance Datum (FAD) and Last Appearance Datum (LAD) as well as the geographical distribution of the taxa were obtained from the relevant published sources (see 28 and 4 for a review). We considered nine ancestral areas according to the regions that fossils were found (Supplementary Table S1). Most areas are so distant from each other that it is not necessary to determine discontinuity between them. Addi- www.nature.com/scientificreports/ tionally, almost all of them have endemic taxa. We performed a Bayesian Binary Markov Chain Monte Carlo (BBM; 68 ) analysis with the calibrated cladogram, considering 100,000 cycles, 100 chains, and the maximum number of areas (nine) per node. We conducted an optimization analysis to test the possible ancestral area of relevant nodes using the Weighted Ancestral Area Analysis (WAAA; 44 ). The number of weighted gain steps (GSW), weighted loss steps (LSW) and the probability index (PI = GSW/LSW) were calculated manually. The PI of each area indicates the probability of that area as a part of the ancestral area. It is important to bear in mind that the incomplete fossil record will bias any quantitative analysis. Hence, the results presented here depend on the known record of the taxa considered and changes are expected as our knowledge improves.

Data availability
The holotype and only known specimen of Tessellatia bonapartei is deposited at the Museo de Ciencias Naturales, Universidad Nacional de La Rioja under the collection number PULR-V121. The 3D model of PULR-V121 is freely accessible in MorphoMuseuM (http:// morph omuse um. com 70 ).