Best practices for justifying fossil calibrations

Species limits and hybridization in Andean leaf-eared mice (Phyllotis)

Leaf-eared mice (genus Phyllotis) are among the most widespread and abundant small mammals in the Andean Altiplano, but species boundaries and distributional limits are often poorly delineated due to sparse survey data from remote mountains and high-elevation deserts. Here we report a combined analysis of mitochondrial DNA variation and whole-genome sequence (WGS) variation in Phyllotis mice to delimit species boundaries, to assess the timescale of diversification of the group, and to examine evidence for interspecific hybridization. Estimates of divergence dates suggest that most diversification of Phyllotis occurred during the past 3 million years. Consistent with the Pleistocene Aridification hypothesis, our results suggest that diversification of Phyllotis largely coincided with climatically induced environmental changes in the mid- to late Pleistocene. Contrary to the Montane Uplift hypothesis, most diversification in the group occurred well after the major phase of uplift of the Central Andean Plateau. Species delimitation analyses revealed surprising patterns of cryptic diversity within several nominal forms, suggesting the presence of much undescribed alpha diversity in the genus. Results of genomic analyses revealed evidence of ongoing hybridization between the sister species Phyllotis limatus and P. vaccarum and suggest that the contemporary zone of range overlap between the two species represents an active hybrid zone.

A global blueberry phylogeny: Evolution, diversification, and biogeography of Vaccinieae (Ericaceae)

Vaccinieae is a morphologically diverse and species-rich (∼1430 species) tribe in Ericaceae. Although the majority of diversity is tropical, Vaccinieae are best known for temperate crops (i.e., blueberries, cranberries, and lingonberries) in Vaccinium. Vaccinium itself (∼500 species) has been previously suggested as highly polyphyletic and taxonomic boundaries among many of the other genera in the tribe remain uncertain. We assessed the evolutionary history of Vaccinieae with phylogenomic analyses based on a target-enrichment dataset containing 256 low-copy nuclear loci and 210 species representing 30 of the 35 genera in the tribe and 25 of the 29 sections of Vaccinium. We conducted time-calibrated biogeographic analyses and diversification analyses to explore the area of origin and global dispersal history of the tribe. The analysis recovered a temperate North American origin for Vaccinieae approximately 30 million years ago. Tropical diversity of Vaccinieae was inferred to result from multiple, independent movements into the tropics from north-temperate ancestors. Diversification rate increases corresponded to radiation into the Andes and SE Asia. The pseudo-10-locular ovary evolved once in the tribe from the five-locular state, coinciding with the diversification of a major clade that includes most Asian Vaccinium and the group from which commercial blueberries are derived (V. sect. Cyanococcus). A reconstruction from available chromosome counts suggests that a major polyploid event predated the evolution of nearly half the diversity of Vaccinieae. The extent of polyphyly in Vaccinium documented here supports the need for a generic reclassification of the tribe.

Phylogenomics resolves a 100-year-old debate regarding the evolutionary history of caddisflies (Insecta: Trichoptera)

Trichoptera (caddisfly) phylogeny provides an interesting example of aquatic insect evolution, with rich ecological diversification, especially for underwater architecture. Trichoptera provide numerous critical ecosystem services and are also one of the most important groups of aquatic insects for assessing water quality. The phylogenetic relationships of Trichoptera have been debated for nearly a century. In particular, the phylogenetic position of the "cocoon-makers" within Trichoptera has long been contested. Here, we designed a universal single-copy orthologue and sets of ultraconserved element markers specific for Trichoptera for the first time. Simultaneously, we reconstructed the phylogenetic relationship of Trichoptera based on genome data from 111 species, representing 29 families and 71 genera. Our phylogenetic inference clarifies the probable phylogenetic relationships of "cocoon-makers" within Integripalpia. Hydroptilidae is considered as the basal lineage within Integripalpia, and the families Glossosomatidae, Hydrobiosidae, and Rhyacophilidae formed a monophyletic clade, sister to the integripalpian subterorder Phryganides. The resulting divergence time and ancestral state reconstruction suggest that the most recent common ancestor of Trichoptera appeared in the early Permian and that diversification was strongly correlated with habitat change.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Molecular time estimates for the Lagomorpha diversification

Despite their importance as members of the Glires group, lagomorph diversification processes have seldom been studied using molecular data. Notably, only a few phylogenetic studies have included most of the examined lagomorph lineages. Previous studies that included a larger sample of taxa and markers used nonconservative tests to support the branches of their proposed phylogeny. The objective of this study was to test the monophyly of families and genera of lagomorphs and to evaluate the group diversification process. To that end, this work expanded the sampling of markers and taxa in addition to implementing the bootstrap, a more rigorous statistical test to measure branch support; hence, a more robust phylogeny was recovered. Our supermatrix included five mitochondrial genes and 14 nuclear genes for eighty-eight taxa, including three rodent outgroups. Our maximum likelihood tree showed that all tested genera and both families, Leporidae and Ochotonidae, were recovered as monophyletic. In the Ochotona genus, the subgenera Conothoa and Pika, but not Ochotona, were recovered as monophyletic. Six calibration points based on fossils were used to construct a time tree. A calibration test was performed (via jackknife) by removing one calibration at a time and estimating divergence times for each set. The diversification of the main groups of lagomorphs indicated that the origin of the order's crown group was dated from the beginning of the Palaeogene. Our diversification time estimates for Lagomorpha were compared with those for the largest mammalian order, i.e., rodent lineages in Muroidea. According to our time-resolved phylogenetic tree, the leporids underwent major radiation by evolving a completely new morphospace-larger bodies and an efficient locomotor system-that enabled them to cover wide foraging areas and outrun predators more easily than rodents and pikas.

Co-evolution of early Earth environments and microbial life

Two records of Earth history capture the evolution of life and its co-evolving ecosystems with interpretable fidelity: the geobiological and geochemical traces preserved in rocks and the evolutionary histories captured within genomes. The earliest vestiges of life are recognized mostly in isotopic fingerprints of specific microbial metabolisms, whereas fossils and organic biomarkers become important later. Molecular biology provides lineages that can be overlayed on geologic and geochemical records of evolving life. All these data lie within a framework of biospheric evolution that is primarily characterized by the transition from an oxygen-poor to an oxygen-rich world. In this Review, we explore the history of microbial life on Earth and the degree to which it shaped, and was shaped by, fundamental transitions in the chemical properties of the oceans, continents and atmosphere. We examine the diversity and evolution of early metabolic processes, their couplings with biogeochemical cycles and their links to the oxygenation of the early biosphere. We discuss the distinction between the beginnings of metabolisms and their subsequent proliferation and their capacity to shape surface environments on a planetary scale. The evolution of microbial life and its ecological impacts directly mirror the Earth's chemical and physical evolution through cause-and-effect relationships.

