Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia

Mitochondrial Genome Characteristics Reveal Evolution of Acanthopsetta nadeshnyi (Jordan and Starks, 1904) and Phylogenetic Relationships

In the present study, the mitochondrial genomic characteristics of Acanthopsetta nadeshnyi have been reported and have depicted the phylogenetic relationship among Pleuronectidae. Combined with a comparative analysis of 13 PCGs, the TN93 model was used to review the neutral evolution and habitat evolution catalysis of the mitogenome to verify the distancing and purification selectivity of the mitogenome in Pleuronectidae. At the same time, a species differentiation and classification model based on mitogenome analysis data was established. This study is expected to provide a new perspective on the phylogenetic relationship and taxonomic status of A. nadeshnyi and lay a foundation for further exploration of environmental and biological evolutionary mechanisms.

Colonization of the ocean floor by jawless vertebrates across three mass extinctions

BACKGROUND

The deep (> 200 m) ocean floor is often considered to be a refugium of biodiversity; many benthic marine animals appear to share ancient common ancestry with nearshore and terrestrial relatives. Whether this pattern holds for vertebrates is obscured by a poor understanding of the evolutionary history of the oldest marine vertebrate clades. Hagfishes are jawless vertebrates that are either the living sister to all vertebrates or form a clade with lampreys, the only other surviving jawless fishes.

RESULTS

We use the hagfish fossil record and molecular data for all recognized genera to construct a novel hypothesis for hagfish relationships and diversification. We find that crown hagfishes persisted through three mass extinctions after appearing in the Permian ~ 275 Ma, making them one of the oldest living vertebrate lineages. In contrast to most other deep marine vertebrates, we consistently infer a deep origin of continental slope occupation by hagfishes that dates to the Paleozoic. Yet, we show that hagfishes have experienced marked body size diversification over the last hundred million years, contrasting with a view of this clade as morphologically stagnant.

CONCLUSION

Our results establish hagfishes as ancient members of demersal continental slope faunas and suggest a prolonged accumulation of deep sea jawless vertebrate biodiversity.

Bayesian Total-Evidence Dating Revisits Sloth Phylogeny and Biogeography: A Cautionary Tale on Morphological Clock Analyses

Combining morphological and molecular characters through Bayesian total-evidence dating allows inferring the phylogenetic and timescale framework of both extant and fossil taxa, while accounting for the stochasticity and incompleteness of the fossil record. Such an integrative approach is particularly needed when dealing with clades such as sloths (Mammalia: Folivora), for which developmental and biomechanical studies have shown high levels of morphological convergence whereas molecular data can only account for a limited percentage of their total species richness. Here, we propose an alternative hypothesis of sloth evolution that emphasizes the pervasiveness of morphological convergence and the importance of considering the fossil record and an adequate taxon sampling in both phylogenetic and biogeographic inferences. Regardless of different clock models and morphological datasets, the extant sloth Bradypus is consistently recovered as a megatherioid, and Choloepus as a mylodontoid, in agreement with molecular-only analyses. The recently extinct Caribbean sloths (Megalocnoidea) are found to be a monophyletic sister-clade of Megatherioidea, in contrast to previous phylogenetic hypotheses. Our results contradict previous morphological analyses and further support the polyphyly of "Megalonychidae," whose members were found in five different clades. Regardless of taxon sampling and clock models, the Caribbean colonization of sloths is compatible with the exhumation of islands along Aves Ridge and its geological time frame. Overall, our total-evidence analysis illustrates the difficulty of positioning highly incomplete fossils, although a robust phylogenetic framework was recovered by an a posteriori removal of taxa with high percentages of missing characters. Elimination of these taxa improved topological resolution by reducing polytomies and increasing node support. However, it introduced a systematic and geographic bias because most of these incomplete specimens are from northern South America. This is evident in biogeographic reconstructions, which suggest Patagonia as the area of origin of many clades when taxa are underrepresented, but Amazonia and/or Central and Southern Andes when all taxa are included. More generally, our analyses demonstrate the instability of topology and divergence time estimates when using different morphological datasets and clock models and thus caution against making macroevolutionary inferences when node support is weak or when uncertainties in the fossil record are not considered.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Dating the evolution of the complex thalloid liverworts (Marchantiopsida): total-evidence dating analysis supports a Late Silurian-Early Devonian origin and post-Mesozoic morphological stasis

Divergence times based on molecular clock analyses often differ from those derived from total-evidence dating (TED) approaches. For bryophytes, fossils have been excluded from previous assessments of divergence times, and thus, their utility in dating analyses remains unexplored. Here, we conduct the first TED analyses of the complex thalloid liverworts (Marchantiopsida) that include fossils and evaluate macroevolutionary trends in morphological 'diversity' (disparity) and rates. Phylogenetic analyses were performed on a combined dataset of 130 discrete characters and 11 molecular markers (sampled from nuclear, plastid and mitochondrial genomes). Taxon sampling spanned 56 extant species - representing all the orders within Marchantiophyta and extant genera within Marchantiales - and eight fossil taxa. Total-evidence dating analyses support the radiation of Marchantiopsida during Late Silurian-Early Devonian (or Middle Ordovician when the outgroup is excluded) and that of Ricciaceae in the Middle Jurassic. Morphological change rate was high early in the history of the group, but it barely increased after Late Cretaceous. Disparity-through-time analyses support a fast increase in diversity until the Middle Triassic (c. 250 Ma), after which phenotypic evolution slows down considerably. Incorporating fossils in analyses challenges previous assumptions on the affinities of extinct taxa and indicates that complex thalloid liverworts radiated c. 125 Ma earlier than previously inferred.

