Bayesian molecular clock dating of species divergences in the genomics era

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.

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families

Negevirus is a recently proposed taxon of arthropod-infecting virus, which is associated with plant viruses of two families (Virgaviridae and Kitaviridae). Nevertheless, the evolutionary history of negevirus-host and its relationship with plant viruses remain poorly understood. Endogenous nege-like viral elements (ENVEs) are ancient nege-like viral sequences integrated into the arthropod genomes, which can serve as the molecular fossil records of previous viral infection. In this study, 292 ENVEs were identified in 150 published arthropod genomes, revealing the evolutionary history of nege-like viruses and two related plant virus families. We discovered three novel and eight strains of nege-like viruses in 11 aphid species. Further analysis indicated that 10 ENVEs were detected in six aphid genomes, and they were divided into four types (ENVE1-ENVE4). Orthologous integration and phylogenetic analyses revealed that nege-like viruses had a history of infection of over 60 My and coexisted with aphid ancestors throughout the Cenozoic Era. Moreover, two nege-like viral proteins (CP and SP24) were highly homologous to those of plant viruses in the families Virgaviridae and Kitaviridae. CP- and SP24-derived ENVEs were widely integrated into numerous arthropod genomes. These results demonstrate that nege-like viruses have a long-term coexistence with arthropod hosts and plant viruses of the two families, Virgaviridae and Kitaviridae, which may have evolved from the nege-like virus ancestor through horizontal virus transfer events. These findings broaden our perspective on the history of viral infection in arthropods and the origins of plant viruses.

IMPORTANCE

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

Challenges in Assembling the Dated Tree of Life

The assembly of a comprehensive and dated Tree of Life (ToL) remains one of the most formidable challenges in evolutionary biology. The complexity of life's history, involving both vertical and horizontal transmission of genetic information, defies its representation by a simple bifurcating phylogeny. With the advent of genome and metagenome sequencing, vast amounts of data have become available. However, employing this information for phylogeny and divergence time inference has introduced significant theoretical and computational hurdles. This perspective addresses some key methodological challenges in assembling the dated ToL, namely, the identification and classification of homologous genes, accounting for gene tree-species tree mismatch due to population-level processes along with duplication, loss, and horizontal gene transfer, and the accurate dating of evolutionary events. Ultimately, the success of this endeavor requires new approaches that integrate knowledge databases with optimized phylogenetic algorithms capable of managing complex evolutionary models.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

Relaxing the molecular clock using models of how substitution rates change across lineages has become essential for addressing evolutionary problems. The diversity of rate evolution models and their implementations are substantial, and studies have demonstrated their impact on divergence time estimates can be as significant as that of calibration information. In this review, we trace the development of rate evolution models from the proposal of the molecular clock concept to the development of sophisticated Bayesian and non-Bayesian methods that handle rate variation in phylogenies. We discuss the various approaches to modeling rate evolution, provide a comprehensive list of available software, and examine the challenges and advancements of the prevalent Bayesian framework, contrasting them to faster non-Bayesian methods. Lastly, we offer insights into potential advancements in the field in the era of big data.

Proterozoic microfossils continue to provide new insights into the rise of complex eukaryotic life

Eukaryotes have evolved to dominate the biosphere today, accounting for most documented living species and the vast majority of the Earth's biomass. Consequently, understanding how these biologically complex organisms initially diversified in the Proterozoic Eon over 539 million years ago is a foundational question in evolutionary biology. Over the last 70 years, palaeontologists have sought to document the rise of eukaryotes with fossil evidence. However, the delicate and microscopic nature of their sub-cellular features affords early eukaryotes diminished preservation potential. Chemical biomarker signatures of eukaryotes and the genetics of living eukaryotes have emerged as complementary tools for reconstructing eukaryote ancestry. In this review, we argue that exceptionally preserved Proterozoic microfossils are critical to interpreting these complementary tools, providing crucial calibrations to molecular clocks and testing hypotheses of palaeoecology. We highlight recent research on their preservation and biomolecular composition that offers new ways to enhance their utility.

A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history

The genera Cephalotrichum and Microascus contain ecologically, morphologically and lifestyle diverse fungi in Microascaceae (Microascales, Sordariomycetes) with a world-wide distribution. Despite previous studies having elucidated that Cephalotrichum and Microascus are highly polyphyletic, the DNA phylogeny of many traditionally morphology-defined species is still poorly resolved, and a comprehensive taxonomic overview of the two genera is lacking. To resolve this issue, we integrate broad taxon sampling strategies and the most comprehensive multi-gene (ITS, LSU, tef1 and tub2) datasets to date, with fossil calibrations to address the phylogenetic relationships and divergence times among major lineages of Microascaceae. Two previously recognised main clades, Cephalotrichum (24 species) and Microascus (49 species), were re-affirmed based on our phylogenetic analyses, as well as the phylogenetic position of 15 genera within Microascaceae. In this study, we provide an up-to-date overview on the taxonomy and phylogeny of species belonging to Cephalotrichum and Microascus, as well as detailed descriptions and illustrations of 21 species of which eight are newly described. Furthermore, the divergence time estimates indicate that the crown age of Microascaceae was around 210.37 Mya (95 % HPD: 177.18-246.96 Mya) in the Late Triassic, and that Cephalotrichum and Microascus began to diversify approximately 27.07 Mya (95 % HPD: 20.47-34.37 Mya) and 70.46 Mya (95 % HPD: 56.96-86.24 Mya), respectively. Our results also demonstrate that multigene sequence data coupled with broad taxon sampling can help elucidate previously unresolved clade relationships. Citation: Wei TP, Wu YM, Zhang X, et al. 2024. A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history. Persoonia 52: 119-160. https://doi.org/10.3767/persoonia.2024.52.05 .

