# The first iguanian lizard from the Mesozoic of Africa

Sebastián Apesteguía, Juan D. Daza, Tiago R. Simões, Jean Claude Rage

► To cite this version:

Sebastián Apesteguía, Juan D. Daza, Tiago R. Simões, Jean Claude Rage. The first iguanian lizard from the Mesozoic of Africa. Royal Society Open Science, 2016, 3 (9), pp.160462. 10.1098/rsos.160462. hal-01426066

HAL Id: hal-01426066

https://hal.sorbonne-universite.fr/hal-01426066

Submitted on 4 Jan 2017

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
ROYAL SOCIETY
OPEN SCIENCE

rsos.royalsocietypublishing.org

Research

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

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

CrossMark
*click for updates*

**Cite this article:** Apesteguía S, Daza JD, Simões TR, Rage JC. 2016 The first iguanian lizard from the Mesozoic of Africa. *R. Soc. open sci. 3:* 160462.

http://dx.doi.org/10.1098/rsos.160462

Received: 28 June 2016

Accepted: 22 August 2016

**Subject Category:**

Earth science

**Subject Areas:**

palaeontology/taxonomy and systematics

**Keywords:**

Acrodonta, biogeography, Cretaceous,
Gondwana, phylogeny, Squamata

**Author for correspondence:**

Sebastián Apesteguía
e-mail: sebapesteguia@gmail.com

Electronic supplementary material is available
at http://dx.doi.org/10.1098/rsos.160462 or via
http://rsos.royalsocietypublishing.org.

THE ROYAL SOCIETY
PUBLISHING

# The first iguanian lizard
from the Mesozoic of Africa

Sebastián Apesteguía<fnref n="1" />, Juan D. Daza<fnref n="2" />, Tiago R. Simões<fnref n="3" /> and Jean Claude Rage<fnref n="4" />

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

SA, 0000-0002-0414-0524

The fossil record shows that iguanian lizards were widely distributed during the Late Cretaceous. However, the biogeographic history and early evolution of one of its most diverse and peculiar clades (acrodontans) remain poorly known. Here, we present the first Mesozoic acrodontan from Africa, which also represents the oldest iguanian lizard from that continent. The new taxon comes from the Kem Kem Beds in Morocco (Cenomanian, Late Cretaceous) and is based on a partial lower jaw. The new taxon presents a number of features that are found only among acrodontan lizards and shares greatest similarities with uromastycines, specifically. In a combined evidence phylogenetic dataset comprehensive of all major acrodontan lineages using multiple tree inference methods (traditional and implied weighting maximum-parsimony, and Bayesian inference), we found support for the placement of the new species within uromastycines, along with *Gueragama sulamericana* (Late Cretaceous of Brazil). The new fossil supports the previously hypothesized widespread geographical distribution of acrodontans in Gondwana during the Mesozoic. Additionally, it provides the first fossil evidence of uromastycines in the Cretaceous, and the ancestry of acrodontan iguanians in Africa.

## 1. Introduction

The phylogenetic position of Iguania (chameleons, agamas, iguanas and New World lizards such as anoles) within Squamata is conflictive between morphological [1,2] and molecular or combined datasets [3–7], being recovered as a nearly diverging branch of the Squamata tree in the former, and nested with snakes and anguimorph lizards (monitors, gila monsters, glass

© 2016 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.

<footnote n="1">¹CEBBAD (CONICET), Fundación de Historia Natural 'Félix de Azara', Universidad Maimónides, Hidalgo 775, 7°p (1405), Buenos Aires, Argentina</footnote>

<footnote n="2">²Department of Biological Sciences, Sam Houston State University, 1900 Avenue I Lee Drain Building Suite 300, Huntsville, TX 77341-2116, USA</footnote>

<footnote n="3">³Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada T6G2E9</footnote>

<footnote n="4">⁴CR2P, Sorbonne Universités, UMR 7207 CNRS, CNRS, Muséum National d'Histoire Naturelle, Université Paris 6, CP 38, rue Cuvier, 75231, Paris cedex 05, France</footnote>
lizards and alligator lizards) in the latter. However, there is a consensus agreement among these studies that acrodontans represent a monophyletic group composed by chameleons and agamas. Acrodontans represent a vast palaeotropical radiation, with extant species living in Africa, Asia, Oceania and parts of southern Europe [8–10]. However, the fossil record of non-acrodontan iguanians dates back to the Early Cretaceous [11] and diversifies into widely distant continental areas by the Late Cretaceous (e.g. [12–17]). The fossil record of acrodontans is considerably more limited.

The oldest putative record of acrodontans (*sensu* Conrad [1] and Gauthier *et al.* [2]) is from the Early–Middle Jurassic of India [18] (but see considerations by Jones *et al.* [19]). Later Mesozoic records come from the Early Cretaceous of China [20], Late Cretaceous of Brazil [21] and the Mid-Cretaceous of Myanmar [22]. Adding to this sparse early fossil record regarding the evolution of the group, phylogenetic and biogeographic studies have also come to conflicting conclusions regarding acrodontan biogeographic origins. Proposed centres of origin and/or dispersal of the crown clades include Eastern Laurasia [23], or Gondwana [18,24], or distinct centres in Eastern Gondwana and Eastern Laurasia [25]. Given that acrodontans alone comprise almost 40% of the extant iguanian species [10], a better understanding of the currently poor knowledge of early acrodontan evolution is paramount.

Africa is a key region of extant acrodontan biodiversity, thus being of major relevance to the study of acrodontan evolution. African acrodontans include large species assemblages of both agamids and chameleons that occupy diverse environments throughout the continent [9,10,25,26]. The oldest unquestionable published record of acrodontans in Africa is from the Oligocene of Egypt [27], with a potential record from the Palaeocene of Morocco [28], and unpublished acrodontans with possible uromastycine affinities from the Eocene of Algeria and Tunisia (J.C.R. and co-workers 2015, unpublished data). Therefore, the available fossil evidence suggests that African acrodontans reached that continent during the Cenozoic, possibly by means of dispersal from Asia [23,25].

Here, we describe an acrodontan fossil species from the early Late Cretaceous Kem Kem beds of Morocco, Northwest Africa. This is the earliest record of any iguanian from the entire Mesozoic of Africa, and brings the record of acrodontans in Africa to the Cretaceous. The new taxon provides valuable clues towards the evolution of acrodontans in Africa, and contributes to our knowledge of early acrodontan evolution worldwide.

## 2. Material and methods

### 2.1. Digital images

Photomicrography and measurements were done using a Nikon D800 36.3 MP digital camera (Tokyo, Japan) at the Muséum National d'Histoire Naturelle (MNHM) in Paris. The specimen is catalogued in the MNHN collection and is registered as collected by René Lavocat. Interpretative illustrations were traced on the photographs using Adobe® Illustrator® CC, 2014.1.1 Release.

### 2.2. Specimens used for comparison

*Agama agama* (FMNH 22190), *Trioceros jacksonii* (AMNH 99984, AMNH 84559), *Pogona barbata* (FMNH 51648), *Physignathus cocincinus* (FMNH 255017), *Uromastyx acanthinura* (AMNH 71836, MCZ 27382), *Uromastyx aegyptia* (AMNH 73160, FMNH 78661), *Saara hardwickii* (UCA.5). The approximate position of the fossil in the skull was estimated using a CT scan of *Uromastyx aegyptia* available from http://digimorph.org.

