# Oxygen isotopes suggest elevated thermometabolism within multiple Permo-Triassic therapsid clades

Kévin Rey, Romain Amiot, François Fourel, Fernando Abdala, Frédéric Fluteau, Nour-Eddine Jalil, Jun Liu, Bruce Rubidge, Roger Mh Smith, Jean-Sébastien Steyer, et al.

► To cite this version:

Kévin Rey, Romain Amiot, François Fourel, Fernando Abdala, Frédéric Fluteau, et al. Oxygen isotopes suggest elevated thermometabolism within multiple Permo-Triassic therapsid clades. eLife, 2017, 6, pp.e28589. 10.7554/eLife.28589. hal-01692941

HAL Id: hal-01692941

https://sde.hal.science/hal-01692941v1

Submitted on 23 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Oxygen isotopes suggest elevated thermometabolism within multiple Permo-Triassic therapsid clades

Kévin Rey¹,²*, Romain Amiot¹, François Fourel¹,³, Fernando Abdala², Frédéric Fluteau⁴, Nour-Eddine Jalil⁵, Jun Liu⁶, Bruce S Rubidge², Roger MH Smith²,⁷, J Sébastien Steyer⁵, Pia A Viglietti², Xu Wang⁸, Christophe Lécuyer¹,⁹

¹Univ Lyon, Université Lyon 1, Ens de Lyon, CNRS, UMR 5276 LGL-TPE, Villeurbanne, France, France; ²Evolutionary Studies Institute and School of Geosciences, University of the Witwatersrand, Johannesburg, South Africa; ³Univ Lyon, Université Lyon 1, CNRS, UMR 5023 LEHNA, Villeurbanne, France, France; ⁴Institut de Physique du Globe de Paris, Paris, France; ⁵Centre de Recherches en Paléobiodiversité et Paléoenvironnements, UMR 7207 CNRS-MNHN-UPMC, Museum National d'Histoire Naturelle, Paris, France; ⁶Key Laboratory of Vertebrate Evolution and Human Origins of Chinese Academy of Sciences, Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing, China; ⁷Iziko South African Museum, Cape Town, South Africa; ⁸Key Laboratory of Cenozoic Geology and Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China; ⁹Institut Universitaire de France, Paris, France

**Abstract** The only true living endothermic vertebrates are birds and mammals, which produce and regulate their internal temperature quite independently from their surroundings. For mammal ancestors, anatomical clues suggest that endothermy originated during the Permian or Triassic. Here we investigate the origin of mammalian thermoregulation by analysing apatite stable oxygen isotope compositions (δ¹⁸Oₚ) of some of their Permo-Triassic therapsid relatives. Comparing of the δ¹⁸Oₚ values of therapsid bone and tooth apatites to those of co-existing non-therapsid tetrapods, demonstrates different body temperatures and thermoregulatory strategies. It is proposed that cynodonts and dicynodonts independently acquired constant elevated thermometabolism, respectively within the Eucynodontia and Lystrosauridae + Kannemeyeriiformes clades. We conclude that mammalian endothermy originated in the Epicynodontia during the middle-late Permian. Major global climatic and environmental fluctuations were the most likely selective pressures on the success of such elevated thermometabolism.

DOI: 10.7554/eLife.28589.001

## Introduction

One key adaptation enabling tetrapods to cope with fluctuating climatic conditions was the acquisition of endothermy (Paaijmans et al., 2013). This character is defined here as the ability to actively produce body heat through metabolic activity (Cannon and Nedergaard, 2004). Its development and anchoring in populations constitutes a major step in vertebrate evolution because it modified the energy relationships between organisms and their environments. By actively raising and maintaining body temperature within a narrow range that allows optimal physiological and biochemical

*For correspondence: kevin.rey@wits.ac.za

**Competing interests:** The authors declare that no competing interests exist.

**Funding:** See page 14

**Received:** 12 May 2017

**Accepted:** 20 June 2017

**Published:** 18 July 2017

**Reviewing editor:** Diethard Tautz, Max-Planck Institute for Evolutionary Biology, Germany

© Copyright Rey et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
**eLife digest** School textbooks often refer to "cold-blooded" and "warm-blooded" animals, but these terms are misleading. Rather than being cold, animals like reptiles have body temperatures that are mostly determined by their external environment and can actually achieve high body temperatures, for example, by basking in the sun. By contrast, "warm-blooded" mammals produce their own heat and typically maintain a body temperature that is warmer than their environment. As such, so-called warm-blooded animals are more accurately referred to as "endotherms" and cold-blooded animals as "ectotherms".

Endothermic animals share several characteristics, including insulating layers – like fur or feathers – that keep the body warm, and a secondary palate that separates the mouth and nose for continuous breathing, even while eating. Many of these traits are seen in fossils belonging to a group of animals called the therapsids. Also known as the "mammal-like reptiles", these animals are descended from ectothermic reptiles but are the ancestors of the endothermic mammals. They dominated the land between 270 and 220 million years ago, during periods of time called the Permian and the Triassic. They also survived two major mass extinction events, including the most devastating mass extinction in all of Earth's history. However, when the ancestors of mammals became truly endothermic remains an open question. Previous studies that have tried to determine this by focusing on the physical characteristics of therapsids have not yet given a consistent date.

Rey et al. took a new approach to answer when endothermy first evolved in the mammal-like reptiles, and instead looked at the chemical makeup of minerals in over 100 fossils. Oxygen can exist in different forms called stable isotopes: oxygen-16 and the rarer and heavier oxygen-18. The ratio of these two isotopes in a fossil will depend on, among other things, where the animal lived and, importantly, its body temperature. Therefore, Rey et al. compared oxygen-containing minerals in the bones and teeth of therapsids to those of other animals that lived alongside them to look for signatures that indicated differences in body temperature and how it was regulated.

It appears that two different branches of the therapsid's family tree independently became endothermic. One branch includes the mammals and their direct ancestors, while the second is more distantly related to mammals. Both became endothermic towards the end of the Permian Period, between about 259 and 252 million years ago. Based on these findings, Rey et al. suggest that endothermy allowed these animals to better cope with fluctuating climates, which helped them to be among the few species that survived the mass extinction event at the end of the Permian.

Going forward, these new findings can help scientists to understand which physical characteristics were necessary for endothermy to first develop and which helped to optimize it afterwards. Furthermore, they also suggest that endothermic animals are more able to survive fluctuations in climate, which could guide efforts to protect modern-day endangered species that are most at risk from the ongoing effects of climate change.

DOI: 10.7554/eLife.28589.002

functioning, endothermic vertebrates are able to colonise environments with extreme thermal conditions, for example freezing at high latitudes and altitudes (Day et al., 2015). Endothermy is commonly associated with homeothermy, being the capacity to regulate the body heat through metabolic activity as well. This combination corresponds to one end of a gradient of thermoregulatory strategies observed in living animals. The other end of the spectrum is ectothermy combined with poikilothermy which animals use as a thermoregulatory strategy to increase their body temperature toward optimal levels by using external heat sources. Their body temperature therefore traces that of their surroundings and is the most commonly occurring energy saving strategy. Amongst modern vertebrates, various thermoregulatory strategies have been adopted between these two end-members, such as regional endothermy (Bernal et al., 2001; Katz, 2002) or inertial homeothermy (McNab and Auffenberg, 1976), and only birds and mammals fall within the endothermy end of the spectrum. It has been proposed that bird thermoregulation originated within non-avian dinosaurs (Seebacher, 2003; Amiot et al., 2006; Grady et al., 2014), or even earlier within basal archosauriforms (Farmer and Carrier, 2000; de Ricqlès et al., 2003; Seymour et al., 2004; Summers, 2005; Gower et al., 2014; Legendre et al., 2016). Various approaches have been tried by
many researchers to assess the origin of mammalian endothermy (McNab, 1978; Bennett and Ruben, 1986; Farmer, 2000; Hillenius, 1992, 1994; Kemp, 2006a; Khaliq et al., 2014; Owerkowicz et al., 2015; Benoit et al., 2016b; Crompton et al., 2017). Some consider the appearance of endothermy to have either occurred during the transition from basal synapsid 'pelycosaurs' to therapsids, and to be either due to a shift in foraging ecology (Hopson, 2012) or due to a response to the availability of a seasonally arid, savanna-like biome by the end of the Early Permian (Kemp, 2006b).

How and why endothermy evolved in mammals remains a contentious issue, mostly because of the very low fossilization potential of anatomical and behavioural features associated with thermo-regulation. Amongst the latter features, the presence of hair as an insulating integument is unequivocally associated with endothermy in all extant mammals. The oldest synapsid fossils preserved with fur imprints are Castorocauda (Ji et al., 2006) and Megaconus (Zhou et al., 2013). These early relatives of mammals were recovered from the Middle-Late Jurassic of China, implying that hair and fur appeared before ~165 Ma. The occurrence of retracted, fully ossified and non-ramified infraorbital canals (a structure associated with the presence of maxillary vibrissae) within non-mammaliaform Prozostrodontia, implies an older age of approximately 240 to 246 Ma for the occurrence of fur and hair (Benoit et al., 2016b).

Another anatomical character interpreted as associated with endothermy is the bony secondary palate. This is a feature associated with efficient respiratory capabilities considered to be linked to the high energy required for elevated metabolic rates. In some Triassic non-mammaliamorph therapsids, bauriid therocephalians and cynodonts, a bony secondary palate is fully developed (Abdala et al., 2014). It is noteworthy that a complete bony secondary palate is also present in dicynodonts, however it is primarily formed by the premaxilla (King, 1988) and not the maxilla as documented in therocephalians, cynodonts and extant mammals. Although a secondary osseous palate is ubiquitous in mammals, it also occurs in a few ectotherms (crocodiles, scincid lizards), thus questioning its direct link to endothermy (Bennett and Ruben, 1986).

Almost all extant endotherms possess nasal turbinate bones covered with mucosa that reduce heat loss and moisten air during respiration (Owerkowicz et al., 2015). This feature, absent in extant ectotherms (Witmer, 1995), may have been present in therocephalian, cynodont and dicynodont therapsids, as postulated from bony ridges in the nasal cavities interpreted as supports for the turbinate complex (Hillenius, 1992, 1994; Crompton et al., 2017).

A peculiar histological structure of fast-growing bone associated with highly vascularised woven-fibred matrix and primary osteons known as fibrolamellar bone (FLB), is another feature often used as evidence of a high metabolic activity (Montes et al., 2010; Legendre et al., 2016). Accordingly, several bone palaeohistological studies have addressed the quest for the presence of FLB in therapsids (de Ricqlès, 1972, 1979; Botha, 2003; Botha and Chinsamy, 2001, 2004; Ray et al., 2004; Olivier et al., 2017). Ray et al. (2004) and Olivier et al. (2017) analysed several therapsid groups (anomodont, gorgonopsian, therocephalian, cynodont) and found FLB in some genera (Aelurognathus, Pristerognathus, Trinitylodon, Oudenodon, Lystrosaurus, Moghreberia), suggesting sustained fast growth, and thus elevated metabolic activity. The presence of FLB has also been demonstrated in some earlier non-therapsid synapsids such as Sphenacodon, Dimetrodon or even Ophiacodon (Huttenlocker et al., 2010; Shelton et al., 2012; Shelton and Sander, 2017). However, FLB also occurs in a few ectotherms such as in some turtles and crocodilians, and is absent in small mammals and passerine birds (Bouvier, 1977), showing that FLB is mostly correlated with high growth rates, which does not always correlate to high metabolic rates. Therefore, these characters alone cannot be considered as definitive evidence of endothermy, leaving the question of the timing and selection pressure for the origin of mammal endothermy still heavily debated.

