StarBEAST2 Brings Faster Species Tree Inference and Accurate Estimates of Substitution Rates

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.

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Traditional multilocus phylogeny fails to fully resolve Palearctic ground squirrels (Spermophilus) relationships but reveals a new species endemic to West Siberia

Previous efforts to reconstruct evolutionary history of Palearctic ground squirrels within the genus Spermophilus have primarily relied on a single mitochondrial marker for phylogenetic data. In this study, we present the first phylogeny with comprehensive taxon sampling of Spermophilus via a conventional multilocus approach utilizing five mitochondrial and five nuclear markers. Through application of the multispecies coalescent model, we constructed a species tree revealing four distinct clades that diverged during the Late Miocene. These clades are 1) S. alaschanicus and S. dauricus from East Asia; 2) S. musicus and S. pygmaeus from East Europe and northwestern Central Asia; 3) the subgenus Colobotis found across Central Asia and its adjacent regions and encompassing S. brevicauda, S. erythrogenys, S. fulvus, S. major, S. pallidicauda, S. ralli, S. relictus, S. selevini, and S. vorontsovi sp. nov.; and 4) a Central/Eastern Europe and Asia Minor clade comprising S. citellus, S. taurensis, S. xanthoprymnus, S. suslicus, and S. odessanus. The latter clade lacked strong support owing to uncertainty of taxonomic placement of S. odessanus and S. suslicus. Resolving relationships within the subgenus Colobotis, which radiated rapidly, remains challenging likely because of incomplete lineage sorting and introgressive hybridization. Most of modern Spermophilus species diversified during the Early-Middle Pleistocene (2.2-1.0 million years ago). We propose a revised taxonomic classification for the genus Spermophilus by recognizing 18 species including a newly identified one (S. vorontsovi sp. nov.), which is found only in a limited area in the southeast of West Siberia. Employing genome-wide single-nucleotide polymorphism genotyping, we substantiated the role of the Ob River as a major barrier ensuring robust isolation of this taxon from S. erythrogenys. Despite its inherent limitations, the traditional multilocus approach remains a valuable tool for resolving relationships and can provide important insights into otherwise poorly understood groups. It is imperative to recognize that additional efforts are needed to definitively determine phylogenetic relationships between certain species of Palearctic ground squirrels.

Branch Length Transforms using Optimal Tree Metric Matching

The abundant discordance between evolutionary relationships across the genome has rekindled interest in ways of comparing and averaging trees on a shared leaf set. However, most attempts at reconciling trees have focused on tree topology, producing metrics for comparing topologies and methods for computing median tree topologies. Using branch lengths, however, has been more elusive, due to several challenges. Species tree branch lengths can be measured in many units, often different from gene trees. Moreover, rates of evolution change across the genome, the species tree, and specific branches of gene trees. These factors compound the stochasticity of coalescence times. Thus, branch lengths are highly heterogeneous across both the genome and the tree. For many downstream applications in phylogenomic analyses, branch lengths are as important as the topology, and yet, existing tools to compare and combine weighted trees are limited. In this paper, we make progress on the question of mapping one tree to another, incorporating both topology and branch length. We define a series of computational problems to formalize finding the best transformation of one tree to another while maintaining its topology and other constraints. We show that all these problems can be solved in quadratic time and memory using a linear algebraic formulation coupled with dynamic programming preprocessing. Our formulations lead to convex optimization problems, with efficient and theoretically optimal solutions. While many applications can be imagined for this framework, we apply it to measure species tree branch lengths in the unit of the expected number of substitutions per site while allowing divergence from ultrametricity across the tree. In these applications, our method matches or surpasses other methods designed directly for solving those problems. Thus, our approach provides a versatile toolkit that finds applications in similar evolutionary questions.

Code availability

The software is available at https://github.com/shayesteh99/TCMM.git .

Data availability

Data are available on Github https://github.com/shayesteh99/TCMM-Data.git .

Multiple transitions between realms shape relict lineages of Proteus cave salamanders

In comparison to biodiversity on Earth's surface, subterranean biodiversity has largely remained concealed. The olm (Proteus anguinus) is one of the most enigmatic extant cave inhabitants, and until now little was known regarding its genetic structure and evolutionary history. Olms inhabit subterranean waters throughout the Dinaric Karst of the western Balkans, with a seemingly uniform phenotypic appearance of cave-specialized traits: an elongate body, snout and limbs, degenerated eyes and loss of pigmentation ("white olm"). Only a single small region in southeastern Slovenia harbours olms with a phenotype typical of surface animals: pigmented skin, eyes, a blunt snout and short limbs ("black olm"). We used a combination of mitochondrial DNA and genome-wide single nucleotide polymorphism data to investigate the molecular diversity, evolutionary history and biogeography of olms along the Dinaric Karst. We found nine deeply divergent species-level lineages that separated between 17 and 4 million years ago, while molecular diversity within lineages was low. We detected no signal of recent admixture between lineages and only limited historical gene flow. Biogeographically, the contemporaneous distribution of lineages mostly mirrors hydrologically separated subterranean environments, while the historical separation of olm lineages follows microtectonic and climatic changes in the area. The reconstructed phylogeny suggests at least four independent transitions to the cave phenotype. Two of the species-level lineages have miniscule ranges and may represent Europe's rarest amphibians. Their rarity and the decline in other lineages call for protection of their subterranean habitats.

Japanese wolves are most closely related to dogs and share DNA with East Eurasian dogs

Although the domestic dog's origin is still unclear, this lineage is believed to have been domesticated from an extinct population of gray wolves, which is expected to be more closely related to dogs than to other populations of gray wolves. Here, we sequence the whole genomes of nine Japanese wolves (7.5-100x: Edo to Meiji periods) and 11 modern Japanese dogs and analyze them together with those from other populations of dogs and wolves. A phylogenomic tree shows that, among the gray wolves, Japanese wolves are closest to the dog, suggesting that the ancestor of dogs is closely related to the ancestor of the Japanese wolf. Based on phylogenetic and geographic relationships, the dog lineage has most likely originated in East Asia, where it diverged from a common ancestor with the Japanese wolf. Since East Eurasian dogs possess Japanese wolf ancestry, we estimate an introgression event from the ancestor of the Japanese wolf to the ancestor of the East Eurasian dog that occurred before the dog's arrival in the Japanese archipelago.

A revision of the trichostrongylid nematode Cooperia Ransom, 1907, from deer game: recent integrative research confirms the existence of the ancient host-specific species Cooperia ventricosa (Rudolphi, 1809)

