Spatiotemporal Distributions of Non-ophidian Ophidiomorphs, With Implications for Their Origin, Radiation, and Extinction

Campbell Mekarski, Michelle; Pierce, Stephanie E.; Caldwell, Michael W.

doi:10.3389/feart.2019.00245

ORIGINAL RESEARCH article

Front. Earth Sci., 24 September 2019

Sec. Paleontology

Volume 7 - 2019 | https://doi.org/10.3389/feart.2019.00245

Spatiotemporal Distributions of Non-ophidian Ophidiomorphs, With Implications for Their Origin, Radiation, and Extinction

Updated

A correction has been applied to this article in:

Corrigendum: Spatiotemporal Distributions of Non-ophidian Ophidiomorphs, With Implications for Their Origin, Radiation, and Extinction
1. Read correction

$\r\nMichelle Campbell Mekarski*$ Michelle Campbell Mekarski^1*

Stephanie E. Pierce²

Michael W. Caldwell^1,3

¹Department of Biological Sciences, University of Alberta, Edmonton, AB, Canada
²Museum of Comparative Zoology and Department of Organismic Biology, Harvard University, Cambridge, MA, United States
³Department of Earth and Atmospheric Sciences, University of Alberta, Edmonton, AB, Canada

Non-ophidian ophidiomorphs, colloquially referred to as ‘dolichosaurs,’ are small-bodied aquatic lizards that lived in shallow seaways, rivers, and reef environments during the Late Cretaceous. Preservational, geographic, and taphonomic biases in this group make trends in biodiversity difficult to assess. This is exemplified by the fact that the majority of the described species are monotypic and known only from single specimens, imparting very little information on morphological or spatial variation. Here we present a revision of the spatial and temporal distributions of non-ophidian ophidiomorph lizards (‘dolichosaurs’) from Cretaceous sediments worldwide. The fossil record of dolichosaurs begins in the Valanginian (Early Cretaceous). The late Early Cretaceous records are sparse but suggest a wide geographic distribution spanning the Tethys and Western Pacific. This is followed by a dense Cenomanian record from Tethyan and British deposits, and rarer specimens from North America. Though there is a substantial drop in the number of specimens recorded from the Turonian–Maastrichtian, these rare occurrences represent the largest geographical distribution of dolichosaurs: spanning Europe, North America, and South America before going extinct during the end-Cretaceous mass extinction. These occurrences indicate that ophidiomorphs most likely originated in the Jurassic Tethys and continued to radiate spatially and phylogenetically until the end of the Mesozoic, showing much more temporally and environmentally diverse patterns than previously indicated.

Introduction

The early Late Cretaceous is undoubtedly one of the most interesting periods in the evolution of squamates. Primarily terrestrial throughout their evolutionary history, Squamata underwent a unique and spectacular aquatic radiation into the marine realm during the Late Cretaceous. This radiation was accomplished primarily by two groups: Mosasauroidea (aigialosaurs and mosasaurs; see Caldwell, 2012) and Ophidiomorpha (snakes and dolichosaurs; see Palci and Caldwell, 2010). Dolichosaurs (Figure 1) represent a phylogenetically contentious group of elongate, limb reduced, marine lizards closely aligned to mosasaurs and snakes (Lee and Caldwell, 2000; Pierce and Caldwell, 2004). They are small—less than a meter in length—and are generally considered to be semi-aquatic piscivores (Camp, 1923; Caldwell, 1999, 2000; Caldwell and Cooper, 1999; Lee and Caldwell, 2000; Caldwell and Dal Sasso, 2004; Palci and Caldwell, 2007, 2010). Dolichosaurs are diagnosable by their elongate necks (>10 cervical vertebrae), elongate trunks (>35 presacral vertebrae), zygosphene-zygantra articulations throughout the entire presacral region, and a reduction in forelimb elements. Historically grouped into a monophyletic family, Dolichosauridae, which was established for the species Dolichosaurus longicollis (Kramberger, 1892), there are currently 15 species (11 genera) that are generally ascribed to the family.

FIGURE 1

Figure 1. A selection of dolichosaur fossils. (A) Adriosaurus skrbinensis (SMNH 2158); (B) Conisasaurus crassidens (BMB 007155); (C) Pontosaurus kornhuberi (MSNM V3662); (D) Dolichosaurus longicollis (BMNH R 49002); (E) Primitivus manduriensis (MPUR NS 161).

Nopcsa (1908, 1923) was the first to propose that dolichosaurs were more closely related to snakes than to other Cretaceous aquatic squamates (the mosasaurs and aigialosaurs). From the start, this hypothesis was contentious, supported by some (e.g., McDowell and Bogert, 1954) but rejected by others (e.g., Féjérváry, 1918; Camp, 1923). The advent of computer generated cladistic analyses in the 1980’s intensified the debate on the relationship of snakes and dolichosaurs. Some recent phylogenetic analyses (Caldwell, 2000; Lee and Caldwell, 2000; Lee and Scanlon, 2002; Pierce and Caldwell, 2004; Palci and Caldwell, 2007, 2010) support Nopcsa’s dolichosaur-snake hypothesis, even supporting the establishment of Ophidiomorpha, to contain exclusively dolichosaurs and ophidians (Palci and Caldwell, 2010). However, instead of finding the Dolichosauridae and the Ophidia to be sister groups, these analyses instead recover ‘dolichosaurs’ as successive sister groups to ophidians. According to this interpretation, dolichosaurs are a paraphyletic assemblage: an evolutionary ‘grade’ along the road to snakes. Since dolichosaurs are potentially critical to the understanding of snake evolution, it is important to recognize aspects of their palaeoecology as a means of understanding the palaeoecology of snake origins. The logic behind such comparisons is based on the idea that if snakes and dolichosaurs are sister taxa then they must share some features of their common ancestor. Marine adaptations might well be on this list.

The dolichosaur fossil record is spotty, with most material coming from the productive Cenomanian-Turonian localities of Europe and the Western Interior Seaway (Pierce, 2003; Jacobs et al., 2005a, b). This bias is primarily due to commercially quarried Adriatic limestone and a culture of collecting ‘curios’ in 18th and 19th century Europe which increased chances of discovery in western and southern Europe (Dadić, 1982). However, it is evident that some combination of preservational, geographic, geological, taphonomic, and depositional biases contributes to the poor record of dolichosaurs worldwide. This is exemplified by the fact that many of the described genera are monotypic and many of the described species are known only from a single specimen. Records of these animals consist of articulated specimens on limestone slabs, or isolated elements (usually vertebrae). Non-deposition, high-energy and/or biologically active environments, erosion of major rock units, and the small size of the animals themselves diminish the likelihood of preservation and discovery, and causes underestimations of diversity and abundance (Brown et al., 2013; Cleary et al., 2015).

Despite these factors, recent finds have dramatically altered traditional views of the stratigraphic distribution and global range of the dolichosaurs (e.g., Haber and Polcyn, 2005; Hontecillas et al., 2015; Paparella et al., 2018). These new occurrence data indicate that our understanding of the evolution of dolichosaurs needs to be reevaluated. Here, we review the stratigraphy and depositional setting of dolichosaur fossils in order to identify patterns in their distribution, generate a hypothesis of dolichosaur spatiotemporal history and create a paleoenvironmental interpretation that would help to explain their evolutionary history and distribution pattern.

Note on definitions: ‘Dolichosaurs’—Several phylogenetic analyses have found dolichosaurs to form a Hennigian Comb as successive sister taxa to snakes (Lee and Caldwell, 2000; Pierce and Caldwell, 2004; Caldwell and Palci, 2010; Palci and Caldwell, 2010). Therefore, the term ‘dolichosaur’ will be used in this study to refer to non-ophidian ophidiomorphs (elongate, limbed—but limb-reduced—Cretaceous marine squamates) as the term is well understood even though it does not include snakes. In a similar fashion, the term ‘aigialosaur’ will be used in this study to refer to non-mosasaurid (plesiopedal and plesiopelvic) mosasauroids.

Spatiotemporal Distribution of Dolichosaurs

Overview

Historically described as a Cenomanian-Turonian aged group from the Tethys and Western Interior Seaway, dolichosaur fossils have recently been reported from well beyond this spatial and temporal range (Figure 2). These records range from articulated specimens to isolated elements (primarily teeth and vertebrae). Most unarticulated fossils originate from outside the Mediterranean Tethys, and are thus more difficult to reliably assign to a genus or family. Uncertainty is expressed where each report is discussed.

FIGURE 2