Because the oxygen isotope fractionation between bone or tooth phosphate and body fluids is temperature dependent, and phosphate has a strong resistance to diagenetic alteration, oxygen isotope compositions of therapsid apatite phosphate (δ¹⁸Oₚ) has been used in this pilot study to investigate the origin of mammalian endothermy. Indeed the δ¹⁸Oₚ value of vertebrate apatite (bone, tooth) reflects both the oxygen isotope composition of the animal body water (δ¹⁸Oᵣₗ) and its body temperature (Tᵣ). Body water derives mainly from drinking meteoric water or plant water (D'Angela and Longinelli, 1990; Kohn, 1996a), and the δ¹⁸O value of this water in turn depends on climatic parameters such as air temperature, hygrometry, and amount of precipitation (Dansgaard, 1964; von Grafenstein et al., 1996; Fricke and O'Neil, 1999).
Variations in the δ¹⁸O values of ectotherm apatite, along with increasing latitude, are expected to reflect decreasing air temperatures as their body temperatures follow those of the environment. In contrast endotherms, which have a constant body temperature, should not be affected by environmental temperatures changes. Moreover, physiological adaptation to specific habitat use (aquatic, semi-aquatic or terrestrial) affects the δ¹⁸O<sub>bw</sub> value by controlling the magnitude of body input and output oxygen fluxes, some of them being associated with oxygen isotopic fractionations (*Amiot et al., 2010*). Therefore, co-existing endotherms and ectotherms should have distinct apatite δ¹⁸O<sub>p</sub> values reflecting their body temperature and ecological differences. By comparing apatite δ¹⁸O<sub>p</sub> values of therapsids with those of co-existing ectotherms of known ecologies at various palaeolatitudes, it should be possible to infer therapsid thermophysiology, a methodology that has previously been applied to non-avian dinosaurs (*Fricke and Rogers, 2000; Amiot et al., 2006*).

Following the protocol of previous research undertaken to establish the Permo-Triassic climatic conditions that prevailed during which South African tetrapods, including therapsids, radiated (*Rey et al., 2016*), this study aims to investigate thermophysiological strategies developed by various Permo-Triassic therapsid groups using the stable oxygen isotope compositions of their phosphatic remains. Our results add new data to the discussion of the origin of mammalian endothermy and its link to global climatic change.

## Results

### Permian therapsids

The 13 sampled South African Permian therapsids come from three different assemblage zones (AZ) of the Beaufort Group: the lower *Tapinocephalus* AZ, the *Tropidostoma* AZ and the lower *Daptocephalus* AZ (*Viglietti et al., 2016*).

Oxygen isotope compositions of three therapsid genera (*Dicynodon*, *Diictodon* and *Oudenodon*) from the youngest assemblage zone (lower *Daptocephalus* AZ) and seven therapsid genera (*Aelurosaurus*, *Diictodon*, *Ictidosuchoides*, *Oudenodon*, *Rhachiocephalus*, *Tropidostoma* and a basal cynodont) from the *Tropidostoma* AZ were respectively compared with one co-occuring *Rhinesuchus* and one co-occuring rhinesuchid. Differences between all the Permian therapsids and ectothermic temnospondyl range from +1.1 ± 0.6‰ to +8.0 ± 0.9‰ (*Figure 1A*), encompassing the expected range for which therapsids are considered ectothermic.

In addition, δ¹⁸O<sub>p</sub> values of the therapsids *Dicynodon*, *Diictodon* and *Oudenodon* are compared to those of the supposedly semi-aquatic parareptile *Pareiasaurus* (*Ivakhnenko, 2001; Kriloff et al., 2008*) (see Appendix 1), with an observed range of +4.3 ± 0.4‰ to +6.8 ± 0.5‰ (*Figure 1A*) which is similar to that measured between therapsids and amphibians. This also supports the ectothermic status of *Dicynodon*, *Diictodon* and *Oudenodon*.

Anteosaurus, *Criocephalosaurus*, *Struthiocephalus*, *Glanosuchus* and a titanosuchid from the lower *Tapinocephalus* AZ have also been compared to two co-occuring basal pareiasaurs which are attributed to either *Embrithosaurus*, *Nochelosaurus* or *Bradysaurus* (*Lee, 1997*) and are considered to have been terrestrial animals (*Canoville et al., 2014*). From *Figure 1B*, the δ¹⁸O<sub>p</sub> differences range from -1.4 ± 0.6‰ to 0.7 ± 1.0‰, also supporting the ectothermic status of these therapsids.

From the Middle Permian of China, one anteaurid *Sinophoneus yumenensis* from the low palaeolatitude locality of Dashankou has a δ¹⁸O<sub>p</sub> value 4.4 ± 0.3‰ lower than the co-existing bolosaurid parareptile *Belebey chengi*, which is considered to have been a terrestrial ectotherm (*Berman et al., 2000; Müller et al., 2008*). This difference between only two values would suggest that *Sinophoneus* was endothermic, but it is also very close to the expected ranges for ectothermic therapsids (*Figure 1B*). Considering *Sinophoneus* as semi-aquatic, as has been suggested for some anteaurids (*Boonstra, 1955, Boonstra, 1962*), the δ¹⁸O<sub>p</sub> difference would imply an ectothermic thermophysiology for this therapsid. This hypothesis needs to be tested with a larger number of samples, which are not yet available.

### Early to Middle Triassic therapsids

From the *Cynognathus* AZ (subzone B) of South Africa, differences between the therapsids *Kanne-meyeria*, *Cynognathus* and *Diademodon* and the temnospondyl amphibians *Xenotosuchus* and *Microposaurus* range from -1.5 ± 1.1‰ to +0.9 ± 1.5‰ (*Figure 2A*), which fit within the range
<figure><img src="image_3.png" /><figcaption>**Figure 1.** δ¹⁸Oₚ differences between Permian therapsids and other tetrapods. Differences in δ¹⁸Oₚ values between therapsids and stereospondyls (white symbols) and between therapsids and parareptiles (black symbols) from the same localities are plotted against their corresponding palaeolatitude. A theoretical framework based on modern temperature gradient (0.6 ± 0.1°C/°Lat; see Appendix 1) and phosphate-water-temperature Figure 1 continued on next page</figcaption></figure>
Figure 1 continued

oxygen isotope fractionation (Lécuyer et al., 2013) predicts various δ¹⁸Oₚ value differences. The lighter orange and red areas correspond to the uncertainty of the temperature gradient. Dap.: Daptocephalus; Tap.: Tapinocephalus.

DOI: 10.7554/eLife.28589.003

predicting endothermic therapsids. Interestingly, these therapsids have values ranging from 0.0 ± 1.6‰ to +1.8 ± 1.6‰ higher than the coexisting terrestrial archosauriform *Erythrosuchus* (*Botha-Brink and Smith*, 2011), a range suggesting that they shared a similar thermophysiology (Figure 2B). Therefore δ¹⁸Oₚ values imply that, as in the case of the therapsids, *Erythrosuchus* was also endothermic which is consistent with the elevated growth rates implied by its palaeohistology (de Ricqlès et al., 2008; Botha-Brink and Angielczyk, 2010).

Also from South Africa, five *Lystrosaurus* specimens from the lower *Lystrosaurus* AZ have δ¹⁸Oₚ values similar to those of the co-existing semi-aquatic stereospondyl *Lydekkerina* (Schoch, 2008; Canoville and Chinsamy, 2015). In addition, an indeterminate lystrosaurid from the Induan Jiu-caiyuan Formation of the Xinjiang Province has a δ¹⁸Oₚ value similar (with a difference of -0.1 ± 0.6‰; *Figure 2B*) to that of the proterosuchid 'Chasmatosaurus' yuani, a basal archosauri-form considered terrestrial and possessing an intermediate thermometabolsim based on a palaeo-histological study (Botha-Brink and Smith, 2011). The combined results from South Africa and China suggest that the analysed lystrosaurids were terrestrial endotherms (*Figure 2*; see Appendix 1).

From the Ermaying Formation of the Shanxi Province (China), the therapsids *Shansiodon wangi* and *Parakannemeyeria youngi* have respectively δ¹⁸Oₚ values of 2.0 ± 0.7‰ and 1.7 ± 0.7‰. These are both lower than the sampled erythrosuchid archosauriform *Shansisuchus shansisuchus*, which fall within two theoretical overlapping ranges (*Figure 2B*). As for the South African erythrosuchids, if we consider *Shansisuchus* as a terrestrial endotherm-like animal and the low palaeolatitude of this part of China region, then the two therapsids also fall within the range of endotherms.

## Middle to Late Triassic therapsids

The late Anisian cynodont *Diademodon* and the kannemeyeriform *Kannemeyeria*, from the Cynog-nathus AZ (subzone C), have both lower δ¹⁸Oₚ values than those of the contemporary semi-aquatic stereospondyls *Paracyclotosaurus* and *Xenotosuchus* with differences ranging from -3.9 ± 2.7‰ to -0.5 ± 0.6‰ (*Figure 3A*). This pattern fits within the main range predicting endothermic therapsids.

The Moroccan kannemeyeriform *Moghreberia nmachouensis* from the early middle Carnian of the Argana Basin has a mean δ¹⁸Oₚ value 2.0 ± 0.5‰ higher than the co-existing aquatic stereospondyl *Almasaurus habbazi* (*Figure 3A*), thus implying that *Moghreberia nmachouensis* was also endothermic.

