Fossils, Molecules, Divergence Times, and the Origin of Lissamphibians

Evolutionary Insights into the Relationship of Frogs, Salamanders, and Caecilians and Their Adaptive Traits, with an Emphasis on Salamander Regeneration and Longevity

The extant amphibians have developed uncanny abilities to adapt to their environment. I compared the genes of amphibians to those of other vertebrates to investigate the genetic changes underlying their unique traits, especially salamanders' regeneration and longevity. Using the well-supported Batrachia tree, I found that salamander genomes have undergone accelerated adaptive evolution, especially for development-related genes. The group-based comparison showed that several genes are under positive selection, rapid evolution, and unexpected parallel evolution with traits shared by distantly related species, such as the tail-regenerative lizard and the longer-lived naked mole rat. The genes, such as EEF1E1, PAFAH1B1, and OGFR, may be involved in salamander regeneration, as they are involved in the apoptotic process, blastema formation, and cell proliferation, respectively. The genes PCNA and SIRT1 may be involved in extending lifespan, as they are involved in DNA repair and histone modification, respectively. Some genes, such as PCNA and OGFR, have dual roles in regeneration and aging, which suggests that these two processes are interconnected. My experiment validated the time course differential expression pattern of SERPINI1 and OGFR, two genes that have evolved in parallel in salamanders and lizards during the regeneration process of salamander limbs. In addition, I found several candidate genes responsible for frogs' frequent vocalization and caecilians' degenerative vision. This study provides much-needed insights into the processes of regeneration and aging, and the discovery of the critical genes paves the way for further functional analysis, which could open up new avenues for exploiting the genetic potential of humans and improving human well-being.

Spatio-Temporal Evolutionary Patterns of the Pieridae Butterflies (Lepidoptera: Papilionoidea) Inferred from Mitogenomic Data

Pieridae is one of the largest and almost cosmopolitan groups of butterflies, which plays an important role in natural ecosystems; however, to date, its phylogeny and evolutionary history have not been fully resolved. In this study, we obtained the complete or nearly complete mitochondrial genomes of 100 pierid taxa (six newly sequenced, sixty extracted from the whole-genome data, and thirty-four directly available from GenBank). At the same time, for the first time, we conducted comparative mitogenomic and phylogenetic analyses based on these mitogenomic data, to further clarify their spatio-temporal evolutionary patterns. Comparative mitogenomic analysis showed that, except for cox2, the GC content of each of the 13 protein-coding genes (PCGs) in the rapidly diverging subfamily Pierinae was higher than in its sister group Coliadinae. Moreover, the dN/dS values of nine genes (atp6, atp8, cox1, cox3, cob, nad1, nad3, nad5, and nad6) in Pierinae were also relatively higher than those in its sister group, Coliadinae. Phylogenetic analysis showed that all the resultant phylogenetic trees were generally in agreement with those of previous studies. The Pierinae family contained six clades in total with the relationship of (Leptosiaini + (((Nepheroniini + Arthocharidini) + Teracolini) + (Pierini + Elodini))). The Pieridae originated in the Palearctic region approximately 72.3 million years ago in the late Cretaceous, and the subfamily Pierinae diverged from this family around 57.9 million years ago in the Oriental region, shortly after the K-Pg mass extinction event; in addition, the spatio-temporal evolutionary patterns of Pierinae were closely correlated with geological events and environmental changes, as well as the host plant coevolutionary scenario in Earth's history. However, some incongruencies were observed between our results and those of previous studies in terms of shallow phylogenies for a few taxa, and should be further investigated.

Modulation of the cholinergic pathway induced by skin secretion of Phyllomedusa iheringii Boulenger, 1885 in a vertebrate model

Amphibians represent one of the main natural sources of bioactive molecules of interest to biotechnological research. The Phyllomedusidae family has several species occurring in Brazil and some studies demonstrate the biological potential of poisons of these species, however many still need to be characterized. Phyllomedusa iheringii is endemic in Brazilian and Uruguayan Pampa Biome and has little data in the literature regarding the action of its poison on experimental organisms. Thus, the present work evaluates the biological activity of P. iheringii secretion on the central and peripheral nervous system of a vertebrate model. The skin secretions of P. iheringii (SSPI) were collected through manual compression and electrical stimulation of the animal's bodies. The resulting content was used in neurobiological tests searching for modulatory effects on the main pathways involved in the neurotoxicity mechanism of vertebrates. SSPI affected the contraction force of the chick biventer cervicis muscle (Gallus gallus domesticus) at some concentrations used (5, 10, and 12 μg/mL). In slices from the cerebral cortex of G. gallus domesticus an increase in cell viability was observed after treatment with SSPI (10 μg/mL) and a neuroprotective effect when treated simultaneously with hydrogen peroxide (H₂O₂), Neostigmine (NEO) and Trichlorfon (TRI). The cholinergic pathway is possibly the main pathway modulated by SSPI since assays with the cerebral cortex and biventer cervicis muscle demonstrated the increased activity of the enzyme acetylcholinesterase (AChE) (SSPI 10 μg/mL and 12 μg/mL, respectively). SSPI (10 μg/mL) also prevented the modulation of NEO and TRI, two recognized anticholinesterase agents, in AChE activity in slices of the cerebral cortex. Therefore, our results have demonstrated the unpublished biotechnological potential of P. iheringii over the vertebrate model and its modulation on the nervous system, with apparent action on the cholinergic pathway.

Pseudogenized Amelogenin Reveals Early Tooth Loss in True Toads (Anura: Bufonidae)

Extant anurans (frogs and toads) exhibit reduced dentition, ranging from a lack of mandibular teeth to complete edentulation, as observed in the true toads of the family Bufonidae. The evolutionary time line of these reductions remains vague due to a poor fossil record. Previous studies have demonstrated an association between the lack of teeth in edentulous vertebrates and the pseudogenization of the major tooth enamel gene amelogenin (AMEL) through accumulation of deleterious mutations and the disruption of its coding sequence. In this study, we have harnessed the pseudogenization of AMEL as a molecular dating tool to correlate loss of dentition with genomic mutation patterns during the rise of the family Bufonidae. Specifically, we have utilized AMEL pseudogenes in three members of the family as a tool to estimate the putative date of edentulation in true toads. Comparison of AMEL sequences from Rhinella marina, Bufo gargarizans and Bufo bufo, with nine extant, dentulous frogs, revealed mutations confirming AMEL inactivation in Bufonidae. AMEL pseudogenes in modern bufonids also exhibited remarkably high 86-93% sequence identity among each other, with only a slight increase in substitution rate and relaxation of selective pressure, in comparison with functional copies in other anurans. Moreover, using selection intensity estimates and synonymous substitution rates, analysis of functional and pseudogenized AMEL resulted in an estimated inactivation window of 46-60 million years ago in the lineage leading to modern true toads, a time line that coincides with the rise of the family Bufonidae.

First assessment of genetic diversity, population structure and historical population dynamics of Myocastor coypus (Rodentia: Echimyidae) in the centre of its native range