Figure 2. Paleogeographic distributions of dolichosaurs according to stage of Cretaceous. (A) Hauterivian (132.9-129.4 My); (B) Barremian (129.4-125.0 My); (C) Albian (113.0-100.5 My); (D) Cenomanian (100.5-93.9 My); (E) Turonian (93.9-89.8 Ma); (F) Coniacian (89.8-86.3 Ma); (G) Santonian (86.3-83.6 Ma); (H) Campanian/Maastrichtian (83.6-66.0 My). A, Adriosaurus sp.; Am, Adriosaurus microbrachis; Ask, Adriosaurus skrbinensis; Asu, Adriosaurus suessi; At, Acteosaurus tommasinii; C, Coniasaurus sp.; Cc, Coniasaurus crassidens; Cg, Coniasaurus gracilodens; Cm, Carentonosaurus mineaui; d, dolichosaur; Dl, Dolichosaurus longicollis; Et, Eidolosaurus trauthi; Hg, Haasiasaurus gittelmani; Jt, Judeasaurus tchernovi; Kh, Kaganaias hakusanensis; M, Mesoleptos sp.; Mz, Mesoleptos zendrinii; Pk, Pontosaurus kornhuberi; Pl, Pontosaurus lesinensis; Pm, Primitivus manduriensis; Pr, Pontosaurus ‘from the reef’ - unnamed; p, pythonomorph lizard (because of their proximity to each other, all the localities in France, and the United States are not shown). Paleogeographic maps from Blakey, R., ©2016 Colorado Plateau Geosystems Inc., used with permission.

Hauterivian (Lower Cretaceous); Figure 2A

Japan

The geologically oldest proposed dolichosaur is Kaganaias hakusanensis, described from the Kuwajima Formation (Tetori Group) of central Honshu, Japan (Evans et al., 2006). The Kuwajima Formation has been variously dated from the Neocomian (Lower Cretaceous, 146-130 mya), though an upper Hauterivian age is best supported by biostratigraphy and zircon fission-track dating (see Evans et al., 2006 for discussion). The formation alternates between beds of fine-grained sandstones, mudstones, and coarse sandstones (Isaji et al., 2005). The coarse sandstone layers are found with pebbles, gravels and driftwood, and are interpreted to be channel deposits. The fine sandstones and mudstones are the fossil bearing units, containing in situ tree trunks, and interpreted as inter-channel deposits. These layers are dominated by three facies (Isaji et al., 2005): one representing a peat marsh, one a shallow freshwater lake, and the third a subaerial swamp. The fossil lizards originate from two facies: one containing freshwater aquatic invertebrates, fish, turtles, and choristoderes which is interpreted as a freshwater lake; the second with lower concentrations of freshwater aquatic invertebrates, fish, turtles, and choristoderes, and dominated by terrestrial fauna (lizards, mammals, dinosaurs). The latter depositional environment would have been at least 100 km from the sea, and is interpreted as a vegetated freshwater swamp that was subject to flooding such as could be found on a delta plain or river floodplain.

Barremian (Lower Cretaceous); Figure 2B

Spain

The next oldest record attributed to a dolichosaur comes from the Barremian of Spain (biostratigraphy; Fregenal-Martínez et al., 2017). These isolated vertebrae were initially described as a snake (Rage and Richter, 1994), but additional data on Cretaceous squamates caused the original authors to rescind this referral (Rage and Werner, 1999). It has since been suggested that it is more likely a dolichosaur (Scanlon and Hocknull, 2008). The vertebrae were recovered from near Uña (Cuenca Province, Castilla-La Mancha) in eastern Spain. The site is part of the La Huérguina Formation which preserves a freshwater lacustrine environment (Buscalioni and Fregenal-Martinez, 2010). This rich lithographic limestone locality produces plants, invertebrates, and complete, fully articulated vertebrates including birds, dinosaurs, pterosaurs, amphibians, turtles, crocodiles, and lizards (Sanz et al., 1988; Evans and Barbadillo, 1997, 1999; Ortego et al., 1999; Bolet and Evans, 2012). The locality is interpreted as a seasonal subtropical wetland ecosystem (Buscalioni and Fregenal-Martinez, 2010).

Albian (Lower Cretaceous); Figure 2C

Australia

A single, incomplete vertebra identified as a dolichosaur (Incertae sedis) comes from the latest Albian of Queensland, Australia (Scanlon and Hocknull, 2008). The specimen originates from the Winton Formation, considered to be latest Albian in age (pollen biostratigraphy; Dettmann and Clifford, 2000). The fossiliferous, silty claystone associated with this find has produced extensive fossil remains including plants, insects, mollusks, lungfish, turtles, crocodilians, pterosaurs, and dinosaurs. The formation is an extensive fluvial sequence composed of several fossiliferous horizons. The larger elements associated with the dolichosaur remains include bones that have been reworked and transported, indicating an active system at the time of the dolichosaur’s deposition (Exon and Senior, 1976; Scanlon and Hocknull, 2008).

Cenomanian (Upper Cretaceous); Figure 2D

Overview of Cenomanian Distribution

By far the highest concentrations of discovered dolichosaurs are Cenomanian in origin. All originate from two seaways: the majority from the Tethys and surrounding regions, and the remainder from the Western Interior Seaway.

In the Late Cretaceous, a transgressive phase of the Tethys Sea covered much of Europe (including Western Europe, the Mediterranean region, and Eastern Europe), the Middle East, and northeastern Africa (Hancock and Kauffman, 1979; Jacobs et al., 2005a; Rage, 2013). The flooded continental margins formed widespread shallow seas dotted with islands. This supported the development of extensive patch-reefs, lagoons, and barrier reefs. In Western Europe, from Ireland to Germany, the warmer waters of the Tethys mixed with the cooler waters of the Atlantic in the Anglo-Paris and North Sea Basins, forming a temperate zone (Diedrich, 2010). This caused the deposition of chalks and marls, rather than the limestones that dominate the rest of the Tethyan region (Robinson, 1986; Owen, 1987).

The Western Interior Seaway was an inland sea that split the western from the eastern side of North America, and stretched from the Arctic Ocean to the Gulf of Mexico. It was formed by the subduction and convergence of the Farallon and North American plates (Robinson Roberts and Kirschbaum, 1995). This caused a basin to form down the center of the continent which was subsequently flooded from the north and south beginning in the Barremian-Middle Albian. The two arms connected in the Late Albian, and the seaway persisted until the Mid-Maastrichtian, close to the end of the Cretaceous (Hattin, 1982; Kauffman and Caldwell, 1993).

Tethys: Adriatic Region

In the Late Cretaceous, the Adriatic region was a massive shallow marine carbonate platform that supported rudist and coral reefs. For many centuries, the yellowish-gray limestone has been quarried for building material, which produced huge numbers of fossil plants, invertebrates, fish, and—more rarely—lizards (Stur, 1891).

Slovenia

Komen (or Comen) is located on the Istrian Peninsula. It is one of the richest sites for Late Cretaceous shallow marine squamates, producing both dolichosaurs and aigialosaurs. The Komen limestones are dark gray and thinly laminated mudstone to wackestone. They contain abundant invertebrates, fish, plants, rudists, and squamates, and represent a very shallow, highly restricted shore comprising lagoons separated by rudist reefs (Cavin and Kolar-Jurkovšek, 2000; Jurkovšek and Kolar-Jurkovšek, 2007).

Though most of the outcrops that produced these lizards can no longer be located, the fish and lizard bearing beds around Komen have been dated to the middle to late Cenomanian using fish fossils and other biomarkers such as foraminiferans (Parona, 1926; Langer, 1961; Polšak et al., 1982; Fuček et al., 1990; Jurkovšek and Kolar-Jurkovšek, 2007).

Besides the Lower Chalk of United Kingdom, the Komen beds produced the first described dolichosaurs: Mesoleptos zendrinii (Cornalia, 1852), Acteosaurus tommasinii (von Meyer, 1860), Adriosaurus suessi (Seeley, 1881), Eidolosaurus trauthi (Nopcsa, 1923), Adriosaurus microbrachis (Palci and Caldwell, 2007), and Adriosaurus skrbinensis (Caldwell and Palci, 2010). The same beds have also produced the basal mosasauroids Carsosaurus marchesetti (Kornhuber, 1893) and Komensaurus carrolli (Caldwell and Palci, 2007). The volume of specimens and species from the same locality (likely the Pivor Formation) indicate an unprecedented concentration of small marine squamates sharing this region.

Croatia

Croatia–Roughly 300 km south of Komen, on the Dalmatian Coast of Croatia, is the Island of Hvar (Isola di Lesina). On the north side of this island, between the towns of Stari Grad and Jelsa is an outcrop of fossiliferous Late Cenomanian carbonates spanning the late Cenomanian to early Turonian (biostratigraphy; see Radovčić, 1973, 1975; Herak et al., 1976; Sarı and Özer, 2009). This outcrop is the most likely origin of a myriad of shallow marine squamates (Herak, 1959; Langer, 1961; Herak et al., 1976). Like the fossils from Komen, the dolichosaurs and aigialosaurs discovered here are among the first described representatives of their respective families. Adriosaurus suessi (Seeley, 1881), Pontosaurus lesinensis (Kornhuber, 1873; Kramberger, 1892), Mesoleptos zendrinii (Cornalia, 1852), a new species of undescribed Pontosaurus (Campbell Mekarski and Caldwell, personal observation), and several indeterminate dolichosaur remains make up the described dolichosaur fanua. Among the basal mosasauroids, the monotypic holotypes of Aigialosaurus dalmaticus, Aigialosaurus bucchichi, and Aigialosaurus novaki originate from Hvar (Kramberger, 1892; Kornhuber, 1901). Interestingly, Hvar has produced the single known specimen of P. lesinensis, Pontosaurus sp. nov., A. dalmaticus, A. bucchichi, and A. novaki, indicating a different ecosystem structure than the paleoenvironment at Komen even though the fossiliferous layers containing the lizards are also dated to the late Cenomanian (Starigrad Formation; Marinčić, 1997; Diedrich et al., 2011). The depositional setting at Hvar was a highly restricted shallow marine environment, most likely in very shallow lagoons surrounded by rudist reefs on an inner platform close to the shore (Radovčić et al., 1983; Fuček et al., 1990; Diedrich et al., 2011). The facies containing abundant rudist fossils is also rich in benthic organisms including mussels, clams, worms, echinoids, oysters and other sparsely preserved organisms of various environmental origin including land plants, fishes of the shallow shelf, and cephalopods of the deeper shelf (Radovčić, 1975; Radovčić et al., 1983; Diedrich et al., 2011).

