Exploring Phylogenetic Relationships within Myriapoda and the Effects of Matrix Composition and Occupancy on Phylogenomic Reconstruction

The venom and telopodal defence systems of the centipede Lithobius forficatus are functionally convergent serial homologues

BACKGROUND

Evolution of novelty is a central theme in evolutionary biology, yet studying the origins of traits with an apparently discontinuous origin remains a major challenge. Venom systems are a well-suited model for the study of this phenomenon because they capture several aspects of novelty across multiple levels of biological complexity. However, while there is some knowledge on the evolution of individual toxins, not much is known about the evolution of venom systems as a whole. One way of shedding light on the evolution of new traits is to investigate less specialised serial homologues, i.e. repeated traits in an organism that share a developmental origin. This approach can be particularly informative in animals with repetitive body segments, such as centipedes.

RESULTS

Here, we investigate morphological and biochemical aspects of the defensive telopodal glandular organs borne on the posterior legs of venomous stone centipedes (Lithobiomorpha), using a multimethod approach, including behavioural observations, comparative morphology, proteomics, comparative transcriptomics and molecular phylogenetics. We show that the anterior venom system and posterior telopodal defence system are functionally convergent serial homologues, where one (telopodal defence) represents a model for the putative early evolutionary state of the other (venom). Venom glands and telopodal glandular organs appear to have evolved from the same type of epidermal gland (four-cell recto-canal type) and while the telopodal defensive secretion shares a great degree of compositional overlap with centipede venoms in general, these similarities arose predominantly through convergent recruitment of distantly related toxin-like components. Both systems are composed of elements predisposed to functional innovation across levels of biological complexity that range from proteins to glands, demonstrating clear parallels between molecular and morphological traits in the properties that facilitate the evolution of novelty.

CONCLUSIONS

The evolution of the lithobiomorph telopodal defence system provides indirect empirical support for the plausibility of the hypothesised evolutionary origin of the centipede venom system, which occurred through functional innovation and gradual specialisation of existing epidermal glands. Our results thus exemplify how continuous transformation and functional innovation can drive the apparent discontinuous emergence of novelties on higher levels of biological complexity.

The structural and functional modularity of ovarian follicle epithelium in the pill-millipede Hyleoglomeris japonica Verhoeff, 1936 (Diplopoda: Glomerida: Glomeridae)

Ovarian somatic tissues typically surround developing oocytes and play a crucial role in oogenesis across various metazoans, often displaying structural properties specific to their functions. However, there is an absence of evident structural modularity in the follicle epithelium of Myriapoda. We report here two structurally and developmentally distinct domains within the follicle epithelium of the Japanese pill millipede, Hyleoglomeris japonica. The follicle epithelium of H. japonica exhibits a thick cell mass at the apex of the follicle. These cells harbor abundant rough endoplasmic reticulum, mitochondria, Golgi complexes, and numerous microvilli, indicative of synthetic/secretory activities. Moreover, their height increases as oogenesis progresses. In contrast, another region of the epithelium lacks these features. Our findings highlight the presence of structural and functional modularity in the follicle epithelium of H. japonica. We suggest classifying the follicle epithelium of Myriapoda into three types: homogenous epithelia with enhanced synthetic activities, homogenous epithelia with diminished such activities, and heterogeneous epithelia with varying synthetic activities. These findings prompt a reevaluation of the nature of ovarian somatic tissues in Myriapoda as well as in Arthropoda.

A new genus of scutigerid centipede from southern South America with the description of two new species and an updated molecular phylogeny of the myriapod order Scutigeromorpha (Myriapoda: Chilopoda)

Scutigeromorph centipedes are conspicuous, yet often ignored myriapods for which little work has been conducted in southern South America. After examining recent and museum collections from Chile and Argentina, two new species of generic uncertainty were identified. A new genus of scutigerid centipede, Edgethreua , is therefore described with two new species, E. chilensis from Central Chile (type species of the genus) and E. goloboffi from Argentinian Patagonia. The new genus is characterised by the presence of scattered setiform bristles with short paired spines and the absence of simple spinulae and spines on all stomatotergites, the presence of a single spine-bristle in the prefemur of the second maxilla, a patch of cuticular ridges and pores surrounding the sensilla of the proximal labral portion of the epipharynx, the morphology of the sensilla of the distal patch of the hypopharynx and the morphology of the female gonopods. A phylogenetic analysis of the new species using two nuclear ribosomal RNA genes (18S and 28S rRNA), two mitochondrial ribosomal RNA genes (12S and 16S rRNA) and the mitochondrial protein-encoding gene cytochrome c oxidase subunit I show that the new genus does not cluster with any other described genus of scutigeromorph represented in molecular phylogenies. The data indicate that the new genus is probably sister group to a clade including the genera Lassophora , Ballonema and the subfamily Thereuoneminae, although one analysis suggests a position as sister group to Scutigerinae. ZooBank: urn:lsid:zoobank.org:pub:A4D453F3-9031-4E21-84C7-87F16C07AD51.

Inference of phylogenetic trees directly from raw sequencing reads using Read2Tree

Current methods for inference of phylogenetic trees require running complex pipelines at substantial computational and labor costs, with additional constraints in sequencing coverage, assembly and annotation quality, especially for large datasets. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes and bypasses traditional steps in phylogeny inference, such as genome assembly, annotation and all-versus-all sequence comparisons, while retaining accuracy. In a benchmark encompassing a broad variety of datasets, Read2Tree is 10-100 times faster than assembly-based approaches and in most cases more accurate-the exception being when sequencing coverage is high and reference species very distant. Here, to illustrate the broad applicability of the tool, we reconstruct a yeast tree of life of 435 species spanning 590 million years of evolution. We also apply Read2Tree to >10,000 Coronaviridae samples, accurately classifying highly diverse animal samples and near-identical severe acute respiratory syndrome coronavirus 2 sequences on a single tree. The speed, accuracy and versatility of Read2Tree enable comparative genomics at scale.