The trichostrongylid roundworms of the genus Cooperia, which are important in veterinary medicine, currently comprise 19 valid species that parasitize the small intestine of both free-living and domestic ruminants. Only four Cooperia spp. have been reported in Europe, namely C. oncophora, C. punctata, C. curticei and C. pectinata. In 2018-2022, 25 red deer (Cervus elaphus) and 30 sika deer (Cervus nippon) of both sexes and various ages from several remote locations in the Czech Republic were parasitologically examined. Intestinal nematodes of the genus Cooperia were found only in two northern regions. Using the globally recognized key book on trichostrongylid nematodes, they were preliminarily identified as C. pectinata. However, a molecular analysis of cox2 and ITS rDNA gene sequences revealed that Cooperia sp. parasitizing Czech deer is a separate taxon that is more closely related to C. oncophora than to C. pectinata. A subsequent morphological analysis and literature survey confirmed the independence of deer Cooperia sp., which is similar but not identical to bovid C. pectinata. Previous long-term correct identifications of bovid C. pectinata and misidentifications of deer Cooperia species were caused by a fundamental error in the key book mentioned above. Interestingly, the ancient trichostrongylid nematode Strongylus ventricosus from the type host red deer (Cervus elaphus) shot near Greifswald (Germany) was described by Rudolphi in 1809. Rudolphi's type material (one male and four females) was deposited in the Museum für Naturkunde (Berlin). Later, the ancient species S. ventricosus was taken as a synonym for various Cooperia spp. Our current re-examination of the type male indicated that there is a relatively good agreement with our new material from Czech deer regarding the most important characteristics of S. ventricosus (i.e., the shape and size of the male spicules); however, Rudolphi's type material is in rather poor condition. The suggested resurrection of the deer Cooperia sp. in this study as Cooperia ventricosa (Rudolphi, 1809) requires verification by collecting and analyzing new nematode material from the type locality near Greifswald.

Ranked Subtree Prune and Regraft

Phylogenetic trees are a mathematical formalisation of evolutionary histories between organisms, species, genes, cancer cells, etc. For many applications, e.g. when analysing virus transmission trees or cancer evolution, (phylogenetic) time trees are of interest, where branch lengths represent times. Computational methods for reconstructing time trees from (typically molecular) sequence data, for example Bayesian phylogenetic inference using Markov Chain Monte Carlo (MCMC) methods, rely on algorithms that sample the treespace. They employ tree rearrangement operations such as [Formula: see text] (Subtree Prune and Regraft) and [Formula: see text] (Nearest Neighbour Interchange) or, in the case of time tree inference, versions of these that take times of internal nodes into account. While the classic [Formula: see text] tree rearrangement is well-studied, its variants for time trees are less understood, limiting comparative analysis for time tree methods. In this paper we consider a modification of the classical [Formula: see text] rearrangement on the space of ranked phylogenetic trees, which are trees equipped with a ranking of all internal nodes. This modification results in two novel treespaces, which we propose to study. We begin this study by discussing algorithmic properties of these treespaces, focusing on those relating to the complexity of computing distances under the ranked [Formula: see text] operations as well as similarities and differences to known tree rearrangement based treespaces. Surprisingly, we show the counterintuitive result that adding leaves to trees can actually decrease their ranked [Formula: see text] distance, which may have an impact on the results of time tree sampling algorithms given uncertain "rogue taxa".

Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas

Despite their limited vagility and pronounced habitat heterogeneity in the tropics, many anuran species have unexpectedly extensive geographic ranges. One prominent example of this phenomenon is Pithecopus hypochondrialis, which is found in the Cerrado, Guianan savanna, and Llanos domains, as well as isolated tracts of savanna and open habitat within the Amazon Forest. The present study employs an integrative species delimitation approach to test the hypothesis that P. hypochondrialis is in fact a species complex. We also reconstruct the relationships among the lineages delimited here and other Pithecopus species. In this study, we employ Ecological Niche Modelling (ENM) and spatiotemporal phylogeographic reconstruction approaches to evaluate a multitude of scenarios of connectivity across the Neotropical savannas. We identified three divergent lineages, two of which have been described previously. The lineages were allocated to a lowland Pithecopus clade, although the relationships among these lineages are weakly supported. Both the ENM and the phylogeographic reconstruction highlight the occurrence of periods of connectivity among the Neotropical savannas over the course of the Pliocene and Pleistocene epochs. These processes extended from eastern Amazonia to the northern coast of Brazil. The findings of the present study highlight the presence of hidden diversity within P. hypochondrialis, and reinforce the need for a comprehensive taxonomic review. These findings also indicate intricate and highly dynamic patterns of connectivity across the Neotropical savannas that date back to the Pliocene.

Delimiting Species with Single-Locus DNA Sequences

DNA sequences are increasingly used for large-scale biodiversity inventories. Because these genetic data avoid the time-consuming initial sorting of specimens based on their phenotypic attributes, they have been recently incorporated into taxonomic workflows for overlooked and diverse taxa. Major statistical developments have accompanied this new practice, and several models have been proposed to delimit species with single-locus DNA sequences. However, proposed approaches to date make different assumptions regarding taxon lineage history, leading to strong discordance whenever comparisons are made among methods. Distance-based methods, such as Automatic Barcode Gap Discovery (ABGD) and Assemble Species by Automatic Partitioning (ASAP), rely on the detection of a barcode gap (i.e., the lack of overlap in the distributions of intraspecific and interspecific genetic distances) and the associated threshold in genetic distances. Network-based methods, as exemplified by the REfined Single Linkage (RESL) algorithm for the generation of Barcode Index Numbers (BINs), use connectivity statistics to hierarchically cluster-related haplotypes into molecular operational taxonomic units (MOTUs) which serve as species proxies. Tree-based methods, including Poisson Tree Processes (PTP) and the General Mixed Yule Coalescent (GMYC), fit statistical models to phylogenetic trees by maximum likelihood or Bayesian frameworks.Multiple webservers and stand-alone versions of these methods are now available, complicating decision-making regarding the most appropriate approach to use for a given taxon of interest. For instance, tree-based methods require an initial phylogenetic reconstruction, and multiple options are now available for this purpose such as RAxML and BEAST. Across all examined species delimitation methods, judicious parameter setting is paramount, as different model parameterizations can lead to differing conclusions. The objective of this chapter is to guide users step-by-step through all the procedures involved for each of these methods, while aggregating all necessary information required to conduct these analyses. The "Materials" section details how to prepare and format input files, including options to align sequences and conduct tree reconstruction with Maximum Likelihood and Bayesian inference. The Methods section presents the procedure and options available to conduct species delimitation analyses, including distance-, network-, and tree-based models. Finally, limits and future developments are discussed in the Notes section. Most importantly, species delimitation methods discussed herein are categorized based on five indicators: reliability, availability, scalability, understandability, and usability, all of which are fundamental properties needed for any approach to gain unanimous adoption within the DNA barcoding community moving forward.

Evolutionary history and systematics of European blind mole rats (Rodentia: Spalacidae: Nannospalax): Multilocus phylogeny and species delimitation in a puzzling group

