Expectation-Maximization enables Phylogenetic Dating under a Categorical Rate Model.

Dating phylogenetic trees to obtain branch lengths in time unit is essential for many downstream applications but has remained challenging. Dating requires inferring substitution rates that can change across the tree. While we can assume to have information about a small subset of nodes from the fossil record or sampling times (for fast-evolving organisms), inferring the ages of the other nodes essentially requires extrapolation and interpolation. Assuming a distribution of branch rates, we can formulate dating as a constrained maximum likelihood (ML) estimation problem. While ML dating methods exist, their accuracy degrades in the face of model misspecification where the assumed parametric statistical distribution of branch rates vastly differs from the true distribution. Notably, most existing methods assume rigid, often unimodal, branch rate distributions. A second challenge is that the likelihood function involves an integral over the continuous domain of the rates and often leads to difficult non-convex optimization problems. To tackle these two challenges, we propose a new method called Molecular Dating using Categorical-models (MD-Cat). MD-Cat uses a categorical model of rates inspired by non-parametric statistics and can approximate a large family of models by discretizing the rate distribution into k categories. Under this model, we can use the Expectation- Maximization (EM) algorithm to co-estimate rate categories and branch lengths in time units. Our model has fewer assumptions about the true distribution of branch rates than parametric models such as Gamma or LogNormal distribution. Our results on two simulated and real datasets of Angiosperms and HIV and a wide selection of rate distributions show that MD-Cat is often more accurate than the alternatives, especially on datasets with exponential or multimodal rate distributions.

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.

The rediscovery of a relict unlocks the first global phylogeny of whip spiders (Amblypygi).

Asymmetrical rates of cladogenesis and extinction abound in the Tree of Life, resulting in numerous minute clades that are dwarfed by larger sister groups. Such taxa are commonly regarded as phylogenetic relicts or "living fossils" when they exhibit an ancient first appearance in the fossil record and prolonged external morphological stasis, particularly in comparison to their more diversified sister groups. Due to their special status, various phylogenetic relicts tend to be well-studied and prioritized for conservation. A notable exception to this trend is found within Amblypygi ("whip spiders"), a visually striking order of functionally hexapodous arachnids that are notable for their antenniform first walking leg pair (the eponymous "whips"). Paleoamblypygi, the putative sister group to the remaining Amblypygi, is known from Late Carboniferous and Eocene deposits, but is survived by a single living species, Paracharon caecus Hansen, 1921, that was last collected in 1899. Due to the absence of genomic sequence-grade tissue for this vital taxon, there is no global molecular phylogeny for Amblypygi to date, nor a fossil-calibrated estimation of divergences within the group. Here, we report a previously unknown species of Paleoamblypygi from a cave site in Colombia. Capitalizing upon this discovery, we generated the first molecular phylogeny of Amblypygi, integrating ultraconserved element sequencing with legacy Sanger datasets and including described extant genera. To quantify the impact of sampling Paleoamblypygi on divergence time estimation, we performed in silico experiments with pruning of Paracharon. We demonstrate that the omission of relicts has a significant impact on the accuracy of node dating approaches that outweighs the impact of excluding ingroup fossils, which bears upon the ancestral range reconstruction for the group. Our results underscore the imperative for biodiversity discovery efforts in elucidating the phylogenetic relationships of "dark taxa", and especially phylogenetic relicts in tropical and subtropical habitats. The lack of reciprocal monophyly for Charontidae and Charinidae leads us to subsume them into one family, Charontidae, new synonymy.

The Neotropical endemic liverwort subfamily Micropterygioideae had circum-Antarctic links to the rest of the Lepidoziaceae during the early Cretaceous.

Lepidoziaceae are the third-largest family of liverworts, with about 860 species distributed on all continents. The evolutionary history of this family has not been satisfactorily resolved, with taxa such as Micropterygioideae yet to be included in phylogenetic analyses. We inferred a dated phylogeny of Lepidoziaceae using a data set consisting of 13 genetic markers, sampled from 147 species. Based on our phylogenetic estimate, we used statistical dispersal-vicariance analysis to reconstruct the biogeographic history of the family. We inferred a crown age of 197 Ma (95% credible interval 157-240 Ma) for the family in the Australian region, with most major lineages also originating in the same region. Micropterygioideae are placed as the sister group to Lembidioideae, with these two lineages diverging from each other about 132 Ma in the South American-Australian region. With South America and Australia being connected through Antarctica at the time, our results suggest a circum-Antarctic link between Micropterygioideae and the rest of the family. Crown Micropterygioideae were inferred to have arisen 45 Ma in South America before the continent separated from Antarctica. Extinction from southern temperate regions might explain the present-day restriction of Micropterygioideae to the Neotropics. Our study reveals the influence of past geological events, such as continental drift, on the evolution and distribution of a widespread and diverse family of liverworts.