The wood frog (Rana sylvatica): An emerging comparative model for anuran immunity and host-ranavirus interactions

The wood frog (Rana sylvatica) is widely distributed across North America and is the only amphibian found north of the Arctic Circle due to its remarkable ability to tolerate whole-body freezing. Recent mass mortalities attributable to Ranavirus spp. (family Iridoviridae) in wild juvenile wood frogs, coupled with the apparent high susceptibility of wood frogs to experimental infection with frog virus 3 (FV3), the type species of the Ranavirus genus, or FV3-like isolates underscore the serious threat ranaviruses poses to wood frog populations. Despite the ecological relevance and unique life history of wood frogs, our understanding of the wood frog immune system and antiviral response to ranaviral infections is in its infancy. Here we aim to (1) synthesize the limited knowledge of wood frog immune defences, (2) review recent progress in establishing the wood frog as a study system for ranavirus infection, and (3) highlight the future use of wood frogs as a model anuran to provide insight into the evolution of anuran immune systems and antiviral responses.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Handling Logical Character Dependency in Phylogenetic Inference: Extensive Performance Testing of Assumptions and Solutions Using Simulated and Empirical Data

Logical character dependency is a major conceptual and methodological problem in phylogenetic inference of morphological data sets, as it violates the assumption of character independence that is common to all phylogenetic methods. It is more frequently observed in higher-level phylogenies or in data sets characterizing major evolutionary transitions, as these represent parts of the tree of life where (primary) anatomical characters either originate or disappear entirely. As a result, secondary traits related to these primary characters become "inapplicable" across all sampled taxa in which that character is absent. Various solutions have been explored over the last three decades to handle character dependency, such as alternative character coding schemes and, more recently, new algorithmic implementations. However, the accuracy of the proposed solutions, or the impact of character dependency across distinct optimality criteria, has never been directly tested using standard performance measures. Here, we utilize simple and complex simulated morphological data sets analyzed under different maximum parsimony optimization procedures and Bayesian inference to test the accuracy of various coding and algorithmic solutions to character dependency. This is complemented by empirical analyses using a recoded data set on palaeognathid birds. We find that in small, simulated data sets, absent coding performs better than other popular coding strategies available (contingent and multistate), whereas in more complex simulations (larger data sets controlled for different tree structure and character distribution models) contingent coding is favored more frequently. Under contingent coding, a recently proposed weighting algorithm produces the most accurate results for maximum parsimony. However, Bayesian inference outperforms all parsimony-based solutions to handle character dependency due to fundamental differences in their optimization procedures-a simple alternative that has been long overlooked. Yet, we show that the more primary characters bearing secondary (dependent) traits there are in a data set, the harder it is to estimate the true phylogenetic tree, regardless of the optimality criterion, owing to a considerable expansion of the tree parameter space. [Bayesian inference, character dependency, character coding, distance metrics, morphological phylogenetics, maximum parsimony, performance, phylogenetic accuracy.].

Estimating the Age of Poorly Dated Fossil Specimens and Deposits Using a Total-Evidence Approach and the Fossilized Birth-Death Process

Bayesian total-evidence approaches under the fossilized birth-death model enable biologists to combine fossil and extant data while accounting for uncertainty in the ages of fossil specimens, in an integrative phylogenetic analysis. Fossil age uncertainty is a key feature of the fossil record as many empirical data sets may contain a mix of precisely dated and poorly dated fossil specimens or deposits. In this study, we explore whether reliable age estimates for fossil specimens can be obtained from Bayesian total-evidence phylogenetic analyses under the fossilized birth-death model. Through simulations based on the example of the Baltic amber deposit, we show that estimates of fossil ages obtained through such an analysis are accurate, particularly when the proportion of poorly dated specimens remains low and the majority of fossil specimens have precise dates. We confirm our results using an empirical data set of living and fossil penguins by artificially increasing the age uncertainty around some fossil specimens and showing that the resulting age estimates overlap with the recorded age ranges. Our results are applicable to many empirical data sets where classical methods of establishing fossil ages have failed, such as the Baltic amber and the Gobi Desert deposits. [Bayesian phylogenetic inference; fossil age estimates; fossilized birth-death; Lagerstätte; total-evidence.].

Syngnathoid Evolutionary History and the Conundrum of Fossil Misplacement

Seahorses, pipefishes, trumpetfishes, shrimpfishes, and allies are a speciose, globally distributed clade of fishes that have evolved a large number of unusual body plans. The clade that includes all these forms, Syngnathoidei, has become a model for the study of life history evolution, population biology, and biogeography. Yet, the timeline of syngnathoid evolution has remained highly contentious. This debate is largely attributable to the nature of the syngnathoid fossil record, which is both poorly described and patchy for several major lineages. Although fossil syngnathoids have been used to calibrate molecular phylogenies, the interrelationships of extinct species and their affinities to major living syngnathoid clades have scarcely been quantitatively tested. Here, I use an expanded morphological dataset to reconstruct the evolutionary relationships and clade ages of fossil and extant syngnathoids. Phylogenies generated using different analytical methodologies are largely congruent with molecular phylogenetic trees of Syngnathoidei but consistently find novel placements for several key taxa used as fossil calibrators in phylogenomic studies. Tip-dating of the syngnathoid phylogeny finds a timeline for their evolution that differs slightly from the one inferred using molecular trees but is generally congruent with a post-Cretaceous diversification event. These results emphasize the importance of quantitatively testing the relationships of fossil species, particularly when they are critical to assessing divergence times.

