Rates of phenotypic and genomic evolution during the Cambrian explosion

Total-evidence phylogenetic analysis resolves the evolutionary timescale of mantis shrimps (Stomatopoda) and provides insights into their molecular and morphological evolutionary rates

The crustacean order Stomatopoda comprises approximately 500 species of mantis shrimps. These marine predators, common in tropical and subtropical waters, possess sophisticated visual systems and specialized hunting appendages. In this study, we infer the evolutionary relationships within Stomatopoda using a combined data set of 77 morphological characters, whole mitochondrial genomes, and three nuclear markers. Our data set includes representatives from all seven stomatopod superfamilies, including the first sequence data from Erythrosquilloidea. Using a Bayesian relaxed molecular clock with fossil-based calibration priors, we estimate that crown-group unipeltatan stomatopods appeared ∼ 143 (95 % credible interval 199-98) million years ago in the Mesozoic. Additionally, our results support the hypothesis that specialized smashing and spearing appendages appeared early in the evolutionary history of Unipeltata. We found no evidence of a correlation between rates of morphological and molecular evolution across the phylogeny, but identified very high levels of among-lineage rate variation in the morphological characters. Our total-evidence analysis recovered evolutionary signals from both molecular and morphological data sets, demonstrating the merit in combining these sources of information for phylogenetic inference and evolutionary analysis.

Copepod phylogenomics supports Canuelloida as a valid order separate from Harpacticoida

Copepods are small crustaceans that are ubiquitous in aquatic environments. They are particularly abundant in marine and freshwater plankton, marine sediments, and as parasites or commensals of other aquatic organisms. Despite their abundance and importance, phylogenetic relationships among copepods are poorly resolved. The validity of higher-level taxa, including several orders, has continued to be controversial throughout the 21st century. This study has two main goals: first, to use phylogenomic data to assess relationships among the four major copepod orders: Calanoida, Cyclopoida, Harpacticoida, and Siphonostomatoida, which together include more than 98 % of copepod species diversity, and second, to test the validity of the recently proposed order Canuelloida. Towards these goals, we sampled 28 copepod transcriptomes and genomes spanning 20 families and 5 orders, including the first transcriptome of a representative of Canuelloida. We identified 2,527 single copy protein coding genes comprising 939,460 amino acid (aa) positions and 530,269 informative sites. All phylogenetic analyses support a monophyletic Podoplea (i.e., the superorder comprising all copepod orders except for Calanoida and Platycopioida) with Calanoida as its sister taxon. We find robust support across all methods for Canuelloida as a distinct order separate from the traditionally recognized Harpacticoida (Oligoarthra). Contrary to several recent studies of smaller sets of nuclear genes or mitochondrial genomes, we recover Cyclopoida and Harpacticoida as sister taxa and find that gene tree discordance analysis rejects the alternative topologies. Transcriptomic data are promising for resolving the backbone of the copepod phylogeny but collecting and sequencing the nearly 15,000 species of copepods, many of which are infrequently encountered and less than 1 mm in size, remains a major hurdle.

Deciphering the biosynthetic potential of microbial genomes using a BGC language processing neural network model

Biosynthetic gene clusters (BGCs), key in synthesizing microbial secondary metabolites, are mostly hidden in microbial genomes and metagenomes. To unearth this vast potential, we present BGC-Prophet, a transformer-based language model for BGC prediction and classification. Leveraging the transformer encoder, BGC-Prophet captures location-dependent relationships between genes. As one of the pioneering ultrahigh-throughput tools, BGC-Prophet significantly surpasses existing methods in efficiency and fidelity, enabling comprehensive pan-phylogenetic and whole-metagenome BGC screening. Through the analysis of 85 203 genomes and 9428 metagenomes, BGC-Prophet has profiled an extensive array of sub-million BGCs. It highlights notable enrichment in phyla like Actinomycetota and the widespread distribution of polyketide, NRP, and RiPP BGCs across diverse lineages. It reveals enrichment patterns of BGCs following important geological events, suggesting environmental influences on BGC evolution. BGC-Prophet's capabilities in detection of BGCs and evolutionary patterns offer contributions to deeper understanding of microbial secondary metabolites and application in synthetic biology.

Exotic cuticular specializations in a Cambrian scalidophoran

Scalidophora, the ecdysozoan group including priapulids, kinorhynchs and loriciferans, comprises some of the most abundant and ecologically important Cambrian animals. However, reconstructions of the morphology and lifestyles of fossil scalidophorans are often hampered by poor preservation of their submillimetre-scale cuticular specializations. Based on exceptionally preserved small carbonaceous fossils (SCFs), we describe a new scalidophoran-grade animal, Scalidodendron crypticum gen. et sp. nov., from the Early to Middle Cambrian Hess River Formation of northern Canada. The Hess River SCFs comprise pharyngeal teeth, coniform sclerites and hook-like sclerites, all closely comparable to known scalidophoran counterparts. The coniform and hook-like sclerites recurrently associate with arborescent cuticular projections that show multiple orders of branching, morphologically unlike those of any known living or fossil scalidophoran. The fine splintering and inferred post-pharyngeal position of these structures argue against locomotory, feeding and defensive roles with direct analogues in extant counterparts. As such, the arborescent structures of Scalidodendron denote a previously cryptic range of morphological variation in Cambrian scalidophorans, paralleling that of coeval panarthropods but expressed at a fundamentally different level of anatomical organization.

Impact of ontogeny and spines on the hydrodynamic performance of the Cambrian arthropod Isoxys