Phylogeny, biogeography and ecological diversification of New Caledonian palms (Arecaceae).

BACKGROUND AND AIMS: The geographical origin and evolutionary mechanisms underpinning the rich and distinctive New Caledonian flora remain poorly understood. This is attributable to the complex geological past of the island and to the scarcity of well-resolved species-level phylogenies. Here, we infer phylogenetic relationships and divergence times of New Caledonian palms, which comprise 40 species. We use this framework to elucidate the biogeography of New Caledonian palm lineages and to explore how extant species might have formed. METHODS: A phylogenetic tree including 37 New Caledonian palm species and 77 relatives from tribe Areceae was inferred from 151 nuclear genes obtained by targeted sequencing. Fossil-calibrated divergence times were estimated and ancestral ranges inferred. Ancestral and extant ecological preferences in terms of elevation, precipitation and substrate were compared between New Caledonian sister species to explore their possible roles as drivers of speciation. KEY RESULTS: New Caledonian palms form four well-supported clades, inside which relationships are well resolved. Our results support the current classification but suggest that Veillonia and Campecarpus should be resurrected and fail to clarify whether Rhopalostylidinae is sister to or nested in Basseliniinae. New Caledonian palm lineages are derived from New Guinean and Australian ancestors, which reached the island through at least three independent dispersal events between the Eocene and Miocene. Palms then dispersed out of New Caledonia at least five times, mainly towards Pacific islands. Geographical and ecological transitions associated with speciation events differed across time and genera. Substrate transitions were more frequently associated with older events than with younger ones. CONCLUSIONS: Neighbouring areas and a mosaic of local habitats shaped the palm flora of New Caledonia, and the island played a significant role in generating palm diversity across the Pacific region. This new spatio-temporal framework will enable population-level ecological and genetic studies to unpick the mechanisms underpinning New Caledonian palm endemism.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates.

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analysed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all thirteen major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Under the hood: Phylogenomics of hooded tick spiders (Arachnida, Ricinulei) uncovers discordance between morphology and molecules.

Ricinulei or hooded tick-spiders are a cryptic and ancient group of arachnids. The order consists of around 100 highly endemic extant species restricted to the Afrotropics and the Neotropics along with 22 fossil species. Their antiquity and low vagility make them an excellent group with which to interrogate biogeographic questions. To date, only four molecular analyses have been conducted on the group and they failed to resolve the relationships of the main lineages and even recovering the non-monophyly of the three genera. These studies were limited to a few Sanger loci or phylogenomic analyses with at most seven ingroup samples. To increase phylogenetic resolution in this little-understood and poorly studied group, we present the most comprehensive phylogenomic study of Ricinulei to date leveraging the Arachnida ultra-conserved element probe set. With a data set of 473 loci across 96 ingroup samples, analyses resolved a monophyletic Neotropical clade consisting of four main lineages. Two of them correspond to the current genera Cryptocellus and Pseudocellus while topology testing revealed one lineage to likely be a phylogenetic reconstruction artefact. The fourth lineage, restricted to Northwestern, Andean South America, is consistent with the Cryptocellus magnus group, likely corresponding to the historical genus Heteroricinoides. Since we did not sample the type species for this old genus, we do not formally re-erect Heteroricinoides but our data suggest the need for a thorough morphological re-examination of Neotropical Ricinulei.

Molecular phylogeny and divergence time of Harpalyce (Leguminosae, Papilionoideae), a lineage with amphitropical diversification in seasonally dry forests and savannas.