### 2.3. Standard symbolic codes for institutional resource collections

Abbreviations used follow [29,30] and the following unlisted collection: UCA, University College Anatomy Collection and University of Central Arkansas, USA.

### 2.4. Anatomical nomenclature

Following [31–33], we term sub-dental shelf as the medial projection of the dentary bone that lies ventral to its teeth, and sub-dental ridge as the medial, usually thickened, margin of the sub-dental shelf [32]. In acrodontans (and most rhynchocephalians), the sub-dental shelf is commonly absent, but the dentary bears a ventrally projecting sub-dental ridge which can be quite deep and partially enclose the dorsal half of the Meckelian canal, as in *Uromastyx*, *Chamaeleo* and *Agama*.
## 2.5. Phylogenetic procedures

Character scores from the new taxon and *Gueragama sulamericana* [21] were added to a dataset comprised of morphological and molecular characters [34], and analysed under traditional maximum-parsimony (MP), implied weighting maximum-parsimony (IWMP), and Bayesian inference. Although a larger morphological dataset is available [25], it contains far less mandibular characters than [34], and only one character would have been scored for the new taxon described herein. The resulting dataset includes 30 taxa, with 22 mandibular morphological characters and 1968 molecular characters of mitochondrial data from ND2 and adjacent loci. Out of these, 544 molecular characters were identified as invariable sites (27% of total number of sites) and were removed for the MP analyses. Therefore, 22 morphological and 1424 molecular characters (total of 1446 characters) were analysed under MP. The MP analyses were run in TNT v. 1.1 64-bit (no taxon limit) [35] using 1000 RAS and tree search with TBR (no gains in the sampling of local optima were obtained using the New Technologies algorithms). IWMP is designed to downweight highly homoplastic characters [36,37] and was also used to analyse the dataset in TNT, using a $K=3.0$. A third analysis using Bayesian inference was conducted in MrBayes v. 3.2.5 [38]. Morphological data were analysed under two distinct model parameters for morphological data: $Mk+F81+\Gamma$ and $Mk+F81+LN$, whereas the molecular data were analysed under a $GTR+I+\Gamma$ model (the best substitution model for this molecular dataset [34]). All remaining parameters were set with the default options of MrBayes. The number of generations was set to 1000 000, with six chains and three swaps per cycle (sample frequency = 1000 and burn in fraction of 0.25). The Bayes factor based on harmonic means of the model likelihoods [39,40] was calculated in order to evaluate whether a gamma ($\Gamma$) or a lognormal (LN) distribution had a better fit to the morphological data. As recent studies indicate a LN distribution of rate variation might be preferred to the most widely used gamma distribution [41,42]. However, no significant differences were observed—2 $log_e(B_{10}) < 2$—and therefore only the analysis using a gamma distribution (default option) is reported herein.

In datasets with relatively small numbers of taxa (less than 50 taxa—such as the one herein), Bayesian inference has been consistently demonstrated to provide greater accuracy and efficiency in finding the correct trees when evolutionary rates vary among branches for molecular characters [43–45], which is a biologically sound assumption for both morphological and molecular characters [46–50]. More recently, it has also been demonstrated that Bayesian outperforms traditional parsimony for discrete morphological characters at variable evolutionary rates, with and without missing data [51,52], although with some potential loss of resolution [52]. Additionally, Bayesian analyses allow explicit considerations regarding variations in evolutionary rates and branch lengths among characters, and between the morphological and molecular partitions of the dataset [38,53–55]. For the above reasons, we illustrate the result of the Bayesian inference analysis as our preferred result (based on methodological reasoning). However, for comparative purposes, we also report on the results obtained from both the MP analyses (see also the electronic supplementary material).

# 3. Results

## 3.1. Systematic palaeontology

Squamata Oppel, 1811

Iguanomorpha Sukhanov, 1961

Iguania, Cope, 1864

Acrodonta, Cope, 1864

†*Jeddaherdan* gen. nov.

Type species †*Jeddaherdan aleadonta* sp. nov.

Figures 1 and 3.

Included species †*Jeddaherdan aleadonta*

### 3.1.1. Etymology

*Jeddaherdan* refers to its close relationships with *Uromastyx*, meaning in the Amazigh berber language from Morocco, grandfather (jeddi) of *Uromastyx* (aherdan); ‘aleadonta’, meaning ‘dice teeth’ in reference to the cube-like dentition.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** The estimated silhouette of the skull of MNHN.F.MRS51.1 is based on *Uromastyx aegyptia* (FMNH 78661). Scale bar equals 5 mm.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Position of the type locality (Gara Tabroumit) indicated on a satellite image of Northwest Africa.</figcaption></figure>

### 3.1.2. Holotype

MNHN.F.MRS51.1: isolated left partial mandible with teeth (figures 1 and 3).

### 3.1.3. Type locality and stratigraphy

Cenomanian (Late Cretaceous) beds from the Kem Kem region of Southeastern Morocco, Gara Tabroumit (figure 2). These deposits have superbly preserved fossils, including fishes, amphibians, crocodyliforms and dinosaurs [56–59]. The expeditions carried on by René Lavocat (1909–2007) resulted in hundreds of specimens collected from 1948 to 1951, preserved in the collections of the MNHN. They belong to the three Kem Kem localities worked by Lavocat, (Gara Tabroumit, Kouah Trick and Gara Sbaa). The material described herein was among the fossils collected from the first locality, and had been originally misinterpreted as a fish jaw in the collection. No specific data for the age of Gara Tabroumit is available, so the age for the specimen can range between the base of the Cenomanian to the Early Turonian (ca 90–100.5 MYA) [59].
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Holotype of *Jeddaherdan aleadonta* (MNHN.F.MRS51.1) in (a) dorsal; (b) lateral; (c) medial; (d) ventral views. Scale bar equals 5 mm. Abbreviations follow Evans [31]. An.ft, angular facet; A.pr, alveolar process; Co.ft, coronoid facet; Co.pr, coronoid process; D.su, dentary sulcus; m.f, mental foramen; Mk.g, Meckel's groove; pv.pr, posteroventral process; San.ft, surangular facet; sd.cr, sub-dental crest; sd.sh, sub-dental shelf; si.cr, sinusoidal crest; tt.w, tooth wear; Part.ft, prearticular facet.</figcaption></figure>

### 3.1.4. Diagnosis

*Jeddaherdan aleadonta* can be distinguished from other acrodontan lizards by the following combination of features: posterodorsally ascending coronoid process of dentary; sub-dental ridge projecting ventrally to mid-height of the Meckelian canal posteriorly; posterior dentary teeth very closely spaced, labiolingually compressed and mesiolabially oriented; dentary teeth with vertically straight mesial and distal margins,
and almost flat apices, bearing mesiodistally oriented sinusoidal crests; the dentition also lacks any
ornamentation (striations, flutings or facets) and bears an apical narrow groove extending mesiodistally.

**3.1.5. Description**

The fossil is represented by an incomplete left dentary lacking the symphysial and posteriormost
portions. The fossil includes the facets for articulation with the surangular, angular and coronoid bones.
The dentary is a stout bone with a straight posterior ventral margin; it bears at least four teeth, and has a
thick and well-developed posterodorsal extension of the coronoid process.