Myocastor coypus is a rodent native to South America that is strongly linked to river systems. Past studies on the coypu in its native distribution range were aimed mainly at its ecology and parasitology and genetic studies are very limited. Here, we used sequences of the mitochondrial D-loop control region to study the genetic diversity, population genetic structure and some aspects of historical population dynamics of coypu at the centre of its native range. Our results showed moderate to high levels of genetic diversity and an absence of genetic structure in the study area. Bayesian analysis of population structure (BAPS) showed the existence of two haplogroups distributed in most sampling sites. These results suggest that movement of coypu is independent of the structure of current known river networks and its dispersal patterns are related to flooding events and the presence of lagoons connecting river basins. The demographic expansion patterns in these populations and those of other South American mammals during the Late Pleistocene support the hypothesis that demographic changes in wild populations are related to fluctuations in climate and ecology. The mitochondrial data obtained in this study constitute the first record of these types of sequences in the native range of M. coypus.

Myocastor coypus es un roedor nativo de Sudamérica fuertemente ligado a los sistemas fluviales. En su área de distribución nativa, los estudios llevados a cabo en la especie estuvieron dirigidos principalmente a estudiar su ecología y zoonosis, sin embargo, los estudios genéticos son muy limitados. En el presente trabajo, utilizamos secuencias del d-loop de la Región Control del ADN mitocondrial para estudiar la diversidad genética, la estructura genética de la población y algunos aspectos de la dinámica poblacional histórica del coipo en el centro de su área de distribución. Nuestros resultados develaron niveles moderados y altos de diversidad genética, y ausencia de estructuración genética en el área de estudio. El BAPS exhibió la existencia de dos haplogrupos distribuidos en la mayoría de los sitios de muestreo. Estos resultados sugieren que el movimiento de los coipos sería independiente de la estructura de las redes fluviales actuales conocidas y que sus patrones de dispersión estarían relacionados con eventos de inundación y con la presencia de lagunas que conectan las cuencas fluviales. Los patrones de expansión demográfica observados en esta población de coipos y en otras poblaciones de mamíferos sudamericanos durante el Pleistoceno tardío, apoyan la hipótesis de que las fluctuaciones climáticas y los cambios ecológicos están relacionados con cambios demográficos en las poblaciones silvestres. Los datos mitocondriales obtenidos en este estudio constituyen el primer registro de este tipo de secuencias en el área de distribución nativa de esta especie.

Cooperation behavior of fore- And hindlimbs during jumping in Rana dybowskii and Xenopus laevis

Frogs are characterized by their outstanding jumping ability, depending on the rapid extension of hindlimbs to propel their bodies into air. A typical jumping cycle could be broken into four phases: preparation, takeoff, flight, and landing. Considerable research has been performed to discuss the function of hindlimbs of frogs during takeoff phase, whereas the literature of limbs' motion in jumping between different species was limited. To profile the evolution of locomotion in anurans, it is necessary to investigate on the motion of fore- and hindlimbs of frogs within different taxa. In this work, we put forward a detailed description of jumping behavior of two frog species, Rana dybowskii and Xenopus laevis. High-speed cameras were used to explore the movement of different joints in fore- and hindlimbs of these two animals, and kinematic analysis was operated to identify both homologous behaviors and significant differences between them. We found that the Rana dybowskii's fore- and hindlimbs had good cooperation during jumping, while the Xenopus laevis' uncooperative behavior in limbs may give a functional explanation for the deficiency in terrestrial jumping; besides, the R. dybowskii's landing followed the "hands-belly-feet slap" strategy, and Xenopus laevis had clumsy landing with "belly-flops" sequence. The result gained here clarifies the cooperation behavior of anuran limbs and may supply a new insight into our understanding of the anuran's evolution.

The Making of Calibration Sausage Exemplified by Recalibrating the Transcriptomic Timetree of Jawed Vertebrates

Molecular divergence dating has the potential to overcome the incompleteness of the fossil record in inferring when cladogenetic events (splits, divergences) happened, but needs to be calibrated by the fossil record. Ideally but unrealistically, this would require practitioners to be specialists in molecular evolution, in the phylogeny and the fossil record of all sampled taxa, and in the chronostratigraphy of the sites the fossils were found in. Paleontologists have therefore tried to help by publishing compendia of recommended calibrations, and molecular biologists unfamiliar with the fossil record have made heavy use of such works (in addition to using scattered primary sources and copying from each other). Using a recent example of a large node-dated timetree inferred from molecular data, I reevaluate all 30 calibrations in detail, present the current state of knowledge on them with its various uncertainties, rerun the dating analysis, and conclude that calibration dates cannot be taken from published compendia or other secondary or tertiary sources without risking strong distortions to the results, because all such sources become outdated faster than they are published: 50 of the (primary) sources I cite to constrain calibrations were published in 2019, half of the total of 280 after mid-2016, and 90% after mid-2005. It follows that the present work cannot serve as such a compendium either; in the slightly longer term, it can only highlight known and overlooked problems. Future authors will need to solve each of these problems anew through a thorough search of the primary paleobiological and chronostratigraphic literature on each calibration date every time they infer a new timetree, and that literature is not optimized for that task, but largely has other objectives.

Brain Reconstruction Across the Fish-Tetrapod Transition; Insights From Modern Amphibians

The fish-tetrapod transition (which incorporates the related fin-limb and water-land transitions) is celebrated as one of the most important junctions in vertebrate evolution. Sarcopterygian fishes (the “lobe-fins”) are today represented by lungfishes and coelacanths, but during the Paleozoic they were much more diverse. It was some of these sarcopterygians, a lineage of the tetrapodomorph fishes, that gave rise to tetrapods (terrestrial vertebrates with limbs bearing digits). This spectacular leap took place during the Devonian Period. Due to the nature of preservation, it is the hard parts of an animal’s body that are most likely to fossilize, while soft tissues such as muscular and brain tissues, typically fail to do so. Thus, our understanding of the adaptations of the hard skeletal structures of vertebrates is considerably greater than that of the soft tissue systems. Fortunately, the braincases of early vertebrates are often ossified and thereby have the potential to provide detailed morphological information. However, the correspondence between brain and endocast (an internal mold of the cavity) has historically been considered poor in most “lower” vertebrates and consequently neglected in such studies of brain evolution. Despite this, recent work documenting the spatial relationship in extant basal sarcopterygians (coelacanth, lungfish, axolotl, and salamander) has highlighted that this is not uniformly the case. Herein, we quantify and illustrate the brain-endocast relationship in four additional extant basal tetrapod exemplars: neobatrachian anurans (frogs) Breviceps poweri and Ceratophrys ornata; and gymnophionans (caecilians) Gegeneophis ramaswamii and Rhinatrema bivittatum. We show that anurans and caecilians appear to have brains that fill their endocasts to a similar degree to that of lungfishes and salamanders, but not coelacanth. Ceratophrys has considerably lower correspondence between the brain and endocast in the olfactory tract and mesencephalic regions, while Breviceps has low correspondence along its ventral endocranial margin. The brains of caecilians reflect their endocasts most closely (vol. ∼70%). The telencephalon is tightly fitted within the endocast in all four taxa. Our findings highlight the need to adequately assess the brain-endocast relationship in a broad range of vertebrates, in order to inform neural reconstructions of fossil taxa using the Extant Phylogenetic Bracket approach and future studies of brain evolution.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