A metazoan-dominated biological pump was established early in the Phanerozoic, a time that saw the evolution of the first pelagic euarthropod zooplankton such as some species of the Cambrian bivalved euarthropod Isoxys. Pelagic groups evolved from benthic stock, in many cases through neoteny and retention of characteristics from planktic larval stages. However, Isoxys brooded eggs and did not have a planktic larval stage, precluding this route into the pelagic realm. Computational fluid dynamics was used to quantify hydrodynamic performance through the ontogeny of two hyperbenthic species of Isoxys, Isoxys auritus and Isoxys minor. Coefficients were used to quantify forces for different carapace shapes over a range of biologically relevant sizes and swimming speeds. Streamlining and lift generation were greater for later growth stages, a consequence of carapace asymmetry and elongated anterior and posterior spines. Simulations performed with anterior spines artificially removed demonstrate the importance of this feature for lift generation, with a minimal impact on drag. Elongated spines and faster swimming can also be considered anti-predatory, and the reduction of drag would have reduced the detectability of Isoxys to predators. Taken together, it is likely that pelagic Isoxys species evolved from benthic stock through the co-option of anti-predatory features.

Ediacaran origin and Ediacaran-Cambrian diversification of Metazoa

The timescale of animal diversification has been a focus of debate over how evolutionary history should be calibrated to geologic time. Molecular clock analyses have invariably estimated a Cryogenian or Tonian origin of animals while unequivocal animal fossils first occur in the Ediacaran. However, redating of key Ediacaran biotas and the discovery of several Ediacaran crown-Metazoa prompt recalibration of molecular clock analyses. We present revised fossil calibrations and use them in molecular clock analyses estimating the timescale of metazoan evolutionary history. Integrating across uncertainties including phylogenetic relationships, clock model, and calibration strategy, we estimate Metazoa to have originated in the early Ediacaran, Eumetazoa in the middle Ediacaran, and Bilateria in the upper Ediacaran, with many crown-phyla originating across the Ediacaran-Cambrian interval or elsewise fully within the Cambrian. These results are in much closer accord with the fossil record, coinciding with marine oxygenation, but they reject a literal reading of the fossil record.

Stem Life: A Framework for Understanding the Prebiotic-Biotic Transition

Abiogenesis is frequently envisioned as a linear, ladder-like progression of increasingly complex chemical systems, eventually leading to the ancestors of extant cellular life. This "pre-cladistics" view is in stark contrast to the well-accepted principles of organismal evolutionary biology, as informed by paleontology and phylogenetics. Applying this perspective to origins, I explore the paradigm of "Stem Life," which embeds abiogenesis within a broader continuity of diversification and extinction of both hereditary lineages and chemical systems. In this new paradigm, extant life's ancestral lineage emerged alongside and was dependent upon many other complex prebiotic chemical systems, as part of a diverse and fecund prebiosphere. Drawing from several natural history analogies, I show how this shift in perspective enriches our understanding of Origins and directly informs debates on defining Life, the emergence of the Last Universal Common Ancestor (LUCA), and the implications of prebiotic chemical experiments.

Arthropod Phylotranscriptomics With a Special Focus on the Basal Phylogeny of the Myriapoda

Arthropoda represents the most diverse animal phylum, but clarifying the phylogenetic relationships among arthropod taxa remains challenging given the numerous arthropod lineages that diverged over a short period of time. In order to resolve the most controversial aspects of deep arthropod phylogeny, focusing on the Myriapoda, we conducted phylogenetic analyses based on ten super-matrices comprised of 751 to 1,233 orthologous genes across 64 representative arthropod species, including 28 transcriptomes that were newly generated in this study. Our findings provide unambiguous support for the monophyly of the higher arthropod taxa, Chelicerata, Mandibulata, Myriapoda, Pancrustacea, and Hexapoda, while the Crustacea are paraphyletic, with the class Remipedia supported as the lineage most closely related to hexapods. Within the Hexapoda, our results largely affirm previously proposed phylogenetic relationships among deep hexapod lineages, except that the Paraneoptera (Hemiptera, Thysanoptera, and Psocodea) was recovered as a monophyletic lineage in some analyses. The results corroborated the recently proposed phylogenetic framework of the four myriapod classes, wherein Symphyla and Pauropoda, as well as Chilopoda and Diplopoda, are each proposed to be sister taxa. The findings provide important insights into understanding the phylogeny and evolution of arthropods.

Two Notorious Nodes: A Critical Examination of Relaxed Molecular Clock Age Estimates of the Bilaterian Animals and Placental Mammals

The popularity of relaxed clock Bayesian inference of clade origin timings has generated several recent publications with focal results considerably older than the fossils of the clades in question. Here, we critically examine two such clades: the animals (with a focus on the bilaterians) and the mammals (with a focus on the placentals). Each example displays a set of characteristic pathologies which, although much commented on, are rarely corrected for. We conclude that in neither case does the molecular clock analysis provide any evidence for an origin of the clade deeper than what is suggested by the fossil record. In addition, both these clades have other features (including, in the case of the placental mammals, proximity to a large mass extinction) that allow us to generate precise expectations of the timings of their origins. Thus, in these instances, the fossil record can provide a powerful test of molecular clock methodology, and why it goes astray, and we have every reason to think these problems are general. [Cambrian explosion; mammalian evolution; molecular clocks.].

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24β1, -A24β2, -A34α2, -A34β1, and -A34β2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded β-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34β1/β2 had an additional two-stranded β-sheet. We hypothesize that this second β-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

Profile of a flower: How rates of morphological evolution drive floral diversification in Ericales and angiosperms

PREMISE

Recent studies of floral disparity in the asterid order Ericales have shown that flowers vary strongly among families and that disparity is unequally distributed between the three flower modules (perianth, androecium, gynoecium). However, it remains unknown whether these patterns are driven by heterogeneous rates of morphological evolution or other factors.

METHODS

Here, we compiled a data set of 33 floral characters scored for 414 species of Ericales sampled from 346 genera and all 22 families. We conducted ancestral state reconstructions using an equal-rates Markov model for each character. We estimated rates of morphological evolution for Ericales and for a separate angiosperm-wide data set of 19 characters and 792 species, creating "rate profiles" for Ericales, angiosperms, and major angiosperm subclades. We compared morphological rates among flower modules within each data set separately and between data sets, and we compared rates among angiosperm subclades using the angiosperm data set.

RESULTS