A new species of Polydesmidae (Myriapoda, Diplopoda, Polydesmida) from the mid-Cretaceous Burmese amber

A new species of the flat-backed millipede family Polydesmidae, Propolydesmus cretaceus sp. nov., is described from the mid-Cretaceous Burmese amber. Detailed morphological characters are provided on the basis of 3 well-preserved adults (2 males and 1 females), 2 female subadults, and 14 additional specimens, using confocal laser scanning microscopy (CLSM) and X-ray micro-computed tomography (CT) with computer-aided 3D-reconstructions. The new species can be placed in the extant genus Propolydesmus. It majorly differs from other congeners due to the stronger gonopod acropodite and the presence of a bipartite gonopod exomerite.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Transcriptomic data recover a new superfamily-level phylogeny of Cucujiformia (Coleoptera, Polyphaga)

Cucujiformia, the largest taxon in the order Coleoptera, exhibits extraordinary morphological, ecological, and behavioral diversity. This infraorder is currently divided into seven superfamilies, but considerably incongruent relationships among superfamilies have been reported by recent phylogenomic studies. Here, we combined the 21 newly sequenced transcriptomes representing six superfamilies with nine previously published cucujiform genomes/transcriptomes to elucidate the phylogeny and evolution of Cucujiformia. The monophyly of each of five superfamilies were consistently supported by all phylogenetic analyses based on the twelve datasets (matrix occupancy, amino acid and nucleotide data) and the two analytical methods (maximum likelihood method and Bayesian inference). Both the amino acid datasets and the RY recoded nucleotide datasets recovered the monophyly of Cucujoidea. Topology test results statistically supported the following robust superfamily-level phylogeny in Cucujiformia: (Coccinelloidea, (Cleroidea, (Tenebrionoidea, (Cucujoidea, (Chrysomeloidea, Curculionoidea))))). Our divergence time analyses recovered a Permian origin of Cucujiformia and a Jurassic-Cretaceous origin of most superfamilies. The diversification of phytophagous beetles that occurred in the Cretaceous can be attributed to its co-evolution with angiosperms, supporting the hypothesis of a Cretaceous Terrestrial Revolution.

Re-evaluating and dating myriapod diversification with phylotranscriptomics under a regime of dense taxon sampling

Recent transcriptomic studies of myriapod phylogeny have been based on relatively small datasets with <40 myriapod terminals and variably supported or contradicted the traditional morphological groupings of Progoneata and Dignatha. Here we amassed a large dataset of 104 myriapod terminals, including multiple species for each of the four myriapod classes. Across the tree, most nodes are stable and well supported. Most analyses across a range of gene occupancy levels provide moderate to strong support for a deep split of Myriapoda into Symphyla + Pauropoda (=Edafopoda) and an uncontradicted grouping of Chilopoda + Diplopoda (=Pectinopoda nov.), as in other recent transcriptome-based analyses; no analysis recovers Progoneata or Dignatha as clades. As in all recent multi-locus and phylogenomic studies, chilopod interrelationships resolve with Craterostigmus excluded from Amalpighiata rather than uniting with other centipedes with maternal brood care in Phylactometria. Diplopod ordinal interrelationships are largely congruent with morphology-based classifications. Chilognathan clades that are not invariably advocated by morphologists include Glomerida + Glomeridesmida, such that the volvation-related characters of pill millipedes may be convergent, and Stemmiulida + Polydesmida more closely allied to Juliformia than to Callipodida + Chordeumatida. The latter relationship implies homoplasy in spinnerets and contradicts Nematophora. A time-tree with nodes calibrated by 25 myriapod and six outgroup fossil terminals recovers Cambrian-Ordovician divergences for the deepest splits in Myriapoda, Edafopoda and Pectinopoda, predating the terrestrial fossil record of myriapods as in other published chronograms, whereas age estimates within Chilopoda and Diplopoda overlap with or do not appreciably predate the calibration fossils. The grouping of Chilopoda and Diplopoda is recovered in all our analyses and is formalized as Pectinopoda nov., named for the shared presence of mandibular comb lamellae. New taxonomic proposals for Chilopoda based on uncontradicted clades are Tykhepoda nov. for the three blind families of Scolopendromorpha that share a "sieve-type" gizzard, and Taktikospina nov. for Scolopendromorpha to the exclusion of Mimopidae.

Read2Tree: scalable and accurate phylogenetic trees from raw reads

The inference of phylogenetic trees is foundational to biology. However, state-of-the-art phylogenomics requires running complex pipelines, at significant computational and labour costs, with additional constraints in sequencing coverage, assembly and annotation quality. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes. In a benchmark encompassing a broad variety of datasets, our assembly-free approach was 10-100x faster than conventional approaches, and in most cases more accurate-the exception being when sequencing coverage was high and reference species very distant. To illustrate the broad applicability of the tool, we reconstructed a yeast tree of life of 435 species spanning 590 million years of evolution. Applied to Coronaviridae samples, Read2Tree accurately classified highly diverse animal samples and near-identical SARS-CoV-2 sequences on a single tree-thereby exhibiting remarkable breadth and depth. The speed, accuracy, and versatility of Read2Tree enables comparative genomics at scale.

Was There a Cambrian Explosion on Land? The Case of Arthropod Terrestrialization