A review of the fossil record of caecilians (Lissamphibia: Gymnophionomorpha) with comments on its use to calibrate molecular timetrees

Gymnophiona, popularly known as caecilians, the most poorly known major taxon of extant amphibians, are elongate and limbless tetrapods, with compact ossified skulls and reduced eyes, mainly adapted to fossorial life as adults. Caecilians are poorly represented in the fossil record, but despite the scarcity of fossil specimens described (only four named taxa, in addition to indeterminate fragmentary material), their fossils play a key role in our knowledge of the origin and evolution of Lissamphibia, as well as contribute directly to a better understanding of the phylogeny, taxonomy and biogeography of extant gymnophionan taxa. These records are scattered throughout geological time (from the Jurassic to the sub-Recent) and space (North and South America and Africa). Here, we revisit the caecilian fossil record, providing a brief description of all known extinct taxa described so far, along with general remarks about their impact on systematics, time range, and geographical distribution of the clade, as well as prospects for future research. Possible calibration constraints based on the caecilian fossil record are provided.

Can We Reliably Calibrate Deep Nodes in the Tetrapod Tree? Case Studies in Deep Tetrapod Divergences

Recent efforts have led to the development of extremely sophisticated methods for incorporating tree-wide data and accommodating uncertainty when estimating the temporal patterns of phylogenetic trees, but assignment of prior constraints on node age remains the most important factor. This depends largely on understanding substantive disagreements between specialists (paleontologists, geologists, and comparative anatomists), which are often opaque to phylogeneticists and molecular biologists who rely on these data as downstream users. This often leads to misunderstandings of how the uncertainty associated with node age minima arises, leading to inappropriate treatments of that uncertainty by phylogeneticists. In order to promote dialogue on this subject, we here review factors (phylogeny, preservational megabiases, spatial and temporal patterns in the tetrapod fossil record) that complicate assignment of prior node age constraints for deep divergences in the tetrapod tree, focusing on the origin of crown-group Amniota, crown-group Amphibia, and crown-group Tetrapoda. We find that node priors for amphibians and tetrapods show high phylogenetic lability and different phylogenetic treatments identifying disparate taxa as the earliest representatives of these crown groups. This corresponds partially to the well-known problem of lissamphibian origins but increasingly reflects deeper instabilities in early tetrapod phylogeny. Conversely, differences in phylogenetic treatment do not affect our ability to recognize the earliest crown-group amniotes but do affect how diverse we understand the earliest amniote faunas to be. Preservational megabiases and spatiotemporal heterogeneity of the early tetrapod fossil record present unrecognized challenges in reliably estimating the ages of tetrapod nodes; the tetrapod record throughout the relevant interval is spatially restricted and disrupted by several major intervals of minimal sampling coincident with the emergence of all three crown groups. Going forward, researchers attempting to calibrate the ages for these nodes, and other similar deep nodes in the metazoan fossil record, should consciously consider major phylogenetic uncertainty, preservational megabias, and spatiotemporal heterogeneity, preferably examining the impact of working hypotheses from multiple research groups. We emphasize a need for major tetrapod collection effort outside of classic European and North American sections, particularly from the southern hemisphere, and suggest that such sampling may dramatically change our timelines of tetrapod evolution.

The lacrimal/ectethmoid region of waterfowl (Aves, Anseriformes): Phylogenetic signal and major evolutionary patterns

Waterfowl (Aves, Anseriformes) constitute an ancient global radiation, and understanding the pattern and timing of their evolution requires a well-corroborated phylogeny including extant species and fossils. Following the molecular advances in avian systematics, however, morphology has often been held as misleading, yet congruence with molecular data has been shown to vary considerably among different skeletal parts. Here, we explore phylogenetic signal in discrete characters of the lacrimal/ectethmoid region of waterfowl, which is highly variable among species and constitutes a rich source of data. We do so by combining cladistic and multivariate approaches, and using phylogenetic comparative methods. We quantitatively recognize three major morphological types among lacrimal bones, and discuss homoplasy and potential synapomorphies of major clades using a molecular backbone tree. Our results clearly indicate that the lacrimal bone carries substantial phylogenetic signal and could be of systematic value at different levels of the phylogeny of waterfowl, feeding the exploration of other regions of the skull with this combined approach.

Rates and Rocks: Strengths and Weaknesses of Molecular Dating Methods

I present here an in-depth, although non-exhaustive, review of two topics in molecular dating. Clock models, which describe the evolution of the rate of evolution, are considered first. Some of the shortcomings of popular approaches-uncorrelated clock models in particular-are presented and discussed. Autocorrelated models are shown to be more reasonable from a biological perspective. Some of the most recent autocorrelated models also rely on a coherent treatment of instantaneous and average substitution rates while previous models are based on implicit approximations. Second, I provide a brief overview of the processes involved in collecting and preparing fossil data. I then review the main techniques that use this data for calibrating the molecular clock. I argue that, in its current form, the fossilized birth-death process relies on assumptions about the mechanisms underlying fossilization and the data collection process that may negatively impact the date estimates. Node-dating approaches make better use of the data available, even though they rest on paleontologists' intervention to prepare raw fossil data. Altogether, this study provides indications that may help practitioners in selecting appropriate methods for molecular dating. It will also hopefully participate in defining the contour of future methodological developments in the field.

Quantifying the Error of Secondary vs. Distant Primary Calibrations in a Simulated Environment

Using calibrations to obtain absolute divergence times is standard practice in molecular clock studies. While the use of primary (e.g., fossil) calibrations is preferred, this approach can be limiting because of their rarity in fast-growing datasets. Thus, alternatives need to be explored, such as the use of secondary (molecularly-derived) calibrations that can anchor a timetree in a larger number of nodes. However, the use of secondary calibrations has been discouraged in the past because of concerns in the error rates of the node estimates they produce with an apparent high precision. Here, we quantify the amount of errors in estimates produced by the use of secondary calibrations relative to true times and primary calibrations placed on distant nodes. We find that, overall, the inaccuracies in estimates based on secondary calibrations are predictable and mirror errors associated with primary calibrations and their confidence intervals. Additionally, we find comparable error rates in estimated times from secondary calibrations and distant primary calibrations, although the precision of estimates derived from distant primary calibrations is roughly twice as good as that of estimates derived from secondary calibrations. This suggests that increasing dataset size to include primary calibrations may produce divergence times that are about as accurate as those from secondary calibrations, albeit with a higher precision. Overall, our results suggest that secondary calibrations may be useful to explore the parameter space of plausible evolutionary scenarios when compared to time estimates obtained with distant primary calibrations.

Exact Distribution of Divergence Times from Fossil Ages and Tree Topologies

Being given a phylogenetic tree of both extant and extinct taxa in which the fossil ages are the only temporal information (namely, in which divergence times are considered unknown), we provide a method to compute the exact probability distribution of any divergence time of the tree with regard to any speciation (cladogenesis), extinction, and fossilization rates under the Fossilized Birth–Death model. We use this new method to obtain a probability distribution for the age of Amniota (the synapsid/sauropsid or bird/mammal divergence), one of the most-frequently used dating constraints. Our results suggest an older age (between about 322 and 340 Ma) than has been assumed by most studies that have used this constraint (which typically assumed a best estimate around 310–315 Ma) and provide, for the first time, a method to compute the shape of the probability density for this divergence time. [Divergence times; fossil ages; fossilized birth–death model; probability distribution.]