The androecium exhibits the highest evolutionary rates across most characters, whereas most perianth and gynoecium characters evolve more slowly in both Ericales and angiosperms. Both high and low rates of morphological evolution can result in high floral disparity in Ericales. Analyses of an angiosperm-wide floral data set reveal that this pattern appears to be conserved across most major angiosperm clades.

CONCLUSIONS

Elevated rates of morphological evolution in the androecium of Ericales may explain the higher disparity reported for this floral module. Comparing rates of morphological evolution through rate profiles proves to be a powerful tool in understanding floral evolution.

Expansion and collapse of VEGF diversity in major clades of the animal kingdom

Together with the platelet-derived growth factors (PDGFs), the vascular endothelial growth factors (VEGFs) form the PDGF/VEGF subgroup among cystine knot growth factors. The evolutionary relationships within this subgroup have not been examined thoroughly to date. Here, we comprehensively analyze the PDGF/VEGF growth factors throughout all animal phyla and propose a phylogenetic tree. Vertebrate whole-genome duplications play a role in expanding PDGF/VEGF diversity, but several limited duplications are necessary to account for the temporal pattern of emergence. The phylogenetically oldest PDGF/VEGF-like growth factor likely featured a C-terminus with a BR3P signature, a hallmark of the modern-day lymphangiogenic growth factors VEGF-C and VEGF-D. Some younger VEGF genes, such as VEGFB and PGF, appeared completely absent in important vertebrate clades such as birds and amphibia, respectively. In contrast, individual PDGF/VEGF gene duplications frequently occurred in fish on top of the known fish-specific whole-genome duplications. The lack of precise counterparts for human genes poses limitations but also offers opportunities for research using organisms that diverge considerably from humans. Sources for the graphical abstract: 326 MYA and older [1]; 72-240 MYA [2]; 235-65 MYA [3].

A corazonin G protein-coupled receptor gene in the tick Ixodes scapularis yields two splice variants, each coding for a specific corazonin receptor

We have identified a corazonin G protein-coupled receptor (GPCR) gene in the tick Ixodes scapularis, which likely plays a central role in the physiology and behavior of this ectoparasite. This receptor gene is unusually large (1.133 Mb) and yields two corazonin (CRZ) receptor splice variants, where nearly half of the coding regions are exchanged: CRZ-Ra (containing exon 2, exon 3, and exon 4 of the gene) and CRZ-Rb (containing exon 1, exon 3, and exon 4 of the gene). CRZ-Ra codes for a GPCR with a canonical DRF sequence at the border of the third transmembrane helix and the second intracellular loop. The positively-charged R residue from the DRF sequence is important for coupling of G proteins after activation of a GPCR. CRZ-Rb, in contrast, codes for a GPCR with an unusual DQL sequence at this position, still retaining a negatively-charged D residue, but lacking a positively-charged R residue, suggesting different G protein coupling. Another difference between the two splice variants is that exon 2 from CRZ-Ra codes for an N-terminal signal sequence. Normally, GPCRs do not have N-terminal signal sequences, although a few mammalian GPCRs have. In the tick CRZ-Ra, the signal sequence probably assists with inserting the receptor correctly into the RER membrane. We stably transfected Chinese Hamster Ovary cells with each of the two splice variants and carried out bioluminescence bioassays that also included the use of the human promiscuous G protein G₁₆. CRZ-Ra turned out to be selective for I. scapularis corazonin (EC₅₀ = 10^-8 M) and could not be activated by related neuropeptides like adipokinetic hormone (AKH) and AKH/corazonin-related peptide (ACP). Similarly, also CRZ-Rb could only be activated by corazonin, although about 4-fold higher concentrations were needed to activate it (EC₅₀ = 4 x 10^-8 M). The genomic organization of the tick corazonin GPCR gene is similar to that of the insect AKH and ACP receptor genes. This similar genomic organization can also be found in the human gonadotropin-releasing hormone (GnRH) receptor gene, confirming previous conclusions that the corazonin, AKH, and ACP receptor genes are the true arthropod orthologues of the human GnRH receptor gene.

Ultraconserved elements-based phylogenomic systematics of the snake superfamily Elapoidea, with the description of a new Afro-Asian family

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes.

Reassessing a cryptic history of early trilobite evolution

Trilobites are an iconic Paleozoic group of biomineralizing marine euarthropods that appear abruptly in the fossil record (c. 521 million years ago) during the Cambrian ‘explosion’. This sudden appearance has proven controversial ever since Darwin puzzled over the lack of pre-trilobitic fossils in the Origin of Species, and it has generally been assumed that trilobites must have an unobserved cryptic evolutionary history reaching back into the Precambrian. Here we review the assumptions behind this model, and suggest that a cryptic history creates significant difficulties, including the invocation of rampant convergent evolution of biomineralized structures and the abandonment of the synapomorphies uniting the clade. We show that a vicariance explanation for early Cambrian trilobite palaeobiogeographic patterns is inconsistent with factors controlling extant marine invertebrate distributions, including the increasingly-recognized importance of long-distance dispersal. We suggest that survivorship bias may explain the initial rapid diversification of trilobites, and conclude that the group’s appearance at c. 521 Ma closely reflects their evolutionary origins.

A reassessment of early trilobite phylogenetic relationships and palaeobiogeographic patterns suggests that a cryptic evolutionary history is unlikely for this group. The abrupt appearance of trilobites is likely to closely reflect their evolutionary origins, and may be explained by survivorship biases inherent in the fossil record.

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.

The origin and early evolution of arthropods