An indeterminate cynodont from the Rhaetian Lower Elliot Formation of Lesotho has a δ¹⁸Oₚ value 2.1 ± 0.3‰, higher than that of an indeterminate basal sauropodomorph. The suspected endothermy and terrestriality of both dinosaurs (*Amiot et al., 2006; D'Emic, 2015*) and cynodonts are in agreement with their δ¹⁸Oₚ difference that falls within the expected range predicting similar thermophysiology between the two (*Figure 3B*).

## Discussion

According to the δ¹⁸Oₚ value differences observed between therapsids and co-existing non-therap-id tetrapods, elevated thermometabolism seems to have been acquired by at least two therapsid clades: the unnamed dicynodont clade comprising *Lystrosauridae* + *Kannemeyeriformes*, abbreviated the 'L+K' clade, and the *Eucynodontia* (*Figure 4*).

Among the interpreted endothermic therapsids, six belong to the L+K clade (Figure 4): *Moghreberia nmachouensis*, *Parakannemeyeria youngi*, *Kannemeyeria simocephalus* and *Shansiodon wangi* belong to the *Kannemeyeriformes* clade, whereas the *Lystrosauridae* clade comprises *Lystrosaurus* and the Chinese indeterminate *lystrosaurid*. The interpretation of these taxa as endothermic-like animals is better supported for the African *Kannemeyeriformes*, *M. nmachouensis* and *K. simocephalus*, where more individuals were sampled, than for the Chinese species, *S. wangi* and *P. youngi*. Concerning the *lystrosaurids*, *Viglietti et al. (2013)* demonstrated aggregating behaviour in the Early
<figure><img src="image_4.png" /><figcaption>**Figure 2.** δ¹⁸Oₚ differences between Early to Middle Triassic therapsids and other tetrapods. Differences in δ¹⁸Oₚ values between therapsids and stereospondyls (white symbols) and between therapsids and archosauriforms (black symbols) from the same localities are plotted against their corresponding palaeolatitude. A theoretical framework based on a lower-than-today thermal gradient (0.4 ± 0.1°C/°Lat; see Appendix 1) and *Figure 2 continued on next page*</figcaption></figure>
Figure 2 continued

phosphate-water-temperature oxygen isotope fractionation (Lécuyer et al., 2013) predicts various $\delta^{18}$O$_p$ value differences. The lighter orange and red areas correspond to the uncertainty of the temperature gradient.

DOI: 10.7554/eLife.28589.004

Triassic *L. declivis* and interpreted this as a means to keep warm under extreme climatic conditions. This is in agreement with our endothermic interpretation for the genus.

Based on the above interpretations, *Dicynodon* seems to have been ectothermic, a fact which would suggests the rise of endothermy amongst the dicynodont L+K clade during the latest Permian (Lopingian). This can be investigated in the future through the isotopic study of a basal dicynodon-toid such as *Daptocephalus* from the Lopingian of South Africa (Kammerer et al., 2011; Viglietti et al., 2016) or *Peramodon* from the Salarevo Formation of Russia (Sushkin, 1926).

Based on our interpretations, the monophyletic group Eucynodontia (Ruta et al., 2013) (including *Cynognathus*, *Diademodon*, an indeterminate cynodont from Lesotho, and extant mammals) possessed endothermic-like metabolism. A parsimonious interpretation would imply rooting the origin of mammal endothermy within non-Eucynodontia Epicynodontia, between the end-Permian and the earliest Triassic. According to our results, the closest sampled relative of Eucynodontia, the late Permian basal cynodont (SAM-PK-K05339), was probably ectothermic. Therefore the origin of mammal endothermy could have taken place among 'intermediate' groups belonging to the Epicynodontia (such as *Cynosaurus* from the latest Permian or *Thrinaxodon* from the Early Triassic of South Africa) that have in the past been considered to have been endothermic, based on anatomical features (Hillenius and Ruben, 2004; Benoit et al., 2015, 2016a).

In agreement with recent phylogenies (Kemp, 2012; Ruta et al., 2013), endothermic-like body temperature regulations seem to root sometime during the late Permian (Lopingian) independently within the L+K and Epicynodontia clades, the latter being at the origin of mammal endothermy.

An alternative hypothesis considers that both the L+K and Epicynodontia clades possessed a homologous endothermy inherited from their direct common ancestors, the basal Neotherapsida (Figure 4). This suggests that biochemical and physiological mechanisms enabling mammal endothermy, appeared at the base of the neotherapsids during the middle Permian, which is a conclusion recently published based on the paleohistology of some dicynodonts (Olivier et al., 2017). In our case, the absence of an endothermic signal in the sampled Permian therapsids could be due to an endothermy being only seasonal, linked to the presence of a cold season or to the reproduction period (as observed today in some reptile species; Tattersall et al., 2016), which would not be visible in a bulk signal. Therefore, effective acquisition of mammal 'true endothermy' was expressed independently within these two lineages, possibly as a result of extrinsic factors.

Global and regional palaeoclimate reconstructions show a cooling trend toward the end-Permian, followed by an abrupt and intense warming at the Permian-Triassic Boundary (Chen et al., 2013; Rey et al., 2016). Interestingly, most of the therapsid clades which survived the end-Permian mass extinction were supposedly endothermic. It thus appears that climatic fluctuations may have acted as selective pressures which favoured or 'activated' elevated thermometabolic capabilities within therapsids, at the origins of mammal endothermy. A possible explanation could be the acquisition of a fast growth rate due to the high metabolic rate of the endothermy. According to a recent palaeohistology study (Botha-Brink et al., 2016), Early Triassic therapsids, such as *Lystrosaurus* or even therocephalians and cynodonts, had a high growth rate allowing them to reach reproductive maturity within a few seasons and compensate their shortened life expectancies. This adaptation might have enabled certain therapsids to survive the intense climatic change of that time and conquer the newly vacant niches.

## Concluding remarks

In order to investigate the origin of mammal endothermy amongst the Permo-Triassic therapsids, stable oxygen isotope compositions of apatite phosphate and carbonate from therapsids and associated taxa recovered from several palaeolatitudes were analysed. The following results are highlighted:
<figure><img src="image_5.png" /><figcaption>**Figure 3.** δ¹⁸Oₚ differences between Middle to latest Triassic therapsids and other tetrapods. Differences in δ¹⁸Oₚ values between therapsids and stereospondyls (white symbols) and between therapsids and archosauriforms (black symbols) from the same localities are plotted against their corresponding palaeolatitude. A theoretical framework based on a lower-than-today thermal gradient (0.5 ± 0.1°C/°Lat; see Appendix 1) and *Figure 3 continued on next page*</figcaption></figure>
Figure 3 continued

phosphate-water-temperature oxygen isotope fractionation (Lécuyer et al., 2013) predicts various δ¹⁸Oₚ value differences. The lighter orange and red areas correspond to the uncertainty of the temperature gradient.

DOI: 10.7554/eLife.28589.005

* Assuming that analysed samples have preserved their original isotope composition of phosphate, all the Permian therapsids analysed appear to have ectotherm-like thermoregulation and representatives of two Triassic therapsid clades are considered to have had endotherm-like thermoregulation: the Lystrosauridae + Kannemeyeriformes and the Eucynodontia.

* It is proposed that constant elevated thermometabolism appeared independently, at least twice during therapsid evolution. Following the principles of parsimony and phylogenetic systematics, both evolutionary events occured during the late Permian.

* It seems that the timing of the acquisition of elevated thermometabolism among amniotes coincides with major global climatic and environmental fluctuations at the Permo-Triassic boundary and may had a selective advantage to survive the extinction event and result ultimately in mammalian endothermy.

# Material and methods

## Sample collection

Nineteen new fossil apatite samples were analysed to determine stable oxygen isotope compositions of apatite phosphate and carbonate, along with 89 samples for which oxygen isotope compositions have already been published (Rey et al., 2016; Supplementary file 1). This sample total comprises 41 teeth and 65 bones of 90 individual tetrapods (Therapsida, Archosauriformes, Parareptilia and Stereospondyli) recovered from Permian and Triassic deposits of South Africa, Lesotho, Morocco and China. All the sample localities are correlated to the marine biostratigraphic stages using the absolute ages accepted by the International Commission on Stratigraphy (Cohen et al., 2013); updated 12/2016), with the Permo-Triassic and Guadalupian-Lopingian boundaries now respectively considered to be at 251.90 ± 0.02 Ma (Burgess et al., 2014) and 259.1 ± 0.5 (Zhong et al., 2014) Ma.

South African samples comprise Permian and Triassic bones and teeth of therapsids, pareiasaurs, archosauriforms and stereospondyls recovered from 10 localities in the Beaufort Group (Karoo Supergroup), and housed in the collections of the Iziko South African Museum, Cape Town (SAM, Supplementary file 1) and at the Evolutionary Studies Institute, University of the Witwatersrand, Johannesburg (ESI, Supplementary file 1). Permian biozone ages of South African localities were taken from (Rubidge et al., 2013; Day et al., 2015), whereas Triassic age determination has been achieved by biostratigraphic correlation with Laurasian sequences (Hancox et al., 1995; Rubidge, 2005; Abdala and Ribeiro, 2010).

Lesotho samples comprise a cynodont therapsid and a basal sauropodomorph dinosaur from a Triassic locality near the town of Pokane, and are part of the Paul Ellenberger Collection at the Institut des Sciences de l'Evolution, University Montpellier, France (ISEM, Supplementary file 1). The locality belongs to the 'Red Beds inférieurs a or b' of the lower Elliot Formation which is currently regarded as latest Triassic (late Rhaetian) (Knoll, 2004).

Moroccan samples comprise therapsid and stereospondyl bones recovered from the 'Locality 11' of the Argana Group (Jalil, 1999) near the village of Alma, and housed at the Museum National d'Histoire Naturelle, Paris, France (MNHN, Supplementary file 1). The locality is biostratigraphically correlated to the upper Timezgadiouine Formation, considered to be Middle to early Late Carnian (Jalil, 1999).

Chinese samples are from Permian and Triassic localities situated in Gansu, Shanxi and Xinjiang provinces and comprise therapsids found in association with archosauriforms or parareptiles. These remains are curated at the Institute of Vertebrate Paleontology and Paleoanthropology in Beijing, China (IVPP, Supplementary file 1). The Dashankou locality, from Gansu Province, is biostratigraphically dated as Early Roadian (Liu et al., 2009; Liu, 2010). From Shanxi Province, sampled fossils originate from three localities in the Ermaying Formation which is considered to be Anisian (Liu et al.,
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Phylogeny of sampled therapsids. Phylogeny of the sampled therapsids plotted alongside a stratigraphic scale, based on proposed therapsid phylogenies (Ruben and Jones, 2000; Hillenius and Ruben, 2004; Gebauer, 2007; Cisneros et al., 2012; Liu, 2013; Ruta et al., 2013) and their biostratigraphic ranges (Kammerer et al., 2011; Kemp, 2012; Ruta et al., 2013; Huttenlocker, 2014; Day et al., 2015; Viglietti et al., 2016). The thickest parts of the bold lines represent the age range uncertainty of the localities where the samples come from. Species identified as endotherm-like *Figure 4 continued on next page*</figcaption></figure>
Figure 4 continued

are written in bold and red. Node numbers refer to clades quoted in the text: 1: Neotherapsida; 2: Dicynodontoidea; 3: Lystrosauridae; 4: Kannemeyeriformes; 5: Epicynodontia; 6: Eucynodontia.

DOI: 10.7554/eLife.28589.006

2013). From Xinjiang Province, two localities in the Jiucaiyuan Formation have been sampled and are considered Early Triassic (Metcalfe et al., 2009).