The teeth have long and wide tooth bases and look cubic in lateral view with a flat apical surface. In dorsal view, the contour of the teeth is parallelogramic, thus having an overall great similarity to the teeth of *Uromastyx* and *Saara* (uromastycines). Despite this the similarity, the relative size of the teeth is much greater than in extant uromastycines (by comparing with the depth of the dentary and size of the coronoid process). Therefore, the total number of dentary teeth is expected to be less than usually observed in extant uromastycines (between 14 and 17 in the adults of the observed specimens). Additionally, the consecutive teeth are interlocked, with the distolingual side of each tooth in parallel to the mesiolabial side of the tooth immediately posterior to it (providing a mesiolabially inclined configuration of the teeth). This interlocking is present for almost half the length of each tooth mesiodistal length, and defines an oblique and minimal interdental space. The teeth have a well-defined apical sinusoidal crest that is directed mesiolabially. Among extant acrodontan lizards, an apically positioned sinusoidal crest also occurs in uromastycines, but it is directed mesiolingually. The posterior dentition is apicolingually placed on the dentary labial wall, thus the dentary has a more prominent alveolar process on the labial rather than on the lingual side, and the tooth bases extend lingually to the dentary dorsal margin (a condition partially comparable in terms of tooth placement on the dentary to the terms 'pleuroacrodont' [18] or 'sub-pleurodont' [60]). The latter trait is seen in most agamids, whereas a strictly apically placed posterior dentition is most common (and perhaps exclusive) to chamaeleonids ([34] and T.R.S. 2015, personal observation). The posterior dentition is ankylosed to the dentary sub-dental shelf, and a comparatively reduced degree of fusion also occurs to the labial wall of the dentary. No resorption pits are visible, which indicates tooth replacement was absent, at least in the ontogenetic stage represented by this specimen (see below).

The labial surface of the teeth is inclined and has subtle furrows, which may be caused by diagenetic
processes, whereas the lingual side is vertical and smooth. Wearing is evident on the alveolar processes
on the labial side, marking where the lingual surface of the maxillary teeth occluded. According to the
classification of teeth by Jones [61], the general morphology of the teeth in *J. aleadonta* is intermediate
between the cut and slice, and between the grind and shred morphotypes, with a tooth shape similar
to *Uromastyx* and *Saara*. This dental morphotype might be indicative of an omnivorous diet with an
important herbivorous component, as indicated by the elongate mesiodistal axis of the teeth, and an
apical sharp crest that could be used for shredding plant material. However, the lack of extensive tooth
wear as that observed in old *Uromastyx* individuals suggests that the holotype of *J. aleadonta* represents a
young specimen and/or that it fed on plant material much softer than that consumed by *Uromastyx*.

In the labial side of the dentary, there are at least two mental foramina. Also, on the labial side, the facets for the angular and surangular bones are preserved, with the angular facet extending anteriorly to the level of the midpoint of the penultimate tooth. The preserved posterior portion of the dentary has a medially open Meckel's canal. The ventral surface of the mandible is broken close to the base of the coronoid process, and at the anteriormost preserved region of the dentary. Therefore, it cannot be determined with precision if the posteroventral margin of the dentary curved inwards posteriorly (as in chamaeleonids and most agamines) in these regions. However, at the level of the angular facet of the dentary, the dentary is better preserved and the labial wall is vertically straight (as in leiolepidines). This suggests that, if any, the degree of inward inflection would be minor.

Attaching lingually to the sub-dental crest is a facet for the anteromedial process of the coronoid bone.
This facet indicates that the coronoid abutted the sub-dental crest, right below a faint sulcus, or dental
groove [62,63], that extends below the tooth bearing section of the dentary. The facet for the anteromedial
process of the coronoid extends anteriorly to the level of the distal portion of the last dentary tooth, as
also seen in taxa such as *Agama agama*. We discard the possibility of a large splenial in *J. aleadonta*, because
there is no indication of a splenial facet on the sub-dental crest; therefore, the splenial in *J. aleadonta* was
absent or reduced. Had this bone been present, it would have been positioned mostly ventral to the sub-
dental crest, having a limited or no suturing contact with the latter structure, as is observed in many
extant acrodontans (e.g. *Agama*, *Leiolepis* and *Uromastyx*.).
# 4. Discussion

## 4.1. Ontogenetic stage

In acrodontans, tooth replacement does not occur as observed in most other lizards. Some replacement occurs early during ontogeny in the anterior tooth series. However, in the posterior series, teeth are continuously added posteriorly to the back of the tooth row without replacement [64–67]. In the holotype of *J. aleadonta*, the observed degree of ankylosis of the posteriormost teeth to the sub-dental shelf, the amount of tooth wearing and the absence of teeth erupting posteriorly, leads to the conclusion that this specimen probably represents an individual in an advanced ontogenetic stage [66,67]. Nevertheless, the reduced degree of fusion of the posteriormost teeth to the labial wall compared to adults of extant acrodontans (and also *G. sulamericana*) raises the possibility that full skeletal maturity was still not attained.

## 4.2. Osteological comparison of *Jeddaherdan* with other lepidosaurs

*Jeddaherdan* shares with most rhynchocephalians and acrodontans some features such as lack of tooth replacement on the posterior tooth series (although replacement is observed in early rhynchocephalians, [68,69]), a coronoid process of the dentary that ascends posterodorsally, and a reduction or loss of the splenial [25,31,70]. However, *J. aleadonta* differs from all known rhynchocephalians and it is similar to most agamid acrodontans in having the following combination of features: apicolingually positioned posterior dentary teeth (tooth bases are placed apically, but also extend lingually to some degree), whereas they are placed entirely apically in most rhynchocephalians and all chamaeleonids; the presence of a very deep Meckel's groove; the absence of a secondary deposition of dentine along the tooth row; the presence of an angular bone that extends well anteriorly relative to the last tooth in medial aspect (revealed by the angular facet, which extends to the level of the anteriormost preserved tooth) and a dentary with a deep facet for the surangular [18,21,31].

*Jeddaherdan* shares some particularly informative features with extant acrodontans, such as the posterodorsally ascending coronoid process of the dentary, the squared-like shape of the posterior dentary teeth in lateral aspect and the almost horizontal apical tooth surfaces. The latter features are all common to *Uromastyx* and closely related taxa (uromastycines, e.g. AMNH 73160, FMNH 78661, MCZ 27382 and UCA.5), but uncommon to *Leiolepis* and other acrodontans (e.g. [25,31,71] and FMNH 22190, AMNH 99984, AMNH 84559, FMNH 51648, FMNH 255017). Among extant uromastycines, *Jeddaherdan* differs from *Uromastyx* and *Saara* by the mesiolabial orientation of its teeth (which are oblique, but mesiolingually oriented in *Uromastyx*), having a shallow sub-dental ridge, and a very elongate anteromedial process of the coronoid. The dorsolateral excavation of the dentary, producing a crest posterodorsally ascending towards the coronoid process, such as observed in *Uromastyx* and *Gueragama sulamericana*, is also absent in the new taxon. Furthermore, when compared with *Gueragama sulamericana*, the only other known acrodontan from the entire Cretaceous of Gondwana, *J. aleadonta* further differs from it in the form of tooth implantation (*J. aleadonta* has a much reduced lingual contribution of its tooth bases), absence of posterior tooth replacement, absence of a wide apically placed groove on the posterior tooth series, and by the presence of sinusoidal tooth apices.