The rise of arthropods is a decisive event in the history of life. Likely the first animals to have established themselves on land and in the air, arthropods have pervaded nearly all ecosystems and have become pillars of the planet's ecological networks. Forerunners of this saga, exceptionally well-preserved Palaeozoic fossils recently discovered or re-discovered using new approaches and techniques have elucidated the precocious appearance of extant lineages at the onset of the Cambrian explosion, and pointed to the critical role of the plankton and hard integuments in early arthropod diversification. The notion put forward at the beginning of the century that the acquisition of extant arthropod characters was stepwise and represented by the majority of Cambrian fossil taxa is being rewritten. Although some key traits leading to Euarthropoda are indeed well documented along a diversified phylogenetic stem, this stem led to several speciose and ecologically diverse radiations leaving descendants late into the Palaeozoic, and a large part, if not all of the Cambrian euarthropods can now be placed on either of the two extant lineages: Mandibulata and Chelicerata. These new observations and discoveries have altered our view on the nature and timing of the Cambrian explosion and clarified diagnostic characters at the origin of extant arthropods, but also raised new questions, especially with respect to cephalic plasticity. There is now strong evidence that early arthropods shared a homologous frontalmost appendage, coined here the cheira, which likely evolved into antennules and chelicerae, but other aspects, such as brain and labrum evolution, are still subject to active debate. The early evolution of panarthropods was generally driven by increased mastication and predation efficiency and sophistication, but a wealth of recent studies have also highlighted the prevalent role of suspension-feeding, for which early panarthropods developed their own adaptive feedback through both specialized appendages and the diversification of small, morphologically differentiated larvae. In a context of general integumental differentiation and hardening across Cambrian metazoans, arthrodization of body and limbs notably prompted two diverging strategies of basipod differentiation, which arguably became founding criteria in the divergence of total-groups Mandibulata and Chelicerata. The kinship of trilobites and their relatives remains a source of disagreement, but a recent topological solution, termed the 'deep split', could embed Artiopoda as sister taxa to chelicerates and constitute definitive support for Arachnomorpha. Although Cambrian fossils have been critical to all these findings, data of exceptional quality have also been accumulating from other Palaeozoic Konservat-Lagerstätten, and a better integration of this information promises a much more complete and elaborate picture of early arthropod evolution in the near future. From the broader perspective of a total-evidence approach to the understanding of life's history, and despite persisting systematic debates and new interpretative challenges, various advances based on palaeontological evidence open the prospect of finally using the full potential of the most diverse animal phylum to investigate macroevolutionary patterns and processes.

Phylogenomic Analysis of Velvet Worms (Onychophora) Uncovers an Evolutionary Radiation in the Neotropics

Onychophora ("velvet worms") are charismatic soil invertebrates known for their status as a "living fossil," their phylogenetic affiliation to arthropods, and their distinctive biogeographic patterns. However, several aspects of their internal phylogenetic relationships remain unresolved, limiting our understanding of the group's evolutionary history, particularly with regard to changes in reproductive mode and dispersal ability. To address these gaps, we used RNA sequencing and phylogenomic analysis of transcriptomes to reconstruct the evolutionary relationships and infer divergence times within the phylum. We recovered a fully resolved and well-supported phylogeny for the circum-Antarctic family Peripatopsidae, which retains signals of Gondwanan vicariance and showcases the evolutionary lability of reproductive mode in the family. Within the Neotropical clade of Peripatidae, though, we found that amino acid-translated sequence data masked nearly all phylogenetic signal, resulting in highly unstable and poorly supported relationships. Analyses using nucleotide sequence data were able to resolve many more relationships, though we still saw discordant phylogenetic signal between genes, probably indicative of a rapid, mid-Cretaceous radiation in the group. Finally, we hypothesize that the unique reproductive mode of placentotrophic viviparity found in all Neotropical peripatids may have facilitated the multiple inferred instances of over-water dispersal and establishment on oceanic islands.

The Continuity Principle and the Evolution of Replication Fidelity

Evolution in modern life requires high replication fidelity to allow for natural selection. A simulation model utilizing simulated phenotype data on cellular probability of survival was developed to determine how self-replication fidelity could evolve in early life. The results indicate that initial survivability and replication fidelity both contribute to overall fitness as measured by growth rates of the cell population. Survival probability was the more dominant feature, and evolution was possible even with zero replication fidelity. A derived formula for the relationship of survival probability and replication fidelity with growth rate was consistent with the simulated empirical data. Quantitative assessment of continuity and other evidence was obtained for a saltation (non-continuous) evolutionary process starting from low to moderate levels of survival probability and self-replication fidelity to reach the high levels seen in modern life forms.

Multiple paths to morphological diversification during the origin of amniotes

Early terrestrial vertebrates (amniotes) provide a classic example of diversification following adaptive zone invasion. The initial terrestrialization of vertebrates was closely followed by dietary diversification, as evidenced by a proliferation of craniomandibular and dental adaptations. However, morphological evolution of early amniotes has received limited study, in analyses with restricted taxonomic scope, leaving substantial questions about the dynamics of this important terrestrial radiation. We use novel analyses of discrete characters to quantify variation in evolutionary rates and constraints during diversification of the amniote feeding apparatus. We find evidence for an early burst, comprising high rates of anatomical change that decelerated through time, giving way to a background of saturated morphological evolution. Subsequent expansions of phenotypic diversity were not associated with increased evolutionary rates. Instead, variation in the mode of evolution became important, with groups representing independent origins of herbivory evolving distinctive, group-specific morphologies and thereby exploring novel character-state spaces. Our findings indicate the importance of plant-animal interactions in structuring the earliest radiation of amniotes and demonstrate the importance of variation in modes of phenotypic divergence during a major evolutionary radiation.

Respiratory Heme A-Containing Oxidases Originated in the Ancestors of Iron-Oxidizing Bacteria

Respiration is a major trait shaping the biology of many environments. Cytochrome oxidase containing heme A (COX) is a common terminal oxidase in aerobic bacteria and is the only one in mammalian mitochondria. The synthesis of heme A is catalyzed by heme A synthase (CtaA/Cox15), an enzyme that most likely coevolved with COX. The evolutionary origin of COX in bacteria has remained unknown. Using extensive sequence and phylogenetic analysis, we show that the ancestral type of heme A synthases is present in iron-oxidizing Proteobacteria such as Acidithiobacillus spp. These bacteria also contain a deep branching form of the major COX subunit (COX1) and an ancestral variant of CtaG, a protein that is specifically required for COX biogenesis. Our work thus suggests that the ancestors of extant iron-oxidizers were the first to evolve COX. Consistent with this conclusion, acidophilic iron-oxidizing prokaryotes lived on emerged land around the time for which there is the earliest geochemical evidence of aerobic respiration on earth. Hence, ecological niches of iron oxidation have apparently promoted the evolution of aerobic respiration.