Characterization of immunoglobulin loci in the gigantic genome of Ambystoma mexicanum

Background

The axolotl, Ambystoma mexicanum is a unique biological model for complete tissue regeneration. Is a neotenic endangered species and is highly susceptible to environmental stress, including infectious disease. In contrast to other amphibians, the axolotl is particularly vulnerable to certain viral infections. Like other salamanders, the axolotl genome is one of the largest (32 Gb) and the impact of genome size on Ig loci architecture is unknown. To better understand the immune response in axolotl, we aimed to characterize the immunoglobulin loci of A. mexicanum and compare it with other model vertebrates.

Methods

The most recently published genome sequence of A. mexicanum (V6) was used for alignment-based annotation and manual curation using previously described axolotl Ig sequences or reference sequences from other vertebrates. Gene models were further curated using A. mexicanum spleen RNA-seq data. Human, Xenopus tropicalis, Danio rerio (zebrafish), and eight tetrapod reference genomes were used for comparison.

Results

Canonical A. mexicanum heavy chain (IGH), lambda (IGL), sigma (IGS), and the putative surrogate light chain (SLC) loci were identified. No kappa locus was found. More than half of the IGHV genes and the IGHF gene are pseudogenes and there is no clan I IGHV genes. Although the IGH locus size is proportional to genome size, we found local size restriction in the IGHM gene and the V gene intergenic distances. In addition, there were V genes with abnormally large V-intron sizes, which correlated with loss of gene functionality.

Conclusion

The A. mexicanum immunoglobulin loci share the same general genome architecture as most studied tetrapods. Consistent with its large genome, Ig loci are larger; however, local size restrictions indicate evolutionary constraints likely to be imposed by high transcriptional demand of certain Ig genes, as well as the V(D)J recombination over very long genomic distance ranges. The A. mexicanum has undergone an extensive process of Ig gene loss which partially explains a reduced potential repertoire diversity that may contribute to its impaired antibody response.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

Relative performance of Bayesian morphological clock and parsimony methods for phylogenetic reconstructions: Insights from the case of Myomiminae and Dryomyinae glirid rodents

Extinct organisms provide crucial information about the origin and time of origination of extant groups. The importance of morphological phylogenetics for rigorously dating the tree of life is now widely recognized and has been revitalized by methodological developments such as the application of tip-dating Bayesian approaches. Traditionally, molecular clocks have been node calibrated. However, node calibrations are often unsatisfactory because they do not allow the fossil age to inform about phylogenetic hypothesis. The introduction of tip calibrations allow fossil species to be included alongside their living relatives, and the absence of molecular sequence data for these taxa to be remedied by supplementing the sequence alignments for living taxa with phenotype character matrices for both living and fossil taxa. Therefore, only phylogenetic analyses that take into account morphological characters can incorporate both fossil and extant species. Herein we present an unprecedented morphological dataset for a vast group of glirid rodents, to which different phylogenetic methodologies have been applied. We have compared the tree topologies resulting from traditional parsimony and Bayesian phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods seem to outperform parsimony with our dataset. The strict consensus tree recovered by tip dating invalidates the classic classification and allows dates to be proposed for the divergence and origin of the different clades.

The Genetic Diversity and the Divergence Time in Extant Primitive Mayfly, Siphluriscus chinensis Ulmer, 1920 Using the Mitochondrial Genome

In this study, the mitochondrial (mt) genomes of Siphluriscus chinensis (Ephemeroptera: Siphluriscidae) were evaluated in specimens collected from two sites in China: Niutou Mountain, Zhejiang Province (S. chinensis NTS) and Leigong Mountain, Guizhou Province (S. chinensis LGS) and were successfully sequenced. The lengths of the mt genomes of S. chinensis NTS and S. chinensis LGS were 15,904 bp (ON729390) and 15,212 bp (ON729391), respectively. However, an in-depth comparison of the two mt genomes showed significant differences between the specimens collected from the two sites. A detailed analysis of the genetic distance between S. chinensis NTS and S. chinensis LGS was undertaken to further achieve an accurate delimitation of S. chinensis. The genetic distance between S. chinensis NTS and the other three species within Siphluriscidae was a high value, above 12.2%. The two mt genomes were used to reconstruct phylogenetic relationships and estimate divergence time. The results demonstrated robust differences between S. chinensis NTS and S. chinensis LGS, which revealed that a kind of cryptic species existed. Maximum likelihood (ML) and Bayesian inference (BI) analyses produced well-supported phylogenetic trees that showed evolutionary relationships between Siphluriscidae (((S. chinensis HQ875717 + S. chinensis MF352165) + S. chinensis LGS) + S. chinensis NTS). The most recent common ancestor (MRCA) of four species within Siphluriscidae began to diversify during the Neogene [11.80 million years ago (Mya); 95% highest posterior densities (HPD) = 6.17–19.28 Mya], and S. chinensis NTS was first to diverge from the branches of S. chinensis LGS. In short, based on mitochondrial genomes, our results showed that the specimens collected from Leigong Mountain, Guizhou Province (S. chinensis LGS) belonged to S. chinensis, and the specimens collected from Niutou Mountain, Zhejiang Province (S. chinensis NTS) were a cryptic species of S. chinensis.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Beetle bioluminescence outshines extant aerial predators