Species delimitation is a powerful approach to assist taxonomic decisions in challenging taxa where species boundaries are hard to establish. European taxa of the blind mole rats (genus Nannospalax) display small morphological differences and complex chromosomal evolution at a shallow evolutionary divergence level. Previous analyses led to the recognition of 25 'forms' in their distribution area. We provide a comprehensive framework to improve knowledge on the evolutionary history and revise the taxonomy of European blind mole rats based on samples from all but three of the 25 forms. We sequenced two nuclear-encoded genetic regions and the whole mitochondrial cytochrome b gene for phylogenetic tree reconstructions using concatenation and coalescence-based species-tree estimations. The phylogenetic analyses confirmed that Aegean N. insularis belongs to N. superspecies xanthodon, and that it represents the second known species of this superspecies in Europe. Mainland taxa reached Europe from Asia Minor in two colonisation events corresponding to two superspecies-level taxa: N. superspecies monticola (taxon established herewith) reached Europe c. 2.1 million years ago (Mya) and was followed by N. superspecies leucodon (re-defined herewith) c. 1.5 Mya. Species delimitation allowed the clarification of the taxonomic contents of the above superspecies. N. superspecies monticola contains three species geographically confined to the western periphery of the distribution of blind mole rats, whereas N. superspecies leucodon is more speciose with six species and several additional subspecies. The observed geographic pattern hints at a robust peripatric speciation process and rapid chromosomal evolution. The present treatment is thus regarded as the minimum taxonomic content of each lineage, which can be further refined based on other sources of information such as karyological traits, crossbreeding experiments, etc. The species delimitation models also allowed the recognition of a hitherto unnamed blind mole rat taxon from Albania, described here as a new subspecies.

Phylogenetic review of the comb-tooth blenny genus Hypleurochilus in the northwest Atlantic and Gulf of Mexico

As some of the smallest vertebrates, yet largest producers of consumed reef biomass, cryptobenthic reef fishes serve a disproportionate role in reef ecosystems and are one of the most poorly understood groups of fish. The blenny genera Hypleurochilus and Parablennius are currently considered paraphyletic and the interrelationships of Parablennius have been the focus of recent phylogenetic studies. However, the interrelationships of Hypleurochilus remain understudied. This genus is transatlantically distributed and comprises 11 species with a convoluted taxonomic history. In this study, relationships for ten Hypleurochilus species are resolved using multi-locus nuclear and mtDNA sequence data, morphological data, and mined COI barcode data.  Mitochondrial and nuclear sequence data from 61 individuals collected from the western Atlantic and northern Gulf of Mexico (N. GoM) delimit seven species into a temperate clade, a tropical clade, and a third distinct lineage. This lineage, herein referred to as H. cf. aequipinnis, may represent a species of Hypleurochilus whose range has expanded into the N. GoM. Inclusion of publicly available COI sequence for an additional three species provides further phylogenetic resolution. H. bananensis forms a new eastern Atlantic clade with H. cf. aequipinnis, providing further evidence for a western Atlantic range expansion. Single marker COI delimitation was unable to elucidate the relationships between H. springeri/H. pseudoaequipinnis and between H. multifilis/H. caudovittatus due to incomplete lineage sorting. Mitochondrial data are also unable to accurately resolve the placement of H. bermudensis. However, a comprehensive approach using multi-locus phylogenetic and species delimitation methods was able to resolve these relationships. While mining publicly available sequence data allowed for the inclusion of an increased number of species in the analysis and a more comprehensive phylogeny, it was not without drawbacks, as a handful of sequences are potentially mis-identified. Overall, we find that the recent divergence of some species within this genus and potential introgression events confound the results of single locus delimitation methods, yet a combination of single and multi-locus analyses has allowed for insights into the biogeography of this genus and uncovered a potential transatlantic range expansion.

Genomic evidence sheds new light on phylogeny of Rhabdophis nuchalis (sensu lato) complex (Serpentes: Natricidae)

Stable taxonomy and robust phylogeny are essential for the evolution and conservation of organisms. The Rhabdophis nuchalis (sensu lato) complex presently contains three species (R. nuchalis, R. chiwen, R. pentasupralabialis). Although several studies have explored the diversity and phylogeography of this group, certain issues related to systematics and taxonomy remain unresolved. Here, based on genome-wide data, including single nucleotide polymorphisms (SNPs) generated from ddRAD-seq and mitochondrial DNA (mtDNA), we re-evaluated the phylogenetic relationships and cryptic diversity of this species group. Our results are generally consistent with previous studies but provide some new insights. Phylogenetic relationship reconstruction based on SNPs and mtDNA revealed that three species in the R. nuchalis (sensu lato) complex did not form a monophyly but each species is well supported as monophyletic lineage in SNP-based analyses. Population structure analyses showed genetic admixture between several species pairs. Additionally, the population in eastern Yunnan, China, was identified as a potential cryptic species and thus described as a new species based on morphological data. From our results and previous studies, we redefined the distribution boundary for each species in the R. nuchalis (sensu lato) species complex.

Patterns of phylogenetic diversification in the Dollarbird (Eurystomus orientalis) and Azure Roller (Eurystomus azureus) complex

Genetic isolation and morphological differentiation are two important factors in the speciation process that not always act in concert. A rapid morphological change in a lineage can hide its close relationship to another lineage, while slight morphological differentiation between two taxa can give the appearance of a closer relationship than is actually the case. The Dollarbird (Eurystomus orientalis) and the Azure Roller (Eurystomus azureus) is such an example. Today the Dollarbird and the Azure Roller are unanimously considered to constitute two distinct species, but in a recent genetic study it has been shown that the latter taxon, despite being larger and having a distinctly different coloration, is phylogenetically nested within the former. Its precise placement within this complex has not been determined, however. In this study, we investigate the phylogenetic relationships within the Dollarbird/Azure Roller complex. We estimate divergence times and infer phylogenetic relationships using sequence data from 6,475 genome-wide intronic regions, as well as complete mitochondrial genomes, using both concatenation and multispecies coalescence approaches. We find that within the Dollarbird/Azure Roller complex there are several examples of discrepancies between genetic and morphological differentiation. The Dollarbird is currently divided into between nine to twelve subspecies. Some of these subspecies are poorly differentiated, whereas others are morphologically more clearly discernable. Our data suggest that the complex consist of at least seven distinct genetic lineages that do not entirely match the morphological variation within the group. For instance, our results show that the subspecies solomonensis from the Solomon Islands, despite being morphologically very similar to its geographically closest neighbors, in fact is a highly distinct lineage that became isolated more than 700,000 years ago. In contrast, the morphologically distinct Azure Roller, which is currently treated as a distinct species, is nested within the Dollarbird and forms a slightly younger lineage than solomonensis and is the sister group to a clade with Australian and New Guinean Dollarbirds. Our results also show a deep genetic split within the Dollarbirds on the Asian mainland. This stands in contrast to the apparent clinal morphological variation reported for the birds on the Asian mainland. We also find support for the presence of a genetically distinct clade in the Wallacea region. The birds from the Wallacea region has previously been recognized as a distinct subspecies, connectens, but is currently placed in synonymy of other subspecies. Our results are thus at odds with the current division of the Dollarbird/Azure Roller complex into two species. Given that the species status of azureus is undisputed, the apparent genetic isolation of solomonensis and its clear separation from the other lineages suggests that this taxon also warrants species status. Based on the genetic and morphological variation observed within the Dollarbird/Azure Roller complex there is little doubt that even more taxa should regarded as species, but this require further examination.

QR-STAR: A Polynomial-Time Statistically Consistent Method for Rooting Species Trees Under the Coalescent