Selective Sampling of Species and Fossils Influences Age Estimates Under the Fossilized Birth-Death Model

The fossilized birth-death (FBD) model allows the estimation of species divergence times from molecular and fossil information in a coherent framework of diversification and fossil sampling. Some assumptions of the FBD model, however, are difficult to meet in phylogenetic analyses of highly diverse groups. Here, I use simulations to assess the impact of extreme model violations, including diversified sampling of species and the exclusive use of the oldest fossils per clade, on divergence times estimated with the FBD model. My results demonstrate that selective sampling of fossils can produce dramatically overestimated divergence times when the FBD model is used for inference, due to an interplay of underestimates for the model parameters net diversification rate, turnover, and fossil-sampling proportion. In contrast, divergence times estimated with CladeAge, a method that uses information about the oldest fossils per clade together with estimates of sampling and diversification rates, are accurate under these conditions. Practitioners of Bayesian divergence-time estimation should therefore ensure that the dataset conforms to the expectations of the FBD model, or estimates of sampling and diversification rates should be obtained a priori so that CladeAge can be used for the inference.

Historical museum collections clarify the evolutionary history of cryptic species radiation in the world's largest amphibians

Inaccurate taxonomic assessment of threatened populations can hinder conservation prioritization and management, with human-mediated population movements obscuring biogeographic patterns and confounding reconstructions of evolutionary history. Giant salamanders were formerly distributed widely across China, and are interpreted as a single species, Andrias davidianus. Previous phylogenetic studies have identified distinct Chinese giant salamander lineages but were unable to associate these consistently with different landscapes, probably because population structure has been modified by human-mediated translocations for recent commercial farming. We investigated the evolutionary history and relationships of allopatric Chinese giant salamander populations with Next-Generation Sequencing methods, using historical museum specimens and late 20th-century samples, and retrieved partial or near-complete mitogenomes for 17 individuals. Samples from populations unlikely to have been affected by translocations form three clades from separate regions of China, spatially congruent with isolation by either major river drainages or mountain ranges. Pliocene-Pleistocene divergences for these clades are consistent with topographic modification of southern China associated with uplift of the Qinghai-Tibet Plateau. General Mixed Yule Coalescent model analysis indicates that these clades represent separate species: Andrias davidianus (Blanchard, 1871) (northern Yangtze/Sichuan), Andrias sligoi (Boulenger, 1924) (Pearl/Nanling), and an undescribed species (Huangshan). Andrias sligoi is possibly the world's largest amphibian. Inclusion of additional reportedly wild samples from areas of known giant salamander exploitation and movement leads to increasing loss of biogeographic signal. Wild Chinese giant salamander populations are now critically depleted or extirpated, and conservation actions should be updated to recognize the existence of multiple species.

Phylogeny of Paleozoic limbed vertebrates reassessed through revision and expansion of the largest published relevant data matrix

The largest published phylogenetic analysis of early limbed vertebrates (Ruta M, Coates MI. 2007. Journal of Systematic Palaeontology 5:69-122) recovered, for example, Seymouriamorpha, Diadectomorpha and (in some trees) Caudata as paraphyletic and found the "temnospondyl hypothesis" on the origin of Lissamphibia (TH) to be more parsimonious than the "lepospondyl hypothesis" (LH)-though only, as we show, by one step. We report 4,200 misscored cells, over half of them due to typographic and similar accidental errors. Further, some characters were duplicated; some had only one described state; for one, most taxa were scored after presumed relatives. Even potentially continuous characters were unordered, the effects of ontogeny were not sufficiently taken into account, and data published after 2001 were mostly excluded. After these issues are improved-we document and justify all changes to the matrix-but no characters are added, we find (Analysis R1) much longer trees with, for example, monophyletic Caudata, Diadectomorpha and (in some trees) Seymouriamorpha; Ichthyostega either crownward or rootward of Acanthostega; and Anthracosauria either crownward or rootward of Temnospondyli. The LH is nine steps shorter than the TH (R2; constrained) and 12 steps shorter than the "polyphyly hypothesis" (PH-R3; constrained). Brachydectes (Lysorophia) is not found next to Lissamphibia; instead, a large clade that includes the adelogyrinids, urocordylid "nectrideans" and aïstopods occupies that position. As expected from the taxon/character ratio, most bootstrap values are low. Adding 56 terminal taxa to the original 102 increases the resolution (and decreases most bootstrap values). The added taxa range in completeness from complete articulated skeletons to an incomplete lower jaw. Even though the lissamphibian-like temnospondyls Gerobatrachus, Micropholis and Tungussogyrinus and the extremely peramorphic salamander Chelotriton are added, the difference between LH (R4; unconstrained) and TH (R5) rises to 10 steps, that between LH and PH (R6) to 15; the TH also requires several more regains of lost bones than the LH. Casineria, in which we tentatively identify a postbranchial lamina, emerges rather far from amniote origins in a gephyrostegid-chroniosuchian grade. Bayesian inference (Analysis EB, settings as in R4) mostly agrees with R4. High posterior probabilities are found for Lissamphibia (1.00) and the LH (0.92); however, many branches remain weakly supported, and most are short, as expected from the small character sample. We discuss phylogeny, approaches to coding, methods of phylogenetics (Bayesian inference vs. equally weighted vs. reweighted parsimony), some character complexes (e.g. preaxial/postaxial polarity in limb development), and prospects for further improvement of this matrix. Even in its revised state, the matrix cannot provide a robust assessment of the phylogeny of early limbed vertebrates. Sufficient improvement will be laborious-but not difficult.

Effects of amphibian phylogeny, climate and human impact on the occurrence of the amphibian-killing chytrid fungus

Chytridiomycosis, due to the fungus Batrachochytrium dendrobatidis (Bd), has been associated with the alarming decline and extinction crisis of amphibians worldwide. Because conservation programs are implemented locally, it is essential to understand how the complex interactions among host species, climate and human activities contribute to Bd occurrence at regional scales. Using weighted phylogenetic regressions and model selection, we investigated geographic patterns of Bd occurrence along a latitudinal gradient of 1500 km within a biodiversity hot spot in Chile (1845 individuals sampled from 253 sites and representing 24 species), and its association with climatic, socio-demographic and economic variables. Analyses show that Bd prevalence decreases with latitude although it has increased by almost 10% between 2008 and 2013, possibly reflecting an ongoing spread of Bd following the introduction of Xenopus laevis. Occurrence of Bd was higher in regions with high gross domestic product (particularly near developed centers) and with a high variability in rainfall regimes, whereas models including other bioclimatic or geographic variables, including temperature, exhibited substantially lower fit and virtually no support based on Akaike weights. In addition, Bd prevalence exhibited a strong phylogenetic signal, with five species having high numbers of infected individuals and higher prevalence than the average of 13.3% across all species. Taken together, our results highlight that Bd in Chile might still be spreading south, facilitated by a subset of species that seem to play an important epidemiological role maintaining this pathogen in the communities, in combination with climatic and human factors affecting the availability and quality of amphibian breeding sites. This information may be employed to design conservation strategies and mitigate the impacts of Bd in the biodiversity hot spot of southern Chile, and similar studies may prove useful to disentangle the role of different factors contributing to the emergence and spread of this catastrophic disease.