Diversification dynamics of total-, stem-, and crown-groups are compatible with molecular clock estimates of divergence times

Molecular evolutionary time scales are expected to predate the fossil evidence, but, particularly for major evolutionary radiations, they can imply extremely protracted stem lineages predating the origin of living clades, leading to claims of systematic overestimation of divergence times. We use macroevolutionary birth-death models to describe the range of total-group and crown-group ages expected under constant rates of speciation and extinction. We extend current predictions on origination times for crown- and total-groups, and extinction of stem-groups, demonstrating that there is broad variance in these predictions. Under constant rates of speciation and extinction, we show that the distribution of expected arthropod total-group ages is consistent with molecular clock estimates. The fossil record cannot be read literally, and our results preclude attempts to interpret the antiquity of clades based on the co-occurrence of stem- and crown-representatives.

On the Evolution of the Biological Framework for Insight

The details of abiogenesis, to date, remain a matter of debate and constitute a key mystery in science and philosophy. The prevailing scientific hypothesis implies an evolutionary process of increasing complexity on Earth starting from (self-) replicating polymers. Defining the cut-off point where life begins is another moot point beyond the scope of this article. We will instead walk through the known evolutionary steps that led from these first exceptional polymers to the vast network of living biomatter that spans our world today, focusing in particular on perception, from simple biological feedback mechanisms to the complexity that allows for abstract thought. We will then project from the well-known to the unknown to gain a glimpse into what the universe aims to accomplish with living matter, just to find that if the universe had ever planned to be comprehended, evolution still has a long way to go.

Fossils improve phylogenetic analyses of morphological characters

Fossils provide our only direct window into evolutionary events in the distant past. Incorporating them into phylogenetic hypotheses of living clades can help time-calibrate divergences, as well as elucidate macroevolutionary dynamics. However, the effect fossils have on phylogenetic reconstruction from morphology remains controversial. The consequences of explicitly incorporating the stratigraphic ages of fossils using tip-dated inference are also unclear. Here, we use simulations to evaluate the performance of inference methods across different levels of fossil sampling and missing data. Our results show that fossil taxa improve phylogenetic analysis of morphological datasets, even when highly fragmentary. Irrespective of inference method, fossils improve the accuracy of phylogenies and increase the number of resolved nodes. They also induce the collapse of ancient and highly uncertain relationships that tend to be incorrectly resolved when sampling only extant taxa. Furthermore, tip-dated analyses under the fossilized birth-death process outperform undated methods of inference, demonstrating that the stratigraphic ages of fossils contain vital phylogenetic information. Fossils help to extract true phylogenetic signals from morphology, an effect that is mediated by both their distinctive morphology and their temporal information, and their incorporation in total-evidence phylogenetics is necessary to faithfully reconstruct evolutionary history.

The origin and evolution of the euarthropod labrum

A widely (although not universally) accepted model of arthropod head evolution postulates that the labrum, a structure seen in almost all living euarthropods, evolved from an anterior pair of appendages homologous to the frontal appendages of onychophorans. However, the implications of this model for the interpretation of fossil arthropods have not been fully integrated into reconstructions of the euarthropod stem group, which remains in a state of some disorder. Here I review the evidence for the nature and evolution of the labrum from living taxa, and reconsider how fossils should be interpreted in the light of this. Identification of the segmental identity of head appendage in fossil arthropods remains problematic, and often rests ultimately on unproven assertions. New evidence from the Cambrian stem-group euarthropod Parapeytoia is presented to suggest that an originally protocerebral appendage persisted well up into the upper stem-group of the euarthropods, which prompts a re-evaluation of widely-accepted segmental homologies and the interpretation of fossil central nervous systems. Only a protocerebral brain was implicitly present in a large part of the euarthropod stem group, and the deutocerebrum must have been a relatively late addition.

Deciphering trophic interactions in a mid-Cambrian assemblage

Exceptionally preserved fossil sites have allowed specimen-based identification of trophic interactions to which network analyses have been applied. However, network analyses of the fossil record suffer from incomplete and indirect data, time averaging that obscures species coexistence, and biases in preservation. Here, we present a high-resolution fossil data set from Raymond Quarry member of the mid-Cambrian Burgess Shale (7,549 specimens, 61 taxa, ∼510 Mya) and formulate a measure of "preservation bias" that aids identification of assemblage subsets to which network analyses can be reliably applied. For these sections, abundance correlation network analyses predicted longitudinally consistent trophic and competitive interactions. Our analyses predicted previously postulated trophic interactions with 83.5% accuracy and demonstrated a shift from specialist interaction-dominated assemblages to ones dominated by generalist and competitive interactions. This approach provides a robust, taphonomically corrected framework to explore and predict in detail the existence and ecological character of putative interactions in fossil data sets.

An interplay of resource availability, population size and mutation rate potentiates the evolution of metabolic signaling

Background

Asexually reproducing populations of single cells evolve through mutation, natural selection, and genetic drift. Environmental conditions in which the evolution takes place define the emergent fitness landscapes. In this work, we used Avida—a digital evolution framework—to uncover a hitherto unexplored interaction between mutation rates, population size, and the relative abundance of metabolizable resources, and its effect on evolutionary outcomes in small populations of digital organisms.

Results

Over each simulation, the population evolved to one of several states, each associated with a single dominant phenotype with its associated fitness and genotype. For a low mutation rate, acquisition of fitness by organisms was accompanied with, and dependent on, an increase in rate of genomic replication. At an increased mutation rate, phenotypes with high fitness values were similarly achieved through enhanced genome replication rates. In addition, we also observed the frequent emergence of suboptimal fitness phenotype, wherein neighboring organisms signaled to each other information relevant to performing metabolic tasks. This metabolic signaling was vital to fitness acquisition and was correlated with greater genotypic and phenotypic heterogeneity in the population. The frequency of appearance of signaling populations increased with population size and with resource abundance.