Phylogeny and biogeography of harmochirine jumping spiders (Araneae: Salticidae)

We use ultraconserved elements (UCE) and Sanger data to study the phylogeny, age, and biogeographical history of harmochirine jumping spiders, a group that includes the species-rich genus Habronattus, whose remarkable courtship has made it the focus of studies of behaviour, sexual selection, and diversification. We recovered 1947 UCE loci from 43 harmochirine taxa and 4 outgroups, yielding a core dataset of 193 UCEs with at least 50 % occupancy. Concatenated likelihood and ASTRAL analyses confirmed the separation of harmochirines into two major clades, here designated the infratribes Harmochirita and Pellenita. Most are African or Eurasian with the notable exception of a clade of pellenites containing Habronattus and Pellenattus of the Americas and Havaika and Hivanua of the Pacific Islands. Biogeographical analysis using the DEC model favours a dispersal of the clade's ancestor from Eurasia to the Americas, from which Havaika's ancestor dispersed to Hawaii and Hivanua's ancestor to the Marquesas Islands. Divergence time analysis on 32 loci with 85 % occupancy, calibrated by fossils and island age, dates the dispersal to the Americas at approximately 4 to 6 million years ago. The explosive radiation of Habronattus perhaps began only about 4 mya. The phylogeny clarifies both the evolution of sexual traits (e.g., the terminal apophyses was enlarged in Pellenes and not subsequently lost) and the taxonomy. Habronattus is confirmed as monophyletic. Pellenattus is raised to the status of genus, and 13 species moved into it as new combinations. Bianor stepposus Logunov, 1991 is transferred to Sibianor, and Pellenes bulawayoensis Wesołowska, 1999 is transferred to Neaetha. A molecular clock rate estimate for spider UCEs is presented and its utility to inform prior distributions is discussed.

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.].

Dating Ammonia-Oxidizing Bacteria with Abundant Eukaryotic Fossils

Evolution of a complete nitrogen (N) cycle relies on the onset of ammonia oxidation, which aerobically converts ammonia to nitrogen oxides. However, accurate estimation of the antiquity of ammonia-oxidizing bacteria (AOB) remains challenging because AOB-specific fossils are absent and bacterial fossils amenable to calibrate molecular clocks are rare. Leveraging the ancient endosymbiosis of mitochondria and plastid, as well as using state-of-the-art Bayesian sequential dating approach, we obtained a timeline of AOB evolution calibrated largely by eukaryotic fossils. We show that the first AOB evolved in marine Gammaproteobacteria (Gamma-AOB) and emerged between 2.1 and 1.9 billion years ago (Ga), thus postdating the Great Oxidation Event (GOE; 2.4 to 2.32 Ga). To reconcile the sedimentary N isotopic signatures of ammonia oxidation occurring near the GOE, we propose that ammonia oxidation likely occurred at the common ancestor of Gamma-AOB and Gammaproteobacterial methanotrophs, or the actinobacterial/verrucomicrobial methanotrophs which are known to have ammonia oxidation activities. It is also likely that nitrite was transported from the terrestrial habitats where ammonia oxidation by archaea took place. Further, we show that the Gamma-AOB predated the anaerobic ammonia-oxidizing (anammox) bacteria, implying that the emergence of anammox was constrained by the availability of dedicated ammonia oxidizers which produce nitrite to fuel anammox. Our work supports a new hypothesis that N redox cycle involving nitrogen oxides evolved rather late in the ocean.

ClockstaRX: Testing Molecular Clock Hypotheses With Genomic Data

Phylogenomic data provide valuable opportunities for studying evolutionary rates and timescales. These analyses require theoretical and statistical tools based on molecular clocks. We present ClockstaRX, a flexible platform for exploring and testing evolutionary rate signals in phylogenomic data. Here, information about evolutionary rates in branches across gene trees is placed in Euclidean space, allowing data transformation, visualization, and hypothesis testing. ClockstaRX implements formal tests for identifying groups of loci and branches that make a large contribution to patterns of rate variation. This information can then be used to test for drivers of genomic evolutionary rates or to inform models for molecular dating. Drawing on the results of a simulation study, we recommend forms of data exploration and filtering that might be useful prior to molecular-clock analyses.

The Sphinx and the egg: Evolutionary enigmas of the (glyco)sphingolipid biosynthetic pathway

In eukaryotes, the de novo synthesis of sphingolipids (SLs) consists of multiple sequential steps which are compartmentalized between the endoplasmic reticulum and the Golgi apparatus. Studies over many decades have identified the enzymes in the pathway, their localization, topology and an array of regulatory mechanisms. However, little is known about the evolutionary forces that underly the generation of this complex pathway or of its anteome, i.e., the metabolic pathways that converge on the SL biosynthetic pathway and are essential for its activity. After briefly describing the pathway, we discuss the mechanisms by which the enzymes of the SL biosynthetic pathway are targeted to their different subcellular locations, how the pathway per se may have evolved, including its compartmentalization, and the relationship of the pathway to eukaryogenesis. We discuss the circular interdependence of the evolution of the SL pathway, and comment on whether current Darwinian evolutionary models are able to provide genuine mechanistic insight into how the pathway came into being.