Dated phylogeny and dispersal history of the butterfly subfamily Nymphalinae (Lepidoptera: Nymphalidae)

The origin and dispersal history of the large butterfly subfamily Nymphalinae are not fully understood, due to internal phylogenetic and time calibration issues. We conducted phylogenetic and dating analyses using mitochondrial and nuclear genes of biogeographically diverse groups of the Nymphalinae in order to resolve some controversial relationships and the paleobiogeographic pattern of the subfamily. Our results support the sister relationship of Vanessa (Tribe Nymphalini) and the Nymphalis-group, and the grouping of the three old-world genera (Rhinopalpa, Kallimoides and Vanessula) within Tribe Victorinini. Molecular dating analyses invoking two additional calibrations under the butterfly-host plant coevolutionary scenarios result in a relatively deeper divergence of the subfamily's two major clades (Nymphalini and the Kallimoids), compatible with the Cretaceous floral turnover scenario during the so-called Cretaceous Terrestrial Revolution. Phylobiogeographic analyses reveal that the Oriental region is probably the center of early divergences for Nymphalinae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by repeated dispersals into the rest of the Old World and the New World during various periods beginning in Eocene. The biogeographic history indicates that temperature changes and host-plant diversification may have facilitated the dispersals of this butterfly subfamily, with accelerated global colonization during the middle to late Miocene.

Stem caecilian from the Triassic of Colorado sheds light on the origins of Lissamphibia

The origin of the limbless caecilians remains a lasting question in vertebrate evolution. Molecular phylogenies and morphology support that caecilians are the sister taxon of batrachians (frogs and salamanders), from which they diverged no later than the early Permian. Although recent efforts have discovered new, early members of the batrachian lineage, the record of pre-Cretaceous caecilians is limited to a single species, Eocaecilia micropodia The position of Eocaecilia within tetrapod phylogeny is controversial, as it already acquired the specialized morphology that characterizes modern caecilians by the Jurassic. Here, we report on a small amphibian from the Upper Triassic of Colorado, United States, with a mélange of caecilian synapomorphies and general lissamphibian plesiomorphies. We evaluated its relationships by designing an inclusive phylogenetic analysis that broadly incorporates definitive members of the modern lissamphibian orders and a diversity of extinct temnospondyl amphibians, including stereospondyls. Our results place the taxon confidently within lissamphibians but demonstrate that the diversity of Permian and Triassic stereospondyls also falls within this group. This hypothesis of caecilian origins closes a substantial morphologic and temporal gap and explains the appeal of morphology-based polyphyly hypotheses for the origins of Lissamphibia while reconciling molecular support for the group's monophyly. Stem caecilian morphology reveals a previously unrecognized stepwise acquisition of typical caecilian cranial apomorphies during the Triassic. A major implication is that many Paleozoic total group lissamphibians (i.e., higher temnospondyls, including the stereospondyl subclade) fall within crown Lissamphibia, which must have originated before 315 million years ago.

Biotic and environmental dynamics through the Late Jurassic-Early Cretaceous transition: evidence for protracted faunal and ecological turnover

The Late Jurassic to Early Cretaceous interval represents a time of environmental upheaval and cataclysmic events, combined with disruptions to terrestrial and marine ecosystems. Historically, the Jurassic/Cretaceous (J/K) boundary was classified as one of eight mass extinctions. However, more recent research has largely overturned this view, revealing a much more complex pattern of biotic and abiotic dynamics than has previously been appreciated. Here, we present a synthesis of our current knowledge of Late Jurassic-Early Cretaceous events, focusing particularly on events closest to the J/K boundary. We find evidence for a combination of short-term catastrophic events, large-scale tectonic processes and environmental perturbations, and major clade interactions that led to a seemingly dramatic faunal and ecological turnover in both the marine and terrestrial realms. This is coupled with a great reduction in global biodiversity which might in part be explained by poor sampling. Very few groups appear to have been entirely resilient to this J/K boundary 'event', which hints at a 'cascade model' of ecosystem changes driving faunal dynamics. Within terrestrial ecosystems, larger, more-specialised organisms, such as saurischian dinosaurs, appear to have suffered the most. Medium-sized tetanuran theropods declined, and were replaced by larger-bodied groups, and basal eusauropods were replaced by neosauropod faunas. The ascent of paravian theropods is emphasised by escalated competition with contemporary pterosaur groups, culminating in the explosive radiation of birds, although the timing of this is obfuscated by biases in sampling. Smaller, more ecologically diverse terrestrial non-archosaurs, such as lissamphibians and mammaliaforms, were comparatively resilient to extinctions, instead documenting the origination of many extant groups around the J/K boundary. In the marine realm, extinctions were focused on low-latitude, shallow marine shelf-dwelling faunas, corresponding to a significant eustatic sea-level fall in the latest Jurassic. More mobile and ecologically plastic marine groups, such as ichthyosaurs, survived the boundary relatively unscathed. High rates of extinction and turnover in other macropredaceous marine groups, including plesiosaurs, are accompanied by the origin of most major lineages of extant sharks. Groups which occupied both marine and terrestrial ecosystems, including crocodylomorphs, document a selective extinction in shallow marine forms, whereas turtles appear to have diversified. These patterns suggest that different extinction selectivity and ecological processes were operating between marine and terrestrial ecosystems, which were ultimately important in determining the fates of many key groups, as well as the origins of many major extant lineages. We identify a series of potential abiotic candidates for driving these patterns, including multiple bolide impacts, several episodes of flood basalt eruptions, dramatic climate change, and major disruptions to oceanic systems. The J/K transition therefore, although not a mass extinction, represents an important transitional period in the co-evolutionary history of life on Earth.

Small RNAs from a Big Genome: The piRNA Pathway and Transposable Elements in the Salamander Species Desmognathus fuscus

Most of the largest vertebrate genomes are found in salamanders, a clade of amphibians that includes 686 species. Salamander genomes range in size from 14 to 120 Gb, reflecting the accumulation of large numbers of transposable element (TE) sequences from all three TE classes. Although DNA loss rates are slow in salamanders relative to other vertebrates, high levels of TE insertion are also likely required to explain such high TE loads. Across the Tree of Life, novel TE insertions are suppressed by several pathways involving small RNA molecules. In most known animals, TE activity in the germline is primarily regulated by the Piwi-interacting RNA (piRNA) pathway. In this study, we test the hypothesis that salamanders' unusually high TE loads reflect the loss of the ancestral piRNA-mediated TE-silencing machinery. We characterized the small RNA pool in the female and male adult gonads, testing for the presence of small RNA molecules that bear the characteristics of TE-targeting piRNAs. We also analyzed the amino acid sequences of piRNA pathway proteins from salamanders and other vertebrates, testing whether the overall patterns of sequence divergence are consistent with conserved pathway function across the vertebrate clade. Our results do not support the hypothesis of piRNA pathway loss; instead, they suggest that the piRNA pathway is expressed in salamanders. Given these results, we propose hypotheses to explain how the extraordinary TE loads in salamander genomes could have accumulated, despite the expression of TE-silencing machinery.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Recent assembly of the global herbaceous flora: evidence from the paper daisies (Asteraceae: Gnaphalieae)