Conclusions

Our results reveal a minimal set of environment–genotype interactions that lead to the emergence of metabolic signaling within evolving populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01782-0.

Morphological Phylogenetics Evaluated Using Novel Evolutionary Simulations

Evolutionary inferences require reliable phylogenies. Morphological data have traditionally been analyzed using maximum parsimony, but recent simulation studies have suggested that Bayesian analyses yield more accurate trees. This debate is ongoing, in part, because of ambiguity over modes of morphological evolution and a lack of appropriate models. Here, we investigate phylogenetic methods using two novel simulation models-one in which morphological characters evolve stochastically along lineages and another in which individuals undergo selection. Both models generate character data and lineage splitting simultaneously: the resulting trees are an emergent property, rather than a fixed parameter. Standard consensus methods for Bayesian searches (Mki) yield fewer incorrect nodes and quartets than the standard consensus trees recovered using equal weighting and implied weighting parsimony searches. Distances between the pool of derived trees (most parsimonious or posterior distribution) and the true trees-measured using Robinson-Foulds (RF), subtree prune and regraft (SPR), and tree bisection reconnection (TBR) metrics-demonstrate that this is related to the search strategy and consensus method of each technique. The amount and structure of homoplasy in character data differ between models. Morphological coherence, which has previously not been considered in this context, proves to be a more important factor for phylogenetic accuracy than homoplasy. Selection-based models exhibit relatively lower homoplasy, lower morphological coherence, and higher inaccuracy in inferred trees. Selection is a dominant driver of morphological evolution, but we demonstrate that it has a confounding effect on numerous character properties which are fundamental to phylogenetic inference. We suggest that the current debate should move beyond considerations of parsimony versus Bayesian, toward identifying modes of morphological evolution and using these to build models for probabilistic search methods. [Bayesian; evolution; morphology; parsimony; phylogenetics; selection; simulation.].

Survival and selection biases in early animal evolution and a source of systematic overestimation in molecular clocks

Important evolutionary events such as the Cambrian Explosion have inspired many attempts at explanation: why do they happen when they do? What shapes them, and why do they eventually come to an end? However, much less attention has been paid to the idea of a 'null hypothesis'-that certain features of such diversifications arise simply through their statistical structure. Such statistical features also appear to influence our perception of the timing of these events. Here, we show in particular that study of unusually large clades leads to systematic overestimates of clade ages from some types of molecular clocks, and that the size of this effect may be enough to account for the puzzling mismatches seen between these molecular clocks and the fossil record. Our analysis of the fossil record of the late Ediacaran to Cambrian suggests that it is likely to be recording a true evolutionary radiation of the bilaterians at this time, and that explanations involving various sorts of cryptic origins for the bilaterians do not seem to be necessary.

Megaevolutionary dynamics and the timing of evolutionary innovation in reptiles

The origin of phenotypic diversity among higher clades is one of the most fundamental topics in evolutionary biology. However, due to methodological challenges, few studies have assessed rates of evolution and phenotypic disparity across broad scales of time to understand the evolutionary dynamics behind the origin and early evolution of new clades. Here, we provide a total-evidence dating approach to this problem in diapsid reptiles. We find major chronological gaps between periods of high evolutionary rates (phenotypic and molecular) and expansion in phenotypic disparity in reptile evolution. Importantly, many instances of accelerated phenotypic evolution are detected at the origin of major clades and body plans, but not concurrent with previously proposed periods of adaptive radiation. Furthermore, strongly heterogenic rates of evolution mark the acquisition of similarly adapted functional types, and the origin of snakes is marked by the highest rates of phenotypic evolution in diapsid history.

Gene content evolution in the arthropods

BACKGROUND

Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods.

RESULTS

Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception.

CONCLUSIONS

These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Homoplasy-Based Partitioning Outperforms Alternatives in Bayesian Analysis of Discrete Morphological Data

Bayesian analysis of morphological data is becoming increasingly popular mainly (but not only) because it allows for time-calibrated phylogenetic inference using relaxed morphological clocks and tip dating whenever fossils are available. As with molecular data, recent studies have shown that modeling among-character rate variation (ACRV) in morphological matrices greatly improves phylogenetic inference. In a likelihood framework this may be accomplished, for instance, by employing a hidden Markov model to assign characters to rate categories drawn from a (discretized) $\Gamma$ distribution and/or by partitioning data sets according to rate heterogeneity and estimating per-partition branch lengths, conditioned on a single topology. While the first approach is available in many phylogenetic analysis software, there is still no clear consensus on how to partition data, except perhaps in the simplest cases (e.g., "by codon" partitioning of coding sequences). Additionally, there is a trade-off between improvement in likelihood scores and the number of free parameters in the analysis, which rises quickly with the number of partitions. This trade-off may be dealt with by employing statistics that penalize overfitting of complex models, such as Akaike or Bayesian information criteria, or the more recently introduced stepping-stone method for marginal likelihood approximation. We applied the latter to three distinct matrices of discrete morphological data and demonstrated that sorting characters by homoplasy scores (obtained from implied weighting parsimony analysis) outperformed other partitioning strategies (anatomically-based and PartitionFinder2). The method was in fact so efficient in segregating characters by rates of evolution that no within-partition ACRV modeling was necessary, while among-partition rate variation was adequately accommodated by rate multipliers. We conclude that partitioning by homoplasy is a powerful and easy-to-implement strategy to address ACRV in complex data sets. We provide some guidelines focusing on morphological matrices, although this approach may be also applicable to molecular data sets.

Evo-Devo Algorithms: Gene-Regulation for Digital Architecture