Calculation of palaeogeographic coordinates of the sampling sites was performed after careful selection of the magnetic poles of West Gondwana (Muttoni et al., 2001), North China (et al., 1992), South Jungar (Choulet et al., 2013) and the Alashan terrane (Meng, 1992; Yuan and Yang, 2015). The Apparent Polar Wander Path (APWP) of South Africa (Torsvik et al., 2012) was used to

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Isotopic preservation assessment. δ¹⁸O_c-δ¹⁸O_p differences between teeth and bones plotted against the structural carbonate content (wt%) of apatite. Samples that have δ¹⁸O_c-δ¹⁸O_p differences higher than 14.7‰ are considered doubtful as regards potential diagenetic alteration (see text). For carbonate contents (wt%) higher than 13.4%, the δ¹⁸O_c values are considered to be inherited from inorganic diagenetic processes. A high difference between δ¹⁸O_c and δ¹⁸O_p is interpreted as the result of a microbially-mediated alteration of the apatite phosphate or too high δ¹⁸O_c values resulting from the addition of inorganic carbonate or isotopic exchange with an external source of inorganic carbon. The grey crosses refer to previously published South African bone and tooth samples (Rey et al., 2016).</figcaption></figure>

DOI: 10.7554/eLife.28589.007
constrain the palaeolatitudinal position of the South African and Lesotho fossil sites. Palaeolatitudes and associated uncertainties (A95) are shown in **Supplementary file 1**.

## Analytical techniques

To measure the oxygen isotope composition of the apatite phosphate group, the phosphate ions were isolated using acid dissolution and anion-exchange resin applying a standard protocol (Lécuyer, 2004). Silver phosphate was quantitatively precipitated in a thermostatic bath set at a temperature of 70°C. After filtration, washing with double deionized water and drying at 50°C, an aliquot of 300 µg of Ag₃PO₄ was mixed with 300 µg of nickelised carbon in a silver reaction capsule. Silver phosphate was then reduced into CO to measure its ¹⁸O/¹⁶O ratio (Lécuyer et al., 2007; Fourel et al., 2011). Each sample was heated at 1450°C by pyrolysis using a VarioPYROcube EA system (Elementar) interfaced to an IsoPrime isotope ratio mass spectrometer working in continuous flow mode at the UMR CNRS 5276 LGLTPE, University Claude Bernard Lyon 1.

Isotopic compositions are quoted in the standard δ notation relative to V-SMOW. Silver phosphate precipitated from standard NBS120c (natural Miocene phosphorite from Florida) was repeatedly analysed (δ¹⁸O = 21.71 ± 0.20‰; n = 30) along with the silver phosphate samples derived from the tetrapod remains. For the oxygen isotope analysis of apatite carbonate, about 10 mg of tooth or bone powder was pre-treated (Koch et al., 1997). Powders were washed with a 2% NaOCl solution to remove organic matter, then rinsed five times with double deionized water and air-dried at 40°C for 24 hr. Potential secondary carbonate was removed by adding 0.1 M acetic acid and leaving for 24 hr, after which the powder was again rinsed five times with double deionized water and air-dried at 40°C for 24 hr. The powder/solution ratio was kept constant at 0.04 g mL⁻¹ for both treatments. Stable isotope ratios were determined by using a Thermo Finnigan Gasbench II at the geochemistry laboratory of the Institute of Geology and Geophysics (Chinese Academy of Sciences, China). For each sample, an aliquot of 2 mg of pre-treated apatite was reacted with 5 drops of supersaturated orthophosphoric acid at 72°C for one hour under a He atmosphere before starting 10 measurement cycles of the isotopic composition of the CO₂ produced with a Finnigan MAT 253 continuous flow isotope ratio mass spectrometer. The measured oxygen isotopic compositions were normalized relative to the NBS-19 calcite standard and have a reproducibility index better than ±0.2‰. Isotopic compositions are quoted in the standard δ notation relative to V-SMOW.

## Robustness of the stable isotope record

Analysed materials consist of bone or tooth dentine, which is more porous than enamel with small and less densely inter-grown apatite crystals (Mills, 1967). Thus, their original stable isotope compositions are more prone to diagenetic alteration that may have taken place through precipitation of secondary minerals within and at the surface of bioapatite crystals, adsorption of ions on the surface of apatite crystals, or dissolution and recrystallization with isotopic exchange. The samples from South Africa have been previously tested for primary preservation through comparison between their δ¹⁸Oₚ values, δ¹⁸O_c values and carbonate content on the basis of the following considerations: (1) the carbonate content in apatite of modern vertebrates typically ranges from less than 1% up to 13.4%. Thus, samples that have a carbonate content exceeding 13.4 wt% likely contain additional inorganic carbonate precipitated from diagenetic fluids, and would result in potentially biased δ¹⁸O_c values of apatite carbonate (Figure 5); (2) In modern vertebrates, the oxygen isotope composition of apatite carbonate is higher than that of co-occurring apatite phosphate (7–9 ‰ in mammals), and up to 14.7‰ in sharks (Vennemann et al., 2001). Experimental (et al., 1967) and empirical studies (Zazzo et al., 2004b) have shown that microbially-mediated diagenetic alteration of apatite phosphate results in a greater difference between δ¹⁸O_c and δ¹⁸Oₚ values. Therefore, fossil samples exhibiting δ¹⁸O_c-δ¹⁸Oₚ differences larger than 14.7‰ are most likely altered and can be disregarded (Figure 5). Inorganic alteration at low temperature has little effect on the δ¹⁸Oₚ values of phosphates, even at geological time scales (Lecuyer et al., 1999), so samples affected by inorganic diagenetic alteration of carbonates, (resulting either in a high overall carbonate content or anomalous δ¹⁸O_c-δ¹⁸Oₚ differences), may still preserve the original oxygen isotope composition of their phosphate (Figure 5). Using these two assessments, newly measured δ¹⁸Oₚ values are considered to have preserved their original isotopic signatures and can be interpreted in terms of ecologies and physiologies.
## Assessment of therapsid thermophysiology

For all localities, average δ¹⁸Oₚ values were calculated for each tetrapod species. Differences in δ¹⁸Oₚ values between therapsid species and co-occurring non-therapsid tetrapods (amphibians, parareptiles or archosauriforms) were calculated and plotted against their corresponding palaeolatitude for three time intervals: the middle to late Permian (Figure 2), the Early to Middle Triassic (Figure 3) and the Middle Triassic to latest Triassic (Figure 4). These differences were compared to the following four theoretical areas of values represented as coloured areas in Figures 1–3. To construct those theoretical areas, both the phosphate-water temperature scale from Lécuyer et al., 2013 and the differences of stable oxygen compositions between mammals of various ecologies from Cerling et al. (2008) have been used. (see Appendix 1 for their construction details).

Orange and green areas in Figures 1A, 2A and 3A represent expected δ¹⁸Oₚ value differences between terrestrial therapsids and semi-aquatic stereospondyls (white symbols) or parareptiles (black symbols); red and blue areas in Figures 1B, 2B and 3B represent expected δ¹⁸Oₚ value differences between terrestrial therapsids and terrestrial Permian parareptiles or Triassic archosauriforms (black symbols). Oblique orange and red areas in Figures 1–3 represent expected δ¹⁸Oₚ value differences between an endotherm and an ectotherm. Vertical green and blue areas in Figures 1–3 represent expected δ¹⁸Oₚ value differences between animals having similar thermophysiology.

## Acknowledgements

The authors thank the MNHN, the IVPP, Iziko SA Museum and the ESI for granting access to the fossils, and the South African Heritage Resources Agency (SAHRA) for their authorizations to sample and export fossils for stable isotope analysis (PermitID: 1858). We also acknowledge S Jiquel (ISEM) for the access to the Ellenberger's collection of the University Montpellier 2, J Falconnet (MNHN) for his identification, L. Cui (IGGCAS) for laboratory help and J Cubo (UPMC) for constructed discussions that improved the manuscript. This work was supported by a French project INSU of the CNRS, the Palaeontological Scientific Trust (PAST) and its Scatterlings of Africa programmes, National Research Foundation (NRF) and DST/NRF Centre of Excellence in Palaeosciences. AR and WX were supported by the National Basic Research Program of China (2012CB821900), and LJ by the National Basic Research Program of China (2012CB821902). We thank the four reviewers, K Angielczyk, J Eiler and two anonymous ones, for their constructive comments.

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Institut National des Sciences de l'Univers, Centre National de la Recherche Scientifique</td><td></td><td>Romain Amiot</td></tr><tr><td>Palaeontological Scientific Trust</td><td></td><td>Bruce S Rubidge</td></tr><tr><td>National Research Foundation</td><td></td><td>Bruce S Rubidge</td></tr><tr><td>Centre of Excellence in Palaeosciences</td><td></td><td>Bruce S Rubidge</td></tr><tr><td>National Basic Research Program of China</td><td>2012CB821900</td><td>Romain Amiot<br>Xu Wang</td></tr><tr><td>National Basic Research Program of China</td><td>2012CB821902</td><td>Jun Liu</td></tr><tr><td>Institut Universitaire de France</td><td></td><td>Christophe Lécuyer</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
## Author contributions

KR, Conceptualization, Data curation, Formal analysis, Investigation, Writing—original draft, Project administration, Writing—review and editing; RA, Conceptualization, Data curation, Formal analysis, Supervision, Project administration, Writing—review and editing; FFo, Conceptualization, Data curation, Project administration, Writing—review and editing; FA, N-EJ, JL, RMHS, JSS, PAV, Resources, Writing—review and editing; FFI, Formal analysis, Investigation, Writing—review and editing; BSR, Resources, Funding acquisition, Project administration, Writing—review and editing; XW, Resources, Data curation, Formal analysis, Writing—review and editing; CL, Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Project administration, Writing—review and editing

## Author ORCIDs

Kévin Rey, http://orcid.org/0000-0002-6788-2453

## Additional files

### Supplementary files

• Supplementary file 1. Stable oxygen isotope compositions of phosphate (δ¹⁸O~P~) and carbonate (δ¹⁸O~C~) of Permo-Triassic tetrapod teeth and bones reported along with their stratigraphic position, estimated age, palaeolatitudes and their carbonate content. Asterisks represent diagenetically altered values.

DOI: 10.7554/eLife.28589.008

## References

Abdala F, Jashashvili T, Rubidge BS, Heever J. 2014. New Material of *Microgomphodon* *oligocynus* (Eutherapsida, Therocephalia) and the Taxonomy of Southern African Bauriidae. In: Kammerer C. F, Angielczyk K. D, Fröbisch J (Eds). *Early Evolutionary History of the Synapsida, Vertebrate Paleobiology and Paleoanthropology*. Netherlands: Springer. p. 209–231.

Abdala F, Ribeiro AM. 2010. Distribution and diversity patterns of Triassic cynodonts (Therapsida, Cynodontia) in Gondwana. *Palaeogeography, Palaeoclimatology, Palaeoecology* **286**:202–217. doi: 10.1016/j.palaeo.2010.01.011

Amiot R, Buffetaut E, Lecuyer C, Wang X, Boudad L, Ding Z, Fourel F, Hutt S, Martineau F, Medeiros MA, Mo J, Simon L, Suteethorn V, Sweetman S, Tong H, Zhang F, Zhou Z. 2010. Oxygen isotope evidence for semi-aquatic habits among spinosaurid theropods. *Geology* **38**:139–142. doi: 10.1130/G30402.1