The global flora is thought to contain a large proportion of herbs, and understanding the general spatiotemporal processes that shaped the global distribution of these communities is one of the most difficult issues in biogeography. We explored patterns of world-wide biogeography in a species-rich herbaceous group, the paper daisy tribe Gnaphalieae (Asteraceae), based on the hitherto largest taxon sampling, a total of 835 terminal accessions representing 80% of the genera, and encompassing the global geographic range of the tribe, with nuclear internal transcribed spacer (ITS) and external transcribed spacer (ETS) sequences. Biogeographic analyses indicate that Gnaphalieae originated in southern Africa during the Oligocene, followed by repeated migrations into the rest of Africa and the Mediterranean region, with subsequent entries into other continents during various periods starting in the Miocene. Expansions in the late Miocene to Pliocene appear to have been the driving force that shaped the global distribution of the tribe as forests were progressively broken up by the mid-continent aridification and savannas and grasslands expanded into the interior of the major continents. This pattern of recent colonizations may explain the world-wide distribution of many other organisms in open ecosystems and it is highlighted here as an emerging pattern in the evolution of the global flora.

Geomolecular Dating and the Origin of Placental Mammals

In modern evolutionary divergence analysis the role of geological information extends beyond providing a timescale, to informing molecular rate variation across the tree. Here I consider the implications of this development. I use fossil calibrations to test the accuracy of models of molecular rate evolution for placental mammals, and reveal substantial misspecification associated with life history rate correlates. Adding further calibrations to reduce dating errors at specific nodes unfortunately tends to transfer underlying rate errors to adjacent branches. Thus, tight calibration across the tree is vital to buffer against rate model errors. I argue that this must include allowing maximum bounds to be tight when good fossil records permit, otherwise divergences deep in the tree will tend to be inflated by the interaction of rate errors and asymmetric confidence in minimum and maximum bounds. In the case of placental mammals I sought to reduce the potential for transferring calibration and rate model errors across the tree by focusing on well-supported calibrations with appropriately conservative maximum bounds. The resulting divergence estimates are younger than others published recently, and provide the long-anticipated molecular signature for the placental mammal radiation observed in the fossil record near the 66 Ma Cretaceous-Paleogene extinction event.

An Extremely Peramorphic Newt (Urodela: Salamandridae: Pleurodelini) from the Latest Oligocene of Germany, and a New Phylogenetic Analysis of Extant and Extinct Salamandrids

We describe an Oligocene newt specimen from western Germany that has gone practically unnoticed in the literature despite having been housed in the Museum für Naturkunde (Berlin) for a century. It is referable to the coeval Chelotriton, but is unusually peramorphic; for many characters it is more peramorphic than all other caudates or even all other lissamphibians. Most noticeable are the position of the jaw joints far caudal to the occiput, the honeycombed sculpture on the maxilla, and the possible presence of a septomaxilla (which would be unique among salamandrids). Referral to a species would require a revision of the genus, but the specimen likely does not belong to the type species. A phylogenetic analysis of nonmolecular characters of Salamandridae, far larger than all predecessors, confirms the referral to Chelotriton. It further loosely associates the Oligocene Archaeotriton and the Miocene Carpathotriton with the extant Lissotriton, though the former may alternatively lie outside Pleurodelinae altogether. The Miocene? I. randeckensis may not belong to the extant Ichthyosaura. The Miocene “Triturus” roehrsi is found neither with the extant Ommatotriton nor with Lissotriton, but inside an Asian/aquatic clade or, when geographic distribution is included as a character, as the sister-group to all other European molgins. The main cause for discrepancies between the results and the molecular consensus is not heterochrony, but adaptations to a life in mountain streams; this is the most likely reason why the Paleocene Koalliella from western Europe forms the sister-group to some or all of the most aquatic extant newts in different analyses. We would like to urge neontologists working on salamandrids to pay renewed attention to the skeleton, not limited to the skull, as a source of diagnostic and phylogenetically informative characters.

Potential for bias and low precision in molecular divergence time estimation of the Canopy of Life: an example from aquatic bird families

Uncertainty in divergence time estimation is frequently studied from many angles but rarely from the perspective of phylogenetic node age. If appropriate molecular models and fossil priors are used, a multi-locus, partitioned analysis is expected to equally minimize error in accuracy and precision across all nodes of a given phylogeny. In contrast, if available models fail to completely account for rate heterogeneity, substitution saturation and incompleteness of the fossil record, uncertainty in divergence time estimation may increase with node age. While many studies have stressed this concern with regard to deep nodes in the Tree of Life, the inference that molecular divergence time estimation of shallow nodes is less sensitive to erroneous model choice has not been tested explicitly in a Bayesian framework. Because of available divergence time estimation methods that permit fossil priors across any phylogenetic node and the present increase in efficient, cheap collection of species-level genomic data, insight is needed into the performance of divergence time estimation of shallow (<10 MY) nodes. Here, we performed multiple sensitivity analyses in a multi-locus data set of aquatic birds with six fossil constraints. Comparison across divergence time analyses that varied taxon and locus sampling, number and position of fossil constraint and shape of prior distribution showed various insights. Deviation from node ages obtained from a reference analysis was generally highest for the shallowest nodes but determined more by temporal placement than number of fossil constraints. Calibration with only the shallowest nodes significantly underestimated the aquatic bird fossil record, indicating the presence of saturation. Although joint calibration with all six priors yielded ages most consistent with the fossil record, ages of shallow nodes were overestimated. This bias was found in both mtDNA and nDNA regions. Thus, divergence time estimation of shallow nodes may suffer from bias and low precision, even when appropriate fossil priors and best available substitution models are chosen. Much care must be taken to address the possible ramifications of substitution saturation across the entire Tree of Life.

SCPP genes in the coelacanth: tissue mineralization genes shared by sarcopterygians

The coelacanth is the basal-most extant sarcopterygian that has teeth and tooth-like structures, comprising bone, dentin, and enamel or enameloid. Formation of these tissues involves many members of the secretory calcium-binding protein (SCPP) family. In tetrapods, acidic-residue-rich SCPPs are used in mineralization of bone and dentin, whereas Pro/Gln-rich SCPPs participate in enamel formation. Teleosts also employ many SCPPs for tissue mineralization. Nevertheless, the repertoire of SCPPs is largely different in teleosts and tetrapods; hence, filling this gap would be critical to elucidate early evolution of mineralized tissues in osteichthyans. In the present study, we searched for SCPP genes in the coelacanth genome and identified 11, of which two have clear orthologs in both tetrapods and teleosts, seven only in tetrapods, and two in neither of them. Given the divergence times of these vertebrate lineages, our discovery of this many SCPP genes shared between the coelacanth and tetrapods, but not with teleosts, suggests a complicated evolutionary scheme of SCPP genes in early osteichthyans. Our investigation also revealed both conserved and derived characteristics of SCPPs in the coelacanth and other vertebrates. Notably, acidic SCPPs independently evolved various acidic repeats in different lineages, while maintaining high acidity, presumably important for interactions with calcium. Furthermore, the three Pro/Gln-rich SCPP genes, required for mineralizing enamel matrix and confirmed only in tetrapods, were all identified in the coelacanth, strongly suggesting that enamel is equivalent in the coelacanth and tetrapods. This finding corroborates the previous proposition that true enamel evolved much earlier than the origin of tetrapods.