Tethys: Eastern Reaches

Kazakhstan

Far to the east of Komen and Hvar is Kazakhstan, which in the Cenomanian also lay on the northern border of the Tethys. A single isolated vertebra was described from the upper mid-Cenomanian upper phosphatic horizon of the Mangyshlak Plateau (Averianov, 2001). The indeterminate dolichosaur is associated with a diverse assemblage of sharks, cephalopods, gastropods, echinoderms, and bivalves. The depositional environment was shallow marine with terrestrial input and can be correlated to the Anglo-Paris Basin based on ammonites and inoceramid bivalves (Gale et al., 1999).

Tethys: Middle Eastern Region

Lebanon

The ‘Fish Beds’ of the Sannine Formation, Lebanon are the origin of several taxa of Cenomanian small marine squamates. The Sannine Formation was laid down in a low energy, shallow carbonate platform thought to represent a system of lagoons, rudist patch reefs, and small islands (Dal Sasso and Pinna, 1997; Dalla Vecchia et al., 2002). The appropriately named “fish beds” produce mainly teleosts and crustaceans, but also bivalves, plants, and rare tetrapods including turtles, lizards, and limbed snakes (Dal Sasso and Renesto, 1999; Dalla Vecchia et al., 2002). Several Sannine outcrops in Lebanon have produced squamate remains: Nammoura, Hakel, and Hadjula. Al Nammoura (or En Nammoura) is a mid-Cenomanian Lagerstätten located in the Garbour Valley in the north of Lebanon that produced the holotype of Pontosaurus kornhuberi (Dal Sasso and Renesto, 1999; Caldwell and Dal Sasso, 2004; Caldwell, 2006). This specimen is the most complete pythonomorph lizard described to date, preserving the entire skeleton, scalation, and cartilage. This same locality also produced the holotype of the hind-limbed snake Eupodophis (= Podophis) descouensi (Rage and Escuillié, 2000). Dalla Vecchia et al. (2002) proposed that this locality developed in an anoxic, possibly hypersaline lagoon, leading to the stunning preservation of soft tissues. It is dated to the latest middle Cenomanian based on foraminiferan biostratigraphy, the same age or slightly older than the locality at Hakel and Hadjula (Dalla Vecchia et al., 2002). Hakel (or Hâqel) is associated with a small basin on the outer margin of the continental shelf, and has also produced four specimens of Eupodophis descouensi (Rieppel and Head, 2004).

Aphanizocnemus libanensis is a small marine squamate of uncertain affinities from the lower-middle Cenomanian of Lebanon (Dal Sasso and Pinna, 1997; Dal Sasso and Renesto, 1999). The exact locality of its discovery is uncertain, but it is almost certainly from one of the three Sannine outcrops mentioned above.

Palestine

The ‘Ain Yabrud (or ‘Ein Jabrud) locality of Palestine is an early to middle Cenomanian outcrop that has produced several small marine squamates. These include the hind-limbed snakes Pachyrhachis and Haasiophis (Caldwell and Lee, 1997; Tchernov et al., 2000), and the pythonomorph lizard Haasiasaurus (= Haasia) gittelmani (Polcyn et al., 1999, 2003). Haasiasaurus was described as a mosasauroid, but may represent a more basal pythonomorph lineage. An articulated specimen of Mesoleptos can also be attributed to this locality (Lee and Scanlon, 2002). The locality is alternatively assigned to the Amminidav Formation (Chalifa, 1985, 1989) or the Bet-Meir Formation (Lee and Caldwell, 1998; Lee and Scanlon, 2002), of the latest early Cenomanian (ammonite dating; Lewy and Raab, 1976). The depositional environment of ‘Ain Yabrud is a low energy, shallow marine platform interpreted as a bay close to a river mouth or an inter-reef lagoon similar to the depositional environment in Al Nammoura, Hvar, and Komen (Lee and Caldwell, 1998; Polcyn et al., 1999; Jacobs et al., 2005b). The fossiliferous interval preserves terrestrial plants, squamates, turtles, teleost fishes, elasmobranchs, ammonites, bivalves, crustaceans, and echinoderms, often including soft tissue structures (Polcyn et al., 1999).

Palestine and Israel

The Judean Hills near Jerusalem have produced another basal pythonomorph lizard of uncertain affinities: Judeasaurus tchernovi (Haber and Polcyn, 2005). Lithological analyses indicate that this specimen was deposited in a shallow, low-energy marine environment with little to no terrigenous input. The most likely origin of this specimen is the mid-Cenomanian Kefar Sha’ul or Bina Formation, part of the Judea Group (biostratigraphy; Lewy, 1989, 1996). The depositional environment was most likely a restricted basin on a continental shelf with input from the open ocean, and was rich in fish and ammonites (Haber and Polcyn, 2005).

North Sea Basin and Anglo-Paris Basins

United Kingdom

In the early late Cenomanian, much of United Kingdom was submerged under the Tethys and the Northern Atlantic, forming a broad continental shelf under a shallow sea which produced the Chalk Group: massive soft chalk deposits interspersed with thinner marl and harder, lithified chalk layers. The soft chalks represent deeper water settings, while the harder chalks and marls are thought to have formed in a shallower environment nearer to shore (Robinson, 1986; Owen, 1987). The fossils of the Chalk are diverse and well-studied. A variety of invertebrate (bivalves, cephalopods, arthropods, sponges, bryozoans, echinoderms) and vertebrate (teleost fishes, sharks, sauropterygians, ichthyopterygians, dinosaurs, and pterosaurs) remains show a relatively complete picture of the relatively shallow, continental shelf ecosystem (Owen, 1842, 1850). Dolichosaur fossils are rarer, and are found only in Chalk outcrops in Sussex and Kent (southeast of United Kingdom). Coniasaurus crassidens and Dolichosaurus longicollis are both represented by multiple fragmentary remains spanning the West Melbury Marly Chalk Formation (lower-middle Cenomanian) and the Zig-zag Chalk Formation (Upper Cenomanian), to the Plenus Marl Member (upper Cenomanian) (Owen, 1850, 1851). This indicates that both species occupied marine shallow water environments in the Cenomanian. Coniasaurus gracilodens is represented by a single specimen from the Zig Zag Chalk Formation (Lower Cenomanian) of Sussex (Caldwell, 1999).

Germany

Cenomanian Germany sat northeast of the Anglo-Paris Basin, at the southern edge of the North Sea Basin, which was also covered by a shallow sea (Diedrich, 2010). The proximity makes it unsurprising that similar taxa make up the fossil fauna. Fragmentary Coniasaurus crassidens and Dolichosaurus longicollis have been found in Halle/Westphalia in northwest Germany (Diedrich, 1997, 1999) in the Schwarzbunte Wechselfolge Formation which are marked by a series of large ammonite accumulations. Biostratigraphic correlations and isotope anomalies have been used to date this formation and correlate it to sections in United Kingdom, France, and the northern Tethys including the Plenus Marl of United Kingdom (Lehmann, 1999; Diedrich, 2010). Remains of Coniasaurus and Dolichosaurus were recovered from ammonite scour troughs indicating that the lizards may have been using the accumulation of large ammonite (60–80 cm) shells for protection or hunting opportunities (Diedrich, 2010). In between the ammonites, hundreds of marine fossils are reported, including chondrichthyans, cephalopods, gastropods, bivalves, echinoids, asteroids, brachiopods, annelids, and crustaceans. The depositional environment is interpreted as an offshore submarine swell on a rapidly shallowing slope (Diedrich, 1999, 2010).

Atlantic Margin

France