Our knowledge of the systematics of the papilionoid legume tribe Brongniartieae has greatly benefitted from recent advances in molecular phylogenetics. The tribe was initially described to include species marked by a strongly bilabiate calyx and an embryo with a straight radicle, but recent research has placed taxa from the distantly related core Sophoreae and Millettieae within it. Despite these advances, the most species-rich genera within the Brongniartieae are still not well studied, and their morphological and biogeographical evolution remains poorly understood. Comprising 35 species, Harpalyce is one of these poorly studied genera. In this study, we present a comprehensive, multi-locus molecular phylogeny of the Brongniartieae, with an increased sampling of Harpalyce, to investigate morphological and biogeographical evolution within the group. Our results confirm the monophyly of Harpalyce and indicate that peltate glandular trichomes and a strongly bilabiate calyx with a carinal lip and three fused lobes are synapomorphies for the genus, which is internally divided into three distinct ecologically and geographically divergent lineages, corresponding to the previously recognized sections. Our biogeographical reconstructions demonstrate that Brongniartieae originated in South America during the Eocene, with subsequent pulses of diversification in South America, Mesoamerica, and Australia. Harpalyce also originated in South America during the Miocene at around 20 Ma, with almost synchronous later diversification in South America and Mexico/Mesoamerica beginning 10 Ma, but mostly during the Pliocene. Migration of Harpalyce from South to North America was accompanied by a biome and ecological shift from savanna to seasonally dry forest.

The Global Evolutionary History of Orf Virus in Sheep and Goats Revealed by Whole Genomes Data.

Orf virus (ORFV) belongs to the genus Parapoxvirus (Poxviridae family). It is the causative agent of contagious ecthyma (CE) that is an economically detrimental disease affecting small ruminants globally. Contagious ecthyma outbreaks are usually reported in intensive breeding of sheep and goats but they have also been reported in wildlife species. Notably, ORFV can infect humans, leading to a zoonotic disease. This study aims to elucidate the global evolutionary history of ORFV genomes in sheep and goats, including the first genomes from Central America in the analyses. In comparison to the last study on ORFV whole genomes, the database now includes 11 more sheep and goat genomes, representing an increase of 42%. The analysis of such a broader database made it possible to obtain a fine molecular dating of the coalescent time for ORFV S and G genomes, further highlighting the genetic structuring between sheep and goat genomes and corroborating their emergence in the latter half of 20th century.

The replication-competent HIV reservoir is a genetically restricted, younger subset of the overall pool of HIV proviruses persisting during therapy, which is highly genetically stable over time.

Within-host HIV populations continually diversify during untreated infection, and this diversity persists within infected cell reservoirs during antiretroviral therapy (ART). Achieving a better understanding of on-ART proviral evolutionary dynamics, and a better appreciation of how the overall persisting pool of (largely genetically defective) proviruses differs from the much smaller replication-competent HIV reservoir, is critical to HIV cure efforts. We reconstructed within-host HIV evolutionary histories in blood from seven participants of the Women's Interagency HIV Study who experienced HIV seroconversion, and used these data to characterize the diversity, lineage origins, and ages of proviral env-gp120 sequences sampled longitudinally up to 12 years on ART. We also studied HIV sequences emerging from the reservoir in two participants. We observed that proviral clonality generally increased over time on ART, with clones frequently persisting long term. While on-ART proviral integration dates generally spanned the duration of untreated infection, HIV emerging in plasma was exclusively younger (i.e., dated to the years immediately pre-ART). The genetic and age distributions of distinct proviral sequences remained stable during ART in all but one participant, in whom there was evidence that younger proviruses had been preferentially eliminated after 12 years on ART. Analysis of the gag region in three participants corroborated our env-gp120-based observations, indicating that our observations are not influenced by the HIV region studied. Our results underscore the remarkable genetic stability of the distinct proviral sequences that persist in blood during ART. Our results also suggest that the replication-competent HIV reservoir is a genetically restricted, younger subset of this overall proviral pool.IMPORTANCECharacterizing the genetically diverse HIV sequences that persist in the reservoir despite antiretroviral therapy (ART) is critical to cure efforts. Our observations confirm that proviruses persisting in blood on ART, which are largely genetically defective, broadly reflect the extent of within-host HIV evolution pre-ART. Moreover, on-ART clonal expansion is not appreciably accompanied by the loss of distinct proviral lineages. In fact, on-ART proviral genetic composition remained stable in all but one participant, in whom, after 12 years on ART, proviruses dating to around near ART initiation had been preferentially eliminated. We also identified recombinant proviruses between parental sequence fragments of different ages. Though rare, such sequences suggest that reservoir cells can be superinfected with HIV from another infection era. Overall, our finding that the replication-competent reservoir in blood is a genetically restricted, younger subset of all persisting proviruses suggests that HIV cure strategies will need to eliminate a reservoir that differs in key respects from the overall proviral pool.