Arthropods, the most diverse form of macroscopic life in the history of the Earth, originated in the sea. Since the early Cambrian, at least ~518 million years ago, these animals have dominated the oceans of the world. By the Silurian-Devonian, the fossil record attests to arthropods becoming the first animals to colonize land, However, a growing body of molecular dating and palaeontological evidence suggests that the three major terrestrial arthropod groups (myriapods, hexapods, and arachnids), as well as vascular plants, may have invaded land as early as the Cambrian-Ordovician. These dates precede the oldest fossil evidence of those groups and suggest an unrecorded continental "Cambrian explosion" a hundred million years prior to the formation of early complex terrestrial ecosystems in the Silurian-Devonian. We review the palaeontological, phylogenomic, and molecular clock evidence pertaining to the proposed Cambrian terrestrialization of the arthropods. We argue that despite the challenges posed by incomplete preservation and the scarcity of early Palaeozoic terrestrial deposits, the discrepancy between molecular clock estimates and the fossil record is narrower than is often claimed. We discuss strategies for closing the gap between molecular clock estimates and fossil data in the evolution of early ecosystems on land.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Two complete mitochondrial genomes in Scolopendra and a comparative analysis of tRNA rearrangements in centipedes

BACKGROUND

Centipedes are one of the oldest terrestrial arthropods belonging to the sub phylum Myriapoda. With the expansion of our understanding of the application of the two centipedes Scolopendra morsitans and Scolopendra hainanum, belonging to the order Scolopendromorpha, an exhaustive classification was required. Although consensus has been reached on the phylogeny of Chilopoda based on morphological traits, recent analyses based on molecular data exhibited differences in results.

METHODS AND RESULTS

The mitochondrial genome sequences of S. morsitans and S. hainanum were obtained by next-generation sequencing. S. morsitans contains 13 PCGs, two rRNAs, 11 tRNAs, and one CR. whereas S. hainanum contains 12 PCGs, of which ATP8 remains unpredicted, two rRNAs, 14 tRNAs, and one CR. An obvious tRNA rearrangement was found in the genus Scolopendra. S. morsitans exhibited a loss of trnW, trnC, trnI, trnK, trnD, trnA, trnN, trnQ, trnF, trnT, trnS, trnL, and trnV, and a repeat of trnR and trnL. S. hainanum exhibited a loss of trnQ, trnC, trnW, trnI, trnD, trnQ, trnP, and trnV. Phylogenetic analyses of centipedes based on 12 PCGs supported the sister relationship between the orders Geophilomorpha and Lithobiomorpha and a close relationship between Scolopendra dehaani and S. hainanum.

CONCLUSIONS

The new mitogenomes determined in this study provide new genomic resources for gene rearrangements and contribute to the understanding of the evolution of gene rearrangement in Chilopoda.