Southwest of the Anglo-Paris basin, lying roughly on the border of the Cenomanian Atlantic Ocean and the Tethys Seaway, is the Aquitaine Basin. This basin currently occupies the south west quadrant of modern France, and in the Cenomanian, was overtaken by the transgression of the Atlantic Ocean. Deposits from this western edge of the Aquitaine Basin have produced isolated material (mostly vertebrae) of Carentonosaurus mineaui, a pythonomorph lizard of uncertain affinities (Rage and Néraudeau, 2004). The specimens come from at least two localities in Western France: the lower upper Cenomanian of Île Madame (Charente Maritime), the upper Cenomanian of La Couronne (Charente), and possibly other locations including the lower Cenomanian of Les Renardičres (Charente Maritime), the basal upper Cenomanian of Roullet-Saint-Estčphe (Charente), and the upper Cenomanian of Douce Amie (Sarthe). The specimens from Île Madame come from two units, both Cenomanian in age (biostratigraphy; Néraudeau et al., 1997; Rage and Néraudeau, 2004). The majority originate from the ‘Dm’ unit, a silty limestone corresponding to a deeper marine environment and containing adundant echnoids and cephalopods (Rage and Néraudeau, 2004; Csiki-Sava et al., 2015). The rest of the specimens are from the ‘C4ms’ unit, a silty unit indicative of very shallow, high-energy marine environment abundant in rudist and stromatoporoid reefal material (Rage and Néraudeau, 2004; Csiki-Sava et al., 2015). Nearby localities with Cenomanian-aged deposits have also produced specimens referred to Simoliophis rochebrunei, a hind-limbed snake (Bardet et al., 2008).

Spain

The La Cabaña Formation (middle?–late Cenomanian) of Asturias, Spain is a microsite that preserves vertebrate remains of a shallow coastal platform and lagoon with continental input (Vullo et al., 2009). The vertebrate fauna includes fish, turtles, crocodilians, plesiosaurs, dinosaurs, pterosaurs, and marine squamates. The mixture of terrestrial and marine fauna is believed to result from a transgressive lag likely containing reworked remains from older middle Cenomanian deposits (dating based on sequence- and biostratigraphy; Wilmsen, 1997; Vullo et al., 2009; Csiki-Sava et al., 2015). The assemblage of La Cabaña compares favorably to that of Charentes (France), which was also part of the Bay of Biscay (Wilmsen, 1997; Csiki-Sava et al., 2015). Vullo et al. (2009) suggest that the Bay formed an area of endemism during the mid-Cretaceous. Among the endemic taxa supporting this hypothesis is the marine squamate Carentonosaurus (Rage and Néraudeau, 2004). Carentonosaurus is the most abundant tetrapod taxon in Charentes and Asturias, but isolated Coniasaurus-like vertebra and Simoliophis-like vertebra have also been identified (Vullo et al., 2009).

Western Interior Seaway

The Western Interior Sea was a continental flooding of the foreland basin which runs longitudinally down the center of North America. The basin, which was caused by the uplift of the Rocky Mountains to the west, was inundated from the north by the Arctic Ocean, and from the south by the Atlantic/Tethys. Little to no polar ice, and higher rates of sea floor spreading caused a sea level that was much higher than it was today. This period lasted from the late Albian (∼100 mya) to the late Maastrichtian (∼66 mya) and was marked by multiple transgressive/regressive events (Kauffman and Caldwell, 1993; Robinson Roberts and Kirschbaum, 1995). This shallow epicontinental seaway, along with the Tethys Sea, has the densest concentration of Cretaceous small marine squamates, though most of the occurrences in the Western Interior are assigned to a single genus: Coniasaurus.

United States

Several Cenomanian formations from the southern United States have produced disarticulated and incomplete remains of the dolichosaur Coniasaurus. The Eagle Ford Group, Greenhorn Formation, Graneros Shale Formation, and Lake Waco Formation all represent marine environments with fragmentary Coniasaurus remains.

The Eagle Ford is the only formation in the Western Interior that has so far produced more than one type of small marine squamate: the dolichosaur Coniasaurus crassidens (Bell et al., 1982; Jacobs et al., 2005a) and the basal mosasauroid Dallasaurus turneri (Bell and Polcyn, 2005). The Eagle Ford spans the early Cenomanian to the latest Turonian, and contacts the Coniacian Austin Chalk (biostratigraphy, correlation, geochemical dating; Denne et al., 2016). The only Cenomanian marine squamate represented is Coniasaurus, which is found in the Tarrant and Britton Formations of Texas (Bell et al., 1982; Jacobs et al., 2005a). The paleoenvironment was shallow, low-energy marine in association with a river delta (Kauffman and Caldwell, 1993; Denne et al., 2016).

The Greenhorn Formation is extensively exposed throughout the central states: Wyoming, South Dakota, Colorado, Nebraska, Kansas, and New Mexico. The Formation spans the mid Cenomanian to mid Turonian (biostratigraphy; Hattin, 1975; Von Loh and Bell, 1998). The higher beds yield derived marine squamates (mosasaurs), while the lower beds (Orman Lake Member) produce Coniasaurus until just below the Cenomanian-Turonian boundary in Kansas (Liggett et al., 2005), Colorado (Shimada et al., 2006), and South Dakota (Von Loh and Bell, 1998). Like the Eagle Ford, the depositional environment is offshore, deep-water, low-energy marine (Von Loh and Bell, 1998). The abundance of Coniasaurus material from this Formation is particularly notable; for example, Shimada et al. (2006) recovered almost 400 isolated teeth from a deposit in southeastern Colorado.

The middle-Cenomanian Graneros Shale Formation contacts the base of the Greenhorn Formation in Kansas. The depositional environment was shallower, and higher energy than the Greenhorn Formation. Coniasaurus remains (including teeth and a partial maxilla) have been found here, but only at the top of the formation, at the contact between the Graneros Shale and Greenhorn (Liggett et al., 2005).

Unlike the previous formations bearing Coniasaurus fossils, the mid-Cenomanian-Early Turonian (dated biostratigraphically, geochemically) Lake Waco Formation of Texas is higher energy, with consistent current action (Denne et al., 2016). Plesiosaur, chondrichthyan, ammonite, and inoceramid fossils have been recovered from its strata (Thurmond, 1968; Denne et al., 2016). The Bluebonnet Member had a relatively shallow depositional environment, likely a foreshore or bar with abundant plant fossils likely representing marsh-growing reeds (Silver, 1963).

Canada

Cumbaa et al. (2006), in a description of a Cenomanian bonebed fauna from Saskatchewan, Canada, report fragmentary vertebral and dental elements that they tentatively identified as?Coniasaurus crassidens. The locality at the Pasquia Hills is a middle-Cenomanian deposit on the northernmost upland of the Manitoba Escarpment, part of the Belle Fourche Member of the Ashville Formation which is dated using biostratigraphic and regional chronostratigraphic (bentonite) markers (Cumbaa et al., 2006). If these remains do represent a coniasaur, it would be the first and only occurrence from Canada. However, Shimada and Ystesund (2007) express their doubt that the specimens represent a dolichosaur, citing personal communications with the original authors. The bonebed also preserves a wide variety of taxa, including chondrichthyans and actinopterygians of several orders, turtles, sauropterygians, and a surprising diversity of birds (Cumbaa et al., 2006). The bonebed appears to have been deposited in shallow, inshore marine waters with a strong tidal influence and may have represented a lagoonal environment (Cumbaa et al., 2006).

Turonian (Upper Cretaceous); Figure 2E

Though the Turonian is nowhere near as rich in dolichosaur species or density as the Cenomanian, it does reflect a wider geographical distribution (Europe, North America, and South America). The drop in diversity could be attributed to sampling bias, or to an extinction event: the Cenomanian-Turonian Boundary Event (CTBE), which caused significant biotic and geologic changes between these two time periods (Lamolda et al., 1994; Caldwell and Cooper, 1999).

Other marine squamates achieve a broader diversity and distribution during this time. Limbed snakes, which were exclusively known from Northern and Southern margins of the Mediterranean Tethys, became known intercontinentally in the Turonian (Gardner and Cifelli, 1999; Rage and Werner, 1999; Caldwell and Albino, 2001). Mosasauroids, also exclusively known from Europe in the Cenomanian, increase in size and distribution and by the mid-Turonian could be found in Africa, Europe, and North and South America (Von Loh and Bell, 1998; Bardet et al., 2003, 2008; Buchy et al., 2005; Schulp et al., 2008; Smith and Buchy, 2008; Longrich, 2016).

Colombia

The Villeta Group near Yaguará, Colombia is an Aptian-Coniacian deposit which contains the La Frontera Formation; a Turonian deposit (based on biostratigraphy) of laminated limestones with well preserved vertebrate remains (Blanco et al., 2004; García Gonzalez et al., 2009). The majority comprised fish, but marine squamates include indeterminate ‘varanoid’ and ‘pythonomorph’ lizard remains, and the basal mosasauroid Yaguarasaurus colombianus (Páramo-Fonseca, 1994, 1997). The Turonian of Colombia represented a maximal transgression flooding the plains of Colombia to the Andes. The depositional environment at Yaguará was shallow, low-energy, and anoxic, likely deposited in an open marine platform to submarine fan setting (Blanco et al., 2004).

Mexico