Amiot R, Lécuyer C, Buffetaut E, Escarguel G, Fluteau F, Martineau F. 2006. Oxygen isotopes from biogenic apatites suggest widespread endothermy in cretaceous dinosaurs. *Earth and Planetary Science Letters* **246**:41–54. doi: 10.1016/j.epsl.2006.04.018

Amiot R, Lécuyer C, Buffetaut E, Fluteau F. Legendre S, Martineau F. 2004. Latitudinal temperature gradient during the Cretaceous Upper Campanian–Middle Maastrichtian: δ18O record of continental vertebrates. *Earth and Planetary Science Letters* **226**:255–272. doi: 10.1016/j.epsl.2004.07.015

Bennett AF, Ruben JA. 1986. *The Metabolic and Thermoregulatory Status of Therapsids. the Ecology and Biology of Mammal-Like Reptiles*. Washington, DC: Smithsonian Institution Press. p. 207–218.

Benoit J, Abdala F, Manger P, Rubidge B. 2016a. The sixth sense in mammalians forerunners: variability of the parietal foramen and the evolution of the pineal eye in south african Permo-Triassic eutheriodont therapsids. *Acta Palaeontologica Polonica* **61**:777–789. doi: 10.4202/app.00219.2015

Benoit J, Abdala F, Van den Brandt MJ, Manger PR, Rubidge BS. 2015. Physiological implications of the abnormal absence of the parietal foramen in a late Permian cynodont (Therapsida). *The Science of Nature* **102**: 69. doi: 10.1007/s00114-015-1321-4, PMID: 26538062

Benoit J, Manger PR, Rubidge BS. 2016b. Palaeoneurological clues to the evolution of defining mammalian soft tissue traits. *Scientific Reports* **6**:25604. doi: 10.1038/srep25604, PMID: 27157809

Berman DS, Reisz RR, Scott D, Henrici AC, Sumida SS, Martens T. 2000. Early Permian bipedal reptile. *Science* **290**:969–972. doi: 10.1126/science.290.5493.969, PMID: 11062126

Bernal D, Dickson KA, Shadwick RE, Graham JB. 2001. Review: analysis of the evolutionary convergence for high performance swimming in lamnid sharks and tunas. *Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology* **129**:695–726. doi: 10.1016/S1095-6433(01)00333-6, PMID: 11423338

Boonstra LD. 1955. The girdles and limbs of the South African Deinocephalia. *Annals of the South African Museum* **42**:185–326.

Boonstra LD. 1962. The dentition of the titanosuchian dinocephalians. *Annals of the South African Museum* **46**: 57–112.
Botha J, Chinsamy A. 2001. Growth patterns deduced from the bone histology of the cynodonts *diademodon* and *Cynognathus*. *Journal of Vertebrate Paleontology* **20**:705–711. doi: 10.1671/0272-4634(2000)020[0705:GPDFTB]2.0.CO;2

Botha J, Chinsamy A. 2004. Growth and life habits of the triassic cynodont *Trirachodon*, inferred from bone histology. *Acta Palaeontologica Polonica* **49**:619–627.

Botha J, Smith RMH. 2006. Rapid vertebrate recuperation in the Karoo Basin of South Africa following the End-Permian extinction. *Journal of African Earth Sciences* **45**:502–514. doi: 10.1016/j.jafrearsci.2006.04.006

Botha J. 2003. Biological aspects of the Permian dicynodont *Oudenodon* (Therapsida: dicynodontia) deduced from bone histology and cross-sectional geometry. *Palaeontologia Africana* **39**:37–44.

Botha-Brink J, Angielczyk KD. 2010. Do extraordinarily high growth rates in Permo-Triassic dicynodonts (Therapsida, Anomodontia) explain their success before and after the end-Permian extinction? *Zoological Journal of the Linnean Society* **160**:341–365. doi: 10.1111/j.1096-3642.2009.00601.x

Botha-Brink J, Codron D, Huttenlocker AK, Angielczyk KD, Ruta M. 2016. Breeding Young as a Survival Strategy during Earth's Greatest Mass Extinction. *Scientific Reports* **6**:24053. doi: 10.1038/srep24053, PMID: 27044713

Botha-Brink J, Smith RMH. 2011. Osteohistology of the Triassic archosauromorphs *Prolacerta, Proterosuchus, Euparkeria*, and *Erythrosuchus* from the Karoo Basin of South Africa. *Journal of Vertebrate Paleontology* **31**: 1238–1254. doi: 10.1080/02724634.2011.621797

Bouvier M. 1977. Dinosaur haversian bone and endothermy. *Evolution* **31**:449–450. doi: 10.1111/j.1558-5646. 1977.tb01028.x, PMID: 28563235

Burgess SD, Bowring S, Shen SZ. 2014. High-precision timeline for Earth's most severe extinction. *PNAS* **111**: 3316–3321. doi: 10.1073/pnas.1317692111, PMID: 24516148

Cannon B, Nedergaard J. 2004. Brown adipose tissue: function and physiological significance. *Physiological Reviews* **84**:277–359. doi: 10.1152/physrev.00015.2003, PMID: 14715917

Canoville A, Chinsamy A. 2015. Bone microstructure of the Stereospondyl Lydekkerina Huxleyi reveals adaptive strategies to the Harsh Post Permian-Extinction Environment. *The Anatomical Record* **298**:1237–1254. doi: 10.1002/ar.23160, PMID: 25857487

Canoville A, Laurin M. 2010. Evolution of humeral microanatomy and lifestyle in Amniotes, and some comments on palaeobiological inferences. *Biological Journal of the Linnean Society* **100**:384–406. doi: 10.1111/j.1095-8312.2010.01431.x

Canoville A, Thomas DB, Chinsamy A. 2014. Insights into the habitat of Middle Permian pareiasaurs (Parareptilia) from preliminary isotopic analyses. *Lethaia* **47**:266–274. doi: 10.1111/let.12056

Cerling TE, Harris JM, Hart JA, Kaleme P, Klingel H, Leakey MG, Levin NE, Lewison RL, Passey BH. 2008. Stable isotope ecology of the common hippopotamus. *Journal of Zoology* **276**:204–212. doi: 10.1111/j.1469-7998. 2008.00450.x

Chen B, Joachimski MM, Shen S-zhong, Lambert LL, Lai X-long, Wang X-dong, Chen J, Yuan D-xun. 2013. Permian ice volume and palaeoclimate history: Oxygen isotope proxies revisited. *Gondwana Research* **24**:77–89. doi: 10.1016/j.gr.2012.07.007

Choulet F, Chen Y, Cogné J-P, Rabillard A, Wang B, Lin W, Faure M, Cluzel D. 2013. First Triassic palaeomagnetic constraints from junggar (NW China) and their implications for the mesozoic tectonics in Central Asia. *Journal of Asian Earth Sciences* **78**:371–394. doi: 10.1016/j.jseaes.2013.01.023

Cisneros JC, Abdala F, Atayman-Güven S, Rubidge BS, Şengörc AM, Schultz CL. 2012. Carnivorous dinocephalian from the Middle Permian of Brazil and tetrapod dispersal in Pangaea. *PNAS* **109**:1584–1588. doi: 10.1073/pnas.1115975109, PMID: 22307615

Cohen KM, Finney SC, Gibbard PL, Fan J-X. 2013. The ICS International Chronostratigraphic Chart. *Episodes* **36**: 199–204.

Crompton AW, Owerkowicz T, Bhullar B-AS, Musinsky C. 2017. Structure of the nasal region of non-mammalian cynodonts and mammaliaforms: speculations on the evolution of mammalian endothermy. *Journal of Vertebrate Paleontology* **37**:e1269116. doi: 10.1080/02724634.2017.1269116

D'Angela D, Longinelli A. 1990. Oxygen isotopes in living mammal's bone phosphate: Further results. *Chemical Geology: Isotope Geoscience Section* **86**:75–82.

D'Emic MD. 2015. Comment on "Evidence for mesothermy in dinosaurs". *Science* **348**:982. doi: 10.1126/science. 1260061, PMID: 26023130

Daniel Bryant J, Froelich PN. 1995. A model of oxygen isotope fractionation in body water of large mammals. *Geochimica Et Cosmochimica Acta* **59**:4523–4537. doi: 10.1016/0016-7037(95)00250-4

Dansgaard W. 1964. Stable isotopes in precipitation. *Tellus* **16**:436–468. doi: 10.3402/tellusa.v16i4.8993

Day MO, Ramezani J, Bowring SA, Sadler PM, Erwin DH, Abdala F, Rubidge BS. 2015. When and how did the terrestrial mid-Permian mass extinction occur? evidence from the tetrapod record of the Karoo Basin, South Africa. *Proceedings of the Royal Society B: Biological Sciences* **282**:20150834. doi: 10.1098/rspb.2015.0834, PMID: 26156768

de Ricqlès A, Padian K, Knoll F, Horner JR. 2008. On the origin of high growth rates in archosaurs and their ancient relatives: complementary histological studies on triassic archosauriforms and the problem of a "phylogenetic signal" in bone histology. *Annales De Paléontologie* **94**:57–76. doi: 10.1016/j.annpal.2008.03.002

de Ricqlès A. 1972. Recherches Paléohistologiques Sur Les Os Longs Des Tétrapodes: Titanosuchiens, Dinocéphales Et Dicynodontes. Masson.

de Ricqlès A. 1979. Quelques remarques sur l'histoire évolutive des tissus squelettiques chez les Vertébrés et plus particulièrement chez les Tétrapodes. *Année Biologique* **18**:1–35.
de Ricqlès AJ, Padian K, Horner JR. 2003. On the bone histology of some triassic pseudosuchian archosaurs and related taxa. *Annales De Paléontologie* **89**:67-101. doi: 10.1016/S0753-3969(03)00005-3

Farmer CG, Carrier DR. 2000. Ventilation and gas exchange during treadmill locomotion in the american Alligator (Alligator mississippiensis). *The Journal of Experimental Biology* **203**:1671-1678. PMID: 10804157

Farmer CG. 2000. Parental care: the key to understanding endothermy and other convergent features in Birds and Mammals. *The American Naturalist* **155**:326-334. doi: 10.1086/303323, PMID: 10718729

Fourel F, Martineau F, Lécuyer C, Kupka HJ, Lange L, Ojeimi C, Seed M. 2011. 18o/16O ratio measurements of inorganic and organic materials by elemental analysis-pyrolysis-isotope ratio mass spectrometry continuous-flow techniques. *Rapid Communications in Mass Spectrometry* **25**:2691-2696. doi: 10.1002/rcm.5056, PMID: 21913245

Fricke HC, O'Neil JR. 1999. The correlation between 18O/16O ratios of meteoric water and surface temperature: its use in investigating terrestrial climate change over geologic time. *Earth and Planetary Science Letters* **170**: 181-196. doi: 10.1016/S0012-821X(99)00105-3