A phylogenomic backbone for gastropod molluscs

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from twelve taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades (i) Vetigastropoda s.l. (including Pleurotomariida + Neomphaliones) and (ii) Neritimorpha + (Caenogastropoda + Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia.

Myriapod genomes reveal ancestral horizontal gene transfer and hormonal gene loss in millipedes

Animals display a fascinating diversity of body plans. Correspondingly, genomic analyses have revealed dynamic evolution of gene gains and losses among animal lineages. Here we sequence six new myriapod genomes (three millipedes, three centipedes) at key phylogenetic positions within this major but understudied arthropod lineage. We combine these with existing genomic resources to conduct a comparative analysis across all available myriapod genomes. We find that millipedes generally have considerably smaller genomes than centipedes, with the repeatome being a major contributor to genome size, driven by independent large gains of transposons in three centipede species. In contrast to millipedes, centipedes gained a large number of gene families after the subphyla diverged, with gains contributing to sensory and locomotory adaptations that facilitated their ecological shift to predation. We identify distinct horizontal gene transfer (HGT) events from bacteria to millipedes and centipedes, with no identifiable HGTs shared among all myriapods. Loss of juvenile hormone O-methyltransferase, a key enzyme in catalysing sesquiterpenoid hormone production in arthropods, was also revealed in all millipede lineages. Our findings suggest that the rapid evolution of distinct genomic pathways in centipede and millipede lineages following their divergence from the myriapod ancestor, was shaped by differing ecological pressures.

Myriapods play an important ecological role in soil and forest ecosystems. Here the authors analyse nine myriapod genomes, showing rapid evolution of distinct genomic pathways in centipede and millipede lineages, shaped by differing ecological pressures.

Novel mitochondrial gene rearrangements pattern in the millipede Polydesmus sp. GZCS-2019 and phylogenetic analysis of the Myriapoda

The subphylum Myriapoda included four extant classes (Chilopoda, Symphyla, Diplopoda, and Pauropoda). Due to the limitation of taxon sampling, the phylogenetic relationships within Myriapoda remained contentious, especially for Diplopoda. Herein, we determined the complete mitochondrial genome of Polydesmus sp. GZCS-2019 (Myriapoda: Polydesmida) and the mitochondrial genomes are circular molecules of 15,036 bp, with all genes encoded on + strand. The A+T content is 66.1%, making the chain asymmetric, and exhibits negative AT-skew (-0.236). Several genes rearrangements were detected and we propose a new rearrangement model: "TD (N\R) L + C" based on the genome-scale duplication + (non-random/random) loss + recombination. Phylogenetic analyses demonstrated that Chilopoda and Symphyla both were monophyletic group, whereas Pauropoda was embedded in Diplopoda to form the Dignatha. Divergence time showed the first split of Myriapoda occurred between the Chilopoda and other classes (Wenlock period of Silurian). We combine phylogenetic analysis, divergence time, and gene arrangement to yield valuable insights into the evolutionary history and classification relationship of Myriapoda and these results support a monophyletic Progoneata and the relationship (Chilopoda + (Symphyla + (Diplopoda + Pauropoda))) within myriapod. Our results help to better explain the gene rearrangement events of the invertebrate mitogenome and lay the foundation for further phylogenetic study of Myriapoda.

Integrated phylogenomics and fossil data illuminate the evolution of beetles

Beetles constitute the most biodiverse animal order with over 380 000 described species and possibly several million more yet unnamed. Recent phylogenomic studies have arrived at considerably incongruent topologies and widely varying estimates of divergence dates for major beetle clades. Here, we use a dataset of 68 single-copy nuclear protein-coding (NPC) genes sampling 129 out of the 193 recognized extant families as well as the first comprehensive set of fully justified fossil calibrations to recover a refined timescale of beetle evolution. Using phylogenetic methods that counter the effects of compositional and rate heterogeneity, we recover a topology congruent with morphological studies, which we use, combined with other recent phylogenomic studies, to propose several formal changes in the classification of Coleoptera: Scirtiformia and Scirtoidea sensu nov., Clambiformia ser. nov. and Clamboidea sensu nov., Rhinorhipiformia ser. nov., Byrrhoidea sensu nov., Dryopoidea stat. res., Nosodendriformia ser. nov. and Staphyliniformia sensu nov., and Erotyloidea stat. nov., Nitiduloidea stat. nov. and Cucujoidea sensu nov., alongside changes below the superfamily level. Our divergence time analyses recovered a late Carboniferous origin of Coleoptera, a late Palaeozoic origin of all modern beetle suborders and a Triassic-Jurassic origin of most extant families, while fundamental divergences within beetle phylogeny did not coincide with the hypothesis of a Cretaceous Terrestrial Revolution.

Reinvestigating the phylogeny of Myriapoda with more extensive taxon sampling and novel genetic perspective

Background

There have been extensive debates on the interrelationships among the four major classes of Myriapoda-Chilopoda, Symphyla, Diplopoda, and Pauropoda. The core controversy is the position of Pauropoda; that is, whether it should be grouped with Symphyla or Diplopoda as a sister group. Two recent phylogenomic studies separately investigated transcriptomic data from 14 and 29 Myriapoda species covering all four groups along with outgroups, and proposed two different topologies of phylogenetic relationships.

Methods

Building on these studies, we extended the taxon sampling by investigating 39 myriapods and integrating the previously available data with three new transcriptomic datasets generated in this study. Our analyses present the phylogenetic relationships among the four major classes of Myriapoda with a more abundant taxon sampling and provide a new perspective to investigate the above-mentioned question, where visual genes' identification were conducted. We compared the appearance pattern of genes, grouping them according to their classes and the visual pathways involved. Positive selection was detected for all identified visual genes between every pair of 39 myriapods, and 14 genes showed positive selection among 27 pairs.

Results

From the results of phylogenomic analyses, we propose that Symphyla is a sister group of Pauropoda. This stance has also received strong support from tree inference and topology tests.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Incorporating Hierarchical Characters into Phylogenetic Analysis

Popular optimality criteria for phylogenetic trees focus on sequences of characters that are applicable to all the taxa. As studies grow in breadth, it can be the case that some characters are applicable for a portion of the taxa and inapplicable for others. Past work has explored the limitations of treating inapplicable characters as missing data, noting that this strategy may favor trees where interval nodes are assigned impossible states, where the arrangement of taxa within subclades is unduly influenced by variation in distant parts of the tree, and/or where taxa that otherwise share most primary characters are grouped distantly. Approaches that avoid the first two problems have recently been proposed. Here, we propose an alternative approach which avoids all three problems. We focus on data matrices that use reductive coding of traits, that is, explicitly incorporate the innate hierarchy induced by inapplicability, and as such our approach extend to hierarchical characters, in general. In the spirit of maximum parsimony, the proposed criterion seeks the phylogenetic tree with the minimal changes across any tree branch, but where changes are defined in terms of dissimilarity metrics that weigh the affects of inapplicable characters. The approach can accommodate binary, multistate, ordered, unordered, and polymorphic characters. We give a polynomial-time algorithm, inspired by Fitch's algorithm, to score trees under a family of dissimilarity metrics, and prove its correctness. We show that the resulting optimality criteria is computationally hard, by reduction to the NP-hardness of the maximum parsimony optimality criteria. We demonstrate our approach using synthetic and empirical data sets and compare the results with other recently proposed methods for choosing optimal phylogenetic trees when the data includes hierarchical characters.

Phylogenetic analyses suggest centipede venom arsenals were repeatedly stocked by horizontal gene transfer

Venoms have evolved over a hundred times in animals. Venom toxins are thought to evolve mostly by recruitment of endogenous proteins with physiological functions. Here we report phylogenetic analyses of venom proteome-annotated venom gland transcriptome data, assisted by genomic analyses, to show that centipede venoms have recruited at least five gene families from bacterial and fungal donors, involving at least eight horizontal gene transfer events. These results establish centipedes as currently the only known animals with venoms used in predation and defence that contain multiple gene families derived from horizontal gene transfer. The results also provide the first evidence for the implication of horizontal gene transfer in the evolutionary origin of venom in an animal lineage. Three of the bacterial gene families encode virulence factors, suggesting that horizontal gene transfer can provide a fast track channel for the evolution of novelty by the exaptation of bacterial weapons into animal venoms.

The Development of Arthropod Segmentation Across the Embryonic/Post-embryonic Divide – An Evolutionary Perspective

In many arthropods, the appearance of new segments and their differentiation are not completed by the end of embryogenesis but continue, in different form and degree, well after hatching, in some cases up to the last post-embryonic molt. Focusing on the segmentation process currently described as post-embryonic segment addition (or, anamorphosis), we revise here the current knowledge and discuss it in an evolutionary framework which involves data from fossils, comparative morphology of extant taxa and gene expression. We advise that for a better understanding of the developmental changes underlying the evolution of arthropod segmentation, some key concepts should be applied in a critical way. These include the notion of the segment as a body block and the idea that hatching represents a well-defined divide, shared by all arthropods, between two contrasting developmental phases, embryonic vs. post-embryonic. This eventually reveals the complexity of the developmental processes occurring across hatching, which can evolve in different directions and with a different pace, creating the observed vagueness of the embryonic/post-embryonic divide.

The tracheal system of scutigeromorph centipedes and the evolution of respiratory systems of myriapods

The tracheal system of scutigeromorph centipedes (Chilopoda) is special, as it consists of dorsally arranged unpaired spiracles. In this study, we investigate the tracheal systems of five different scutigeromorph species. They are strikingly similar to each other but depict unique characters compared to the tracheal systems of pleurostigmophoran centipedes, which has engendered an ongoing debate over a single versus independent origin of tracheal systems in Chilopoda. Up to now, only the respiratory system of Scutigera coleoptrata was investigated intensively using LM-, TEM-, and SEM-techniques. We supplement this with data for species from all three families of Scutigeromorpha. These reveal interspecific differences in atrial width and the shape and branching pattern of the tracheal tubules. Further, we investigated the tracheal system of Scutigera coleoptrata with three additional techniques: light sheet microscopy, microCT and synchrotron radiation based microCT analysis. This set of techniques allows a comparison between fresh versus fixed and dried material. The question of a unique vs. multiple origin of tracheal systems in centipedes and in Myriapoda as a whole is discussed with regard to their structural similarities and differences and the presence of hemocyanin as an oxygen carrier. We used morphological and molecular data and the fossil record to evaluate the alternative hypotheses.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

Life after the mother's hug: Late post-embryonic development of Cryptops parisi (Chilopoda: Scolopendromorpha: Cryptopidae)

Here we describe in detail the late post-embryonic development of the common European scolopendromorph centipede Cryptops parisi. Canonical variate analyses of two groups of external morphological characters, viz., cephalic capsule characters (head length, length of the anterior and posterior paramedian cephalic sutures) and coxopleuron surface characters (number of pores in the coxal pore-field, number of setae on the posterior coxopleuron edge, their number on the coxal pore-field, and their number posterior to the coxal pore-field) were conducted on a large sample of specimens collected from two localities in Serbia. Ten free-living stages are recognized: three pre-adult stages (adolescens I, II, and III) and seven adult stages (one maturus junior stage, four maturus, and two maturus senior stages). The fourth late post-embryonic stage is the first mature stage in both sexes. Sexual dimorphism in the aforementioned characters was not observed. Morphological variation of coxopleuron characters was more informative for the discrimination of developmental stages in Cryptops than the morphological variation of cephalic capsule characters.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola)