We address the problem of rooting an unrooted species tree given a set of unrooted gene trees, under the assumption that gene trees evolve within the model species tree under the multispecies coalescent (MSC) model. Quintet Rooting (QR) is a polynomial time algorithm that was recently proposed for this problem, which is based on the theory developed by Allman, Degnan, and Rhodes that proves the identifiability of rooted 5-taxon trees from unrooted gene trees under the MSC. However, although QR had good accuracy in simulations, its statistical consistency was left as an open problem. We present QR-STAR, a variant of QR with an additional step and a different cost function, and prove that it is statistically consistent under the MSC. Moreover, we derive sample complexity bounds for QR-STAR and show that a particular variant of it based on "short quintets" has polynomial sample complexity. Finally, our simulation study under a variety of model conditions shows that QR-STAR matches or improves on the accuracy of QR. QR-STAR is available in open-source form on github.

Phylogenomic and population genomic analyses reveal the spatial-temporal dynamics of diversification of the Nigella arvensis complex (Ranunculaceae) in the Aegean archipelago

The continental-shelf islands of the Aegean Sea provide an ideal geographical setting for evolutionary-biogeographical studies but disentangling the relationships between palaeogeographical history and the times, orders of modes of taxon divergence is not straightforward. Here, we used phylogenomic and population genomic approaches, based on orthologous gene sequences and transcriptome-derived SNP data, to reconstruct the spatial-temporal evolution of the Aegean Nigella arvensis complex (Ranunculaceae; 11 out of 12 taxa). The group's early diversification in the Early/Mid-Pliocene (c. 3.77 Mya) resulted in three main lineages (Greek mainland vs. central Aegean + Turkish mainland/eastern Aegean islands), while all extant taxa are of Late Plio-/Early Pleistocene origin (c. 3.30-1.59 Mya). Demographic modelling of the outcrossing taxa uncovered disparate modes of (sub)speciation, including divergence with gene flow on the Greek mainland, para- or peripatric diversification across eastern Aegean islands, and a 'mixing-isolation-mixing (MIM)' mode of subspeciation in the Cyclades. The two selfing species (N. stricta, N. doerfleri) evolved independently from the outcrossers. Present-day island configurations are clearly insufficient to explain the spatial-temporal history of lineage diversification and modes of (sub)speciation in Aegean Nigella. Moreover, our identification of positively selected genes in almost all taxa calls into question that this plant group represents a case of 'non-adaptive' radiation. Our study revealed an episodic diversification history of the N. arvensis complex, giving new insight into the modes and drivers of island speciation and adaption across multiple spatiotemporal scales.

Resolving a nearly 90-year-old enigma: The rare Fagus chienii is conspecific with F. hayatae based on molecular and morphological evidence

Taxonomic uncertainties of rare species often hinder effective prioritization for conservation. One such taxonomic uncertainty is the 90-year-old enigma of Fagus chienii. F. chienii was previously only known from the type specimens collected in 1935 in Pingwu County of Sichuan Province, China, and has long been thought to be on the verge of extinction. However, morphological similarities to closely related Fagus species have led many to question the taxonomic status of F. chienii. To clarify this taxonomic uncertainty, we used the newly collected samples to reconstruct a molecular phylogeny of Chinese Fagus species against the phylogenetic backbone of the whole genus using seven nuclear genes. In addition, we examined nine morphological characters to determine whether F. chienii is morphologically distinct from its putatively closest relatives (F. hayatae, F.longipetiolata, and F.lucida). Both morphological and phylogenetic analyses indicated that F. chienii is conspecific with F. hayatae. We recommended that F. chienii should not be treated as a separate species in conservation management. However, conservation strategies such as in situ protection and ex situ germplasm preservation should be adopted to prevent the peculiar "F. chienii" population from extinction.

PhyloAcc-GT: A Bayesian Method for Inferring Patterns of Substitution Rate Shifts on Targeted Lineages Accounting for Gene Tree Discordance

An important goal of evolutionary genomics is to identify genomic regions whose substitution rates differ among lineages. For example, genomic regions experiencing accelerated molecular evolution in some lineages may provide insight into links between genotype and phenotype. Several comparative genomics methods have been developed to identify genomic accelerations between species, including a Bayesian method called PhyloAcc, which models shifts in substitution rate in multiple target lineages on a phylogeny. However, few methods consider the possibility of discordance between the trees of individual loci and the species tree due to incomplete lineage sorting, which might cause false positives. Here, we present PhyloAcc-GT, which extends PhyloAcc by modeling gene tree heterogeneity. Given a species tree, we adopt the multispecies coalescent model as the prior distribution of gene trees, use Markov chain Monte Carlo (MCMC) for inference, and design novel MCMC moves to sample gene trees efficiently. Through extensive simulations, we show that PhyloAcc-GT outperforms PhyloAcc and other methods in identifying target lineage-specific accelerations and detecting complex patterns of rate shifts, and is robust to specification of population size parameters. PhyloAcc-GT is usually more conservative than PhyloAcc in calling convergent rate shifts because it identifies more accelerations on ancestral than on terminal branches. We apply PhyloAcc-GT to two examples of convergent evolution: flightlessness in ratites and marine mammal adaptations, and show that PhyloAcc-GT is a robust tool to identify shifts in substitution rate associated with specific target lineages while accounting for incomplete lineage sorting.

Their fates intertwined: diversification patterns of the Asian gliding vertebrates may have been forged by dipterocarp trees

The repeated evolution of gliding in diverse Asian vertebrate lineages is hypothesized to have been triggered by the dominance of tall dipterocarp trees in the tropical forests of Southeast Asia. These dipterocarp forests have acted as both centres of diversification and climatic refugia for gliding vertebrates, and support most of their extant diversity. We predict similarities in the diversification patterns of dipterocarp trees and gliding vertebrates, and specifically test whether episodic diversification events such as rate shifts and/or mass extinctions were temporally congruent in these groups. We analysed diversification patterns in reconstructed timetrees of Asian dipterocarps, the most speciose gliding vertebrates from different classes (Draco lizards, gliding frogs and Pteromyini squirrels) and compared them with similar-sized clades of non-gliding relatives (Diploderma lizards, Philautus frogs and Callosciurinae squirrels) from Southeast Asia. We found significant declines in net-diversification rates of dipterocarps and the gliding vertebrates during the Pliocene-Pleistocene, but not in the non-gliding groups. We conclude that the homogeneity and temporal coincidence of these rate declines point to a viable ecological correlation between dipterocarps and the gliding vertebrates. Further, we suggest that while the diversification decay in dipterocarps was precipitated by post-Miocene aridification of Asia, the crises in the gliding vertebrates were induced by both events concomitantly.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Estimation of species divergence times in presence of cross-species gene flow

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent].

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Weighted ASTRID: fast and accurate species trees from weighted internode distances

BACKGROUND