We understand very little about the timing and origins of bioluminescence, particularly as a predator avoidance strategy. Understanding the timing of its origins, however, can help elucidate the evolution of this ecologically important signal. Using fireflies, a prevalent bioluminescent group where bioluminescence primarily functions as aposematic and sexual signals, we explore the origins of this signal in the context of their potential predators. Divergence time estimations were performed using genomic-scale datasets providing a robust estimate for the origin of firefly bioluminescence as both a terrestrial and as an aerial signal. Our results recover the origin of terrestrial beetle bioluminescence at 141.17 (122.63-161.17) Ma and firefly aerial bioluminescence at 133.18 (117.86-152.47) Ma using a large dataset focused on Lampyridae; and terrestrial bioluminescence at 148.03 (130.12-166.80) Ma, with the age of aerial bioluminescence at 104.97 (99.00-120.90) Ma using a complementary Elateroidea dataset. These ages pre-date the origins of all known extant aerial predators (i.e. bats and birds) and support much older terrestrial predators (assassin bugs, frogs, ground beetles, lizards, snakes, hunting spiders and harvestmen) as the drivers of terrestrial bioluminescence in beetles. These ages also support the hypothesis that sexual signalling was probably the original function of this signal in aerial fireflies.

Middle Jurassic fossils document an early stage in salamander evolution

Salamanders are an important group of living amphibians and model organisms for understanding locomotion, development, regeneration, feeding, and toxicity in tetrapods. However, their origin and early radiation remain poorly understood, with early fossil stem-salamanders so far represented by larval or incompletely known taxa. This poor record also limits understanding of the origin of Lissamphibia (i.e., frogs, salamanders, and caecilians). We report fossils from the Middle Jurassic of Scotland representing almost the entire skeleton of the enigmatic stem-salamander Marmorerpeton. We use computed tomography to visualize high-resolution three-dimensional anatomy, describing morphologies that were poorly characterized in early salamanders, including the braincase, scapulocoracoid, and lower jaw. We use these data in the context of a phylogenetic analysis intended to resolve the relationships of early and stem-salamanders, including representation of important outgroups alongside data from high-resolution imaging of extant species. Marmorerpeton is united with Karaurus, Kokartus, and others from the Middle Jurassic-Lower Cretaceous of Asia, providing evidence for an early radiation of robustly built neotenous stem-salamanders. These taxa display morphological specializations similar to the extant cryptobranchid "giant" salamanders. Our analysis also demonstrates stem-group affinities for a larger sample of Jurassic species than previously recognized, highlighting an unappreciated diversity of stem-salamanders and cautioning against the use of single species (e.g., Karaurus) as exemplars for stem-salamander anatomy. These phylogenetic findings, combined with knowledge of the near-complete skeletal anatomy of Mamorerpeton, advance our understanding of evolutionary changes on the salamander stem-lineage and provide important data on early salamanders and the origins of Batrachia and Lissamphibia.

The role of paleontological data in bryophyte systematics

Systematics reconstructs tempo and mode in biological evolution by resolving the phylogenetic fabric of biodiversity. The staggering duration and complexity of evolution, coupled with loss of information (extinction), render exhaustive reconstruction of the evolutionary history of life unattainable. Instead, we sample its products-phenotypes and genotypes-to generate phylogenetic hypotheses, which we sequentially reassess and update against new data. Current consensus in evolutionary biology emphasizes fossil integration in total-evidence analyses, requiring in-depth understanding of fossils-age, phenotypes, and systematic affinities-and a detailed morphological framework uniting fossil and extant taxa. Bryophytes present a special case: deep evolutionary history but sparse fossil record and phenotypic diversity encompassing small dimensional scales. We review how these peculiarities shape fossil inclusion in bryophyte systematics. Paucity of the bryophyte fossil record, driven primarily by phenotypic (small plant size) and ecological constraints (patchy substrate-hugging populations), and incomplete exploration, results in many morphologically isolated, taxonomically ambiguous fossil taxa. Nevertheless, instances of exquisite preservation and pioneering studies demonstrate the feasibility of including bryophyte fossils in evolutionary inference. Further progress will arise from developing extensive morphological matrices for bryophytes, continued exploration of the fossil record, re-evaluation of previously described fossils, and training specialists in identification and characterization of bryophyte fossils, and in bryophyte morphology.

What is the age of flowering plants?