The early Turonian limestones at Vallecillo, northeast Mexico are part of the Agua Neuva Formation (Buchy et al., 2005) and are dated based on biostratigraphy (Blanco et al., 2001). The Formation does not typically produce macrofossils, except near Vallecillo, which corresponds to a low-energy, outer shelf environment with localized anoxic episodes (Blanco et al., 2001). These restricted preservational conditions allowed a diverse vertebrate fauna—including marine turtles, fishes, sauropterygians, and squamates—to be preserved complete and often with soft tissue impressions and gastric contents (Blanco-Piñón et al., 2002). This is the state of preservation of Vallecillosaurus donrobertoi (Smith and Buchy, 2008). Though this animal is described as an aigialosaur, the designation is based on caudal neural spine length. Otherwise, the incompleteness of the specimen (pelvis, hind limbs, and anterior caudals) and the lack of diagnostic elements allow that this specimen could be an aigialosaur or a dolichosaur.

United States

Several formations in the United States preserve Turonian-aged Coniasaurus material. These include groups that span the Cenomanian–Santonian, indicating a long history of Coniasaurus occupation in the Western Interior Basin.

As mentioned above, the Eagle Ford (spanning the early Cenomanian to latest Turonian of Texas) is the only formation in the Western Interior to produce more than one small marine squamate: the dolichosaur Coniasaurus crassidens (Bell et al., 1982; Jacobs et al., 2005a) and the basal mosasauroid Dallasaurus turneri (Bell and Polcyn, 2005). At the top of the Group is the Arcadia Park Formation, which contacts the Austin Chalk (Coniacian). Like the Britton Formation, the Arcadia Park Shale was deposited on a shallow marine shelf in association with a freshwater delta (Denne et al., 2016). The fauna is represented by the mosasaur Russellosaurus (Bell et al., 1982; Jacobs et al., 2005a; Polcyn and Bell, 2005), chondrichthyans, osteichthyans, ammonites, clams, oysters, and foraminifera (Moreman, 1942; Denne et al., 2016).

The Boquillas Formation (upper middle Turonian, southeastern Texas) is correlated with the Britton and Arcadia Park Formations of the Eagle Ford Group in Central Texas (Udden, 1907; Huffman, 1960) and has produced Coniasaurus material, as well as indeterminate Mosasauroidea and Russellosaurina material, sharks, and bony fishes (Cicimurri and Bell, 1996). Invertebrate fossils found in the formation include ammonites, crinoids, clams, urchins, and oysters (Denne et al., 2016). The Ernst Member, which preserves the fossil lizards, was deposited on the warm, shallow open shelf of a transgressive sea (Huffman, 1960; Cicimurri and Bell, 1996).

The Carlile Shale (which overlies the Greenhorn Formation; Shimada and Bell, 2006) of South Dakota and Kansas has also yielded isolated Coniasaurus remains (Bell and Polcyn, 1996; Shimada and Ystesund, 2007). The Turner Sandy Member of South Dakota (Bell and Polcyn, 1996), and the Fairport Chalk Member of Kansas (Everhart and Darnell, 2004; Shimada and Ystesund, 2007) are the source of the coniasaur material. The Shales were deposited in a very shallow marine environment close to the shoreline and are Upper Turonian in age (biostratigraphy; Hattin, 1982). The vertebrate fauna includes cartilaginous fishes, bony fishes, turtles, mosasaurs, plesiosaurs, crocodyliforms, and a bird (Shimada and Ystesund, 2007).

France

An isolated vertebra of an indeterminate ‘dolichosaurid’ was described by Rage (1989). The vertebra originates from Sainte-Maure, Touraine, in western France. It is derived from the Saumur Tuffeau, dated to the middle Turonian (based on ammonite dating), and representing relatively shallow outer to inner shelf environments (Rage, 1989; Bardet et al., 1998). Isolated aigialosaurid and mosasaurid vertebrae, and remains from a cheloniid turtle have also been found in these deposits (Bardet et al., 1998). The specimen was later compared more favorably to indeterminate pythonomorph vertebrae from France (Houssaye, 2010), thus altering its status to Pythonomorpha incertae sedis.

A set of 25 isolated vertebrae were discovered in a quarry south of Le Paluau, Indre-et-Loire, western France corresponding to the stratigraphical unit ‘Falun de Continvoir’ (Houssaye, 2010): a shallow marine environment rich in bryozoan and bivalve fragment of late Turonian age (Alcaydé, 1975). The unique shape of the vertebrae was used to suggest that these disarticulated elements may represent a new pythonomorph taxon, though due to a lack of clearly distinctive characters, the taxon was not named (Houssaye, 2010). The different size classes recovered suggest that several individuals are represented.

Coniacian (Upper Cretaceous); Figure 2F

United States

Shimada et al. (2007) described what appears to be an unusually large dolichosaurid tooth similar in anatomy to Coniasaurus crassidens. The find originates from an upper Coniacian unit (based on lithographic marker units) in the Smoky Hill Chalk Member of the Niobrara Chalk Formation in Trego County, Kansas. The Smoky Hill Chalk was deposited hundreds of kilometers from shore under warm (temperate to subtropical), calm, slightly hypoxic conditions (Hattin, 1982). This horizon is rich in marine vertebrate fossils, including cartilaginous and bony fish, plesiosaurs, mosasaurs, and marine turtles (see Shimada et al., 2007, Table 1). Due to the isolated nature of the tooth, it is impossible to tell if it is indeed a dolichosaur, but if it is, it would represent the largest dolichosaur described (based on tooth scaling; Shimada et al., 2007). It would also represent the only Coniacian record of dolichosaurs worldwide.

Santonian (Upper Cretaceous); Figure 2G

United States

The Coniacian-aged Coniasaurus discovery of Shimada et al. (2007) is not the youngest reported Coniasaurus. Another occurrence was described by Shimada and Bell (2006) from the same Member and Formation: the Smoky Hill Member of the Niobrara Chalk (western Kansas). A single isolated vertebra diagnosed as Coniasaurus sp. was recovered from layers established as middle Santonian in age based on lithostratigraphic markers and biostratigraphy (Kauffman and Caldwell, 1993). Mentioned above, the Smoky Hill Member represents the Coniacian-Campanian shallow carbonate platform flooded by the Western Interior Sea (Shimada and Bell, 2006). Unlike most other Coniasaurus-bearing units, the Formation was deposited a few hundred kilometers from shore. Within the same horizon are other marine reptiles including plesiosaurs, marine turtles, and mosasaurs, as well as non-marine reptiles including nodosaurid dinosaurs, pterosaurs, and birds (Everhart, 2005; Shimada and Bell, 2006; Shimada and Ystesund, 2007).

Spain

A collection of 10 isolated vertebrae and 16 teeth was discovered in the Calvero de la Higuera complex, near the village of Pinilla del Valle, northwest of Madrid, Spain (Hontecillas et al., 2015). The complex is a series of cave deposits known for Pleistocene mammals and hominids. The 26 specimens originate from karstic structures formed from Santonian aged marine marls, sandstones, and dolomites. During the Quaternary, the area was karstified, and eventually infilled with fluvial sediments and with sediments formed during the degradation of the karst structures (Peréz-González et al., 2010). It is these reworked Quaternary sediments which produce both Cretaceous and Quaternary vertebrate remains, including the squamate material described. The reptile remains (which are found in association with the teeth of chondrichthyans and osteichthyans) are within Late Cretaceous dolomites and sandstones with dolomitic cement (Peréz-González et al., 2010). Based on correlation with other Late Cretaceous outcrops in the surrounding area, they are assigned to the Santonian (Alonso, 1981; Hontecillas et al., 2015). The squamate material was assigned to Carentonosaurus mineaui, which greatly extends the range of this genus, which was previously only known from the Cenomanian.

Campanian-Maastrichtian (Upper Cretaceous); Figure 2H

Argentina

Albino (2000) reported on an isolated vertebra from the La Colonia Formation, on the southeastern margin of the Somún Curá Plateau, Argentina. The age of this Formation is not well established, but based on micropaleontological data and stratigraphic relations, it seems to span the Campanian-Maastrichtian (Pascual et al., 2000). The deposit was suggested to be an estuary, coastal plain, or tidal flat, since marine, fluvial, and terrestrial taxa are all present (Pascual et al., 2000; Gasparini et al., 2015). Among these was a single vertebra that Albino (2000) assigned to ‘Serpentes incertae sedis.’ However, this assignment has since been questioned, with Scanlon and Hocknull (2008) comparing it more favorably to a dolichosaur.

Spain

The “Lo Hueco” site of Cuenca, central Spain, is a late Campanian-early Maastrichtian site dated based on relative stratigraphic position and paleontological content (Pereda-Suberbiola et al., 2002; Barroso-Barcenilla et al., 2009). It has a diverse vertebrate fauna including terrestrial squamates, turtles, crocodiles, fishes, pterosaurs, and dinosaurs; notably an abundance of titanosaurs (Barroso-Barcenilla et al., 2009). This is in addition to a diverse plant flora, and rarer invertebrates including molluscs and gastropods. The site is located in the upper part of the “Margas, Arcillas y Yesos de Villalba de la Sierra” Formation, and likely represents a muddy flood plain near the coast, swampy, crossed by sandy delta channels, and subjected to both brackish and fresh water influence (Barroso-Barcenilla et al., 2009; Csiki-Sava et al., 2015). Two isolated vertebrae were identified as a novel type of pythonomorph lizard, though the taxon was left unnamed due to the minimal amount of represented material (Houssaye et al., 2013).