BACKGROUND

The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations.

RESULTS

Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera.

CONCLUSION

Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

Shedding light: a phylotranscriptomic perspective illuminates the origin of photosymbiosis in marine bivalves

BACKGROUND

Photosymbiotic associations between metazoan hosts and photosynthetic dinoflagellates are crucial to the trophic and structural integrity of many marine ecosystems, including coral reefs. Although extensive efforts have been devoted to study the short-term ecological interactions between coral hosts and their symbionts, long-term evolutionary dynamics of photosymbiosis in many marine animals are not well understood. Within Bivalvia, the second largest class of mollusks, obligate photosymbiosis is found in two marine lineages: the giant clams (subfamily Tridacninae) and the heart cockles (subfamily Fraginae), both in the family Cardiidae. Morphologically, giant clams show relatively conservative shell forms whereas photosymbiotic fragines exhibit a diverse suite of anatomical adaptations including flattened shells, leafy mantle extensions, and lens-like microstructural structures. To date, the phylogenetic relationships between these two subfamilies remain poorly resolved, and it is unclear whether photosymbiosis in cardiids originated once or twice.

RESULTS

In this study, we establish a backbone phylogeny for Cardiidae utilizing RNASeq-based transcriptomic data from Tridacninae, Fraginae and other cardiids. A variety of phylogenomic approaches were used to infer the relationship between the two groups. Our analyses found conflicting gene signals and potential rapid divergence among the lineages. Overall, results support a sister group relationship between Tridacninae and Fraginae, which diverged during the Cretaceous. Although a sister group relationship is recovered, ancestral state reconstruction using maximum likelihood-based methods reveals two independent origins of photosymbiosis, one at the base of Tridacninae and the other within a symbiotic Fraginae clade.

CONCLUSIONS

The newly revealed common ancestry between Tridacninae and Fraginae brings a possibility that certain genetic, metabolic, and/or anatomical exaptations existed in their last common ancestor, which promoted both lineages to independently establish photosymbiosis, possibly in response to the modern expansion of reef habitats.

The complete mitochondrial genome sequence of Scolopendra mutilans L. Koch, 1878 (Scolopendromorpha, Scolopendridae), with a comparative analysis of other centipede genomes