There are marked differences between *Jeddaherdan* and the Palaeogene fossil acrodont *Barbaturex morrisoni* from Myanmar, which also bears uromastycine affinities [72]. Although both taxa have incomplete dentaries, *Barbaturex* has larger dentaries and an estimated SVL of about 1 m [72], whereas our size estimate for *Jeddaherdan* is no more than 10 cm in SVL (based on jaw/SVL ratio of *Uromastyx acanthinura* AMNH-R 71836, and assuming isometric scaling). Some additional morphological differences include: the lack of prominent ridges on the ventral surface of the dentary (very conspicuous in *Barbaturex*); the shape of the posterior dentary teeth, which have flat apical surfaces in *Jeddaherdan* versus triangular in *Barbaturex*; and consecutive posterior teeth arranged in an interlocking pattern in *Jeddaherdan*, instead of juxtaposed teeth as in *Barbaturex*. Additionally, the presence of an angular facet on the jaw of *Jeddaherdan* indicates that element was present, but not fused to the dentary as in *Barbaturex*. Finally, it is very likely that the splenial bone of *Jeddaherdan* was not as developed as in *Barbaturex* (see our inference about this bone above).

Additional important differences are also found between *Jeddaherdan* and other Palaeogene uromastycines, such as *Khaichinsaurus*, *Lentisaurus Lavatisaurus*, *Graminisaurus Agamimus* and *Acrodontopsis* from the Eocene of Mongolia [73], *Creberidentat* from the Eocene of China [74], and *Qianshanosaurus* from the Palaeocene of China [73,75], grouped under a clade termed
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Maximum-parsimony analyses of the combined morphological and molecular datasets. (a) Strict consensus tree (8313 steps) inferred from 53 most parsimonious trees (7620 steps; consistency index = 0.357; retention index = 0.409 each) obtained from the traditional MP analysis. (b) Strict consensus (fit = 632.00542) of the five best fit trees (7631 steps, fit = 625.290572, CI = 0.357, RI = 0.408 in each tree) inferred from the MP analysis under implied weighting. Blue box = Uromastycinae.</figcaption></figure>

Changjiangosauridae by Alifanov [73]. *Qianshanosaurus*, *Khaichinsaurus*, *Lentisaurus* and *Creberidentat* differ from *Jeddaherdan* by having considerably more teeth in the tooth row (at least 18, in *Qianshanosaurus*, 19 in *Creberidentat* and more than 20 in *Khaichinsaurus* and *Lentisaurus*); a very deep dentary posteriorly; a large splenial (in *Qianshanosaurus*), and labiolingually expanded tooth bases (in *Qianshanosaurus* and *Khaichinsaurus*, also possibly in *Creberidentat*) [73,75].

Some of the most similar Palaeogene taxa to *Jeddaherdan* include *Lavatisaurus*, *Graminisaurus Agamimus* and *Acrodontopsis* [73], an indeterminate uromastycine from the Eocene of Kyrgyzstan [60], *Uromastyx europaeus* from the Oligocene of France [76], and another indeterminate uromastycine from the Oligocene of Egypt [27]. All these taxa share with *Jeddaherdan* and most other uromastycines a posterodorsally ascending coronoid process of the dentary, teeth apiclingually placed on the dentary/maxillary labial margin, lower tooth counts (e.g. less than 14–17 as in extant *Uromastyx*), closely spaced tooth bases, teeth at least slightly elongated mesiodistally, and intensive apical tooth wear in most taxa. However, *Jeddaherdan* differs from all of these taxa by its very distinct interlocking and mesiolabially oriented teeth, and the teeth apical sinusoidal crest.

## 4.3. Phylogenetic position

*Jeddaherdan aleadonta* and *Gueragama sulamericana* are based on a single fragmentary dentary specimen each. Therefore, despite the relatively well-sampled amount of mandibular characters for acrodontans in the present dataset, the strict consensus of the traditional MP analysis could not find resolution at the early divergence of acrodontans (figure 4a). It can be determined, however, that *J. aleadonta* (and *G. sulamericana*) are nested within *Acrodonta*.

In the IWMP analysis (figure 4b), the early dichotomy between chamaeleonids and agamids, as obtained by Hutchinson *et al.* [34], is not retrieved, nor a monophyletic *Leiolepidinae sensu Frost & Etheridge* [77]. Instead, *Leiolepis* falls as the sister taxon to *Uromastycinae* (here being inclusive of *Uromastyx*, *Jeddaherdan* and *Gueragama*), and a clade formed by *Chamaeleontidae + Agaminae sensu Frost & Etheridge* [77].

Under Bayesian inference (figure 5), the analysis we consider most appropriate to our dataset (see Material and methods), *Jeddaherdan* and *Gueragama* are nested with *Uromastyx*, as they also are in the IWMP tree (although forming a monophyletic *Leiolepidinae* with *Leiolepis*), and *leiolepidines* are the sister clade to *Chamaeleontidae + ‘agamines’*. The early divergence of *leiolepidines* (either forming a monophyletic group or a Hennigean comb), followed by the divergence of *chamaeleontids* and *agamines*, is a similar result to the one obtained using MP of morphological data by Gauthier *et al.* [2], and
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Bayesian consensus tree drawn from the 1073 trees obtained from the Bayesian inference of the combined datasets. Values on branches indicate clade posterior probabilities, and branches are proportional to their lengths. Green box = *Leiolepidinae*, blue box = *Uromastycinae*.</figcaption></figure>

both ordered and unordered modified versions of that dataset in Simões *et al.* [78]. This overall topology was also obtained using molecular data under maximum-likelihood [4], and in combined dataset analyses of Wiens *et al.* [7] and Reeder *et al.* [3] using MP. All these results differ from another commonly inferred hypothesis of acrodontan relationships, in which the dichotomy between chamaeleontids and all other acrodontans (Agamidae) occurs early in acrodontan evolution, as previously inferred using both morphological and/or molecular data, under multiple tree inference methods [3,6,79–85], including the analysis of the present dataset without *Jeddaherdan* and *Gueragama* [34].

Importantly, under all of the phylogenetic analyses just discussed (and regardless of the position of chamaeleonids), the lineage leading to *Uromastyx* is always found as one of the earliest diverging ones in the phylogenetic history of acrodontans whenever *Uromastyx* or other uromastycines are included [1–4,6,7,79–85]. Therefore, it is implicit that the divergence of uromastycines might have been an ancient one in acrodontan evolution given our current knowledge of the clade. Considering the fossil evidence of acrodontans in the Cretaceous [20,21]—possibly even in the Jurassic [18]—and most molecular clock estimates for the divergence time for acrodontans [5,86–88] and uromastycines [19,85] (also estimated to be in the Cretaceous; but see [72]), the presence of Cretaceous uromastycines should not be unexpected. However, previous fossil findings only provided evidence for uromastycines, or closely related taxa, past the K-Pg boundary [27,60,72,75] with the oldest known uromastycines being from the Palaeocene of China [75]. *Jeddaherdan*, therefore, provides the first fossil evidence for the phylogenetically predicted presence of uromastycines in the Cretaceous.

