Michael S. Y. Lee

Calibration Choice, Rate Smoothing, and the Pattern of Tetrapod Diversification According to the Long Nuclear Gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.

Snake phylogeny based on osteology, soft anatomy and ecology.

Relationships between the major lineages of snakes are assessed based on a phylogenetic analysis of the most extensive phenotypic data set to date (212 osteological, 48 soft anatomical, and three ecological characters). The marine, limbed Cretaceous snakes Pachyrhachis and Haasiophis emerge as the most primitive snakes: characters proposed to unite them with advanced snakes (macrostomatans) are based on unlikely interpretations of contentious elements or are highly variable within snakes. Other basal snakes include madtsoiids and Dinilysia– both large, presumably non‐burrowing forms. The inferred relationships within extant snakes are broadly similar to currently accepted views, with scolecophidians (blindsnakes) being the most basal living forms, followed by anilioids (pipesnakes), booids and booid‐like groups, acrochordids (filesnakes), and finally colubroids. Important new conclusions include strong support for the monophyly of large constricting snakes (erycines, boines, pythonines), and moderate support for the non‐monophyly of the ‘trophidophiids’ (dwarf boas). These phylogenetic results are obtained whether varanoid lizards, or amphisbaenians and dibamids, are assumed to be the nearest relatives (outgroups) of snakes, and whether multistate characters are treated as ordered or unordered. Identification of large marine forms, and large surface‐active terrestrial forms, as the most primitive snakes contradicts with the widespread view that snakes arose via minute, burrowing ancestors. Furthermore, these basal fossil snakes all have long flexible jaw elements adapted for ingesting large prey (‘macrostomy’), suggesting that large gape was primitive for snakes and secondarily reduced in the most basal living foms (scolecophidians and anilioids) in connection with burrowing. This challenges the widespread view that snake evolution has involved progressive, directional elaboration of the jaw apparatus to feed on larger prey.

ADRIOSAURUS AND THE AFFINITIES OF MOSASAURS, DOLICHOSAURS, AND SNAKES

Abstract The poorly-known, long bodied, limb-reduced marine lizard Adriosaurus suessi Seeley, 1881, is reassessed. Adriosaurus and a number of other marine lizards are known from Upper Cretaceous (Upper Cenomanian-Lower Turonian) marine carbonate rocks exposed along the Dalmatian coast of the Adriatic Sea, from Komen, Slovenia, to Hvar Island, Croatia. A revised vertebral count reveals 10 cervical, 29 dorsal, and at least 65 caudal vertebrae. The projections previously interpreted as hypapophyses are instead transverse processes. Openings on the anterior part of the skull, previously described as external nares, are probably internal nares. Important features not noted previously include accessory articulations on all presacral vertebrae, pachyostosis of dorsal vertebrae and ribs, and the presence of two pygal vertebrae. Phylogenetic analysis of 258 osteological characters and all the major squamate lineages suggests that Adriosaurus and dolichosaurs are successive sister-taxa to snakes. This is consistent with their long-bodied, limb-reduced morphology being intermediate between typical marine squamates (e.g., mosasaurs) and primitive marine snakes (pachyophiids). The analysis further reveals that up to five successive outgroups to living snakes (pachyophiids, Adriosaurus, dolichosaurs, Aphanizocnemus, and mosasauroids) are all marine, suggesting a marine (or at least, semi-aquatic) phase in snake origins. These phylogenetic results are robust whether multistate characters are ordered or unordered, thus refuting recent suggestions that snakes cluster with amphisbaenians and dibamids (rather than aquatic lizards) if multistate characters are left unordered. Also, the recent suggestion that Pachyrhachis shares synapomorphies with advanced snakes (macrostomatans) is shown to be poorly supported, because the reinterpretations of the relevant skull elements are unlikely and, even if accepted, the character states proposed to unite Pachyrhachis and advanced snakes are also present in more basal snakes and/or the nearest lizard outgroups, and are consequently primitive for snakes.

Historical Burden In Systematics And The Interrelationships Of ‘Parareptiles’