The origin of flowering plants (angiosperms) was one of the most transformative events in the history of our planet. Despite considerable interest from multiple research fields, numerous questions remain, including the age of the group as a whole. Recent studies have reported a perplexing range of estimates for the crown-group age of angiosperms, from ~140 million years (Ma; Early Cretaceous) to 270 Ma (Permian). Both ends of the spectrum are now supported by both macroevolutionary analyses of the fossil record and fossil-calibrated molecular dating analyses. Here, we first clarify and distinguish among the three ages of angiosperms: the age of their divergence with acrogymnosperms (stem age); the age(s) of emergence of their unique, distinctive features including flowers (morphological age); and the age of the most recent common ancestor of all their living species (crown age). We then demonstrate, based on recent studies, that fossil-calibrated molecular dating estimates of the crown-group age of angiosperms have little to do with either the amount of molecular data or the number of internal fossil calibrations included. Instead, we argue that this age is almost entirely conditioned by its own prior distribution (typically a calibration density set by the user in Bayesian analyses). Lastly, we discuss which future discoveries or novel types of analyses are most likely to bring more definitive answers. In the meantime, we propose that the age of angiosperms is best described as largely unknown (140-270 Ma) and that contrasting age estimates in the literature mostly reflect conflicting prior distributions. We also suggest that future work that depends on the time scale of flowering plant diversification be designed to integrate over this vexing uncertainty.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Combined phylogenetic analysis of Pleurodema (Anura: Leptodactylidae: Leiuperinae)

The genus Pleurodema comprises 15 species distributed through the Neotropical region, from sea level up to 5000 m.a.s.l. A total-evidence analysis of Pleurodema is provided based on the parsimony criterion. The combined dataset included morphometric, phenotypic, and DNA evidence (34 taxa, 4441 characters). The parsimony analysis yielded one most-parsimonious tree. Pleurodema was recovered as a well-supported clade composed of two major subclades. One subclade has an identical topology to that of previous analyses, the P. brachyops Clade (P. alium, P. borellii, P. brachyops, P. cinereum, P. diplolister, and P. tucumanum). The other subclade includes the remaining nine species of the genus, exhibiting a topology different from that of previous studies. According to the present phylogeny, this second lineage is formed by the P. nebulosum Clade (P. guayapae + P. nebulosum), P. marmoratum, the re-defined P. thaul Clade (P. bufoninum, P. somuncurense, P. thaul) and the P. bibroni Clade (P. bibroni, P. cordobae, P. kriegi). The reproductive modes of Pleurodema represent a unique combination of features within Leiuperinae, including three egg-clutch structures, two types of amplexus, and lack of vocalization. Also, some species of Pleurodema have been considered fossorial, because they are capable of digging with their hind-limbs and remaining in self-made burrows during dry seasons. The evolution of characters associated with reproductive biology and fossoriality is discussed in light of the obtained results.

Phylogeny indicates polyphyly in Cnodocentron (Trichoptera: Xiphocentronidae): biogeography and revision of New World species (Caenocentron)

Cnodocentron has a disjunct Laurasian distribution, with species in South-East Asia and the New World. It is divided into two subgenera: Cnodocentron and Caenocentron. Here, we infer the evolutionary history of the genus through phylogenetic and biogeographic data combining COI and 46 morphological characters. Phylogenetic relationships and divergence-time estimation were simultaneously inferred through Bayesian inference. The dating analysis was performed through relaxed morphological and molecular clocks. The historical biogeography was investigated using the dispersal–extinction–cladogenesis model. Our results indicate the polyphyly of the genus Cnodocentron, with each subgenus being more related to other genera than to each other. Therefore, the subgenus Caenocentron is elevated to genus status here. The biogeographical analysis showed that the Oriental Cnodocentron diverged in the Indian subcontinent during the Middle Eocene, while the New World Caenocentron stat. nov. originated in the Chortis Block in the Late Eocene. The dispersal of Caenocentron to South America occurred only after the Late Miocene, around 10 Mya. Additionally, we provide a revision of Caenocentron, with an identification key and description of the male and female of two new species from Costa Rica: Caenocentron carlosdelarosai sp. nov. and Caenocentron rafamoralesi sp. nov..

Returning to the roots: resolution, reproducibility, and robusticity in the phylogenetic inference of Dissorophidae (Amphibia: Temnospondyli)

The phylogenetic relationships of most Paleozoic tetrapod clades remain poorly resolved, which is variably attributed to a lack of study, the limitations of inference from phenotypic data, and constant revision of best practices. While refinement of phylogenetic methods continues to be important, any phylogenetic analysis is inherently constrained by the underlying dataset that it analyzes. Therefore, it becomes equally important to assess the accuracy of these datasets, especially when a select few are repeatedly propagated. While repeat analyses of these datasets may appear to constitute a working consensus, they are not in fact independent, and it becomes especially important to evaluate the accuracy of these datasets in order to assess whether a seeming consensus is robust. Here I address the phylogeny of the Dissorophidae, a speciose clade of Paleozoic temnospondyls. This group is an ideal case study among temnospondyls for exploring phylogenetic methods and datasets because it has been extensively studied (eight phylogenetic studies to date) but with most (six studies) using a single matrix that has been propagated with very little modification. In spite of the conserved nature of the matrix, dissorophid studies have produced anything but a conserved topology. Therefore, I analyzed an independently designed matrix, which recovered less resolution and some disparate nodes compared to previous studies. In order to reconcile these differences, I carefully examined previous matrices and analyses. While some differences are a matter of personal preference (e.g., analytical software), others relate to discrepancies with respect to what are currently considered as best practices. The most concerning discovery was the identification of pervasive dubious scorings that extend back to the origins of the widely propagated matrix. These include scores for skeletal features that are entirely unknown in a given taxon (e.g., postcrania in Cacops woehri) and characters for which there appear to be unstated working assumptions to scoring that are incompatible with the character definitions (e.g., scoring of taxa with incomplete skulls for characters based on skull length). Correction of these scores and other pervasive errors recovered a distinctly less resolved topology than previous studies, more in agreement with my own matrix. This suggests that previous analyses may have been compromised, and that the only real consensus of dissorophid phylogeny is the lack of one.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Relative time constraints improve molecular dating