Scolopendramutilans L. Koch, 1878 is an important Chinese animal with thousands of years of medicinal history. However, the genomic information of this species is limited, which hinders its further application. Here, the complete mitochondrial genome (mitogenome) of S.mutilans was sequenced and assembled by next-generation sequencing. The genome is 15,011 bp in length, consisting of 13 protein-coding genes (PCGs), 14 tRNA genes, and two rRNA genes. Most PCGs start with the ATN initiation codon, and all PCGs have the conventional stop codons TAA and TAG. The S.mutilans mitogenome revealed nine simple sequence repeats (SSRs), and an obviously lower GC content compared with other seven centipede mitogenomes previously sequenced. After analysis of homologous regions between the eight centipede mitogenomes, the S.mutilans mitogenome further showed clear genomic rearrangements. The phylogenetic analysis of eight centipedes using 13 conserved PCG genes was finally performed. The phylogenetic reconstructions showed Scutigeromorpha as a separate group, and Scolopendromorpha in a sister-group relationship with Lithobiomorpha and Geophilomorpha. Collectively, the S.mutilans mitogenome provided new genomic resources, which will improve its medicinal research and applications in the future.

Trait differentiation and modular toxin expression in palm-pitvipers

Background

Modularity is the tendency for systems to organize into semi-independent units and can be a key to the evolution and diversification of complex biological systems. Snake venoms are highly variable modular systems that exhibit extreme diversification even across very short time scales. One well-studied venom phenotype dichotomy is a trade-off between neurotoxicity versus hemotoxicity that occurs through the high expression of a heterodimeric neurotoxic phospholipase A₂ (PLA₂) or snake venom metalloproteinases (SVMPs). We tested whether the variation in these venom phenotypes could occur via variation in regulatory sub-modules through comparative venom gland transcriptomics of representative Black-Speckled Palm-Pitvipers (Bothriechis nigroviridis) and Talamancan Palm-Pitvipers (B. nubestris).

Results

We assembled 1517 coding sequences, including 43 toxins for B. nigroviridis and 1787 coding sequences including 42 toxins for B. nubestris. The venom gland transcriptomes were extremely divergent between these two species with one B. nigroviridis exhibiting a primarily neurotoxic pattern of expression, both B. nubestris expressing primarily hemorrhagic toxins, and a second B. nigroviridis exhibiting a mixed expression phenotype. Weighted gene coexpression analyses identified six submodules of transcript expression variation, one of which was highly associated with SVMPs and a second which contained both subunits of the neurotoxic PLA₂ complex. The sub-module association of these toxins suggest common regulatory pathways underlie the variation in their expression and is consistent with known patterns of inheritance of similar haplotypes in other species. We also find evidence that module associated toxin families show fewer gene duplications and transcript losses between species, but module association did not appear to affect sequence diversification.

Conclusion

Sub-modular regulation of expression likely contributes to the diversification of venom phenotypes within and among species and underscores the role of modularity in facilitating rapid evolution of complex traits.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Transcriptomics illuminate the phylogenetic backbone of tiger beetles

Phylogenomics is progressing rapidly, allowing large strides forward into our understanding of the tree of life. In this study, we generated transcriptomes from ethanol-preserved specimens of 13 tiger beetle species (Coleoptera: Cicindelinae) and one Scaritinae outgroup. From these 14 transcriptomes and seven publicly available transcriptomes, we recovered an average of 2538 loci for phylogenetic analysis. We constructed an evolutionary tree of tiger beetles to examine deep-level relationships and examined the extent to which the composition of the dataset, missing data, gene tree inconsistency and codon position saturation impacted phylogenetic accuracy. Ethanol-preserved specimens yielded similar numbers of loci to specimens originally preserved in costly reagents, showcasing more flexibility in transcriptomics than anticipated. The number of loci and gene tree inconsistency had less impact on downstream results than third codon position saturation and missing data. Our results recovered tiger beetles as sister to Carabidae with strong support, confirming their taxonomic status as an independent family within Adephaga. Within tiger beetles, phylogenetic relationships were robust across all nodes. This new phylogenomic backbone represents a useful framework for future endeavours in tiger beetle systematics and serves as a starting point for the development of less costly target capture toolkits to expand the taxonomic breadth of the future tiger beetle tree of life.

Fossils from South China redefine the ancestral euarthropod body plan

BACKGROUND

Early Cambrian Lagerstätten from China have greatly enriched our perspective on the early evolution of animals, particularly arthropods. However, recent studies have shown that many of these early fossil arthropods were more derived than previously thought, casting uncertainty on the ancestral euarthropod body plan. In addition, evidence from fossilized neural tissues conflicts with external morphology, in particular regarding the homology of the frontalmost appendage.

RESULTS

Here we redescribe the multisegmented megacheirans Fortiforceps and Jianfengia and describe Sklerolibyon maomima gen. et sp. nov., which we place in Jianfengiidae, fam. nov. (in Megacheira, emended). We find that jianfengiids show high morphological diversity among megacheirans, both in trunk ornamentation and head anatomy, which encompasses from 2 to 4 post-frontal appendage pairs. These taxa are also characterized by elongate podomeres likely forming seven-segmented endopods, which were misinterpreted in their original descriptions. Plesiomorphic traits also clarify their connection with more ancestral taxa. The structure and position of the "great appendages" relative to likely sensory antero-medial protrusions, as well as the presence of optic peduncles and sclerites, point to an overall homology with the anterior head of radiodontans. This is confirmed by our Bayesian phylogeny, which places jianfengiids as the basalmost euarthropods, paraphyletic with other megacheirans, and in contiguity with isoxyids and radiodontans.

CONCLUSIONS

Sklerolibyon and other jianfengiids expand the disparity of megacheirans and suggest that the common euarthropod ancestor possessed a remarkable phenotypic variability associated with the externalized cephalon, as well as endopods that were already heptopodomerous, which differs from previous hypotheses and observations. These animals also demonstrate that the frontalmost pair of arthrodized appendage is homologous between radiodontans and megacheirans, refuting the claim that the radiodontan frontal appendages evolved into the euarthropod labrum, and questioning its protocerebral identity. This evidence based on external anatomy now constitutes a solid benchmark upon which we should address issues of homology, with the help of carefully examined palaeoneurological data.