Paleopteran molecular clock: Time drift and recent acceleration

Applying BEAST v1.10.4, we constructed a Bayesian Inference tree comprising 322 taxa, primarily representing Paleoptera (Odonata and Ephemeroptera; Pterygota), Zygentoma and Archaeognatha (Apterygota; paraphyly), and Neoptera (Plecoptera; Pterygota), based on a 2685 bp sequence dataset. Our analyses revealed that robust dating required the incorporation of both Quaternary and pre-Quaternary dates. To achieve this, our dating incorporated a 1.55 Ma (Quaternary) geological event (the formation of the Ryukyu Islands) and a set of chronologically well-founded fossil dates, spanning from up to 400 Ma (Devonian) for the stem Archaeognatha, 320 Ma (Carboniferous) for the crown of Paleoptera, 300 Ma (Carboniferous) for the crown Ephemeroptera, and 280 Ma (Permian) for the crown Odonata, down to 1.76 Ma (Quaternary) for Calopteryx japonica, encompassing a total of 22 calibration points (events: 6, fossils: 16; Quaternary: 7, pre-Quaternary: 15). The resulting dated tree aligns with previous research, albeit with some dates being overestimated. This overestimation was mainly due to the lack of Quaternary calibration and the exclusive dependence on pre-Quaternary calibration, though the application of maximum age constraints also played a role. Our minimum age dating demonstrates that the molecular clock did not uniformly progress, rendering rate dating an inapplicable approach. We observed that the base substitution rate is time-dependent, with an exponential increase evident from around 20 Ma (Miocene) to the present time, exceeding an order of magnitude. The extensive radiation and speciation of Insecta and Paleoptera, potentially resulting from the severe climatic changes associated with the Quaternary, including the commencement of glacial and interglacial cycles, may have significantly contributed to this increase in base substitution rates. Additionally, we identified a potential peak in base substitution rates during the Carboniferous period, around 320 million years ago, possibly corresponding to the Late Paleozoic Ice Age.

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

The Cretaceous Diatom Database: A tool for investigating early diatom evolution

The Cretaceous period is the time of the first appearance of the diatoms in the fossil record. These fossils give us direct evidence of the age and early evolution of the diatom lineage. The fossil record, however, is incomplete and therefore often extrapolated through time-calibrated phylogenies. These two approaches offer different perspectives on the early evolution of diatoms, which is still poorly understood. We compiled the first comprehensive Cretaceous Diatom Database, a tool to investigate the taxonomy, diversity, and occurrence of the earliest known diatom lineages. To further aid the integration and use of the oldest diatom fossils in molecular clock analyses, we present a set of well-documented Cretaceous fossils that can be placed onto molecular phylogenetic trees of extant and extinct species, making them ideal candidates for the calibration of molecular clocks. The analysis of the fossil record and the Cretaceous Diatom Database revealed Cretaceous diversity is substantially greater than previously thought, yet considerable taxonomic work is still needed. The Cretaceous Diatom Database and the list of Cretaceous fossils for calibrating molecular clocks represent valuable resources for future evolutionary and taxonomic studies of modern and fossil diatoms.

Bayesian tip-dated phylogeny and biogeography of Cissampelideae (Menispermaceae): Mitigating the effects of homoplastic morphological characters

The integration of morphological and molecular data is essential to understand the affinities of fossil taxa and spatio-temporal evolutionary processes of organisms. However, homoplastic morphological characters can mislead the placement of fossil taxa and impact downstream analyses. Here, we provide an example of how to mitigate effectively the effect of morphological homoplasy on the placement of fossil taxa and biogeographic inferences of Cissampelideae. We assembled three data types, morphological data only, morphological data with a molecular scaffold and combined morphological and molecular data. By removing high-level homoplastic morphological data or reweighting the morphological characters, we conducted 15 parsimony, 12 undated Bayesian and four dated Bayesian analyses. Our results show that the 14 selected Cissampelideae fossil taxa are placed poorly when based only on morphological data, but the addition of molecular scaffold and combination of morphological and molecular data greatly improve the resolution of fossil nodes. We raise the monotypic Stephania subg. Botryodiscia to generic status and discover that three fossils previously assigned to Stephania should be members of Diploclisia. The Bayesian tip-dated tree recovered by removing homoplastic morphological characters with a Rescaled Consistency Index <0.25 has the highest stratigraphic fit and consequently generates more reasonable biogeographic reconstruction for Cissampelideae. Cissampelideae began to diversify in Asia in the latest Cretaceous and subsequently dispersed to South America around the Cretaceous-Palaeogene boundary. Two dispersal events from Asia to Africa occurred in the Early Eocene and the Late Eocene-Late Oligocene, respectively. These findings provide guidelines and practical methods for mitigating the effects of homoplastic morphological characters on fossil placements and Bayesian tip-dating, as well as insights into the past tropical floristic exchanges among different continents.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Sequential trait evolution did not drive deep-time diversification in sharks

Estimating how traits evolved and impacted diversification across the tree of life represents a critical topic in ecology and evolution. Although there has been considerable research in comparative biology, large parts of the tree of life remain underexplored. Sharks are an iconic clade of marine vertebrates, and key components of marine ecosystems since the early Mesozoic. However, few studies have addressed how traits evolved or whether they impacted their extant diversity patterns. Our study aimed to fill this gap by reconstructing the largest time-calibrated species-level phylogeny of sharks and compiling an exhaustive database for ecological (diet, habitat) and biological (reproduction, maximum body length) traits. Using state-of-the-art models of evolution and diversification, we outlined the major character shifts and modes of trait evolution across shark species. We found support for sequential models of trait evolution and estimated a small to medium-sized lecithotrophic and coastal-dwelling most recent common ancestor for extant sharks. However, our exhaustive hidden traits analyses do not support trait-dependent diversification for any examined traits, challenging previous works. This suggests that the role of traits in shaping sharks' diversification dynamics might have been previously overestimated and should motivate future macroevolutionary studies to investigate other drivers of diversification in this clade.

Convergent Adaptation of True Crabs (Decapoda: Brachyura) to a Gradient of Terrestrial Environments