The phylogenetic position inferred herein for *G. sulamericana*, previously inferred to lie at the stem of Acrodonta [21], is intriguing and of great relevance. *Gueragama* is nested within uromastycines
along with *Jeddaherdan* (supported by features shared with uromastycines regarding the shape of the dentary bone). Different from *Jeddaherdan*, *G. sulamericana* has a very distinct dentition from other known uromastycines (including the presence of resorption pits), and has been inferred to represent plesiomorphic features of that taxon among other acrodontans [21]. The present tree topology (figure 5) indicates that the dentition and tooth replacement mode of *G. sulamericana* (similar to non-acrodontan iguanians), would be the result of a reversal within acrodontans and not a plesiomorphic trait. Nevertheless, the result for *G. sulamericana* could be affected by the reduced number of non-acrodontan iguanians in the dataset used herein, and thus increasing the sampling of the latter could help to better assess the transition between acrodontans and other iguanians. A future updated morphological dataset of acrodontan iguanians containing a diverse sample of well-constructed morphological characters (and with a more complete sampling of non-acrodontan iguanians) shall provide important insights into early acrodontan evolution, and the position of *Gueragama* within the group.

One important consideration in every phylogenetic analysis of very fragmentary fossil taxa is the potential effects of missing data upon their placement on the tree, which may cause fossils to be found more rootward on the tree [89-91]. However, as some studies have previously indicated, even very incomplete fossils can be accurately placed if they preserve key synapomorphies that allow their correct systematic placement by phylogenetic inference methods [92,93]. In the case of *Jeddaherdan* and *Gueragama* in the analyses herein, the potential rootward slippage caused by missing data does not seem to be an issue. Both are found in a crown group (uromastycines) along with an extant taxon (*Uromastyx*), and do not lie at the leiolepidines or acrodontan stem. This might be a consequence of the unique set of traits on the shape of the dentary that are quite conspicuous to uromastycines, and sometimes also in *Leiolepis* (characters 9 and 12 herein), and their tooth shape too (character 4 herein), which is quite different from the condition observed in other acrodontan lineages.

## 4.4. Biogeographic implications

The presence of a second acrodontan in West Gondwana during the Late Cretaceous provides important clues towards the early evolution of the group. *Jeddaherdan* represents the earliest iguanian known from Africa, and one of the few known acrodontans *sensu strictu* from anywhere in the world during the Cretaceous (considering that priscagamids have most recently been proposed to be the sister taxon to acrodontans, and not within the latter [1-3]). The acrodontans *Xianglong* from the Early Cretaceous of China [20] and *Bharatagama* from the Early-Middle Jurassic of India [18] represent the oldest members of the group (but see considerations concerning *Bharatagama* in [19]). *Jeddaherdan* thus represents a biogeographic link between the Asian Cretaceous forms and the recently discovered South American record (*Gueragama*).

Uromastycines and other agamids are currently widespread in northern Africa [8,26,94], with at least 11 uromastycine species endemic to that region ranging from the Atlantic coast of the Sahara Desert to the western margin of the Red Sea [10]. Previous records indicated the presence of acrodontans with affinities to uromastycines in the Oligocene of Egypt [27] and Eocene of Algeria and Tunisia (see above), but no record was known from the Cretaceous. Therefore, *Jeddaherdan* indicates uromastycines occurred in Africa at least as far back as the early Late Cretaceous, about 40-45 million years before the currently known Cenozoic records. Additionally, *Jeddaherdan* fills a biogeographical gap and confirms previous predictions of the widespread distribution of acrodontan iguanians (along with other iguanians) during the Cretaceous [21]. Given our phylogenetic hypothesis presented above (using our preferred topology, from Bayesian inference—see Material and methods), the oldest record of uromastycines, and of acrodontans in western Gondwana, would be close to the time of the final separation between South America and Africa (Albian-Cenomanian) [95,96]. Whether the distribution of Mesozoic acrodontans on both sides of the early Atlantic Ocean represents a dispersal event or the result of vicariance is currently unknown due to the paucity of fossil data.

Most of the currently known terrestrial limbed lizards from the Aptian/Albian to the Maastrichtian ('middle' to Late Cretaceous) of South America come from lower latitudes when compared with similar aged rhynchocephalians [21], which, in turn, were incredibly more abundant than terrestrial lizards in higher latitudes [97,98] in a ratio of 200:1. The discovery of *Jeddaherdan* indicates lizards also inhabited lower latitudes during the Late Cretaceous of Africa, where no sphenodontians are currently known from. Furthermore, the Cenomanian fossil record shows abundant herbivorous rhynchocephalians southwards from the Caiuá Desert in South America [97] and the presence of at least partially herbivorous lizards and no rhynchocephalians in northern Africa. If eastern South American and western African faunas had a similar terrestrial lepidosaurian composition during and immediately
after the break-up of West Gondwana, then it is possible that further corresponding patterns of faunal distribution between Africa and South America may be revealed with the discovery of additional African material.

Data accessibility. The datasets supporting this article have been uploaded as part of the supplementary material. They are also available from the Dryad Digital Repository [99] http://dx.doi.org/10.5061/dryad.pn35j.

Authors' contributions. S.A. designed the project, S.A., J.D.D. and T.R.S. performed the description; T.R.S. performed the phylogenetic analyses; all authors contributed to discussions, writing and production of the figures.

Competing interests. We declare that we have no competing interests.

Funding. The material was found under the project (to S.A.) 'Recherche pour la identité gondwanienne entre les tétrapodes du Crétacé de France: Une comparaison avec la faune du Amérique du Sud', postdoctorant program CNRS UMR 8538—ENS (2009), pour chercheurs étrangères du Mairie de Paris. T.R.S. thanks the Izaak Walton Killam Memorial Scholarship program for a PhD scholarship.

Acknowledgements. We thank Philippe Loubry (MNHN) for taking the photographs of the holotype of *Jeddaherdan*. We thank Elizabeth Glynne for English corrections. We thank Ronan Allain (MNHN), Alan Resetar (FMNH), Lauren Vonnahme, Margareth Arnold, David Kizirian and Christopher Raxworthy (AMNH), Jonathan Losos and José Rosado (MCZ) for specimen access. We thank Noureddine Jalil and Mehdi Mesmoudi Padinha for discussions on the etymology of new species name. We appreciate the information shared by the colleagues of one of the authors (J.C.R.) from Montpellier, Tlemcen, Oran and Tunis which permitted us to cite the unpublished presence of acrodontans in the Eocene of Africa. We thank the Willi Hennig Society for the free availability of the software TNT. The fossil specimen under study herein is currently housed at the *Muséum National d'Histoire Naturelle*, Paris, France.

## References

1. Conrad JL. 2008 Phylogeny and systematics of Squamata (Reptilia) based on morphology. *Bull. Am. Mus. Nat. Hist.* **310**, 1–182. (doi:10.1206/310.1)

2. Gauthier JA, Kearney M, Maisano JA, Rieppel O, Behlke ADB. 2012 Assembling the squamate tree of life: perspectives from the phenotype and the fossil record. *Bull. Peabody Mus. Nat. Hist.* **53**, 3–308. (doi:10.3374/014.053.0101)

3. Reeder TW, Townsend TM, Mulcahy DG, Noonan BP, Wood Jr PL, Sites JW Jr, Wiens JJ. 2015 Integrated analyses resolve conflicts over squamate reptile phylogeny and reveal unexpected placements for fossil taxa. *PLoS ONE* **10**, e0118199. (doi:10.1371/journal.pone.0118199)

4. Townsend T, Larson A, Louis E, Macey RJ. 2004 Molecular phylogenetics of Squamata: the position of snakes, amphisbaenians, and dibamids, and the root of the squamate tree. *Syst. Biol.* **53**, 735–757. (doi:10.1080/10635150490522340)