A timeline of bacterial and archaeal diversification in the ocean.

Microbial plankton play a central role in marine biogeochemical cycles, but the timing in which abundant lineages diversified into ocean environments remains unclear. Here, we reconstructed the timeline in which major clades of bacteria and archaea colonized the ocean using a high-resolution benchmarked phylogenetic tree that allows for simultaneous and direct comparison of the ages of multiple divergent lineages. Our findings show that the diversification of the most prevalent marine clades spans throughout a period of 2.2 Ga, with most clades colonizing the ocean during the last 800 million years. The oldest clades - SAR202, SAR324, Ca. Marinimicrobia, and Marine Group II - diversified around the time of the Great Oxidation Event, during which oxygen concentration increased but remained at microaerophilic levels throughout the Mid-Proterozoic, consistent with the prevalence of some clades within these groups in oxygen minimum zones today. We found the diversification of the prevalent heterotrophic marine clades SAR11, SAR116, SAR92, SAR86, and Roseobacter as well as the Marine Group I to occur near to the Neoproterozoic Oxygenation Event (0.8-0.4 Ga). The diversification of these clades is concomitant with an overall increase of oxygen and nutrients in the ocean at this time, as well as the diversification of eukaryotic algae, consistent with the previous hypothesis that the diversification of heterotrophic bacteria is linked to the emergence of large eukaryotic phytoplankton. The youngest clades correspond to the widespread phototrophic clades Prochlorococcus, Synechococcus, and Crocosphaera, whose diversification happened after the Phanerozoic Oxidation Event (0.45-0.4 Ga), in which oxygen concentrations had already reached their modern levels in the atmosphere and the ocean. Our work clarifies the timing at which abundant lineages of bacteria and archaea colonized the ocean, thereby providing key insights into the evolutionary history of lineages that comprise the majority of prokaryotic biomass in the modern ocean.

HIV reservoirs are dominated by genetically younger and clonally enriched proviruses.

IMPORTANCE: Characterizing the human immunodeficiency virus (HIV) reservoir that endures despite antiretroviral therapy (ART) is critical to cure efforts. We observed that the oldest proviruses persisting during ART were exclusively defective, while intact proviruses (and rebound HIV) dated to nearer ART initiation. This helps explain why studies that sampled sub-genomic proviruses on-ART (which are largely defective) routinely found sequences dating to early infection, whereas those that sampled replication-competent HIV found almost none. Together with our findings that intact proviruses were more likely to be clonal, and that on-ART low-level/isolated viremia originated from proviruses of varying ages (including possibly defective ones), our observations indicate that (i) on-ART and rebound viremia can have distinct within-host origins, (ii) intact proviruses have shorter lifespans than grossly defective ones and thus depend more heavily on clonal expansion for persistence, and (iii) an HIV reservoir predominantly "dating" to near ART initiation will be substantially adapted to within-host pressures, complicating immune-based cure strategies.

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.

Early Detection of the Recombinant SARS-CoV-2 XAN Variant in Bulgaria: Initial Genomic Insights into Yet Another Piece of the Growing Puzzle of Recombinant Clades.

The first recombinant SARS-CoV-2 variants were identified in 2022, causing public health concerns. The importance of recombinant variants has increased especially since the WHO designated the recombinant variant XBB and its lineages as subvariants that require monitoring on 20 November 2022. In this study, we provide the first insights into the new SARS-CoV-2 variant named XAN, a recombinant composed of Omicron sub-lineages BA.2 and BA.5. To our knowledge, this is the first report on the recombinant SARS-CoV-2 XAN variant identified in Bulgaria.

Inferring Human Immunodeficiency Virus 1 Proviral Integration Dates With Bayesian Inference.