Species tree estimation is a basic step in many biological research projects, but is complicated by the fact that gene trees can differ from the species tree due to processes such as incomplete lineage sorting (ILS), gene duplication and loss (GDL), and horizontal gene transfer (HGT), which can cause different regions within the genome to have different evolutionary histories (i.e., "gene tree heterogeneity"). One approach to estimating species trees in the presence of gene tree heterogeneity resulting from ILS operates by computing trees on each genomic region (i.e., computing "gene trees") and then using these gene trees to define a matrix of average internode distances, where the internode distance in a tree T between two species x and y is the number of nodes in T between the leaves corresponding to x and y. Given such a matrix, a tree can then be computed using methods such as neighbor joining. Methods such as ASTRID and NJst (which use this basic approach) are provably statistically consistent, very fast (low degree polynomial time) and have had high accuracy under many conditions that makes them competitive with other popular species tree estimation methods. In this study, inspired by the very recent work of weighted ASTRAL, we present weighted ASTRID, a variant of ASTRID that takes the branch uncertainty on the gene trees into account in the internode distance.

RESULTS

Our experimental study evaluating weighted ASTRID typically shows improvements in accuracy compared to the original (unweighted) ASTRID, and shows competitive accuracy against weighted ASTRAL, the state of the art. Our re-implementation of ASTRID also improves the runtime, with marked improvements on large datasets.

CONCLUSIONS

Weighted ASTRID is a new and very fast method for species tree estimation that typically improves upon ASTRID and has comparable accuracy to weighted ASTRAL, while remaining much faster. Weighted ASTRID is available at https://github.com/RuneBlaze/internode .

Phylogenomic branch length estimation using quartets

MOTIVATION

Branch lengths and topology of a species tree are essential in most downstream analyses, including estimation of diversification dates, characterization of selection, understanding adaptation, and comparative genomics. Modern phylogenomic analyses often use methods that account for the heterogeneity of evolutionary histories across the genome due to processes such as incomplete lineage sorting. However, these methods typically do not generate branch lengths in units that are usable by downstream applications, forcing phylogenomic analyses to resort to alternative shortcuts such as estimating branch lengths by concatenating gene alignments into a supermatrix. Yet, concatenation and other available approaches for estimating branch lengths fail to address heterogeneity across the genome.

RESULTS

In this article, we derive expected values of gene tree branch lengths in substitution units under an extension of the multispecies coalescent (MSC) model that allows substitutions with varying rates across the species tree. We present CASTLES, a new technique for estimating branch lengths on the species tree from estimated gene trees that uses these expected values, and our study shows that CASTLES improves on the most accurate prior methods with respect to both speed and accuracy.

AVAILABILITY AND IMPLEMENTATION

CASTLES is available at https://github.com/ytabatabaee/CASTLES.

Global phylogenomic assessment of Leptoseris and Agaricia reveals substantial undescribed diversity at mesophotic depths

BACKGROUND

Mesophotic coral communities are increasingly gaining attention for the unique biological diversity they host, exemplified by the numerous mesophotic fish species that continue to be discovered. In contrast, many of the photosynthetic scleractinian corals observed at mesophotic depths are assumed to be depth-generalists, with very few species characterised as mesophotic-specialists. This presumed lack of a specialised community remains largely untested, as phylogenetic studies on corals have rarely included mesophotic samples and have long suffered from resolution issues associated with traditional sequence markers.

RESULTS

Here, we used reduced-representation genome sequencing to conduct a phylogenomic assessment of the two dominant mesophotic genera of plating corals in the Indo-Pacific and Western Atlantic, respectively, Leptoseris and Agaricia. While these genome-wide phylogenies broadly corroborated the morphological taxonomy, they also exposed deep divergences within the two genera and undescribed diversity across the current taxonomic species. Five of the eight focal species consisted of at least two sympatric and genetically distinct lineages, which were consistently detected across different methods.

CONCLUSIONS

The repeated observation of genetically divergent lineages associated with mesophotic depths highlights that there may be many more mesophotic-specialist coral species than currently acknowledged and that an urgent assessment of this largely unstudied biological diversity is warranted.

Power of Bayesian and Heuristic Tests to Detect Cross-Species Introgression with Reference to Gene Flow in the Tamias quadrivittatus Group of North American Chipmunks

In the past two decades, genomic data have been widely used to detect historical gene flow between species in a variety of plants and animals. The Tamias quadrivittatus group of North America chipmunks, which originated through a series of rapid speciation events, are known to undergo massive amounts of mitochondrial introgression. Yet in a recent analysis of targeted nuclear loci from the group, no evidence for cross-species introgression was detected, indicating widespread cytonuclear discordance. The study used the heuristic method HYDE to detect gene flow, which may suffer from low power. Here we use the Bayesian method implemented in the program BPP to re-analyze these data. We develop a Bayesian test of introgression, calculating the Bayes factor via the Savage-Dickey density ratio using the Markov chain Monte Carlo (MCMC) sample under the model of introgression. We take a stepwise approach to constructing an introgression model by adding introgression events onto a well-supported binary species tree. The analysis detected robust evidence for multiple ancient introgression events affecting the nuclear genome, with introgression probabilities reaching 63%. We estimate population parameters and highlight the fact that species divergence times may be seriously underestimated if ancient cross-species gene flow is ignored in the analysis. We examine the assumptions and performance of HYDE and demonstrate that it lacks power if gene flow occurs between sister lineages or if the mode of gene flow does not match the assumed hybrid-speciation model with symmetrical population sizes. Our analyses highlight the power of likelihood-based inference of cross-species gene flow using genomic sequence data. [Bayesian test; BPP; chipmunks; introgression; MSci; multispecies coalescent; Savage-Dickey density ratio.].

Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance

Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture

BACKGROUND

Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.

RESULTS

Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia.

CONCLUSIONS

We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

A New Assessment of Robust Capuchin Monkey (Sapajus) Evolutionary History Using Genome-Wide SNP Marker Data and a Bayesian Approach to Species Delimitation

Robust capuchin monkeys, Sapajus genus, are among the most phenotypically diverse and widespread groups of primates in South America, with one of the most confusing and often shifting taxonomies. We used a ddRADseq approach to generate genome-wide SNP markers for 171 individuals from all putative extant species of Sapajus to access their evolutionary history. Using maximum likelihood, multispecies coalescent phylogenetic inference, and a Bayes Factor method to test for alternative hypotheses of species delimitation, we inferred the phylogenetic history of the Sapajus radiation, evaluating the number of discrete species supported. Our results support the recognition of three species from the Atlantic Forest south of the São Francisco River, with these species being the first splits in the robust capuchin radiation. Our results were congruent in recovering the Pantanal and Amazonian Sapajus as structured into three monophyletic clades, though new morphological assessments are necessary, as the Amazonian clades do not agree with previous morphology-based taxonomic distributions. Phylogenetic reconstructions for Sapajus occurring in the Cerrado, Caatinga, and northeastern Atlantic Forest were less congruent with morphology-based phylogenetic reconstructions, as the bearded capuchin was recovered as a paraphyletic clade, with samples from the Caatinga biome being either a monophyletic clade or nested with the blond capuchin monkey.

A taxonomic revision of Boiga multomaculata (Boie, 1827) and B. ochracea (Theobald, 1868), with the description of a new subspecies (Squamata, Serpentes, Colubridae)