5. Vidal N, Hedges BS. 2005 The phylogeny of squamate reptiles (lizards, snakes, and amphisbaenians) inferred from nine nuclear protein-coding genes. *C. R. Biol.* **328**, 1000–1008. (doi:10.1016/j.crvl.2005.10.001)

6. Wiens JJ, Hutter CR, Mulcahy DG, Noonan BP, Townsend TM, Sites JW, Reeder TW. 2012 Resolving the phylogeny of lizards and snakes (Squamata) with extensive sampling of genes and species. *Biol. Lett.* **8**, 1043–1046. (doi:10.1098/rsbl.2012.0703)

7. Wiens JJ, Kuczynski CA, Townsend T, Reeder TW, Mulcahy DG, Sites JW. 2010 Combining phylogenomics and fossils in higher-level squamate reptile phylogeny: molecular data change the placement of fossil taxa. *Syst. Biol.* **59**, 674–688. (doi:10.1093/sysbio/syq048)

8. Zug GR, Vitt LJ, Caldwell JP. 2001 Herpetology—an introductory biology of amphibians and reptiles, 630 p. San Diego, CA: Academic Press.

9. Pianka ER, Vitt LJ. 2003 Lizards: windows to the evolution of diversity, 333 p. Berkeley, CA: University of California Press.

10. Uetz P, Hošek J. 2016 The reptile database. 6 June. See http://www.reptile-database.org

11. Gao K, Nessov L. 1998 Early Cretaceous squamates from the Kyzylkum Desert, Uzbekistan. *N. Jb. Geol. Paläont.*, Abh. **207**, 289–310.

12. Gao K, Fox RC. 1996 Taxonomy and evolution of Late Cretaceous lizards (Reptilia: Squamata) from western Canada. *Bull. Carnegie Mus. Nat. Hist.* **33**, 1–107.

13. Estes R, Price LI. 1973 Iguanid lizard from the Upper Cretaceous of Brazil. *Science* **180**, 748–751. (doi:10.1126/science.180.4087.748)

14. Gao K, Norell MA. 2000 Taxonomic composition and systematics of Late Cretaceous lizard assemblages from Ukhaa Tolgod and adjacent localities, Mongolian Gobi Desert. *Bull. Am. Mus. Nat. Hist.* **249**, 1–118. (doi:10.1206/0003-0090(2000)249:0001:tcasol>2.0.co;2)

15. Sigé B, Buscalioni AD, Duffaud S, Gayet M, Orth B, Rage J-C, Sanz J. 1997 Etat des données sur le gisement crétacé supérieur continental de Champ-Garimond (Gard, Sud de la France). *Münchner Geowiss.* Abh. **34**, 111–130.

16. Blain H-A, Canudo J-I, Cuenca-Boscés G, López-Martínez N. 2010 Amphibians and squamate reptiles from the latest Maastrichtian (Upper Cretaceous) of Blasi 2 (Huesca, Spain). *Cretaceous Res.* **31**, 433–446. (doi:10.1016/j.cretres.2010.06.001)

17. Borsuk-Bialynicka M, Moody SM. 1984 Priscagaminae: a new subfamily of the Agamidae (Sauria) from the Late Cretaceous of the Gobi Desert. *Acta Palaeontol. Pol.* **29**, 51–81.

18. Evans SE, Prasad GVR, Manhas BK. 2002 Fossil lizards from the Jurassic Kota Formation of India. *J. Vert. Paleontol.* **22**, 299–312. (doi:10.1671/0272-4634 (2002)022[0299:FLFTJK]2.0.CO;2)

19. Jones ME, Anderson CL, Hipsley CA, Müller J, Evans SE, Schoch RR. 2013 Integration of molecules and new fossils supports a Triassic origin for Lepidosauria (lizards, snakes, and tuatara). *BMC Evol. Biol.* **13**, 208. (doi:10.1186/1471-2148-13-208)

20. Li P-P, Gao K-Q, Hou L-H, Xu X. 2007 A gliding lizard from the Early Cretaceous of China. *Proc. Natl Acad. Sci. USA* **104**, 5507–5509. (doi:10.1073/pnas. 0609552104)

21. Simões TR, Wilner E, Caldwell MW, Weinschütz LC, Kellner AWA. 2015 A stem acrodont lizard in the Cretaceous of Brazil revises early lizard evolution in Gondwana. *Nat. Commun.* **6**, 9149. (doi:10.1038/ ncomms9149)

22. Daza JD, Stanley EL, Wagner P, Bauer AM, Grimaldi DA. 2016 Mid-Cretaceous amber fossils illuminate the past diversity of tropical lizards. *Sci. Adv.* **2**, e1501080. (doi:10.1126/sciadv.1501080)

23. Honda M, Ota H, Kobayashi M, Nabbitahbata J, Yong H-S, Sengoku S, Hikida T. 2000 Phylogenetic relationships of the family Agamidae (Reptilia: Iguania) inferred from mitochondrial DNA sequences. *Zool. Sci.* **17**, 527–537. (doi:10.2108/0289-0003)

24. Macey JR, Schulte JA, Larson A, Ananjeva NB, Wang Y, Petihayagoda R, Rastegar-Pouyani N, Papenfuss TJ. 2000 Evaluating trans-ethys migration: an example using acrodont lizard phylogenetics. *Syst. Biol.* **49**, 233–256. (doi:10.1093/sysbio/49.2.233)

25. Moody SM. 1980 Phylogenetic and historical biogeographical relationships of the genera in the family Agamidae (Reptilia: Lacertilia). *Ann. Arbor. Ml.* University of Michigan.

26. Pianka ER. 2003 Angleheads, calotes, dragon lizards, and relatives. In Grizmex's animal life encyclopedia (eds M Hutchins, JB Murphy, N Schlager), pp. 209–222, 2nd edn. Farmington Hills, MI: Gale Group.

27. Holmes RB, Murray AM, Chatrath P, Attia YS, Simons EL. 2010 Agamid lizard (Agamidae: Uromastycinae) from the lower Oligocene of Egypt. *Hist. Biol.* **22**, 215–223. (doi:10.1080/08912960903302128)

28. Augé M, Rage J-C. 2006 Herpetofaunas from the upper Paleocene and lower Eocene of Morocco. *Ann. Paleontol.* **92**, 235–253. (doi:10.1016/j.annpal.2005.09.001)

29. Leviton AE, Gibbs Jr RH, Heal E, Dawson C. 1985 Standards in herpetology and ichthyology: Part I.
Standard symbolic codes for institutional resource collections in herpetology and ichthyology. *Copeia* **1985**, 802–832.

30. Sabaj Pérez M. 2014 Standard symbolic codes for institutional resource collections in herpetology and ichthyology: an online reference. Version 5.0. 2016. See http://www.asih.org/

31. Evans SE. 2008 The skull of lizards and tuatara. In *biology of the Reptilia* (eds C Gans, A Gaunt, K Adler), pp. 1–347. Ithaca, NY: Academic Press.

32. Klembara J, Böhme M, Rummel M. 2010 Revision of the anguine lizard *Pseudopus laurillardi* (Squamata, Anguidae) from the Miocene of Europe, with comments on paleoecology. *J. Paleontol.* **84**, 159–196. (doi:10.1666/09-033R1.1)

33. Nydam RL, Eaton JG, Sankey J. 2007 New taxa of transversely-toothed lizards (Squamata: Scincomorpha) and new information on the evolutionary history of 'Teids'. *J. Paleontol.* **81**, 538–549. (doi:10.1666/03097.1)