France

Another isolated vertebra was described from the late Maastrichtian Cassagnau locality (Auzas Marls Formation) of Haute-Garonne, south-west France (Laurent et al., 2002). The authors mention that the undescribed form is known from Campanian and Early Maastrichtian localities in Southern France, in addition to the Late Maastrichtian locality (dated using micropaleontological evidence) described in the paper (Laurent et al., 1999, 2002). The vertebra is incomplete, and as mentioned by Houssaye et al. (2013), it cannot be definitively determined whether it represents a pythonomorph or not. The Late Maastrichtian locality at Cassagnau preserves terrestrial and marine vertebrates including fishes, amphibians, chelonians, crocodiles, dinosaurs, and birds, as well as bivalves, ostracods, and gastropods (Laurent et al., 2002). The depositional environment was brackish, and marginally marine to estuarine (Bilotte, 1980; Lepicard, 1985).

Italy

Primitivus manduriensis is a newly reported pythonomorph lizard originating from upper Campanian-lower Maastrichtian deposits in Apulia, Italy (Paparella et al., 2018). It is the first record of dolichosaurs from the southern Italian Carbonate Platform, and the youngest reliable occurrence of dolichosaurs worldwide (dating based on nanofossils; Sorbini, 1978; Medizza and Sorbini, 1980). From the Late Jurassic to the Late Cretaceous, the Apulian Platform was a shallow sea abundant in fish and small invertebrates, and containing an archipelago-like arrangement of small islands (Citton et al., 2015). The locality from which the fossil was recovered (near Nardò, Lecce, Puglia) is known for its fossiliferous limestones, which were deposited in a shallow inner lagoon on the Apulian Platform (Bossio et al., 2006). The lithology, lack of microfossils, and lack of trace fossils associated with the specimen suggest that it was deposited in warm, anoxic to dysoxic waters (Paparella et al., 2018) in potentially a sabkha-like environment similar to that described from Hvar, Croatia by Diedrich et al. (2011).

Discussion

Oceanographic Context

The radiation of aquatic squamates (including dolichosaurs) is obviously very closely tied to the seaways and oceans they were colonizing. The beginning of the Cretaceous (∼145 mya) is roughly 55 million years after the break up of Pangea, and about 35 my after Gondwana started to separate. The continents were still relatively close together (or connected), global sea levels were relatively low, and the Atlantic was only just maturing into an ocean and beginning to form between South America and southern Africa (Haq et al., 1988; Polcyn et al., 1999; Blakey, 2008; Bardet et al., 2014; Scotese, 2014; Csiki-Sava et al., 2015). The earliest Cretaceous represented a greenhouse period, with high global CO₂ concentrations and ocean temperatures exceeding 30°C in the tropics (Littler et al., 2011). Cold sea-surface temperatures and a global anoxic event dominated the Aptian-early Albian, until increased oceanic crust production and hydrothermal activity through the Albian-early Turonian drove a long-term warming trend caused by CO₂, and a global eustatic sea level rise (Haq et al., 1988; Clarke and Jenkyns, 1999; Davey and Jenkyns, 1999; Leckie et al., 2002; Jenkyns et al., 2012; Csiki-Sava et al., 2015). The resulting changes in ocean circulation, oceanic stratification, and nutrient partitioning lead to a reorganization of planktonic communities that subsequently caused widespread carbonate deposition in the Late Cretaceous (Leckie et al., 2002). These conditions drove the development of extensive foraminiferan chalk deposits in temperate waters (e.g., the Anglo-Paris Basin and the Western Interior Sea) and of extensive rudistid carbonate platforms in more tropical waters (e.g., Caribbean Tethys and Mediterranean Tethys) (Fuček et al., 1991; Davey and Jenkyns, 1999).

By the mid-Cretaceous, the separation of Gondwana was nearly complete. The rapid spreading of the seafloor and the fragmentation of the continents resulted in shallow epicontinental seas bordering many of the major landmasses that acted as corridors between larger water bodies such as the Tethys, Atlantic, and Pacific (see Bardet et al., 2014; Csiki-Sava et al., 2015 for discussion). In Europe, the complex tectonics created a very dynamic archipelago within the Tethys, unlike any other ecosystem of the time. Blocks of continental crust moving into the spreading area of the Tethys created a series of moving and fluctuating islands that altered population dynamics and geological features in their wake (Csiki-Sava et al., 2015). This unique region formed the basis of many unique faunas and was arguably the driving force behind many distinctive radiations (Bardet et al., 2008; Bardet, 2012; Csiki-Sava et al., 2015).

Reconstruction of Cenomanian current patterns shows a circumglobal current flowed westward through the Tethys Seaway and across the Atlantic, forming the ‘SuperTethys’ (Föllmi and Delamette, 1991; Caldwell and Cooper, 1999). This current is thought to be the primary method of dispersal of certain invertebrate taxa. For example, bivalves go through a planktonic stage, during which time their larvae are free floating and subject to dispersal along ocean currents. Rudist bivalves are known from the Lower Cretaceous (Berriasian–Valanginian) of the European Tethys but are not recorded in the Caribbean region of the SuperTethys until the Barremian (Johnson and Kauffman, 1990). This dispersal is thought to have been a result of the westward prevailing current carrying the rudist larvae across the Atlantic. This current is one of the factors believed to have facilitated the explosive radiation of rudist reefs that replaced coral-algal paleocommunities during the Barremian–Cenomanian (Johnson and Kauffman, 1990). In addition to rudist-dominated communities being ubiquitous throughout the SuperTethys, invertebrate faunas from the mid-Cretaceous of Australia and New Zealand have also been shown to correspond to European faunas, implying an unrestricted migration route between these two ecosystems during or before the Cenomanian which would have been accessible to dolichosaurs (Stilwell and Henderson, 2002).

The mid-Turonian represents a maximum transgression, with the highest water levels of any time during the Cretaceous (Haq et al., 1988; Davey and Jenkyns, 1999). After the mid-Turonian, mid-Cretaceous sea levels and temperatures drop (Clarke and Jenkyns, 1999; Jenkyns et al., 2004). Even though cyclic transgressions continued throughout the Late Cretaceous (early Turonian, early Coniacian, middle Santonian, late Campanian; Hancock and Kauffman, 1979), there was an overall trend toward lower sea levels. This period also represents the closing of the Tethys: the center of dolichosaur diversity and evolution. As the Atlantic spread and Africa rotated, several microcontinents (Central Turkey, Western Iran, Apulia) drifted northward (Polcyn et al., 1999; Blakey, 2008). The eventual result was the closing of the Tethys, the formation of the Alps, and the assimilation of southern Europe. Thus, lower sea levels, tectonic activity, fewer carbonate shelves, and lower sea surface temperatures likely negatively influenced the evolution and dispersion of dolichosaurs during the latest Cretaceous and could have been a key factor in their extinction.

Paleoenvironment

Dolichosaurs, though widespread across the world (on every continent except Africa and Antarctica), were primarily distributed within a paleolatitudinal belt between 25°N and 50°N. Within this belt, three major seaways have produced the majority of dolichosaur records: the southern Western Interior Sea, the northern Tethys, and the eastern margin of the Atlantic where it meets the Tethys. No dolichosaur remains have been found north of this paleolatitudinal belt. The southernmost occurrences are isolated remains from Colombia (∼5°N), Argentina (∼40°S), and Australia (∼45°S). Within these paleolatitudes, dolichosaurs are found in a range of depositional environments: from freshwater and marine deposits (though the vast majority are marine) representing tidal flats, lagoons, sabhkas, swamps, reefs, deltas, basins, and open shelf. Their paleoenvironments were generally warm, but still ranged from the more temperate chalks of United Kingdom to the more tropical limestones from southern Europe (Polšak and Leskovšek, 1975). Assuming that the depositional environments reflect ecological occurrence, this indicates a wide range of ecological niches within relatively shallow, warm water. This pattern agrees with hypothesized paleobiological niche of dolichosaurs as shallow, warm water, marine ambush predators (Caldwell, 1999, 2000, 2006; Lee and Caldwell, 2000; Pierce and Caldwell, 2004; Caldwell and Palci, 2010). Optimal preservation seems to have taken place not on the reefs themselves, but in intra-reef lagoons with localized anoxic conditions (Radovčić et al., 1983; Diedrich et al., 2011).

It has been suggested that, like similarly sized modern marine reptiles, dolichosaurs may have been subject to thermophysiological constraints (due to small body size) that limited their habitats to water between 20° and 30°C (Jacobs et al., 2005b). For example, in the Western Interior Sea, Coniasaurus-bearing rocks were deposited in waters ranging from 20 to 30°C, which corresponds well to the temperature range inhabited by modern marine squamates: sea snakes and the marine iguana (Jacobs et al., 2005a). The warm waters of the mid-Cretaceous (late Albian-early Campanian), which are estimated to have reached 20°C in the Arctic (Jenkyns et al., 2004) and 35°C equatorially (Norris et al., 2002), would have provided little barrier to migration, even into high latitudes. Selection for slightly warmer water could explain why dolichosaurs seem to be restricted to the paleolatitudinal belt between 25°N and 50°N.