For much of terrestrial biodiversity, the evolutionary pathways of adaptation from marine ancestors are poorly understood and have usually been viewed as a binary trait. True crabs, the decapod crustacean infraorder Brachyura, comprise over 7600 species representing a striking diversity of morphology and ecology, including repeated adaptation to non-marine habitats. Here, we reconstruct the evolutionary history of Brachyura using new and published sequences of 10 genes for 344 tips spanning 88 of 109 brachyuran families. Using 36 newly vetted fossil calibrations, we infer that brachyurans most likely diverged in the Triassic, with family-level splits in the late Cretaceous and early Paleogene. By contrast, the root age is underestimated with automated sampling of 328 fossil occurrences explicitly incorporated into the tree prior, suggesting such models are a poor fit under heterogeneous fossil preservation. We apply recently defined trait-by-environment associations to classify a gradient of transitions from marine to terrestrial lifestyles. We estimate that crabs left the marine environment at least 7 and up to 17 times convergently, and returned to the sea from non-marine environments at least twice. Although the most highly terrestrial- and many freshwater-adapted crabs are concentrated in Thoracotremata, Bayesian threshold models of ancestral state reconstruction fail to identify shifts to higher terrestrial grades due to the degree of underlying change required. Lineages throughout our tree inhabit intertidal and marginal marine environments, corroborating the inference that the early stages of terrestrial adaptation have a lower threshold to evolve. Our framework and extensive new fossil and natural history datasets will enable future comparisons of non-marine adaptation at the morphological and molecular level. Crabs provide an important window into the early processes of adaptation to novel environments, and different degrees of evolutionary constraint that might help predict these pathways. [Brachyura; convergent evolution; crustaceans; divergence times; fossil calibration; molecular phylogeny; terrestrialization; threshold model.].

The skull of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae)

We provide a detailed and first description of the skull, hyoid apparatus, and trachea of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae). Cyclura is a radiation of iguanas restricted to islands of the Caribbean Sea. Species of Cyclura have high rates of endemism, and all species are severely threatened with extinction. Our anatomical description of this threatened iguana is based on high-resolution computed tomography scans of one adult, one putative adult or near adult, and one juvenile specimen, and includes three-dimensional segmented renderings and visualizations. We discuss some observations of intraspecific and ontogenetic variation, and provide a brief comparison with specimens of another species of Cyclura and published descriptions of other iguanas. Our study provides a cranial osteological framework for Cyclura and augments the body of knowledge on iguana anatomy generally. Finally, we posit that our description and future studies may facilitate identification of fossil Cyclura, which could help understand the paleobiogeography of the genus.

Compound osteoderms preserved in amber reveal the oldest known skink

Scincidae is one of the most species-rich and cosmopolitan clades of squamate reptiles. Abundant disarticulated fossil material has also been attributed to this group, however, no complete pre-Cenozoic crown-scincid specimens have been found. A specimen in Burmite (99 MYA) is the first fossil that can be unambiguously referred to this clade. Our analyses place it as nested within extant skinks, supported by the presence of compound osteoderms formed by articulated small ostedermites. The specimen has a combination of dorsal and ventral compound osteoderms and overlapping cycloid scales that is limited to skinks. We propose that this type of osteoderm evolved as a response to an increased overlap of scales, and to reduced stiffness of the dermal armour. Compound osteoderms could be a key innovation that facilitated diversification in this megadiverse family.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Geomorphic impacts within Red River Fault and island shifting as witnessed by the phylogeography of the largest water strider

Palaeogeological events and climate oscillations profoundly impact the demographics and distributions of small-range species, increasing the extinction risk. The largest water strider worldwide, Gigantometra gigas (Hemiptera: Gerridae), exhibits restricted distributions in Vietnam and southern China. Herein, we generated three genomic datasets (mitogenomes, 146 nuclear protein-coding genes and single nucleotide polymorphisms) with ecological niche modelling (ENM) to explicitly test whether the present-day distribution of G. gigas actually resulted from geographical and climatic effects. We found that the origin of this largest water strider reached the divergence time of the genus within Gerridae, providing a greater opportunity to explore its response to geographic movements. The right-lateral motion of the Red River Fault facilitated the divergence of two phylogeographic lineages, resulting in the "north-south component" genetic pattern in G. gigas. The Hainan and southeast Vietnam populations of the southern linage were completely separated by the Beibu Gulf but exhibited similar genetic compositions, confirming that Hainan had a continental origin and that Hainan Island joined with the Indo-China Peninsula to promote gene exchange among populations. Additionally, we noticed the low genetic diversity but long demographic history of the northern lineage, which displayed population dynamics opposite to those of other organisms. Integrating the demographic changes and ENM findings revealed that suitable habitat contraction and rapid demographic decline during the Last Glacial Maximum (LGM) triggered the low genetic diversity of the northern lineage. Overall, the demographic history of the largest water strider was mainly shaped by geographical features, and first provided evidence from the phylogeographic perspective of aquatic insects to support the hypothesis of Hainan Island shifting.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

A morphological basis for path-dependent evolution of visual systems

Path dependence influences macroevolutionary predictability by constraining potential outcomes after critical evolutionary junctions. Although it has been demonstrated in laboratory experiments, path dependence is difficult to demonstrate in natural systems because of a lack of independent replicates. Here, we show that two types of distributed visual systems recently evolved twice within chitons, demonstrating rapid and path-dependent evolution of a complex trait. The type of visual system that a chiton lineage can evolve is constrained by the number of openings for sensory nerves in its shell plates. Lineages with more openings evolve visual systems with thousands of eyespots, whereas those with fewer openings evolve visual systems with hundreds of shell eyes. These macroevolutionary outcomes shaped by path dependence are both deterministic and stochastic because possibilities are restricted yet not entirely predictable.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Cryogenian Origins of Multicellularity in Archaeplastida

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

Phylogenomics and species delimitation of an abundant and little-studied Amazonian forest spiny rat

Species delimitation studies based on integrating different datasets such as genomic, morphometric, and cytogenetics data are rare in studies focused on Neotropical rodents. As a consequence, the evolutionary history of most of these genera remains poorly understood. Proechimys is a highly diverse and widely distributed genus of Neotropical spiny rats with unique traits like multiple sympatry, micro-habitat segregation, and fuzzy species limits. Here, we applied RAD-Seq to infer the phylogenetic relationships, estimate the species boundaries, and estimate the divergence times for Proechimys, one of the most common and least studied small mammals in the Amazon. We tested whether inferred lineages in the phylogenetic trees could be considered distinct species based on the genomic dataset and morphometric data. Analyses revealed the genus is not monophyletic, with Proechimys hoplomyoides sister to a group of Hoplomys gymnurus + all other Proechimys species, contesting the generic status of Hoplomys. There are five main clades in Proechimys stricto sensu (excluding H. gymnurus and P. hoplomyoides). Species delimitation analyses supported 25 species within the genus Proechimys. The five main clades in Proechimys stricto sensu also showed similar ages for their origins, and two rapid diversification events were identified in the Early Pliocene and in the Early Pleistocene. Most cases of sympatry in Proechimys occur among species from the different main clades, and although Proechimys is an inhabitant of the Amazon, three species occupied the Cerrado biome during the Pleistocene. We could associate available nominal taxon, cytogenetics information, and DNA sequences in Genbank to most of the 25 species we hypothesized from our delimitation analyses. Based on our analyses, we estimate that eight forms represent putative new species that need a taxonomic revision.