Insight into Hyalomma anatolicum biology by comparative genomics analyses

Hyalomma anatolicum is an obligatory blood-sucking ectoparasite and contributes to the transmission of Crimean-Congo haemorrhagic fever (CCHF) virus, Theileria spp. and Babesia spp. Progress in exploring the adaptive strategy of this ectoparasite and developing tools to fight it has been hindered by the lack of a complete genome. Herein, we assembled the genome using diverse sources of data from multiple sequencing platforms and annotated the 1.96 Gb genome of Hy. anatolicum. Comparative genome analyses and the predicted protein encoding genes reveal unique facets of this genome, including gene family expansion associated with blood feeding and digestion, multi-gene families involved in detoxification, a great number of neuropeptides and corresponding receptors regulating tick growth, development, and reproduction, and glutathione S-transferase genes playing roles in insecticide resistance and detoxification of multiple xenobiotic factors. This high quality reference genome provides fundamental data for obtaining insights into a variety of aspects of tick biology and developing novel strategies to fight notorious tick vectors of human and animal pathogens.

Phylogenomic insights into the first multicellular streptophyte

Streptophytes are best known as the clade containing the teeming diversity of embryophytes (land plants).1,2,3,4 Next to embryophytes are however a range of freshwater and terrestrial algae that bear important information on the emergence of key traits of land plants. Among these, the Klebsormidiophyceae stand out. Thriving in diverse environments-from mundane (ubiquitous occurrence on tree barks and rocks) to extreme (from the Atacama Desert to the Antarctic)-Klebsormidiophyceae can exhibit filamentous body plans and display remarkable resilience as colonizers of terrestrial habitats.5,6 Currently, the lack of a robust phylogenetic framework for the Klebsormidiophyceae hampers our understanding of the evolutionary history of these key traits. Here, we conducted a phylogenomic analysis utilizing advanced models that can counteract systematic biases. We sequenced 24 new transcriptomes of Klebsormidiophyceae and combined them with 14 previously published genomic and transcriptomic datasets. Using an analysis built on 845 loci and sophisticated mixture models, we establish a phylogenomic framework, dividing the six distinct genera of Klebsormidiophyceae in a novel three-order system, with a deep divergence more than 830 million years ago. Our reconstructions of ancestral states suggest (1) an evolutionary history of multiple transitions between terrestrial-aquatic habitats, with stem Klebsormidiales having conquered land earlier than embryophytes, and (2) that the body plan of the last common ancestor of Klebsormidiophyceae was multicellular, with a high probability that it was filamentous whereas the sarcinoids and unicells in Klebsormidiophyceae are likely derived states. We provide evidence that the first multicellular streptophytes likely lived about a billion years ago.

The Genome of Plasmodium gonderi: Insights into the Evolution of Human Malaria Parasites

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Dating ancient splits in phylogenetic trees, with application to the human-Neanderthal split

BACKGROUND

We tackle the problem of estimating species TMRCAs (Time to Most Recent Common Ancestor), given a genome sequence from each species and a large known phylogenetic tree with a known structure (typically from one of the species). The number of transitions at each site from the first sequence to the other is assumed to be Poisson distributed, and only the parity of the number of transitions is observed. The detailed phylogenetic tree contains information about the transition rates in each site. We use this formulation to develop and analyze multiple estimators of the species' TMRCA. To test our methods, we use mtDNA substitution statistics from the well-established Phylotree as a baseline for data simulation such that the substitution rate per site mimics the real-world observed rates.

RESULTS

We evaluate our methods using simulated data and compare them to the Bayesian optimizing software BEAST2, showing that our proposed estimators are accurate for a wide range of TMRCAs and significantly outperform BEAST2. We then apply the proposed estimators on Neanderthal, Denisovan, and Chimpanzee mtDNA genomes to better estimate their TMRCA with modern humans and find that their TMRCA is substantially later, compared to values cited recently in the literature.

CONCLUSIONS

Our methods utilize the transition statistics from the entire known human mtDNA phylogenetic tree (Phylotree), eliminating the requirement to reconstruct a tree encompassing the specific sequences of interest. Moreover, they demonstrate notable improvement in both running speed and accuracy compared to BEAST2, particularly for earlier TMRCAs like the human-Chimpanzee split. Our results date the human - Neanderthal TMRCA to be [Formula: see text] years ago, considerably later than values cited in other recent studies.

Leaping through Tree Space: Continuous Phylogenetic Inference for Rooted and Unrooted Trees