34. Hutchinson MN, Skinner A, Lee MSY. 2012 Tikiguania and the antiquity of squamate reptiles (lizards and snakes). *Biol. Lett.* **8**, 665–669. (doi:10.1098/rsbl.2011.1216)

35. Goloboff PA, Farris JS, Nixon KC. 2008 TNT, a free program for phylogenetic analysis. *Cladistics* **24**, 774–786. (doi:10.1111/j.1096-0313.2008.00217.x)

36. Goloboff PA. 1993 Estimating character weights during tree search. *Cladistics* **9**, 83–91. (doi:10.1111/j.1096-0313.1993.tb00209.x)

37. Goloboff PA, Carpenter JM, Arias JS, Esquivel DRM. 2008 Weighting against homoplasy improves phylogenetic analysis of morphological data sets. *Cladistics* **24**, 758–773. (doi:10.1111/j.1096-0313.2008.00209.x)

38. Ronquist F et al. 2012 MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. *Syst. Biol.* **61**, 539–542. (doi:10.1093/sysbio/sys029)

39. Kass RE, Raftery AE. 1995 Bayes factors. *J. Am. Stat. Assoc.* **90**, 773–795. (doi:10.1080/01621459.1995.10476572)

40. Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey J. 2004 Bayesian phylogenetic analysis of combined data. *Syst. Biol.* **53**, 47–67. (doi:10.1080/1063515049264699)

41. Wagner PJ. 2012 Modelling rate distributions using character compatibility: implications for morphological evolution among fossil invertebrates. *Biol. Lett.* **8**, 143–146. (doi:10.1098/rslb.2011.0523)

42. Harrison LB, Larsson HC. 2015 Among-character rate variation distributions in phylogenetic analysis of discrete morphological characters. *Syst. Biol.* **64**, 307–324. (doi:10.1093/sysbio/syu098)

43. Kuhner MK, Felsenstein J. 1994 A simulation comparison of phylogeny algorithms under equal and unequal evolutionary rates. *Mol. Biol. Evol.* **11**, 459–468.

44. Saitou N, Imanishi T. 1989 Relative efficiencies of the Fitch-Margoliash, maximum-parsimony, maximum-likelihood, minimum-evolution, and neighbor-joining methods of phylogenetic tree construction in obtaining the correct tree. *Mol. Biol. Evol.* **6**, 514–525.

45. Huelsenbeck JP, Rannala B. 1997 Phylogenetic methods come of age: testing hypotheses in an evolutionary context. *Science* **276**, 227–232. (doi:10.1126/science.276.5310.227)

46. Yang Z. 1994 Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. *J. Mol. Evol.* **39**, 306–314. (doi:10.1007/BF00160154)

47. Chang JT. 1996 Inconsistency of evolutionary tree topology reconstruction methods when substitution rates vary across characters. *Math. Biosci.* **134**, 189–215. (doi:10.1016/0025-5564 (09)00172-7)

48. Farris JS. 1966 Estimation of conservatism of characters by constancy within biological populations. *Evolution* **20**, 587–591. (doi:10.2307/2406593)

49. Farris JS. 1969 A successive approximations approach to character weighting. *Syst. Biol.* **18**, 374–385. (doi:10.2307/2412182)

50. Felsenstein J. 1973 Maximum likelihood and minimum-steps methods for estimating evolutionary trees from data on discrete characters. *Syst. Zool.* **22**, 240–249. (doi:10.2307/2412304)

51. Wright AM, Hillis DM. 2014 Bayesian analysis using a simple likelihood model outperforms parsimony for estimation of phylogeny from discrete morphological data. *PLoS ONE* **9**, e109210. (doi:10.1371/journal.pone.0109210)

52. O'Reilly JE, Puttick MN, Parry L, Tanner AR, Tarver JE, Fleming J, Pisani D, Donoghue PCJ. 2016 Bayesian methods outperform parsimony but at the expense of precision in the estimation of phylogeny from discrete morphological data. *Biol. Lett.* **12**, 20160081. (doi:10.1098/rsbl.2016.0081)

53. Felsenstein J. 2004 *Inferring phylogenies*, 664 p. Sunderland, MA: Sinauer Associates.

54. Ronquist F, van der Mark P, Huelsenbeck JP. 2009 Bayesian phylogenetic analysis using MrBayes. In *The phylogenetic handbook* (eds P Lemey, M Salemi, A-M Vandamme), 2nd edn. Cambridge, UK: Cambridge University Press.

55. Yang Z, Rannala B. 2012 Molecular phylogenetics: principles and practice. *Nat. Rev. Gen.* **13**, 303–314. (doi:10.1038/nrg3186)

56. Martin M. 1984 Deux Lepidosirenidae (Dipnoi) crétacés du Sahara, *Protopterus humei* (PRIEM) et *Protopterus prototeroides* (TABASTE). *Palaontol. Z.* **58**, 265–277. (doi:10.1007/BF02986065)

57. Rage J-C, Dutheil DB. 2008 Amphibians and squamates from the Cretaceous (Cenomanian) of Morocco. A preliminary study, with description of a new genus of pipid frog. *Palaeontog. Abt. A Palaeozool.-Strat.* **285**, 1–22.

58. Calvo J, Salgado L. 1995 Rebbachisaurus tessonei sp. nov. a new Sauropoda from the Albian-Cenomanian of Argentina; new evidence on the origin of the Diplodocidae. *GAIA* **11**, 13–33.

59. Cavin L et al. 2010 Vertebrate assemblages from the early Late Cretaceous of southeastern Morocco: an overview. *J. Afr. Earth Sci.* **57**, 391–412. (doi:10.1016/j.jafreesci.2009.12.007)

60. Averianov A, Danilov I. 1996 Agamid lizards (Reptilia, Sauria, Agamidae) from the Early Eocene of Kyrgyzstan. *N. Jb. Geol. Paläontl. Mh.* **12**, 739–750.

61. Jones MEH. 2009 Dentary tooth shape in Sphenodon and its fossil relatives (Diapsida: Lepidosauria: Rhynchocephalia). In *Interdisciplinary dental morphology: frontiers in oral biology* (eds T Koppe, G Meyer, KW Alt). Greifswald, Germany: Karger.

62. Čerňanský A. 2010 A revision of chamaeleonids from the Lower Miocene of the Czech Republic with description of a new species of Chamaeleo (Squamata, Chamaeleonidae). *Geobios* **43**, 605–613. (doi:10.1016/j.geobios.2010.04.001)

63. Smith KT. 2009 Eocene Lizards of the Clade Geiseltalieli from Messel and Geiseltal, Germany, and the early radiation of Iguanidae (Reptilia: Squamata). *Bull. Peabody Mus. Nat. Hist.* **50**, 219–306. (doi:10.3374/014.050.0201)

64. Edmund AG. 1960 Tooth replacement phenomena in the lower vertebrates. *Contr. R. Ont. Mus. Life Sci. Div.* **52**, 1–190. (doi:10.5962/bhl title.52196)

65. Edmund AG. 1969 Dentition. In *Biology of the Reptilia* (eds C Gans, ADA Bellairs, TS Parsons), pp. 117–200. London, UK: Academic Press.

66. Cooper JS, Poole DFG, Lawson R. 1970 The dentition of agamid lizards with special reference to tooth replacement. *J. Zool.* **162**, 85–98. (doi:10.1111/j.1469-7998.1970.tb01259.x)