A new perspective on the taxonomy and systematics of Arvicolinae (Gray, 1821) and a new time-calibrated phylogeny for the clade

Background

Arvicoline rodents are one of the most speciose and rapidly evolving mammalian lineages. Fossil arvicolines are also among the most common vertebrate fossils found in sites of Pliocene and Pleistocene age in Eurasia and North America. However, there is no taxonomically robust, well-supported, time-calibrated phylogeny for the group.

Methods

Here we present well-supported hypotheses of arvicoline rodent systematics using maximum likelihood and Bayesian inference of DNA sequences of two mitochondrial genes and three nuclear genes representing 146 (82% coverage) species and 100% of currently recognized arvicoline genera. We elucidate well-supported major clades, reviewed the relationships and taxonomy of many species and genera, and critically compared our resulting molecular phylogenetic hypotheses to previously published hypotheses. We also used five fossil calibrations to generate a time-calibrated phylogeny of Arvicolinae that permitted some reconciliation between paleontological and neontological data.

Results

Our results are largely congruent with previous molecular phylogenies, but we increased the support in many regions of the arvicoline tree that were previously poorly-sampled. Our sampling resulted in a better understanding of relationships within Clethrionomyini, the early-diverging position and close relationship of true lemmings (Lemmus and Myopus) and bog lemmings (Synaptomys), and provided support for recent taxonomic changes within Microtini. Our results indicate an origin of ∼6.4 Ma for crown arvicoline rodents. These results have major implications (e.g., diversification rates, paleobiogeography) for our confidence in the fossil record of arvicolines and their utility as biochronological tools in Eurasia and North America during the Quaternary.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

A Palaeogene stem crotaphytid (Aciprion formosum) and the phylogenetic affinities of early fossil pleurodontan iguanians

Pleurodonta is an ancient, diverse clade of iguanian lizard distributed primarily in the Western Hemisphere. Although the clade is a frequent subject of systematic research, phylogenetic resolution among the major pleurodontan clades is elusive. That uncertainty has complicated the interpretations of many fossil pleurodontans. I describe a fossil skull of a pleurodontan lizard from the Palaeogene of Wyoming that was previously allocated to the puzzling taxon Aciprion formosum, and provide an updated morphological matrix for iguanian lizards. Phylogenetic analyses using Bayesian inference demonstrate that the fossil skull is the oldest and first definitive stem member of Crotaphytidae (collared and leopard lizards), establishing the presence of that clade in North America during the Palaeogene. I also discuss new or revised hypotheses for the relationships of several early pleurodontans. In particular, I examine potential evidence for crown-Pleurodonta in the Cretaceous of Mongolia (Polrussia), stem Pleurodonta in the Cretaceous of North America (Magnuviator) and a stem anole in the Eocene of North America (Afairiguana). I suggest that the placement of the fossil crotaphytid is stable to the uncertain phylogeny of Pleurodonta, but recognize the dynamic nature of fossil diagnosis and the potential for updated systematic hypotheses for the other fossils analysed here.

Nucleotide Substitution Model Selection Is Not Necessary for Bayesian Inference of Phylogeny With Well-Behaved Priors

Model selection aims to choose the most adequate model for the statistical analysis at hand. The model must be complex enough to capture the complexity of the data but should be simple enough not to overfit. In phylogenetics, the most common model selection scenario concerns selecting an adequate substitution and partition model for sequence evolution to infer a phylogenetic tree. Previously, several studies showed that substitution model under-parameterization can bias phylogenetic studies. Here, we explored the impact of substitution model over-parameterization in a Bayesian statistical framework. We performed simulations under the simplest substitution model, the Jukes-Cantor model, and compare posterior estimates of phylogenetic tree topologies and tree length under the true model to the most complex model, the $\text{GTR}+\Gamma+\text{I}$ substitution model, including over-splitting the data into additional subsets (i.e., applying partitioned models). We explored 4 choices of prior distributions: the default substitution model priors of MrBayes, BEAST2, and RevBayes and a newly devised prior choice (Tame). Our results show that Bayesian inference of phylogeny is robust to substitution model over-parameterization and over-partitioning but only under our new prior settings. All 3 current default priors introduced biases for the estimated tree length. We conclude that substitution and partition model selection are superfluous steps in Bayesian phylogenetic inference pipelines if well-behaved prior distributions are applied and more effort should focus on more complex and biologically realistic substitution models.

Unraveling phylogenetic relationships and species boundaries in the arid adapted Gerbillus rodents (Muridae: Gerbillinae) by RAD-seq data

Gerbillus is one of the most speciose genera among rodents, with ca. 51 recognized species. Previous attempts to reconstruct the evolutionary history of Gerbillus mainly relied on the mitochondrial cyt-b marker as a source of phylogenetic information. In this study, we utilize RAD-seq genomic data from 37 specimens representing 11 species to reconstruct the phylogenetic tree for Gerbillus, applying concatenation and coalescence methods. We identified four highly supported clades corresponding to the traditionally recognized subgenera: Dipodillus, Gerbillus, Hendecapleura and Monodia. Only two uncertain branches were detected in the resulting trees, with one leading to diversification of the main lineages in the genus, recognized by quartet sampling analysis as uncertain due to possible introgression. We also examined species boundaries for four pairs of sister taxa, including potentially new species from Morocco, using SNAPP. The results strongly supported a speciation model in which all taxa are treated as separate species. The dating analyses confirmed the Plio-Pleistocene diversification of the genus, with the uncertain branch coinciding with the beginning of aridification of the Sahara at the the Plio-Pleistocene boundary. This study aligns well with the earlier analyses based on the cyt-b marker, reaffirming its suitability as an adequate marker for estimating genetic diversity in Gerbillus.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

ATP synthase evolution on a cross-braced dated tree of life