Phylogenetics is now fundamental in life sciences, providing insights into the earliest branches of life and the origins and spread of epidemics. However, finding suitable phylogenies from the vast space of possible trees remains challenging. To address this problem, for the first time, we perform both tree exploration and inference in a continuous space where the computation of gradients is possible. This continuous relaxation allows for major leaps across tree space in both rooted and unrooted trees, and is less susceptible to convergence to local minima. Our approach outperforms the current best methods for inference on unrooted trees and, in simulation, accurately infers the tree and root in ultrametric cases. The approach is effective in cases of empirical data with negligible amounts of data, which we demonstrate on the phylogeny of jawed vertebrates. Indeed, only a few genes with an ultrametric signal were generally sufficient for resolving the major lineages of vertebrates. Optimization is possible via automatic differentiation and our method presents an effective way forward for exploring the most difficult, data-deficient phylogenetic questions.

A timescale for placental mammal diversification based on Bayesian modeling of the fossil record

The timing of the placental mammal radiation has been the focus of debate over the efficacy of competing methods for establishing evolutionary timescales. Molecular clock analyses estimate that placental mammals originated before the Cretaceous-Paleogene (K-Pg) mass extinction, anywhere from the Late Cretaceous to the Jurassic. However, the absence of definitive fossils of placentals before the K-Pg boundary is compatible with a post-Cretaceous origin. Nevertheless, lineage divergence must occur before it can be manifest phenotypically in descendent lineages. This, combined with the non-uniformity of the rock and fossil records, requires the fossil record to be interpreted rather than read literally. To achieve this, we introduce an extended Bayesian Brownian bridge model that estimates the age of origination and, where applicable, extinction through a probabilistic interpretation of the fossil record. The model estimates the origination of placentals in the Late Cretaceous, with ordinal crown groups originating at or after the K-Pg boundary. The results reduce the plausible interval for placental mammal origination to the younger range of molecular clock estimates. Our findings support both the Long Fuse and Soft Explosive models of placental mammal diversification, indicating that the placentals originated shortly prior to the K-Pg mass extinction. The origination of many modern mammal lineages overlapped with and followed the K-Pg mass extinction.

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

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

Molecular timetrees using relaxed clocks and uncertain phylogenies

A common practice in molecular systematics is to infer phylogeny and then scale it to time by using a relaxed clock method and calibrations. This sequential analysis practice ignores the effect of phylogenetic uncertainty on divergence time estimates and their confidence/credibility intervals. An alternative is to infer phylogeny and times jointly to incorporate phylogenetic errors into molecular dating. We compared the performance of these two alternatives in reconstructing evolutionary timetrees using computer-simulated and empirical datasets. We found sequential and joint analyses to produce similar divergence times and phylogenetic relationships, except for some nodes in particular cases. The joint inference performed better when the phylogeny was not well resolved, situations in which the joint inference should be preferred. However, joint inference can be infeasible for large datasets because available Bayesian methods are computationally burdensome. We present an alternative approach for joint inference that combines the bag of little bootstraps, maximum likelihood, and RelTime approaches for simultaneously inferring evolutionary relationships, divergence times, and confidence intervals, incorporating phylogeny uncertainty. The new method alleviates the high computational burden imposed by Bayesian methods while achieving a similar result.

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.

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Diversification Models Conflate Likelihood and Prior, and Cannot be Compared Using Conventional Model-Comparison Tools

Time-calibrated phylogenetic trees are a tremendously powerful tool for studying evolutionary, ecological, and epidemiological phenomena. Such trees are predominantly inferred in a Bayesian framework, with the phylogeny itself treated as a parameter with a prior distribution (a "tree prior"). However, we show that the tree "parameter" consists, in part, of data, in the form of taxon samples. Treating the tree as a parameter fails to account for these data and compromises our ability to compare among models using standard techniques (e.g., marginal likelihoods estimated using path-sampling and stepping-stone sampling algorithms). Since accuracy of the inferred phylogeny strongly depends on how well the tree prior approximates the true diversification process that gave rise to the tree, the inability to accurately compare competing tree priors has broad implications for applications based on time-calibrated trees. We outline potential remedies to this problem, and provide guidance for researchers interested in assessing the fit of tree models. [Bayes factors; Bayesian model comparison; birth-death models; divergence-time estimation; lineage diversification].

From fossils to mind

Fossil endocasts record features of brains from the past: size, shape, vasculature, and gyrification. These data, alongside experimental and comparative evidence, are needed to resolve questions about brain energetics, cognitive specializations, and developmental plasticity. Through the application of interdisciplinary techniques to the fossil record, paleoneurology has been leading major innovations. Neuroimaging is shedding light on fossil brain organization and behaviors. Inferences about the development and physiology of the brains of extinct species can be experimentally investigated through brain organoids and transgenic models based on ancient DNA. Phylogenetic comparative methods integrate data across species and associate genotypes to phenotypes, and brains to behaviors. Meanwhile, fossil and archeological discoveries continuously contribute new knowledge. Through cooperation, the scientific community can accelerate knowledge acquisition. Sharing digitized museum collections improves the availability of rare fossils and artifacts. Comparative neuroanatomical data are available through online databases, along with tools for their measurement and analysis. In the context of these advances, the paleoneurological record provides ample opportunity for future research. Biomedical and ecological sciences can benefit from paleoneurology's approach to understanding the mind as well as its novel research pipelines that establish connections between neuroanatomy, genes and behavior.

Biogeography and evolution of social parasitism in Australian Myrmecia bulldog ants revealed by phylogenomics