Histaminergic interneurons in the ventral nerve cord: assessment of their value for Euarthropod phylogeny

Despite numerous approaches to the resolution of euarthropod phylogeny, mainly based on modern sequence information and traditional external morphology, the resulting hypotheses are often contradictory and leave many questions about euarthropod evolution unanswered. The comparison of developmental and structural aspects of the nervous system has shown to be a valuable contribution to the assessment of current phylogenetic hypotheses. One promising approach for the generation of new character sets is the morphology of transmitter systems and the discovery of individually identifiable neurons, which allow phylogenetic comparisons on the single cell level. In this context, the serotonin transmitter system has been investigated to a considerable degree. Studies to date have yielded important stimuli to our understanding of euarthropod relationships and the evolution of their nervous systems. However, data on other transmitter systems remain fragmented, and their value with respect to phylogenetic questions remains speculative. The biogenic amine histamine is a promising transmitter; a substantial amount of data has been reported in the literature and the homology of some histaminergic neurons has been suggested. Here, we present a comprehensive review of histaminergic neurons in the ventral nerve cord of Euarthropoda. Using immunocytochemical labeling of histamine combined with confocal laser-scanning microscopy, we investigated the transmitter system in phylogenetically relevant taxa, such as Zygentoma, Remipedia, Diplopoda, and Arachnida. By reconstructing ground patterns, we evaluated the significance of this specific character set for euarthropod phylogeny. With this approach, we identified a set of neurons, which can be considered homologous within the respective major taxon. In conclusion, the histaminergic system contains useful information for our understanding of euarthropod phylogeny, supporting the proposed clades Tetraconata and Mandibulata. Furthermore, this character set has considerable potential to help resolve relationships within the major clades at a deeper level of taxonomy, due to the considerable variability in neurite morphology.

Parallel Evolution of Complex Centipede Venoms Revealed by Comparative Proteotranscriptomic Analyses

Centipedes are among the most ancient groups of venomous predatory arthropods. Extant species belong to five orders, but our understanding of the composition and evolution of centipede venoms is based almost exclusively on one order, Scolopendromorpha. To gain a broader and less biased understanding we performed a comparative proteotranscriptomic analysis of centipede venoms from all five orders, including the first venom profiles for the orders Lithobiomorpha, Craterostigmomorpha, and Geophilomorpha. Our results reveal an astonishing structural diversity of venom components, with 93 phylogenetically distinct protein and peptide families. Proteomically-annotated gene trees of these putative toxin families show that centipede venom composition is highly dynamic across macroevolutionary timescales, with numerous gene duplications as well as functional recruitments and losses of toxin gene families. Strikingly, not a single family is found in the venoms of representatives of all five orders, with 67 families being unique for single orders. Ancestral state reconstructions reveal that centipede venom originated as a simple cocktail comprising just four toxin families, with very little compositional evolution happening during the approximately 50 My before the living orders had diverged. Venom complexity then increased in parallel within the orders, with scolopendromorphs evolving particularly complex venoms. Our results show that even venoms composed of toxins evolving under the strong constraint of negative selection can have striking evolutionary plasticity on the compositional level. We show that the functional recruitments and losses of toxin families that shape centipede venom arsenals are not concentrated early in their evolutionary history, but happen frequently throughout.

Evolutionary emergence of Hairless as a novel component of the Notch signaling pathway

Suppressor of Hairless [Su(H)], the transcription factor at the end of the Notch pathway in Drosophila, utilizes the Hairless protein to recruit two co-repressors, Groucho (Gro) and C-terminal Binding Protein (CtBP), indirectly. Hairless is present only in the Pancrustacea, raising the question of how Su(H) in other protostomes gains repressive function. We show that Su(H) from a wide array of arthropods, molluscs, and annelids includes motifs that directly bind Gro and CtBP; thus, direct co-repressor recruitment is ancestral in the protostomes. How did Hairless come to replace this ancestral paradigm? Our discovery of a protein (S-CAP) in Myriapods and Chelicerates that contains a motif similar to the Su(H)-binding domain in Hairless has revealed a likely evolutionary connection between Hairless and Metastasis-associated (MTA) protein, a component of the NuRD complex. Sequence comparison and widely conserved microsynteny suggest that S-CAP and Hairless arose from a tandem duplication of an ancestral MTA gene.

Combining morphological and genomic evidence to resolve species diversity and study speciation processes of the Pallenopsis patagonica (Pycnogonida) species complex

BACKGROUND

Pallenopsis patagonica (Hoek, 1881) is a morphologically and genetically variable sea spider species whose taxonomic classification is challenging. Currently, it is considered as a species complex including several genetic lineages, many of which have not been formally described as species. Members of this species complex occur on the Patagonian and Antarctic continental shelves as well as around sub-Antarctic islands. These habitats have been strongly influenced by historical large-scale glaciations and previous studies suggested that communities were limited to very few refugia during glacial maxima. Therefore, allopatric speciation in these independent refugia is regarded as a common mechanism leading to high biodiversity of marine benthic taxa in the high-latitude Southern Hemisphere. However, other mechanisms such as ecological speciation have rarely been considered or tested. Therefore, we conducted an integrative morphological and genetic study on the P. patagonica species complex to i) resolve species diversity using a target hybrid enrichment approach to obtain multiple genomic markers, ii) find morphological characters and analyze morphometric measurements to distinguish species, and iii) investigate the speciation processes that led to multiple lineages within the species complex.

RESULTS