The timing of early cellular evolution, from the divergence of Archaea and Bacteria to the origin of eukaryotes, is poorly constrained. The ATP synthase complex is thought to have originated prior to the Last Universal Common Ancestor (LUCA) and analyses of ATP synthase genes, together with ribosomes, have played a key role in inferring and rooting the tree of life. We reconstruct the evolutionary history of ATP synthases using an expanded taxon sampling set and develop a phylogenetic cross-bracing approach, constraining equivalent speciation nodes to be contemporaneous, based on the phylogenetic imprint of endosymbioses and ancient gene duplications. This approach results in a highly resolved, dated species tree and establishes an absolute timeline for ATP synthase evolution. Our analyses show that the divergence of ATP synthase into F- and A/V-type lineages was a very early event in cellular evolution dating back to more than 4 Ga, potentially predating the diversification of Archaea and Bacteria. Our cross-braced, dated tree of life also provides insight into more recent evolutionary transitions including eukaryogenesis, showing that the eukaryotic nuclear and mitochondrial lineages diverged from their closest archaeal (2.67-2.19 Ga) and bacterial (2.58-2.12 Ga) relatives at approximately the same time, with a slightly longer nuclear stem-lineage.

Critical re-examination of known purported fossil Bombycoidea (Lepidoptera)

We critically re-examine 17 records of fossils currently assigned to the lepidopteran superfamily Bombycoidea, which includes the silk moths, emperor moths and hawk moths. These records include subfossils, compression and impression fossils, permineralizations and ichnofossils. We assess whether observable morphological features warrant their confident assignment to the superfamily. None of the examined fossils displays characters that allow unequivocal identification as Sphingidae, but three fossils and a subfossil (Mioclanis shanwangiana Zhang, Sun and Zhang, 1994, two fossil larvae, and a proboscis in asphaltum) have combinations of diagnostic features that support placement in the family. The identification of a fossil pupa as Bunaeini (Saturniidae) is well supported. The other fossils that we evaluate lack definitive bombycoid and, in several cases, even lepidopteran characters. Some of these dubious fossils have been used as calibration points in earlier studies casting doubt on the resulting age estimates. All fossil specimens reliably assigned to Bombycoidea are relatively young, the earliest fossil evidence of the superfamily dating to the middle Miocene.

MicroCT data provide evidence correcting the previous misidentification of an Eocene amber beetle (Coleoptera, Cicindelidae) as an extant species

The fossil record suggests some insect species have a marked longevity. The oldest fossils purported to represent extant insect species are from the Oligocene and Eocene. One of the most cited fossils is the extant tiger beetle Tetracha carolina (Coleoptera: Cicindelidae) that was identified over a century ago by Walther Horn in Eocene Baltic amber. We examined this and compared it to the previously described cincindelid Baltic amber fossil Palaeoiresina cassolai using X-ray microscopy and 3D imaging techniques. We conclude that Horn's fossil tiger beetle specimen is conspecific with the Eocene P. cassolai and is a member of an extinct stem group lineage of Cicindelidae. Based on a review of all the tiger beetle fossils described from Cretaceous and Paleogene deposits, we found that the assignment of these fossil species to extant lineages is not supported. There are currently no synapomorphies known from fossils that can provide evidence for Cretaceous Manticorni or Megacephalini nor is there evidence for Eocene Iresina. We provide evidence that rejects the idea of a recent beetle species persisting since the Eocene period, which is crucial for using the currently known fossil Cicindelidae species to calibrate divergence dating of beetle phylogenies.

Biogeographic diversification of Actaea (Ranunculaceae): Insights into the historical assembly of deciduous broad-leaved forests in the Northern Hemisphere

The deciduous broad-leaved forests (DBLFs) cover large temperate and subtropical high-altitude regions in the Northern Hemisphere. They are home to rich biodiversity, especially to numerous endemic and relict species. However, we know little about how this vegetation in the Northern Hemisphere has developed through time. Here, we used Actaea (Ranunculaceae), an herbaceous genus almost exclusively growing in the understory of the Northern Hemisphere DBLFs, to shed light on the historical assembly of this biome in the Northern Hemisphere. We present a complete species-level phylogenetic analysis of Actaea based on five plastid and nuclear loci. Using the phylogenetic framework, we estimated divergence times, ancestral ranges, and diversification rates. Phylogenetic analyses strongly support Actaea as monophyletic. Sections Podocarpae and Oligocarpae compose a clade, sister to all other Actaea. The sister relationship between sections Chloranthae and Souliea is strongly supported. Section Dichanthera is not monophyletic unless section Cimicifuga is included. Actaea originated in East Asia, likely the Qinghai-Tibet Plateau, in the late Paleocene (c. 57 Ma), and subsequently dispersed into North America in the middle Eocene (c. 43 Ma) via the Thulean bridge. Actaea reached Europe twice, Japan twice, and Taiwan once, and all these five colonization events occurred in the late Miocene-early Pliocene, a period when sea level dropped. Actaea began to diversify at c. 43 Ma. The section-level diversification took place at c. 27-37 Ma and the species-level diversification experienced accelerations twice, which occurred at c. 15 Ma and c. 5 Ma, respectively. Our findings suggest that the Northern Hemisphere DBLFs might have risen in the middle Eocene and further diversified in the late Eocene-Oligocene, middle Miocene and early Pliocene, in association with climatic deterioration during these four periods.

A hybrid-capture approach to reconstruct the phylogeny of Scleractinia (Cnidaria: Hexacorallia)

A well-supported evolutionary tree representing most major lineages of scleractinian corals is in sight with the development and application of phylogenomic approaches. Specifically, hybrid-capture techniques are shedding light on the evolution and systematics of corals. Here, we reconstructed a broad phylogeny of Scleractinia to test previous phylogenetic hypotheses inferred from a few molecular markers, in particular, the relationships among major scleractinian families and genera, and to identify clades that require further research. We analysed 449 nuclear loci from 422 corals, comprising 266 species spanning 26 families, combining data across whole genomes, transcriptomes, hybrid capture and low-coverage sequencing to reconstruct the largest phylogenomic tree of scleractinians to date. Due to the large number of loci and data completeness (<38% missing data), node supports were high across shallow and deep nodes with incongruences observed in only a few shallow nodes. The "Robust" and "Complex" clades were recovered unequivocally, and our analyses confirmed that Micrabaciidae Vaughan, 1905 is sister to the "Robust" clade, transforming our understanding of the "Basal" clade. Several families remain polyphyletic in our phylogeny, including Deltocyathiidae Kitahara, Cairns, Stolarski & Miller, 2012, Caryophylliidae Dana, 1846, and Coscinaraeidae Benzoni, Arrigoni, Stefani & Stolarski, 2012, and we hereby formally proposed the family name Pachyseridae Benzoni & Hoeksema to accommodate Pachyseris Milne Edwards & Haime, 1849, which is phylogenetically distinct from Agariciidae Gray, 1847. Results also revealed species misidentifications and inconsistencies within morphologically complex clades, such as Acropora Oken, 1815 and Platygyra Ehrenberg, 1834, underscoring the need for reference skeletal material and topotypes, as well as the importance of detailed taxonomic work. The approach and findings here provide much promise for further stabilising the topology of the scleractinian tree of life and advancing our understanding of coral evolution.