The analyses of molecular genetic data (mtDNA markers 16S, ND4, CYTB, and the nuclear marker c-mos) provided evidence that the Asian cat snake taxa Boiga multomaculata and B. ochracea actually represent a single species. They form mixed clades of low intraclade genetic differentiation. This evidence for conspecificy is supported by the lack of differentiation in all examined pholidotic and morphometric characters. Therefore, we formally place Dipsas ochracea Theobald, 1868 in the synonymy of Dipsas multomaculata Boie, 1827. We provide a summary of the tangled taxonomic history of the taxa involved in this study. Also, we resurrect Dipsadomorphus stoliczkae Wall, 1909 from synonymy of B. ochracea, for specimens exhibiting 21 midbody dorsal scale rows. Boiga stoliczkae is found in the Himalayas north and west of the Brahmaputra valley. Finally, based on the detection of historical genetic lineages within the newly defined species Boiga multomaculata we recognize three subspecies: Boiga multomaculata multomaculata (Boie, 1827), Boiga multomaculata ochracea (Theobald, 1868), and Boiga multomaculata septentrionalis n. ssp. which is distributed in northern Myanmar and Assam and Nagaland, India. We designate BMNH 1946.1.2.60 (1) as neotype of Dipsas ochracea Theobald, (2) as lectotype of D. ochraceus Günther, and (3) as lectotype of Boiga ochracea walli Smith, thereby making these names objective synonyms. Finally, we designate BMNH 94.12.31.55 as lectotype of Dipsadomorphus stoliczkae Wall.

The use of Anchored Hybrid Enrichment data to resolve higher-level phylogenetic relationships: A proof-of-concept applied to Asterales (Eudicotyledoneae; Angiosperms)

Anchored Hybrid Enrichment (AHE) is a tool for capturing orthologous regions of the nuclear genome shared in low or single copy across lineages. Despite the increasing number of studies using this method, its usefulness to estimate relationships at deeper taxonomic levels in plants has not been fully explored. Here we present a proof of concept about the performance of nuclear loci obtained with AHE to infer phylogenetic relationships and explore the use of gene sampling schemes to estimate divergence times in Asterales. We recovered low-copy nuclear loci using the AHE method from herbarium material and silica-preserved samples. Maximum likelihood, Bayesian inference, and coalescence approaches were used to reconstruct phylogenomic relationships. Dating analyses were conducted under a multispecies coalescent approach by jointly inferring species tree and divergence times with random gene sampling schemes and multiple calibrations. We recovered 403 low-copy nuclear loci for 63 species representing nine out of eleven families of Asterales. Phylogenetic hypotheses were congruent among the applied methods and previously published results. Analyses with concatenated datasets were strongly supported, but coalescence-based analyses showed low support for the phylogenetic position of families Argophyllaceae and Alseuosmiaceae. Estimated family ages were congruent among gene sampling schemes, with the mean age for Asterales around 130 Myr. Our study documents the usefulness of AHE for resolving phylogenetic relationships at deep phylogenetic levels in Asterales. Observed phylogenetic inconsistencies were possibly due to the non-inclusion of families Phellinceae and Pentaphragmataceae. Random gene sampling schemes produced consistent age estimates with coalescence and species tree relaxed clock approaches.

Comparing inference under the multispecies coalescent with and without recombination

Accurate inference of population parameters plays a pivotal role in unravelling evolutionary histories. While recombination has been universally accepted as a fundamental process in the evolution of sexually reproducing organisms, it remains challenging to model it exactly. Thus, existing coalescent-based approaches make different assumptions or approximations to facilitate phylogenetic inference, which can potentially bring about biases in estimates of evolutionary parameters when recombination is present. In this article, we evaluate the performance of population parameter estimation using three methods-StarBEAST2, SNAPP, and diCal2-that represent three different types of inference. We performed whole-genome simulations in which recombination rates, mutation rates, and levels of incomplete lineage sorting were varied. We show that StarBEAST2 using short or medium-sized loci is robust to realistic rates of recombination, which is in agreement with previous studies. SNAPP, as expected, is generally unaffected by recombination events. Most surprisingly, diCal2, a method that is designed to explicitly account for recombination, performs considerably worse than other methods under comparison.

Late Cenozoic History And The Role Of Beringia In Assembling A Holarctic Cestode Species Complex

The dynamic climate history that drove sea level fluctuation during past glacial periods mediated the movement of organisms between Asia and North America via the Bering Land Bridge. Investigations of the biogeographic histories of small mammals and their parasites demonstrate facets of a complex history of episodic geographic colonization and refugial isolation that structured diversity across the Holarctic. We use a large multi-locus nuclear DNA sequence dataset to robustly resolve relationships within the cestode genus Arostrilepis (Cyclophyllidea: Hymenolepididae), a widespread parasite of predominantly arvicoline rodents (voles, lemmings). Using this phylogeny, we confirm that several Asian Arostrilepis lineages colonized North America during up to four distinct glacial periods in association with different rodent hosts, consistent with taxon-pulse dynamics. A previously inferred westward dispersal across the land bridge is rejected. We also refine interpretations of past host colonization, providing evidence for several distinct episodes of expanding host range, which probably contributed to diversification by Arostrilepis. Finally, Arostrilepis is shown to be paraphyletic with respect to Hymenandrya thomomyis, a parasite of pocket gophers, confirming that ancient Arostrilepis species colonized new host lineages upon arriving in North America.

Evaluating species boundaries using coalescent delimitation in pine-killing Monochamus (Coleoptera: Cerambycidae) sawyer beetles

Plant-feeding beetle species are diverse and often individually highly variable. Accurate classifications can be difficult to establish yet are essential for study of evolutionary patterns and processes. Molecular data are key to further characterizing morphologically difficult groups and defining genus and species boundaries. Monochamus Dejean species are ecologically and economically significant, and in coniferous forests they vector the nematode that causes Pine Wilt Disease. This study uses nuclear and mitochondrial genes to test the monophyly and relationships of Monochamus and applies coalescent methods to further delimit the conifer-feeding species. Monochamus has also included approximately 120 Old World species associated with diverse angiosperm tree species. We sample from these additional morphologically diverse species to determine their placement in the Lamiini. Through supermatrix and coalescent methods, the higher-level relationships of Monochamus show that conifer-feeders are a monophyletic group that includes the type species and has split into Nearctic and Palearctic clades. Molecular dating indicates a single dispersal of conifer-feeders to North America over the second Bering Land Bridge circa 5.3 Ma. All other Monochamus sampled fall in different parts of the Lamiini tree. Small-bodied angiosperm-feeding Monochamus group with the monotypic genus Microgoes Casey. The African Monochamus subgenera sampled are distantly related to the conifer-feeding clade. The multispecies coalescent delimitation methods BPP and STACEY delimit 17 conifer-feeding Monochamus species for a total of 18 species, and supports the retention of all current species. An interrogation with nuclear gene allele phasing reveals that unphased data can be unreliable for accurate delimitations and divergence times. The delimited species are discussed with integrative evidence, highlighting real-world challenges in recognizing the completion of speciation trajectories.

Evolutionary origins of the eastern North American-Mesoamerican floristic disjunction: Current status and future prospects