Phylogenomic results support most of the previously reported lineages within the P. patagonica species complex and morphological data show that several lineages are distinct species with diagnostic characters. Two lineages are proposed as new species, P. aulaeturcarum sp. nov. Dömel & Melzer, 2019 and P. obstaculumsuperavit sp. nov. Dömel, 2019, respectively. However, not all lineages could be distinguished morphologically and thus likely represent cryptic species that can only be identified with genetic tools. Further, morphometric data of 135 measurements showed a high amount of variability within and between species without clear support of adaptive divergence in sympatry.

CONCLUSIONS

We generated an unprecedented molecular data set for members of the P. patagonica sea spider species complex with a target hybrid enrichment approach, which we combined with extensive morphological and morphometric analyses to investigate the taxonomy, phylogeny and biogeography of this group. The extensive data set enabled us to delineate species boundaries, on the basis of which we formally described two new species. No consistent evidence for positive selection was found, rendering speciation in allopatric glacial refugia as the most likely model of speciation.

OMA standalone: orthology inference among public and custom genomes and transcriptomes

Genomes and transcriptomes are now typically sequenced by individual laboratories but analyzing them often remains challenging. One essential step in many analyses lies in identifying orthologs—corresponding genes across multiple species—but this is far from trivial. The Orthologous MAtrix (OMA) database is a leading resource for identifying orthologs among publicly available, complete genomes. Here, we describe the OMA pipeline available as a standalone program for Linux and Mac. When run on a cluster, it has native support for the LSF, SGE, PBS Pro, and Slurm job schedulers and can scale up to thousands of parallel processes. Another key feature of OMA standalone is that users can combine their own data with existing public data by exporting genomes and precomputed alignments from the OMA database, which currently contains over 2100 complete genomes. We compare OMA standalone to other methods in the context of phylogenetic tree inference, by inferring a phylogeny of Lophotrochozoa, a challenging clade within the protostomes. We also discuss other potential applications of OMA standalone, including identifying gene families having undergone duplications/losses in specific clades, and identifying potential drug targets in nonmodel organisms. OMA standalone is available under the permissive open source Mozilla Public License Version 2.0.

Solovyeva E. On the taxonomic position of the enigmatic genus Tonkinodentus Schileyko, 1992 (Chilopoda, Scolopendromorpha): the first molecular data

The taxonomic position of the monotypic Vietnamese genus Tonkinodentus Schileyko, 1992 (for T.lestes Schileyko, 1992) has been considered in the light of the first obtained molecular data. Both molecular (28S rRNA) and morphological data support the position of this extraordinary eye-less genus within the family Scolopendridae Leach, 1814, a sighted clade, and thus suggests the polyphyly of blind scolopendromorphs. The species diagnosis has been amended and color images of T.lestes provided for the first time.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.

Evolutionary biogeography of the centipede genus Ethmostigmus from Peninsular India: testing an ancient vicariance hypothesis for Old World tropical diversity

BACKGROUND

Understanding the relative influence of vicariance and dispersal in shaping Old World tropical biodiversity remains a challenge. We aimed to infer the roles of these alternative biogeographic processes using a species time-tree for the centipede genus Ethmostigmus from the Old World tropics. Additionally, we explored fine-scale biogeographic patterns for an endemic radiation of Ethmostigmus from the peninsular Indian Plate (PIP), an area with complex geological and climatic history.

RESULTS

Divergence time estimates suggest that Ethmostigmus began diversifying in the Late Cretaceous, 99 (± 25) million years ago (Ma), its early biogeographic history shaped by vicariance. Members of Ethmostigmus in PIP form a monophyletic group that underwent endemic radiation in the Late Cretaceous, 72 (± 25) Ma. In contrast, a new species of Ethmostigmus from north-east India formed a clade with African/Australian species. Fine-scale biogeographic analyses in PIP predict that Indian Ethmostigmus had an ancestor in southern-central parts of the Western Ghats. This was followed by four independent dispersal events from the southern-central Western Ghats to the Eastern Ghats, and between different parts of the Western Ghats in the Cenozoic.

CONCLUSIONS

Our results are consistent with Gondwanan break-up driving the early evolutionary history of the genus Ethmostigmus. Multiple dispersal events coinciding with geo-climatic events throughout the Cenozoic shaped diversification in PIP. Ethmostigmus species in PIP are restricted to wet forests and have retained that niche throughout their diversification.

Blind scolopendrid centipedes of the genus Cormocephalus from subterranean habitats in Western Australia (Myriapoda: Scolopendromorpha: Scolopendridae)

Only a single blind species is known in the centipede family Scolopendridae, representing the monotypic genus Tonkinodentus Schileyko, 1992, from Vietnam. All of more than 400 other species have four ocelli on each side of the cephalic plate. A complex of three new blind species of the genus Cormocephalus Newport, 1844, is described from the subterranean fauna of the central Pilbara region of Western Australia. Phylogenies based on sequence data for the barcode region of COI and a concatenated matrix that also includes 12S rRNA, 28S rRNA and ITS2 unite the blind Pilbara species as a monophyletic group, albeit with moderate bootstrap support, informally named the C. sagmus species group. Cormocephalus sagmus, C. pyropygus and C. delta spp. nov. supplement 17 epigean congeners previously described from Australia. The new species are all morphologically similar, but can be distinguished using the shape and spinulation of the ultimate leg prefemur. Two additional genetically distinct lineages were recovered that are not described, owing to the specimens being immature or lacking diagnostic morphological characters. The subterranean radiation in the Pilbara is more closely related to species from forests in the south-west of Western Australia than to congeners from the arid zone. http://zoobank.org/urn:lsid:zoobank.org:pub:6F67FD31-A373-4DC5-A5FD-374D32DEE02C