Studying the historical biogeography and life history transitions from eusocial colony life to social parasitism contributes to our understanding of the evolutionary mechanisms generating biodiversity in eusocial insects. The ants in the genus Myrmecia are a well-suited system for testing evolutionary hypotheses about how their species diversity was assembled through time because the genus is endemic to Australia with the single exception of the species M. apicalis inhabiting the Pacific Island of New Caledonia, and because at least one social parasite species exists in the genus. However, the evolutionary mechanisms underlying the disjunct biogeographic distribution of M. apicalis and the life history transition(s) to social parasitism remain unexplored. To study the biogeographic origin of the isolated, oceanic species M. apicalis and to reveal the origin and evolution of social parasitism in the genus, we reconstructed a comprehensive phylogeny of the ant subfamily Myrmeciinae. We utilized Ultra Conserved Elements (UCEs) as molecular markers to generate a comprehensive molecular genetic dataset consisting of 2,287 loci per taxon on average for 66 out of the 93 known Myrmecia species as well as for the sister lineage Nothomyrmecia macrops and selected outgroups. Our time-calibrated phylogeny inferred that: (i) stem Myrmeciinae originated during the Paleocene ∼58 Ma ago; (ii) the current disjunct biogeographic distribution of M. apicalis was driven by long-distance dispersal from Australia to New Caledonia during the Miocene ∼14 Ma ago; (iii) the single social parasite species, M. inquilina, evolved directly from one of the two known host species, M. nigriceps, in sympatry via the intraspecific route of social parasite evolution; and (iv) 5 of the 9 previously established taxonomic species groups are non-monophyletic. We suggest minor changes to reconcile the molecular phylogenetic results with the taxonomic classification. Our study enhances our understanding of the evolution and biogeography of Australian bulldog ants, contributes to our knowledge about the evolution of social parasitism in ants, and provides a solid phylogenetic foundation for future inquiries into the biology, taxonomy, and classification of Myrmeciinae.

Comparative genomics reveals unique features of two Babesia motasi subspecies: Babesia motasi lintanensis and Babesia motasi hebeiensis

Parasites of the Babesia genus are prevalent worldwide and infect a wide diversity of domestic animals and humans. Herein, using Oxford Nanopore Technology and Illumina sequencing technologies, we sequenced two Babesia sub-species, Babesia motasi lintanensis and Babesia motasi hebeiensis. We identified 3,815 one-to-one ortholog genes that are specific to ovine Babesia spp. Phylogenetic analysis reveals that the two B. motasi subspecies form a distinct clade from other Piroplasma spp. Consistent with their phylogenetic position, comparative genomic analysis reveals that these two ovine Babesia spp. share higher colinearity with Babesia bovis than with Babesia microti. Concerning the speciation date, B. m. lintanensis split from B. m. hebeiensis approximately 17 million years ago. Genes correlated to transcription, translation, protein modification and degradation, as well as differential/specialized gene family expansions in these two subspecies may favor adaptation to vertebrate and tick hosts. The close relationship between B. m. lintanensis and B. m. hebeiensis is underlined by a high degree of genomic synteny. Compositions of most invasion, virulence, development, and gene transcript regulation-related multigene families, including spherical body protein, variant erythrocyte surface antigen, glycosylphosphatidylinositol anchored proteins, and transcription factor Apetala 2 genes, is largely conserved, but in contrast to this conserved situation, we observe major differences in species-specific genes that may be involved in multiple functions in parasite biology. For the first time in Babesia spp., we find abundant fragments of long terminal repeat-retrotransposons in these two species. We provide fundamental information to characterize the genomes of B. m. lintanensis and B. m. hebeiensis, providing insights into the evolution of B. motasi group parasites.

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.

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

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

Renewal of planktonic foraminifera diversity after the Cretaceous Paleogene mass extinction by benthic colonizers

The biotic crisis following the end-Cretaceous asteroid impact resulted in a dramatic renewal of pelagic biodiversity. Considering the severe and immediate effect of the asteroid impact on the pelagic environment, it is remarkable that some of the most affected pelagic groups, like the planktonic foraminifera, survived at all. Here we queried a surface ocean metabarcoding dataset to show that calcareous benthic foraminifera of the clade Globothalamea are able to disperse actively in the plankton, and we show using molecular clock phylogeny that the modern planktonic clades originated from different benthic ancestors that colonized the plankton after the end-Cretaceous crisis. We conclude that the diversity of planktonic foraminifera has been the result of a constant leakage of benthic foraminifera diversity into the plankton, continuously refueling the planktonic niche, and challenge the classical interpretation of the fossil record that suggests that Mesozoic planktonic foraminifera gave rise to the modern communities.

Unveiling new perspective of phylogeography, genetic diversity, and population dynamics of Southeast Asian and Pacific chickens