67. Cooper JS, Poole DFG. 1973 The dentition and dental tissues of the agamid lizard, *Urofastyx. J. Zool.* **169**, 85–100. (doi:10.1111/j.1469-7998.1973.tb04654.x)

68. Evans SE. 1985 Tooth replacement in the Lower Jurassic leptidosaur *Gephyrosaurus bridensis*. *N. Jb. Geol. Paläontl. Mh.* **7**, 411–420.

69. Whiteside DI. 1986 The head skeleton of the Rhaetian Sphenodontid *Diphydontosaurus avonis* gen. et sp. nov. and the modernization of a living fossil. *Phil. Trans. R. Soc. Lond.* **B 312**, 379–430. (doi:10.1098/rstb.1986.0014)

70. Jollie MT. 1960 The head skeleton of the lizard. *Acta Zool.* **41**, 1–64. (doi:10.1111/j.1463-6395.1960.tb00474.x)

71. The Deep Scaly Project. 2010 'Leiolepis bellandii' (on-line), digital morphology. See http://digimorph.org/specimens/Leiolepis_belliana/ (accessed 10 June 2016).

72. Head JJ, Gunnell GF, Holroyd PA, Hutchison JH, Ciochon RL. 2013 Giant lizards occupied herbivorous mammalian ecospace during the Paleogene greenhouse in Southeast Asia. *Proc. R. Soc. B* **280**, 20130665. (doi:10.1098/rspb.2013.0665)

73. Alifanov VR. 2009 New acrodon lizards (Lacertilia) from the Middle Eocene of southern Mongolia. *Paleontol. J.* **43**, 675–685. (doi:10.1134/S0031030109060124)

74. Li J-L. 1991 Fossil reptiles from Hetaoyuan Formation, Xichuan, Henan. *Vert. Palas.* **29**, 199–203.

75. Dong L-P, Evans SE, Wang Y. 2016 Taxonomic revision of lizards from the Paleocene deposits of the Qianshan Basin, Anhui, China. *Vert. Palas.* **54**, 243–268.

76. Rage J-C, Augé M. 2015 Valbro: A new site of vertebrates from the early Oligocene (MP 22) of France (Quercy). III-Imphibians and squamates. *Ann. Paleontol.* **101**, 29–41. (doi:10.1016/j.annpal.2014.10.002)

77. Frost DR, Etheridge R. 1989 A phylogenetic analysis and taxonomy of iguanian lizards (Reptilia: Squamata). *Univ. Kansas Nat. Hist. Mus. Misc. Pub.* **81**, 1–65.

78. Simões TR, Caldwell MW, Palci A, Nydam RL. 2016 Giant taxon-character matrices: quality of character constructions remains critical regardless of size. *Cladistics*, early view. (doi:10.1111/cla.12163)
79. Lee MSY, Caldwell MW. 2000 Adriosaurus and the affinities of mosasaurs, dolichosaurs and snakes. *J. Paleontol.* **74**, 915–937. (doi:10.1017/S00223360000 33102)

80. Amer SA, Kumazawa Y. 2005 Mitochondrial genome of *Pogona vitticeps* (Reptilia; Agamidae): control region duplication and the origin of Australasian agamids. *Gene* **346**, 249–256. (doi:10.1016/j.gene. 2004.11.014)

81. Lee MSY. 2005 Molecular evidence and marine snake origins. *Biol. Lett.* **1**, 227–230. (doi:10.1098/ rsbl.2004.0282)

82. Lee MSY. 2005 Squamate phylogeny, taxon sampling, and data congruence. *Org. Divers. Evol.* **5**, 25–45. (doi:10.1016/j.ode.2004.05.003)

83. Pyron RA, Burbrink FT, Wiens JJ. 2013 A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. *BMC Evol. Biol.* **13**, 93. (doi:10.1186/1471-2148-13-93)

84. Streicher JW, Schulte JA, Wiens JJ. 2016 How should genes and taxa be sampled for phylogenomic analyses with missing data? An empirical study in iguanian lizards. *Syst. Biol.* **65**, 128–145. (doi:10.1093/sysbio/syv058).

85. Zheng Y, Wiens JJ. 2016 Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species. *Mol. Phylogen. Evol.* **94**, 537–547. (doi:10.1016/j.ympev. 2015.10.009)

86. Kumazawa Y. 2007 Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations. *Gene* **388**, 19–26. (doi:10.1016/j.gene.2006. 09.026)

87. Wiens JJ, Brandley MC, Reeder TW. 2006 Why does a trait evolve multiple times within a clade? Repeated evolution of snakeline body form in squamate reptiles. *Evolution* **60**, 123–141.

88. Hugall AF, Foster R, Lee MSY. 2007 Calibration choice, rate smoothing, and the pattern of tetrapod diversification according to the long nuclear gene RAG-1. *Syst. Biol.* **56**, 543–563. (doi:10.1080/ 10635150701477825)

89. Sansom RS. 2014 Bias and sensitivity in the placement of fossil taxa resulting from interpretations of missing data. *Syst. Biol.* **64**, 256–266. (doi:10.1093/sysbio/syu093)

90. Sansom RS, Wills MA. 2013 Fossilization causes organisms to appear erroneously primitive by distorting evolutionary trees. *Sci. Rep.* **3**, 2545. (doi:10.1038/srep02545)

91. Sansom RS, Gabbott SE, Purnell MA. 2010 Non-random decay of chordate characters causes bias in fossil interpretation. *Nature* **463**, 797–800. (doi:10.1038/nature08745)

92. Wiens JJ. 2003 Missing data, incomplete taxa, and phylogenetic accuracy. *Syst. Biol.* **52**, 528–538. (doi:10.1080/10635150390218330)

93. Wiens JJ. 2006 Missing data and the design of phylogenetic analyses. *J. Biomed. Inf.* **39**, 34–42. (doi:10.1016/j.jbi.2005.04.001)

94. Wilms TM, Böhme W, Wagner P, Lutzmann N, Schmitz A. 2009 On the phylogeny and taxonomy of the genus *Uromastyx* Merrem, 1820 (Reptilia: Squamata: Agamidae: Uromastycinae): resurrection of the genus Saara Gray, 1845. *Bonn. Zool. Beitr.* **56**, 55–99.

95. Valença LMM, Neumann VH, Mabesoone JM. 2003 An overview on Callovian-Cenomanian intracratonic basins of northeast Brazil: onshore stratigraphic record of the opening of the southern Atlantic. *Geol. Acta* **1**, 261–275.

96. Blakey RC. 2008 Gondwana paleogeography from assembly to breakup—A 500 m.y. odyssey. *Geol. Soc. Am. Spec. Pap.* **441**, 1–28. (doi:10.1130/2008. 2441(01))

97. Apesteguía S, Novas FE. 2003 Large Cretaceous sphenodontian from Patagonia provides insight into lepidosaur evolution in Gondwana. *Nature* **425**, 609–612. (doi:10.1038/nature01995)

98. Apesteguía S, Agnolin FL, Lio GL. 2005 An early Late Cretaceous lizard from Patagonia, Argentina. *C. R. Palevol* **4**, 311–315. (doi:10.1016/j.crpv.2005. 03.003)

99. Apesteguía S, Daza J, Simões TR, Rage J. 2016 Data from: The first iguanian lizard from the Mesozoic of Africa. Dryad Digital Repository. (doi:10.5061/dryad.pn35j)