In comparison to other Cretaceous marine squamates, dolichosaurs occupy a broader range of paleoenvironments than limbed snakes, which exclusively occupied the carbonate platforms of the Mediterranean Tethys, especially the more tropical southern margin (Rage and Escuillié, 2003; Bardet et al., 2008). Spatially, however, dolichosaurs are less diverse than mosasauroids. Though the smaller aigialosaurs were mostly restricted to the Tethys and Western Interior Seaway, derived mosasaurids were cosmopolitan in their distribution, with a diverse range of habitats including open ocean, near shore, and even freshwater environments (Bardet et al., 2008; Makádi et al., 2012). This supports the hypothesis that dolichosaurs may have been thermophysiologically constrained, as the larger mosasaurs (by virtue of their size), would likely have been less restricted by water temperature.

Origin, Radiation, Extinction

Robust hypotheses on the spatiotemporal and environmental origin of dolichosaurs have been precluded by the poor fossil record in the Early Cretaceous and in fact in the Upper Jurassic as well. Three records (an articulated specimen and two isolated vertebrae) from three distant localities on three continents (Japan, Spain, and Australia; Rage and Richter, 1994; Evans et al., 2006; Scanlon and Hocknull, 2008) have justifiably prohibited much spatiotemporal analysis. These early occurrences are equally confounding because all three come from freshwater deposits, while most Cenomanian–Campanian dolichosaurs originate from marine beds. Another confounding factor is the origin of ophidians, the closest sister group of the dolichosaurs. The oldest known snake is Middle Jurassic in age (Caldwell et al., 2015), while dolichosaurs are first documented from Lower Cretaceous sediments. This means that if dolichosaurs do represent a paraphyletic assemblage from which the monophyletic Ophidia originated, then the origin of non-ophidian ophidiomorphs happened during the Middle Jurassic, or perhaps even earlier (Caldwell, 2019). Currently, two models have been proposed for the origin of non-ophidian ophidiomorphs: a Tethyan origin, or a Western Pacific origin, though the possibility of multiple independent lineages cannot yet be ruled out, and thus multiple ancient areas of origin (see for example: Evans et al., 2006; Bardet et al., 2008; Scanlon and Hocknull, 2008). For the sake of argument, we examine the two existing origin models based on our assemblage of all existing data.

Tethyan Origin

Early pythonomorph diversity is highest in the Western Tethys (southern Europe and the Middle East; Bardet et al., 2008). Therefore, the most parsimonious explanation is that this area was the area of origin for all pythonomorphs, including dolichosaurs, sometime during the Jurassic. Shallow carbonate shelves bordering southern Laurasia and eastern Gondwana from the Berriasian to the Barremian (Lower Cretaceous) would have provided an easy route by which small pythonomorph lizards, adapted to shallow marine environments, could have spread to Japan and Australia by the Early Cretaceous (Kear et al., 2017). These shallow epicontinental shelves kept the animals comfortably within their tropical to temperate belt, and likely provided very little obstruction to dispersal. Once established on coastlines bordering the Pacific, independent populations could have adapted to freshwater environments several times, eventually colonizing the widely distant river and lake deposits. High sea surface temperatures in the Southern Hemisphere in the late Albian could explain the presence of a dolichosaur in Australia, at the highest latitudinal occurrence of any platynotan during the mid-Cretaceous (Jacobs et al., 2005a; Scanlon and Hocknull, 2008). A native dolichosaur population remaining within the Tethys would have formed the basis of the explosive radiation in the Cenomanian. An Early Jurassic, circa-Tethyan origin of the Ophidiomorpha fits well with these findings, and time-wise would accommodate the occurrence of dolichosaurs many hundreds of kilometers away in Japan and Australia.

Similar models for the origin and distribution of other marine tetrapods, such as mosasaurs, support this hypothesis. Prior to the Turonian, mosasauroids (aigialosaurs) are exclusively found in Tethyan deposits from the Mediterranean and the Middle East (Kramberger, 1892; Kornhuber, 1893, 1901; Polcyn et al., 1999; Dutchak and Caldwell, 2006, 2009; Caldwell and Palci, 2007). After an initial westward dispersal to the Western Interior, they quickly developed more specialized adaptations to the marine environment, dispersed globally, and by the mid-Turonian could be found in the Western Interior, the Atlantic, the Tethys, and possibly the Pacific; and by the end of the Turonian, they were essentially cosmopolitan (Bardet et al., 2003, 2008, 2014; Kear et al., 2005; Houssaye and Bardet, 2013). Since aigialosaurs are here interpreted to have had similar paleoecological niches as dolichosaurs, it stands to reason that their distribution would have followed a similar route.

The absence of dolichosaur material in Tethyan deposits prior to the Cenomanian is the primary challenge to this model. The lack of pythonomorph lizards in the Early Cretaceous could be simply explained by a sampling bias caused by the limited exposure of certain ages (Csiki-Sava et al., 2015).

Western Pacific Origin

The alternative to the Tethyan origin model, is that of a Western Pacific origin. Evans et al. (2006) hypothesized that the ancestral pythonomorph lizards initially colonized freshwater, before adapting to brackish and finally marine environments. They suggest that Kaganaias could represent an early, freshwater stage of that transition, implying that the high concentration of Cenomanian Tethyan dolichosaurs represents a lineage already well established and adapted to the marine environment, having successfully radiated and migrated from their freshwater origins.

Following an initial invasion of the marine environment somewhere along the Western Pacific margin, a rapid radiation in the upper Early Cretaceous would have had basal dolichosaurs expanding their range westwards to the Tethys, or eastwards across North America and into Europe while making multiple independent freshwater reinvasions in Australia and Spain. The North American radiation would have likely happened later, in the early Late Cretaceous. Several terrestrial groups—including squamates (e.g., Nydam et al., 2007)—may have followed this migration pattern in the Late Cretaceous, indicating that such a model is certainly within the realm of possibility (see Csiki-Sava et al., 2015 for discussion).

Late Cretaceous Radiations and Extinctions

As of the Cenomanian, the distribution pattern of dolichosaurs appears relatively straightforward. In the early to mid-Cenomanian, a rapid radiation and diversification within the European Tethys is reflected in the diversity of the fossil record. A parallel radiation likely had Coniasaurus-like forms expanding their range westward across the Tethys, crossing the northern margin of the narrow Atlantic, and radiating throughout the interior of North America during the maximum Cretaceous transgression (Jacobs et al., 2005a). The speed of their migration was such that by the middle Cenomanian, they were well established up and down the Western Interior Seaway. The speed and route of their migration closely matches that of the rudist bivalves, that spread westward through the SuperTethys during the Barremian–Cenomanian (Johnson and Kauffman, 1990; Caldwell and Cooper, 1999). Whereas the bivalve larvae would have floated on prevailing currents, it is more likely that the actively swimming dolichosaurs followed either the southern coastline of the Atlantic along western Africa to South America, or the northern coastline through the Baltics and along eastern North America (Caldwell and Cooper, 1999). The migration pattern also compares well with that proposed for the Turonian distribution of russellosaurines and aigialosaurs (Bell, 1997; Bardet et al., 2008; Kear et al., 2013).

Interestingly, after the Cenomanian, there is no record of dolichosaurs in Asia or Australasia, both of which had possible pre-Cenomanian records. Several hypotheses exist to explain this: first, that these tentative, fragmentary records are not dolichosaurs at all. Second, it is possible that dolichosaur populations did exist on these continents but have not been found due to sampling bias and their relative rarity. Finally, that these animals suffered from the mass-extinction event associated with the end of the Cenomanian, the Cenomanian-Turonian boundary event (CTBE; Lamolda et al., 1994). The CTBE was a short (∼250 ka) global oceanic anoxic event possibly linked to rapid oceanic plateau formation and/or increased rates of ridge crest volcanism (Lamolda et al., 1994; Kerr, 1998; Davey and Jenkyns, 1999; Leckie et al., 2002; Wan et al., 2003; Turgeon and Creaser, 2008). The event lead to an extensive faunal turnover, especially in the marine realm where it caused the extinction of roughly 27% of marine invertebrates (Turgeon and Creaser, 2008), and an estimated 90% of benthic foraminiferans, and 50–70% of planktonic foraminiferans (Wang et al., 2001). These invertebrate extinctions were particularly important since these taxa were the builders of the carbonate platforms that formed the basis of the dolichosaurs’ ecosystems. The CTBE, which restructured marine ecosystems globally, could easily have affected dolichosaur ecosystems and thus populations, causing an extinction or near-extinction of dolichosaurs outside of their population centers: the Tethys and the Western Interior. The radiation of mosasauroids in the Tethys during the Turonian (Bardet et al., 2008; Kear et al., 2013) likely exacerbated the situation, possibly putting competition and predation pressures on any surviving dolichosaur populations.