The complex geographic and temporal origins of chicken domestication have attracted wide interest in molecular phylogeny and phylogeographic studies as they continue to be debated up to this day. In particular, the population dynamics and lineage-specific divergence time estimates of chickens in Southeast Asia (SEA) and the Pacific region are not well studied. Here, we analyzed 519 complete mitochondrial DNA control region sequences and identified 133 haplotypes with 70 variable sites. We documented 82.7% geographically unique haplotypes distributed across major haplogroups except for haplogroup C, suggesting high polymorphism among studied individuals. Mainland SEA (MSEA) chickens have higher overall genetic diversity than island SEA (ISEA) chickens. Phylogenetic trees and median-joining network revealed evidence of a new divergent matrilineage (i.e., haplogroup V) as a sister-clade of haplogroup C. The maximum clade credibility tree estimated the earlier coalescence age of ancestral D-lineage (i.e., sub-haplogroup D2) of continental chickens (3.7 kya; 95% HPD 1985-4835 years) while island populations diverged later at 2.1 kya (95% HPD 1467-2815 years). This evidence of earlier coalescence age of haplogroup D ancestral matriline exemplified dispersal patterns to the ISEA, and thereafter the island clade diversified as a distinct group.

Phylogenomic Evidence for the Origin of Obligate Anaerobic Anammox Bacteria Around the Great Oxidation Event

The anaerobic ammonium oxidation (anammox) bacteria can transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remain enigmatic. Here, we performed a comprehensive phylogenomic and molecular clock analysis of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which include implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based molecular dating approach, we estimated a consistent origin of anammox bacteria at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32-2.5 Ga) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress adaptation, bioenergetics, and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.

Bayesian phylogenetic inference using relaxed-clocks and the multispecies coalescent

The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.

Common Ground between Biological Rhythms and Forensics

Simple Summary

Biological clocks regulate the timing of numerous body functions in adaption to daily repeating cycles in the environment, such as the sleep–wake phases that are trained by the cycling changes of night and day light. The identification of a deceased victim is a critical component in a forensic investigation, but it can be significantly hampered by the condition of the dead body and the lack of personal records and documents. This review links current knowledge on the molecular mechanisms of biological rhythms to forensically relevant aspects, including the time period since death, cause of death, the use of insects for forensics, sex and age of a person, ethnic background and development. Putting these findings in context demonstrates how the analysis of molecular clock analysis could be used as tool for future personal identification in forensic investigations.

Abstract

Biological clocks set the timing for a large number of essential processes in the living human organism. After death, scientific evidence is required in forensic investigations in order to collect as much information as possible on the death circumstances and personal identifiers of the deceased victim. We summarize the associations between the molecular mechanisms of biological rhythms and forensically relevant aspects, including post-mortem interval and cause of death, entomological findings, sex, age, ethnicity and development. Given their importance during lifetime, biological rhythms could be potential tools to draw conclusions on the death circumstances and the identity of a deceased person by mechanistic investigations of the different biological clocks in a forensic context. This review puts the known effects of biological rhythms on the functions of the human organism in context with potential applications in forensic fields of interest, such as personal identification, entomology as well as the determination of the post-mortem interval and cause of death.

Comparative genomic analysis of Babesia duncani responsible for human babesiosis

Background

Human babesiosis, caused by parasites of the genus Babesia, is an emerging and re-emerging tick-borne disease that is mainly transmitted by tick bites and infected blood transfusion. Babesia duncani has caused majority of human babesiosis in Canada; however, limited data are available to correlate its genomic information and biological features.

Results

We generated a B. duncani reference genome using Oxford Nanopore Technology (ONT) and Illumina sequencing technology and uncovered its biological features and phylogenetic relationship with other Apicomplexa parasites. Phylogenetic analyses revealed that B. duncani form a clade distinct from B. microti, Babesia spp. infective to bovine and ovine species, and Theileria spp. infective to bovines. We identified the largest species-specific gene family that could be applied as diagnostic markers for this pathogen. In addition, two gene families show signals of significant expansion and several genes that present signatures of positive selection in B. duncani, suggesting their possible roles in the capability of this parasite to infect humans or tick vectors.

Conclusions

Using ONT sequencing and Illumina sequencing technologies, we provide the first B. duncani reference genome and confirm that B. duncani forms a phylogenetically distinct clade from other Piroplasm parasites. Comparative genomic analyses show that two gene families are significantly expanded in B. duncani and may play important roles in host cell invasion and virulence of B. duncani. Our study provides basic information for further exploring B. duncani features, such as host-parasite and tick-parasite interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01361-9.

Studying speciation and extinction dynamics from phylogenies: addressing identifiability issues

A lot of what we know about past speciation and extinction dynamics is based on statistically fitting birth-death processes to phylogenies of extant species. Despite their wide use, the reliability of these tools is regularly questioned. It was recently demonstrated that vast 'congruent' sets of alternative diversification histories cannot be distinguished (i.e., are not identifiable) using extant phylogenies alone, reanimating the debate about the limits of phylogenetic diversification analysis. Here, we summarize what we know about the identifiability of the birth-death process and how identifiability issues can be addressed. We conclude that extant phylogenies, when combined with appropriate prior hypotheses and regularization techniques, can still tell us a lot about past diversification dynamics.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Why Plasmodium vivax and Plasmodium falciparum are so different? A tale of two clades and their species diversities