Dating the tree of life is central to understanding the evolution of life on Earth. Molecular clocks calibrated with fossils represent the state of the art for inferring the ages of major groups. Yet, other information on the timing of species diversification can be used to date the tree of life. For example, horizontal gene transfer events and ancient coevolutionary interactions such as (endo)symbioses occur between contemporaneous species and thus can imply temporal relationships between two nodes in a phylogeny. Temporal constraints from these alternative sources can be particularly helpful when the geological record is sparse, for example, for microorganisms, which represent the majority of extant and extinct biodiversity. Here, we present a new method to combine fossil calibrations and relative age constraints to estimate chronograms. We provide an implementation of relative age constraints in RevBayes that can be combined in a modular manner with the wide range of molecular dating methods available in the software. We use both realistic simulations and empirical datasets of 40 Cyanobacteria and 62 Archaea to evaluate our method. We show that the combination of relative age constraints with fossil calibrations significantly improves the estimation of node ages. [Archaea, Bayesian analysis, cyanobacteria, dating, endosymbiosis, lateral gene transfer, MCMC, molecular clock, phylogenetic dating, relaxed molecular clock, revbayes, tree of life.]

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Molecular and morphological clocks for estimating evolutionary divergence times

BACKGROUND

Matrices of morphological characters are frequently used for dating species divergence times in systematics. In some studies, morphological and molecular character data from living taxa are combined, whereas others use morphological characters from extinct taxa as well. We investigated whether morphological data produce time estimates that are concordant with molecular data. If true, it will justify the use of morphological characters alongside molecular data in divergence time inference.

RESULTS

We systematically analyzed three empirical datasets from different species groups to test the concordance of species divergence dates inferred using molecular and discrete morphological data from extant taxa as test cases. We found a high correlation between their divergence time estimates, despite a poor linear relationship between branch lengths for morphological and molecular data mapped onto the same phylogeny. This was because node-to-tip distances showed a much higher correlation than branch lengths due to an averaging effect over multiple branches. We found that nodes with a large number of taxa often benefit from such averaging. However, considerable discordance between time estimates from molecules and morphology may still occur as  some intermediate nodes may show large time differences between these two types of data.

CONCLUSIONS

Our findings suggest that node- and tip-calibration approaches may be better suited for nodes with many taxa. Nevertheless, we highlight the importance of evaluating the concordance of intrinsic time structure in morphological and molecular data before any dating analysis using combined datasets.

Divergence-time estimates for hominins provide insight into encephalization and body mass trends in human evolution

Quantifying speciation times during human evolution is fundamental as it provides a timescale to test for the correlation between key evolutionary transitions and extrinsic factors such as climatic or environmental change. Here, we applied a total evidence dating approach to a hominin phylogeny to estimate divergence times under different topological hypotheses. The time-scaled phylogenies were subsequently used to perform ancestral state reconstructions of body mass and phylogenetic encephalization quotient (PEQ). Our divergence-time estimates are consistent with other recent studies that analysed extant species. We show that the origin of the genus Homo probably occurred between 4.30 and 2.56 million years ago. The ancestral state reconstructions show a general trend towards a smaller body mass before the emergence of Homo, followed by a trend towards a greater body mass. PEQ estimations display a general trend of gradual but accelerating encephalization evolution. The obtained results provide a rigorous temporal framework for human evolution.

Morphology and Bayesian tip-dating recover deep Cretaceous-age divergences among major chrysidid lineages (Hymenoptera: Chrysididae)

We integrated phylogenetic, biogeographic and palaeontological data to reconstruct the evolutionary history of the cuckoo wasps. We propose a phylogenetic hypothesis based on a comprehensive morphological study resulting in 300 characters coded for both living and extinct species. Phylogenetic relationships and divergence time estimation were simultaneously inferred in a Bayesian tip-dating framework, applying a relaxed morphological clock. Results unequivocally indicate Chrysididae to be monophyletic, as well as all traditionally recognized subfamilies and tribes. Within the Chrysidinae, Elampini was placed as the sister-group of the other three chrysidine tribes, with Parnopini as sister to the clade including Allocoeliini and Chrysidini. Dating analysis indicates that the major lineages started to differentiate around 130 Mya during the Early Cretaceous. The clades recognized as subfamilies started differentiating during the Palaeogene and the Neogene. Our results reveal an intricate process on the geographic evolution of chrysidid wasps and dispute previous ideas that Cretaceous-old splits in their early history could be associated with vicariant events related to the breakup between Africa and South America. The present-day southern disjunctions of some groups are interpreted as the outcome of more recent dispersals and extinctions of representatives from Nearctic and Palaearctic faunas during the Neogene, when northern continents became significantly colder.

Phylogeny and time scale of diversification in the fossil-rich sunfishes and black basses (Teleostei: Percomorpha: Centrarchidae)