Phylogeny and diversification of planthoppers (Hemiptera Fulgoromorpha) based on a comprehensive molecular dataset and large taxon sampling

Our understanding of the evolution of Fulgoromorpha (Insects, Hemiptera) has relied on molecular studies that have only considered either a limited number of taxa where all the families were not represented simultaneously, or a reduced number of genes.The absence of a global analysis comparing all the available data has thus led to significant biases in the analyzes, as evidenced by the incongruence of the results reported for planthopper phylogeny. Here we provide a phylogenetic and dating analysis of the Fulgoromorpha with a large sampling of 531 ingroup taxa, representing about 80% of the currently described suprageneric taxonomic diversity in this group. This study is based on most of the molecular sequences available to date and duly verified, for a set of nuclear and mitochondrial genes from a taxonomic sampling as complete as possible. The most significant results of our study are: (1) the unexpected paraphyly of Delphacidae whose Protodelphacida seem more related to Cixiidae than to other Delphacidae;(2) the group Meenoplidae-Kinnaridae recovered sister to the remaining Fulgoroidea families; (3) the early branching node of Tettigometridae sister of all the other families;(4) the Achilidae-Derbidae clade with Achilidae Plectoderini including Achilixiidae recovered as monophyletic as well as theFulgoridae-Dictyopharidae clade; and (5) the Tropiduchidae placed sister to the other so called 'higher' families (sec. Shcherbakov, 2006).Our divergence times analysis, calibrated with a set of duly verified fossils, suggests that the first diversification of planthoppers occurred in the Early Triassic around 240 Mya and those of the superfamilies Delphacoidea and Fulgoroidea in the Middle-Late Triassic around 210 Mya and 230 Mya, respectively. By the end of the Jurassic, all major planthopper lineages were originated, and all families, around 125 Mya, might havebeen driven in their distribution and evolution (in their first subfamilial divisions) by the geographical constraints of the Gondwanan break-up.Rapid evolutionary radiations occurred particularly in Fulgoridae around 125-130 Mya. Our results stress the importance of the good quality of the sequences used in the molecular analyzes and the primordial importance of a large sampling when analyzing the phylogeny of the group.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

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.

Phylogeny and biogeography of Tiliacoreae (Menispermaceae), a tribe restricted to tropical rainforests

BACKGROUND AND AIMS

Modern tropical rainforests house the highest biodiversity of Earth's terrestrial biomes and are distributed in three low-latitude areas. However, the biogeographical patterns and processes underlying the distribution of biodiversity among these three areas are still poorly known. Here, we used Tiliacoreae, a tribe of pantropical lianas with a high level of regional endemism, to provide new insights into the biogeographical relationships of tropical rainforests among different continents.

METHODS

Based on seven plastid and two nuclear DNA regions, we reconstructed a phylogeny for Tiliacoreae with the most comprehensive sampling ever. Within the phylogenetic framework, we then estimated divergence times and investigated the spatiotemporal evolution of the tribe.

KEY RESULTS

The monophyletic Tiliacoreae contain three major clades, which correspond to Neotropical, Afrotropical and Indo-Malesian/Australasian areas, respectively. Both Albertisia and Anisocycla are not monophyletic. The most recent common ancestor of Tiliacoreae occurred in Indo-Malesia, the Afrotropics and Neotropics in the early Eocene, then rapidly diverged into three major clades between 48 and 46 Ma. Three dispersals from Indo-Malesia to Australasia were inferred, one in the middle Eocene and two in the late Oligocene-late Miocene, and two dispersals from the Afrotropics to Indo-Malesia occurred in the late Eocene-Oligocene.

CONCLUSIONS

The three main clades of Anisocycla correspond to three distinct genera [i.e. Anisocycla sensu stricto and two new genera (Georgesia and Macrophragma)]. Epinetrum is a member of Albertisia. Our findings highlight that sea-level fluctuations and climate changes in the Cenozoic have played important roles in shaping the current distribution and endemism of Tiliacoreae, hence contributing to the knowledge on the historical biogeography of tropical rainforests on a global scale.

Small population of the largest water strider after the late Pleistocene and the implications for its conservation

Climate oscillation and its synergistic impacts on habitat fragmentation have been identified as threatening the survival of some extant species. However, the mechanisms by which semi-aquatic insects impacted by such events remain poorly understood. Herein, we studied the largest water strider in the world, Gigantometra gigas, to explore the effect of these two factors on its evolutionary history. The sequences of mitogenomic and nrDNA cluster were utilized to reconstruct phylogenetic relationship among G. gigas populations and its demographic history. Mitochondrial genes were separately reconstructed topologies of that populations and detected remarkable differences. We found that G. gigas populations conform to the isolation-by-distance model, and decline occurred at about 120 ka, which was probably influenced by the climate change during the late Pleistocene and eventually maintained a small effective population size (Ne) around 85,717. The populations in Guangdong Province of China are worthy of note in that they exhibit low genetic diversity, a small Ne around 18,899 individuals, and occupy an area with little suitable future habitat for G. gigas. This work recommends that conservation efforts are implemented to ensure the long-term survival of small G. gigas populations, and notes that further evaluation of their extinction risk under the impacts of human activities is required.

A complete genus-level phylogeny reveals the Cretaceous biogeographic diversification of the poppy family