Fricke HC, Rogers RR. 2000. Multiple taxon-multiple locality approach to providing oxygen isotope evidence for warm-blooded theropod dinosaurs. *Geology* **28**:799. doi: 10.1130/0091-7613(2000)28<799:MTLATP>2.0.CO;2

Gebauer EV. 2007. Phylogeny and evolution of the Gorgonopsia with a Special Reference to the Skull and Skeleton of GPIT/RE/7113 ('Aelurognathus?'parringtoni). *Universität Tübingen*.

Germain D, Laurin M. 2005. Microanatomy of the radius and lifestyle in amniotes (Vertebrata, Tetrapoda). *Zoologica Scripta* **34**:335-350. doi: 10.1111/j.1463-6409.2005.00198.x

Gower DJ, Hancox PJ, Botha-Brink J, Sennikov AG, Butler RJ. 2014. A new species of *Garjainia* Ochev, 1958 (Diapsida: Archosauriformes: Erythrosuchidae) from the Early Triassic of South Africa. *PLoS One* **9**:e111154. doi: 10.1371/journal.pone.0111154, PMID: 25386937

Grady JM, Enquist BJ, Dettweiler-Robinson E, Wright NA, Smith FA. 2014. Evidence for mesothermy in dinosaurs. *Science* **344**:1268-1272. doi: 10.1126/science.1253143

Hancox PJ, Shishkin MA, Rubidge BS, Kitching JW. 1995. A threefold subdivision of the *Cynognathus* Assemblage Zone (Beaufort Group, South Africa) and its palaeogeographical implications. *South African Journal of Science* **91**:143-144.

Hillenius WJ, Ruben JA. 2004. The evolution of endothermy in terrestrial vertebrates: who? when? why? *Physiological and Biochemical Zoology* **77**:1019-1042. doi: 10.1086/425185, PMID: 15674773

Hillenius WJ. 1992. The evolution of nasal turbinates and mammalian endothermy. *Paleobiology* **18**:17-29. doi: 10.1017/S0094837300012197

Hillenius WJ. 1994. Turbinates in therapsids: evidence for late permian origins of mammalian endothermy. *Evolution* **48**:207-229. doi: 10.1111/j.1558-5646.1994.tb01308.x, PMID: 28568303

Hopson JA. 2012. The role of foraging Mode in the origin of Therapsids: implications for the origin of mammalian endothermy. *Fieldiana Life and Earth Sciences* **5**:126-148. doi: 10.3158/2158-5520-5.1.126

Huttenlocker AK, Rega E, Sumida SS. 2010. Comparative anatomy and osteohistology of hyperelongate neural spines in the sphenacodontids *Sphenacodon* and *Dimetrodon* (Amniota: synapsida). *Journal of Morphology* **271**:1407-1421. doi: 10.1002/jmor.10876, PMID: 20886514

Huttenlocker AK. 2014. Body size reductions in nonmammalian eutheriodont therapsids (Synapsida) during the end-Permian mass extinction. *PLoS One* **9**:e87553. doi: 10.1371/journal.pone.0087553, PMID: 24498335

Ivakhnenko MF. 2001. Tetrapods from the east european placket—Late Paleozoic natural territorial complex. *Tr. Paleontol. Inst. Ross. Akad. Nauk* **283**:1-200.

Jalil N-E. 1999. Continental Permian and Triassic vertebrate localities from Algeria and Morocco and their stratigraphical correlations. *Journal of African Earth Sciences* **29**:219-226. doi: 10.1016/S0899-5362(99)00091-3

Ji Q, Luo ZX, Yuan CX, Tabrum AR. 2006. A swimming mammaliaform from the Middle jurassic and ecomorphological diversification of early mammals. *Science* **311**:1123-1127. doi: 10.1126/science.1123026, PMID: 16497926

Kammerer CF, Angielczyk KD, Fröbisch J. 2011. A comprehensive taxonomic revision of *Dicynodon* (Therapsida, Anomodontia) and its implications for dicynodont phylogeny, biogeography, and biostratigraphy. *Journal of Vertebrate Paleontology* **31**:1-158. doi: 10.1080/02724634.2011.627074

Katz SL. 2002. Design of heterothermic muscle in fish. *The Journal of Experimental Biology* **205**:2251-2266. PMID: 12110659

Kemp TS. 2006a. The origin of mammalian endothermy: a paradigm for the evolution of complex biological structure. *Zoological Journal of the Linnean Society* **147**:473-488. doi: 10.1111/j.1096-3642.2006.00226.x

Kemp TS. 2006b. The origin and early radiation of the therapsid mammal-like reptiles: a palaeobiological hypothesis. *Journal of Evolutionary Biology* **19**:1231-1247. doi: 10.1111/j.1420-9101.2005.01076.x, PMID: 16780524

Kemp TS. 2012. The origin and radiation of Therapsids. In: Chinsamy-Turan A (Ed). *Forerunners of Mammals: Radiation, Histology, Biology*. Bloomington: Indiana University Press. p. 3-28.

Khaliq I, Hof C, Prinzinger R, Böhning-Gaese K, Pfenninger M. 2014. Global variation in thermal tolerances and vulnerability of endotherms to climate change. *Proceedings of the Royal Society B: Biological Sciences* **281**: 20141097. doi: 10.1098/rspb.2014.1097, PMID: 25009066

King GM. 1988. Anomodontia. Stuttgart: Gustav Fischer Verlag GmbH & Company KG.

Knoll F. 2004. Review of the tetrapod fauna of the "Lower Stormberg Group" of the main Karoo Basin (southern Africa): implication for the age of the Lower Elliot Formation. *Bulletin De La Societe Geologique De France* **175**:73-83. doi: 10.2113/175.1.73
Koch PL, Tuross N, Fogel ML. 1997. The effects of Sample treatment and diagenesis on the Isotopic integrity of Carbonate in Biogenic Hydroxylapatite. *Journal of Archaeological Science* **24**:417–429. doi: 10.1006/jasc.1996.0126

Kohn MJ. 1996a. Predicting animal δ18O: Accounting for diet and physiological adaptation. *Geochimica Et Cosmochimica Acta* **60**:4811–4829. doi: 10.1016/S0016-7037(96)00240-2

Kriloff A, Germain D, Canoville A, Vincent P, Sache M, Laurin M. 2008. Evolution of bone microanatomy of the tetrapod tibia and its use in palaeobiological inference. *Journal of Evolutionary Biology* **21**:807–826. doi: 10.1111/j.1420-9101.2008.01512.x, PMID: 18312321

Lecuyer C, Grandjean P, Sheppard SMF. 1999. Oxygen isotope exchange between dissolved phosphate and water at temperatures ≤135°C: inorganic versus biological fractionations. *Geochimica Et Cosmochimica Acta* **63**:855–862. doi: 10.1016/S0016-7037(99)00096-4

Lee MSY. 1997. A taxonomic revision of pareiasaurian reptiles: implications for Permian terrestrial palaeoecology. *Modern Geology* **21**:231–298.

Legendre LJ, Guénard G, Botha-Brink J, Cubo J. 2016. Palaeohistological evidence for ancestral high metabolic rate in Archosaurs. *Systematic Biology* **65**:989–996. doi: 10.1093/sysbio/syw033, PMID: 27073251

Lécuyer C, Amiot R, Touzeau A, Trotter J. 2013. Calibration of the phosphate δ18o thermometer with carbonate-water oxygen isotope fractionation equations. *Chemical Geology* **347**:217–226. doi: 10.1016/j.chemgeo.2013.03.008

Lécuyer C, Fourel F, Martineau F, Amiot R, Bernard A, Daux V, Escarguel G, Morrison J. 2007. High-precision determination of 18O/16O ratios of silver phosphate by EA-pyrolysis-IRMS continuous flow technique. *Journal of Mass Spectrometry* **42**:36–41. doi: 10.1002/jms.1130, PMID: 17149791

Lécuyer C. 2004. Oxygen isotope analysis of phosphate. *Handbook of Stable Isotope Analytical Techniques* **1**: 482–499.

Liu J, Li L, Li XW. 2013. SHRIMP U-Pb zircon dating of the Triassic Ermaying and Tongchuan formations in Shanxi, China and its stratigraphic implications. *Vert PalAs* **51**:162–168.

Liu J, Rubidge B, Li J. 2009. New basal synapsid supports laurasian origin for therapsids. *Acta Palaeontologica Polonica* **54**:393–400. doi: 10.4202/app.2008.0071

Liu J. 2010. Dashankou Fauna: a unique window on the early evolution of therapsids. *Bulletin of the Chinese Academy of Sciences* **24**:83–85.

Liu J. 2013. Osteology, ontogeny, and phylogenetic position of *Sinophoneus yumenensis* (Therapsida, Dinocephalia) from the Middle Permian Dashankou Fauna of China. *Journal of Vertebrate Paleontology* **33**: 1394–1407. doi: 10.1080/02724634.2013.781505

Luz B, Kolodny Y, Horowitz M. 1984. Fractionation of oxygen isotopes between mammalian bone-phosphate and environmental drinking water. *Geochimica Et Cosmochimica Acta* **48**:1689–1693. doi: 10.1016/0016-7037(84)90338-7

McHugh JB. 2014. Paleohistology and histovariability of the Permian stereospondyl *Rhinesuchus*. *Journal of Vertebrate Paleontology* **34**:59–68. doi: 10.1080/02724634.2013.787429

McNab BK, Auffenberg W. 1976. The effect of large body size on the temperature regulation of the Komodo Dragon, *Varanus komodoensis*. *Comparative Biochemistry and Physiology Part A: Physiology* **55**:345–350. doi: 10.1016/0300-9629(76)90058-X, PMID: 9248

McNab BK. 1978. The evolution of Endothermy in the Phylogeny of Mammals. *The American Naturalist* **112**:1–21. doi: 10.1086/283249

Meng ZF. 1992. Paleomagnetic results from the Late Permian redbeds along the Hexi Corridor, China. *Chinese Science Bulletin* = *Kexue Tongbao*:1452–1456.

Metcalfe I, Foster CB, Afonin SA, Nicoll RS, Mundil R, Xiaofeng W, Lucas SG. 2009. Stratigraphy, biostratigraphy and C-isotopes of the Permian–Triassic non-marine sequence at Dalongkou and Lucaogou, Xinjiang Province, China. *Journal of Asian Earth Sciences* **36**:503–520. doi: 10.1016/j.jseaes.2008.06.005

Mills A. 1967. *Structural and Chemical Organization of Teeth*. Elsevier.

Montes L, Castanet J, Cubo J. 2010. Relationship between bone growth rate and bone tissue organization in amniotes: first test of Amprino's rule in a phylogenetic context. *Animal Biology* **60**:25–41. doi: 10.1163/157075610X12610595764093

Muttoni G, Garzanti E, Alfonsi L, Cirilli S, Germani D, Lowrie W. 2001. Motion of Africa and Adria since the Permian: paleomagnetic and paleoclimatic constraints from northern Libya. *Earth and Planetary Science Letters* **192**:159–174. doi: 10.1016/S0012-821X(01)00439-3

