Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies

Highly diversified mitochondrial genomes provide new evidence for interordinal relationships in the Arachnida

Arachnida is an exceptionally diverse class in the Arthropoda, consisting of 20 orders and playing crucial roles in the terrestrial ecosystems. However, their interordinal relationships have been debated for over a century. Rearranged or highly rearranged mitochondrial genomes (mitogenomes) were consistently found in this class, but their various extent in different lineages and efficiency for resolving arachnid phylogenies are unclear. Here, we reconstructed phylogenetic trees using mitogenome sequences of 290 arachnid species to decipher interordinal relationships as well as diversification through time. Our results recovered monophyly of ten orders (i.e. Amblypygi, Araneae, Ixodida, Mesostigmata, Opiliones, Pseudoscorpiones, Ricinulei, Sarcoptiformes, Scorpiones and Solifugae), while rejecting monophyly of the Trombidiformes due to the unstable position of the Eriophyoidea. The monophyly of Acari (subclass) was rejected, possibly due to the long-branch attraction of the Pseudoscorpiones. The monophyly of Arachnida was further rejected because the Xiphosura nested within arachnid orders with unstable positions. Mitogenomes that are highly rearranged in mites but less rearranged or conserved in the remaining lineages point to their exceptional diversification in mite orders; however, shared derived mitochondrial (mt) gene clusters were found within superfamilies rather than interorders, confusing phylogenetic signals in arachnid interordinal relationships. Molecular dating results show that arachnid orders have ancient origins, ranging from the Ordovician to the Carboniferous, yet have significantly diversified since the Cretaceous in orders Araneae, Mesostigmata, Sarcoptiformes, and Trombidiformes. By summarizing previously resolved key positions of some orders, we propose a plausible arachnid tree of life. Our results underline a more precise framework for interordinal phylogeny in the Arachnida and provide new insights into their ancient evolution.

Excluding loci with substitution saturation improves inferences from phylogenomic data

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.]

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

Multilocus phylogeny and cryptic diversity of white-toothed shrews (Mammalia, Eulipotyphla, Crocidura) in China

BACKGROUND

Crocidura, the most speciose mammalian genus, occurs across much of Asia, Europe and Africa. The taxonomy of Chinese representatives has been studied primarily based on cursory morphological comparisons and their molecular phylogenetic analyses remain unexplored. In order to understand the phylogeny of this group in China, we estimated the first multilocus phylogeny and conducted species delimitation, including taxon sampling throughout their distribution range.

RESULTS

We obtained one mitochondrial gene (cytb) (~ 1, 134 bp) and three nuclear genes (ApoB, BRCA1, RAG1) (~ 2, 170 bp) for 132 samples from 57 localities. Molecular analyses identified at least 14 putative species that occur within two major well-supported groups in China. Polyphyletic C. wuchihensis appears to be composed of two putative species. Two subspecies, C. rapax rapax and C. rapax kurodai should be elevated to full species status. A phylogenetic tree based on mitochondrial gene from Asian Crocidura species showed that the C. rapax rapax is embedded within C. attenuata, making the latter a paraphyletic group. Three strongly supported undescribed species (C. sp.1, C. sp.2 and C. sp.3) are revealed from Zada County of Tibet (Western China), Hongjiang County of Hunan Province (Central China) and Dongyang County of Zhejiang Province (Eastern China), Motuo County of Tibet, respectively. The divergence time estimation suggested that China's Crocidura species began to diversify during the late Pliocene (3.66 Ma) and the Early Pleistocene (2.29 Ma), followed by a series of diversifications through the Pleistocene.

CONCLUSIONS

The cryptic diversity found in this study indicated that the number of species is strongly underestimated under the current taxonomy. We propose that the three undescribed species should be evaluated using extensive taxon sampling and comprehensive morphological and morphometric approaches. Climate change since the late Pliocene and the uplift of the Qinghai-Tibet Plateau may result in the diversification and speciation of China's Crocidura species. In short, the underestimated diversity underlines the need for a taxonomic revision of Chinese Crocidura species.

The anatomy of an unstable node: a Levantine relict precipitates phylogenomic dissolution of higher-level relationships of the armoured harvestmen (Arachnida: Opiliones: Laniatores)