Species of the North American freshwater fish lineage Centrarchidae are apex predators in their habitats and are among the world's most popular sport fishes. Centrarchids boast a rich fossil record that extends from the latest Eocene to the Pleistocene. To investigate the phylogeny and timing of diversification of Centrarchidae, we deploy a dataset of DNA sequences of 16 nuclear genes sampled from nearly all of the recognized and undescribed species. We also utilize previously published morphological datasets to assess the phylogenetic placement of one of the oldest known centrarchid fossils, †Plioplarchus whitei. A Bayesian multispecies coalescent species tree analysis provides insight on relationships that evaded resolution in earlier studies, such as the relationships of Acantharchus pomotis, the resolution of a clade consisting of species previously synonymized under the Spotted Bass, Micropterus punctulatus, and a clade of recently described species previously considered populations of the Redeye Bass, Micropterus coosae. This new molecular phylogeny and the inclusion of †P. whitei and other centrarchid fossils in the tip-dated fossilized birth-death analysis results in a new hypothesis of the timing of diversification in Centrarchidae that contextualizes the ages of centrarchid fossils to the timing of speciation among the extant species. In addition to providing new temporal perspectives on the diversification of freshwater fishes in North America, this study may close of the chapter of centrarchid phylogeny inferred using Sanger-sequenced genes, as the use of genomic-scale datasets becomes mainstream in the phylogenetics of fishes.

Redescription of the skull of the Australian flatback sea turtle, Natator depressus, provides new morphological evidence for phylogenetic relationships among sea turtles (Chelonioidea)

Chelonioidea (sea turtles) are a group where available morphological evidence for crown-group relationships are incongruent with those established using molecular data. However, morphological surveys of crown-group taxa tend to focus on a recurring subset of the extant species. The Australian flatback sea turtle, Natator depressus, is often excluded from comparisons and it is the most poorly known of the seven extant species of Chelonioidea. Previous descriptions of its skull morphology are limited and conflict. Here we describe three skulls of adult N. depressus and re-examine the phylogenetic relationships according to morphological character data. Using X-ray micro Computed Tomography we describe internal structures of the braincase and identify new phylogenetically informative characters not previously reported. Phylogenetic analysis using a Bayesian approach strongly supports a sister-group relationship between Chelonia mydas and N. depressus, a topology that was not supported by previous analyses of morphological data but one that matches the topology supported by analysis of molecular data. Our results highlight the general need to sample the morphological anatomy of crown-group taxa more thoroughly before concluding that morphological and molecular evidence are incongruous.

Mind the Outgroup and Bare Branches in Total-Evidence Dating: a Case Study of Pimpliform Darwin Wasps (Hymenoptera, Ichneumonidae)

Taxon sampling is a central aspect of phylogenetic study design, but it has received limited attention in the context of total-evidence dating, a widely used dating approach that directly integrates molecular and morphological information from extant and fossil taxa. We here assess the impact of commonly employed outgroup sampling schemes and missing morphological data in extant taxa on age estimates in a total-evidence dating analysis under the uniform tree prior. Our study group is Pimpliformes, a highly diverse, rapidly radiating group of parasitoid wasps of the family Ichneumonidae. We analyze a data set comprising 201 extant and 79 fossil taxa, including the oldest fossils of the family from the Early Cretaceous and the first unequivocal representatives of extant subfamilies from the mid-Paleogene. Based on newly compiled molecular data from ten nuclear genes and a morphological matrix that includes 222 characters, we show that age estimates become both older and less precise with the inclusion of more distant and more poorly sampled outgroups. These outgroups not only lack morphological and temporal information but also sit on long terminal branches and considerably increase the evolutionary rate heterogeneity. In addition, we discover an artifact that might be detrimental for total-evidence dating: "bare-branch attraction," namely high attachment probabilities of certain fossils to terminal branches for which morphological data are missing. Using computer simulations, we confirm the generality of this phenomenon and show that a large phylogenetic distance to any of the extant taxa, rather than just older age, increases the risk of a fossil being misplaced due to bare-branch attraction. After restricting outgroup sampling and adding morphological data for the previously attracting, bare branches, we recover a Jurassic origin for Pimpliformes and Ichneumonidae. This first age estimate for the group not only suggests an older origin than previously thought but also that diversification of the crown group happened well before the Cretaceous-Paleogene boundary. Our case study demonstrates that in order to obtain robust age estimates, total-evidence dating studies need to be based on a thorough and balanced sampling of both extant and fossil taxa, with the aim of minimizing evolutionary rate heterogeneity and missing morphological information. [Bare-branch attraction; ichneumonids; fossils; morphological matrix; phylogeny; RoguePlots.].

Investigating Morphological Complexes Using Informational Dissonance and Bayes Factors: A Case Study in Corbiculate Bees