Müller J, Li JL, Reisz RR. 2008. A new bolosaurid parareptile, *Belebey chengi* sp. nov., from the Middle Permian of China and its paleogeographic significance. *Naturwissenschaften* **95**:1169–1174. doi: 10.1007/s00114-008-0438-0, PMID: 18726080

Olivier C, Houssaye A, Jalil N-E, Cubo J. 2017. First palaeohistological inference of resting metabolic rate in an extinct synapsid, *moghreberia nmachouensis* (Therapsida: anomodontia). *Biological Journal of the Linnean Society* **121**:409–419. doi: 10.1093/biolinnean/blw044

Owerkowicz T, Musinsky C, Middleton KM, Crompton AW. 2015. Respiratory turbinates and the evolution of endothermy in mammals and birds. In: *Great Transformations in Vertebrate Evolution*. Chicago: University of Chicago Press. p. 143–165.

Paaijmans KP, Heinig RL, Seliga RA, Blanford JI, Blanford S, Murdock CC, Thomas MB. 2013. Temperature variation makes ectotherms more sensitive to climate change. *Global Change Biology* **19**:2373–2380. doi: 10.1111/gcb.12240, PMID: 23630036
Ray S, Botha J, Chinsamy A. 2004. Bone histology and growth patterns of some nonmammalian therapsids. *Journal of Vertebrate Paleontology* **24**:634–648. doi: 10.1671/0272-4634(2004)024[0634:BHAGPO]2.0.CO;2

Ray S, Chinsamy A, Bandyopadhyay S. 2005. *Lystrosaurus murrayi* (THERAPSIDA, dicynodontia): BONE histology, growth and lifestyle adaptations. *Palaeontology* **48**:1169–1185. doi: 10.1111/j.1475-4983.2005.00513.x

Rey K, Amiot R, Fourel F, Rigaudier T, Abdala F, Day MO, Fernandez V, Fluteau F, France-Lanord C, Rubidge BS, Smith RM, Viglietti PA, Zipfel B, L?cuyer C, Rey K, Fourel F, Fluteau F. 2016. Global climate perturbations during the Permo-Triassic mass extinctions recorded by continental tetrapods from South Africa. *Gondwana Research* **37**:384–396. doi: 10.1016/j.gr.2015.09.008

Ruben J, Hillenius WJ, Kemp TS, Quick DE, The evolution of mammalian endothermy. 2012. Chinsamy-Turan A (Ed). *Forerunners of Mammals: Radiation, Histology, Biology*. Bloomington: Indiana University Press. p. 273–286.

Ruben JA, Jones TD. 2000. Selective factors associated with the origin of Fur and Feathers. *American Zoologist* **40**:585–596. doi: 10.1093/icb/40.4.585

Rubidge BS, Erwin DH, Ramezani J, Bowring SA, de Klerk WJ. 2013. High-precision temporal calibration of late Permian vertebrate biostratigraphy: u-pb zircon constraints from the Karoo Supergroup, South Africa. *Geology* **41**:363–366. doi: 10.1130/G33622.1

Rubidge BS. 2005. 27th Du Toit Memorial Lecture: Re-uniting lost continents - Fossil reptiles from the ancient Karoo and their wanderlust. *South African Journal of Geology* **108**:135–172. doi: 10.2113/108.1.135

Ruta M, Botha-Brink J, Mitchell SA, Benton MJ. 2013. The radiation of cynodonts and the ground plan of mammalian morphological diversity. *Proceedings of the Royal Society B: Biological Sciences* **280**:20131865. doi: 10.1098/rspb.2013.1865, PMID: 23986112

Schoch RR. 2008. The Capitosauria (Amphibia): characters, phylogeny, and stratigraphy. *Palaeodiversity* **1**:189–226.

Seebacher F. 2003. Dinosaur body temperatures: the occurrence of endothermy and ectothermy. *Paleobiology* **29**:105–122. doi: 10.1666/0094-8373(2003)029<0105:DBTTOO>2.0.CO;2

Seymour RS, Bennett-Stamper CL, Johnston SD, Carrier DR, Grigg GC. 2004. Evidence for endothermic ancestors of crocodiles at the stem of archosaur evolution. *Physiological and Biochemical Zoology* **77**:1051–1067. doi: 10.1086/422766, PMID: 15674775

Shelton CD, Sander PM, Stein K, Winkelhorst H. 2012. Long bone histology indicates sympatric species of *Dimetrodon* (Lower Permian, Sphenacodontidae). *Earth and Environmental Science Transactions of the Royal Society of Edinburgh* **103**:217–236. doi: 10.1017/S175569101300025X

Shelton CD, Sander PM. 2017. Long bone histology of *Ophiacodon* reveals the geologically earliest occurrence of fibrolamellar bone in the mammalian stem lineage. *Comptes Rendus Palevol* **16**:397–424. doi: 10.1016/j.crpv.2017.02.002

Summers AP. 2005. Evolution: warm-hearted crocs. *Nature* **434**:833–834. doi: 10.1038/434833a, PMID: 1582 9945

Sushkin PP. 1926. Notes on the pre-Jurassic tetrapoda from Russia. *Palaeontologia Hungarica* **1**:323–344.

Tattersall GJ, Leite CA, Sanders CE, Cadena V, Andrade DV, Abe AS, Milsom WK. 2016. Seasonal reproductive endothermy in tegu lizards. *Science Advances* **2**:e1500951. doi: 10.1126/sciadv.1500951, PMID: 26844295

Torsvik TH, Van der Voo R, Preeden U, Mac Niocaill C, Steinberger B, Doubrovine PV, van Hinsbergen DJJ, Domeier M, Gaina C, Tohver E, Meert JG, McCausland PJA, Cocks LRM. 2012. Phanerozoic polar wander, palaeogeography and dynamics. *Earth-Science Reviews* **114**:325–368. doi: 10.1016/j.earscirev.2012.06.007

Trotter JA, Williams IS, Nicora A, Mazza M, Rigo M. 2015. Long-term cycles of triassic climate change: a new δ18o record from conodont apatite. *Earth and Planetary Science Letters* **415**:165–174. doi: 10.1016/j.epsl.2015.01.038

Vennemann TW, Hegner E, Cliff G, Benz GW. 2001. Isotopic composition of recent shark teeth as a proxy for environmental conditions. *Geochimica Et Cosmochimica Acta* **65**:1583–1599. doi: 10.1016/S0016-7037(00) 00629-3

Viglietti PA, Smith RMH, Angielczyk KD, Kammerer CF, Fröbisch J, Rubidge BS. 2016. The *Daptocephalus* Assemblage Zone (Lopingian), South Africa: a proposed biostratigraphy based on a new compilation of stratigraphic ranges. *Journal of African Earth Sciences* **113**:153–164. doi: 10.1016/j.jafrearsci.2015.10.011

Viglietti PA, Smith RMH, Compton JS. 2013. Origin and palaeoenvironmental significance of *lystrosaurus* bonebeds in the earliest Triassic Karoo Basin, South Africa. *Palaeogeography, Palaeoclimatology, Palaeoecology* **392**:9–21. doi: 10.1016/j.palaeo.2013.08.015

Voigt S, Hminna A, Saber H, Schneider JW, Klein H. 2010. Tetrapod footprints from the uppermost level of the Permian Ikakern formation (Argana Basin, Western High Atlas, Morocco). *Journal of African Earth Sciences* **57**: 470–478. doi: 10.1016/j.jafrearsci.2009.12.003

von Grafenstein U, Erlenkeuser H, Müller J, Trimborn P, Alefs J. 1996. A 200 year mid-European air temperature record preserved in lake sediments: an extension of the δ18Op-air temperature relation into the past. *Geochimica Et Cosmochimica Acta* **60**:4025–4036. doi: 10.1016/S0016-7037(96)00230-X

Watson JM, Graves JAM. 2013. Monotreme Cell-Cycles and the evolution of Homeothermy. *Australian Journal of Zoology* **36**:573–584. doi: 10.1071/ZO9880573

Williams M, Haywood AM, Gregory FJ, Schmidt DN. 2007. *Deep-Time Perspectives on Climate Change: Marrying the Signal From Computer Models and Biological Proxies*. Geological Society of London.

Witmer LM. 1995. Homology of facial structures in extant archosaurs (birds and crocodilians), with special reference to paranasal pneumaticity and nasal conchae. *Journal of Morphology* **225**:269–327. doi: 10.1002/ jmor.1052250304
Yang Z, Besse J. 2001. New mesozoic apparent polar wander path for south China: tectonic consequences. *Journal of Geophysical Research: Solid Earth* **106**:8493–8520. doi: 10.1029/2000JB900338

Yuan W, Yang Z. 2015. The Alashan terrane did not amalgamate with North China block by the Late Permian: evidence from Carboniferous and Permian Paleomagnetic results. *Journal of Asian Earth Sciences* **104**:145–159. doi: 10.1016/j.jseaes.2014.02.010

Zazzo A, Lécuyer C, Mariotti A. 2004a. Experimentally-controlled carbon and oxygen isotope exchange between bioapatites and water under inorganic and microbially-mediated conditions. *Geochimica Et Cosmochimica Acta* **68**:1–12. doi: 10.1016/S0016-7037(03)00278-3

Zazzo A, Lécuyer C, Sheppard SMF, Grandjean P, Mariotti A. 2004b. Diagenesis and the reconstruction of paleoenvironments: a method to restore original δ18O values of carbonate and phosphate from fossil tooth enamel. *Geochimica Et Cosmochimica Acta* **68**:2245–2258. doi: 10.1016/j.gca.2003.11.009

Zhong Y-T, He B, Mundil R, Xu Y-G. 2014. CA-TIMS zircon U–Pb dating of felsic ignimbrite from the Binchuan section: Implications for the termination age of Emeishan large igneous province. *Lithos, Special Issue Permian Large Igneous Provinces: Characteristics, Mineralization and Paleo-Environment Effects* **204**:14–19.

Zhou CF, Wu S, Martin T, Luo ZX. 2013. A jurassic mammaliaform and the earliest mammalian evolutionary adaptations. *Nature* **500**:163–167. doi: 10.1038/nature12429, PMID: 23925238
# Appendix 1

## Theoretical curves construction

In order to infer therapsid thermophysiologies based on the $\delta^{18}O_p$ differences between therapsids and the associated fauna, the following theoretical framework is proposed and represented in **Figures 1, 2** and **3** as colored areas. These areas are based on the two main factors influencing animal $\delta^{18}O_p$ values: their thermometabolism (here simplified as ectotherm and endotherm) and their lifestyle (here terrestrial and semi-aquatic).

The $\delta^{18}O_p$ value ranges are estimated based on the phosphate-water temperature scale (Lécuyer et al., 2013):

$$T_b = 117.4 - 4.5(\delta^{18}O_b - \delta^{18}O_{bw}) \quad (1)$$