After tumultuous revisions to the family-level systematics of Laniatores (the armored harvestmen), the basally branching family Phalangodidae presently bears a disjunct and irregular distribution, attributed to the fragmentation of Pangea. One of the curious lineages assigned to Phalangodidae is the monotypic Israeli genus Haasus, the only Laniatores species that occurs in Israel, and whose presence in the Levant has been inferred to result from biogeographic connectivity with Eurasia. Recent surveys of Israeli caves have also yielded a new troglobitic morphospecies of Haasus. Here, we describe this new species as Haasus naasane sp. nov. So as to test the biogeographic affinity of Haasus, we sequenced DNA from both species and RNA from Haasus naasane sp. nov., to assess their phylogenetic placement. Our results showed that the new species is clearly closely related to Haasus judaeus, but Haasus itself is unambiguously nested within the largely Afrotropical family Pyramidopidae. In addition, the Japanese ‘phalangodid’ Proscotolemon sauteri was recovered as nested within the Southeast Asian family Petrobunidae. Phylogenomic placement of Haasus naasane sp. nov. in a 1550-locus matrix indicates that Pyramidopidae has an unstable position in the tree of Laniatores, with alternative partitioning of the matrix recovering high nodal support for mutually exclusive tree topologies. Exploration of phylogenetic signal showed the cause of this instability to be a considerable conflict between partitions, suggesting that the basal phylogeny of Laniatores may not yet be stable to addition of taxa. We transfer Haasus to Pyramidopidae (new familial assignment). Additionally, we transfer Proscotolemon to the family Petrobunidae (new familial assignment). Future studies on basal Laniatores phylogeny should emphasise the investigation of small-bodied and obscure groups that superficially resemble Phalangodidae.

The Opiliones tree of life: shedding light on harvestmen relationships through transcriptomics

Opiliones are iconic arachnids with a Palaeozoic origin and a diversity that reflects ancient biogeographic patterns dating back at least to the times of Pangea. Owing to interest in harvestman diversity, evolution and biogeography, their relationships have been thoroughly studied using morphology and PCR-based Sanger approaches to infer their systematic relationships. More recently, two studies utilized transcriptomics-based phylogenomics to explore their basal relationships and diversification, but sampling was limiting for understanding deep evolutionary patterns, as they lacked good taxon representation at the family level. Here, we analysed a set of the 14 existing transcriptomes with 40 additional ones generated for this study, representing approximately 80% of the extant familial diversity in Opiliones. Our phylogenetic analyses, including a set of data matrices with different gene occupancy and evolutionary rates, and using a multitude of methods correcting for a diversity of factors affecting phylogenomic data matrices, provide a robust and stable Opiliones tree of life, where most families and higher taxa are precisely placed. Our dating analyses using alternative calibration points, methods and analytical parameters provide well-resolved old divergences, consistent with ancient regionalization in Pangea in some groups, and Pangean vicariance in others. The integration of state-of-the-art molecular techniques and analyses, together with the broadest taxonomic sampling to date presented in a phylogenomic study of harvestmen, provide new insights into harvestmen interrelationships, as well as an overview of the general biogeographic patterns of this ancient arthropod group.

Ten genes and two topologies: an exploration of higher relationships in skipper butterflies (Hesperiidae)

Despite multiple attempts to infer the higher-level phylogenetic relationships of skipper butterflies (Family Hesperiidae), uncertainties in the deep clade relationships persist. The most recent phylogenetic analysis included fewer than 30% of known genera and data from three gene markers. Here we reconstruct the higher-level relationships with a rich sampling of ten nuclear and mitochondrial markers (7,726 bp) from 270 genera and find two distinct but equally plausible topologies among subfamilies at the base of the tree. In one set of analyses, the nuclear markers suggest two contrasting topologies, one of which is supported by the mitochondrial dataset. However, another set of analyses suggests mito-nuclear conflict as the reason for topological incongruence. Neither topology is strongly supported, and we conclude that there is insufficient phylogenetic evidence in the molecular dataset to resolve these relationships. Nevertheless, taking morphological characters into consideration, we suggest that one of the topologies is more likely.