Angiosperms, a trigger for the Cretaceous Terrestrial Revolution (KTR), underwent a rapid expansion and occupied all the environments during the Mid-Upper Cretaceous. Yet, Cretaceous biogeographic patterns and processes underlying the distribution of angiosperm diversity in the Northern Hemisphere are still poorly known. Here, we elucidated the biogeographic diversification of the angiosperm family Papaveraceae, an ancient Northern Hemisphere clade characterized by poor dispersal ability and high level of regional endemism. Based on both plastome and multi-locus datasets, we reconstructed a robust time-calibrated phylogeny that includes all currently recognized 45 genera of this family. Within the time-calibrated phylogenetic framework, we conducted 72 biogeographic analyses by testing the sensitivity of uncertainties of area delimitation, maxarea constraints, and the parameters of the model, i.e., j (describing jump-dispersal events) and w (modifying dispersal multiplier matrices), to ancestral range estimations. We also inferred ancestral habitat and ecological niches. Phylogenetic analyses strongly support Papaveraceae as monophyletic. Pteridophylloideae is strongly supported as sister to Hypecoideae-Fumarioideae. Our results indicate that the j parameter and number of predefined areas strongly affect ancestral range estimates, generating questionable ancestral ranges, whereas maxarea constraint and w parameter have no effect and improve model fit. After accounting for these uncertainties, our results indicate that Papaveraceae differentiated in Asian wet forests during the Lower Cretaceous and subsequently occupied the Asian and western North American arid and open areas. Three dispersals from Asia to western North America via the Bering land bridge occurred in the Mid-Upper Cretaceous, largely in agreement with the KTR. Habitat shift and ecological niche divergence resulted in the subsequent disjunctions between Asia and western North America. These findings suggest that the interplay of range expansion and niche divergence-driven vicariance might have shaped Cretaceous biogeographic patterns of angiosperms with Papaveraceae-like ecological requirements and dispersal abilities in the Northern Hemisphere, hence contributing to the knowledge on the geographic expansion of angiosperms during the KTR.

The Australian gulf snapping turtle Elseya lavarackorum (Testudines: Chelidae) revisited—Is the late Pleistocene fossil species extant?

Disagreement exists on the taxonomic identity of the extant populations of the Australian Elseya referred to in 1992 as the gulf Elseya (= Elseya sp. aff. dentata [Nicholson]). The extant form has since 1997 been considered conspecific with the late Pleistocene fossil Elseya lavarackorum (White and Archer, 1994). Recently it has been considered a new species, Elseya oneirosJoseph-Ouni et al., 2020, conspecific with another fossil found in the same site and stratum as Elseya lavarackorum. Here we re-examine the fossil material and reassess the characters used by previous authors in an attempt to decide the issue. We find that the anterior bridge suture with the carapace of the fossil Elseya lavarackorum is associated with extensive and prominent plastral elements, which has led to misinterpretation of characters associated with this structure. We furthermore show that interindividual variation in sulci patterns is so great as to render them of little taxonomic value. On the basis of (a) deviation of the anterior shape of the carapace from ovoid such that, in aged individuals, the most anterior point of the carapace occurs at marginal scutes M2 (a resultant nuchal bay occurs in such individuals); (b) the typical absence of a cervical scute; (c) no evidence of a medial constriction in the anterior bridge strut suture; and (d) absence of evidence of any other informative variation of taxonomic value; we conclude that the decision to consider the late Pleistocene (ca 23 kyr old) fossil and the extant Elseya sp. aff. dentata [Nicholson] as Elseya lavarackorum (White and Archer, 1994) as conspecific should stand.

Phylogenomic analyses of Sapindales support new family relationships, rapid Mid-Cretaceous Hothouse diversification, and heterogeneous histories of gene duplication

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order's spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research.

Key innovations and the diversification of Hymenoptera

The order Hymenoptera (wasps, ants, sawflies, and bees) represents one of the most diverse animal lineages, but whether specific key innovations have contributed to its diversification is still unknown. We assembled the largest time-calibrated phylogeny of Hymenoptera to date and investigated the origin and possible correlation of particular morphological and behavioral innovations with diversification in the order: the wasp waist of Apocrita; the stinger of Aculeata; parasitoidism, a specialized form of carnivory; and secondary phytophagy, a reversal to plant-feeding. Here, we show that parasitoidism has been the dominant strategy since the Late Triassic in Hymenoptera, but was not an immediate driver of diversification. Instead, transitions to secondary phytophagy (from parasitoidism) had a major influence on diversification rate in Hymenoptera. Support for the stinger and the wasp waist as key innovations remains equivocal, but these traits may have laid the anatomical and behavioral foundations for adaptations more directly associated with diversification.

Phylogenomics reconciles molecular data with the rich fossil record on the origin of living turtles

Although a consensus exists that all living turtles fall within either Pleurodira or Cryptodira clades, estimating when these lineages split is still under debate. Most molecular studies date the split in the Triassic Period, whereas a Jurassic age is unanimous among morphological studies. Each hypothesis implies different paleobiogeographical scenarios to explain early turtle evolution. Here we explored the rich turtle fossil record with the Fossilized Birth-Death (FBD) and the traditional node dating (ND) methods using complete mitochondrial genomes (147 taxa) and a set of nuclear orthologs with over 10 million bp (25 taxa) to date the major splits in Testudines. Our results support an Early Jurassic split (191-182 Ma) for the crown Testudines with great consistency across different dating methods and datasets, with a narrow confidence interval. This result is independently supported by the oldest fossils of Testudines that postdate the Middle Jurassic (174 Ma), which were not used for calibration in this study. This age coincides with the Pangaea fragmentation and the formation of saltwater barriers such as the Atlantic Ocean and the Turgai Strait, supporting that diversification in Testudines was triggered by vicariance. Our ages of the splits in Pleurodira coincide with the geologic events of the Late Jurassic and Early Cretaceous. Conversely, the early Cryptodira radiation remained in Laurasia, and its diversification ensued as all its major lineages expanded their distribution into every continent during the Cenozoic. We provide the first detailed hypothesis of the evolution of Cryptodira in the Southern Hemisphere, in which our time estimates are correlated with each contact between landmasses derived from Gondwana and Laurasia. Although most South American Cryptodira arrived through the Great American Biotic Interchange, our results indicate that the Chelonoidis ancestor probably arrived from Africa through the chain islands of the South Atlantic during the Paleogene. Together, the presence of ancient turtle diversity and the vital role that turtles occupy in marine and terrestrial ecosystems underline South America as a chief area for conservation.

Systematics and Phylogenetic Interrelationships of the Enigmatic Late Jurassic Shark Protospinax annectans Woodward, 1918 with Comments on the Shark-Ray Sister Group Relationship