More surprising is the absence of any dolichosaurs in the Eastern Tethys after the CTBE, which had been the center of known dolichosaur diversity in the Cenomanian. This absence is likely an artifact of collection bias combined with the rarity of these fossils. However, it is also possible that the Eastern European dolichosaurs suffered extinction associated with the CTBE. They were living and foraging within the extensive carbonate platforms that were eventually submerged below the euphotic zone and inundated by anoxic water during the transgression of the paleo-seaways (Jenkyns, 1991; Gušić and Jelaska, 1993). In Europe and the Middle East, dolichosaurs rarely appear in the fossil record after the CTBE, and related marine squamates, such as aigialosaurs and limbed snakes (i.e., Pachyrhachis) that were occupying a similar ecological niche, also disappear from the fossil record during this period. The combination of the CTBE and the Turonian radiation of mosasauroids could explain why from the Turonian–Santonian, the only records of dolichosaurs are from the Eastern Atlantic Margin/Western Tethys and the Americas (see Polcyn et al., 2014, Figures 1, 3).

The Turonian through Santonian has a sparse though consistent record of dolichosaurs from the Western Interior and Western Tethys, indicating that though these animals persisted, they generally remained a small part of the fauna in their established biogeographic centers. A peak transgression or maximum flooding during the mid-Turonian (Davey and Jenkyns, 1999) could have facilitated a possible South American radiation represented by two occurrences from South America. This first, indeterminate pythonomorph remains from Colombia in the Turonian hint at the possibility of a southward migration starting from the Western Interior, moving down along the island system that was Central America, to the shallow epicontinental sea that covered much of the northern half of South America in the mid-Turonian. A second possible occurrence, a single vertebra from Argentina, would imply the continuation of the southward radiation, which reached the southern coast of South America by the Campanian–Maastrichtian.

A recently described pythonomorph Primitivus manduriensis from Apulia, Italy in the Campanian–Maastrichtian (Paparella et al., 2018) implies that even though the diversity of these animals decreased into the Late Cretaceous, they continued to persist along with the more diverse lineages of Late Cretaceous aquatic reptiles including mosasaurs and sauropterygians. The Maastrichtian finds of Italy and France [see section Campanian–Maastrichtian (Upper Cretaceous); Figure 2H] are the youngest records of non-ophidian ophidiomorphs, implicating the End Cretaceous Mass Extinction as the event that caused the extinction of this group.

Species Distributions and Endemism

The internal relationships of dolichosaurs are not well resolved, therefore it is difficult to discuss the distribution patterns of separate lineages. However, some interesting patterns do emerge surrounding some distinct paleocommunities.

In the late Cenomanian of Europe, four distinct dolichosaur communities appear to exist despite the fact that the Mediterranean Tethys was continuous with the Northern and Western Tethys regions, and the Atlantic ocean. The most diverse dolichosaur fauna is in the Adriatic region, with multiple genera from Slovenia and Croatia (Calligaris, 1988; Bardet et al., 2008). A second, Middle Eastern fauna appears to have some unique genera such as Judeasaurus and Haasiasaursu, but also appears to have some crossover with the Adriatic region, as both faunas have produced Mesoleptos and Pontosaurus fossils. The Anglo-Paris basin has produced abundant fragmentary Coniasaurus and Dolichosaurus fossils originating from the Northern Tethys margin, but which do not appear in the southern or eastern Tethys. Finally, Carentonosaurus seems to be endemic to the Bay of Biscany/Charentes region at the Eastern margin of the Atlantic. The spatial distribution of these marine lizards must have been impacted by their paleoenvironment and paleoclimate. Climatically, it seems that Coniasaurus and Dolichosaurus must have been adapted to the more temperate climate of the Northern Tethys margin, while the Middle Eastern and Adriatic faunas were better suited for warmer, more tropical areas on the southern Tethys. Environmentally, Carentonosaurus seems to have been the only dolichosaur adapted to the Atlantic coastal terrigenous landscape of the Bay of Biscay. The broad carbonate platforms of the Adriatic and Middle East supported the highest diversity, while the cooler chalk deposits of the Anglo-Paris basin were intermediate.

In contrast to most dolichosaur species, specimens attributed to Coniasaurus are more pandemic: being spatially diverse and numerous. Coniasaurus crassidens is recognized from the Cenomanian Chalk of United Kingdom, spanning the early to late Cenomanian. Fossils with similar vertebrae and characteristically swollen teeth are also found in late Cenomanian deposits in Germany, and late Cenomanian to Santonian rocks from Kansas, Colorado, South Dakota, and Texas (see above for specific localities and references). These occurrences are all from temperate waters, which indicates a cooler environmental tolerance for this genus. The fossil record indicates an early Cenomanian origin for the genus, which likely spread to the Western Interior following the prevailing westward current of the SuperTethys throughout the Cenomanian (Caldwell and Cooper, 1999). Coniasaurus fossils have not been recovered from European deposits younger than the CTBE, implying that the European lineage may have been driven to extinction by this event. However, numerous Turonian records in the Western Interior indicate an established lineage that may have persisted to at least the Santonian (based on fragmentary remains from the Coniacian and Santonian of Kansas). A single indeterminate Colombian record in the Turonian hints that the North American lineage even radiated after the boundary event, migrating and establishing new populations in South America. The tentatively assigned dolichosaur vertebra from the Campanian–Maastrichtian of Argentina (Albino, 2000) may attest to the continued evolution and radiation of the North American Coniasaurus lineage.

Overall, the fossil record of dolichosaurs shows a high degree of endemism. This is surprising, as no known physical barriers would have obstructed their potential migration. And yet, dolichosaurs are not the only group that show a surprising amount of characteristic localized faunas. An interesting correlation exists between the multiple endemic populations of dolichosaurs and that of Late Cretaceous pan-chelonioids, certain Tethyan mosasaurs, and the latest Cretaceous plesiosaurs, who all—despite a suite of marine dispersal routes—showed a high degree of specific endemism, even though they were distributed worldwide in shallow marine environments (Hirayama, 1997; Vincent et al., 2011; Bardet, 2012; Kear et al., 2013; Bardet et al., 2014). It is possible that this is related to island biogeography, and that these populations were specialized to the particular ecology surrounding an island or island system. This alludes to the danger posed to dolichosaurs by the CTBE: if these faunas were so specialized to a particular region of the Tethys, then their disappearance could have been caused by changes to, or the loss of, their particularly specialized and unique habitats.

Summary

Previously thought to be limited in space and time to the Tethys and Western Interior Seaway of the mid-Cretaceous, dolichosaurs are now understood to have a more complex history than previously reported. Current knowledge of dolichosaurs indicates a large diversity of species, despite the low number of identifiable specimens. In some areas (Mediterranean region), dolichosaurs are preserved primarily as articulated skeletons, while in others (western Europe, North America) fragmentary remains provide little diagnostic information. Though biogeographic patterns and processes remain obscured, some general patterns are clear. Early pythonomorph diversity is highest in the Western Tethys and this is the likely area of origin of ophidiomorphs, including dolichosaurs. A radiation in the late Early Cretaceous likely had early dolichosaurs expanding their range across the Tethys and making multiple independent freshwater incursions in areas as far away as Japan and Australia. The North American radiation likely happened later, in the early Late Cretaceous. An extinction event at the end of the Cenomanian may have caused the extinction or near-extinction of most dolichosaur lineages, save for an eastern Atlantic fauna, and the North American fauna. These lineages persisted and radiated in the Late Cretaceous, eventually going extinct at the Cretaceous-Tertiary boundary. We expect that future discoveries—especially in areas with tentative and fragmentary dolichosaur remains—will provide new data about the evolutionary and paleogeographical history of pythonomorph lizards. Understanding the occurrence of this group through space and time is critical to the study of their evolution and adaptation, including evolutionary drivers and their aquatic radiation.

Author Contributions

MCM wrote the manuscript and performed the literature review. SP performed the literature review, provided edits on the manuscript, and wrote the part of the sections “Introduction” and “Discussion.” MC critically revised the manuscript, provided funding for the study and longer-term research projects, and supervised the graduate research of MCM and SP that contributed new insights to this manuscript.

Funding

This study was funded in part by an NSERC Discovery Grant (#238458), an NSERC Accelerator Grant (#412275), and a Chair’s Research Allowance to MC. The authors also thank the Vanier CGS for a scholarship to MCM.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank the reviewers JS and SL for their constructive comments, and the editor KB for his support during the review process. Thanks to Ron Blakey for the use of his maps (www.deeptimemaps.com), which were used in the production of Figure 2. This publication is derived from a chapter from the Ph.D. thesis of MCM (Campbell Mekarski, 2017).

References

Albino, A. M. (2000). New record of snakes from the cretaceous of patagonia (Argentina). Geodiversitas 22, 247–253.

PubMed Abstract | Google Scholar

Alcaydé, G. (1975). Notice explicative, carte géologique de la France à 1/50 000, feuille Chinon (486), XVII-23. Orléans: Bureau de recherches géologiques et miničres.