Human immunodeficiency virus 1 (HIV) proviruses archived in the persistent reservoir currently pose the greatest obstacle to HIV cure due to their evasion of combined antiretroviral therapy and ability to reseed HIV infection. Understanding the dynamics of the HIV persistent reservoir is imperative for discovering a durable HIV cure. Here, we explore Bayesian methods using the software BEAST2 to estimate HIV proviral integration dates. We started with within-host longitudinal HIV sequences collected prior to therapy, along with sequences collected from the persistent reservoir during suppressive therapy. We built a BEAST2 model to estimate integration dates of proviral sequences collected during suppressive therapy, implementing a tip date random walker to adjust the sequence tip dates and a latency-specific prior to inform the dates. To validate our method, we implemented it on both simulated and empirical data sets. Consistent with previous studies, we found that proviral integration dates were spread throughout active infection. Path sampling to select an alternative prior for date estimation in place of the latency-specific prior produced unrealistic results in one empirical data set, whereas on another data set, the latency-specific prior was selected as best fitting. Our Bayesian method outperforms current date estimation techniques with a root mean squared error of 0.89 years on simulated data relative to 1.23-1.89 years with previously developed methods. Bayesian methods offer an adaptable framework for inferring proviral integration dates.

HIV reservoirs are dominated by genetically younger and clonally enriched proviruses.

In order to cure HIV, we need to better understand the within-host evolutionary origins of the small reservoir of genome-intact proviruses that persists within infected cells during antiretroviral therapy (ART). Most prior studies on reservoir evolutionary dynamics however did not discriminate genome-intact proviruses from the vast background of defective ones. We reconstructed within-host pre-ART HIV evolutionary histories in six individuals and leveraged this information to infer the ages of intact and defective proviruses sampled after an average >9 years on ART, along with the ages of rebound and low-level/isolated viremia occurring during this time. We observed that the longest-lived proviruses persisting on ART were exclusively defective, usually due to large deletions. In contrast, intact proviruses and rebound HIV exclusively dated to the years immediately preceding ART. These observations are consistent with genome-intact proviruses having shorter lifespans, likely due to the cumulative risk of elimination following viral reactivation and protein production. Consistent with this, intact proviruses (and those with packaging signal defects) were three times more likely to be genetically identical compared to other proviral types, highlighting clonal expansion as particularly important in ensuring their survival. By contrast, low-level/isolated viremia sequences were genetically heterogeneous and sometimes ancestral, where viremia may have originated from defective proviruses. Results reveal that the HIV reservoir is dominated by clonally-enriched and genetically younger sequences that date to the untreated infection period when viral populations had been under within-host selection pressures for the longest duration. Knowledge of these qualities may help focus strategies for reservoir elimination.

Genomic characterization and molecular dating of the novel bacterium Permianibacter aggregans HW001T, which originated from Permian ground water.

The Permian Basin is a unique ecosystem located in the southwest of the USA. An unanswered question is whether the bacteria in the Permian Basin adapted to the changing paleomarine environment and survived in the remnants of Permian groundwater. In our previous study, a novel bacterial strain, Permianibacter aggregans HW001T, was isolated from microalgae cultures incubated with Permian Basin waters, and was shown to originate from the Permian Ocean. In this study, strain HW001T was shown to be the representative strain of a novel family, classified as 'Permianibacteraceae'. The results of molecular dating suggested that the strain HW001T diverged ~ 447 million years ago (mya), which is the early Permian period (~ 250 mya). Genome analysis was used to access its potential energy utilization and biosynthesis capacity. A large number of transporters, carbohydrate-active enzymes and protein-degradation related genes have been annotated in the genome of strain HW001T. In addition, a series of important metabolic pathways, such as peptidoglycan biosynthesis, osmotic stress response system and multifunctional quorum sensing were annotated, which may confer the ability to adapt to various unfavorable environmental conditions. Finally, the evolutionary history of strain HW001T was reconstructed and the horizontal transfer of genes was predicted, indicating that the adaptation of P. aggregans to a changing marine environment depends on the evolution of their metabolic capabilities, especially in signal transmission. In conclusion, the results of this study provide genomic information for revealing the adaptive mechanism of strain HW001T to the changing ancient oceans. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00164-3.

Is it inappropriate to ask for your age? Evaluating parameter impact on tree dating in a challenging clade (Macroscelidea).