The Late Jurassic elasmobranch Protospinax annectans is often regarded as a key species to our understanding of crown group elasmobranch interrelationships and the evolutionary history of this group. However, since its first description more than 100 years ago, its phylogenetic position within the Elasmobranchii (sharks and rays) has proven controversial, and a closer relationship between Protospinax and each of the posited superorders (Batomorphii, Squalomorphii, and Galeomorphii) has been proposed over the time. Here we revise this controversial taxon based on new holomorphic specimens from the Late Jurassic Konservat-Lagerstätte of the Solnhofen Archipelago in Bavaria (Germany) and review its skeletal morphology, systematics, and phylogenetic interrelationships. A data matrix with 224 morphological characters was compiled and analyzed under a molecular backbone constraint. Our results indicate a close relationship between Protospinax, angel sharks (Squatiniformes), and saw sharks (Pristiophoriformes). However, the revision of our morphological data matrix within a molecular framework highlights the lack of morphological characters defining certain groups, especially sharks of the order Squaliformes, hampering the phylogenetic resolution of Protospinax annectans with certainty. Furthermore, the monophyly of modern sharks retrieved by molecular studies is only weakly supported by morphological data, stressing the need for more characters to align morphological and molecular studies in the future.

Basal Anseriformes from the Early Paleogene of North America and Europe

We describe nearly complete skeletons of basal Anseriformes from the Latest Paleocene to the early Eocene of North America and Europe. Collectively, these birds appear to be representative of anseriforms near the divergence of Anhimae and Anseres, but their exact positions relative to these clades remains uncertain. A new family, Anachronornithidae nov. fam., is erected on the basis of one of these, Anachronornis anhimops nov. gen., nov. gen. et sp., to which the others cannot be confidently assigned. The new fossils augment a growing collection of early Pan-Anseriformes, which in their diversity do not paint an unambiguous picture of phylogeny or character state evolution on the path to or within crown-Anseriformes. Anachronornis nov. gen. is similar in some aspects of both cranial and postcranial anatomy to other well-represented early Paleogene Anseriformes and members of Anseres, such as Presbyornis Wetmore, 1926. However, it exhibits a more landfowl-like bill, like that of Anhimae and unlike the spatulate bill of Anseres. Additional specimens of similar basal Anseriformes of uncertain affinities from the early Eocene of North America and Europe further complicate interpretation of character state polarity due to the mosaicism of primitive and derived characters they exhibit.

Deriving calibrations for Arawakan using archaeological evidence

This paper identifies time calibration points for accurately rooting and dating the phylogeny of Arawakan, the largest Indigenous linguistic family of the Americas. We present and model a methodology for extracting calibration points from the archaeological record, based on principles of geographical overlap between archaeological sites and Arawakan peoples, and on continuity in material culture between archaeological finds and modern Arawakan practices. Based on a consensus model of the expansion of the Arawakan family from Central Amazonia, we focus on archaeological finds in Arawakan expansion zones, where Arawakan material culture abruptly appears in a given region, and where only a single major Arawakan subgroup/clade is present. We find 12 calibration points from archaeological sites in Arawakan expansion zones and also identify more recent calibration points from the historical record based on first mentions of ethnonyms and early sources of lexical data.

No phylogenomic support for a Cenozoic origin of the "living fossil" Isoetes

PREMISE

The isoetalean lineage has a rich fossil record that extends to the Devonian, but the age of the living clade is unclear. Recent results indicate that it is young, from the Cenozoic, whereas earlier work based on less data from a denser taxon sampling yielded Mesozoic median ages.

METHODS

We analyzed node ages in Isoetes using two genomic data sets (plastome and nuclear ribosomal cistron), three clock models implemented in MrBayes (ILN, WN, and TK02 models), and a conservative approach to calibration.

RESULTS

While topological results were consistently resolved in Isoetes estimated crown group ages range from the latest Paleozoic (mid-Permian) to the Mesozoic depending on data type and clock model. The oldest estimates were retrieved using the autocorrelated TK02 clock model. An (early) Cenozoic age was only obtained under one specific condition (plastome data analyzed with the uncorrelated ILN clock model). That same plastome data set also yielded the oldest (mid-Permian) age estimate when analyzed with the autocorrelated TK02 clock model. Adding the highly divergent, recently established sister species Isoetes wormaldii to the data set approximately doubled the average median node depth to the Isoetes crown group.

CONCLUSIONS

There is no consistent support for a Cenozoic origin of the living clade Isoetes. We obtained seemingly well-founded, yet strongly deviating results depending on data type and clock model. The single most important future improvement is probably to add calibration points, which requires an improved understanding of the isoetalean fossil record or alternative bases for calibration.

Constraining Whole-Genome Duplication Events in Geological Time

The timing of whole-genome duplication (WGD) events is crucial to understanding their role in evolution and underpins many hypotheses linking WGD to increased diversity and complexity. As such, means of estimating the timing of the WGD events relative to their macroevolutionary outcomes are of considerable importance. Molecular clock methods facilitate direct estimation of the absolute timing of WGD events, integrating information on the rate of sequence evolution between species while accommodating the uncertainty inherent to the fossil record. We present an explanation of the best practice for constructing fossil calibrations and estimating the age of WGD events via molecular clock methods in the program MCMCtree, with an example dataset based on a well-characterized WGD event within the flowering dogwoods (Cornus). The approach presented herein allows for the estimation of the age of WGD events and subsequent speciation events, allowing the relationship between WGD and the macroevolutionary outcomes to be explored. In our example, we show that in the case of flowering dogwoods, the WGD event long predates the end-Cretaceous mass extinction and that the two events may be independent.

The Jasmine (Jasminum sambac) Genome Provides Insight into the Biosynthesis of Flower Fragrances and Jasmonates

Jasminum sambac (jasmine flower), a world-renowned plant appreciated for its exceptional flower fragrance, is of cultural and economic importance. However, the genetic basis of its fragrance is largely unknown. Here, we present the first de novo genome of J. sambac with 550.12 Mb (scaffold N50 = 40.10 Mb) assembled into 13 pseudochromosomes. Terpene synthase genes associated with flower fragrance are significantly amplified in the form of gene clusters through tandem duplications in the genome. Gene clusters within the salicylic acid/benzoic acid/theobromine (SABATH) and BAHD superfamilies were identified as related to the biosynthesis of phenylpropanoid/benzenoid compounds. Several key genes involved in jasmonate biosynthesis were duplicated, causing increased copy numbers. In addition, multi-omics analyses identified various aromatic compounds and many genes involved in fragrance biosynthesis pathways. Furthermore, the roles of JsTPS3 in β-ocimene biosynthesis, as well as JsAOC1 and JsAOS in jasmonic acid biosynthesis, were functionally validated. The genome assembled in this study for J. sambac offers a basic genetic resource for studying floral scent and jasmonate biosynthesis and provides a foundation for functional genomic research and variety improvements in Jasminum.