Where $T_b$ corresponds to body temperature, $\delta^{18}O_p$ correspond to the oxygen isotope composition of apatite phosphate, and $\delta^{18}O_{bw}$ correspond to the oxygen isotope composition of body fluids. For endothermic vertebrates, $T_b$ is assumed to be at 37°C, the average $T_b$ of most extant placental and possibly of the common ancestor of all extant mammals (Watson and Graves, 2013). For ectothermic vertebrates, $T_b$ is assumed here to reflect immediate environmental temperature ($T_e$). According to Equation 1, the $\delta^{18}O_p$ difference between an endotherm and an ectotherm can be expressed as:

$$\delta^{18}O_{p-endotherm} - \delta^{18}O_{p-ectotherm} = (T_e - 37)/4.5 + (\delta^{18}O_{bw-endotherm} - \delta^{18}O_{bw-ectotherm})$$

Because vertebrate $\delta^{18}O_{bw}$ value depend on their ecology that affects the input-output balance of body water (Luz et al., 1984; Bryant and Froelich, 1995; Kohn, 1996a), the difference $\delta^{18}O_{bw-endotherm} - \delta^{18}O_{bw-ectotherm}$ will mainly reflects that of their lifestyle (terrestrial, semi aquatic or aquatic), as well as their dependency on the surface water they ingest. Cerling et al., 2008, the $\delta^{18}O_p$ difference between water-dependent and water-independent terrestrial mammals can be up to 4‰. This is illustrated by the ranges 1 and 3 in **Appendix 1—Figure 1** and by the red and blue ranges in the main text **Figures 1A, 2A and 3A**.
<figure><img src="image_8.png" /></figure>

**Appendix 1—figure 1.** Expected latitudinal variation of the δ¹⁸Oₚ difference between vertebrate taxa of various physiologies and ecologies. Based on modern relationships between climate and phosphate-water-temperature oxygen isotope fractionation (Lécuyer et al., 2013) the following δ¹⁸Oₚ values differences are predicted: (1) corresponds to a terrestrial endotherm compared to a terrestrial ectotherm; (2) corresponds to a terrestrial endotherm compared to a semi-aquatic ectotherm; The vertical range (3) corresponds to terrestrial animals having similar thermometabolism; (4) corresponds to the difference between a terrestrial and a semi-aquatic animals having a similar thermometabolism.

DOI: 10.7554/eLife.28589.009

This variation in δ¹⁸Oₚ values can be overprinted by the animal lifestyle, a semi aquatic animal having δ¹⁸Oₚ value 3‰ lower than that of co-existing water-dependent terrestrial species, and 7‰ lower than that of co-existing water-independent terrestrial one (Appendix 1—figure 1, ranges 2 and 4, and main text Figures 1B, 2B and 3B, ranges orange and green).

The environmental temperature (T_e) used to approximate ectotherm body temperature (T_b) is estimated based on the present-day relationship between mean air temperature and latitude. Because this latitude-temperature relationship is not valid for low latitudes below about 10° (corresponding to the thermal equators of present-day Earth), we assume a constant temperature between 0° and 10° of latitude.

As a simplification, three periods are considered:
1. The Middle to Late Permian with equatorial sea surface temperatures close to the modern ones (Chen et al., 2013) (~25°C) for which we assume a present-day thermal gradient of 0.6 °C/°Latitude (Amiot et al., 2004; Williams et al., 2007) (Main text Figure 1).

2. The Early to Middle Triassic having globally warmer temperatures and a flatter temperature gradient. Based on (Trotter et al., 2015), mean equatorial sea surface temperatures were about 10°C higher than Late Permian ones (Trotter et al., 2015) (~35°C). For this time period, we assume a flatter thermal gradient that we arbitrarily set at 0.4 °C/°Latitude (Main text Figure 2).

3. The Middle Triassic to Early Jurassic having globally intermediates temperatures (Trotter et al., 2015) that we set to 30°C with an intermediate gradient of 0.5 °C/°Latitude (Main text Figure 3).

In order to take into account the possible variations of the thermal gradient within the selected periods, an interval of ±0.1 °C/°L is added for each Main text figures.

Consequently, latitudinal thermal gradients will lead to δ¹⁸Oₚ-endotherm - δ¹⁸Oₚ-ectotherm differences varying along with latitude (Appendix 1—figure 1, ranges 1 and 2). This simplified framework is used in this study to predict the following scenarios based on δ¹⁸Oₚ differences between therapsids and co-existing other tetrapods (Appendix 1—figure 1): Endothermic and terrestrial therapsid vs ectothermic and terrestrial tetrapod (range 1); Endothermic and terrestrial therapsid vs ectothermic and semi-aquatic tetrapod (range 2); Therapsids and other tetrapods having similar thermophysiologies and lifestyle (range 3); Terrestrial therapsid and semi-aquatic tetrapods having similar thermophysiologies (range 4).

## Ecology of sampled tetrapods and its impact on stable oxygen isotope compositions

Prior to interpreting differences between therapsid and associated non-therapsid tetrapod δ¹⁸Oₚ values in terms of differences in thermophysiologies, ecological traits must be considered as they also affect the oxygen isotope compositions of apatite phosphate. Indeed, δ¹⁸Oₚ values recorded in phosphatic tissues depend on the animal body temperature, as well as on the oxygen isotope composition of body water, the latter being affected by water turnover rate and isotopic fractionations associated with water loss (Bryant and Froelich, 1995). Indeed, water loss through transcutaneous evaporation, sweat and exhaled water vapour tends to ¹⁸O-enrich the remaining body water (Kohn, 1996a). This isotopic enrichment is amplified or mitigated depending on the rate of water turnover, which depends itself on the animal ecology. Aquatic and semi-aquatic animals regularly ingest and release large amounts of their environmental water compared to terrestrial ones, which in turn reduces the magnitude of body water ¹⁸O-enrichment relative to that of environmental water (Kohn, 1996a).

## Stereospondyl amphibians

The sampled stereospondyl amphibian clade (Almasaurus habbazi, Lydekkerina, Microposaurus, Paracyclotosaurus, Rhinesuchus and Xenotosuchus) includes semi-aquatic to aquatic animals (Schoch, 2008). A study dedicated to Rhinesuchus palaeohistology concluded that it had a fully aquatic lifestyle (McHugh, 2014), whereas Lydekkerina, a basal stereospondyl, was amphibious with a tendency to be terrestrial (Canoville and Chinsamy, 2015).

## Pareiasaurid and bolosaurid parareptiles

The sampled parareptiles include the Chinese Bolosauridae, Belebey chengi, and the South African Pareiasauria, Pareiasaurus and basal pareiasaurs (Embrithosaurus, Nochelosaurus or
Bradyssaurus; (Lee, 1997). The Chinese taxon is considered to have been terrestrial (Berman et al., 2000; Müller et al., 2008), but interpretations of the life habits and habitats of the South Africa pareiasaurs still lack a consensus. According to authors, they are variously considered fully aquatic (Ivakhnenko, 2001), semi-aquatic (Kriloff et al., 2008) or fully terrestrial (Voigt et al., 2010; Canoville et al., 2014). The low δ¹⁸Oₚ values measured in both *Pareiasaurus* and *Rhinesuchus* are similar to each other, *Pareiasaurus* having slightly higher δ¹⁸Oₚ values than *Rhinesuchus* (~1‰), but significantly lower than those measured in therapsids (about 4‰ to 8‰). Considering *Rhinesuchus* as an aquatic ectothermic amphibian (McHugh, 2014), it is suggested that the ectothermic *Pareiasaurus* was also aquatic or semi-aquatic as previously suggested (Ivakhnenko, 2001; Kriloff et al., 2008). In contrast to the semi-aquatic *Pareiasaurus* from the Lower Daptocephalus AZ, the *Tapinocephalus* AZ pareiasaurs may have been terrestrial as recently proposed in a study on pareiasaur ecology and based on oxygen and carbon isotope compositions of apatite carbonate (Canoville et al., 2014). Moreover, those authors found no significant differences between therocephalians and pareiasaurs, agreeing with our dataset, and a ~ 3‰ overlap between pareiasaur and dinocephalian δ¹⁸O_c values, even though most dinocephalian values are several per mil lower than those of pareiasaurs.

## Archosauriformes

Archosauriformes are represented by the Erythrosuchidae *Erythrosuchus* and *Shansisuchus shansisuchus*, by the Proterosuchidae *'Chasmatosaurus'* *yuani*, as well as by an indeterminate basal sauropodomorph dinosaur. Proterosuchid lifestyle is still unresolved, but *Botha-Brink and Smith, 2011* favoured a terrestrial rather than an aquatic lifestyle, whereas the erythrosuchids are considered as the largest terrestrial predators of their time (*Botha-Brink and Smith, 2011*).

## Lystrosaurid therapsids

Therapsids are generally considered as terrestrial dwellers (Kemp, 2012) with a few uncertainties such as *Anteosaurus* considered as riparian (Boonstra, 1955, Boonstra, 1962) and *Lystrosaurus* considered by some as semi-aquatic (Germain and Laurin, 2005; Ray et al., 2005; Botha-Brink and Angielczyk, 2010; Canoville and Laurin, 2010).

Because the *Lystrosaurus-Lydekkerina* δ¹⁸Oₚ difference falls within the same range as the *Kannemeyeria-Xenotosuchus* and *Kannemeyeria-Microposaurus* values from the *Cynognathus* subzone B assemblage, *Lystrosaurus* could be considered as having the same thermophysiology and lifestyle as *Kannemeyeria*, i.e. a terrestrial endotherm (Figure 2A). Because *Lystrosaurus* may have either been terrestrial (*Botha and Smith, 2006*) or semi-aquatic (*Canoville and Laurin, 2010*), it cannot be totally excluded that *Lystrosaurus* shared a similar lifestyle and thermometabolism with the amphibian *Lydekkerina*. If *Lystrosaurus* was indeed semi aquatic, then it could also be predicted to have been an ectotherm.

Given that the thermophysiology of the terrestrial archosauriform *'Chasmatosaurus'* remains unclear, the two similar δ¹⁸Oₚ values observed in the lower Triassic Jiucaiyuan Fm. of Xinjiang province of China may indicate either that both the lystrosaurid and the proterosuchid shared similar thermometabolism and lifestyle (terrestrial endotherms or terrestrial ectotherms; Figure 2), or that one of them may have had a lower body temperature. Indeed, compared to low latitudes where ectotherms have slightly lower δ¹⁸Oₚ values than co-existing endotherms due to differences in metabolic activity and body temperatures, mid latitude ectotherms have even lower body temperatures (reflecting the environmental thermal gradient) resulting in higher δ¹⁸Oₚ values that mimic those of co-existing endotherms (*Amiot et al., 2004*). It is worth noting that based on histological features, *Proterosuchus*, had intermediate growth rates, suggesting an intermediate
thermometabolism (Botha-Brink and Smith, 2011) that could apply to the Chinese 'Chasmatosaurus'. Considering a terrestrial lifestyle for the Chinese lystrosaurids (see above), and the conflicting hypotheses of Lystrosaurus from South Africa (terrestrial endotherm vs. semi-aquatic ectotherm), the hypothesis for terrestrial endothermy for the lystrosaurids is favoured.