Sengis (order Macroscelidea) are small mammals endemic to Africa. The taxonomy and phylogeny of sengis has been difficult to resolve due to a lack of clear morphological apomorphies. Molecular phylogenies have already significantly revised sengi systematics, but until now no molecular phylogeny has included all 20 extant species. In addition, the age of origin of the sengi crown clade and the divergence age of its two extant families remain unclear. Two recently published studies based on different datasets and age-calibration parameters (DNA type, outgroup selection, fossil calibration points) proposed highly different divergent age estimates and evolutionary scenarios. We obtained nuclear and mitochondrial DNA from mainly museum specimens using target enrichment of single-stranded DNA libraries to generate the first phylogeny of all extant macroscelidean species. We then explored the effects of different parameters (type of DNA, ratio of ingroup to outgroup sampling, number and type of fossil calibration points) and their resulting impacts on age estimates for the origin and initial diversification of Macroscelidea. We show that, even after correcting for substitution saturation, both using mitochondrial DNA in conjunction with nuclear DNA or alone results in much older ages and different branch lengths than when using nuclear DNA alone. We further show that the former effect can be attributed to insufficient amounts of nuclear data. If multiple calibration points are included, the age of the sengi crown group fossil prior has minimal impact on the estimated time frame of sengi evolution. In contrast, the inclusion or exclusion of outgroup fossil priors has a major effect on the resulting node ages. We also find that a reduced sampling of ingroup species does not significantly affect overall age estimates and that terminal specific substitution rates can serve as a means to evaluate the biological likeliness of the produced temporal estimates. Our study demonstrates how commonly varied parameters in temporal calibration of phylogenies affect age estimates. Dated phylogenies should therefore always be seen in the context of the dataset which was used to produce them.

Inter- and intra-island speciation and their morphological and ecological correlates in Aeonium (Crassulaceae), a species-rich Macaronesian radiation.

BACKGROUND AND AIMS: The most species-rich and ecologically diverse plant radiation on the Canary Islands is the Aeonium alliance (Crassulaceae). In island radiations like this, speciation can take place either within islands or following dispersal between islands. Aiming at quantifying intra- and inter-island speciation events in the evolution of Aeonium, and exploring their consequences, we hypothesized that (1) intra-island diversification resulted in stronger ecological divergence of sister lineages, and that (2) taxa on islands with a longer history of habitation by Aeonium show stronger ecological differentiation and produce fewer natural hybrids. METHODS: We studied the biogeographical and ecological setting of diversification processes in Aeonium with a fully sampled and dated phylogeny inferred using a ddRADseq approach. Ancestral areas and biogeographical events were reconstructed in BioGeoBEARS. Eleven morphological characters and three habitat characteristics were taken into account to quantify the morphological and ecological divergence between sister lineages. A co-occurrence matrix of all Aeonium taxa is presented to assess the spatial separation of taxa on each island. KEY RESULTS: We found intra- and inter-island diversification events in almost equal numbers. In lineages that diversified within single islands, morphological and ecological divergence was more pronounced than in lineages derived from inter-island diversification, but only the difference in morphological divergence was significant. Those islands with the longest history of habitation by Aeonium had the lowest percentages of co-occurring and hybridizing taxon pairs compared with islands where Aeonium arrived later. CONCLUSIONS: Our findings illustrate the importance of both inter- and intra-island speciation, the latter of which is potentially sympatric speciation. Speciation on the same island entailed significantly higher levels of morphological divergence compared with inter-island speciation, but ecological divergence was not significantly different. Longer periods of shared island habitation resulted in the evolution of a higher degree of spatial separation and stronger reproductive barriers.

Analysis of variola virus molecular evolution suggests an old origin of the virus consistent with historical records.

Archaeovirology efforts provided a rich portrait of the evolutionary history of variola virus (VARV, the cause of smallpox), which was characterized by lineage extinctions and a relatively recent origin of the virus as a human pathogen (~1700 years ago, ya). This contrasts with historical records suggesting the presence of smallpox as early as 3500 ya. By performing an analysis of ancestry components in modern, historic, and ancient genomes, we unveil the progressive drifting of VARV lineages from a common ancestral population and we show that a small proportion of Viking Age ancestry persisted until the 18th century. After the split of the P-I and P-II lineages, the former experienced a severe bottleneck. With respect to the emergence of VARV as a human pathogen, we revise time estimates by accounting for the time-dependent rate phenomenon. We thus estimate that VARV emerged earlier than 3800 ya, supporting its presence in ancient societies, as pockmarked Egyptian mummies suggest.