Biogeographic disjunction patterns, where multiple taxa are shared between isolated geographic areas, represent excellent systems for investigating the historical assembly of modern biotas and fundamental biological processes such as speciation, diversification, niche evolution, and evolutionary responses to climate change. Studies on plant genera disjunct across the northern hemisphere, particularly between eastern North America (ENA) and eastern Asia (EAS), have yielded tremendous insight on the geologic history and assembly of rich temperate floras. However, one of the most prevalent disjunction patterns involving ENA forests has been largely overlooked: that of taxa disjunct between ENA and cloud forests of Mesoamerica (MAM), with examples including Acer saccharum, Liquidambar styraciflua, Cercis canadensis, Fagus grandifolia, and Epifagus virginiana. Despite the remarkable nature of this disjunction pattern, which has been recognized for over 75 years, there have been few recent efforts to empirically examine its evolutionary and ecological origins. Here I synthesize previous systematic, paleobotanical, phylogenetic, and phylogeographic studies to establish what is known about this disjunction pattern to provide a roadmap for future research. I argue that this disjunction pattern, and the evolution and fossil record of the Mexican flora more broadly, represents a key missing piece in the broader puzzle of northern hemisphere biogeography. I also suggest that the ENA-MAM disjunction represents an excellent system for examining fundamental questions about how traits and life history strategies mediate plant evolutionary responses to climate change and for predicting how broadleaf temperate forests will respond to the ongoing climatic pressures of the Anthropocene.

Historical demography and climatic niches of the Natal multimammate mouse (Mastomys natalensis) in the Zambezian region

The Natal multimammate mouse (Mastomys natalensis) is the most widespread rodent species in sub-Saharan Africa, often studied as an agricultural pest and reservoir of viruses. Its mitochondrial (Mt) phylogeny revealed six major lineages parapatrically distributed across open habitats of sub-Saharan Africa. In this study we used 1949 sequences of the mitochondrial cytochrome b gene to elaborate on distribution and evolutionary history of three Mt lineages inhabiting the open habitats of the Zambezian region (corresponding roughly to the African savannas south of the Equator). We describe in more detail contact zones between the lineages—their location and extent of co-occurrence within localities—and infer past population trends. The estimates are interpreted in the light of climatic niche models. The lineages underwent reduction in effective population size during the last glacial, but they spread widely after that: two of them after the last glacial maximum and the last one in mid-Holocene. The centers of expansion, i.e., possible long-term savanna refugia, were estimated to lie close to the Eastern Arc Mountains and lakes of the Great African Rift, geomorphological structures likely to have had long-term influence on geographical distribution of the lineages. Environmental niche modeling shows climate could also affect the broad scale distribution of the lineages but is unlikely to explain the narrow width of the contact zones. The intraspecific Mt differentiation of M. natalensis echoes phylogeographic patterns observed in multiple co-distributed mammal species, which suggests the mammal communities in the region are shaped by the same long-term processes.

DISCO+QR: rooting species trees in the presence of GDL and ILS

Motivation

Genes evolve under processes such as gene duplication and loss (GDL), so that gene family trees are multi-copy, as well as incomplete lineage sorting (ILS); both processes produce gene trees that differ from the species tree. The estimation of species trees from sets of gene family trees is challenging, and the estimation of rooted species trees presents additional analytical challenges. Two of the methods developed for this problem are STRIDE, which roots species trees by considering GDL events, and Quintet Rooting (QR), which roots species trees by considering ILS.

Results

We present DISCO+QR, a new approach to rooting species trees that first uses DISCO to address GDL and then uses QR to perform rooting in the presence of ILS. DISCO+QR operates by taking the input gene family trees and decomposing them into single-copy trees using DISCO and then roots the given species tree using the information in the single-copy gene trees using QR. We show that the relative accuracy of STRIDE and DISCO+QR depend on the properties of the dataset (number of species, genes, rate of gene duplication, degree of ILS and gene tree estimation error), and that each provides advantages over the other under some conditions.

Availability and implementation

DISCO and QR are available in github.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Initiation of speciation across multiple dimensions in a rock-restricted, tropical lizard

Population isolation and concomitant genetic divergence, resulting in strong phylogeographical structure, is a core aspect of speciation initiation. If and how speciation then proceeds and ultimately completes depends on multiple factors that mediate reproductive isolation, including divergence in genomes, ecology and mating traits. Here we explored these multiple dimensions in two young (Plio-Pleistocene) species complexes of gekkonid lizards (Heteronotia) from the Kimberley-Victoria River regions of tropical Australia. Using mitochondrial DNA screening and exon capture phylogenomics, we show that the rock-restricted Heteronotia planiceps exhibits exceptional fine-scale phylogeographical structure compared to the codistributed habitat generalist Heteronotia binoei. This indicates pervasive population isolation and persistence in the rock-specialist, and thus a high rate of speciation initiation across this geographically complex region, with levels of genomic divergence spanning the "grey zone" of speciation. Proximal lineages of H. planiceps were often separated by different rock substrates, suggesting a potential role for ecological isolation; however, phylogenetic incongruence and historical introgression were inferred between one such pair. Ecomorphological divergence among lineages within both H. planiceps and H. binoei was limited, except that limestone-restricted lineages of H. planiceps tended to be larger than rock-generalists. By contrast, among-lineage divergence in the chemical composition of epidermal pore secretions (putative mating trait) exceeded ecomorphology in both complexes, but with less trait overlap among lineages in H. planiceps. This system-particularly the rock-specialist H. planiceps-highlights the role of multidimensional divergence during incipient speciation, with divergence in genomes, ecomorphology and chemical signals all at play at very fine spatial scales.

On the effects of selection and mutation on species tree inference

The effect of selection acting on regions of the genome on the accuracy of species-level phylogenetic inference using methods that do not explicitly model selection is an open question that is relevant to most, if not all, phylogenomic studies. To address this, we derive a mathematical approximation to the Wright-Fisher model with mutation and selection in the limit as the population size becomes large. In contrast to previous approximations based on diffusion processes, our approximation can be used to study the distribution of coalescent times for an arbitrary number of lineages, allowing calculation of the probability distribution of gene genealogies under the coalescent model. We use these calculations to show that direct selection at strengths typically encountered in practice has only a small effect on the distribution of coalescent times, and hence on the distribution of gene trees. This implies that many coalescent-based methods for estimating the species tree topology will be robust to the presence of selection in a subset of the underlying genes. Selection will, however, bias the estimation of speciation times, causing them to underestimate the true speciation times. Our model captures the effects of selection on the genealogies that generate the observed sequence data, but does not model selective pressures that act only on the subsequent sequences or that negatively impact gene tree estimation.

Phylogeny of Saxifraga section Saxifraga subsection Arachnoideae (Saxifragaceae) and the origin of low elevation shade-dwelling species

Saxifraga section Saxifraga subsection Arachnoideae is a lineage of 12 species distributed mainly in the European Alps. It is unusual in terms of ecological diversification by containing both high elevation species from exposed alpine habitats and low elevation species from shady habitats such as overhanging rocks and cave entrances. Our aims are to explore which of these habitat types is ancestral, and to identify the possible drivers of this remarkable ecological diversification. Using a Hybseq DNA-sequencing approach and a complete species sample we reconstructed and dated the phylogeny of subsection Arachnoideae. Using Landolt indicator values, this phylogenetic tree was used for the reconstruction of the evolution of temperature, light and soil pH requirements in this lineage. Diversification of subsection Arachnoideae started in the late Pliocene and continued through the Pleistocene. Both diversification among and within clades was largely allopatric, and species from shady habitats with low light requirements are distributed in well-known refugia. We hypothesize that low light requirements evolved when species persisting in cold-stage refugia were forced into marginal habitats by more competitive warm-stage vegetation. While we do not claim that such competition resulted in speciation, it very likely resulted in adaptive evolution.