The global malaria burden sometimes obscures that the genus Plasmodium comprises diverse clades with lineages that independently gave origin to the extant human parasites. Indeed, the differences between the human malaria parasites were highlighted in the classical taxonomy by dividing them into two subgenera, the subgenus Plasmodium, which included all the human parasites but Plasmodium falciparum that was placed in its separate subgenus, Laverania. Here, the evolution of Plasmodium in primates will be discussed in terms of their species diversity and some of their distinct phenotypes, putative molecular adaptations, and host–parasite biocenosis. Thus, in addition to a current phylogeny using genome-level data, some specific molecular features will be discussed as examples of how these parasites have diverged. The two subgenera of malaria parasites found in primates, Plasmodium and Laverania, reflect extant monophyletic groups that originated in Africa. However, the subgenus Plasmodium involves species in Southeast Asia that were likely the result of adaptive radiation. Such events led to the Plasmodium vivax lineage. Although the Laverania species, including P. falciparum, has been considered to share “avian characteristics,” molecular traits that were likely in the common ancestor of primate and avian parasites are sometimes kept in the Plasmodium subgenus while being lost in Laverania. Assessing how molecular traits in the primate malaria clades originated is a fundamental science problem that will likely provide new targets for interventions. However, given that the genus Plasmodium is paraphyletic (some descendant groups are in other genera), understanding the evolution of malaria parasites will benefit from studying “non-Plasmodium” Haemosporida.

Predominantly Eastward Long-Distance Dispersal in Pantropical Ochnaceae Inferred From Ancestral Range Estimation and Phylogenomics

Ochnaceae is a pantropical family with multiple transoceanic disjunctions at deep and shallow levels. Earlier attempts to unravel the processes that led to such biogeographic patterns suffered from insufficient phylogenetic resolution and unclear delimitation of some of the genera. In the present study, we estimated divergence time and ancestral ranges based on a phylogenomic framework with a well-resolved phylogenetic backbone to tackle issues of the timing and direction of dispersal that may explain the modern global distribution of Ochnaceae. The nuclear data provided the more robust framework for divergence time estimation compared to the plastome-scale data, although differences in the inferred clade ages were mostly small. While Ochnaceae most likely originated in West Gondwana during the Late Cretaceous, all crown-group disjunctions are inferred as dispersal-based, most of them as transoceanic long-distance dispersal (LDD) during the Cenozoic. All LDDs occurred in an eastward direction except for the SE Asian clade of Sauvagesieae, which was founded by trans-Pacific dispersal from South America. The most species-rich clade by far, Ochninae, originated from either a widespread neotropical-African ancestor or a solely neotropical ancestor which then dispersed to Africa. The ancestors of this clade then diversified in Africa, followed by subsequent dispersal to the Malagasy region and tropical Asia on multiple instances in three genera during the Miocene-Pliocene. In particular, Ochna might have used the South Arabian land corridor to reach South Asia. Thus, the pantropical distribution of Ochnaceae is the result of LDD either transoceanic or via land bridges/corridors, whereas vicariance might have played a role only along the stem of the family.

A species-level timeline of mammal evolution integrating phylogenomic data

High-throughput sequencing projects generate genome-scale sequence data for species-level phylogenies1-3. However, state-of-the-art Bayesian methods for inferring timetrees are computationally limited to small datasets and cannot exploit the growing number of available genomes4. In the case of mammals, molecular-clock analyses of limited datasets have produced conflicting estimates of clade ages with large uncertainties5,6, and thus the timescale of placental mammal evolution remains contentious7-10. Here we develop a Bayesian molecular-clock dating approach to estimate a timetree of 4,705 mammal species integrating information from 72 mammal genomes. We show that increasingly larger phylogenomic datasets produce diversification time estimates with progressively smaller uncertainties, facilitating precise tests of macroevolutionary hypotheses. For example, we confidently reject an explosive model of placental mammal origination in the Palaeogene8 and show that crown Placentalia originated in the Late Cretaceous with unambiguous ordinal diversification in the Palaeocene/Eocene. Our Bayesian methodology facilitates analysis of complete genomes and thousands of species within an integrated framework, making it possible to address hitherto intractable research questions on species diversifications. This approach can be used to address other contentious cases of animal and plant diversifications that require analysis of species-level phylogenomic datasets.

No link between population isolation and speciation rate in squamate reptiles

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Rooting Species Trees Using Gene Tree-Species Tree Reconciliation

Interpreting phylogenetic trees requires a root, which provides the direction of evolution and polarizes ancestor-descendant relationships. But inferring the root using genetic data is difficult, particularly in cases where the closest available outgroup is only distantly related, which are common for microbes. In this chapter, we present a workflow for estimating rooted species trees and the evolutionary history of the gene families that evolve within them using probabilistic gene tree-species tree reconciliation. We illustrate the pipeline using a small dataset of prokaryotic genomes, for which the example scripts can be run using modest computer resources. We describe the rooting method used in this work in the context or other rooting strategies and discuss some of the limitations and opportunities presented by probabilistic gene tree-species tree reconciliation methods.

Dating Microbial Evolution with MCMCtree

This protocol explains how to use the program MCMCtree to estimate divergence times in microbial phylogenies. The main advantage of MCMCtree is the implementation of an approximation to the molecular data likelihood that dramatically speeds up computation during Bayesian MCMC sampling of divergence times and evolutionary rates. The approximation allows the analysis of large phylogenies with hundreds of taxa and molecular alignments with thousands or millions of sites. Two examples are used to illustrate Bayesian clock dating with MCMCtree. The first is a phylogeny of (mostly) microbial eukaryotes and prokaryotes encompassing the major groups of life on Earth, and for which fossil information, to calibrate the nodes of the phylogeny, is available. The second is a phylogeny of influenza viruses with known sampling times. An overview of Bayesian MCMC sampling is given as well as practical advice on issues such as construction of the time and rate prior and assessment of convergence of MCMC chains. Strategies for estimating times in microbial phylogenies for which neither fossil information nor sampling times are known are discussed.