Hellbender genome sequences shed light on genomic expansion at the base of crown salamanders

Among animals, genome sizes range from 20 Mb to 130 Gb, with 380-fold variation across vertebrates. Most of the largest vertebrate genomes are found in salamanders, an amphibian clade of 660 species. Thus, salamanders are an important system for studying causes and consequences of genomic gigantism. Previously, we showed that plethodontid salamander genomes accumulate higher levels of long terminal repeat (LTR) retrotransposons than do other vertebrates, although the evolutionary origins of such sequences remained unexplored. We also showed that some salamanders in the family Plethodontidae have relatively slow rates of DNA loss through small insertions and deletions. Here, we present new data from Cryptobranchus alleganiensis, the hellbender. Cryptobranchus and Plethodontidae span the basal phylogenetic split within salamanders; thus, analyses incorporating these taxa can shed light on the genome of the ancestral crown salamander lineage, which underwent expansion. We show that high levels of LTR retrotransposons likely characterize all crown salamanders, suggesting that disproportionate expansion of this transposable element (TE) class contributed to genomic expansion. Phylogenetic and age distribution analyses of salamander LTR retrotransposons indicate that salamanders' high TE levels reflect persistence and diversification of ancestral TEs rather than horizontal transfer events. Finally, we show that relatively slow DNA loss rates through small indels likely characterize all crown salamanders, suggesting that a decreased DNA loss rate contributed to genomic expansion at the clade's base. Our identification of shared genomic features across phylogenetically distant salamanders is a first step toward identifying the evolutionary processes underlying accumulation and persistence of high levels of repetitive sequence in salamander genomes.

Low levels of LTR retrotransposon deletion by ectopic recombination in the gigantic genomes of salamanders

Across the tree of life, species vary dramatically in nuclear genome size. Mutations that add or remove sequences from genomes-insertions or deletions, or indels-are the ultimate source of this variation. Differences in the tempo and mode of insertion and deletion across taxa have been proposed to contribute to evolutionary diversity in genome size. Among vertebrates, most of the largest genomes are found within the salamanders, an amphibian clade with genome sizes ranging from ~14 to ~120 Gb. Salamander genomes have been shown to experience slower rates of DNA loss through small (i.e., <30 bp) deletions than do other vertebrate genomes. However, no studies have addressed DNA loss from salamander genomes resulting from larger deletions. Here, we focus on one type of large deletion-ectopic-recombination-mediated removal of LTR retrotransposon sequences. In ectopic recombination, double-strand breaks are repaired using a "wrong" (i.e., ectopic, or non-allelic) template sequence-typically another locus of similar sequence. When breaks occur within the LTR portions of LTR retrotransposons, ectopic-recombination-mediated repair can produce deletions that remove the internal transposon sequence and the equivalent of one of the two LTR sequences. These deletions leave a signature in the genome-a solo LTR sequence. We compared levels of solo LTRs in the genomes of four salamander species with levels present in five vertebrates with smaller genomes. Our results demonstrate that salamanders have low levels of solo LTRs, suggesting that ectopic-recombination-mediated deletion of LTR retrotransposons occurs more slowly than in other vertebrates with smaller genomes.

Methods for the quantitative comparison of molecular estimates of clade age and the fossil record

Approaches quantifying the relative congruence, or incongruence, of molecular divergence estimates and the fossil record have been limited. Previously proposed methods are largely node specific, assessing incongruence at particular nodes for which both fossil data and molecular divergence estimates are available. These existing metrics, and other methods that quantify incongruence across topologies including entirely extinct clades, have so far not taken into account uncertainty surrounding both the divergence estimates and the ages of fossils. They have also treated molecular divergence estimates younger than previously assessed fossil minimum estimates of clade age as if they were the same as cases in which they were older. However, these cases are not the same. Recovered divergence dates younger than compared oldest known occurrences require prior hypotheses regarding the phylogenetic position of the compared fossil record and standard assumptions about the relative timing of morphological and molecular change to be incorrect. Older molecular dates, by contrast, are consistent with an incomplete fossil record and do not require prior assessments of the fossil record to be unreliable in some way. Here, we compare previous approaches and introduce two new descriptive metrics. Both metrics explicitly incorporate information on uncertainty by utilizing the 95% confidence intervals on estimated divergence dates and data on stratigraphic uncertainty concerning the age of the compared fossils. Metric scores are maximized when these ranges are overlapping. MDI (minimum divergence incongruence) discriminates between situations where molecular estimates are younger or older than known fossils reporting both absolute fit values and a number score for incompatible nodes. DIG range (divergence implied gap range) allows quantification of the minimum increase in implied missing fossil record induced by enforcing a given set of molecular-based estimates. These metrics are used together to describe the relationship between time trees and a set of fossil data, which we recommend be phylogenetically vetted and referred on the basis of apomorphy. Differences from previously proposed metrics and the utility of MDI and DIG range are illustrated in three empirical case studies from angiosperms, ostracods, and birds. These case studies also illustrate the ways in which MDI and DIG range may be used to assess time trees resultant from analyses varying in calibration regime, divergence dating approach or molecular sequence data analyzed.

Beyond fossil calibrations: realities of molecular clock practices in evolutionary biology

Molecular-based divergence dating methods, or molecular clocks, are the primary neontological tool for estimating the temporal origins of clades. While the appropriate use of vertebrate fossils as external clock calibrations has stimulated heated discussions in the paleontological community, less attention has been given to the quality and implementation of other calibration types. In lieu of appropriate fossils, many studies rely on alternative sources of age constraints based on geological events, substitution rates and heterochronous sampling, as well as dates secondarily derived from previous analyses. To illustrate the breadth and frequency of calibration types currently employed, we conducted a literature survey of over 600 articles published from 2007 to 2013. Over half of all analyses implemented one or more fossil dates as constraints, followed by geological events and secondary calibrations (15% each). Vertebrate taxa were subjects in nearly half of all studies, while invertebrates and plants together accounted for 43%, followed by viruses, protists and fungi (3% each). Current patterns in calibration practices were disproportionate to the number of discussions on their proper use, particularly regarding plants and secondarily derived dates, which are both relatively neglected in methodological evaluations. Based on our survey, we provide a comprehensive overview of the latest approaches in clock calibration, and outline strengths and weaknesses associated with each. This critique should serve as a call to action for researchers across multiple communities, particularly those working on clades for which fossil records are poor, to develop their own guidelines regarding selection and implementation of alternative calibration types. This issue is particularly relevant now, as time-calibrated phylogenies are used for more than dating evolutionary origins, but often serve as the backbone of investigations into biogeography, diversity dynamics and rates of phenotypic evolution.

A comparative cluster analysis of nicotinamide adenine dinucleotide phosphate (NADPH)-diaphorase histochemistry in the brains of amphibians