Contrasting continental patterns of adaptive population divergence in the holarctic ectomycorrhizal fungus Boletus edulis

In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome-wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis, a globally distributed species complex with broad ecological amplitude. By whole-genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphisms to characterize patterns of genome-wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62-2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene-Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome-wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation.

Considerations on the Systematics of Neotropical Pierina, with the Description of Two New Species of Phulia Herrich-Schäffer from the Peruvian Andes (Lepidoptera: Pieridae, Pierinae, Pierini)

A comparative analysis of high-Andean Pierina was carried out, including a total of 25 species. Based on morphological evidence, with an emphasis on venation and genitalia and molecular data, using three genetic markers, we confirm the recent subjective synonymy of the generic names Tatochila Butler, 1870, Piercolias, Staudinger, 1894, Hypsochila Ureta, 1955, Infraphulia Field, 1958, Pierphulia Field, 1958, and Theochila Field, 1958 with Phulia Herrich-Schäffer, 1867. Two new species are described, namely Phulia stoddardi Pyrcz & Cerdeña n. sp., from the Andes of Central Peru, which occurs at an unusually high altitude of close to 5000 m a.s.l. in dry puna habitat, and Phulia phantasma Lamas, Willmott & Boyer n. sp., from dry montane forests in northern Peru and southern Ecuador. An overview of high-elevation butterflies is presented, with some discussion on adaptations to this environment.

Insights from the genomes of 4 diploid Camelina spp

Plant evolution has been a complex process involving hybridization and polyploidization making understanding the origin and evolution of a plant's genome challenging even once a published genome is available. The oilseed crop, Camelina sativa (Brassicaceae), has a fully sequenced allohexaploid genome with 3 unknown ancestors. To better understand which extant species best represent the ancestral genomes that contributed to C. sativa's formation, we sequenced and assembled chromosome level draft genomes for 4 diploid members of Camelina: C. neglecta C. hispida var. hispida, C. hispida var. grandiflora, and C. laxa using long and short read data scaffolded with proximity data. We then conducted phylogenetic analyses on regions of synteny and on genes described for Arabidopsis thaliana, from across each nuclear genome and the chloroplasts to examine evolutionary relationships within Camelina and Camelineae. We conclude that C. neglecta is closely related to C. sativa's sub-genome 1 and that C. hispida var. hispida and C. hispida var. grandiflora are most closely related to C. sativa's sub-genome 3. Further, the abundance and density of transposable elements, specifically Helitrons, suggest that the progenitor genome that contributed C. sativa's sub-genome 3 maybe more similar to the genome of C. hispida var. hispida than that of C. hispida var. grandiflora. These diploid genomes show few structural differences when compared to C. sativa's genome indicating little change to chromosome structure following allopolyploidization. This work also indicates that C. neglecta and C. hispida are important resources for understanding the genetics of C. sativa and potential resources for crop improvement.

Cryptic diversity impacts model selection and macroevolutionary inferences in diversification analyses

Species persist in landscapes through ecological dynamics but proliferate at wider spatial scales through evolutionary mechanisms. Disentangling the contribution of each dynamic is challenging, but the increasing use of dated molecular phylogenies opened new perspectives. First, the increasing use of DNA sequences in biodiversity inventory shed light on a substantial amount of cryptic diversity in species-rich ecosystems. Second, explicit diversification models accounting for various eco-evolutionary models are now available. Integrating both advances, we explored diversification trajectories among 10 lineages of freshwater fishes in Sundaland, for which time-calibrated and taxonomically rich phylogenies are available. By fitting diversification models to dated phylogenies and incorporating DNA-based species delimitation methods, the impact of cryptic diversity on diversification model selection and related inferences is explored. Eight clades display constant speciation rate model as the most likely if cryptic diversity is accounted, but nine display a signature of diversification slowdowns when cryptic diversity is ignored. Cryptic diversification occurs during the last 5 Myr for most groups, and palaeoecological models received little support. Most cryptic lineages display restricted range distribution, supporting geographical isolation across homogeneous landscapes as the main driver of diversification. These patterns question the persistence of cryptic diversity and its role during species proliferation.

Estimation of speciation times under the multispecies coalescent

MOTIVATION

The multispecies coalescent model is now widely accepted as an effective model for incorporating variation in the evolutionary histories of individual genes into methods for phylogenetic inference from genome-scale data. However, because model-based analysis under the coalescent can be computationally expensive for large datasets, a variety of inferential frameworks and corresponding algorithms have been proposed for estimation of species-level phylogenies and associated parameters, including speciation times and effective population sizes.

RESULTS

We consider the problem of estimating the timing of speciation events along a phylogeny in a coalescent framework. We propose a maximum a posteriori estimator based on composite likelihood (MAPCL) for inferring these speciation times under a model of DNA sequence evolution for which exact site-pattern probabilities can be computed under the assumption of a constant θ throughout the species tree. We demonstrate that the MAPCL estimates are statistically consistent and asymptotically normally distributed, and we show how this result can be used to estimate their asymptotic variance. We also provide a more computationally efficient estimator of the asymptotic variance based on the non-parametric bootstrap. We evaluate the performance of our method using simulation and by application to an empirical dataset for gibbons.

AVAILABILITY AND IMPLEMENTATION

The method has been implemented in the PAUP* program, freely available at https://paup.phylosolutions.com for Macintosh, Windows and Linux operating systems.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent progress on methods for estimating and updating large phylogenies

With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Weighting by Gene Tree Uncertainty Improves Accuracy of Quartet-based Species Trees

Phylogenomic analyses routinely estimate species trees using methods that account for gene tree discordance. However, the most scalable species tree inference methods, which summarize independently inferred gene trees to obtain a species tree, are sensitive to hard-to-avoid errors introduced in the gene tree estimation step. This dilemma has created much debate on the merits of concatenation versus summary methods and practical obstacles to using summary methods more widely and to the exclusion of concatenation. The most successful attempt at making summary methods resilient to noisy gene trees has been contracting low support branches from the gene trees. Unfortunately, this approach requires arbitrary thresholds and poses new challenges. Here, we introduce threshold-free weighting schemes for the quartet-based species tree inference, the metric used in the popular method ASTRAL. By reducing the impact of quartets with low support or long terminal branches (or both), weighting provides stronger theoretical guarantees and better empirical performance than the unweighted ASTRAL. Our simulations show that weighting improves accuracy across many conditions and reduces the gap with concatenation in conditions with low gene tree discordance and high noise. On empirical data, weighting improves congruence with concatenation and increases support. Together, our results show that weighting, enabled by a new optimization algorithm we introduce, improves the utility of summary methods and can reduce the incongruence often observed across analytical pipelines.