The majority of current visual-algorithmic architecture is constricted to specific parameters that are gradient related, keeping their parts’ relation fixed within the algorithm, far away from a truly parametric modeling with a flexible topology. Recent findings around genetics and certain genes capable of shape conditioning (development) have succeeded in recovering the science of embryology as a valid field that connects and affects the evolutionary ecosystem, showing the existence of universal mechanisms that are present in living species, thus describing powerful strategies for generation and emergence. Therefore, a new dual discipline is justified: Evolutionary developmental biology science. Authors propose the convergence of genetics algorithms and simulated features from evolutionary developmental biology into a single data-flow that will prove itself capable of generating great diversity through a simple and flexible structure of data, commands, and polygonal geometry. For that matter, a case study through visual-algorithmic software deals with the hypothesis that for obtaining a greater emergence and design space, a simpler and more flexible approach might only be required, prioritizing hierarchical levels over complex and detailed operations.

The Temporal and Environmental Context of Early Animal Evolution: Considering All the Ingredients of an "Explosion"

Animals originated and evolved during a unique time in Earth history-the Neoproterozoic Era. This paper aims to discuss (1) when landmark events in early animal evolution occurred, and (2) the environmental context of these evolutionary milestones, and how such factors may have affected ecosystems and body plans. With respect to timing, molecular clock studies-utilizing a diversity of methodologies-agree that animal multicellularity had arisen by ∼800 million years ago (Ma) (Tonian period), the bilaterian body plan by ∼650 Ma (Cryogenian), and divergences between sister phyla occurred ∼560-540 Ma (late Ediacaran). Most purported Tonian and Cryogenian animal body fossils are unlikely to be correctly identified, but independent support for the presence of pre-Ediacaran animals is recorded by organic geochemical biomarkers produced by demosponges. This view of animal origins contrasts with data from the fossil record, and the taphonomic question of why animals were not preserved (if present) remains unresolved. Neoproterozoic environments demanding small, thin, body plans, and lower abundance/rarity in populations may have played a role. Considering environmental conditions, geochemical data suggest that animals evolved in a relatively low-oxygen ocean. Here, we present new analyses of sedimentary total organic carbon contents in shales suggesting that the Neoproterozoic ocean may also have had lower primary productivity-or at least lower quantities of organic carbon reaching the seafloor-compared with the Phanerozoic. Indeed, recent modeling efforts suggest that low primary productivity is an expected corollary of a low-O2 world. Combined with an inability to inhabit productive regions in a low-O2 ocean, earliest animal communities would likely have been more food limited than generally appreciated, impacting both ecosystem structure and organismal behavior. In light of this, we propose the "fire triangle" metaphor for environmental influences on early animal evolution. Moving toward consideration of all environmental aspects of the Cambrian radiation (fuel, heat, and oxidant) will ultimately lead to a more holistic view of the event.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

Bayesian and parsimony approaches reconstruct informative trees from simulated morphological datasets

Phylogenetic analysis aims to establish the true relationships between taxa. Different analytical methods, however, can reach different conclusions. In order to establish which approach best reconstructs true relationships, previous studies have simulated datasets from known tree topologies, and identified the method that reconstructs the generative tree most accurately. On this basis, researchers have argued that morphological datasets should be analysed by Bayesian approaches, which employ an explicit probabilistic model of evolution, rather than parsimony methods-with implied weights parsimony sometimes identified as particularly inaccurate. Accuracy alone, however, is an inadequate measure of a tree's utility: a fully unresolved tree is perfectly accurate, yet contains no phylogenetic information. The highly resolved trees recovered by implied weights parsimony in fact contain as much useful information as the more accurate, but less resolved, trees recovered by Bayesian methods. By collapsing poorly supported groups, this superior resolution can be traded for accuracy, resulting in trees as accurate as those obtained by a Bayesian approach. By contrast, equally weighted parsimony analysis produces trees that are less resolved and less accurate, leading to less reliable evolutionary conclusions.

Trilobite evolutionary rates constrain the duration of the Cambrian explosion

Trilobites are often considered exemplary for understanding the Cambrian explosion of animal life, due to their unsurpassed diversity and abundance. These biomineralized arthropods appear abruptly in the fossil record with an established diversity, phylogenetic disparity, and provincialism at the beginning of Cambrian Series 2 (∼521 Ma), suggesting a protracted but cryptic earlier history that possibly extends into the Precambrian. However, recent analyses indicate elevated rates of phenotypic and genomic evolution for arthropods during the early Cambrian, thereby shortening the phylogenetic fuse. Furthermore, comparatively little research has been devoted to understanding the duration of the Cambrian explosion, after which normal Phanerozoic evolutionary rates were established. We test these hypotheses by applying Bayesian tip-dating methods to a comprehensive dataset of Cambrian trilobites. We show that trilobites have a Cambrian origin, as supported by the trace fossil record and molecular clocks. Surprisingly, they exhibit constant evolutionary rates across the entire Cambrian, for all aspects of the preserved phenotype: discrete, meristic, and continuous morphological traits. Our data therefore provide robust, quantitative evidence that by the time the typical Cambrian fossil record begins (∼521 Ma), the Cambrian explosion had already largely concluded. This suggests that a modern-style marine biosphere had rapidly emerged during the latest Ediacaran and earliest Cambrian (∼20 million years), followed by broad-scale evolutionary stasis throughout the remainder of the Cambrian.

Global Rate Variation in Bony Vertebrates

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

A general theory of genital homologies for the Hexapoda (Pancrustacea) derived from skeletomuscular correspondences, with emphasis on the Endopterygota

No consensus exists for the homology and terminology of the male genitalia of the Hexapoda despite over a century of debate. Based on dissections and the literature, genital skeletomusculature was compared across the Hexapoda and contrasted with the Remipedia, the closest pancrustacean outgroup. The pattern of origin and insertion for extrinsic and intrinsic genitalic musculature was found to be consistent among the Ectognatha, Protura, and the Remipedia, allowing for the inference of homologies given recent phylogenomic studies. The penis of the Hexapoda is inferred to be derived from medially-fused primary gonopods (gonopore-bearing limbs), while the genitalia of the Ectognatha are inferred to include both the tenth-segmental penis and the ninth-segmental secondary gonopods, similar to the genitalia of female insects which comprise gonopods of the eighth and ninth segments. A new nomenclatural system for hexapodan genitalic musculature is presented and applied, and a general list of anatomical concepts is provided. Novel and refined homologies are proposed for all hexapodan orders, and a series of groundplans are postulated. Emphasis is placed on the Endopterygota, for which fine-grained transition series are hypothesized given observed skeletomuscular correspondences.