Estimating the Divergence Times of Alphaproteobacteria Based on Mitochondrial Endosymbiosis and Eukaryotic Fossils

Alphaproteobacteria is one of the most abundant bacterial lineages that successfully colonize diverse marine and terrestrial environments on Earth. In addition, many alphaproteobacterial lineages have established close association with eukaryotes. This makes Alphaproteobacteria a promising system to test the link between the emergence of ecologically important bacteria and related geological events and the co-evolution between symbiotic bacteria and their hosts. Understanding the timescale of evolution of Alphaproteobacteria is key to testing these hypotheses, which is limited by the scarcity of bacterial fossils, however. Based on the mitochondrial endosymbiosis which posits that the mitochondrion originated from an alphaproteobacterial lineage, we propose a new strategy to estimate the divergence times of lineages within the Alphaproteobacteria by leveraging the fossil records of eukaryotes. In this chapter, we describe the workflow of the mitochondria-based method to date Alphaproteobacteria evolution by detailing the software, methods, and commands used for each step. Visualization of data and results is also described. We also provide related notes with background information and alternative options. All codes used to build this protocol are made available to the public, and we strive to make this protocol user-friendly in particular to microbiologists with limited practical skills in bioinformatics.

A Computational Protocol for Dating the Evolution of Cyanobacteria

Cyanobacteria are known to play important roles in driving biological and geochemical innovations in ancient Earth. The origin of Cyanobacteria is the key to understanding these evolutionary events and thus has gained much interest to biologists and geobiologists. Recent development of the molecular dating approaches provides us an opportunity to assess the timeline of Cyanobacteria evolution based on relaxed clock models. The implementation of Bayesian phylogenetic approaches accommodates the uncertainties from different sources, such as fossil calibrations and topological structure of the phylogenomic tree, and provides us converged estimates of posterior mean ages. In this chapter, by taking Cyanobacteria as an example, we introduce a refined strategy to perform molecular dating analysis, as well as a practical method to evaluate the precision of dating analysis.

Progress and Prospects in Epigenetic Studies of Ancient DNA

Development of technologies for high-throughput whole-genome sequencing and improvement of sample preparation techniques made it possible to study ancient DNA (aDNA) from archaeological samples over a million year old. The studies of aDNA have shed light on the history of human migration, replacement of populations, interbreeding of Cro-Magnons with Neanderthals and Denisovans, evolution of human pathogens, etc. Equally important is the possibility to investigate epigenetic modifications of ancient genomes, which has allowed to obtain previously inaccessible information on gene expression, nucleosome positioning, and DNA methylation. Analysis of methylation status of certain genomic sites can predict an individual's age at death and reconstruct some phenotypic features, as it was done for the Denisovan genome, and even to elucidate unfavorable environmental factors that had affected this archaic individual. In this review, we discuss current progress in epigenetic studies of aDNA, including methodological approaches and promising research directions in this field.

Unexpected cryptic species among streptophyte algae most distant to land plants

Streptophytes are one of the major groups of the green lineage (Chloroplastida or Viridiplantae). During one billion years of evolution, streptophytes have radiated into an astounding diversity of uni- and multicellular green algae as well as land plants. Most divergent from land plants is a clade formed by Mesostigmatophyceae, Spirotaenia spp. and Chlorokybophyceae. All three lineages are species-poor and the Chlorokybophyceae consist of a single described species, Chlorokybus atmophyticus. In this study, we used phylogenomic analyses to shed light into the diversity within Chlorokybus using a sampling of isolates across its known distribution. We uncovered a consistent deep genetic structure within the Chlorokybus isolates, which prompted us to formally extend the Chlorokybophyceae by describing four new species. Gene expression differences among Chlorokybus species suggest certain constitutive variability that might influence their response to environmental factors. Failure to account for this diversity can hamper comparative genomic studies aiming to understand the evolution of stress response across streptophytes. Our data highlight that future studies on the evolution of plant form and function can tap into an unknown diversity at key deep branches of the streptophytes.

Assessing rapid relaxed-clock methods for phylogenomic dating

Rapid relaxed-clock dating methods are frequently applied to analyze phylogenomic data sets containing hundreds to thousands of sequences because of their accuracy and computational efficiency. However, the relative performance of different rapid dating methods is yet to be compared on the same data sets, and, thus, the power and pitfalls of selecting among these approaches remain unclear. We compared the accuracy, bias, and coverage probabilities of RelTime, treePL, and least-squares dating time estimates by applying them to analyze computer-simulated data sets in which evolutionary rates varied extensively among branches in the phylogeny. RelTime estimates were consistently more accurate than the other two, particularly when evolutionary rates were autocorrelated or shifted convergently among lineages. The 95% confidence intervals (CIs) around RelTime dates showed appropriate coverage probabilities (95% on average), but other methods produced rather low coverage probabilities because of overly narrow CIs of time estimates. Overall, RelTime appears to be a more efficient method for estimating divergence times for large phylogenies.

Relative time constraints improve molecular dating

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