The profile is currenly unclaimed by the seller. All information is provided by CB Insights.

## Stage

Leveraged Buyout | Alive

Provider of cementum products.

Netherlands

## Predict your next investment

The CB Insights tech market intelligence platform analyzes millions of data points on venture capital, startups, patents , partnerships and news mentions to help you see tomorrow's opportunities, today.

## Latest Cementum News

Reptile-like physiology in Early Jurassic stem-mammals

Oct 12, 2020

Abstract Despite considerable advances in knowledge of the anatomy, ecology and evolution of early mammals, far less is known about their physiology. Evidence is contradictory concerning the timing and fossil groups in which mammalian endothermy arose. To determine the state of metabolic evolution in two of the earliest stem-mammals, the Early Jurassic Morganucodon and Kuehneotherium, we use separate proxies for basal and maximum metabolic rate. Here we report, using synchrotron X-ray tomographic imaging of incremental tooth cementum, that they had maximum lifespans considerably longer than comparably sized living mammals, but similar to those of reptiles, and so they likely had reptilian-level basal metabolic rates. Measurements of femoral nutrient foramina show Morganucodon had blood flow rates intermediate between living mammals and reptiles, suggesting maximum metabolic rates increased evolutionarily before basal metabolic rates. Stem mammals lacked the elevated endothermic metabolism of living mammals, highlighting the mosaic nature of mammalian physiological evolution. Introduction Recent discoveries and analyses have revolutionized our knowledge of Mesozoic mammals, revealing novel aspects of their ecology 1 , 2 , development 2 , 3 , systematics 2 , 3 and macroevolution 4 , 5 . However, details of physiology are more difficult to determine from fossils, and our knowledge of physiological evolution remains comparatively poor. Living mammals are endotherms, possessing the ability to control and maintain metabolically produced heat and have a substantially higher capacity for sustained aerobic activity than ectothermic animals 6 , 7 , 8 . The origin of endothermy is an important event in mammalian evolution, often noted as key to their success 6 , 7 , 8 . There are a number of competing evolutionary hypotheses for the origin of endothermy: (a) selection for higher maximum metabolic rates (MMRs) enhanced sustained aerobic activity 6 , 9 , 10 , (b) selection for higher basal metabolic rates (BMRs) enhanced thermoregulatory control 11 , 12 , or (c) MMRs and BMRs evolved in lockstep with each other 7 , 8 . Direct evidence from living mammals to support these hypotheses is equivocal 7 . Recent analyses find no long-term evolutionary trend in BMR 13 contradicting earlier suggestions of increasing BMR throughout the Cenozoic 12 and so implying that the Middle Jurassic (∼170 Ma) most recent common ancestor (MRCA) of living mammals 13 possessed a BMR within the range of present-day mammals. Several indirect indicators of metabolic physiology in fossil synapsids have been suggested but provide contradictory evidence for the timing of the origin of endothermy and its evolutionary tempo. These include: the presence of fibrolamellar long-bone histology 14 , 15 , first seen in the Early Permian synapsid Ophiacodon ∼300 million years (Ma) ago 16 ; the presence of an infraorbital canal and lack of parietal foramen, used to infer facial whiskers, fur, lactation and endothermy in Early Triassic (∼245 Ma) cynodonts 17 ; inferred maxillary nasal turbinates in the Late Permian (∼255 Ma) therapsid Glanosuchus, used to suggest that mammalian levels of endothermy evolved by the Late Triassic (∼210 Ma) 18 ; a trend towards increased relative brain size initiated in Late Triassic non-mammaliaform cynodonts 19 and the mammaliaform (stem mammal sensu Rowe 20 ) Morganucodon 21 , 22 ; and acquisition of a parasagittal gait in the Early Cretaceous (∼125 Ma) therian mammals Eomaia and Sinodelphys 23 . Several recent studies provide more quantitative links to physiological parameters. Oxygen isotopes were used to infer elevated thermometabolism in Middle–Late Permian (∼270–255 Ma) eucynodonts 24 , red blood cell size diminution in Late Permian (∼255 Ma) eutheriodontid therapsids was linked via two proxies to increased MMR 25 and osteocyte lacuna shape correlations suggested “mammalian” resting metabolic rates in Permo-Triassic (∼250 Ma) dicynodonts 26 . However, the inconsistency of these characters, in time and with respect to phylogeny 27 , 28 , along with re-assessments of function in relation to endothermy 7 , 29 , 30 , limit their use as conclusive indicators of modern mammalian levels of endothermy in fossil taxa. Such temporal and phylogenetic heterogeneity suggests that the evolution of mammalian endothermy followed a complex, mosaic pattern with different physiological aspects likely evolving independently, and at separate rates, towards current mammalian levels. Additionally, few of these physiological proxies are directly related to measurable aspects of metabolic rate. To address these issues, we use two proxies to improve understanding of physiology at one of the most important nodes along this transition. We do so by estimating BMR and growth rate, and calculating a known proxy for MMR, for two of the earliest known mammaliaforms, Morganucodon and Kuehneotherium 1 , 31 . Using cementochronology to estimate maximum lifespan by counting growth increments in synchrotron radiation-based micro-computed tomographic (µCT) data of fossil dental cementum, we estimate that both taxa had significantly longer lifespans than extant mammals of comparable size. By regressing lifespan against BMR and postnatal growth rate in extant mammals and reptiles, we in turn estimate significantly lower values for both of these metrics for the earliest mammaliaforms. However, when we compare the blood flow index (the ratio between femoral nutrient foramina area and femur length that serves as a proxy for MMR) of Morganucodon with those of extant taxa, we find that Morganucodon had an intermediate value between living mammals and reptiles. These results suggest that basal mammaliaforms occupied a metabolic grade similar to living reptiles and had yet to achieve the endothermic physiology of living mammals. Results Lifespan as a proxy for mammaliaform physiology We used maximum lifespan (i.e. the single longest known lifespan of a taxon) estimates for fossil mammaliaform taxa as a proxy for both BMR 32 and postnatal growth rate 33 . In extant tetrapods, negative correlations exist between maximum lifespan and BMR 32 and between maximum lifespan and growth rate 33 , 34 . In general, the longer a mammal’s lifespan, the lower its size-adjusted BMR and growth rate. Growth rates have been shown to correlate strongly with metabolic power in extant vertebrates, with endotherms growing an order of magnitude faster than ectotherms 34 , 35 . Maximum lifespan is an applicable value for fossil samples, as, unlike other metrics (e.g. 10% most long lived or mean lifespan of a cohort), it does not rely on cohort- or population-based statistics that fossil samples cannot fulfil 36 . This value is also less susceptible to extrinsic population-level factors on lifespan, such as disease or predation, and relates most closely to the physiological limit of lifespan of an organism. An accurate assessment of maximum lifespan in fossil mammals can therefore be used to estimate their metabolic potential. To estimate mammaliaform lifespans, we used cementochronology. This well-established technique, which counts annual growth increments in tooth-root cementum, has been used to record lifespans in extant mammals 37 with >70 species aged using this technique 38 . Cementum is a mineralized dental tissue surrounding the tooth root (Fig. 1a, b ), attaching it to the periodontal ligament and anchoring the tooth within the alveolus 37 . Growth of cementum is continuous throughout life in extant mammals and seasonally appositional in nature, forming a series of increments of differing thickness and opacity when viewed in histological thin sections under light microscopy. The correlation between increment count and chronological age is well documented, with one thick and one thin increment deposited every year 38 . It has been shown that the thin, hyper-mineralized opaque increments record growth rate reduction in less favourable seasons 39 , 40 . Fig. 1: Cementum of Morganucodon and Kuehneotherium. a, b Three-dimensional reconstructions of a Morganucodon right lower molar tooth NHMUK PV M 104134 (voxel size 2 μm, μCT) and b Kuehneotherium right lower molar tooth NHMUK PV M 21095 (voxel size 1.2 μm, PPC-SRµCT). Green = cementum. c, d Transverse PPC-SRµCT virtual thin sections (0.33 μm voxel size) of roots of c NHMUK PV M 104134 and d NHMUK M 27436. Red bracketed line highlights extent of cementum surrounding dentine. e, f Close-ups of boxes in c, d, with five and four circumferential light/dark increment pairs highlighted by red arrows, respectively. g Synchrotron nanotomographic virtual thin section of NHMUK PV M 104134 (30 nm voxel size) provides close-up of region close to box in e. Vertical red arrows = cementum increments; horizontal blue arrows, dashed blue lines and Sf = radial bands of Sharpey’s fibres; yellow dashed line and glT = granular layer of Tomes; green dashed line and hlH-S = hyaline layer of Hopewell-Smith. Scale bars represent 500 µm in a, b, 100 µm in c, d, 30 µm in e, f and 10 µm in g. We estimated a Qi of 3.829e−7 mm3 for Morganucodon and compared this with published and new data (“Methods”) for extant mammals (n = 69) and reptiles (n = 30). The latter includes varanids (n = 8), which in the absence of mammalian predators fill an active hunting niche and tend to have mammalian MMR levels while retaining reptilian BMR levels 51 (Supplementary Table  2 ). One-way ANCOVA comparisons show that means of GLS regression slopes for extant mammals and non-varanid reptiles are significantly different (p < 0.01) while the slopes are similar (p = 0.16). Log10 GLS regression of body mass and Qi shows that Morganucodon is further above (higher Qi for its mass) the non-varanid reptile mean than all non-varanid reptiles (phylogenetically informed statistical comparisons were not used here due to the non-significant lambda values showing no phylogenetic signal in the taxa used, see “Methods”). However, Morganucodon is also slightly further from the mammalian mean than the non-varanid reptile mean and considerably closer to small non-varanid reptile species data points than those of small mammalian species (Fig. 7b ). This intermediate Qi, and so inferred intermediate MMR, suggests that, while retaining typical reptilian BMR and growth rates, Morganucodon had MMR above non-varanid reptiles but not as high as mammals or actively foraging varanid reptiles. Discussion We have used two quantitative proxies to determine the metabolic status of early mammaliaforms. Relatively long lifespans for both Morganucodon and Kuehneotherium result in SMR/BMR and growth rate estimates equivalent to modern reptiles of comparable size and indeed at the higher lifespan/lower BMR/slower growth end of the reptile scale for Morganucodon. This is true whether we compare our fossil estimates to wild lifespans of extant taxa or estimate fossil “captive” lifespans and compare them to captive values for extant taxa. In contrast, femoral blood flow estimates (Qi) suggest that the MMR of Morganucodon was intermediate between extant non-varanid reptiles and mammals. We therefore infer that in Morganucodon increased MMR (and so also absolute aerobic capacity (AAC) = MMR – BMR) was initially selected for before BMR and that the MMR-first hypothesis 6 is the best-supported model for the evolution of mammalian endothermy. We suggest that at least Morganucodon, if not also Kuehneotherium, occupied a metabolic grade approaching extant varanids: able to undergo longer bouts of aerobically demanding activity than non-varanid reptiles but not capable of sustaining either mammalian levels of aerobic activity or the elevated thermometabolism exhibited by living endotherms. Evidence from non-mammalian synapsids (including changes in gait 7 , long bone histology 15 and development of secondary osteological features correlated with increased metabolic rate 17 , 18 ) indicate unquestionable changes in physiology from pelycosaur- to mammaliaform-grade taxa. Determinate growth 47 and reduction of dental replacement (diphyodonty) in basal mammaliaforms permitted more precise occlusion 52 , which has been considered a key innovation in the development of mammalian endothermy by enabling increased assimilation and higher metabolism 53 . However, determinate growth and diphyodonty appear to have preceded the appearance of modern mammalian levels of endothermy, at least in Morganucodon and Kuehneotherium. We therefore suggest that the development of precise occlusion in basal mammaliaforms 52 may be more associated with dietary specialization and niche partitioning 1 . Comparison of our results to those of other recent studies of physiology in fossil synapsids supports the hypothesis of a complex, mosaic pattern for the evolution of endothermy, with different characters being selected for at different rates through time, and with respect to phylogeny. For example, the size diminution associated with the cynodont–mammaliaform transition 54 may have reversed the evolutionary trajectory of some previous histological proxies for endothermy 55 , contributing to the complex, contradictory patterns observed. Our study also suggests that more work is needed to compare fossil and extant ectothermic and endothermic taxa directly in order to better understand their relative metabolic properties. Many previous studies rely on simple binary divisions, such as the presence/absence of fibrolamellar bone and/or respiratory nasal turbinates. These proxies cannot represent accurately the complex series of physiological characteristics that range between “ectothermy” and “endothermy” and are frequently distributed homoplastically across the synapsid phylogeny, individually and with respect to each other. Other studies provide relative data such as preserved apatite oxygen isotopes 24 that allow comparisons with co-habiting ectothermic taxa but cannot be directly compared to extant data and so do not suggest where the studied fossil taxa fall in the metabolic spectrum of extant vertebrates. However, our results are compatible with recent work on living mammals, suggesting that the BMR of the Middle Jurassic (∼170 Ma) mammalian MRCA was comparable to present-day values 13 . This indicates that evolution towards modern-day mammalian endothermy occurred during the 25 million year-long Early Jurassic and suggests that the mammalian mid-Jurassic adaptive radiation 4 , 5 was driven by this or vice versa. In conclusion, our data offer a direct link to measurable aspects of endothermy, such as BMR and MMR, at a key point in mammalian evolution. Further work applying these methods to additional Mesozoic mammaliaforms and mammals, and comparison with evidence from other physiological characteristics, will allow the evolutionary tempo and mode of multiple aspects of mammalian physiology to be determined. The early mammaliaforms Morganucodon and Kuehneotherium possessed surprisingly low, reptile-like metabolic rates plus a mixture of plesiomorphic and derived characters 7 relating to life history and physiology. Ultimately, we can no longer assume that the endothermic metabolism of living mammals had evolved in the earliest mammaliaforms. Methods Choice of fossil specimens All Morganucodon and Kuehneotherium specimens are from the Early Jurassic St Brides Island fissure suite, from Glamorgan, South Wales (UK) 1 , 56 , 57 . Morganucodon specimens used for the cementum analysis are from the collections of the Natural History Museum, London, UK (NHMUK), Kuehneotherium specimens are from the collections of the NHMUK and University Museum of Zoology Cambridge, UK (UMZC) and the Morganucodon femora specimens are from the collections of UMZC. The Morganucodon specimens are all one species, M. watsoni, from Pontalun 3 fissure, excavated in 1962–1963 in Pontalun quarry (now known as Lithalun) 45 , 56 . Kuehneotherium is not abundant in Pontalun 3 fissure and the specimens used here for cementum analysis are from three Glamorgan fissures in Pontalun and Pant quarries: Pontalun 3, Pant 2, and Pant 4. Pontalun 3 and Pant 2 fissures have a relatively impoverished fauna, and Pant 4 has a more diverse biota 45 (Supplementary Data  2 ). Kuehneotherium from Pontalun 3 fissure is Kuehneotherium praecursoris 57 , but those from the Pant fissures are considered to be different species, based on small molar differences 58 . All specimens from Pontalun 3 fissure in the NHMUK (Morganucodon and Kuehneotherium) were prepared by immersing dried blocks of clay matrix in hot tap water, with the addition of dilute hydrogen peroxide or sodium hexametaphosphate (Calgon) if required, but those in the UMZC collection were prepared with 10% acetic acid. Kuehneotherium specimens from Pant quarry (Pant 2 and Pant 4) were collected by the Kermack team from University College London, between 1955 and 1978, and are from a harder matrix, which was prepared with 15% acetic acid 58 . The Morganucodon and Kuehneotherium specimens used for cementum analysis are either isolated teeth or dentary specimens with a range of teeth in situ. The lower dental formula of Morganucodon is i4:c1:p4:m4–5 and for Kuehneotherium it is i? :c1:p6:m6 1 . In Morganucodon, the majority of isolated teeth measured were lower second molars (m2) since they are easily identified, are relatively large, have robust separated roots and, as anterior diphyodont molars, they should erupt relatively early and are not replaced, therefore offering a near-complete record of life history. However, i4, c, p3, p4, m1, m3 and m4 teeth were also studied in dentulous Morganucodon specimens. Dentulous Kuehneotherium specimens are extremely rare and so isolated teeth were scanned in almost all cases. For Kuehneotherium, all appropriate teeth were chosen; the distinctive p5 and p6, and a range of lower and upper molars, which can be identified to anterior, mid or posterior tooth row on the degree of triangulation 58 . Tomographic imaging of cementum Pilot scans of two Morganucodon lower second molars (NHMUK PV M 104131 and NHMUK PV M 104132) were carried out in 2011 on the nanotomography imaging beamline ID22 at the European Synchrotron Radiation Facility (ESRF), Grenoble, France (project EC 1064). For ID22, we used the following experimental settings for computed tomographic (CT) imaging: X-ray energy of 29.6 keV, 1999 projections over a 180° rotation, 0.5 ms exposure time, 321 nm voxel size, 405 mm sample-to-detector distance, diamond window, a 20-µm-thick LSO scintillator doped with Tb, and 1.5 mm Al filter, single propagation distance tomography. During a 4-day experiment at the ID19 beamline of the ESRF (18/04/2014-22/04/2014), 71 additional Morganucodon specimens (52 isolated teeth and 19 dentaries) and 2 pilot Kuehneotherium specimens were scanned (project ES 152). A single harmonic U13 undulator was used as the X-ray source, delivering a pink X-ray beam with peak energy at 26.5 keV, with a 1.4 mm Al filter used to cut the background lower energies. The detector was a microscope optic system coupled to a sCMOS sensor (PCO edge 5.5), mounted with a 10-µm thick GGG:Eu scintillator. Scans were performed using single propagation distance tomography (15 mm sample-to-detector propagation distance), an exposure time of 250–300 ms, 2499 angular projections over a 360° scan and at voxel sizes of 280, 347 and 700 nm. Subsequently, 117 additional Kuehneotherium specimens (116 isolated teeth and one dentary) and 12 additional Morganucodon dentary specimens were scanned during a 3-day experiment at the TOMCAT tomographic beamline of the Swiss light Source (SLS), Villigen, Switzerland (13/04/2015-16/04/2015). The beam was set at an energy of 20 keV using a double multilayer monochromator, a LSO:Tb scintillator and a pco.EDGE 5.5 detector. Samples were scanned using single propagation distance tomography (14 mm sample-to-detector propagation distance), an exposure time of 150 ms and 1500 angular projections over a 180° scan at a voxel size of 330 nm. A Kuehneotherium lower molar (UMZC Sy 141) was imaged at 1.2 μm voxel size (with an exposure time of 150 ms and 1500 angular projections over 180°) to provide the 3D volume presented in Fig. 1b . Three juvenile Morganucodon dentary specimens, with roots from final deciduous premolars (NHMUK PV M 27312, NHMUK PV M 27474 and NHMUK PV M 27475), and an older individual with extensive molar wear (NHMUK PV M 27465), were scanned during a 3-day experiment at the TOMCAT beamline of the SLS (07/03/2016-10/03/2016). The beam energy was set at 21 keV using a double multi-layer monochromator. Samples were scanned using single propagation distance tomography (14 mm sample-to-detector propagation distance), an exposure time of 200 ms and 1601 angular projections over a 180° scan at a voxel size of 330 nm. CT reconstructions of the above tomographic data were generated using a filtered back-projection algorithm coupled with “Paganin-style” single distance phase retrieval 59 algorithms developed in-house at the respective beamlines 60 , 61 . For data from ID19, β = 8.1 × 10−8, δ = 9.8 × 10−9. For data from TOMCAT, β = 3.7 × 10−8, δ = 1.7 × 10−10. Of the 71 Morganucodon molar specimens imaged at beamline ID19 in 2014, 4 were additionally imaged at the nano-imaging beamline ID16A of the ESRF synchrotron (project ES 152). These were imaged using holotomography 62 from four propagation distances. Holograms were recorded using a charge-coupled device detector with an effective pixel size of 3 µm and a 23-µm-thick GGG:Eu scintillator at both 17 keV and 33.6 keV. The selected voxel sizes were 10, 25, 30 and 130 nm. The number of angular projections recorded over 180° varied between 1200 and 2000 and the exposure times were set at 250–800 ms. To generate the image containing the virtual section in Fig. 1g with 30 nm voxel size, the four focus-to-sample distances were 2.65, 13.19, 15.36 and 19.87 mm and the sample-to-detector distance was 1.2648 m. To produce the 3D model presented in Fig. 1a , a Morganucodon lower molar (NHMUK PV M 104134) was imaged using μCT at the University of Helsinki in March 2013. μCT was performed using a Nanotom 180 NF (phoenix X-ray Systems & Services GmbH) with a CMOS detector (Hamamatsu Photonics) and a high-power transmission-type X-ray nanofocus source with a tungsten anode. A total of 900 angular projections were collected for a 180° rotation, at an exposure time of 1840 ms and a voxel size of 2 × 2 × 2 μm3. The raw projection data were reconstructed using filtered back-projection by the reconstruction software datos|x rec, supplied by the system manufacturer. Increment counting and creation of virtual thin sections Cementum increments were counted in CT data using modifications to the techniques suggested by the Cementochronology Research Program 38 to take into account the 3D nature of the PPC-SRµCT cementum data. First, the cementum was visually inspected throughout the entire volume of each scan, in transverse PPC-SRµCT slices using ImageJ/Fiji 63 to distinguish between specimens that could be confidently interpreted as preserving cementum increments or those that were too badly affected by diagenesis for increment counting. Phase-contrast imaging of incremental features is understood to be prone to recurrent destructive interference patterns from Fresnel diffraction that create periodic blurring at differing frequencies when they are scanned using inappropriate experimental parameters (principally X-ray energy, sample-to-detector propagation distance and voxel size). However, our parameters produce blurring frequencies that are too narrow (approximately 500–900 nm) to significantly affect the contrast between cementum increments (1–3 μm radial thickness) (Tafforeau, personal observation). Therefore, no significant masking of increments from Fresnel diffraction blurring should be expected in our data. In specimens that preserved increments, volumes were inspected by eye to identify regions of highest increment contrast with no lensing and/or coalescence between increments. Increments identified in these regions were followed by eye throughout the entire cementum tissue surrounding these regions, both longitudinally and transversely through the root, in order to distinguish between principal increments and accessory increments formed by lensing and coalescence of primary increments in discrete portions of the tissue (Supplementary Fig. 1 ). Primary increments were distinguished as those that persisted vertically through the entire scanned region of cementum, whereas accessory increments lasted only for short periods before coalescing into neighbouring increments (Supplementary Fig. 1 ). Once regions of highly contrasting primary cementum increments had been identified, virtual thin sections of these regions were created. This was performed by isolating ten transverse PPC-SRµCT slices through each region and summing their greyscale values using the “Sum slices” option of the “Z projection” tool in ImageJ/Fiji to create a new image of increased contrast between dark and light cementum increments and reduced image noise. Between three and five virtual thin sections were created for all specimens with readable cementum increments. For each virtual thin section, increments were counted manually by three different observers: Observer One (E.N.) had considerable experience in counting cementum increments (>100 specimens studied); Observer Two (K.W.) had training in counting cementum increments (30 specimens studied under guidance from Observer One) and experience in studying growth patterns in PPC-SRµCT data of long-bones; Observer Three (C.N.) had no prior experience in counting increments or studying growth patterns. Each observer studied virtual thin sections blind, after collections of virtual thin sections were numbered and randomized between specimens using the RAND function in Microsoft Excel. For each observer, the final increment count for every specimen was determined as the maximum number counted in all of its virtual thin sections. Once each observer had counted increments in every virtual thin section (Supplementary Fig. 3 and (Supplementary Data  1 ) and the precision between their increment counts was compared by calculating the CV (Eq. 1 ; Supplementary Table  1 ): $${\mathrm{CV}} = \left[ {\frac{{{\mathrm{standard}}\,{\mathrm{deviation}}}}{{{\mathrm{mean}}}}} \right]^\ast 100$$ (1) 3D modelling of cementum increments 3D modelling was performed, using the Avizo image analysis software (Avizo 8.0; Thermo Fisher Scientific), on a subsample of teeth comprising the first tooth imaged using PPC-SRµCT (NHMUK PV M 104131) and the specimens that provided the highest cementum increment counts for each fossil taxon (NHMUK M 104127 for Morganucodon and UMCZ Sy 141 for Kuehneotherium). Original tomographic data were downsampled in each axis by a factor of two to decrease manual processing time while retaining sufficient spatial resolution to preserve the cementum increments. Principal increments originally identified by eye were manually traced in each PPC-SRµCT slice and assigned to different materials in the “label field” tool-kit. Models were subsequently created and analysed using the “surface view” feature of Avizo. This allowed the pattern of incrementation to be viewed with a 3D perspective, to test whether increments defined in two-dimensional slices were true principal increments or diffuse accessory increments (Supplementary Fig. 1 ). Fossil body mass estimation and choice of extant taxa Body mass for our fossil taxa was estimated using two techniques, based on scaling between single cranial dimensions and measured body mass in extant mammals. Maximum body mass estimates were made using the scaling relationship (Eq. 2 ) between dentary length (mm) and body mass (g) published for extant marsupial mammals of small body mass and subsequently used to estimate body mass in several Mesozoic fossil taxa from the Late Jurassic Morrison Formation and elsewhere 48 : $${\mathrm{ln}}\,{\mathrm{body}}\,{\mathrm{mass}}\left( {\mathrm{g}} \right) = 2.9677^\ast \left( {{\mathrm{ln}}\,{\mathrm{dentary}}\,{\mathrm{length}}\left( {{\mathrm{mm}}} \right)} \right)-5.6712$$ (2) Dentary lengths are from published CT reconstructions 1 (20 mm for Morganucodon and 21.9 mm for Kuehneotherium), resulting in body mass estimates of 25.0 g and 32.7 g respectively. However, this estimate may be an overestimate for Kuehneotherium as it has a longer, more gracile, dentary relative to Morganucodon, due to the differing feeding ecologies of the two taxa 1 . Second, minimum body mass estimates were calculated using the scaling relationship (Eq. 3 ) between skull length (mm) and body mass (g) found for 64 extant species of small “lipotyphlan” insectivores 64 and used to estimate the body masses of several Mesozoic mammaliaforms 21 , 64 , 65 . $${\mathrm{ln}}\left( {{\mathrm{body}}\,{\mathrm{mass}}\left( {\mathrm{g}} \right)} \right) = 3.68 \ast \left( {{\mathrm{ln}}\left( {{\mathrm{skull}}\,{\mathrm{length}}\left( {{\mathrm{mm}}} \right)} \right)} \right.-3.83$$ (3) Due to a lack of complete, diagnostic cranial material for UK samples of Morganucodon and Kuehneotherium, their skull lengths were estimated using a scaling relationship between dentary length and skull length calculated for Morganucodon oehleri from Rowe et al. Figure S1 21 (skull length (mm) = 1.0458 × jaw length (mm)). Following this relationship, we used dentary length 1 , 21 to estimate a skull length of 21.0 mm for M. watsoni and 22.9 mm for Kuehneotherium, which results in body mass estimates of 10.7 and 14.9 g, respectively (as above, this may be an overestimate for Kuehneotherium). Body masses presented in Figs. 5 and 7 are mean values of the two estimates (17.9 g for Morganucodon and 23.8 g for Kuehneotherium). Extant data Information on maximum wild lifespan and mean body mass was obtained for 278 extant terrestrial mammal species (body mass: mean = 70.2 kg, standard deviation = 395 kg; wild lifespan: mean = 10.8 years, standard deviation = 11.2 years) and 256 extant terrestrial non-avian reptile species (body mass: mean = 3.6 kg, standard deviation = 20 kg; wild lifespan: mean = 12 years, standard deviation = 14 years) (Supplementary Data File  3 ). Information on maximum captive lifespan was obtained for 644 extant terrestrial mammal species (body mass: mean = 62.2 kg, standard deviation = 324 kg; captive lifespan: mean = 15.5 years, standard deviation = 11.5 years) and 866 extant reptile species (body mass: mean = 2.4 kg, standard deviation = 15.9 kg; captive lifespan: mean = 15.3 years, standard deviation = 11.6 years). Information for flying or gliding taxa (including birds) was not included as their body masses are known to be secondarily reduced to aid their lifestyle and so may distort observed overall relationships between lifespan and body mass. In particular, birds and bats live on average three-to-four times as long as terrestrial mammals of similar body mass 66 . Similarly, marine mammals were not included as their environment positively affects body mass, allowing significantly higher body masses than terrestrial taxa 67 . Mammalian wild and captive lifespan data were obtained from the primary literature. The majority of these data originated from a download from the online databases of the Max Planck Institute ( https://www.demogr.mpg.de/longevityrecords/0203.htm ) on 10/03/2019. The majority of reptile wild and captive lifespans were obtained from the supplementary information of Scharf et al. 68 . Body mass estimates were obtained from an online Ecological Archives database 69 ( http://www.esapubs.org/archive/ecol/E084/094/metadata.htm ), the AnAge database ( https://genomics.senescence.info/species/ ) and published data 50 , 68 . Although the Max Planck Institute data are entitled “longevity records”, it notes in the introduction that this is equivalent to “highest documented age”. The authors of the database believe that, despite any analytical and conceptual shortcomings, the records are useful in a number of comparative and disciplinary contexts, including demographic, gerontological, ecological and evolutionary. Similarly, the AnAge database states that maximum longevity is synonymous with maximum lifespan. Data on SMR for extant mammals (n = 117; mean = 915; standard deviation = 3018) and reptiles (n = 55; mean = 62.3; standard deviation = 159) was obtained from a published electronic appendix 50 . SMRs for mammals were included in this data set only if measured under basal metabolic conditions 50 and so are synonymous with BMRs. Data on the postnatal growth rate constant K was obtained upon request for extant mammals 32 (n = 115, mean = 0.03, standard deviation = 0.03) and reptiles (n = 33, mean = 0.002, standard deviation = 0.003) from previous publications 70 , 71 . The growth rate constant K is a magnitude free measure of growth rate, measured in days−1, that is considered the best measure for comparing growth rates between multiple species 72 . Physiological metrics (lifespan/body mass/SMR/K) were compared using PGLS regression in the “R” statistical environment with the “ape”, “geiger”, “nlme” and “phytools” packages installed. For each regression, phylogenetic subsets of squamate and mammal taxa were downloaded from https://vertlife.org (100 trees per subset, birth-death node-dated completed tree distribution for mammals) representing the phylogenetic relationships of every taxon in the respective subsample. For crocodilian taxa, a phylogeny was manually constructed in R following the time-calibrated phylogeny of Oaks 73 (using node mean age values from the species tree/90 My maximum analysis) and added to the base of the squamate phylogeny to produce a reptile clade. Each subset was investigated to find the tree that produced the highest PGLS F values and this tree was then input into the “corPagel” covariance structure for the “gls()” function to produce a phylogenetically informed regression model between the respective metrics. The relationships between body mass and maximum lifespan were compared between extant mammals and extant reptiles using phylogenetic ANCOVA following Smaers et al. 74 using the “gls.ancova” function in “R” with the same libraries installed as for PGLS. The individual phylogenetic subsets representing the mammal and reptile samples were combined with connecting branches connecting the Last Common Ancestors (LCAs) of the crocodilian, squamate and mammal clades. The LCA node connecting squamates and crocodilians was dated at 275.9 Ma using the mean of minimum and maximum node ages from Benton et al. 75 ; similarly, the LCA node connecting reptiles (squamates and crocodiles) and mammals was dated at 325.5 Ma 75 . The mammal and reptile clades were then distinguished using two factor elements, and the influence of these factors upon the fit of the data to a series of models of increasing complexity was tested. First, a baseline model was produced consisting of one slope and one intercept for the entire data set. Three models were then compared to this baseline; the first varied in slope between mammals/reptiles, keeping their intercepts constant; the second allowed the intercept to vary between mammals/reptiles, keeping their slopes constant; and the final model allowed both slope and intercept to vary between mammals/reptiles. Comparisons between these models allowed the determination of their improvements upon the baseline model. Phylogenetic ANCOVA identified a lambda value of <0.001 between body mass and Qi for reptiles and mammals, and so non-phylogenetically informed GLS and ANCOVA were used to compare these metrics. μCT study of femoral nutrient foramina Nutrient foramina in the femora of Morganucodon and 11 extant mammal taxa of comparable size were imaged using μCT in the X-Ray Micro-Imaging Laboratory, University of Helsinki (Supplementary Table  2 ). We used μCT data of femoral foramina rather than photographic images of the foramina 51 because many of the Morganucodon foramina were filled with sediment, making their photographic measurement difficult or impossible, whereas the differential density of sediment and fossilized bone in μCT data allowed visualization and analysis. Six Morganucodon femur specimens were selected for analysis from the collections of the Cambridge University Museum of Zoology based on their relative completeness. All femora (both fossil and extant) were scanned using a Bruker Skyscan 1272 μCT scanner at 70 kVp source tension and with a 0.5 mm Al filter. For each scan, 1125–1800 angular projections were collected during a 180° rotation, at an exposure time of 1344 ms. One fossil specimen was scanned at an isotropic voxel size of 4 μm (Fig. 7a ) and five at 5 μm. All femora from extant taxa were scanned at 4 μm voxel sizes. The raw projection data were reconstructed using filtered back-projection with the Feldkamp algorithm by the Bruker reconstruction software “NRecon” (Version 1.7.1.0). Analysis of μCT femora data was conducted using the Avizo image analysis software (Version 9.3.0; Thermo Fisher Scientific). Foramina were located using the “Orthoslice” tool, to scroll through transverse CT slices of each femur. Once located, the foramen was imaged in 3D using the “Volume rendering” tool. This allowed the minimum diameter, and so minimum radius (in cm), of the foramen aperture to be assessed and measured using the “3D line measurement” tool. Minimum diameter was used by Seymour et al. 51 , in the study of foramina in small mammals, as it is not always possible to see the direction of nutrient vessel penetration. Resulting foramina radii were used to generate estimates of blood flow index (Qi) following the method outlined in Seymour et al. 51 . If multiple foramina were found in a single specimen, then the radii were summed following the assumption that this represents the total entry/exit potential of nutrient circulation through the femur 51 . This sum radius, or the single radius in specimens with single foramina (both abbreviated to r; measured in cm), was then used along with the length of the femur (L; measured in cm) to estimate Qi as follows (Eq. 4 ): $$Q_{\rm{i}} = r^4/L$$ (4) Although Morganucodon femoral elements are relatively common in the Hirmeriella fissure suite, almost all are incomplete 76 . We here chose the six most complete femoral specimens known from the Pontalun 3 fissure, from where the Morganucodon teeth imaged using PPC-SRµCT originated. Specimens chosen (UMZC EoPC 19_1 to EoPC 19_6) preserve at least two-thirds of the femoral shaft. Length was estimated from the minimum mid-shaft diaphysial width of each element, using the scaling relationship between the two measurements established for a reconstructed Morganucodon femur created by concatenating three incomplete femora by P.G.G. (Eq. 5 ): $${\mathrm{Femur}}\,{\mathrm{length}} = 10.3^\ast {\mathrm{maximum}}\,{\mathrm{mid}} - {\mathrm{diaphysial}}\,{\mathrm{diameter}}.$$ (5) Statistics One-way ANOVA comparison of intra-observer CV between cementochronological studies; Shapiro–Wilk normality test: PPC-SRµCT data W = 1, p = 1, histological data W = 0.93, p = 0.41; test statistics F = 11.12, degrees of freedom (df) = 10, Cohen’s effect size d = 3.13, p = 0.00728. Phylogenetic ANCOVA comparison of PGLS regression slopes for lifespan against body mass in mammals (log10 lifespan = 0.26(log10 body mass) + 0.16; 95% confidence interval (CI) = 0.05; r2 = 0.69) and reptiles (log10 lifespan = 0.26(log10 body mass) + 0.60; 95% CI = 0.08); r2 = 0.46); slopes are statistically similar (F = 0.868, p = 0.352) while means are significantly separated (F = 4.44, df = 529, partial eta squared effect size = 0.32, p = 0.036). PGLS regression of mammalian lifespan against msSMR: log10 msSMR = −0.237(log10 lifespan) – 0.083; 95% CI = 0.07; r2 = 0.59, p < 0.001. PGLS regression of reptilian lifespan against msSMR: log10 msSMR = −0.83(log10 lifespan) – 0.31; 95% CI = 0.255; r2 = 0.43, p < 0.01. PGLS regression of mammalian lifespan against growth constant K; log10 K = −0.692(log10 lifespan) – 1.171; 95% CI = 0.101; r2 = 0.66, p < 0.01. PGLS regression of reptilian lifespan against growth constant K; log10 K = −0.69(log10 lifespan) – 2.523; 95% CI = 0.339; r2 = 0.43, p < 0.01. One-way ANCOVA comparison of OLS regression slopes for Qi against body mass in extant mammals (log10(Qi) = 0.513 × log10(body mass) – 6.104) and non-varanid reptiles (log10(Qi) = 0.685 × log10(body mass) – 8.139); slopes are statistically similar (F = 2, p = 0.16) while means are significantly different (F = 87.6, df = 89, partial eta squared effect size = 0.50, p = 7.4E−15). Reporting summary Further information on research design is available in the  Nature Research Reporting Summary linked to this article. Data availability The tomographic data that support the findings of this study are available from the corresponding authors upon reasonable request. Virtual thin sections used for cementum and bone increment counts are available at the University of Southampton’s PURE data repository as data number D1506 ( https://doi.org/10.5258/SOTON/D1506 ). Physiological and phylogenetic data are from online databases of the Max Planck Institute ( https://www.demogr.mpg.de/longevityrecords/0203.htm ), an online Ecological Archives database ( http://www.esapubs.org/archive/ecol/E084/094/metadata.htm ), the AnAge database ( https://genomics.senescence.info/species/ ), the VertLife online project ( https://vertlife.org ) and the literature (references in Supplementary Data file  3 ) and are provided in Supplementary Tables, as Supplementary Data files, and as a part of the Source data provided with this paper. References 1. Gill, P. G. et al. Dietary specializations and diversity in feeding ecology of the earliest stem mammals. Nature 512, 303–305 (2014). 3. Han, G., Mao, F., Bi, S., Wang, Y. & Meng, J. A Jurassic gliding euharamiyidan mammal with an ear of five auditory bones. Nature 551, 451–456 (2017). 4. Wilson, G. P. et al. Adaptive radiation of multituberculate mammals before the extinction of dinosaurs. Nature 483, 457–460 (2012). 5. Close, R. A., Friedman, M., Lloyd, G. T. & Benson, R. B. Evidence for a mid-Jurassic adaptive radiation in mammals. Curr. Biol. 25, 2137–2142 (2015). 7. Kemp, T. S. The origin of mammalian endothermy: a paradigm for the evolution of complex biological structure. Zool. J. Linn. Soc. 147, 473–488 (2006). 8. Clarke, A. & Pörtner, H. O. Temperature, metabolic power and the evolution of endothermy. Biol. Rev. 85, 703–727 (2010). 9. Koteja, P. Energy assimilation, parental care and the evolution of endothermy. Proc. R. Soc. B 267, 479–484 (2000). 10. Hopson, J. A. The role of foraging mode in the origin of therapsids: implications for the origin of mammalian endothermy. Fieldiana Life Earth Sci. 49, 126–149 (2012). 12. Lovegrove, B. G. The evolution of endothermy in Cenozoic mammals: a plesiomorphic-apomorphic continuum. Biol. Rev. 87, 128–162 (2012). 13. Avaria-Llautureo, J., Hernández, C. E., Rodríguez-Serrano, E. & Venditti, C. The decoupled nature of basal metabolic rate and body temperature in endotherm evolution. Nature 572, 651–654 (2019). 14. Gross, W. Die typen des mikroskopischen knochenbaues bei fossilen stegoccphalen und reptilien. Anat. Embryol. 103, 731–764 (1934). 15. Ray, S., Botha, J. & Chinsamy, A. Bone histology and growth patterns of some nonmammalian therapsids. J. Vert. Paleontol. 24, 634–648 (2004). 16. Shelton, C. D. & Sander, P. M. Long bone histology of Ophiacodon reveals the geologically earliest occurrence of fibrolamellar bone in the mammalian stem lineage. C. R. Palevol. 16, 397–424 (2017). 17. Benoit, J., Manger, P. R. & Rubidge, B. S. Palaeoneurological clues to the evolution of defining mammalian soft tissue traits. Sci. Rep. 6, 25604 (2016). 18. Hillenius, W. J. Turbinates in therapsids: evidence for Late Permian origins of mammalian endothermy. Evolution 48, 207–229 (1994). 19. Benoit, J., Fernandez, V., Manger, P. R. & Rubidge, B. S. Endocranial casts of pre-mammalian therapsids reveal an unexpected neurological diversity at the deep evolutionary root of mammals. Brain Behav. Evol. 90, 311–333 (2017). 21. Rowe, T. B., Macrini, T. E. & Luo, Z. X. Fossil evidence on origin of the mammalian brain. Science 332, 955–957 (2011). 22. Hoffman, E. A. & Rowe, T. B. Jurassic stem–mammal perinates and the origin of mammalian reproduction and growth. Nature 561, 104–108 (2018). 23. Kielan-Jaworowska, Z. & Hurum, J. H. Limb posture in early mammals: sprawling or parasagittal. Acta Palaeontol. Pol. 51, 393–406 (2006). 24. Rey, K. et al. Oxygen isotopes suggest elevated thermometabolism within multiple Permo-Triassic therapsid clades. eLife 6, e28589 (2017). 25. Huttenlocker, A. K. & Farmer, C. G. Bone microvasculature tracks red blood cell size diminution in Triassic mammal and dinosaur forerunners. Curr. Biol. 27, 48–54 (2017). 26. Olivier, C., Houssaye, A., Jalil, N. E. & Cubo, J. First palaeohistological inference of resting metabolic rate in an extinct synapsid, Moghreberia nmachouensis (Therapsida: Anomodontia). Biol. J. Linn. Soc. 121, 409–419 (2017). 27. Crompton, A. W., Owerkowicz, T., Bhullar, B. A. & Musinsky, C. Structure of the nasal region of non–mammalian cynodonts and mammaliaforms: speculations on the evolution of mammalian endothermy. J. Vertebrate Paleontol. 37, e1269116 (2017). 28. Rodrigues, P. G. et al. Digital cranial endocast of Riograndia guaibensis (Late Triassic, Brazil) sheds light on the evolution of the brain in non-mammalian cynodonts. Hist. Biol. 31, 1195–1212 (2019). 29. Hayes, J. P. & Garland, T. Jr The evolution of endothermy: testing the aerobic capacity model. Evolution 49, 836–847 (1995). 30. Köhler, M., Marín-Moratalla, N., Jordana, X. & Aanes, R. Seasonal bone growth and physiology in endotherms shed light on dinosaur physiology. Nature 487, 358–361 (2012). 31. Kielan-Jaworowska, Z., Cifelli, R. L. & Luo, Z. X. Mammals from the Age of Dinosaurs: Origins, Evolution, and Structure (Columbia University Press, New York, 2004). 32. Hulbert, A. J., Pamplona, R., Buffenstein, R. & Buttemer, W. A. Life and death: metabolic rate, membrane composition, and life span of animals. Physiol. Rev. 87, 1175–1213 (2007). 33. Magalhães, J. P. D., Costa, J. & Church, G. M. An analysis of the relationship between metabolism, developmental schedules, and lifespan using phylogenetic independent contrasts. J. Gerontol. A Biol. Sci. Med. Sci. 62, 149–160 (2007). 34. Werner, J., Sfakianakis, N., Rendall, A. D. & Griebeler, E. M. Energy intake functions and energy budgets of ectotherms and endotherms derived from their ontogenetic growth in body mass and timing of sexual maturation. J. Theor. Biol. 444, 83–92 (2018). 35. Grady, J. M., Enquist, B. J., Dettweiler-Robinson, E., Wright, N. A. & Smith, F. A. Evidence for mesothermy in dinosaurs. Science 344, 1268–1272 (2014). 36. Dong, X., Milholland, B. & Viig, J. Evidence for a limit to human lifespan. Nature 538, 257–259 (2016). 38. Naji, S. et al. Cementochronology, to cut or not to cut? Int. J. Paleopathol. 15, 113–119 (2016). Stock, S. R. et al. Cementum structure in Beluga whale teeth. Acta Biomater. 48, 289–299 (2017). Acknowledgements We acknowledge the European Synchrotron Radiation Facility, Grenoble, France for provision of synchrotron radiation facilities on beamlines ID19 and ID16A (project ES152) and thank Peter Cloetens for assistance in using beamline ID16A. We also acknowledge the Paul Scherrer Institut, Villigen, Switzerland for provision of synchrotron radiation beamtime at beamline TOMCAT of the Swiss Light Source (project 20141278). The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 312284 (for CALIPSO). We thank Keijo Hämäläinen for his help in the initial stages of the synchrotron imaging. This project was part-funded by a Natural Environmental Research Council studentship and an Engineering and Physical Sciences Research Council studentship, awarded to E.N. and P.S. (Grant number NE/R009783/1), and we also thank the Academy of Finland for part-funding the project. Thank you to the Natural History Museum London for contributing to travel for P.B. via the Departmental Investment Fund for Earth Sciences. We thank Ginko Investments Ltd. for funding for materials and travel and the University of Bristol Bob Savage memorial fund for travel for E.N. Many thanks to the Natural History Museum, London, University Museum of Zoology, Cambridge and the Finnish Museum of Natural History, Helsinki, Finland for loans of specimens. facilitated by Martha Richter, Rob Asher, Matt Lowe and Martti Hildén. For assistance with laboratory work and materials, we thank Tom Davies, Wendy Dirks, Dani Schmidt, Remmert Schouten, Pedro Viegas, John Cunningham & Duncan Murdock. Discussions: Roger Benson, Chris Dean, Wendy Dirks, Jim Hopson, Fabien Lafuma, Thomas Martin, Rachel O’Meara, Stephen Naji, Tanya Smith and Emily Rayfield. Author information Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ .

## Discover the right solution for your team

The CB Insights tech market intelligence platform analyzes millions of data points on vendors, products, partnerships, and patents to help your team find their next technology solution.

Request a demo

CBI websites generally use certain cookies to enable better interactions with our sites and services. Use of these cookies, which may be stored on your device, permits us to improve and customize your experience. You can read more about your cookie choices at our privacy policy here. By continuing to use this site you are consenting to these choices.