History is written by the victors: The effect of the push of the past on the fossil record

Survivorship biases can generate remarkable apparent rate heterogeneities through time in otherwise homogeneous birth‐death models of phylogenies. They are a potential explanation for many striking patterns seen in the fossil record and molecular phylogenies. One such bias is the “push of the past”: clades that survived a substantial length of time are likely to have experienced a high rate of early diversification. This creates the illusion of a secular rate slow‐down through time that is, rather, a reversion to the mean. An extra effect increasing early rates of lineage generation is also seen in large clades. These biases are important but relatively neglected influences on many aspects of diversification patterns in the fossil record and elsewhere, such as diversification spikes after mass extinctions and at the origins of clades; they also influence rates of fossilization, changes in rates of phenotypic evolution and even molecular clocks. These inevitable features of surviving and/or large clades should thus not be generalized to the diversification process as a whole without additional study of small and extinct clades, and raise questions about many of the traditional explanations of the patterns seen in the fossil record.

Tetraconatan phylogeny with special focus on Malacostraca and Branchiopoda: highlighting the strength of taxon-specific matrices in phylogenomics

Understanding the evolution of Tetraconata or Pancrustacea-the clade that includes crustaceans and insects-requires a well-resolved hypothesis regarding the relationships within and among its constituent taxa. Here, we assembled a taxon-rich phylogenomic dataset focusing on crustacean lineages based solely on genomes and new-generation Illumina-generated transcriptomes, including 89 representatives of Tetraconata. This constitutes, to our knowledge, the first phylogenomic study specifically addressing internal relationships of Malacostraca (with 26 species included) and Branchiopoda (36 species). Seven matrices comprising 81-684 orthogroups and 17 690-242 530 amino acid positions were assembled and analysed under five different analytical approaches. To maximize gene occupancy and to improve resolution, taxon-specific matrices were designed for Malacostraca and Branchiopoda. Key tetraconatan taxa (i.e. Oligostraca, Multicrustacea, Branchiopoda, Malacostraca, Thecostraca, Copepoda and Hexapoda) were monophyletic and well supported. Within Branchiopoda, Phyllopoda, Diplostraca, Cladoceromorpha and Cladocera were monophyletic. Within Malacostraca, the clades Eumalacostraca, Decapoda and Reptantia were well supported. Recovery of Caridoida or Peracarida was highly dependent on the analysis for the complete matrix, but it was consistently monophyletic in the malacostracan-specific matrices. From such examples, we demonstrate that taxon-specific matrices and particular evolutionary models and analytical methods, namely CAT-GTR and Dayhoff recoding, outperform other approaches in resolving certain recalcitrant nodes in phylogenomic analyses.

The soft explosive model of placental mammal evolution

BACKGROUND

Recent molecular dating estimates for placental mammals echo fossil inferences for an explosive interordinal diversification, but typically place this event some 10-20 million years earlier than the Paleocene fossils, among apparently more "primitive" mammal faunas.

RESULTS

However, current models of molecular evolution do not adequately account for parallel rate changes, and result in dramatic divergence underestimates for large, long-lived mammals such as whales and hominids. Calibrating among these taxa shifts the rate model errors deeper in the tree, inflating interordinal divergence estimates. We employ simulations based on empirical rate variation, which show that this "error-shift inflation" can explain previous molecular dating overestimates relative to fossil inferences. Molecular dating accuracy is substantially improved in the simulations by focusing on calibrations for taxa that retain plesiomorphic life-history characteristics. Applying this strategy to the empirical data favours the soft explosive model of placental evolution, in line with traditional palaeontological interpretations - a few Cretaceous placental lineages give rise to a rapid interordinal diversification following the 66 Ma Cretaceous-Paleogene boundary mass extinction.

CONCLUSIONS

Our soft explosive model for the diversification of placental mammals brings into agreement previously incongruous molecular, fossil, and ancestral life history estimates, and closely aligns with a growing consensus for a similar model for bird evolution. We show that recent criticism of the soft explosive model relies on ignoring both experimental controls and statistical confidence, as well as misrepresentation, and inconsistent interpretations of morphological phylogeny. More generally, we suggest that the evolutionary properties of adaptive radiations may leave current molecular dating methods susceptible to overestimating the timing of major diversification events.

Early fossil record of Euarthropoda and the Cambrian Explosion

Euarthropoda is one of the best-preserved fossil animal groups and has been the most diverse animal phylum for over 500 million years. Fossil Konservat-Lagerstätten, such as Burgess Shale-type deposits (BSTs), show the evolution of the euarthropod stem lineage during the Cambrian from 518 million years ago (Ma). The stem lineage includes nonbiomineralized groups, such as Radiodonta (e.g., Anomalocaris) that provide insight into the step-by-step construction of euarthropod morphology, including the exoskeleton, biramous limbs, segmentation, and cephalic structures. Trilobites are crown group euarthropods that appear in the fossil record at 521 Ma, before the stem lineage fossils, implying a ghost lineage that needs to be constrained. These constraints come from the trace fossil record, which show the first evidence for total group Euarthropoda (e.g., Cruziana, Rusophycus) at around 537 Ma. A deep Precambrian root to the euarthropod evolutionary lineage is disproven by a comparison of Ediacaran and Cambrian lagerstätten. BSTs from the latest Ediacaran Period (e.g., Miaohe biota, 550 Ma) are abundantly fossiliferous with algae but completely lack animals, which are also missing from other Ediacaran windows, such as phosphate deposits (e.g., Doushantuo, 560 Ma). This constrains the appearance of the euarthropod stem lineage to no older than 550 Ma. While each of the major types of fossil evidence (BSTs, trace fossils, and biomineralized preservation) have their limitations and are incomplete in different ways, when taken together they allow a coherent picture to emerge of the origin and subsequent radiation of total group Euarthropoda during the Cambrian.