The interrelationships within the clade comprised of turtles, pareiasaurs, and procolophonid‐like taxa are investigated via a cladistic analysis incorporating 56 characters. A single most parsimonious tree was found (80 steps, c. i. = 0·8) in which the successive outgroups to turtles are: pareiasaurs, Sclerosaurus, lanthanosuchids, procolophonoids (=Owenetta, Barasaurus and procolophonids), and nyctiphruretians (= nycteroleterids). Thus, as suggested recently by other workers (Reisz, in Fischman, 1993) turtles are the highly modified survivors of a radiation of poorly‐known reptiles commonly called ‘parareptiles’. Pareiasaurs are united with turtles on the basis of twenty unambiguous derived features which are absent in other basal amniotes (=‘primitive reptiles’) and reptiliomorph amphibians: for example, the medially located choana, enlarged foramina palatinum posterius, blunt cultriform process, fully ossified medial wall of the prootic, opisthotic‐squamosal suture, lateral flange of exoccipital, loss of ventral cranial fissure, thickened braincase floor, ‘pleurosphenoid’ ossification, reduced presacral count, acromion process, trochanter major, reduced fifth pedal digit, and presence of transverse processes on most caudals. Recent phylogenetic proposals linking turtles with captorhinids, with dicynodonts, and with procolophonoids are evaluated. None of the proposed traits supporting the first two hypotheses is compelling. The procolophonoid hypotheses is supported by only one synapomorphy (the slender stapes). All other synapomorphies proposed in favour of the above groupings either occur in many other primitive amniotes, or are not primitive for turtles, or are not primitive for the proposed chelonian sister‐group.

TESTING THE RELATIONSHIP BETWEEN MORPHOLOGICAL AND MOLECULAR RATES OF CHANGE ALONG PHYLOGENIES

Abstract.— Molecular evolution has been considered to be essentially a stochastic process, little influenced by the pace of phenotypic change. This assumption was challenged by a study that demonstrated an association between rates of morphological and molecular change estimated for “total‐evidence” phylogenies, a finding that led some researchers to challenge molecular date estimates of major evolutionary radiations. Here we show that Omlands (1997) result is probably due to methodological bias, particularly phylogenetic nonindependence, rather than being indicative of an underlying evolutionary phenomenon. We apply three new methods specifically designed to overcome phylogenetic bias to 13 published phylogenetic datasets for vertebrate taxa, each of which includes both morphological characters and DNA sequence data. We find no evidence of an association between rates of molecular and morphological rates of change.

Molecular phylogeny and divergence dates for Australasian elapids and sea snakes (hydrophiinae): evidence from seven genes for rapid evolutionary radiations

One of the most prolific radiations of venomous snakes, the Australo‐Melanesian Hydrophiinae includes ∼100 species of Australasian terrestrial elapids plus all ∼60 species of viviparous sea snakes. Here, we estimate hydrophiine relationships based on a large data set comprising 5800 bp drawn from seven genes (mitochondrial: ND4, cytb, 12S, 16S; nuclear: rag1, cmos, myh). These data were analysed using parsimony, likelihood and Bayesian methods to better resolve hydrophiine phylogeny and provide a timescale for the terrestrial and marine radiations. Among oviparous forms, Cacophis, Furina and Demansia are basal to other Australian elapids (core oxyuranines). The Melanesian Toxicocalamus and Aspidomorphus group with Demansia, indicating multiple dispersal events between New Guinea and Australia. Oxyuranus and Pseudonaja form a robust clade. The small burrowing taxa form two separate clades, one consisting of Vermicella and Neelaps calanotus, and the other including Simoselaps, Brachyurophis and Neelaps bimaculatus. The viviparous terrestrial elapids form three separate groups: Acanthophis, the Rhinoplocephalus group and the Notechis–Hemiaspis group. True sea snakes (Hydrophiini) are robustly united with the Notechis–Hemiaspis group. Many of the retrieved groupings are consistent with previous molecular and morphological analyses, but the polyphyly of the viviparous and burrowing groups, and of Neelaps, are novel results. Bayesian relaxed clock analyses indicate very recent divergences: the ∼160 species of the core Australian radiation (including sea snakes) arose within the last 10 Myr, with most inter‐generic splits dating to between 10 and 6 Ma. The Hydrophis sea snake lineage is an exceptionally rapid radiation, with > 40 species evolving within the last 5 Myr.

Acute vision in the giant Cambrian predator Anomalocaris and the origin of compound eyes