It is widely recognized that different regions of a genome often have different evolutionary histories and that ignoring this variation when estimating phylogenies can be misleading. However, the extent to which this is also true for morphological data is still largely unknown. Discordance among morphological traits might plausibly arise due to either variable convergent selection pressures or else phenomena such as hemiplasy. Here, we investigate patterns of discordance among 282 morphological characters, which we scored for 50 bee species particularly targeting corbiculate bees, a group that includes the well-known eusocial honeybees and bumblebees. As a starting point for selecting the most meaningful partitions in the data, we grouped characters as morphological modules, highly integrated trait complexes that as a result of developmental constraints or coordinated selection we expect to share an evolutionary history and trajectory. In order to assess conflict and coherence across and within these morphological modules, we used recently developed approaches for computing Bayesian phylogenetic information allied with model comparisons using Bayes factors. We found that despite considerable conflict among morphological complexes, accounting for among-character and among-partition rate variation with individual gamma distributions, rate multipliers, and linked branch lengths can lead to coherent phylogenetic inference using morphological data. We suggest that evaluating information content and dissonance among partitions is a useful step in estimating phylogenies from morphological data, just as it is with molecular data. Furthermore, we argue that adopting emerging approaches for investigating dissonance in genomic datasets may provide new insights into the integration and evolution of anatomical complexes. [Apidae; entropy; morphological modules; phenotypic integration; phylogenetic information.].

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Raising names from the dead: A time-calibrated phylogeny of frog shells (Bursidae, Tonnoidea, Gastropoda) using mitogenomic data

With 59 Recent species, Bursidae, known as «frog shells», are a small but widely distributed group of tropical and subtropical gastropods that are most diverse in the Indo-West Pacific. The present study is aimed at reconstructing phylogenetic relationships of bursid gastropods based on extensive and representative taxon sampling. Five genetic markers (cytochrome c oxidase subunit I (cox1), 16 s and 12 s rRNA mitochondrial genes, 28 s rRNA and Histone H3 nuclear gene) were sequenced for over 30 species in every known genus but Crossata. Furthermore, we sequenced the complete mt-genome of 9 species (10 specimens) (Aspa marginata, Marsupina bufo, Korrigania quirihorai, Korrigania fijiensis, Tutufa rubeta, Bursa lamarckii, Lampasopsis rhodostoma (twice), Bufonaria perelegans and Bursa aff. tuberosissima). Our analysis recovered Bursidae as a monophyletic group, whereas the genus Bursa was found to be polyphyletic. The genera Talisman and Dulcerana are resurrected and the genera Alanbeuella gen. nov. and Korrigania gen. nov. are described. Dating analysis using 21 extinct taxa for node and simplified tip calibrations was performed, showing a diversification of the group in two phases. Diversification may be linked to tectonic events leading to biodiversity relocation from the western Tethys toward the Indo-Pacific.

Drimolen cranium DNH 155 documents microevolution in an early hominin species

Paranthropus robustus is a small-brained extinct hominin from South Africa characterized by derived, robust craniodental morphology. The most complete known skull of this species is DNH 7 from Drimolen Main Quarry, which differs from P. robustus specimens recovered elsewhere in ways attributed to sexual dimorphism. Here, we describe a new fossil specimen from Drimolen Main Quarry, dated from approximately 2.04-1.95 million years ago, that challenges this view. DNH 155 is a well-preserved adult male cranium that shares with DNH 7 a suite of primitive and derived features unlike those seen in adult P. robustus specimens from other chronologically younger deposits. This refutes existing hypotheses linking sexual dimorphism, ontogeny and social behaviour within this taxon, and clarifies hypotheses concerning hominin phylogeny. We document small-scale morphological changes in P. robustus associated with ecological change within a short time frame and restricted geography. This represents the most highly resolved evidence yet of microevolutionary change within an early hominin species.

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.

Morphological Evidence for an Oral Venom System in Caecilian Amphibians

Amphibians are known for their skin rich in glands containing toxins employed in passive chemical defense against predators, different from, for example, snakes that have active chemical defense, injecting their venom into the prey. Caecilians (Amphibia, Gymnophiona) are snake-shaped animals with fossorial habits, considered one of the least known vertebrate groups. We show here that amphibian caecilians, including species from the basal groups, besides having cutaneous poisonous glands as other amphibians do, possess specific glands at the base of the teeth that produce enzymes commonly found in venoms. Our analysis of the origin of these glands shows that they originate from the same tissue that gives rise to teeth, similar to the venom glands in reptiles. We speculate that caecilians might have independently developed mechanisms of production and injection of toxins early in their evolutionary history.

•

Amphibian caecilians have tooth-related glands in both upper and lower jaws

•

The glands have the same origin of reptile venom glands

•

The secretion contains proteins with enzymatic activities commonly found in venoms

•

Caecilians might have developed the ability to inject oral toxins early in evolution

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.

Key Steps in the Evolution of Mammalian Movement: A Prolegomenal Essay

Rich repertoires of movements underlie the complex social interactions of mammals. The building blocks, or syllables, of these movements are produced by spinal cord circuits that are comprised of diverse neuronal types that control musculoskeletal systems comprised of multi-segmented limbs. Together, these systems provide mammals with the evolutionary advantages of power, speed, and endurance. Here, I propose that the key steps in chordate evolution that led to these traits began with the development of the notochord and a proliferative ventricular zone (with associated Notch signalling). This step led to the production of diverse neuronal types that included the development of a sympathetic nervous system that could regulate the evolving cardiovascular system. And the sympathetic nervous system in turn led to the development of homeothermic endothermy, a requirement for motor systems to produce a combination of power, speed, and endurance. Furthermore, the evolution of the continuous structure of the spinal cord led not only to a structure fit for cartesian signalling molecules, but also to one with high processing power in which circuits for effecting movement syllables formed. These syllables are harnessed by higher regions of nervous systems to produce the complex movements required for interactions with others and with the surrounding environment.