Nicotinamide adenine dinucleotide phosphate-diaphorase (NADPH-d) is a key enzyme in the synthesis of the gaseous neurotransmitter nitric oxide. We compare the distribution of NADPH-d in the brain of four species of hylid frogs. NADPH-d-positive fibers are present throughout much of the brain, whereas stained cell groups are distributed in well-defined regions. Whereas most brain areas consistently show positive neurons in all species, in some areas species-specific differences occur. We analyzed our data and those available for other amphibian species to build a matrix on NADPH-d brain distribution for a multivariate analysis. Brain dissimilarities were quantified by using the Jaccard index in a hierarchical clustering procedure. The whole brain dendrogram was compared with that of its main subdivisions by applying the Fowlkes-Mallows index for dendrogram similarity, followed by bootstrap replications and a permutation test. Despite the differences in the distribution map of the NADPH-d system among species, cluster analysis of data from the whole brain and hindbrain faithfully reflected the evolutionary history (framework) of amphibians. Dendrograms from the secondary prosencephalon, diencephalon, mesencephalon, and isthmus showed some deviation from the main scheme. Thus, the present analysis supports the major evolutionary stability of the hindbrain. We provide evidence that the NADPH-d system in main brain subdivisions should be cautiously approached for comparative purposes because specific adaptations of a single species could occur and may affect the NADPH-d distribution pattern in a brain subdivision. The minor differences in staining pattern of particular subdivisions apparently do not affect the general patterns of staining across species.

On the fossil record of the Gekkota

Gekkota is often interpreted as sister to all remaining squamates, exclusive of dibamids, or as sister to Autarchoglossa. It is the only diverse lineage of primarily nocturnal lizards and includes some of the smallest amniotes. The skeleton of geckos has often been interpreted as paedomorphic and/or "primitive" but these lizards also display a wide range of structural specializations of the postcranium, including modifications associated with both scansorial locomotion and limb reduction. Although the concept of "Gekkota" has been variously applied by different authors, we here apply a rigorous apomorphy based definition, recent advances in gekkotan morphology and phylogenetics, and diverse comparative material to provide a comprehensive assessment of 28 known pre-Quaternary geckos, updating the last such review, published three decades ago. Fossils evaluated include both sedimentary fossils and amber-embedded specimens. Known Cretaceous geckos are exclusively Asian and exhibit character combinations not seen in any living forms. Cenozoic gekkotans derive from sites around the world, although Europe is especially well represented. Paleogene geckos are largely known from disarticulated remains and show similarities to Sphaerodactylidae and Diplodactylidae, although resemblances may be plesiomorphic in some cases. Many Neogene gekkotans are referable to living families or even genera, but their geographic occurrences are often extralimital to those of modern groups, as is consistent with paleoclimatic conditions. The phylogenetic placement of fossil gekkotans has important repercusions for timetree calibration, but at present only a small number of fossils can be confidently assigned to even family level groupings, limiting their utility in this regard.

Large-scale phylogenetic analyses reveal the causes of high tropical amphibian diversity

Many groups show higher species richness in tropical regions but the underlying causes remain unclear. Despite many competing hypotheses to explain latitudinal diversity gradients, only three processes can directly change species richness across regions: speciation, extinction and dispersal. These processes can be addressed most powerfully using large-scale phylogenetic approaches, but most previous studies have focused on small groups and recent time scales, or did not separate speciation and extinction rates. We investigate the origins of high tropical diversity in amphibians, applying new phylogenetic comparative methods to a tree of 2871 species. Our results show that high tropical diversity is explained by higher speciation in the tropics, higher extinction in temperate regions and limited dispersal out of the tropics compared with colonization of the tropics from temperate regions. These patterns are strongly associated with climate-related variables such as temperature, precipitation and ecosystem energy. Results from models of diversity dependence in speciation rate suggest that temperate clades may have lower carrying capacities and may be more saturated (closer to carrying capacity) than tropical clades. Furthermore, we estimate strikingly low tropical extinction rates over geological time scales, in stark contrast to the dramatic losses of diversity occurring in tropical regions presently.

Combining fossil and molecular data to date the diversification of New World Primates

Recent methodological advances in molecular dating associated with the growing availability of sequence data have prompted the study of the evolution of New World Anthropoidea in recent years. Motivated by questions regarding historical biogeography or the mode of evolution, these works aimed to obtain a clearer scenario of Platyrrhini origins and diversification. Although some consensus was found, disputed issues, especially those relating to the evolutionary affinities of fossil taxa, remain. The use of fossil taxa for divergence time analysis is traditionally restricted to the provision of calibration priors. However, new analytical approaches have been developed that incorporate fossils as terminals and, thus, directly assign ages to the fossil tips. In this study, we conducted a combined analysis of molecular and morphological data, including fossils, to derive the timescale of New World anthropoids. Differently from previous studies that conducted total-evidence analysis of molecules and morphology, our approach investigated the morphological clock alone. Our results corroborate the hypothesis that living platyrrhines diversified in the last 20 Ma and that Miocene Patagonian fossils compose an independent evolutionary radiation that diversified in the late Oligocene. When compared to the node ages inferred from the molecular timescale, the inclusion of fossils augmented the precision of the estimates for nodes constrained by the fossil tips. We show that morphological data can be analysed using the same methodological framework applied in relaxed molecular clock studies.

Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs

Anura (frogs and toads) constitute over 88% of living amphibian diversity but many important questions about their phylogeny and evolution remain unresolved. For this study, we developed an efficient method for sequencing anuran mitochondrial DNAs (mtDNAs) by amplifying the mitochondrial genome in 12 overlapping fragments using frog-specific universal primer sets. Based on this method, we generated 47 nearly complete, new anuran mitochondrial genomes and discovered nine novel gene arrangements. By combining the new data and published anuran mitochondrial genomes, we assembled a large mitogenomic data set (11,007 nt) including 90 frog species, representing 39 of 53 recognized anuran families, to investigate their phylogenetic relationships and evolutionary history. The resulting tree strongly supported a paraphyletic arrangement of archaeobatrachian (=nonneobatrachian) frogs, with Leiopelmatoidea branching first, followed by Discoglossoidea, Pipoidea, and Pelobatoidea. Within Neobatrachia, the South African Heleophrynidae is the sister-taxon to all other neobatrachian frogs and the Seychelles-endemic Sooglossidae is recovered as the sister-taxon to Ranoidea. These phylogenetic relationships agree with many nuclear gene studies. The chronogram derived from two Bayesian relaxed clock methods (MultiDivTime and BEAST) suggests that modern frogs (Anura) originated in the early Triassic about 244 Ma and the appearance of Neobatrachia took place in the late Jurassic about 163 Ma. The initial diversifications of two species-rich superfamilies Hyloidea and Ranoidea commenced 110 and 133 Ma, respectively. These times are older than some other estimates by approximately 30-40 My. Compared with nuclear data, mtDNA produces compatible time estimates for deep nodes (>150 Ma), but apparently older estimates for more shallow nodes. Our study shows that, although it evolves relatively rapidly and behaves much as a single locus, mtDNA performs well for both phylogenetic and divergence time inferences and will provide important reference hypotheses for the phylogeny and evolution of frogs.