Until recently, intricate details of the optical design of non-biomineralized arthropod eyes remained elusive in Cambrian Burgess-Shale-type deposits, despite exceptional preservation of soft-part anatomy in such Konservat-Lagerstätten. The structure and development of ommatidia in arthropod compound eyes support a single origin some time before the latest common ancestor of crown-group arthropods, but the appearance of compound eyes in the arthropod stem group has been poorly constrained in the absence of adequate fossils. Here we report 2–3-cm paired eyes from the early Cambrian (approximately 515 million years old) Emu Bay Shale of South Australia, assigned to the Cambrian apex predator Anomalocaris. Their preserved visual surfaces are composed of at least 16,000 hexagonally packed ommatidial lenses (in a single eye), rivalling the most acute compound eyes in modern arthropods. The specimens show two distinct taphonomic modes, preserved as iron oxide (after pyrite) and calcium phosphate, demonstrating that disparate styles of early diagenetic mineralization can replicate the same type of extracellular tissue (that is, cuticle) within a single Burgess-Shale-type deposit. These fossils also provide compelling evidence for the arthropod affinities of anomalocaridids, push the origin of compound eyes deeper down the arthropod stem lineage, and indicate that the compound eye evolved before such features as a hardened exoskeleton. The inferred acuity of the anomalocaridid eye is consistent with other evidence that these animals were highly mobile visual predators in the water column. The existence of large, macrophagous nektonic predators possessing sharp vision—such as Anomalocaris—within the early Cambrian ecosystem probably helped to accelerate the escalatory ‘arms race’ that began over half a billion years ago.

Molecules, morphology, and ecology indicate a recent, amphibious ancestry for echidnas

The semiaquatic platypus and terrestrial echidnas (spiny anteaters) are the only living egg-laying mammals (monotremes). The fossil record has provided few clues as to their origins and the evolution of their ecological specializations; however, recent reassignment of the Early Cretaceous Teinolophos and Steropodon to the platypus lineage implies that platypuses and echidnas diverged >112.5 million years ago, reinforcing the notion of monotremes as living fossils. This placement is based primarily on characters related to a single feature, the enlarged mandibular canal, which supplies blood vessels and dense electrosensory receptors to the platypus bill. Our reevaluation of the morphological data instead groups platypus and echidnas to the exclusion of Teinolophos and Steropodon and suggests that an enlarged mandibular canal is ancestral for monotremes (partly reversed in echidnas, in association with general mandibular reduction). A multigene evaluation of the echidna–platypus divergence using both a relaxed molecular clock and direct fossil calibrations reveals a recent split of 19–48 million years ago. Platypus-like monotremes (Monotrematum) predate this divergence, indicating that echidnas had aquatically foraging ancestors that reinvaded terrestrial ecosystems. This ecological shift and the associated radiation of echidnas represent a recent expansion of niche space despite potential competition from marsupials. Monotremes might have survived the invasion of marsupials into Australasia by exploiting ecological niches in which marsupials are restricted by their reproductive mode. Morphology, ecology, and molecular biology together indicate that Teinolophos and Steropodon are basal monotremes rather than platypus relatives, and that living monotremes are a relatively recent radiation.

Unalignable sequences and molecular evolution

Abstract Alignment ambiguity is a widespread problem in molecular evolutionary studies that has received insufficient attention. Most studies ignore such regions by deleting them before analyses, even though alignment-ambiguous regions can contain useful phylogenetic and evolutionary information. The alignment ambiguity might affect only one taxon, the region being readily alignable and phylogenetically informative across all other taxa. Alternatively, all possible alignments can consistently imply certain relationships. Because they are usually the most rapidly evolving regions, alignment-ambiguous regions might be those that are most able to resolve closely spaced divergences and contribute to estimates of branch lengths, evolutionary rates and divergence times. Three methods to incorporate such regions into phylogenetic and evolutionary analyses have been devised. The multiple analysis method evaluates each plausible alignment separately and seeks areas of congruence among the resultant trees, whereas the elision method combines all plausible alignments into a single analysis. Fragment-level alignment (= fixed states, INAASE) treats the entire unalignable section as a single but highly complex multistate character. Although these methods still need refining, they are preferable to discarding large portions of hard-earned and potentially informative sequence data.

Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds

A limiting factor in many molecular dating studies is shortage of reliable calibrations. Current methods for choosing calibrations (e.g. cross-validation) treat them as either correct or incorrect, whereas calibrations probably lie on a continuum from highly accurate to very poor. Bayesian relaxed clock analysis permits inclusion of numerous candidate calibrations as priors: provided most calibrations are reliable, the model appropriate and the data informative, the accuracy of each calibration prior can be evaluated. If a calibration is accurate, then the analysis will support the prior so that the posterior estimate reflects the prior; if a calibration is poor, the posterior will be forced away from the prior. We use this approach to test two fossil dates recently proposed as standard calibrations within vertebrates. The proposed bird–crocodile calibration (approx. 247 Myr ago) appears to be accurate, but the proposed bird–lizard calibration (approx. 255 Myr ago) is substantially too recent.